1、例如高速飞行器在空气动力学方案分析阶段可以采用无粘性流动计算初步确定外形,然后进行粘性计算,将流体粘性和湍流粘性对升力和阻力的影响计入。无粘性流动计算的另一个用途是给复杂的流动提供好的迭代初值。对于特别复杂的问题有时这是唯一能使求解过程进行下去的方法。无粘性流动的计算求解Euler方程。其中质量方程与粘性流动的相同: 其动量方程与粘性流动的相比,没有粘性应力项粘性耗散项能量方程与粘性流动相比,式（2.34）式（2.36）中符号的意义与粘性流动控制方程的相同见（2.1.12.1.3节）。2.6多孔介质模型多孔介质（Porousmedia）模型可用于模拟许多问题,包括流过填充床、滤纸、多孔板、布流
2、器、管排等的流动。多孔介质模型在流体区上定义（见17.2.1节）。此外,一个被称为多孔阶跃面（porousjump）的多孔介质模型的一维简化可用于模拟已知速度?压降特性的薄膜。多孔阶跃面在界面区上定义。多孔阶跃面比多孔介质模型更健壮,收敛性更好.应anSYSFLUEnT参考手册12首选采用。2.6.1基于表观速度的多孔介质动量方程对于单相介质和多相介质,多孔介质模型可以使用表观速度或物理速度形式的公式。基于表观速度的多孔介质模型根据多孔介质区中的体积流量率计算表观相速度或混合物速度。基于表观速度的多孔介质模型能够较好模拟多孔介质区内部的压力损失。但是在多孔介质区与非多孔介质区的交界面处的表观速
3、度与的速度是相同的,不能反映实际速度变化所引起的动量变化,对计算精度不利。多孔介质模型通过在动量方程中增加源项来模拟计算域中多孔性材料对流体的流动阻力。该源项由两部分组成即darcy粘性阻力项和惯性损失项 其中,d和c分别为粘性阻力和惯性损失系数矩阵。这个负的动量源项导致多孔介质单元中的压力降。同时,在全部变量的输运方程和连续性方程中,瞬态项变为,其中为孔隙率。对于简单的均匀多孔介质,分别在系数矩阵d和c中对角线项代入1/和c2,而其它项为零,则有:其中为渗透率c2为惯性阻力系数。也可以用速度大小的幂函数来模拟阻力: 式中c0和c1为经验系数,且c0的单位为Si制。采用幂函数时压力降为各向同性
4、的。2.6.2darcy粘性阻力项多孔介质中流动为层流时,典型情况下压力降与速度成正比,即多孔介质模型简化为darcy定律:于是,在三个坐标方向上的压力降为（）式中1/ij为系数矩阵d的项ni为多孔介质在三个坐标方向上的厚度。2.6.3惯性损失项当速度比较高,或模拟多孔板和管排时,有时可忽略渗透项.只保留惯性损失项,则多孔介质方程简化为或写成三个坐标方向上的压力降:式中,c2,ij为系数矩阵c中的项ni为多孔介质在三个坐标方向上的厚度。2.6.4多孔介质中能量方程的处理对多孔介质修正了扩散项和瞬态项的能量方程为其中Ef为流体总能；Es为多孔介质基体固体总能r为孔隙率；为流体焓的源项，keff为
5、多孔介质的有效导热系数，采用流体导热系数（包括湍流有效导热系数）kf与多孔介质中固体材料的导热系数ks的体积加权平均：采用UdF可以定义各向异性的有效导热系数。孔隙率定义为多孔介质区中流体的体积分数，也就是介质中空的部分所占的比例。孔隙率影响传热计算、输运方程中的非稳态项、以及介质中的化学反应和体积力。如果希望模拟介质为全空（即没有固体介质）的情况，应给定孔隙率等于1。2.6.5多孔介质中湍流的处理在多孔介质中，当介质的渗透性很大且介质的几何尺度与湍流涡的尺度不发生相互作用时，可以认为固体基体对湍流的生成和耗散率没有影响。但其它情况下应降低多孔介质中湍流的影响。当采用湍流模型时（LES除外），
6、可通过将多孔介质指定为层流区（Laminarzone）。而使湍流粘性t为零来抑制多孔介质区中湍流效应。此时，进口湍流量被输运穿过多孔介质区，而其对流体混合及动量的影响被忽略，同时介质中湍流生成被置为零。2.6.6粘性阻力系数和惯性阻力系数阻力系数一般是基于流体在多孔介质中的表观速度定义的。阻力源项的计算可以采用相对速度或绝对速度。选择RelativeVelocityResistanceFormulation （相对速度阻力公式），选项可以更精确计算有动网格和运动参考坐标系时的源项。对于高度各向异性的多孔介质，当使用基于压力求解器时，选择alternativeFormulation非常规公式，选
7、项可以使求解过程更稳定。采用非常规公式时，通过多孔介质的压力损失取决于速度矢量第i个方向分量的大小 计算粘性阻力系数和惯性阻力系数的方法如下：（1）已知压力降，计算基于表观速度的阻力系数使用多孔介质模型时，FLUEnT假定单元中没有多孔介质的固体基体，即单元是100%开孔的（100%open），且所给定的阻力系数值是基于这一假设的。在已知流体流过实际设备中多孔介质的压力降p与速度的关系时，可计算c2。流体流过开孔率为open%的多孔板时，基于实际流动速度的压力损失系数KL定义为式中 V%open为流过多孔板的实际流速。对于100%开孔时的压力损失系数值，有式中V100%为流过开孔率100%多孔
8、板时的流速。而在相同流量下,速度与开孔率成反比,将KL折算为100%开孔时的压力损失系数值阻力系数c2为单位厚度多孔板的压力损失系数 式中n为多孔板厚度。（2）使用Ergun公式计算通过层床的阻力系数在湍流时，层床用渗透率和惯性损失系数模拟。对于多种类型的层床，在较宽的Re数范围内阻力系数可以采用半经验的Ergun公式计算：当层床中为层流时，忽略式（2.51）中的第二项，可得Blake-Kozeny方程： 式中 为粘性系数，dp为平均颗粒直径，L为床厚度， 为孔隙率，其定义为孔隙体积与层床总体积之比。比较式（2.40）、式（2.42）和式（2.51），可得各方向粘性阻力系数和惯性损失系数（3）
9、使用经验公式计算流过多孔板湍流的阻力系数流过锐边孔多孔板的压力损失系数可以采用Vanwinkle等的公式计算（适用于孔呈等边三角形布置的情况）:式中,为通过板的流量；af为孔的总面积；ap为板的总面积；c为适用于不同Re数范围和不同孔径厚度比d/t情况下的系数，t/d1.6且Re4000时（Re数的特征尺寸为孔径，特征速度为孔内的速度）c0.98。利用式（2.55）和式中v为表观速度而非孔内的流速。与式（2.42）比较可得在垂直于板方向的阻力系数c2：（4）用实验数据计算流过纤维状材料层流的阻力系数在已知任意排列的纤维材料的无量纲渗透率B与纤维体积分数之间关系的情况下，粘性阻力系数1/可由无量
10、纲渗透率的定义a为纤维直径确定。（5）用压力降与速度关系实验数据计算阻力系数可以用通过多孔介质的压力降p与速度v关系的实验数据确定阻力系数。设实验数据用二次多项式拟合为式中a1和a2为拟合系数。下页余下全文多孔介质-技术总结动量方程源项为单位长度的压力降， 即式中n为多孔介质厚度。则比较式（2.38）和式（2.58）及式（2.59），可得阻力系数和该方法也可以用于多孔阶跃面。2.6.7基于物理速度的多孔介质模型FLUEnT默认情况下，在多孔介质中使用按体积流量率计算的表观速度。表观速度（SuperficialVelocity）与物理速度（PhysicalVelocity）即真实速度的关系为式中
11、 为介质的孔隙率。由于孔隙率小于1流体流入多孔介质中物理速度会提高，而表观速度不反映出来。为精确模拟多孔介质中的流动，应求解物理速度，而不是表观速度。（1）单相多孔介质模型，单相流动情况下各向同性多孔介质中的通用标量输运控制方程为体积平均质量方程和动量方程为式（2.65）中最后一项代表多孔介质对流体的粘性阻力和惯性阻力。采用物理速度求解时，式（2.65）中的两个阻力系数仍以表观速度计算（见本节2.6.6）FLUEnT将其转换为与物理速度公式相应的值。入口质量流量亦是以表观速度计算的。对于相同的入口质量流量和阻力系数对于表观速度或物理速度均应得到相同的压力降。（2）多相多孔介质模型可以使用物理速
12、度多孔介质公式模拟包含有多孔介质区的多相流。关于多相流理论见第5章。各向同性多孔介质中第q相的通用变量的控制方程取如下形式：其中 为孔隙率；pq为第q相的物理密度；q为第q相流体体积分数；为第q相的速度；为第q相通用扩散系数；为源项。质量方程和动量方程为通用变量控制方程（2.66）适用于Euler多相流模型的所有输运方程。质量方程和动量方程为式中最后一项为多孔介质中的动量阻力源项。该项由两部分组成 粘性损失项和惯性损失项。K为渗透率，c2为惯性阻力系数，二者均为（1?）的函数。能量方程为式中 Qsp为多孔介质中固体表面与第q相的传热量。默认情况下，FLUEnT假定多孔介质的固体与多相流体之间处
13、于热平衡，则且（）但也可以用求解用户定义标量（UdS）的方式单独求解多孔介质固体的导热方程： 这时 如仅考虑对流换热，有式中hq,eff为有效传热系数，Ts为多孔介质固体表面温度。2.6.8多孔介质模型的限制和求解策略多孔介质模型的假定和限制条件多孔介质对湍流影响的模拟是近似的。当在运动坐标系中应用多孔介质模型，多孔介质采用相对速度形式的阻力公式时，动量方程可以采用相对速度形式或绝对速度形式。当多孔介质区中在流动方向上压力降较大时（例如渗透率较小或惯性系数c2较大），收敛速度较慢。解决收敛性问题的最好方法是估算多孔介质压力降的合适的迭代初值，并以分块（Patch）的方式初始化，使多孔介质区上、下游的初始压力差满足该压力降值。另一个方法是暂时停用多孔介质模型求解 获得没有多孔介质的流场的初步解，然后再启用多孔介质模型继续求解（此方法对于高流阻多孔介质不适用）。高度各向异性的多孔介质模型可能有收敛困难的问题 可将各方向多孔介质系数（1/ij和c2,ij）的相差倍数限制在23个数量级以内来解决这一问题。如果某一方向的介质阻力为无限大，只需将其置为主流方向阻力的1000倍。上页
