燃油箱燃油晃动的CFD计算.ppt

July 30, 1999 燃油箱燃油晃动的CFD计算 Advanced Contact Fasteners Training Manual 1.问题的描述 如图所示，油箱内制有四个隔板，油箱的尺寸为0.4*0.3*0.2,液体的初始液面高度为0.08m，油箱的运动规律： T=0；VX=0 T=3; VX=25 T=3.5;VX=0 T=6;VX=0 2. VOF的计算原理 VOF模型通过求解单独的动量方程和处理穿过区域的每一流体的volume fraction来模拟两种或三种不能混合的流体。典型的应用包括预测，射流断裂、流体中大泡的运动、和气液界面的稳态和瞬态处理。 在FLUENT中VOF公式通常用于计算时间依赖解，但是对于只关心稳态解的问题，它也可以执行稳态计算。稳态VOF计算是敏感的只有当你的解是独立于初始时间并且对于单相有明显的流入边界。例如，由于在旋转的杯子中自由表面的形状依赖于流体的出事水平，这样的问题必须使用time-dependent公式。另一方面，渠道内顶部有空气的水的流动和分离的空气入口可以采用steady-state公式求解 2. VOF的计算原理 VOF公式依靠的是两种或多种流体（或相）没有互相穿插（interpenetrating）这一事实。对你增加到模型里的每一附加相，就引进一个变量：即计算单元里的相的容积比率（the volume fraction of the phase）。在每个控制容积内，所有相的volume fraction的和为1。所有变量及其属性的区域被各相共享并且代表了容积平均值（volume-averaged values）,只要每一相的容积比率在每一位置是可知的。这样，在任何给定单元内的变量及其属性或者纯粹代表了一相，或者代表了相的混合，这取决于容积比率值。换句话说，在单元中，如果第q相流体的容积比率记为aq，那么下面的三个条件是可能的： 2. VOF的计算原理 液气分界面最典型的是水与大气的分界面，即自由面。由于自由面本身是运动和变形的，而且其形状常常也是一个需要求解的未知函数，因此就有一个自由面的运动学条件问题。设自由面方程为 并假定在自由面上的流体质点始终保持在自由面上，则流体质点在自由面上一点的法向速度，应该等于自由面本身在这一点的法向速度。经过一系列推导(参见文献[2])，得到自由液面运动学条件： 如果要考虑液气边界上的表面张力，则在界面两侧，两种介质的压强差与表面张力有如下关系： 2. VOF的计算原理 FLUENT中表面张力模型是由Brackbill et al[25]提出的连续表面力模型。用这个模型，VOF计算中附加的表面张力导致了动量方程中的源项。为了理解这个源项的起源，考虑沿着表面表面张力为常数的的特殊情况，那些地方只考虑垂直于界面的力。 3. 问题的求解 3.1 导入模型 3. 问题的求解 3.2 网格划分 3. 问题的求解 3.2 网格划分 3. 问题的求解 3.3 Fluent前处理 3. 问题的求解 3.3 Fluent前处理 3. 问题的求解 3.3 Fluent前处理 3. 问题的求解 3.3 Fluent前处理 3. 问题的求解 3.3 Fluent前处理 3. 问题的求解 3.3 Fluent前处理 在ANSYS FLUENT中，有两种方式用于指定瞬态网格区域及边界条件： （1）与标准profile格式相同的瞬态profile文件 （2）表格格式的瞬态profile文件。 对于这两种方式，网格域及边界条件的变化均只与时间有关。然而如果使用了in-cylinder模型，则用户可以使用crank角代替时间变量。 1、标准瞬态profile文件 一个标准瞬态profile文件格式如下： ((profile-name transient n periodic?) (field_name_1 a1 a2 a3 … an) (field_name_2 b1 b2 b3 … bn) 3. 问题的求解 3.3 Fluent前处理 注意： （1）profile-name：profile名，必须全部为小写字母，少于64个字符。 （2）transient为保留关键字，不可随便更改。 （3）n表示每一个场变量所包含的变量个数。 （4）periodic?用于标识是否使用周期，1为使用时间周期，0表示不使用时间周期。 （5）通常field_name_1为time，后面的变量值为升序排列。 （6）所有的值，包括坐标值，都必须使用SI单位制。fluent读入profile文件时不进行单位转换。 3. 问题的求解 3.3 Fluent前处理 一个实际的例子： ((sampleprofile transient