Fluent-movingzones动网格.ppt

2、域在移动的同时还会有形状的改变(变形),我们可以用动网格(DM)技术解算方程域的位置和形状的改变都可以用时间函数来表示非定常问题,y,运动坐标系域随着坐标系转动,动网格域的形状是时间函数,运动坐标系和动网格对比,建模方法概览,单参考系(SRF)整体计算域都置于运动坐标系下。多参考系(MRF)选择域中的部分区域置于运动坐标系下。忽略相互作用稳态混和面(MPM)在旋转/静止域的分界面上，混合模型需要考虑相邻区域相互之间的影响。忽略流动中的不均匀圆周运动稳态滑移网格(SMM)特定区域的移动使用移动网格算法对滑移分界面进行流动变量的插值非定常问题把边界的相互作用完全考虑进去了，但是计算量比 S

3、RF,MRF,或者MPM都更大。动网格(DM)和滑移网格相似，但是域随着时间变化可以移动和变形。网格变形技术有弹簧压缩，动态铺层，局部重构,大纲,介绍和模型建立方法概览单参考系(SRF)模型多重域和多参考系(MRF)模型混合面模型(MPM)滑移网格模型(SMM)动网格(DM)模型概要附录,建立SRF模型介绍,SRF一般针对单流域，只能是绕着某个特定的坐标点，以恒定的速度进行旋转。为什么使用旋转坐标系?非定常的流域当放在旋转坐标系下来看时，可以认为是定常的情况。定常问题比较容易求解边界条件简单低的计算机资源占用易于进行后处理和分析壁面条件必须符合下列标准：随着流域动的壁面可以是任意形状。固定(

4、至于确定的坐标)的壁面必须是旋转表面。选择旋转周期性边界条件可以提高效率(减少网格数量),SRF模型中的静止壁面,stationary wall,rotor,baffle,Correct,Wrong!,Wall with baffles not a surfaceof revolution!,N-S 方程：旋转坐标系,方程可以在绝对坐标系或者旋转坐标系下求解。相对速度公式(RVF)通过把静止坐标系下的N-S方程转化为旋转坐标系下获得。相对速度和相对总内能为依赖变量。绝对速度公式(AVF)由相对速度方程得到绝对速度和绝对总内能为依赖变量动量方程中的旋转源项,公式比较,相对速度方程：x方向上动量方

5、程,绝对速度方程：x方向上动量方程,Coriolis acceleration,Centripetal acceleration,Coriolis+Centripetal accelerations,SRF 设置：求解器,求解器选项推荐：Pressure Based不可压缩，低速（亚音速）可压缩流Density Based 高速(跨音速，超音速)可压缩流速度方程推荐：当从静止域流入时，使用绝对速度方程(AVF)当为闭域（所有表面都运动）或者从旋转域流入时，使用相对速度方程(RVF)注意：RVF 只在分离求解器中有效梯度选项非结构网格选用节点梯度多面体网格选用最小二乘梯度,SRF 设置：流体边界

6、条件,使用流体边界条件设置旋转坐标系的坐标原点和方向向量。方向向量必须是单位向量，如果不是，Fluent将会自动转化成单位向量。SRF中的Motion Type 选择运动参考坐标系（Moving Reference Frame）输入旋转速度旋转方向使用右手定则负速度表示速度方向相反,SRF 设置：入口/出口边界条件,速度入口：压力入口：通过速度方程来定义总压压力出口：当轴向流问题在出口处有旋转时，可以选择径向平衡选项(radial equilibrium assumption option)指定的压力为中心压力SRF问题的其它边界条件无反射边界条件目标质量流出口（Target mass flo

7、w outlet）,incompressible,RVF,incompressible,AVF,壁面边界条件,在运动参考坐标系下，可以选择壁面是相对坐标还是绝对坐标。推荐在所有的运动参考坐标系下，壁面条件都设为选上“Moving Wall”选项旋转坐标原点和方向的设置和流域中的设置一样对静止的表面(在绝对坐标下)输入零旋转速度，绝对速度（Absolute）对运动的表面，输入零旋转速度，相对于邻近的单元域,SRF问题的求解方法,SRF 问题可能相对较难求解，因为流域的旋转将会带相当来大的流动梯度。对于旋转产生的影响，可能需要很好的网格才能求解大梯度。推荐策略:确定网格质量足够好开始时使用一阶离

8、散方法。缩小松弛因子/克朗特数来保证数值稳定性。当遇到难以开始（hard-to-start）问题时，使用FMG初始化。特别需要注意的是压缩机、抽水机之类的问题。考虑使用瞬态计算来解算问题。相对于静态的方法，可以提供更好的收敛性。在时间上使用一阶离散方法，而且大概每次迭代2-3步。一直运行到静态。,大纲,介绍和模型建立方法概览单参考系(SRF)模型多重域和多参考系(MRF)模型混合面模型(MPM)滑移网格模型(SMM)动网格(DM)模型概要附录,多重域建模,许多旋转机械问题都有固定部分，不能使用表面旋转来描述(SRF 无效)。这样的系统可以分成多重流动域（multiple fluid zone

9、s）有些域是旋转着的，其它是静止的。使用一个或多个分界面（interfaces）把多重域联系起来。多重域模型中分界面的处理方法：多重参考系模型(MRF)定常求解法对分界面进行简化不考虑不同坐标系间的转动相互作用。混合面模型(MPM)定常求解法不同坐标系间的相互作用都近似成在流域分界面(mixing planes)上的周期性平均。滑移网格模型(SMM)非定常求解法精确模拟运动和静止域间的相对运动，但是需要占用很高的计算机资源。,单体和多体模型间的比较,interface,Single Component(blower wheel blade passage),Multiple Component

10、(blower wheel+casing),MRF模型的介绍,计算域被分为定子和转子。分界面把域分隔开。每个流域分别解算各自的流动方程每个域中都假定流动是静止的。在转子中，使用了SRF 方程。在定子和转子之间的分界面上，速度向量和速度梯度适当变化，便可以计算得出质量，动量，能量和其它标量。MRF 不同域之间的相对运动定子和转子之间，不考虑动力学相互作用。因为这个特性，MRF经常称作“冷冻转子”（frozen rotor）方法。,MRF模型的分界面,正交分界面把相邻的流域分开的内部网格面。两边的网格必须一样。非正交分界面网格单元是物理上分开的。分界面由两个重叠的面组成。(type=interfa

11、ce)用户创建分界面的方式是：DefineGrid Interfaces分界面可以是周期性的,Conformal interface,Non-conformal interface,分界面形状,包含旋转流域上分界面的壁面，将和流域一样移动，并且可以是任意形状。旋转表面允许是静止的。两个域之间的分界面必须是旋转表面，而且旋转轴和旋转域是一样的。,Correct,Wrong!,Interface is not a surface or revolution,Consider a mixing impeller inside a rectangular vessel:-Problem can be

12、described with two reference frames or zones.,MRF 设置,为静止和旋转的流域生成合适的的网格。可以为单元域选择正交和非正交类型。为每个旋转流域(流体边界条件),在 Motion Type 选项中选择 Moving Reference Frame，输入旋转速度。SRF中除了多重域，其它都一样静止域选上静止选项（Stationary）设置其它边界条件，求解器设置和 SRF 模型设置相同。,MRF问题的计算方法,和SRF问题相同，因为一个或者多个流域的流动梯度比较大，MRF问题也比较难解。相邻单元域间的作用，需要用比默认值更小的欠松弛因子求解。处理策

13、略：确认网格质量足够好，特别是相邻域间分界面周围的网格。使用FMG 初始化。可以使用非定常时间步长的方法，达到和静止状态下稳定解相同的效果。如果MRF的结果不稳定或者不符合实际，你就可能要用非定常滑移网格模型求解。,大纲,介绍和模型建立方法概览单参考系(SRF)模型多重域和多参考系(MRF)模型混合面模型(MPM)滑移网格模型(SMM)动网格(DM)模型概要附录,混合平面(MPM)模型介绍,MPM模型起初是用来模拟定子/转子，叶轮/叶片和轴心/离心涡轮机组。同样也能用于很多相同种类的问题。在每个域中都包含有恒稳态的SRF解，并且每个域的边界条件都互相关联。混合平面是指每个域间“相互关联”的域

14、。边界条件之间通过周期性平均侧面流动变量来传递的，这些变量每个迭代步都会进行更新。,radial machines,混合平面构成,axial machines,混合平面的分界面由上游域出口和邻近下游域入口组成。入口/出口边界条件必须是下面这几种边界条件中的一种：压力出口/压力入口Pressure-Outlet/Pressure-Inlet压力出口/速度入口Pressure-Outlet/Velocity-Inlet 压力出口/流量入口Pressure-Outlet/Mass-Flow-InletMPM可以同时实现轴向和径向涡轮组叶片流场的模拟。,MPM 设置,设置流域速度。为入口和出口边界条件

15、设置合适的边界条件类型。为混合平面对选择上游和下游域。为轮廓插值设置点的数量。混合平面几何（Mixing Plane Geometry）决定了周期性平均的方法轴向流动机械选择径向（Radial）。径向流动机械选择轴向（Axial）。混合平面控制松弛系数（Under-Relaxation）01,混合模型的计算方法,因为混合平面模型包括修改混合平面分界面的边界条件，流动条件变化过快很可能导致混合平面求解收敛上的困难。试着降低混合平面的欠松弛因子到0.10.5之间，可能对解的稳定性有所帮助。处理策略：确定网格质量足够好(最大网格偏斜 0.9 0.95)。为难以开始（hard-to-start）的问题

16、使用FMG初始化。FMG初始化适合混合平面模型。降低松弛因子或者克朗特数。使用固定的边界条件初始化运行，再使用混合平面法。,大纲,介绍和模型建立方法概览单参考系(SRF)模型多重域和多参考系(MRF)模型混合面模型(MPM)滑移网格模型(SMM)动网格(DM)模型概要附录,滑移网格(SMM)模型介绍,涡轮机组中定子和转子之间的相对运动可能会有非定常相互作用。这些相互作用大体上分为下列几种：位势相互作用(压力波相互作用)尾迹相互作用震动相互作用MRF和MPM都忽略了全部的非定常相互作用，而且限定这些流动中的效果十分弱。如果非定常相互作用不可以忽略，我们可以使用滑移网格模型，把静止和旋转组成中的

17、相对运动考虑进去。,滑移网格模型的工作原理,和 MRF模型相同的是，非正交分界面把域分成移动和静止域。和 MRF模型不同的是，每个移动域的网格会随着一个时间函数更新，因此为非定常。另外和MRF模型不同的是，控制方程有新的移动网格形式，而且在绝对坐标系下求解。(附录中有详细信息)没有使用运动参考坐标系方程(例如：没有科氏力，向心加速度)。方程是普通移动/变形网格的特殊公式。假设刚体网格运动和滑移，非正交分界面。,滑移分界面,MRF问题中的滑移网格必须符合下列规则：界面的平移不能垂直于自身。位于旋转域和它相邻的静止/旋转域必须是个旋转面，这个旋转面的转动轴和旋转子域重合。许多失败的滑移网格模型，很

18、多都是因为网格移动时，分界面变得不相连！滑移分界面可以是部分重叠。可以是：周期性的壁面,Elliptic interface is not a surface of revolution.,滑移网格设置,选择非定常求解器（unsteady solver）。定义滑移域的分界面条件类型。对于移动域，选择移动网格（Moving Mesh）作为移动类型（Motion Type）。为每对分界面域，创建非正交分界面。滑移/旋转是周期性的，选择周期性（Periodic）选项。耦合传热，选择Coupled。其它的边界条件和SRF,MRF模型设置相同。,SMM问题计算方法,设置适当的时间步长和每次时间步最大迭代

19、步，得到每个时间步较好的收敛情况。时间步长不能比移动单元时间长：在流动达到周期性之前不断推进求解通常需要几次网格的旋转。好的初始条件可以减少到达时间周期性需要的时间。,Ds=average cell sizewR=characteristic velocity of the moving zone(R=mean radius of interface surface),大纲,介绍和模型建立方法概览单参考系(SRF)模型多重域和多参考系(MRF)模型混合面模型(MPM)滑移网格模型(SMM)动网格(DM)模型概要附录,什么是动网格(DM)模型?,FLUENT求解中，一种可以适应移动边界条件/物

20、体，并且可以相应的调节网格与之适应的方法。,案例：活塞运动机翼抖动阀门开关动脉扩张和收缩,Volumetric fuel pump,动网格模型：特点,根据用户自定义的边界/物体移动，单元类型和网格形式，FLUENT会自动计算出内部节点位置。弹簧压缩局部重构动态铺层2.5 D，6DOF用户自定义网格移动边界/物体运动可以根据下列运动类型设定：活塞运动（In-cylinder）(RPM,行程长度,曲柄角,)通过预置文件（profiles）或者UDF指定的运动流体动力学耦合运动（6 DOF）。不同的网格运动形式可以运用在不同的域上。相邻域间的连接可以是非正交的。,弹簧压缩,网格的移动就像相连的弹簧，

21、或者是海绵。连通性并没有改变;对结构和非结构网格都有效;,局部网格重构,当超出用户自定义的偏斜和尺度限制范围时，局部网格节点和单元就会增加或减少。当单元增加或减少时，连通性便发生了改变。只有非结构网格有效。这个过程同时也有光滑作用,动态铺层,单元增加或减少随着域的增长或收缩改变。当单元增加或减少时，连通性也相应改变。只有结构网格才有效。,混合运用的方法,需要正确的分解来初始生成网格动态铺层:阀程区域；下汽缸区域。网格重构:上汽缸区域。域间非正交分界面,动网格设置,选择非定常求解器。选择动网格模型 DefineDynamic Mesh。激活需要的网格方法，设置合适的参数。在动网格用户操作界面中定

22、义边界运动。可能需要UDF。其它模型，边界条件和求解器设置和SRF,MRF模型相同。网格运动可以提前查看SolveMesh Motion utility.通过提前查看，如果网格质量有问题，Fluent可以提前发现。,其它动网格选项,2.5 D 网格运动挤压几何（extruded geometries）的特殊网格重构方法。可以允许三角形网格像楔形一样在体积内挤压。优点可以比四面体网格更好的贴合小裂缝。用户自定义网格运动通过UDF控制网格的运动。网格拓扑结构不能发生改变。优点可以提供其它三种方法不能提供的运动模式。6自由度模型根据计算结果来调整网格运动,大纲,介绍和模型建立方法概览单参考系(SR

23、F)模型多重域和多参考系(MRF)模型混合面模型(MPM)滑移网格模型(SMM)动网格(DM)模型概要附录,总结,五种可以模拟流体流动部分的方法。单（旋转）参考系模型多重参考系模型混合平面模型滑移网格模型动网格模型前三种方法是完全静态定常方法，而滑移网格和动网格则是非定常方法。大部分物理模型适合用运动参考系或者动网格（例如：多相，燃烧，传热等）,附录,运动参考系下的N-S方程相对速度方程绝对速度方程动网格问题中的N-S方程,N-S 方程:旋转参考坐标系,Fluent中使用的两种不同的方程相对速度方程(RVF)包括把固定参考系下的N-S方程转换到旋转坐标系下。使用相对速度作为动量方程中的依赖

24、变量。使用相对总内能作为能量方程中的依赖变量。绝对速度方程(AVF)从相对速度方程得来。使用绝对速度作为动量方程中的依赖变量。使用绝对总内能作为能量方程中的依赖变量。注意:N-S方程中RVF and AVF 是等价的！边界条件等价的情况下，两种方程的解应该是相同的。,参考系,x,y,z,z,y,x,stationaryframe,rotatingframe,axis ofrotation,CFD domain,Note:R is perpendicularto axis of rotation,假设没有交换()特定坐标下的定常旋转(w=constant)通过原点旋转坐标定义坐标忽略物体的受力，

25、比如重力和其它效果(像方程中列出的)忽略能量源(像方程中列出的)定义绝对速度()流体在静止（绝对）参考坐标系的速度相对速度()-流体在旋转参考坐标系下的速度3-D 可压缩，层流形式的方程在接下来的幻灯片中将会介绍到(其它的形式相似),假设和定义,相对速度方程,(Continuity),(Momentum),(Energy),相对速度方程(2),(Relative total internal energy),(Fouriers Law),(Viscous stress),旋转坐标系下的相对加速度,Coriolisacceleration,centrifugalacceleration,Roth

26、alpy,考虑在旋转参考系下，流动是定常，绝热，无粘流。能量方程可以简化成：数量便是 rothalpy。从上面可以看出，rothalpy在旋转坐标系下，定常，无粘流将是守恒的。如果粘性流的壁面和参考系一起旋转，Rothalpy也可以是守恒的。(因为这时旋转壁面的相对速度为零)。,绝对速度方程,(Continuity),(Momentum),(Energy),绝对速度方程(2),(Relative total internal energy),(Fouriers Law),(Viscous stress),旋转坐标系下的加速度,Acceleration reduces to single te

27、rm involvingrotational speed and absolute velocity,推荐的速度方程,当入流从静止域进入时，使用绝对速度参考系（AVF）在这种情况下，通常知道的是绝对总压和绝对速度例子：导管风扇系统中的流动，入流处是个静止管封闭域（所有的表面是动态的）或者入流冲旋转域进入时，使用相对速度参考系（RVF）在这种情况下，通常知道的是相对总压力和相对速度例子：圆形洞穴中的旋转流正如先前知道的，RVF和AVF 等价的，而且两种都能在许多问题中很好的运用。如果绝对速度梯度量和相对速度梯度量非常不同，相同的网格结果可能会有差异。使用合适精致的网格时，不同求解方式中的结果应该

28、相似。,运动参考系中的标量方程,运动参考系中普通的标量输运方程形式可以写成下面的形式：,方程中的标量是沿着相对流线（relative streamlines）的标量。梯度，源项都是关于运动控制体积来定义的。例子：湍流模型方程，物质守恒方程等。,N-S 方程：动网格形式,滑移网格（aka动网格）和动网格方程中假设计算域在相对静止参考系下的运动。计算域没有参考系。域中任何运动的点都是关于一个给定的时间变化率的位置向量().可以认为是网格速度刚体的旋转速度是常数方程是积分型的。,动网格说明,x,y,z,stationaryframe,axis ofrotation,Moving CFD domain

29、,N-S 方程：动网格(1),(Continuity),(Momentum),(Energy),N-S 方程：动网格(2),前面所说的方程中，W 和分别是移动控制体积中的体积和边界面。滑移网格模型中，即使W 在网格中没有定义，W 始终是常数。移动/变形模型中，W=W(t)，而且符合几何守恒定律（geometric conservation law）(GCL):时间导数(d/dt)表示的是时间顺着运动控制体积的微分。所有空间梯度都是在静止参考系下计算的。,2.5 D 模型,2.5D网格本质上是2D上的三角形网格沿着特定动域坐标轴拓展而成的，你可能更关注怎样建立这样的模型。刚体运动是随着运动边界

30、域来运动。三角形拓展面是通过网格重构和光滑来进行调整的。拓展三角形网格的对边，然后同时分配变形域，只是使用光滑技术。,用户自定义网格运动,网格运动是通过用户子定义函数（UDF）来定义的。如果使用UDF移动网格，将不能产生连通性的改变。实际运用包括：叶轮泵离心泵轴承旋转压缩机,6自由度耦合运动,物体得到的运动，同时考虑了空气动力和瞬间其它力一起的作用，例如重力，推力或者爆破力（例如：为了避免碰撞，初始时推动物体离开飞机或者火箭的力）。这些情况下，运动和流域是耦合的，我们于是把这种运动叫做耦合运动。Fluent提供可以计算物体运动轨迹的UDF，计算时考虑空气动力/动量，重力和爆破力。这个通常叫做6自由度（6-DOF）求解器，一般和它相关的UDF叫做 6-DOF UDF。The 6-DOF UDF 是完全并行的。,6 DOF 耦合运动(contd),储箱在马赫数1.2时，从delta翼形上掉落(NACA 64A010)爆破力在一个很短的时间起作用。全非结构网格使用尺度函数来光滑和网格重构Fluent的计算结果和风动试验数据相稳和！,