多孔介质

ANSYSFLUENT参考手册112.4.3可压缩流动的求解策略可压缩流动求解中速度、密度、压力和能量的高度耦合以及可能存在的激波导致求解过程不稳定。有助于改善可压缩流动计算过程稳定性的方法有：（仅适用于基于压力求解器）以接近于滞止条件的流动参数进行初始化（即，压力很小但不为零，压力和温度分别等于进口总压和总温）。在迭代过程的最初几十步不求解能量方程。设置能量方程的亚松驰因子等于1，压力的亚松驰因子0.4，动量的亚松驰因子0.3。求解过程稳定后再加入能量方程的求解，并将压力的亚松驰因子提高到0.7。设置合理的温度和压力限制值以避免求解过程发散。必要时，先以较低的进、出口边界压力比进行求解，然后再逐步升高压力比直到预定工况。对于低Mach数流动，也可以先求解不可压缩流动，然后以所得到的解作为可压缩流动的迭代初值。某些情况下，也可以先求解无粘性流动作为迭代初值。2.5无粘性流动在高Re数流动中，惯性力相对于粘性力而言起支配作用，可忽略粘性的影响。例如高速飞行器在空气动力学方案分析阶段可以采用无粘性流动计算初步确定外形，然后进行粘性计算，将流体粘性和湍流粘性对升力和阻力的影响计入。无粘性流动计算的另一个用途是给复杂的流动提供好的迭代初值。对于特别复杂的问题，有时这是唯一能使求解过程进行下去的方法。无粘性流动的计算求解Euler方程。其中，质量方程与粘性流动的相同：mvSt(2.34)其动量方程与粘性流动的相比，没有粘性应力项：vvvpgFt(2.35)能量方程与粘性流动相比，没有导热项effkT和粘性耗散项effv：jjhjEvEphJSt(2.36)式(2.34)~式(2.36)中符号的意义与粘性流动控制方程的相同（见2.1.1~2.1.3节）。2.6多孔介质模型多孔介质（PorousMedia）模型可用于模拟许多问题，包括流过填充床、滤纸、多孔板、布流器、管排等的流动。多孔介质模型在流体区上定义（见17.2.1节）。此外，一个被称为多孔阶跃面（porousjump）的多孔介质模型的一维简化可用于模拟已知速度−压降特性的薄膜。多孔阶跃面在界面区上定义。多孔阶跃面比多孔介质模型更健壮，收敛性更好，应ANSYSFLUENT参考手册12首选采用（见　节）。2.6.1基于表观速度的多孔介质动量方程对于单相介质和多相介质，多孔介质模型可以使用表观速度或物理速度形式的公式。基于表观速度的多孔介质模型根据多孔介质区中的体积流量率计算表观相速度或混合物速度。基于表观速度的多孔介质模型能够较好模拟多孔介质区内部的压力损失。但是，在多孔介质区与非多孔介质区的交界面处的表观速度与的速度是相同的，不能反映实际速度变化所引起的动量变化，对计算精度不利。多孔介质模型通过在动量方程中增加源项来模拟计算域中多孔性材料对流体的流动阻力。该源项由两部分组成，即Darcy粘性阻力项和惯性损失项：33111,,2iijjijjjjSDvCvvixyz(2.37)其中，D和C分别为粘性阻力和惯性损失系数矩阵。这个负的动量源项导致多孔介质单元中的压力降。同时，在全部变量的输运方程和连续性方程中，瞬态项变为t，其中γ为孔隙率。对于简单的均匀多孔介质，分别在系数矩阵D和C中对角线项代入1/α和C2，而其它项为零，则有21,,2iiiSvCvvixyz(2.38)其中，α为渗透率；C2为惯性阻力系数。也可以用速度大小的幂函数来模拟阻力：11100,,CCiiSCvCvvixyz(2.39)式中，C0和C1为经验系数，且C0的单位为SI制。采用幂函数时，压力降为各向同性的。2.6.2Darcy粘性阻力项多孔介质中流动为层流时，典型情况下压力降与速度成正比，即多孔介质模型简化为Darcy定律：pv(2.40)于是，在三个坐标方向上的压力降为31,,ijiijjpvnixyz(2.41)式中，1/αij为系数矩阵D的项；Δni为多孔介质在三个坐标方向上的厚度。ANSYSFLUENT参考手册132.6.3惯性损失项当速度比较高，或模拟多孔板和管排时，有时可忽略渗透项，只保留惯性损失项，则多孔介质方程简化为3211,,2ijijpCvvixyz(2.42)或写成三个坐标方向上的压力降：3211,,2iijijjpCnvvixyz(2.43)式中，C2,ij为系数矩阵C中的项；Δni为多孔介质在三个坐标方向上的厚度。2.6.4多孔介质中能量方程的处理对多孔介质修正了扩散项和瞬态项的能量方程为1hffssffeffjjefffjEEvEpkThJvSt(2.44)其中，Ef为流体总能；Es为多孔介质基体固体总能；γ为孔隙率；hfS为流体焓的源项；keff为多孔介质的有效导热系数，采用流体导热系数（包括湍流有效导热系数）kf与多孔介质中固体材料的导热系数ks的体积加权平均：1efffskkk(2.45)采用UDF可以定义各向异性的有效导热系数。孔隙率γ定义为多孔介质区中流体的体积分数，也就是介质中空的部分所占的比例。孔隙率影响传热计算、输运方程中的非稳态项、以及介质中的化学反应和体积力。如果希望模拟介质为全空（即没有固体介质）的情况，应给定孔隙率等于1。2.6.5多孔介质中湍流的处理在多孔介质中，当介质的渗透性很大且介质的几何尺度与湍流涡的尺度不发生相互作用时，可以认为固体基体对湍流的生成和耗散率没有影响。但其它情况下应降低多孔介质中湍流的影响。当采用湍流模型时（LES除外），可通过将多孔介质指定为层流区（LaminarZone）而使湍流粘性μt为零来抑制多孔介质区中湍流效应。此时，进口湍流量被输运穿过多孔介质区，而其对流体混合及动量的影响被忽略，同时介质中湍流生成被置为零。2.6.6粘性阻力系数和惯性阻力系数阻力系数一般是基于流体在多孔介质中的表观速度定义的。阻力源项的计算可以采用相对速度或绝对速度。选择RelativeVelocityResistanceFormulation（相ANSYSFLUENT参考手册14对速度阻力公式）选项可以更精确计算有动网格和运动参考坐标系时的源项。对于高度各向异性的多孔介质，当使用基于压力求解器时，选择AlternativeFormulation（非常规公式）选项可以使求解过程更稳定。采用非常规公式时，通过多孔介质的压力损失取决于速度矢量第i个方向分量的大小：12iiiiSCvv(2.46)计算粘性阻力系数和惯性阻力系数的方法如下：(1)已知压力降，计算基于表观速度的阻力系数使用多孔介质模型时，FLUENT假定单元中没有多孔介质的固体基体，即单元是100%开孔的（100%open），且所给定的阻力系数值是基于这一假设的。在已知流体流过实际设备中多孔介质的压力降Δp与速度的关系时，可计算C2。流体流过开孔率为open%的多孔板时，基于实际流动速度的压力损失系数KL定义为：2%12LopenpKv(2.47)式中，%openv为流过多孔板的实际流速。对于100%开孔时的压力损失系数值LK，有2100%12LpKv(2.48)式中，100%v为流过开孔率100%多孔板时的流速。而在相同流量下，速度与开孔率成反比，将KL折算为100%开孔时的压力损失系数值LK：22%100%100%%openLLLvKKKvopen(2.49)阻力系数C2为单位厚度多孔板的压力损失系数：2LKCn(2.50)式中，Δn为多孔板厚度。(2)使用Ergun公式计算通过层床的阻力系数在湍流时，层床用渗透率和惯性损失系数模拟。对于多种类型的层床，在较宽的Re数范围内阻力系数可以采用半经验的Ergun公式计算：222331501.7511pppvvLDD¥¥(2.51)当层床中为层流时，忽略式(2.51)中的第二项，可得Blake-Kozeny方程：2231501ppvLD¥(2.52)式中，μ为粘性系数；Dp为平均颗粒直径；L为床厚度；ε为孔隙率，其定义为孔隙体积与层床总体积之比。比较式(2.40)、式(2.42)和式(2.51)，可得各方向粘性阻力系数和惯性损失系数ANSYSFLUENT参考手册1522311150pD(2.53)2313.5pCD(2.54)(3)使用经验公式计算流过多孔板湍流的阻力系数流过锐边孔多孔板的压力损失系数可以采用VanWinkle等的公式计算（适用于孔呈等边三角形布置的情况）：221ffppmCAAA(2.55)式中，m为通过板的流量；Af为孔的总面积；Ap为板的总面积；C为适用于不同Re数范围和不同孔径厚度比D/t情况下的系数，t/D1.6且Re4000时（Re数的特征尺寸为孔径，特征速度为孔内的速度），C0.98。利用式(2.55)和pmvA，可得2221112fpAApvxCt(2.56)式中，v为表观速度而非孔内的流速。与式(2.42)比较可得在垂直于板方向的阻力系数C2：22211fpAACCt(2.57)(4)用实验数据计算流过纤维状材料层流的阻力系数在已知任意排列的纤维材料的无量纲渗透率B与纤维体积分数之间关系的情况下，粘性阻力系数1/α可由无量纲渗透率的定义B=α/a2（a为纤维直径）确定。(5)用压力降与速度关系实验数据计算阻力系数可以用通过多孔介质的压力降Δp与速度v关系的实验数据确定阻力系数。设实验数据用二次多项式拟合为：212pavav(2.58)式中，a1和a2为拟合系数。动量方程源项为单位长度的压力降，即ipSn(2.59)式中，Δn为多孔介质厚度。则比较式(2.38)和式(2.58)及式(2.59)，可得阻力系数2212aCn(2.60)和11an(2.61)该方法也可以用于多孔阶跃面。ANSYSFLUENT参考手册162.6.7基于物理速度的多孔介质模型FLUENT默认情况下，在多孔介质中使用按体积流量率计算的表观速度。表观速度（SuperficialVelocity）与物理速度（PhysicalVelocity，即真实速度）的关系为superficialphysialvv(2.62)式中，γ为介质的孔隙率。由于孔隙率小于1，流体流入多孔介质中，物理速度会提高，而表观速度不反映出来。为精确模拟多孔介质中的流动，应求解物理速度，而不是表观速度。(1)单相多孔介质模型单相流动情况下各向同性多孔介质中的通用标量输运控制方程为vSt(2.63)体积平均质量方程和动量方程为0vt(2.64)22fCvvvpBvvt(2.65)式(2.65)中最后一项代表多孔介质对流体的粘性阻力和惯性阻力。采用物理速度求解时，式(2.65)中的两个阻力系数仍以表观速度计算（见本节2.6.6），FLUENT将其转换为与物理速度公式相应的值。入口质量流量亦是以表观速度计算的。对于相同的入口质量流量和阻力系数，对于表观速度或物理速度均应得到相同的压力降。(2)多相多孔介质模型可以使用物理速度多孔介质公式模拟包含有多孔介质区的多相流。关于多相流理论见第5章。各向同性多孔介质中第q相的通用变量q的控制方程取如下形式：,qqqqqqqqqqvSt(2.66)其中，γ为孔隙率；ρq为第q相的物理密度；αq为第q相流体体积