comsol流体动力学模型油箱的晃动.doc

流体动力学模型：油箱的晃动 模型背景软件可以借助移动网格技术来仿真动态自由表面流本案例模型表明COMSOL Multiphysics 方程来模拟液体的流动。矩形油箱内的液体（甘油）动。模型使用不可压缩Navier-Stokes?4 。原本处于静止状态。运动由来回摆动的重力矢量来驱动，重力矢量的最大摆动角度为ALE(arbitrary 。然而，不规范的计算问题约束由于液体的表面没有，该模型为一个COMSOL Lagrangian-Eulerian)技术可以很好的解决这类问题。不仅仅是可以方便的在）应用模式，而且在移动网格求解域边界，它代表了自中建立移动网格（ALEMultiphysics也使得表面张量分析成这也允许更加精确的估计边界条件特性，比如曲率，由的边界条件。 为可能。本实例忽略了表面张量的影响。 模型几何、控制方程及边界条件 Navier-Stokes方程式来描述流体动力特性：本模型使用不可压缩????u???T????F?u?I?????u??u???u t?0??u? ???????ftsin2F?gsinmaxx ???????ft2cos?F?singmaxy2??Hz1/180,f/g?9.81ms,??4 其中 max 移动网格 为了能够用动网格来跟随模拟运动的流体，必须将移动网格耦合到流体表面的法线上去。结果是，对于这种自由表面运动，不能将移动网格耦合到流体运动的切线方向。如果这样做网格会变形过度从而导致问题不再收敛。自由表面的网格等式边界条件在这种条件下可写为： ??T?n?u?x,yn ( 即u*nx+v*ny) tt??Ty,x 表示网格的移动速度表示边界的法向矢量，其中ntt 模型数据 Parameters Name Expression Value rho 1270[kg/m^3] nu 1.49[Pa*s] phi_max (4*pi/180)[rad] freq 1[Hz] g 9.81[m/s^2] Variables Name Expression Unit rad phi phi_max*sin(2*pi*freq*t) grav_x g*sin(phi) grav_y -g*cos(phi) 分析结果要求：在t=1s，t=1.2s,t=1.4s,t=1.6s时的速度场分布情况。重力矢量的倾斜是由油箱的倾斜引起的。 追踪自由液面在X=0.5m处，时间段t∈[0,6s]的波峰高度。 建模过程 1.将所需的材料参数预先写如TXT格式的记事本文档中(包括全局定义参数Parameters，和变量Variables 2.运行COMSO4.2创建模型，选择2D——Laminar Flow(spf)层流和Moving mesh动网格——Time Dependent瞬态——生成 3.进入几何绘图界面，尺寸参数不用改变，长度单位为m,右键单击Geometry——选择Rectangle绘制矩形尺寸如图 Form a union 生成集合，4. 选择右上角Build Selected键生成 5.定义流体属性 mu 材料动力学粘度rho,定义材料物理参数，如图，材料密度Fluid Properties 1,选择． 6.添加滑动边界条件，右键选择Laminar Flow， 选择Wall(壁)，在几何模型中选择如图所示1、2、4边，在Boundary condition选项中选择Slip(滑移) 7.添加开边界，右键选择Laminar Flow， 选择Open Boundary，在几何模型中选择3边，其余保持默认设置 8.添加体积力，右键选择Laminar Flow， 选择Volume Force，选中体积力作用区域1（如图几何模型区域），在Volume force 中输入体积力Fx=rho*grav_x Fy=rho*grav_y 9.定义移动网格，添加自由变形区域，右键选择Moving Mesh选项，选择Free Deformation，并选中自由变形区域，几何模型1，如图。 10.定义网格位移，右键选择Moving Mesh选项，选择Prescribed mesh displacement,在几何模型中选择1、4边，定义x方向的位移为0（将y方向位移前的 去掉）． 11.同10步，定义2边的位移，x、y方向的位移都为0 12.定义开边界上的网格移动速度，右键选择Moving Mesh选项，选择Prescribed mesh Velocity, 在几何模型中选择3边