1、动网格学习资料汇总动网格学习资料汇总为了提高流体传动与控制研究室科研助理对FLUENT软件的学习效率,现将有关动网格学习资料进行汇总,资料来源于书籍和网上资源,不能以偏概全,希望科研助理在学习的过程中,多搜集相关资料,作为后续补充,逐渐完善,为以后科研助理的学习提供帮助。 1、动网格简介动网格模型可以用来模拟流场形状由于边界运动而随时间改变的问题。边界的运动形式可以是预先定义的运动,即可以在计算前指定其速度或角速度;也可以是预先未做定义的运动,即边界的运动要由前一步的计算结果决定。网格的更新过程由FLUENT根据每个迭代步中边界的变化情况自动完成。在使用动网格模型时,必须首先定义初始网格、边界
2、运动的方式并指定参予运动的区域。可以用边界型函数或者UDF定义边界的运动方式。FLUENT 要求将运动的描述定义在网格面或网格区域上。如果流场中包含运动与不运动两种区域,则需要将它们组合在初始网格中以对它们进行识别。那些由于周围区域运动而发生变形的区域必须被组合到各自的初始网格区域中。不同区域之间的网格不必是正则的,可以在模型设置中用FLUENT软件提供的非正则或者滑动界面功能将各区域连接起来。注:一般来讲,在Fluent中使用动网格,基本上都要使用到UDF,所以你最好具备一定的C语言编程基础。2、动网格更新方法动网格计算中网格的动态变化过程可以用三种模型进行计算,即弹簧近似光滑模型(spri
3、ng-based smoothing)、动态分层模型(dynamic layering)和局部重划模型(local remeshing)。弹簧近似光滑模型原则上弹簧光顺模型可以用于任何一种网格体系,但是在非四面体网格区域(二维非三角形),最好在满足下列条件时使用弹簧光顺方法: (1)移动为单方向。(2)移动方向垂直于边界。如果两个条件不满足,可能使网格畸变率增大。另外,在系统缺省设置中,只有四面体网格(三维)和三角形网格(二维)可以使用弹簧光顺法,如果想在其他网格类型中激活该模型,需要在dynamic-mesh-menu 下使用文字命令spring-on-all-shapes,然后激活该选项即
4、可。动态分层模型动态分层模型的应用有如下限制:(1)与运动边界相邻的网格必须为楔形或者六面体(二维四边形)网格。(2)在滑动网格交界面以外的区域,网格必须被单面网格区域包围。(3)如果网格周围区域中有双侧壁面区域,则必须首先将壁面和阴影区分割开,再用滑动交界面将二者耦合起来。(4)如果动态网格附近包含周期性区域,则只能用FLUENT 的串行版求解,但是如果周期性区域被设置为周期性非正则交界面,则可以用FLUENT的并行版求解。如果移动边界为内部边界,则边界两侧的网格都将作为动态层参与计算。如果在壁面上只有一部分是运动边界,其他部分保持静止,则只需在运动边界上应用动网格技术,但是动网格区与静止网
5、格区之间应该用滑动网格交界面进行连接。局部网格重划模型需要注意的是,局部网格重划模型仅能用于四面体网格和三角形网格。在定义了动边界面以后,如果在动边界面附近同时定义了局部重划模型,则动边界上的表面网格必须满足下列条件:(1)需要进行局部调整的表面网格是三角形(三维)或直线(二维)。(2)将被重新划分的面网格单元必须紧邻动网格节点。(3)表面网格单元必须处于同一个面上并构成一个循环。(4)被调整单元不能是对称面(线)或正则周期性边界的一部分。动网格的实现在FLUENT 中是由系统自动完成的。如果在计算中设置了动边界,则FLUENT 会根据动边界附近的网格类型,自动选择动网格计算模型。如果动边界附
6、近采用的是四面体网格(三维)或三角形网格(二维),则FLUENT 会自动选择弹簧光顺模型和局部重划模型对网格进行调整。如果是棱柱型网格,则会自动选择动态层模型进行网格调整。在静止网格区域则不进行网格调整。动网格问题中对于固体运动的描述,是以固体相对于重心的线速度和角速度为基本参数加以定义的。既可以用型函数定义固体的线速度和角速度,也可以用UDF 来定义这两个参数。同时需要定义的是固体在初始时刻的位置。注:这一小节主要讲述了动网格的更新方法,最好能掌握,尤其是各种方法的适用范围,通常来讲,在一个case中,我们使用的更新方法都是根据网格类型以及和要实现的运动来选择的,很多时候都是几种更新方法搭配
7、起来使用的。总结一下:使用弹簧近似光滑法网格拓扑始终不变,无需插值,保证了计算精度。但弹簧近似光滑法不适用于大变形情况,当计算区域变形较大时,变形后的网格会产生较大的倾斜变形,从而使网格质量变差,严重影响计算精度。动态分层法在生成网格方面具有快速的优势,同时它的应用也受到了一些限制。它要求运动边界附近的网格为六面体或楔形,这对于复杂外形的流场区域是不适合的。使用局部网格重划法要求网格为三角形(二维)或四面体(三维),这对于适应复杂外形是有好处的,局部网格重划法只会对运动边界附近区域的网格起作用。3、动网格问题的建立设置动网格问题的步骤中需要注意的如下:在Solver(求解器)面板中选择非定常流
8、(unsteady)计算;设定边界条件,即设定壁面运动速度;激活动网格模型,并设定相应参数;指定移动网格区域的运动参数;保存算例文件和数据文件;预览动网格设置,菜单操作为:Solve - Mesh Motion.;应用自动保存功能保存计算结果,在动网格计算中,因为每个计算步中网格信息都会改变,而网格信息是储存在算例文件中的,所以必须同时保存算例文件和数据文件;如果想建立网格运动的动画过程,可以在Solution Animation(计算结果动画)面板中进行相关设置。注:在这一步中,需要提醒一下,使用动网格进行正式计算之前,最好养成预览动网格更新的习惯;就是在正式计算前,浏览一下动网格的更新情况
9、,这样可以避免在计算过程中出现动网格更新本身的问题。在预览更新时,很多人都说会出现负体积的警告,更新不成功,出现这样的问题时,最好先把时间步长改的更小点儿试试,一般来讲,排除UDF本身的原因,出现更新出错的原因都与时间步长有关,这需要结合所使用的更新方法多琢磨。4、设定动网格参数为了使用动网格模型,需要在dynamic mesh(动网格)面板中激活Dynamic Mesh(动网格)选项。如果计算的是活塞运动,则同时激活In-Cylinder(活塞)选项。然后选择动网格模型,并设置相关参数。选择网格更新模型在Mesh Methods(网格划分方法)下面选择Smothing(弹簧光顺模型),Lay
10、ering(动态层模型)和(或)Remshing(局部网格重划模型)。设置弹簧光顺参数激活弹簧光顺模型,相关参数设置位于Smoothing(光顺)标签下,可以设置的参数包括Spring Constant Factor(弹簧弹性系数)、Boundary Node Relaxation(边界点松弛因子)、Convergence Tolerance(收敛判据)和Number of Iterations(迭代次数)。弹簧弹性系数应该在0 到1 之间变化,弹性系数等0时,弹簧系统没有耗散过程,在图中算例中,靠近壁面的网格没有被改变,而是保持了原来的网格形状和密度;在弹性系数等于1 时,弹簧系统的耗散过程
11、与缺省设置相同,从图中可以发现壁面发生变形,壁面附近网格因为过度加密而质量下降。因此在实际计算中应该在0 到1 之间选择一个适当的值。边界点松弛因子用于控制动边界上网格点的移动。当这个值为零时,边界节点不发生移动;在这个值为1 时,则边界节点的移动计算中不采用松弛格式。在大多数情况下,这个值应该取为0 到1 之间的一个值,以保证边界节点以合适的移动量发生移动。收敛判据就是网格节点移动计算中,迭代计算的判据。迭代次数是指网格节点移动计算的最大迭代次数。动态分层在Layering(动态分层)标签下,可以设置与动态层模型相关的参数。通过设定Constant Height(常值高度)与Constant
12、 Ratio(常值比例)可以确定分解网格的两种方法。Split Factor(分割因子)和Collapse Factor(合并因子)则分别为上面介绍的alpha_s和alpha_c。局部网格重新划分网格在Remeshing(重划网格)标签下,设置与局部重划模型相关的参数。可以设置的参数包括Maximum Cell Skewness(最大畸变率)、Maximum Cell Volume(最大网格体积)和Minimum Cell Volume(最大网格体积),其含义如前所述,主要用于确定哪些网格需要被重新划分。在缺省设置中,如果重新划分的网格优于原网格,则用新网格代替旧网格;否则,将保持原网格划分
13、不变。如果无论如何都要采用新网格的话,则可以在Options(选项)下面选择Must Improve Skewness(必须改善畸变率)选项。如果Options(选项)下面的Size Function(尺寸函数)被激活,则还可以用网格尺寸分布函数标志需要重新划分的网格。设定活塞运动参数如果在计算中选择使用In-cylinder(活塞)模型,需要指定Crank Shaft Speed(曲柄速度)、Starting Crank Speed(曲柄起始速度)、Crank Period(曲柄周期)以及Crank Angle Step Size(曲柄角度时间步长)。FLUENT 中还提供一个内建函数用于计
14、算活塞位置。这个函数的自变量为曲柄转角,如果选用这个函数,还需要指定Piston Stroke(活塞行程)和Connecting Rod Length(连接杆长度)二个参数。注:以上的一些参数设定对动网格的更新质量影响很大,如果设置不当,可能引起更新不成功,如果不清楚这些设定因子对更新的影响,可以固定一些参数调节一两个参数查看效果。另外,在这些更新方法中比较容易出问题的就是Remeshing方法中的一些参数设定:Remeshing中的参数Minimum length scale和Maximum Length Scale,这两个参数你可以参考mesh scale info中的值,仅是参考,因为m
15、esh scale info中的值是整个网格的评价值,设置的时候看一下动网格附近的网格和整个网格区域的大小比较,然后确定这两个参数,一般来讲,动网格附近的网格较密,这些值都比整体的小,所以在设置时通常设置为比mesh scale info中的Minimum length scale大一点,比Maximum Length Scale小一点。5、定义动网格的运动方式在计算动网格问题时,必须定义动网格区的运动方式。在动网格区为刚体运动时,可以用型函数和UDF 来定义其运动;在动网格区为变形区域时,则需要定义其几何特征及局部网格重划参数;如果动网格区既做刚体运动又有变形发生,则只能用UDF 来定义其几
16、何形状的变化和运动过程。上述定义在Dynamic Zones(动态区域)面板中设置,在这个面板中可以修改动态区域的设置、计算刚体运动区域的重心或删除一个动态区域。方法是首先在Dynamic Zones(动态区域)列表中选择一个动网格区,然后修改其设置参数,或计算其重心,或进行删除操作,最后点击Create(创建)按钮保存设置。对于新加入的区域,需要先从Zone Names(区域名称)下选择相关区域,然后在Type(类型)下选择其运动类型。可供选择的运动类型包括Stationary(静止)、Rigid Body(刚体运动)、Deforming(变形)和User-Defined(用户自定义)四种。
17、补充资料: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中的就是UDFname,从中可看出它已被制定
18、成用于控制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
19、2对应z坐标; FLUENT自动遍历所有的边界节点,因此适用于有规律的可以用函数描述的运动边界。3) void DEFINE_ GRID_MOTION(name,d,dt,time,dtime)。此函数接口也用于控制形体的边界运动。主要用于更加复杂的控制,用户需要自己利用 FLUENT提供的其他函数来遍历运动边界上的节点,并对其位置进行控制,因此 UDF编程比前面两种复杂得多。它甚至可以事先生成好边界数据,在计算中把数据读入,完成复杂形体控制。6、预览动网格在设置好动网格模型及动网格区的运动方式后,可以通过预览的方式检查设置效果。预览功能在Mesh Motion(网格运动)面板中进行设置,启动
20、这个面板的菜单操作为:Solve - Mesh Motion.预览操作步骤如下:(1)在参数设置完毕后,首先保存算例(case)文件。因为与网格设置有关的记录都保存在算例文件中,在预览过程中伴随着网格的更新,与网格有关的记录也不断被刷新,如果不进行保存,则无法恢复原始设置状态,在发现参数设置问题后就无法进行更改了。(2)设置迭代时间步数和时间步长。在计算过程中,当前时间将被显示在Current Mesh Time(当前网格时间)栏中。如果在计算中选用了活塞模型选项,则时间步长用曲柄速度(shaft speed)和曲柄转动步长(crank angel step size)计算得出。(3)为了在图
21、形窗口中预览网格变化过程,需要激活Display Options(显示选项)下的Display Grid(显示网格),并在Display Frequency(显示频率)中设置显示频率,即每分钟显示图幅数量。如果要保存显示的图形,则同时激活Save Hardcopy(保存硬拷贝)选项。(4)点击Preview(预览)按钮开始预览。在定义了活塞运动时,活塞运动的预览是在IC Zone Motion(活塞运动)面板中实现的,激活这个面板的菜单操作为:Display - IC Zone Motion.预览操作的具体步骤为:(1)在Display Grid(显示网格)面板中选择准备预览的网格区域。(2)
22、在IC Zone Motion(网格运动)面板中,设置曲柄角度增量(Increment)和迭代步数(Number of Steps)。(3)点击Preview(预览)按钮开始预览。注:在动网格的建立内容中,已经说到了预览动网格在整个计算之中的重要性,请参考。7、定义活塞事件在计算活塞运动时,可以通过定义活塞事件(events)的方式,在计算中增加问题的复杂度,使计算更接近真实情况。这些事件都是以曲柄角度为自变量的,比如可以将打开排气阀的时间定义为曲柄角度的函数,然后在曲柄转到这个角度时,开始打开排气阀操作,等等。下面介绍定义活塞事件的具体步骤和事件种类,从这些介绍中可以了解活塞事件的基本概念和
23、使用方法。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 Ty
24、pe(改变网格区域类型)、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(删除内部网格层)等。在确定了事件类
25、型后,还需要对所选运动类型做进一步设置,详细情况请参见下一小节中的内容。(6)对其他事件重复(2)到(5)的设置过程。(7)将所有事件设置完毕后,点击Apply(应用)按钮保存所有设置。(8)在Dynamic Mesh Events(动网格事件)面板中点击Preview(预览)按钮,预览上述设置的事件。2)活塞事件在上一小节中讲到的各种事件类型的含义如下:(1)Change Zone Type(改变网格区域类型)。在计算过程中可以改变网格区域的类型,即将网格类型在wall(壁面)、interface(交界面)、fluid(流体)和solid(固体)等类型之间进行转换。在New Zone Typ
26、e(新网格区域类型)中指定一种类型,则当曲柄运动到相应位置时,指定网格区域类型就会发生转换。(2)Copy Zone BC(拷贝网格区域边界条件)。在网格区域类型改变后,原来设置在网格上的边界条件也要发生相应的变化,这种变化可以通过拷贝另一个区域的边界条件实现,即在设定改变网格区类型的同时,设定拷贝边界条件事件,在网格类型发生改变时,随即从另一个边界上拷贝其边界条件作为新类型下网格区域的边界条件。(3)Creating Sliding Interface(创建滑动网格交界面)。在选择了这个事件类型后,需要定义构成交界面的两个面,即两个相邻网格区域Interface Zone 1(交界面区1)和
27、InterfaceZone 2(交界面区2)。在这里不用考虑界面类型转换的问题,系统会自动将定义为交界面区的网格区域转换为内部面类型。(4)Deleting Sliding Interface(删除滑动网格交界面)。在选择了这个事件类型后,需要在Interface Name(界面名称)中指定将被删除的交界面。在删除界面后,界面所在网格区域将被自动转换为壁面,如果不准备采用系统缺省设置的壁面条件的话,可以考虑如(2)中提示的那样,从其他壁面上拷贝边界。(5)Change Motion Attribute(改变运动属性)。可以用这个事件类型修改网格运动类型。选择了这个事件类型后,需要在Attrib
28、ute(运动属性)中指定一种运动类型,可以选择的类型包括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(临近动网格区)。边界层网格将以基础动网格区
29、为起点生成,临近动网格区就是网格即将发生变形的区域。(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 名
30、称)和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 文本菜单下,用文本命令设置活塞的初始位置,例如:/define/models/dynamic-mesh-controls position-starting-meshStart Crank Angle (deg) 0 340即将初始位置设置在曲柄转角为340处。在进行活塞计算的设置中,需要妥善处理动网格的拓扑关系、活塞和阀门的运动属性、进气阀和排气阀的开关以及活塞事件的先后次序,尽量全面、合理地进行设置,以保证计算能顺利完成。8、动网格实例旋转式动网格FLUENT技术基础与应用实例中动网格实例
copyright@ 2008-2022 冰豆网网站版权所有
经营许可证编号:鄂ICP备2022015515号-1