Fluent理论手册3—滑移网格及动网格理论
- 格式:pdf
- 大小:462.42 KB
- 文档页数:9
图1采用滑移网格
图2采用动网格
滑移网格与动网格都可以计算瞬态运动问题。
但是存在以下区别:
(1)滑移网格需要建立多个域,实际上还是计算区域运动,是一个独立区域内所有网格一起运动,而动网格则是真正意义上的网格运动。
(2)滑移网格不会涉及到网格的变形与重生,但是要设计到交界面设置
(3)滑移网格不会造成负体积,而动网格极易形成负网格
(4)滑移网格是一种简化了的模型,最大计算误差出现在交界面位置
(5)动网格误差常出现在运动壁面位置,因此实际应用中往往将边界层与壁面合在一起运动。
(6)SRF,MRF与MP只能计算稳态,在新版本的FLUENT中应用的是坐标系变换,和滑移网格的设置有所区别。
第一章Fluent 软件的介绍fluent 软件的组成:软件功能介绍:GAMBIT 专用的CFD 前置处理器(几何/网格生成) Fluent4.5 基于结构化网格的通用CFD 求解器 Fluent6.0 基于非结构化网格的通用CFD 求解器 Fidap 基于有限元方法的通用CFD 求解器 Polyflow 针对粘弹性流动的专用CFD 求解器 Mixsim 针对搅拌混合问题的专用CFD 软件 Icepak专用的热控分析CFD 软件软件安装步骤:step 1: 首先安装exceed软件,推荐是exceed6.2版本,再装exceed3d,按提示步骤完成即可,提问设定密码等,可忽略或随便填写。
step 2: 点击gambit文件夹的setup.exe,按步骤安装;step 3: FLUENT和GAMBIT需要把相应license.dat文件拷贝到FLUENT.INC/license目录下;step 4:安装完之后,把x:\FLUENT.INC\ntbin\ntx86\gambit.exe命令符拖到桌面(x为安装的盘符);step 5: 点击fluent源文件夹的setup.exe,按步骤安装;step 6: 从程序里找到fluent应用程序,发到桌面上。
注:安装可能出现的几个问题:1.出错信息“unable find/open license.dat",第三步没执行;2.gambit在使用过程中出现非正常退出时可能会产生*.lok文件,下次使用不能打开该工作文件时,进入x:\FLUENT.INC\ntbin\ntx86\,把*.lok文件删除即可;3.安装好FLUENT和GAMBIT最好设置一下用户默认路径,推荐设置办法,在非系统分区建一个目录,如d:\usersa) win2k用户在控制面板-用户和密码-高级-高级,在使用fluent用户的配置文件修改本地路径为d:\users,重起到该用户运行命令提示符,检查用户路径是否修改;b) xp用户,把命令提示符发送到桌面快捷方式,右键单击命令提示符快捷方式在快捷方式-起始位置加入D:\users,重起检查。
滑移网格法滑移网格法作为众多网格方法中的一种,在CFD 计算方面应用十分广泛,特别是在处理有旋转域的相关问题上。
但是“滑移”在中文的传统意义上是属于“动”的一种,使得很多人将滑移网格法错误的认为是动网格方法的一种。
以下将具体介绍动网格的三种实现-方法,以及滑移网格和滑移网格方法的具体实现步骤。
1 动网格方法动网格可以用来模拟流场形状由于边界运动而随时间改变的问题。
其在商用软件中的应用十分广泛。
-动网格可以分为以下三种:弹簧近似光顺法、动态分层法、局部网格分层法。
1.1 弹簧近似光顺法弹簧近似光顺法近似将网格节点间通过弹簧相连,任意一个网格节点的位移均会打破网格系统的力的平衡,通过反复迭代,得到一个新的网格系统,变化过程中网格总量不变。
但在处理计算域较大变形问题时,误差较大。
1.2 动态分层法动态分层法是根据移动边界的运动规律,在变形区域实时增加或减少网格。
局部网格发生增减。
相对弹簧近似光顺法,可以处理计算域变形较大的问题。
1.3 局部网格分层法局部网格分层法是在弹性光顺法的基础上发展得到的。
在弹性光顺法得到新的网格系统的基础上,删除部分网格并从新生成。
该方法有一个网格拉伸度尺寸标准,弹性光顺法得到的网格满足标准则继续使用原网格,不满足则从新生成。
由上可见,动网格方法在数值模拟计算域形状发生变化的相关问题上具有较好的效果。
但在处理计算域旋转但形状不变的问题上有所不足。
而滑移网格方法可以很好的弥补这一不足。
2 滑移网格法滑移网方法是在计算过程中,移动单元区域沿网格分界面滑动,移动网格区域内部网格保持不变。
这一特点使得其在数值模拟带有旋转区域相关问题时,具有较大的优势。
2.1 滑移网格法实现步骤(1)读取各计算计算域网格,识别转/静交接面;(2)对转/静交界面上的网格节点进行外延,构造滑移边界;(3)找到每个滑移点的宿主单元,并计算对应的插值型函数;(4)进行流场的定常数值计算,达到收敛标准;(5)开始非定常计算,第一个物理时刻t=0;(6)对旋转域进行相应旋转,重新构造滑移边界、宿主单元搜索和插值函数计算;(7)进行第t 个物理时刻计算,直到满足内迭代收敛标准;判断非定常计算是否完成,完成则终止计算;否则开始下一物理时刻t=t+1时刻的计算,返回(6)。
第一章Fluent 软件的介绍fluent 软件的组成:软件功能介绍:GAMBIT 专用的CFD 前置处理器(几何/网格生成) Fluent4.5 基于结构化网格的通用CFD 求解器 Fluent6.0 基于非结构化网格的通用CFD 求解器 Fidap 基于有限元方法的通用CFD 求解器 Polyflow 针对粘弹性流动的专用CFD 求解器 Mixsim 针对搅拌混合问题的专用CFD 软件 Icepak专用的热控分析CFD 软件软件安装步骤:step 1: 首先安装exceed软件,推荐是exceed6.2版本,再装exceed3d,按提示步骤完成即可,提问设定密码等,可忽略或随便填写。
step 2: 点击gambit文件夹的setup.exe,按步骤安装;step 3: FLUENT和GAMBIT需要把相应license.dat文件拷贝到FLUENT.INC/license目录下;step 4:安装完之后,把x:\FLUENT.INC\ntbin\ntx86\gambit.exe命令符拖到桌面(x为安装的盘符);step 5: 点击fluent源文件夹的setup.exe,按步骤安装;step 6: 从程序里找到fluent应用程序,发到桌面上。
注:安装可能出现的几个问题:1.出错信息“unable find/open license.dat",第三步没执行;2.gambit在使用过程中出现非正常退出时可能会产生*.lok文件,下次使用不能打开该工作文件时,进入x:\FLUENT.INC\ntbin\ntx86\,把*.lok文件删除即可;3.安装好FLUENT和GAMBIT最好设置一下用户默认路径,推荐设置办法,在非系统分区建一个目录,如d:\usersa) win2k用户在控制面板-用户和密码-高级-高级,在使用fluent用户的配置文件修改本地路径为d:\users,重起到该用户运行命令提示符,检查用户路径是否修改;b) xp用户,把命令提示符发送到桌面快捷方式,右键单击命令提示符快捷方式在快捷方式-起始位置加入D:\users,重起检查。
网格的读入和使用FLUENT可以从输入各种类型,各种来源的网格。
你可以通过各种手段对网格进行修改,如:转换和调解节点坐标系,对并行处理划分单元,在计算区域内对单元重新排序以减少带宽以及合并和分割区域等。
你也可以获取网格的诊断信息,其中包括内存的使用与简化,网格的拓扑结构,解域的信息。
你可以在网格中确定节点、表面以及单元的个数,并决定计算区域内单元体积的最大值和最小值,而且检查每一单元内适当的节点数。
以下详细叙述了FLUENT关于网格的各种功能。
(请参阅网格适应一章以详细了解网格适应的具体内容。
)网格拓扑结构FLUENT是非结构解法器,它使用内部数据结构来为单元和表面网格点分配顺序,以保持临近网格的接触。
因此它不需要i,j,k指数来确定临近单元的位置。
解算器不会要求所有的网格结构和拓扑类型,这使我们能够灵活使用网格拓扑结构来适应特定的问题。
二维问题,可以使用四边形网格和三角形网格,三维问题,可以使用六面体、四面体,金字塔形以及楔形单元,具体形状请看下面的图形。
FLUENT可以接受单块和多块网格,以及二维混合网格和三维混合网格。
另外还接受FLUENT有悬挂节点的网格(即并不是所有单元都共有边和面的顶点),有关悬挂节点的详细信息请参阅“节点适应”一节。
非一致边界的网格也可接受(即具有多重子区域的网格,在这个多重子区域内,内部子区域边界的网格节点并不是同一的)。
详情请参阅非一致网格Figure 1: 单元类型可接受网格拓扑结构的例子正如网格拓扑结构一节所说,FLUENT可以在很多种网格上解决问题。
图1—11所示为FLUENT的有效网格。
O型网格,零厚度壁面网格,C型网格,一致块结构网格,多块结构网格,非一致网格,非结构三角形,四边形和六边型网格都是有效的。
Note that while FLUENT does not require a cyclic branch cut in an O-type grid,it will accept a grid that contains one.Figure 1: 机翼的四边形结构网格Figure 2:非结构四边形网格Figure 3: 多块结构四边形网格Figure 4: O型结构四边形网格Figure 5: 降落伞的零厚度壁面模拟Figure 6: C型结构四边形网格Figure 7:三维多块结构网格Figure 8: Unstructured Triangular Grid for an AirfoilFigure 9:非结构四面体网格Figure 10:具有悬挂节点的混合型三角形/四边形网格Figure 11:非一致混合网格for a Rotor-Stator Geometry选择适当的网格类型FLUENT在二维问题中可以使用由三角形、四边形或混合单元组成的网格,在三维问题中可以使用四面体,六面体,金字塔形以及楔形单元,或者两种单元的混合。
FLUENT动网格技术简介FLUENT动网格简介在固体有限元计算中,网格运动实非什么稀奇事儿。
而且在绝多数固体计算的基本物理量是网格的节点位移,所以,固体计算中,网格节点运动是对的,没有运动反而不正常了。
也可以这么说:正因为计算域内部节点间的相对运动,才导致了内应力的产生。
流体计算与固体完全不同。
其根源在于它们使用的网格类型不同。
当前固体有限元计算采用的是拉格朗日网格,而流体计算则大多数采用的欧拉网格。
如果说把拉格朗日网格中的节点点看作是真实世界的物质原子的话,那么欧拉网格的节点则好比是真实世界中的一个个传感器,它们总是呆在相同的位置,真实的记录着各自位置上的物理量。
正常情况下,欧拉网格系统是这样的:计算域和节点保持位置不变,发生变化的是物理量,网格节点就像一个个布置在计算域中的传感器,记录该位置上的物理量。
这其实是由流体力学研究方法所决定的。
宏观与微观的差异决定了固体力学计算采用拉格朗日网格,流体计算采用欧拉网格。
关于这部分的详细解说,可以参阅任何一本计算流体动力学书籍。
世界是公平的。
有利必有弊。
朗格朗日网格适合计算节点位移,然而对于过大的网格变形却难以处理。
欧拉网格生来可以处理大变形(因为节点不动),然而对于对于节点运动的处理,则是其直接软肋。
然而很不幸的是,现实生活中有太多网格边界运动的实例。
如汽车发动机中的气缸运动、阀门开启与关闭、机翼的运动、飞机投弹等等等等举不胜举。
计算流体动力学计算的基本物理量通常为:速度、温度、压力、组分。
并不计算网格节点位移。
因此要让网格产生运动,通常给节点施加的物理约束是速度。
CFD中的动网格大体分为两类:(1)显式规定的网格节点速度。
配合瞬态时间,即可很方便的得出位移。
当然一些求解器(如FLUENT)也支持稳态动网格,这时候可以直接指定节点位移。
(2)网格节点速度是通过求解得到的。
如6DOF模型基本上都属于此类。
用户将力换算成加速度,然后将其积分成速度。
对于第一类动网格问题,在fluent中通常可以使用profile与UDF 进行网格设置,通过规定节点或区域的速度、角速度或位移等方式来显式确定网格的运动,通常大部分的动网格问题都归于此类。