《FLUENT知识点解析》由会员分享,可在线阅读,更多相关《FLUENT知识点解析(10页珍藏版)》请在金锄头文库上搜索。
1、本文格式为Word版,下载可任意编辑FLUENT知识点解析 一、根本设置 1Double Precision的选择 启动设置如图,这里着重说说Double Precision(双精度)复选框,对于大多数处境,单精度求解器已能很好的得志精度要求,且计算量小,这里我们选择单精度。然而对于以下一些特定的问题,使用双精度求解器可能更有利1。 a. 几何特征包含某些极端的尺度(如分外长且窄的管道),单精度求解器可能不能足够精确地表达各尺度方向的节点信息。 b. 假设几何模型包含多个通过小直径管道相互连接的体,而某一个区域的压力更加大(由于用户只能设定一个总体的参考压力位置),此时,双精度求解器可能更能表
2、达压差带来的滚动(如渐缩渐扩管的无粘与可压缩滚动模拟)。 c. 对于某些高导热系数比或高宽纵比的网格,使用单精度求解器可能会遇到收敛性不佳或精确度缺乏缺乏的问题,此时,使用双精度求解器可能会有所扶助。 2网格光顺化 用光滑和交换的方式改善网格:通过Mesh下的Smooth/Swap来实现,可用来提高网格质量,一般用于三角形或四边形网格,不过质量提高的效果一般般,影响较小,网格质量的提高主要还是在网格生成软件里面实现,所以这里不再用光滑和交换的方式改善网格,其原理可参考FLUENT全攻略(已下载)。 3Pressure-based与Density-based 求解器设置如图。下面说一说Press
3、ure-based和Density-based的识别: Pressure-Based Solver是Fluent的优势,它是基于压力法的求解器,使用的是压力修正算法,求解的操纵方程是标量形式的,擅长求解不成压缩滚动,对于可压滚动也可以求解;Fluent 6.3以前的版本求解器,只有Segregated Solver和Coupled Solver,其实也是Pressure-Based Solver的两种处理方法; Density-Based Solver是Fluent 6.3新进展出来的,它是基于密度法的求解器,求解的操纵方程是矢量形式的,主要离散格式有Roe,AUSM+,该方法的初衷是让Flu
4、ent具有对比 1 李鹏飞,徐敏义,王飞飞.精通CFD工程仿真与案例实战:FLUENT GAMBIT ICEM CFD TecplotM. 北京,人民邮电出版社,2022:114-116 好的求解可压缩滚动才能,但目前格式没有添加任何限制器,因此还不太完善;它只有Coupled的算法;对于低速问题,他们是使用Preconditioning方法来处理,使之也能够计算低速问题。Density-Based Solver下断定是没有SIMPLEC,PISO这些选项的,由于这些都是压力修正算法,不会在这种类型的求解器中展现的;一般还是使用Pressure-Based Solver解决问题。 基于压力的求
5、解器适用于求解不成压缩和中等程度的可压缩流体的滚动问题。而基于密度的求解器最初用于高速可压缩滚动问题的求解。虽然目前两种求解器都适用于各类滚动问题的求解(从不成压缩滚动到高度可压缩滚动),但对于高速可压缩滚动而言,使用基于密度的求解器通常能获得比基于压力的求解器更为精确的结果。 4axisymmetric和axisymmetric swirl 从字面的意思很好理解axisymmetric和axisymmetric swirl的区别: axisymmetric:是轴对称的意思,也就是关于一个坐标轴对称,2D的axisymmetric问题仍为2D问题。而axisymmetric swirl:是轴对
6、称旋转的意思,就是一个区域关于一条坐标轴回转所产生的区域,这产生的将是一个回转体,是3D的问题。在Fluent中使用这个,是将一个3D的问题简化为2D问题,裁减计算量,需要留神的是,在Fluent中,回转轴务必是x轴。 5操作工况参数(Operating Conditions) 操作压力的介绍 关于参考压力的设定,首先需了解有关压力的一些定义。ANSYS FLUENT中有以下几个压力,即Static Pressure(静压)、Dynamic Pressure(动压)与Total Pressure(总压);Absolute Pressure(十足压力)、Relative Pressure(参考压
7、力)与Operating Pressure(操作压力)。 这些压力间的关系为,Total Pressure(总压)=Static Pressure(静压)+Dynamic Pressure(动压);Absolute Pressure(十足压力)=Operating Pressure(操作压力)+Gauge Pressure(表压)。 其中,静压、动压和总压是流体力学中关于压力的概念。静压是测量到的压力,动压是有关速度动能的压力,是滚动速度能量的表达。 而十足压力、操作压力和表压是FLUENT引入的压力参考量,在ANSYS FLUENT中,全体设定的压力都默认为表压。这是考虑到计算精度的问题。
8、操作压力的设定 设定操作压力时需要留神的事项如下: ? 对于不成压缩梦想气体的滚动,操作压力的设定直接影响流体密度的计算,由于对 于梦想气体而言,滚动的密度由梦想气体方程获得,梦想气体方程中的压力为操作压力。 ? 对于低马赫数的可压缩滚动而言,相比十足静压,总压降是很小的,因此其计算精 度很轻易受到数值截断误差的影响。需要采取措施来制止此误差的形成,ANSYS FLUENT通过采用表压(由十足压力减去操作压力)的形式来制止截断误差的形成,操作压力一般等于流场中的平均总压。 ? 对于高马赫数可压缩滚动的求解而言,由于此时的压力比低马赫可压缩滚动的大得 多,所以求解过程中的截断误差的影响不大,可以
9、不设定表压。由于ANSYS FLUENT中全体需输入的压力都为表压,因此此时可以将操作压力设定为0(这样可以最小化由于压力脉动而引起的误差),使表压与十足压力相等。 ? 假设密度设定为常数或者其值由通过温度变化的函数获得,操作压力并没有在计算 密度的过程中被使用。 ? 默认的操作压力为101325Pa。 操作压力的设定主要基于两点考虑,一是滚动马赫数的大小,二是密度计算方法。 表格 1 操作压力的推举设置 密度关系式 梦想气体定律 梦想气体定律 关于温度的函数 常数 不成压缩的梦想气体 大于0.1 小于0.1 不成压缩 不成压缩 不成压缩 马赫数 操作压力 0或约等于流场的平均压力 约等于流场
10、的平均压力 不使用 不使用 约等于流场的平均压力 关于参考压力位置的设定 对于不涉及任何压力边界条件的不成压缩滚动,ANSYS FLUENT在每次迭代后要调整表压值。这个过程通过使用参考压力位置处(或该位置邻近)节点的压力完成。因此,参考压力位置处的表压应一向为0。假设使用了压力边界条件,那么不会使用到上述关系,因此参考压力位置不被使用。 参考压力位置默认为等于或接近(0,0,0)的节点中心位置。实际计算中可能需要设置参考压力位置到十足静压已知的位置处。在Operating Conditions对话框中的Reference Pressure Location选项组中设置新的参考压力位置的x,y
11、,z的坐标即可。 假设要考虑某一方向的加速度,如重力,可以勾选Gravity复选框。 对于VOF计算,应选中择Specified Operating Density,并且在Operating Density 下为最轻相设置密度。这样做摈弃了水力静压的积累,提高了round-off精度为动量平衡。同样需要开启 Implicit Body Force,片面平衡压力梯度和动量方程中体积力,提高解的收敛性。 Reference Pressure Location(参考压强位置)应是位于流体永远是100%的某一相(空气)的区域,光滑和快速收敛是其根本条件。 二、求解模型的设定 1滚动模型的设置 无粘模型
12、 梦想流体是一种设想的没有粘性的流体,在滚动时各层之间没有相互作用的切应力,即没有内摩擦力。特别明显,梦想流体对于切向变形没有任何抗拒才能。理应强调指出,真正的梦想流体在客观实际中是不存在的,它只是实际流体在某些条件下的一种近似模型。 在Inviscid滚动模型应用方面,无粘滚动疏忽了粘性对滚动的影响,这对高雷诺数的滚动是适合的,由于高雷诺数滚动惯性力的作用远大于粘性力的作用,粘性力可以疏忽,所以可以将其考虑成无粘滚动。无粘滚动的求解更快,其激波在某些值上预料的偏高。无粘滚动能对滚动状态和激波位置举行快速预料。 马赫数与激波 马赫数的定义是 M?v a它表示流体的滚动速度与当地声速之比,是一个
13、无量纲的参量。对应于M?1,M?1和M?1这三种处境的滚动分别称为亚声速流、声速流和超声速流。当马赫数很小时,速度的相对变化只能引起很小的密度相对变化,但当马赫数很大时,那么将引起较大的密度相对变化,这也说领略马赫数是流体压缩性的一个表征。 当飞机、炮弹和火箭以超音速飞行时,或者发生强爆炸、强爆震时,气流受到急剧的压缩,压强和密度突然显著增加,这时所产生的压强扰度将比声速大得多的速度传播,波阵面所到之处气流的各种参数都将发生显著变化,参数突跃。这样一个强休止面叫做激波阵面。 渐缩渐扩管的滚动是计算流体力学模拟的经典问题之一。在这类滚动中,激波的展现是滚动中可压缩效应的表达。精确的激波模拟是CF
14、D研究的热点之一。为了更好抓获压力梯度,需要采用较细的网格并结合适合的数值模拟和格式。好多实际模拟中,局部网格的自适应会很有扶助。 层流模型 滚动有层流和湍流之分,判断湍流的标准可以参考2,这里写出内流的判断标准: Re?UD?2300 对于内流而言,一般大多数滚动都是湍流,一般不使用湍流模型。而对一些外流而言(如外掠平板或是外掠障碍物),那么很有可能是层流运动。 湍流模型的评价与选择 a. k?湍流模型 这里我们使用的湍流模型是Standard k?模型,这种模型应用较多,计算量适中,有较多数据积累和对比高的精度,对于曲率较大和压力梯度较强等繁杂滚动模拟效果欠佳。一般工程计算都使用该模型,其
15、收敛性和计算精度能得志一般的工程计算要求,但模拟旋流和绕流时有缺陷。 壁面函数的选择 对于有壁面的滚动,当主流为充分进展湍流时,根据离壁面法线距离不同,可将滚动划分为壁面区(或称内区、近壁区)和核心区(或称外区)。 核心区是完全湍流区,为充分进展的湍流。 在壁面区,由于有壁面的影响,滚动与核心区不同。壁面区可分为3个子层:粘性底层、过渡层和对数率层。 粘性底层是一个紧贴壁面的极薄层,在动量、热量和质量的交换过程中粘性力起主要作用,而湍流切应力可以疏忽,因此滚动几乎可以看成层流滚动,且在平行于壁面方向上的速度呈线性分布。 过渡层处于粘性底层之外,在此层中,粘性力和湍流切应力的作用相当,滚动状况较为繁杂,很难用公式或定律表述。实际工程计算中由于过渡层厚度微小,可不考虑此层,直接以对数率层的方法处理。 对数率层处于近壁区的最外层,粘性力的影响不明显,湍流切应力占主要地位,