动网格技术

题记：在学习使用Fluent的时候，有不少朋友需要使用动网格模型(Dynamic Mesh Model)，因此，本版推出这个专题，进行大讨论，使大家在使用动网格时尽量少走弯路，更快更好地掌握；也欢迎使用过的版友积极参与讨论指导，谢谢！

该专题主要包括以下的主要内容：

##1.动网格的相关知识介绍；

##2.以NACA0012翼型俯仰振荡实例进行讲解动网格的应用过程；

##3. 与动网格应用有关的参考文献；

##4. 使用动网格进行计算的一些例子。

##1.动网格的相关知识介绍

有关动网格基础方面的东西，请具体参考FLUENT User’s Guide或FLUENT全攻略的相关章节，这里只给出一些提要性的知识要点。

1、简介

动网格模型可以用来模拟流场形状由于边界运动而随时间改变的问题。边界的运动形式可以是预先定义的运动，即可以在计算前指定其速度或角速度；也可以是预先未做定义的运动，即边界的运动要由前一步的计算结果决定。网格的更新过程由FLUENT 根据每个迭代步中边界的变化情况自动完成。在使用动网格模型时，必须首先定义初始网格、边界运动的方式并指定参予运动的区域。可以用边界型函数或者UDF 定义边界的运动方式。FLUENT 要求将运动的描述定义在网格面或网格区域上。如果流场中包含运动与不运动两种区域，则需要将它们组合在初始网格中以对它们进行识别。那些由于周围区域运动而发生变形的区域必须被组合到各自的初始网格区域中。不同区域之间的网格不必是正则的，可以在模型设置中用FLUENT软件提供的非正则或者滑动界面功能将各区域连接起来。

注：一般来讲，在Fluent中使用动网格，基本上都要使用到UDF，所以你最好具备一定的C 语言编程基础。

2、动网格更新方法

动网格计算中网格的动态变化过程可以用三种模型进行计算，即弹簧近似光滑模型(spring-based smoothing)、动态分层模型（dynamic layering）和局部重划模型（local remeshing）。

弹簧近似光滑模型

在弹簧近似光滑模型中，网格的边被理想化为节点间相互连接的弹簧。移动前的网格间距相当于边界移动前由弹簧组成的系统处于平衡状态。在网格边界节点发生位移后，会产生与位移成比例的力，力量的大小根据胡克定律计算。边界节点位移形成的力虽然破坏了弹簧系统原有的平衡，但是在外力作用下，弹簧系统经过调整将达到新的平衡，也就是说由弹簧连接在一起的节点，将在新的位置上重新获得力的平衡。从网格划分的角度说，从边界节点的位移出发，采用虎克定律，经过迭代计算，最终可以得到使各节点上的合力等于零的、新的网格节点位置，这就是弹簧光顺法的核心思想。

原则上弹簧光顺模型可以用于任何一种网格体系，但是在非四面体网格区域（二维非三角形），最好在满足下列条件时使用弹簧光顺方法：

（1）移动为单方向。

（2）移动方向垂直于边界。

如果两个条件不满足，可能使网格畸变率增大。另外，在系统缺省设置中，只有四面体网格（三维）和三角形网格（二维）可以使用弹簧光顺法，如果想在其他网格类型中激活该模型，需要在dynamic-mesh-menu 下使用文字命令spring-on-all-shapes?，然后激活该选项即可。

动态层模型

对于棱柱型网格区域（六面体和或者楔形），可以应用动态层模型。动态层模型的中心思想是根据紧邻运动边界网格层高度的变化，添加或者减少动态层，即在边界发生运动时，如果紧邻边界的网格层高度增大到一定程度，就将其划分为两个网格层；如果网格层高度降低到一定程度，就将紧邻边界的两个网格层合并为一个层：

如果网格层j扩大，单元高度的变化有一临界值：

H_min>(1+alpha_s)*h_0

式中h_min为单元的最小高度，h_0为理想单元高度，alpha_s为层的分割因子。在满足上述条件的情况下，就可以对网格单元进行分割，分割网格层可以用常值高度法或常值比例法。在使用常值高度法时，单元分割的结果是产生相同高度的网格。在采用常值比例法时，网格单元分割的结果是产生是比例为alpha_s的网格。

若对第j层进行压缩，压缩极限为：

H_min

式中alpha_c为合并因子。在紧邻动边界的网格层高度满足这个条件时，则将这一层网格与外面一层网格相合并。

动网格模型的应用有如下限制：

（1）与运动边界相邻的网格必须为楔形或者六面体（二维四边形）网格。

（2）在滑动网格交界面以外的区域，网格必须被单面网格区域包围。

（3）如果网格周围区域中有双侧壁面区域，则必须首先将壁面和阴影区分割开，再用

滑动交界面将二者耦合起来。

（4）如果动态网格附近包含周期性区域，则只能用FLUENT 的串行版求解，但是如果周期性区域被设置为周期性非正则交界面，则可以用FLUENT 的并行版求解。

如果移动边界为内部边界，则边界两侧的网格都将作为动态层参与计算。如果在壁面上只有一部分是运动边界，其他部分保持静止，则只需在运动边界上应用动网格技术，但是动网格区与静止网格区之间应该用滑动网格交界面进行连接。

局部重划模型

在使用非结构网格的区域上一般采用弹簧光顺模型进行动网格划分，但是如果运动边界的位移远远大于网格尺寸，则采用弹簧光顺模型可能导致网格质量下降，甚至出现体积为负值的网格，或因网格畸变过大导致计算不收敛。为了解决这个问题，FLUENT 在计算过程中将畸变率过大，或尺寸变化过于剧烈的网格集中在一起进行局部网格的重新划分，如果重新划分后的网格可以满足畸变率要求和尺寸要求，则用新的网格代替原来的网格，如果新的网格仍然无法满足要求，则放弃重新划分的结果。

在重新划分局部网格之前，首先要将需要重新划分的网格识别出来。FLUENT 中识别不合乎要求网格的判据有二个，一个是网格畸变率，一个是网格尺寸，其中网格尺寸又分最大尺寸和最小尺寸。在计算过程中，如果一个网格的尺寸大于最大尺寸，或者小于最小尺寸，或者网格畸变率大于系统畸变率标准，则这个网格就被标志为需要重新划分的网格。在遍历所有动网格之后，再开始重新划分的过程。局部重划模型不仅可以调整体网格，也可以调整动边界上的表面网格。

需要注意的是，局部重划模型仅能用于四面体网格和三角形网格。在定义了动边界面以后，如果在动边界面附近同时定义了局部重划模型，则动边界上的表面网格必须满足下列条件：（1）需要进行局部调整的表面网格是三角形（三维）或直线（二维）。

（2）将被重新划分的面网格单元必须紧邻动网格节点。

（3）表面网格单元必须处于同一个面上并构成一个循环。

（4）被调整单元不能是对称面（线）或正则周期性边界的一部分。

动网格的实现在FLUENT 中是由系统自动完成的。如果在计算中设置了动边界，则FLUENT

会根据动边界附近的网格类型，自动选择动网格计算模型。如果动边界附近采用的是四面体网格（三维）或三角形网格（二维），则FLUENT 会自动选择弹簧光顺模型和局部重划模型对网格进行调整。如果是棱柱型网格，则会自动选择动态层模型进行网格调整。在静止网格区域则不进行网格调整。

动网格问题中对于固体运动的描述，是以固体相对于重心的线速度和角速度为基本参数加以定义的。既可以用型函数定义固体的线速度和角速度，也可以用UDF 来定义这两个参数。同时需要定义的是固体在初始时刻的位置。

注：这一小节主要讲述了动网格的更新方法，最好能掌握，尤其是各种方法的适用范围，通常来讲，在一个case中，我们使用的更新方法都是根据网格类型以及和要实现的运动来选择的，很多时候都是几种更新方法搭配起来使用的。

总结一下：

使用弹簧近似光滑法网格拓扑始终不变，无需插值，保证了计算精度。但弹簧近似光滑法不适用于大变形情况，当计算区域变形较大时，变形后的网格会产生较大的倾斜变形，从而使网格质量变差，严重影响计算精度。动态分层法在生成网格方面具有快速的优势，同时它的应用也受到了一些限制。它要求运动边界附近的网格为六面体或楔形，这对于复杂外形的流场区域是不适合的。使用局部网格重划法要求网格为三角形（二维）或四面体（三维），这对于适应复杂外形是有好处的，局部网格重划法只会对运动边界附近区域的网格起作用。

3、动网格问题的建立

设置动网格问题的步骤如下：

（1）在Solver（求解器）面板中选择非定常流（unsteady）计算。

（2）设定边界条件，即设定壁面运动速度。

（3）激活动网格模型，并设定相应参数，菜单操作如下：

Define -> Dynamic Mesh -> Parameters...

（4）指定移动网格区域的运动参数，菜单操作如下：

Define -> Dynamic Mesh -> Zones...

（5）保存算例文件和数据文件。

（6）预览动网格设置，菜单操作为：

Solve -> Mesh Motion...

（7）在计算活塞问题时，设定活塞计算中的事件：

Define -> Dynamic Mesh -> Events...

并可以通过显示阀与活塞的运动，检查上述设置是否正确：

Display -> IC Zone Motion...

（8）应用自动保存功能保存计算结果。

File -> Write -> Autosave...

在动网格计算中，因为每个计算步中网格信息都会改变，而网格信息是储存在算例文件中的，所以必须同时保存算例文件和数据文件。

（9）如果想建立网格运动的动画过程，可以在Solution Animation（计算结果动画）面板中进行相关设置。

注：在这一步中，需要提醒一下，使用动网格进行正式计算之前，最好养成预览动网格更新的习惯；就是在正式计算前，浏览一下动网格的更新情况，这样可以避免在计算过程中出现动网格更新本身的问题。在预览更新时，很多人都说会出现负体积的警告，更新不成功，出现这样的问题时，最好先把时间步长改的更小点儿试试，一般来讲，排除UDF本身的原因，出现更新出错的原因都与时间步长有关，这需要结合所使用的更新方法多琢磨。

4、设定动网格参数

为了使用动网格模型，需要在dynamic mesh（动网格）面板中激活Dynamic Mesh（动网格）选项。如果计算的是活塞运动，则同时激活In-Cylinder（活塞）选项。然后选择动网格模型，并设置相关参数。

1）选择网格更新模型

在Mesh Methods（网格划分方法）下面选择Smothing（弹簧光顺模型），Layering（动态层模型）和（或）Remshing（局部重划模型）。

2）设置弹簧光顺参数

激活弹簧光顺模型，相关参数设置位于Smoothing（光顺）标签下，可以设置的参数包括Spring Constant Factor（弹簧弹性系数）、Boundary Node Relaxation（边界点松弛因子）、Convergence Tolerance（收敛判据）和Number of Iterations（迭代次数）。

弹簧弹性系数应该在0 到1 之间变化，弹性系数等于0 时，弹簧系统没有耗散过程，在图中算例中，靠近壁面的网格没有被改变，而是保持了原来的网格形状和密度；在弹性系数等于1 时，弹簧系统的耗散过程与缺省设置相同，从图中可以发现壁面发生变形，壁面附近网格因为过度加密而质量下降。因此在实际计算中应该在0 到1 之间选择一个适当的值。边界点松弛因子用于控制动边界上网格点的移动。当这个值为零时，边界节点不发生移动；在这个值为1 时，则边界节点的移动计算中不采用松弛格式。在大多数情况下，这个值应该取为

0 到1 之间的一个值，以保证边界节点以合适的移动量发生移动。

收敛判据就是网格节点移动计算中，迭代计算的判据。迭代次数是指网格节点移动计算的最大迭代次数。

3）动态层

在Layering（动态层）标签下，可以设置与动态层模型相关的参数。通过设定Constant Height（常值高度）与Constant Ratio（常值比例）可以确定分解网格的两种方法。Split Factor （分割因子）和Collapse Factor（合并因子）则分别为上面介绍的alpha_s和alpha_c。

4）局部重新划分网格

在Remeshing（重划网格）标签下，设置与局部重划模型相关的参数。可以设置的参数包括Maximum Cell Skewness（最大畸变率）、Maximum Cell Volume（最大网格体积）和Minimum Cell Volume（最大网格体积），其含义如前所述，主要用于确定哪些网格需要被重新划分。在缺省设置中，如果重新划分的网格优于原网格，则用新网格代替旧网格；否则，将保持原网格划分不变。如果无论如何都要采用新网格的话，则可以在Options（选项）下面选择Must Improve Skewness（必须改善畸变率）选项。如果Options（选项）下面的Size Function（尺寸函数）被激活，则还可以用网格尺寸分布函数标志需要重新划分的网格。假设在某点附近的理想网格尺寸为L ，而某个网格的尺寸为L' ，如果：

L’不属于[0.8*gamma*L，1.25*gamma*L]

则网格被标志为需要重新划分的网格，并在随后的计算中被重新划分。式中的gamma用下面的公式计算：

当alpha>0时，gamma=1+alpha*d_b^(1+2*beta)

当alpha<0时，gamma=1+alpha*d_b^[(1-beta)^-1]

式中d_b为网格到壁面边界的最小距离，alpha和beta就是需要设置的Size Function Variation（尺寸函数增量）和Size Function Rate（尺寸函数变化率）。

以上两式中，alpha以边界网格尺寸为基准控制内部网格的大小，alpha等于0.5表示内部网格的尺寸至少是边界网格的1.5倍，alpha等于-0.5表示内部网格的尺寸等于边界网格的0.5倍，如果alpha等于0则表示内部网格与边界网格大小相同；bata的取值在-1到1之间，用于控制网格从边界到内部区域的变化速率。Beta取正值表示网格变化速率较慢取负值则表示变化速率较快，取0 则表示网格从边界到内部区域呈线性变化。

另外一个选项Size Function Resolution（尺寸函数分辨率）用于定义尺寸分布函数对网格大小的分辨率，其参照值是最小网格尺寸。

点击Use defaults（使用缺省值）按钮，可以恢复系统缺省设置。

5）设定活塞运动参数

如果在计算中选择使用In-cylinder（活塞）模型，需要指定Crank Shaft Speed（曲柄速度）、Starting Crank Speed（曲柄起始速度）、Crank Period（曲柄周期）以及Crank Angle Step Size（曲柄角度时间步长）。

FLUENT 中还提供一个内建函数用于计算活塞位置。这个函数的自变量为曲柄转角，如果选用这个函数，还需要指定Piston Stroke（活塞行程）和Connecting Rod Length（连接杆长度）二个参数。

注：以上的一些参数设定对动网格的更新质量影响很大，如果设置不当，可能引起更新不成功，如果不清楚这些设定因子对更新的影响，可以固定一些参数调节一两个参数查看效果。另外，在这些更新方法中比较容易出问题的就是Remeshing方法中的一些参数设定：Remeshing中的参数Minimum length scale和Maximum Length Scale，这两个参数你可以参考mesh scale info中的值，仅是参考，因为mesh scale info中的值是整个网格的评价值，设置的时候看一下动网格附近的网格和整个网格区域的大小比较，然后确定这两个参数，一般来讲，动网格附近的网格较密，这些值都比整体的小，所以在设置时通常设置为比mesh scale info中的Minimum length scale大一点，比Maximum Length Scale小一点。

5、定义动网格的运动方式

在计算动网格问题时，必须定义动网格区的运动方式。在动网格区为刚体运动时，可以用型函数和UDF 来定义其运动；在动网格区为变形区域时，则需要定义其几何特征及局部网格重划参数；如果动网格区既做刚体运动又有变形发生，则只能用UDF 来定义其几何形状的变化和运动过程。

上述定义在Dynamic Zones（动态区域）面板中设置，启动该面板的菜单操作顺序为：Define -> Dynamic Mesh -> Zones...

在这个面板中可以修改动态区域的设置、计算刚体运动区域的重心或删除一个动态区域。方法是首先在Dynamic Zones（动态区域）列表中选择一个动网格区，然后修改其设置参数，或计算其重心，或进行删除操作，最后点击Create（创建）按钮保存设置。

对于新加入的区域，需要先从Zone Names（区域名称）下选择相关区域，然后在Type（类型）下选择其运动类型。可供选择的运动类型包括Stationary（静止）、Rigid Body（刚体运动）、Deforming（变形）和User-Defined（用户自定义）四种。

1）静止区域设置

如果被指定区域为静止区域，则首先在Zone Names（区域名称）下选择这个区域，然在Type （类型）下选择Stationary（静止），再指定Adjacent Zone（相邻区域）的Cell Height网格高度）用于网格重新划分，最后点击Create（创建）按钮完成设置。

2）刚体运动区域设置

如果被指定区域为刚体运动区域，则其设置过程如下：

（1）在Zone Names（区域名称）下选择这个区域的名称，然后在Type（类型）下选择Rigid Body（刚体）。

（2）在Motion Attributes（运动属性）标签下的Motion UDF/Profile（用UDF 或型函数定义运动）中确定究竟用型函数，还是UDF 来做运动定义。

（3）在C.G. Location（重心位置）中定义刚体重心的初始位置。

（4）在C.G. Orentation（重心方向）中定义重力在惯性系中的方向。

（5）如果计算中包含活塞计算，则需要在Valve/Piston Axis（阀门或活塞轴）中指定阀门或活塞的参考轴。如果在所定义的网格区域中，某种形状的网格单元需要被排除在上述设置之外，则可以在Motion Mask（运动屏蔽）中选择这些单元形状。

（6）如果所定义的区域是面区域，则还需要定义Cell Height（网格高度）。这个参数用于定义局部网格重划时，与面区域相邻的网格的理想高度。

（7）点击Create（创建）按钮，完成设置。

3）变形运动

变形区域的设置过程为：

（1）在Zone Names（区域名称）中选择区域，并在Type（类型）下选择Deforming（变形）。（2）在Geometry Definition（几何定义）下定义变形区的几何特征，共有4 个选项：如果没有合适的几何形状，就在Definition 中选择none；如果变形区为平面，则选择plane（平面）并在Point on Plane（平面上的点）中定义平面上一点，同时在Plane Normal（平面法向）中定义法线方向；如果变形区为圆柱面，则选择Cylinder（圆柱），并同时定义Cylinder Radius（圆柱半径）、Cylinder Origin（圆柱原点）和Cylinder Axis（圆柱轴）；如果变形区几何形状需要用UDF 来定义，则在Definition（定义）中选择user-defined（用户定义），并在Geometry UDF（几何UDF）中选择适当的函数。

（3）在Remeshing Options（重划选项）标签下定义与网格局部重划相关的参数。重划方法在Mesh Methods（网格划分方法）中选择，其中包括Smoothing（弹簧光顺模型）、Layering （动态层模型）和Remeshing（网格重划模型）。如果动网格区域为面域，则需要设置局部

重划模型中的几个参数，包括Height（网格高度）、Height Factor（高度因子）和Maximum Skewness（最大畸变率）。如果动网格区域为体积域时，还可以设置Minimum Volume（最小体积）、Maximum Volume（最大体积）和Maximum Skewness（最大畸变率），这项设置可以覆盖系统缺省设置，从而改变网格重划过程。

（4）点击Create（创建）按钮完成设置。

4）用户定义的运动方式

对于同时存在运动和变形的区域，只能使用UDF 来定义其运动方式，定义步骤如下：（1）在Zone Names（区域名称）中选择需要定义的区域名称，并在Type（类型）下选择User-Defined（用户定义）。

（2）在Motion Attributes（运动属性）标签下，然后在Mesh Motion UDF（网格运动UDF）下选择相应的UDF 函数。

（3）点击Create（创建）按钮完成设置。

关于UDF 函数的编制方法请参见与UDF 相关的章节。

补充资料：FLUENT三种动边界控制实现方法

1) void DEFINE_CG_MOTION (UDFname,Dynamic_Thread * dt,real vel[ ], real omega[ ], real time,real dtime)。

此函数接口用于控制刚体的运动,用户把刚体质心运动速度和角速度分别赋值给vel和omega, FLUENT根据它们的值来自动计算出边界下一步的位置,从而实现动边界的控制; 刚体质心的位置可以在函数接口界面对话框中定义。Dynamic Zones中的dwall就是要控制的动边界,Motion UDF/Profile中的stc1sta010a0ph0就是UDFname,从中可看出它已被制定成用于控制dwall,理论上FLUEN T可以通过这种方式实现无穷多个动边界的控制;

C.G.Location用于设定初始位置的质心,C.G.Orientation用于设定刚体的初始角度。一般适用于刚体本身不变形的运动。

2) void DEFINE_GEOM(char name,Domain * d,Dynamic_Thread * dt,real * position)。此函数接口用于控制变形体的边界运动, position就是运动边界上某网格节点的位置值,用户可以通过对其赋值达到控制效果, position [0]对应边界节点的x坐标, position [1]对应y 坐标, position [2]对应z坐标; FLUENT自动遍历所有的边界节点,因此适用于有规律的可以用函数描述的运动边界。

3) void DEFINE_ GRID_MOTION(name,d,dt,time,dtime)。

此函数接口也用于控制形体的边界运动。主要用于更加复杂的控制,用户需要自己利用

FLUENT提供的其他函数来遍历运动边界上的节点,并对其位置进行控制,因此UDF编程比前面两种复杂得多。它甚至可以事先生成好边界数据,在计算中把数据读入,完成复杂形体控制。

6、预览动网格

在设置好动网格模型及动网格区的运动方式后，可以通过预览的方式检查设置效果。预览功能在Mesh Motion（网格运动）面板中进行设置，启动这个面板的菜单操作为：

Solve -> Mesh Motion...

预览操作步骤如下：

（1）在参数设置完毕后，首先保存算例（case）文件。因为与网格设置有关的记录都保存在算例文件中，在预览过程中伴随着网格的更新，与网格有关的记录也不断被刷新，如果不进行保存，则无法恢复原始设置状态，在发现参数设置问题后就无法进行更改了。

（2）设置迭代时间步数和时间步长。在计算过程中，当前时间将被显示在Current Mesh Time（当前网格时间）栏中。如果在计算中选用了活塞模型选项，则时间步长用曲柄速度（shaft speed）和曲柄转动步长（crank angel step size）计算得出。

（3）为了在图形窗口中预览网格变化过程，需要激活Display Options（显示选项）下的Display Grid（显示网格），并在Display Frequency（显示频率）中设置显示频率，即每分钟显示图幅数量。如果要保存显示的图形，则同时激活Save Hardcopy（保存硬拷贝）选项。（4）点击Preview（预览）按钮开始预览。

在定义了活塞运动时，活塞运动的预览是在IC Zone Motion（活塞运动）面板中实现的，激活这个面板的菜单操作为：

Display -> IC Zone Motion...

预览操作的具体步骤为：

（1）在Display Grid（显示网格）面板中选择准备预览的网格区域。

（2）在IC Zone Motion（网格运动）面板中，设置曲柄角度增量（Increment）和迭代步数（Number of Steps）。

（3）点击Preview（预览）按钮开始预览。

注：在动网格的建立内容中，已经说到了预览动网格在整个计算之中的重要性，请参考。7、定义活塞事件

在计算活塞运动时，可以通过定义活塞事件（events）的方式，在计算中增加问题的复杂度，使计算更接近真实情况。这些事件都是以曲柄角度为自变量的，比如可以将打开排气阀的时

间定义为曲柄角度的函数，然后在曲柄转到这个角度时，开始打开排气阀操作，等等。下面介绍定义活塞事件的具体步骤和事件种类，从这些介绍中可以了解活塞事件的基本概念和使用方法。

1）活塞事件的设置步骤

活塞事件在Dynamic Mesh Events（动网格事件）面板中定义。

启动这个面板的菜单操作次序为：

Define -> Dynamic Mesh -> Events...

活塞事件设置步骤为：

（1）在Number of Events（事件数量）栏中输入需要定义的事件数量。

（2）在Name（名称）下面输入事件名称。

（3）在At Crank Angel（曲柄角度）下输入事件发生的曲柄角度。

（4）点击Define（定义）按钮进入Define Event（定义事件）面板做详细定义

（5）在Type（类型）下面选择动网格运动类型，可供选择的运动类型包括Change Zone Type（改变网格区域类型）、Copy Zone BC（拷贝网格区域边界条件）、Create Sliding Interface（创建滑动网格交界面）、Delete Sliding Interface（删除滑动网格交界面）、Change Motion Attribute（改变运动属性）、Change Time Step Size（改变时间步长）、Insert Boundary Layer（插入边界层）、Remove Boundary Layer（删除边界层）、Insert Interior Layer（插入内部网格层）和Remove Interior Layer（删除内部网格层）等。在确定了事件类型后，还需要对所选运动类型做进一步设置，详细情况请参见下一小节中的内容。

（6）对其他事件重复（2）到（5）的设置过程。

（7）将所有事件设置完毕后，点击Apply（应用）按钮保存所有设置。

（8）在Dynamic Mesh Events（动网格事件）面板中点击Preview（预览）按钮，预览上述设置的事件。

2）活塞事件

在上一小节中讲到的各种事件类型的含义如下：

（1）Change Zone Type（改变网格区域类型）。在计算过程中可以改变网格区域的类型，即将网格类型在wall（壁面）、interface（交界面）、fluid（流体）和solid（固体）等类型之间进行转换。在New Zone Type（新网格区域类型）中指定一种类型，则当曲柄运动到相应位置时，指定网格区域类型就会发生转换。

（2）Copy Zone BC（拷贝网格区域边界条件）。在网格区域类型改变后，原来设置在网格

上的边界条件也要发生相应的变化，这种变化可以通过拷贝另一个区域的边界条件实现，即在设定改变网格区类型的同时，设定拷贝边界条件事件，在网格类型发生改变时，随即从另一个边界上拷贝其边界条件作为新类型下网格区域的边界条件。

（3）Creating Sliding Interface（创建滑动网格交界面）。在选择了这个事件类型后，需要定义构成交界面的两个面，即两个相邻网格区域Interface Zone 1（交界面区1）和Interface Zone 2（交界面区2）。在这里不用考虑界面类型转换的问题，系统会自动将定义为交界面区的网格区域转换为内部面类型。

（4）Deleting Sliding Interface（删除滑动网格交界面）。在选择了这个事件类型后，需要在Interface Name（界面名称）中指定将被删除的交界面。在删除界面后，界面所在网格区域将被自动转换为壁面，如果不准备采用系统缺省设置的壁面条件的话，可以考虑如（2）中提示的那样，从其他壁面上拷贝边界。

（5）Change Motion Attribute（改变运动属性）。可以用这个事件类型修改网格运动类型。选择了这个事件类型后，需要在Attribute（运动属性）中指定一种运动类型，可以选择的类型包括sliding（滑动）、moving（移动）和remesh（网格重划），同时在Status（状态）中选择enable（激活）或disable（关闭）。

（6）Change Time Step（改变时间步长）。在Crank Angel Step Size（曲柄角度步长）中输入新的步长，可以修改迭代的时间步长）。

（7）Insert Boundary Cell Layer（插入边界网格层）。选择这个事件类型后，需要指定Base Dynamic Zone（基础动网格区）和Side Dynamic Zone（临近动网格区）。边界层网格将以基础动网格区为起点生成，临近动网格区就是网格即将发生变形的区域。

（8）Deleting Boundary Cell Layer（删除边界网格层）。这个事件类型选项与（7）中的事件相对应，在选择这个选项后，只要指定Base Dynamic Zone（基础动网格区）就可以将插入的边界网格删除。

（9）Insert Interior Cell Layer（插入内部网格层）。如同插入边界网格层一样，这项操作也需要指定Base Dynamic Zone（基础动网格区）和Side Dynamic Zone（临近动网格区），不同的是还需要为新生成的界面命名，即指定Internal Zone 1 Name（内部网格区1 名称）和Internal Zone 2 Name（内部网格区2 名称）。

（10）Remove Interior Cell Layer（删除内部网格层）。选择这个事件选项后，只要指定被删除的网格区名称，Internal Zone 1 Name（内部网格区1 名称）和Internal Zone 2 Name （内部网格区2 名称）即可。

3）输入输出活塞事件

在Dynamic Mesh Events（动网格事件）面板上，点击Write...（输出）按钮可以将设置好的活塞事件以文件形式输出到磁盘上；点击Read...（输入）按钮可以将活塞事件文件读入计算模型。

4）活塞计算中初始网格的生成

可以在dynamesh-mesh-control 文本菜单下，用文本命令position-starting-mesh 设置活塞的初始位置，例如：

/define/models/dynamic-mesh-controls> position-starting-mesh

Start Crank Angle (deg) [0] 340

即将初始位置设置在曲柄转角为340°处。

在进行活塞计算的设置中，需要妥善处理动网格的拓扑关系、活塞和阀门的运动属性、进气阀和排气阀的开关以及活塞事件的先后次序，尽量全面、合理地进行设置，以保证计算能顺利完成。

动网格的入门专题

题记：在学习使用Fluent的时候，有不少朋友需要使用动网格模型(Dynamic Mesh Model)，因此，本版推出这个专题，进行大讨论，使大家在使用动网格时尽量少走弯路，更快更好地掌握；也欢迎使用过的版友积极参与讨论指导，谢谢！ 该专题主要包括以下的主要内容： ##1. 动网格的相关知识介绍； ##2. 以NACA0012翼型俯仰振荡实例进行讲解动网格的应用过程； ##3. 与动网格应用有关的参考文献； ##4. 使用动网格进行计算的一些例子。 ##1. 动网格的相关知识介绍 有关动网格基础方面的东西，请具体参考FLUENT User’s Guide或FLUENT全攻略的相关章节，这里只给出一些提要性的知识要点。 1、简介 动网格模型可以用来模拟流场形状由于边界运动而随时间改变的问题。边界的运动形式可以是预先定义的运动，即可以在计算前指定其速度或角速度；也可以是预先未做定义的运动，即边界的运动要由前一步的计算结果决定。网格的更新过程由FLUENT 根据每个迭代步中边界的变化情况自动完成。在使用动网格模型时，必须首先定义初始网格、边界运动的方式并指定参予运动的区域。可以用边界型函数或者UDF 定义边界的运动方式。FLUENT 要求将运动的描述定义在网格面或网格区域上。如果流场中包含运动与不运动两种区域，则需要将它们组合在初始网格中以对它们进行识别。那些由于周围区域运动而发生变形的区域必须被组合到各自的初始网格区域中。不同区域之间的网格不必是正则的，可以在模型设置中用FLUENT软件提供的非正则或者滑动界面功能将各区域连接起来。 注：一般来讲，在Fluent中使用动网格，基本上都要使用到UDF，所以你最好具备一定的C语言编程基础。 2、动网格更新方法 动网格计算中网格的动态变化过程可以用三种模型进行计算，即弹簧近似光滑模型(spring-based smoothing)、动态分层模型（dynamic layering）和局部重划模型（local remeshing）。 弹簧近似光滑模型 在弹簧近似光滑模型中，网格的边被理想化为节点间相互连接的弹簧。移动前的网格间距相当于边界移动前由弹簧组成的系统处于平衡状态。在网格边界节点发生位移后，会产生与位移成比例的力，力量的大小根据胡克定律计算。边界节点位移形成的力虽然破坏了弹簧系统原有的平衡，但是在外力作用下，弹簧系统经过调整将达到新的平衡，也就是说由弹簧连接在一起的节点，将在新的位置上重新获得力的平衡。从网格划分的角度说，从边界节点的位移出发，采用虎克定律，经过迭代计算，最终可以得到使各节点上的合力等于零的、新的网格节点位置，这就是弹簧光顺法的核心思想。 原则上弹簧光顺模型可以用于任何一种网格体系，但是在非四面体网格区域（二维非三角形），最好在满足下列条件时使用弹簧光顺方法： （1）移动为单方向。 （2）移动方向垂直于边界。 如果两个条件不满足，可能使网格畸变率增大。另外，在系统缺省设置中，只有四面体网格（三维）和三角形网格（二维）可以使用弹簧光顺法，如果想在其他网格类型中激活该模型，

网格技术发展方向的探讨

网格技术发展方向的探讨 [日期：2006-05-31] 来源：作者：[字体：大中小] 方莉英 四川机电职业技术学院信息工程系攀枝花617000 摘要该文主要介绍国内外一些有代表性的网格研究计划的开发进展和网格计算的起源及发展，探讨了网格技术的未来发展方向。 关键词网格网格计算服务网格； 1 引言 网络技术是一次计算革命，有可能将全球计算机联合起来协同工作，被人们视为21世纪的新型网络基础架构。网格的概念来源于电力网格，如电力网把强大的电力输送到每家每户的每个插座，用户不必过问电是水电、火电，还是核电，也不必过问电站位于何方。网格技术可把分布在各地的计算机联网，将充足的计算资源分配给每个用户，如同个人使用一台虚拟的超级计算机一样。但网格不是万能技术，不是所有问题都能解决，只有那些能并行运算的应用才可能被拆成若干细小任务分配到每个网格计算节点，通过并行处理来提高计算速度。 网格研究源于美国联邦政府过去10年来资助的高性能计算项目。这类研究使用的名词就是“网格”或“计算网格”。这类研究的目标是将跨地域的多台高性能计算机、大型数据库、贵重科研设备、通信设备、可视化设备和各种传感器整合成一个巨大的超级计算机系统，支持科学计算和科学研究。 2 国外网格技术研究动向 2.1 美国的网格技术研究 美国是网格研究起步最早的国家，美国多家研究机构开展了与网格相关的研究工作，制定了很多网格研究计划，如美国国家科学基金会资助的TeraGrid、美国国防部的“全球信息网格”(GIG)、美国政府资助的“大型物理实验网格”(GriPhyN)及美国能源部的ASCIGrid、国家技术网格(NTG)等计划。 (1)TeraGrid 美国国家科学基金会TeraGrid计划是一个需要多年分段实施的，为开放式科学研究而建立和使用的世界上最大、最全面的分布式基础设施。目前，TeraGrid已能提供40TFL OPS的聚合计算及近4千万亿字节的网络数据存储能力，并建成了能融合许多系统结构的

动网格

FLUENT6.1全攻略 图10-16 转子定子模型的静压等值线图 在显示速度矢量时，同样有绝对速度和相对速度两种形式。另外需要注意的是，后处理过程不能在交界区中的壁面、内部、周期等类型的边界上建立数据显示面（surface），但是可以在交界面上建立数据显示面，但结果将是单边的，就是只显示交界面一侧的结果。而且在跨越交接面时，等值线中可能会有细微的不连续。在画三维填充等值线时，图形中可能会出现一些小缝，但是这些缝只是图形显示问题，与解的连续性无关。 10.6 动网格模型 10.6.1 简介 动网格模型可以用来模拟流场形状由于边界运动而随时间改变的问题。边界的运动形式可以是预先定义的运动，即可以在计算前指定其速度或角速度；也可以是预先未做定义的运动，即边界的运动要由前一步的计算结果决定。网格的更新过程由FLUENT根据每个迭代步中边界的变化情况自动完成。在使用移动网格模型时，必须首先定义初始网格、边界运动的方式并指定参予运动的区域。可以用边界型函数或者UDF定义边界的运动方式。 FLUENT要求将运动的描述定义在网格面或网格区域上。如果流场中包含运动与不运动两种区域，则需要将它们组合在初始网格中以对它们进行识别。那些由于周围区域运动而发生变形的区域必须被组合到各自的初始网格区域中。不同区域之间的网格不必是正则 21

FLUENT6.1全攻略 22 的，可以在模型设置中用FLUENT 软件提供的非正则或者滑动界面功能将各区域连接起来。 10.6.2 动网格守恒方程 在任意一个控制体中，广义标量Φ的积分守恒方程为： ()∫∫∫∫??+??=??+V V Φg V V dV S A d ΦΓA d u u ρΦV ρΦd dt d G G G G （10-7） 式中ρ为流体密度，u G 为速度向量，g u G 移动网格的网格速度，Γ为扩散系数，ΦS 为 源项，V ?代表控制体V 的边界。 方程（10-7）中的时间导数项，可以用一阶后向差分格式写成： ()()t V V dV dt d n n V ΔΦ?Φ=Φ+∫ρρρ1 （10-8） 式中n 和n+1代表不同的时间层。n+1层上的V n+1由下式计算： t dt dV V V n n Δ+=+1 （10-9） 式中dV/dt 是控制体的时间导数。为了满足网格守恒定律，控制体的时间导数由下式计算： ∫∑??=?=V n j j j g g f A u A d u dt dV G G G G , （10-10） 式中n f 是控制体积的面网格数，j A G 为面j 的面积向量。点乘j j g A u G G ?,由下式计算： t V A u j j j g Δ=?δG G , （10-11） 式中j δV 为控制体积面j 在时间间隔Δt 中扫过的空间体积。 10.6.3 动网格更新方法 动网格计算中网格的动态变化过程可以用三种模型进行计算，即弹簧光滑模型、动态层模型和局部重划模型。

动网格流沙版完美整理.

动网格 让网格动起来（1）—闲谈动网格 在固体有限元计算中，网格运动实非什么稀奇事儿。而且在绝多数固体计算的基本物理量是网格的节点位移，所以，固体计算中，网格节点运动是对的，没有运动反而不正常了。也可以这么说：正因为计算域内部节点间的相对运动，才导致了内应力的产生。 流体计算与固体完全不同。其根源在于它们使用的网格类型不同。当前固体有限元计算采用的是拉格朗日网格，而流体计算则大多数采用的欧拉网格。如果说把拉格朗日网格中的节点点看作是真实世界的物质原子的话，那么欧拉网格的节点则好比是真实世界中的一个个传感器，它们总是呆在相同的位置，真实的记录着各自位置上的物理量。正常情况下，欧拉网格系统是这样的：计算域和节点保持位置不变，发生变化的是物理量，网格节点就像一个个布置在计算域中的传感器，记录该位置上的物理量。这其实是由流体力学研究方法所决定的。宏观与微观的差异决定了固体力学计算采用拉格朗日网格，流体计算采用欧拉网格。关于这部分的详细解说，可以参阅任何一本计算流体动力学书籍。 世界是公平的。有利必有弊。朗格朗日网格适合计算节点位移，然而对于过大的网格变形却难以处理。欧拉网格生来可以处理大变形（因为节点不动），然而对于对于节点运动的处理，则是其直接软肋。然而很不幸的是，现实生活中有太多网格边界运动的实例。如汽车发动机中的气缸运动、阀门开启与关闭、机翼的运动、飞机投弹等等等等举不胜举。 计算流体动力学计算的基本物理量通常为：速度、温度、压力、组分。并不计算网格节点位移。因此要让网格产生运动，通常给节点施加的物理约束是速度。CFD中的动网格大体分为两类：（1）显式规定的网格节点速度。配合瞬态时间，即可很方便的得出位移。当然一些求解器（如FLUENT）也支持稳态动网格，这时候可以直接指定节点位移。（2）网格节点速度是通过求解得到的。如6DOF模型基本上都属于此类。用户将力换算成加速度，然后将其积分成速度。 对于第一类动网格问题，在fluent中通常可以使用profile与UDF进行网格设置，通过规定节点或区域的速度、角速度或位移等方式来显式确定网格的运动，通常大部分的动网格问题都归于此类。而对于第二类问题，通常涉及到力的计算，力在流体中通常是对压力进行积分而来。将力转换为速度或位移，一般涉及到加速度、转动惯量等物理量的计算。在fluent 中，可以使用6DOF模型进行处理，在CFX中，可以使用刚体模型（13.0以上版本才有）。在FLUENT中，动网格涉及的内容包括： （1）运动的定义。主要是PROFILE文件与UDF中的动网格宏。 （2）网格更新。FLUENT中关于网格更新方法有三种：网格光顺、动态层、网格重构。 需要详细了解这些网格更新方法的运作机理，每个参数所代表的具体含义及设置方法，每种方法的适用范围。 动网格的最在挑战来自于网格更新后的质量，避免负体积是动网格调试的主要目标。在避免负网格的同时，努力提高运动更新后的网格质量。

结构化网格和非结构化网格

1. 什么是结构化网格和非结构化网格 1.1结构化网格 从严格意义上讲，结构化网格是指网格区域内所有的内部点都具有相同的毗邻单元。 它可以很容易地实现区域的边界拟合，适于流体和表面应力集中等方面的计算。它的主要优点是： 网格生成的速度快。 网格生成的质量好。 数据结构简单。 对曲面或空间的拟合大多数采用参数化或样条插值的方法得到，区域光滑，与实际的模型更容易接近。 它的最典型的缺点是适用的范围比较窄，只适用于形状规则的图形。尤其随着近几年的计算机和数值方法的快速发展，人们对求解区域的几何形状的复杂性的要求越来越高，在这种情况下，结构化网格生成技术就显得力不从心了。 1.2非结构化网格 同结构化网格的定义相对应，非结构化网格是指网格区域内的内部点不具有相同的毗邻单元。即与网格剖分区域内的不同内点相连的网格数目不同。从定义上可以看出，结构化网格和非结构化网格有相互重叠的部分，即非结构化网格中可能会包含结构化网格的部分。 2.如果一个几何造型中既有结构化网格，也有非结构化网格，分块完成的，分别生成网格后，也可以直接就调入fluent中计算。 3.在fluent中，对同一个几何造型，如果既可以生成结构化网格，也可生成非结构化网格，当然前者要比后者的生成复杂的多，那么应该选择哪种网格，两者计算结果是否相同，哪个的计算结果更好些呢？ 一般来说，结构网格的计算结果比非结构网格更容易收敛，也更准确。但后者容易做。 影响精度主要是网格质量，和你是用那种网格形式关系并不是很大，如果结构话网格的质量很差，结果同样不可靠，相对而言，结构化网格更有利于计算机存储数据和加快计算速度。

结构化网格据说计算速度快一些，但是网格划分需要技巧和耐心。非结构化网格容易生成，但相对来说速度要差一些。 4.在gambit中，只有map和submap生成的是结构化网格，其余均为非结构化网格。 采用分块网格划分的时候，在两个相邻块之间设置了connected，但是这两个块我要用不同尺寸的网格来划分。比如说我用结构化的六面体网格来划分，一遍的尺寸为2，另一边的尺寸为3，这时候公共边界面该怎么处理？如果采用cooper 的格式来划分这个网格，尺寸就是前面所说的，该怎么来做呢？ 我用单独的两个块试过，就是在公共边界上采用interface的格式，但是由于与这个公共边界相邻的另一个边界也不得不用interface格式，结果导入fluent 的时候就说can not creat a bound loop,也不清楚这是什么问题。 如果中间面两侧的面网格一致，可以直接在fluent中merge,如果不一致，可以设interface 网格的正交性是指三个方向上的网格边之间互相垂直的程度。一般而言，三维网格单元中，三个方向上的网格边之间的夹角越接近90度则质量越好。这一点在规则区域(例如正方形方腔)很容易实现，但对于流动区域比较复杂的问题则非常困难。但一般情况下，应当保证所有的网格单元内的网格边夹角大于10度，否则网格本身就会引入较大的数值误差。 EquiSize Skew（尺寸扭曲率）和EquiAngle Skew（角度扭曲率）是评判网格质量最主要标准，其值越小，网格质量越高 一般来说，Fluent要求扭曲率3D小于0.85，2D小于0.75。 关于复杂模型和gambit中的实体及虚体 模型比较复杂，是在pro/E中建的模，然后用igs导入gambit，不过这样就产生了很多碎线和碎面并且在一些面交界的地方还存在尖角。我曾经做成功过把它们统统merge成一个虚面，中间设置了一个可以容忍尖角的参数，也可以划分网格，但把生成的msh文件导入fluent就会出错，这是virtual geometry的原因还是因为尖角的原因？还有，virtual geometry和普通的真实的几何体到底有什么区别？好像最大的区别是virtual geometry不能进行布尔操作，布尔操作（boolean operation）又是什么？使用virtual geometry需要注意哪些问题？virtual geometry是很头疼的问题。你把它们统统merge成一个虚面 按理说全是虚的也是可以算的。可能是因为尖角的原因，虚实最大差别：是virtual geometry不能进行布尔操作，boolean operation即是并 对于复杂外形的网格生成，不可避免的会用到virtual geometry，virtual face ,和virtual edge等， 1。作网格的时候，把所有的面全部合成一个虚面的做法不好，特别是对于复杂外形的网格生成，你最好在模型变化剧烈的地方多分几个面，这样会更有效的控制网格能够在模型表面曲率比较大的地方能够生成规则的结构或者非结构网格。

Fluent动网格专题讨论

Fluent动网格专题讨论（-） 题记：在学习使用Fluent的时候，有不少朋友需要使用动网格模型(Dynamic Mesh Model)，因此，本版推出这个专题，进行大讨论，使大家在使用动网格时尽量少走弯路，更快更好地掌握；也欢迎使用过的版友积极参与讨论指导，谢谢！ 该专题主要包括以下的主要内容： ##1. 动网格的相关知识介绍； ##2. 以NACA0012翼型俯仰振荡实例进行讲解动网格的应用过程； ##3. 与动网格应用有关的参考文献； ##4. 使用动网格进行计算的一些例子。 ##1. 动网格的相关知识介绍 有关动网格基础方面的东西，请具体参考FLUENT User’s Guide或FLUENT全攻略的相关章节，这里只给出一些提要性的知识要点。 1、简介 动网格模型可以用来模拟流场形状由于边界运动而随时间改变的问题。边界的运动形式可以是预先定义的运动，即可以在计算前指定其速度或角速度；也可以是预先未做定义的运动，即边界的运动要由前一步的计算结果决定。网格的更新过程由FLUENT 根据每个迭代步中边界的变化情况自动完成。在使用动网格模型时，必须首先定义初始网格、边界运动的方式并指定参予运动的区域。可以用边界型函数或者UDF 定义边界的运动方式。FLUENT 要求将运动的描述定义在网格面或网格区域上。如果流场中包含运动与不运动两种区域，则需要将它们组合在初始网格中以对它们进行识别。那些由于周围区域运动而发生变形的区域必须被组合到各自的初始网格区域中。不同区域之间的网格不必是正则的，可以在模型设置中用FLUENT软件提供的非正则或者滑动界面功能将各区域连接起来。 注：一般来讲，在Fluent中使用动网格，基本上都要使用到UDF，所以你最好具备一定的C语言编程基础。 2、动网格更新方法 动网格计算中网格的动态变化过程可以用三种模型进行计算，即弹簧近似光滑模型(spring-based smoothing)、动态分层模型（dynamic layering）和局部重划模型（local remeshing）。 弹簧近似光滑模型 在弹簧近似光滑模型中，网格的边被理想化为节点间相互连接的弹簧。移动前的网格间距相当于边界移动前由弹簧组成的系统处于平衡状态。在网格边界节点发生位移后，会产生与位移成比例的力，力量的大小根据胡克定律计算。边界节点位移形成的力虽然破坏了弹簧系统原有的平衡，但是在外力作用下，弹簧系统经过调整将达到新的平衡，也就是说由弹簧连接在一起的节点，将在新的位置上重新获得力的平衡。从网格划分的角度说，从边界节点的位移出发，采用虎克定律，经过迭代计算，最终可以得到使各节点上的合力等于零的、新的网格节点位置，这就是弹簧光顺法的核心思想。 原则上弹簧光顺模型可以用于任何一种网格体系，但是在非四面体网格区域（二维非三角形），最好在满足下列条件时使用弹簧光顺方法： （1）移动为单方向。 （2）移动方向垂直于边界。 如果两个条件不满足，可能使网格畸变率增大。另外，在系统缺省设置中，只有四面体网格（三维）和三角形网格（二维）可以使用弹簧光顺法，如果想在其他网格类型中激活该模型，需要在dynamic-mesh-menu 下使用文字命令spring-on-all-shapes?，然后激活该选项即可。 动态层模型 对于棱柱型网格区域（六面体和或者楔形），可以应用动态层模型。动态层模型的中心思想是根据紧邻运动边界网格层高度的变化，添加或者减少动态层，即在边界发生运动时，如果紧邻边界的网格层高度增大到一定程度，就将其划分为两个网格层；如果网格层高度降低到一定程度，就将紧邻边界的两个网格层合并为一个层： 如果网格层j扩大，单元高度的变化有一临界值：

FLUENT 动网格教程

FLUENT动网格教程 摘自https://www.sodocs.net/doc/f74210573.html,/dvbbs/dispbbs.asp?boardid=61&id=1396题记：在学习使用Fluent的时候，有不少朋友需要使用动网格模型(Dynamic Me sh Model)，因此，本版推出这个专题，进行大讨论，使大家在使用动网格时尽量少走弯路，更快更好地掌握；也欢迎使用过的版友积极参与讨论指导，谢谢！。 该专题主要包括以下的主要内容： §一、动网格的相关知识介绍； §二、以NACA0012翼型俯仰振荡实例进行讲解动网格的应用过程； §三、与动网格应用有关的参考文献； §四、使用动网格进行计算的一些例子。 §一、动网格的相关知识介绍 有关动网格基础方面的东西，请具体参考FLUENT User’s Guide或FLUENT全攻略的相关章节，这里只给出一些提要性的知识要点。 1、简介 动网格模型可以用来模拟流场形状由于边界运动而随时间改变的问题。边界的运动形式可以是预先定义的运动，即可以在计算前指定其速度或角速度；也可以是预先未做定义的运动，即边界的运动要由前一步的计算结果决定。网格的更新过程由FLUE NT 根据每个迭代步中边界的变化情况自动完成。在使用动网格模型时，必须首先定义初始网格、边界运动的方式并指定参予运动的区域。可以用边界型函数或者UDF定义边界的运动方式。FLUENT 要求将运动的描述定义在网格面或网格区域上。如果流场中包含运动与不运动两种区域，则需要将它们组合在初始网格中以对它们进行识别。那些由于周围区域运动而发生变形的区域必须被组合到各自的初始网格区域中。不同区域之间的网格不必是正则的，可以在模型设置中用FLUENT软件提供的非正则或者滑动界面功能将各区域连接起来。 注：一般来讲，在Fluent中使用动网格，基本上都要使用到UDF，所以你最好具备一定的C语言编程基础。

ICEM CFD动网格学习

FLUENT动网格教程 该专题主要包括以下的主要内容： §一、动网格的相关知识介绍； §二、以NACA0012翼型俯仰振荡实例进行讲解动网格的应用过程； §三、与动网格应用有关的参考文献； §四、使用动网格进行计算的一些例子。 §一、动网格的相关知识介绍 有关动网格基础方面的东西，请具体参考FLUENT User’s Guide或FLUENT全攻略的相关章节，这里只给出一些提要性的知识要点。 1、简介 动网格模型可以用来模拟流场形状由于边界运动而随时间改变的问题。边界的运动形式可以是预先定义的运动，即可以在计算前指定其速度或角速度；也可以是预先未做定义的运动，即边界的运动要由前一步的计算结果决定。网格的更新过程由FLUE NT 根据每个迭代步中边界的变化情况自动完成。在使用动网格模型时，必须首先定义初始网格、边界运动的方式并指定参予运动的区域。可以用边界型函数或者UDF定义边界的运动方式。FLUENT 要求将运动的描述定义在网格面或网格区域上。如果流场中包含运动与不运动两种区域，则需要将它们组合在初始网格中以对它们进行识别。那些由于周围区域运动而发生变形的区域必须被组合到各自的初始网格区域中。不同区域之间的网格不必是正则的，可以在模型设置中用FLUENT软件提供的非正则或者滑动界面功能将各区域连接起来。 注：一般来讲，在Fluent中使用动网格，基本上都要使用到UDF，所以你最好具备一定的C语言编程基础。 2、动网格更新方法 动网格计算中网格的动态变化过程可以用三种模型进行计算，即 弹簧近似光滑模型(spring-based smoothing)、 动态分层模型（dynamic layering） 局部重划模型（local remeshing）

动网格技术的研究

Fluent中被动运动Six DOF (2016-01-30 10:07:51) 转载▼ 利用CFD软件解决动网格问题，通常可分为以下两类： （1）主动型动网格 主动型动网格问题通常指的是边界运动规律及运动状态已知，可由软件使用者通过函数或程序进行描述。在程序计算过程中，求解器调用动网格运动轨迹描述程序实现边界运动。这类动网格例子很多，如各类泵、风扇等。 （2）被动型动网格 被动型动网格问题通常是指其边界运动规律往往是未知的，常常需要通过计算边界上的力或力矩，以此来求取动网格的运动规律。在这类动网格计算设置中，网格变化规律难以预料，导致网格参数经常需要进行多次调整才能达到目的。这类例子在现实中也很多，比如风力发电机的叶轮、水轮机等。 解决主动型动网格问题比较容易，利用CFD软件提供的动网格模拟能力很容易解决。需要关注的地方是边界运动后，网格节点如何重新布置和生成。如在FLUENT软件中，其动网格主要包括三种网格功能：Smoothing、Layering、Remeshing。利用Remeshing 功能几乎可以解决所有主动型动网格问题。 那被动型动网格问题怎么处理呢？一般来说，这类边界的运动都是由于内部流体对其压力所造成的，那么就涉及到力和力矩计算的问题。对于这类问题，在FLUENT软件中可以采用Six DOF模型进行计算。

需要注意的是，以上所有类型动网格计算均以刚性移动边界为前提，即不考虑由于流体作用力导致的边界变形。若要计算边界变形，则需要采用流固耦合方法，利用固体求解器计算被动型动网格中的力和力矩，均是压力对面的积分计算而来。 下面介绍一下Six DOF运用方法： 在FLUENT中利用Six DOF是需要定义UDF宏的。该宏的定义形式如下： DEFINE_SDOF_PROPERTIES(name, properties, dt ,time ,dtime) 函数中： Name：宏名称 Real *properties：存储6DOF属性的数组 Dynamic_Thread *dt：存储制定的动网格属性 Real time：当前时间 Real dtime：时间步长 该UDF宏没有返回值。用户需要定义的变变量为name、properties、dt、dtime。在利用该UDF宏的过程中,需要注意的是properties数组，其包含多种属性：SDOF_MASS //质量，定义方式如properties[SDOF_MASS]=7.0定义质量7kg SDOF_IXX，SDOF_IYY，SOF_IZZ，//X、Y、Z方向惯性矩 SDOF_IXY，SDOF_IXZ，SDOF_IYZ，//惯性积 SDOF_LOAD_LOCAL，//布尔值，FALSE为全局坐标系，TRUE为体坐标系（局部坐标系），默认为FALSE SDOF_LOAD_F_X，SDOF_LOAD_F_Y，SDOF_LOAD_F_Z，//X，Y，Z方向外力

FLUENT动网格技术简介

FLUENT动网格简介 在固体有限元计算中，网格运动实非什么稀奇事儿。而且在绝多数固体计算的基本物理量是网格的节点位移，所以，固体计算中，网格节点运动是对的，没有运动反而不正常了。也可以这么说：正因为计算域内部节点间的相对运动，才导致了内应力的产生。 流体计算与固体完全不同。其根源在于它们使用的网格类型不同。当前固体有限元计算采用的是拉格朗日网格，而流体计算则大多数采用的欧拉网格。如果说把拉格朗日网格中的节点点看作是真实世界的物质原子的话，那么欧拉网格的节点则好比是真实世界中的一个个传感器，它们总是呆在相同的位置，真实的记录着各自位置上的物理量。正常情况下，欧拉网格系统是这样的：计算域和节点保持位置不变，发生变化的是物理量，网格节点就像一个个布置在计算域中的传感器，记录该位置上的物理量。这其实是由流体力学研究方法所决定的。宏观与微观的差异决定了固体力学计算采用拉格朗日网格，流体计算采用欧拉网格。关于这部分的详细解说，可以参阅任何一本计算流体动力学书籍。 世界是公平的。有利必有弊。朗格朗日网格适合计算节点位移，然而对于过大的网格变形却难以处理。欧拉网格生来可以处理大变形（因为节点不动），然而对于对于节点运动的处理，则是其直接软肋。然而很不幸的是，现实生活中有太多网格边界运动的实例。如汽车发动机中的气缸运动、阀门开启与关闭、机翼的运动、飞机投弹等等等等举不胜举。 计算流体动力学计算的基本物理量通常为：速度、温度、压力、组分。并不计算网格节点位移。因此要让网格产生运动，通常给节点施加的物理约束是速度。CFD中的动网格大体分为两类：（1）显式规定的网格节点速度。配合瞬态时间，即可很方便的得出位移。当然一些求解器（如FLUENT）也支持稳态动网格，这时候可以直接指定节点位移。（2）网格节点速度是通过求解得到的。如6DOF 模型基本上都属于此类。用户将力换算成加速度，然后将其积分成速度。 对于第一类动网格问题，在fluent中通常可以使用profile与UDF进行网格设置，通过规定节点或区域的速度、角速度或位移等方式来显式确定网格的运动，通常大部分的动网格问题都归于此类。而对于第二类问题，通常涉及到力的计算，力在流体中通常是对压力进行积分而来。将力转换为速度或位移，一般涉及到加

Fluent中的动网格

Fluent中的动网格 动网格是目前求解计算域变化问题的常用方法。参考Fluent帮助，可以知道动网格技术与一般流动计算设置的主要区别在于网格更新方法和更新域设置。这里就这两方面问题的一点体会作一简单记录。 一、网格更新方法 弹簧近似光滑法将任意两网格节点之间的连线理想地看成一条弹簧，并通过近似弹簧的压缩或拉伸实现网格和计算域的改变。该方法网格拓扑不变，无需网格的插值处理，对结构化（四边形、六面体）和非结构化（三角形、四面体）网格同样适用。但不适合于大变形情况，当计算区域变形较大时，变形后的网格质量变差，严重影响计算精度。 动态分层法在运动边界相邻处根据运动规律动态增加或减少网格层数，以此来更新变形区域的网格。该方法适用于结构化网格，通过设置适当的分层和缩减系数，更新后的网格依然为较为均匀的结构化网格，对计算精度影响较小。对于运动域具有多自由度和任意变形情况，该方法处理起来非常困难。 网格重生方法在整个网格更新区域内依据设定的最大和最小网格尺寸判断需要进行网格重生的网格，并依据设置的更新频率进行网格重生处理。该方法适用于非结构化网格，能够较好的应用于任意变形的计算区域处理。 二、更新域设置 更新域设置是动网格设置中的一项重要工作，最常用的设置是刚体运动域和变形域，这里针对这两种域的设置注意事项和技巧作一简单介绍。 1、域动网格 一般来讲，设置为刚体运动域的区域一般为壁面类边界，通过设置固壁的运动，模拟计算域内物体的运动。 由于固壁边界有时形状较为复杂，壁面附近网格尺度与周围网格尺度存在较大差别，网格更新时变形较大。在这种情况下，可以设置一个包含固壁运动边界的计算域，通过该计算域的整体运动模拟域内物体的运动，在有的地方将这种方法称为域动网格法。在域动网格法中，需要设置包含运动物体的内部计算域、内部计算域界面均为刚体运动域。如下图所示。

径向基函数在动网格中的应用及可并行性研究

4 ·技术?/? TECHNOLOGY ·科研信息化技术与应用2012, 3(5): 4–12 径向基函数在动网格中的应用及可并行性研究 马文鹏，陆忠华，胡晓东 中国科学院计算机网络信息中心 超级计算中心，北京 100190摘 要： 关键词： Applications of Radial Basis Functions in Dynamic Mesh and Its Parallelizability Ma Wenpeng, Lu Zhonghua, Hu Xiaodong Supercomputing Center, Computer Network Information Center , Chinese Academy of Sciences , Beijing 100190, China Abstract: 径向基函数广泛应用于网格变形、气动外形优化设计、网格优化等领域。近年来，基于径向基函数的动网格技术得到了深入的研究和广泛的应用。本文结合计算流体力学和高性能计算的应用背景，从径向基函数对网格的变形质量和变形效率进行了总结和进一步研究：在网格变形方面，重点对比了不同基函数对同一网格运动变形能力和同一基函数对不同网格运动的适应能力；在网格变形效率方面，分析了算法在计算和存储的瓶颈所在，考虑了? OpenMP 和?GPU 这两种共享内存的加速方式，得到较好加速比。最后，分析了当网格规模增大时，动网格在分布式计算和存储模型?(MPI) 下的处理方法。 径向基函数；径向函数；动网格；网格变形；并行 Radial basis functions are widely used for mesh deformation, aerodynamic shape optimization, grid optimization, etc. Recently, dynamic mesh techniques based on radial basis functions are widely used and much more attention is paid to them. This paper summaries the quality and efficiency of mesh deformation with the applications of radial basis functions on both CFD (Computational Fluid Dynamics) and HPC (High Performance Computing). In addition, a further research is conducted on the following two aspects: the abilities to satisfy the same mesh motion among different radial basis functions and the abilities to satisfy different mesh motion using one radial basis function; computation and storage bottleneck are analyzed and then parallel solutions based on shared memory models including OpenMP and GPU are considered 基金项目：国家高技术研究发展计划 (863 计划) (2012AA01A304)；国家自然科学基金 (91130019)

结构化网格与非结构化网格比较

对于连续的物理系统的数学描述，如航天飞机周围的空气的流动，水坝的应力 集中等等,通常是用偏微分方程来完成的。为了在计算机上实现对这些物理系统的行为或状态的模 拟，连续的方程必须离散化，在方程的求解域上（时间和空间）仅仅需要有限个点，通过 计算这些点上的未知变量既而得到整个区域上的物理量的分布。有限差分，有限体积和有限元等数值方法都是通过这种方法来 实现的。这些数值方法的非常重要的一个部分就是实现对求解区域的网格剖分。 网格剖分技术已经有几十年的发展历史了。到目前为止，结构化网格技术发展 得相对比较成熟，而非结构化网格技术由于起步较晚，实现比较困难等方面的 原因，现在正在处于逐渐走向成熟的阶段。下面就简要介绍一些这方面的情况。 从严格意义上讲，结构化网格是指网格区域内所有的内部点都具有相同的毗邻 单元。 结构化网格生成技术有大量的文献资料[1,2,3,4]。结构化网格有很多 优点: 1.它可以很容易地实现区域的边界拟合，适于流体和表面应力集中等方面的计算。 5.对曲面或空间的拟合大多数采用参数化或样条插值的方法得到，区域光滑， 与实际的模型更容易接近。 它的最典型的缺点是适用的范围比较窄。尤其随着近几年的计算机和数值方法 的快速发展，人们对求解区域的复杂性的要求越来越高，在这种情况下，结构 化网格生成技术就显得力不从心了。 代数网格生成方法。主要应用参数化和插值的方法，对处理简单的求解区域十 分有效。 同结构化网格的定义相对应，非结构化网格是指网格区域内的内部点不具有相 同的毗邻单元。即与网格剖分区域内的不同内点相连的网格数目不同。从定义 上可以看出，结构化网格和非结构化网格有相互重叠的部分，即非结构化网格 中可能会包含结构化网格的部分。 非结构化网格技术从六十年代开始得到了发展,主要是弥补结构化网格不能够 解决任意形状和任意连通区域的网格剖分的缺欠.到90年代时,非结构化网格的文献达到了它的高峰时期.由于非结构化网格的生成技术比较复杂,随着人们对 求解区域的复杂性的不断提高,对非结构化网格生成技术的要求越来越高.从现 在的文献调查的情况来看,非结构化网格生成技术中只有平面三角形的自动生

动网格相关知识--1

https://www.sodocs.net/doc/f74210573.html,/dvbbs/dispbbs.asp?boardid=61&Id=1396 题记：在学习使用Fluent的时候，有不少朋友需要使用动网格模型(Dynamic Mesh M odel)，因此，本版推出这个专题，进行大讨论，使大家在使用动网格时尽量少走弯路，更快更好地掌握；也欢迎使用过的版友积极参与讨论指导，谢谢！ 该专题主要包括以下的主要内容： §一、动网格的相关知识介绍； §二、以NACA0012翼型俯仰振荡实例进行讲解动网格的应用过程； §三、与动网格应用有关的参考文献； §四、使用动网格进行计算的一些例子。 §一、动网格的相关知识介绍 有关动网格基础方面的东西，请具体参考FLUENT User’s Guide或FLUENT全攻略的相关章节，这里只给出一些提要性的知识要点。 1、简介 动网格模型可以用来模拟流场形状由于边界运动而随时间改变的问题。边界的运动形式可以是预先定义的运动，即可以在计算前指定其速度或角速度；也可以是预先未做定义的运动，即边界的运动要由前一步的计算结果决定。网格的更新过程由FLUENT 根据每个迭代步中边界的变化情况自动完成。在使用动网格模型时，必须首先定义初始网格、边界运动的方式并指定参予运动的区域。可以用边界型函数或者UDF 定义边界的运动方式。FLUENT 要求将运动的描述定义在网格面或网格区域上。如果流场中包含运动与不运动两种区域，则需要将它们组合在初始网格中以对它们进行识别。那些由于周围区域运动而发生变形的区域必须被组合到各自的初始网格区域中。不同区域之间的网格不必是正则的，可以在模型设置中用FLUENT软件提供的非正则或者滑动界面功能将各区域连接起来。 注：一般来讲，在Fluent中使用动网格，基本上都要使用到UDF，所以你最好具备一定的C语言编程基础。 2、动网格更新方法 动网格计算中网格的动态变化过程可以用三种模型进行计算，即弹簧近似光滑模型(spr ing-based smoothing)、动态分层模型（dynamic layering）和局部重划模型（local reme shing）。 弹簧近似光滑模型 原则上弹簧光顺模型可以用于任何一种网格体系，但是在非四面体网格区域（二维非三

Fluent动网格----layering个一个简单实例(作者Snow)

Fluent动网格----layering个一个简单实例我这几天看了点动网格技术方面的东西，在学习过程中发现这方面的例子很少，自己也走了一些弯路。现在还好，弄明白了一些，能够应付现在我的工作。为了让更多学习者快速了解动网格，我打算尽量把我学习心得在这里和大家分享，这里给出一个layering的一个简单例子。 1.Gambit画网格 本例很简单，在Gambit里画一个10*10的矩形，网格间隔为1，也就是有100个网格，具体见下图。都学动网格的人了，不至于这个不会做！ 这里需要注意一个问题：设置边界条件的时候，一定要把要移动的边单独设定，本例中一右边界作为移动的边，设成wall就可以，这里再后面需要制定。 2.编写UDF #include "udf.h" #include "unsteady.h" #include "stdio.h" #include "stdlib.h" /************************************************************/ real current_time = 0.0 ; Domain * domain ; Thread * thread ; real NV_VEC( origin ),NV_VEC( force ),NV_VEC( moment ) ; /************************************************************/ DEFINE_CG_MOTION(throttle,dt,vel,omega,time,dtime) { current_time = CURRENT_TIME ; vel[0] = 30; Message("time=%f omega=%f\n",current_time) ; }

fluent 动网格

Remeshing方法中的一些参数设定：Remeshing中的参数Minimum length scale和Maximum Length Scale，这两个参数你可以参考mesh scale info中的值，仅是参考，因为mesh scale info中的值是整个网格的评价值，设置的时候看一下动网格附近的网格和整个网格区域的大小比较，然后确定这两个参数，一般来讲，动网格附近的网格较密，这些值都比整体的小，所以在设置时通常设置为比mesh scale info中的Minimum length scale大一点，比Maximum Length Scale小一点。 以上是一般来讲的设置思路。下面是我在NACA0012翼型动网格例子中的设置： Remeshing中的参数设定： 为了得到较好的网格更新，本例在使用局部网格重新划分方法时，使用尺寸函数，也就是Remeshing+Must Improve Skewness+Size Function的策略。 将Minimum Length Scale及Maximum Length Scale均设置为0，为了使所有的区域都被标记重新划分； Maximum Cell Skewness（最大单元畸变），参考Mesh Scale Info…中的参考值0.51，将其设定为0.4，以保证更新后的单元质量； Size Remesh Interval（依照尺寸标准重新划分的间隔），将这个值设定为1，在FLUENT，不满足最大网格畸变的网格在每个时间步都会被标记，而后重新划分，而不满足最小，最大及尺寸函数的网格，只有在Current Time=（Size Remesh Interval）*delta t的时候，才根据这些尺寸的标准标记不合格的单元进行重新划分，为了保证每步的更新质量，将其修改为1，就是每个时间都根据尺寸的标准标记及更新网格。 Size Function Resolution（尺寸函数分辨率），保持默认的3； Size Function Variation（尺寸函数变量）：建议使用一个小值，在0.1到0.5之间，本例将其设置为0.3；Size Function Rate（尺寸函数变化率），保持默认的0.3。 动网格(dynamic mesh)是CFD中专有的概念。由于当前流体计算多采用欧拉坐标系，该坐标系区别于拉格朗日坐标系的一个最直观特点是：计算过程中网格保持静止。因此，在CFD计算中应用动网格，具有其特别的难处。 1、动网格控制方式 最主要的困难在于边界运动后的网格质量控制。由于边界的运动，不可避免的导致网格变形。我们知道，求解器对于网格质量的容忍是有限度的。当网格扭曲过大引起网格质量的急剧下降，可能导致计算发散、形成负网格，进而终止计算。因此，在边界运动过程中，对网格质量进行控制尤为重要。在fluent软件中动网格主要有三种控制方式：smoothing，layering，remeshing。其中layering主要应用与四边形网格及六面体网格，remeshing主要应用于三角形网格及四面体网格等费结构网格中，至于smoothing方法则在各类网格中均可应用。 layering方法应用于结构网格也是有条件限制的：边界运动最好是沿着某单一方向。如若是旋转，最好还是采用非结构网格配合remeshing方式。 非结构网格是最适合应用动网格模型的，但是网格质量不好控制，通常需要仔细调节。结构网格采用layering 方法，能够很好的控制网格质量，但是几何适应性差。具体采用何种网格类型以及何种控制方式，还是要