求解器设置IntroductoryFLUENTTraining

上传人:go****e 文档编号:137233798 上传时间:2020-07-06 格式:PPT 页数:46 大小:964.01KB
返回 下载 相关 举报
求解器设置IntroductoryFLUENTTraining_第1页
第1页 / 共46页
求解器设置IntroductoryFLUENTTraining_第2页
第2页 / 共46页
求解器设置IntroductoryFLUENTTraining_第3页
第3页 / 共46页
求解器设置IntroductoryFLUENTTraining_第4页
第4页 / 共46页
求解器设置IntroductoryFLUENTTraining_第5页
第5页 / 共46页
点击查看更多>>
资源描述

《求解器设置IntroductoryFLUENTTraining》由会员分享,可在线阅读,更多相关《求解器设置IntroductoryFLUENTTraining(46页珍藏版)》请在金锄头文库上搜索。

1、求解器的设置,Introductory FLUENT Training,大纲,使用求解器 设置求解器参数 收敛 定义 监视 稳定性 加速收敛 准确度 网格独立性 网格自适应 非定常流动模型 非定常流动问题设置 非定常流动模型选项 摘要 附录,大纲,使用求解器(大致求解过程) 设置求解器参数 收敛 定义 监视 稳定性 加速收敛 准确度 网格独立性 网格自适应 非定常流动模型 非定常流动问题设置 非定常流动模型选项 摘要 附录,否,设置求解参数,初始化,激活感兴趣的监视,检查收敛,检查准确性,停,求解过程概览,求解参数 选择求解器 离散化方案 初始化 收敛 监视收敛 稳定性 设置低松弛 设置柯朗数

2、 加速收敛 精确度 网格独立性 自适应,是,是,否,可用的求解器,在FLUENT中有两个可用的求解器 基于压力的求解器 基于密度的求解器 pressure-based求解器把动量和压力(或者压力修正)作为主要变量 压力速度耦合算法则起源于从连续性方程出发推导一个压力方程 基于压力的求解器中两个运算法则可用 分离求解器按顺序求解算压力修正和动量方程 压力耦合求解器同时解算压力和动量方程,Segregated,PBCS,Solve Turbulence Equation(s),Solve Species,Solve Energy,DBCS,Solve Other Transport Equatio

3、ns as required,可用求解器,基于密度的耦合求解器以矢量形式同时求解连续性方程,动量方程,能量方程,如果需要还可求解组分方程。压力由状态方程得到。其他的标量方程用分离求解器求解。 基于密度的求解器既能用于显式也能用于隐式求解方式 隐式使用高斯/赛德尔对称块方法来解出变量 显式使用龙格库塔显式时间积分格式 注意:基于压力的求解器只有隐式格式,选择一个求解器,pressure-based求解器在从低速不可压流到高速可压流这样大范围的流动体系中都适用 需要较小内存 在求解过程中弹性较大 基于压力的耦合求解器对大多数的单相流动是适用的,而且它基于压力求解器效果好。对多相流,周期质量流和NI

4、TA型不适用。 比分离式求解器要多用1.52倍内存。 (DBCS)密度耦合解算器一般适用于有强烈耦合或者在密度,能量,动量或组分之间关联性较强的情况 举例:伴随燃烧,高超音速流,相互冲击的高速可压流 隐式解法通常更受喜欢,因为显式方式在时间步长上有严格限制 显式方式用在流动的特征时间比例与声学时间比例处于相同量级的情况下(比如高马赫数冲击波传播),离散化(插值方法),场变量(存储于单元体中心)必须内插于控制体的面上 对流形式的内插方案 First-Order Upwind 一阶迎风最容易收敛,只有一阶精度 Power Law 幂率定理当雷诺数小于5(低雷诺数)时比一阶精度更高 Second-O

5、rder Upwind 二阶迎风对一阶精度来说使用更小的截断误差,适用三角/四面体网格或者当流动与网格不在一条直线上;收敛可能变慢 Monotone Upstream-Centered Schemes for Conservation Laws (MUSCL) 当地3阶对流离散化格式,用于非结构化网格;在预测二次流,漩涡,力等等时更精确 Quadratic Upwind Interpolation (QUICK) 二次迎风插值法,插值方法(梯度),对于更高阶数的离散化方案,求解扩散通量和速度派生需要先解出变量的梯度。 单元体中心处结果变量的梯度可以用三种方法得到 Green-Gauss Cel

6、l-Based 格林高斯基于单元体默认方法;解有伪扩散(求解域的拖尾现象) Green-Gauss Node-Based 格林高斯基于节点更精确;最小化伪扩散;推荐用在三角网格上 Least-Squares Cell-Based 基于单元体的最小二乘法推荐用于多面体网格;与上一种有相同精度和属性 在面上的结果变量梯度用多维泰勒级数展开,对面压力的插值方法,下面是FLUENT中可用的使用分离求解器的用来计算单元体面上压力的插值格式: Standard 标准缺省格式;对大表面边界层附近的曲面法线压力梯度流动求解精度会降低(但是不能够用于流动中压力急剧变化的地方应该用PRESTO!格式 代替) PR

7、ESTO! 用于高回旋流,压力急剧变化流(多孔介质,涡扇模型,等等),或者剧烈弯曲的区域 Linear 线性的当其它选项导致收敛困难或者出现非物理现象时使用 Second-Order 二阶用于可压流;不能用于多孔介质,跳跃,风扇,等等。或者VOF/混合多相流。 Body Force Weighted 体积力当体积力很大时,比如高雷诺数自然对流或者高回旋流。,压力速度耦合,当使用基于压力的求解器时,压力速度耦合是指在离散连续性方程时,利用连续性方程和动量方程联合推导出一个压力方程(或者压力修正)的数值运算法则 在FLUENT中有下面四种耦合方法 Semi-Implicit Method for

8、Pressure-Linked Equations (SIMPLE) 半隐式连接压力方程方法 缺省格式 SIMPLE-Consistent (SIMPLEC) 允许对简单问题快速收敛(比如,没有物理上可用模型的层流) Pressure-Implicit with Splitting of Operators (PISO) 对非定常流动问题或者包含比平均网格倾斜度更高的网格适用 Fractional Step Method (FSM) 对非定常流的分步方法 用于NITA格式;与PISO特征相同,初始化,迭代计算之前需要给所有变量初始化 和实际情况近似的假设可改善解的稳定性和加速收敛 对有些问题,

9、一个好的初始假设是必需的。 举例:高温区域引起化学反应 修补值用于特定区域的独立变量 自由射流(高速射流) 燃烧问题(高温区域来初始化反应) 单元寄存器(在自适应面板上给单元体作记号)能够用于补缀不同领域的不同地方,全部多重网格初始化,全部多重网格初始化能够用来创建一个更好的流场初始化 TUI command: TUI 命令: /solve/init/fmg-initialization 在计算上又省又快。在粗网格上用一阶精度的欧拉方程计算 在基于压力和密度的求解器中都能使用,但是只能用于定常状态。 在把解转到实际网格上之前, FMG 使用完全近似存储(FAS)多重网格方法用在一系列粗网格上来

10、解决流动问题。 设置能够用TUI命令存取/solve/init/set-fmg-initialization FMG初始化方法对很多复杂问题。包括在一个大的计算域上存在很大压力和速度梯度的复杂流动问题都比较有用(比如:旋转喷灌器,扩展螺旋管),检查Case,Case Check是FLUENT中的一种非常有用的功能,用于寻找常规的设置错误和在选择case参数和模型时提供指导 使用法则和最优方法 Case 检查主要依次检查以下内容 网格 模型选择 边界条件 材料属性 求解器设置 列表选项部分包含推荐 自动建议功能将自动改变错误设置 手动推荐:手动改变,Solve,Case Check,大纲,使用求

11、解器(大致求解过程) 设置求解器参数 收敛 定义 监视 稳定性 加速收敛 准确度 网格独立性 网格自适应 非定常流动模型 非定常流动问题设置 非定常流动模型选项 摘要 附录,收敛,收敛必须满足下面条件: 所有离散化的守恒方程(动量,能量,等等)在所有单元体中都达到设定的残差或者解算在之后的迭代中不再改变 全部的质量,动量,能量,和标量平衡都要达到 使用残差记录来监视收敛 一般来说,残差减少了三个数量级就意味着定性的收敛。在这个时候,主要的流动特征就应该已经呈现出来了。 能量残差大小必须减少到 10-6 (对基于压力的求解器而言) 组分残差大小需要减少到10-5 以达到介质平衡。 监视定量收敛:

12、 监视其它相关关键变量/物理量来确认 确保全部质量/热量/物质守恒得到满足,收敛监视:残差,残差图显示直到残差值达到了指定的容差范围,All equations converged.,10-3,10-6,Solve,Monitors,Residual,收敛监视:力/面,除了残差之外,你还能够监视: 升力,阻力或力矩 相关变量或函数(比如,面积分)在边界上或者定义的面上,Solve,Monitors,Force,Solve,Monitors,Surface,校核通量守恒,除了监视残差和变量历史之外,你还能检查整体的热量和质量平衡 净不平衡量应该小于区域边界上最小通量的1%,Report,Flux

13、es,固定收敛容差,如果解的监视器指示出解是收敛的,但是解依然跳动或者在质量/热量上不平衡,这也清楚地表明解还没有收敛 在这种情况下,你需要: 减少收敛性判据的值或者在残差监视面板关掉残差监视面板 继续迭代直到解收敛 在收敛性判据下选择none将会指示FLUENT不给任何方程检查收敛。,Solve,Monitors,Solve,Iterate,Residual,收敛困难,收敛性差,数值计算稳定性差,可能是因为提出的问题本身就存在问题或是由于网格的质量差,或者不适当的求解器设置出现的。 出现增长的(分叉的)或者“停滞的”残差 分离的残差意味着在守恒方程中不平衡一直在增长 不收敛的结果是错误的,C

14、ontinuity equation convergence trouble affects convergence of all equations.,解决方案 确保问题建立是正确的 用一阶离散化格式计算一个初始解 对有收敛问题的方程减少低松弛因子(基于压力的求解器) 减少柯朗数(基于密度的求解器) 重画网格和修改有大纵横比和大斜度的单元体,修饰欠松弛因子,内嵌欠松弛因子 用来稳定压力求解器的迭代过程 使用默认欠松弛因子来开始计算 减少动量的欠松弛因子经常能够帮助收敛 默认设置是对较大范围的问题适用的,如果需要的话,你可以减少这个值 适当的设置最好从试验中得到,对基于密度的求解器来说,欠松弛

15、因子对于外部方程连接的设置修正用在基于压力的求解器,Solve,Controls,Solution,修正柯朗数,即使是定常问题也包括衰减项 柯朗数定义时间步 对基于密度的显式求解器 稳定性约束给柯朗数一个最大的限制 不能大于2(默认值为1) 当收敛困难时减小柯朗数 基于密度的隐式算法 柯朗数不受稳定性约束限制 默认值是5,Solve,Controls,Solution,Mesh size,Appropriate velocity scale,加速收敛,收敛可以由以下方式加速: 提供更好的初始条件 从一个已有的解开始(如果需要的话可以使用文件或者插值) 逐渐增加的欠松弛因子或者柯朗数 过高的值会

16、导致不稳定或者收敛问题 建议在继续迭代之前保存case和data文件 控制多重网格求解器设定(但是默认设置提供了一个精确的网格设置而且一般不需要改动),从已有的解开始,当改变用于问题说明时已有的解可以作为初始条件 使用解的插值来初始化一个运行(当粗糙网格解是可用的时候对开始一个细网格case特别有用) 一旦这个解初始化,附加的迭代经常使用现在的数据作为起始点 对一些实际问题如何提供初始解的一些建议,大纲,使用求解器(大致求解过程) 设置求解器参数 收敛 定义 监视 稳定性 加速收敛 准确度 网格独立性 网格自适应 非定常流动模型 非定常流动问题设置 非定常流动模型选项 摘要 附录,解的精确性,收敛解不一定是正确解 经常使用测量数据,物理定理来检查和评价结果 对最终结果使用二阶迎风离散化格式 确保结果是独立于网格的 使用自适应来修正网格或者创建增加的网格数使结果独立于网格If 如果流动特征看起来不合理: 重新考虑物理模型和边界条件 检查网格质量而且可能要重画 重新考虑边界位置的选择(或者计算域的位置):计算域的选择不适当(特别是出流边界)会显著影响解的精确性,网

展开阅读全文
相关资源
相关搜索

当前位置:首页 > 幼儿/小学教育 > 其它小学文档

电脑版 |金锄头文库版权所有
经营许可证:蜀ICP备13022795号 | 川公网安备 51140202000112号