ansys12-FLUENT04求解器设置

检查Case
Case Check 功能发现一些常 见的错误设置和不一致性
– 提供选择参数和模型的指导
用于可以直接应用或忽略这 些建议
A Pera Global Company © PERA China
收敛性
计算收敛时应该满足:
– 所有离散的守恒方程（动量、能量等）在所有的单元中满足指定 的误差或者结果随计算不再改变 – 全局的质量、动量、能量和标量达到平衡
A Pera Global Company © PERA China
求解器选择
密度基耦合求解器
– 以矢量方式求解连续性方程、动 量方程、能量方程和组分方程
– 通过状态方程得到压力
– 其他标量方程按照分离方式求解
DBCS 可以显式或隐式方式求解
– 隐式 – 使用高斯赛德尔方法求解 所有变量 – 显式: 用多步龙格库塔显式时间积 分法。
No
Stop
A Pera Global Company © PERA China
求解器选择
FLUENT中有两种求解器 – 压 力基和密度基。 压力基求解器以动量和压力为 基本变量 – 通过连续性方程导出压力和 速度的耦合算法 压力基求解器有两种算法 – 分离求解器 – 压力修正和动 量方程顺序求解。 – 耦合求解器 (PBCS) –压力和 动量方程同时求解
ANSYS FLUENT 培训教材 第四节：求解器设置
安世亚太科技（北京）有限公司
A Pera Global Company © PERA China
概要
使用求解器（求解过程概览） – 设置求解器参数 – 收敛 • 定义 • 监测
• 稳定性
• 加速收敛 – 精度 • 网格无关性 • 网格自适应 – 非稳态流模拟（后续章节中介绍） • 非稳态流问题设置
A Pera Global Company © PERA China
离散化（插值方法）
存储在单元中心的流场变量必须插值到控制体面上
对流项的插值方法有:
– First-Order Upwind – 易收敛，一阶精度。
– Power Law –对低雷诺数流动 （ Recell < 5 ）比一阶格式更精确
– 表现为残差曲线上扬（发散）或不下降 – 发散意味守恒方程的不平衡增加 – 没收敛的结果会误导使用者
解决方法 – 确保问题是物理合理的 – 用一阶离散格式计算一个初场 – 对压力基求解器，减少发散方程的 松弛因子 – 对密度基求解器，减少Courant 数 – 重新生成网格或加密质量差的网格
• 注意网格自适应不能提高扭曲度大 的网格质量
单元中心的变量梯度由以下三种方法得到:
– Green-Gauss Cell-Based – 可能会引起伪扩散 – Green-Gauss Node-Based – 更精确，更少伪扩散，建议对三角 形/四面体网格采用 – Least-Squares Cell-Based – 建议对多面体网格采用，精度和属 性同Node-based
Enabling pressurebased coupled solver (PBCS)
A Pera Global Company © PERA China
如何选择求解器
压力基求解器应用范围覆盖从低压不可压缩流到高速压缩流 – 需要的内存少 – 求解过程灵活 压力基耦合求解器 (PBCS) 适用于大多数单相流，比分离求解器性能更 好 – 不能用于多相流（欧拉）、周期质量流和 NITA – 比分离求解器多用1.5–2倍内存 密度基耦合求解器 (DBCS)适用于密度、能量、动量、组分间强耦合的 现象 – 例如: 伴有燃烧的高速可压缩流动，超高音速流动、激波干扰 隐式方法一般优于显式，因为其对时间步有严格的限制 显式方法一般用于流动时间尺度和声学时间尺度相当的情况（如高马 赫激波的传播）
– 压力基求解器: /solve/init/fmg-initialization – 密度基求解器: 当选择密度基求解器后在 GUI里可见
FMG 在粗网格上用多重网格求解
– 通过 TUI 命令来设置
/solve/init/set-fmg-initialization
A Pera Global Company © PERA China
使用残差历史曲线来监测收敛：
– 一般地，残差下降三个量级表示至少达到定性的收敛，流场的主 要特征已经形成。 – 压力基求解器的能量残差应下降到10-6 – 组分残差应下降到10-5
监测定量的收敛:
– 监测其他关键的物理量 – 确保全局的质量、能量、组分守恒。
A Pera Global Company © PERA China
– 缺省值对大多数问题都适用， 需要时你可以改变这些值
– 合适的设置最好通过经验获得
对密度基求解器，对耦合方程组 外的方程，松弛因子同样有用
A Pera Global Company © PERA China
修改 Courant 数
对密度基求解器，即使稳态问题， 也存在瞬态项
– 用Courant 定义时间步长
– 设置更好的初场
• 从前次的计算结果开始（如需要，可以使用file/interpolation）
– 渐次增加松弛因子或Courant数
• 过高的值容易引起发散 • 继续迭代是应保存case和date文件
– 控制多重网格求解器设置 （一般不推荐）
• 缺省设置一般足够稳定，不建议修改
A Pera Global Company © PERA China
FMG 初始化
Full MultiGrid (FMG) 能用来创建更好的初场。
– FMG 初始化对包括大的压力梯度和速度梯度的复杂流动有用 – 在粗级别网格上求解一阶欧拉方程 – 可用于压力基或密度基求解器，但限于稳态问题
启动 FMG 初始化
面上的梯度用多级泰勒级数展开求得
A Pera Global Company © PERA China
压力的插值方法
使用分离算法时，计算面上压力的插值方法有:
– Standard – 默认格式，对于近边界的沿面法向存在大压力梯度流 动，精度下降（如果存在压力突变，建议改用 PRESTO! ） – PRESTO! – 用于高度旋流，包括压力梯度突变（多孔介质，风 扇模型等）或者计算域存在大曲率的面 – Linear – 当其他格式导致收敛问题或非物理解时使用 – Second-Order – 用于压缩流，不适用多孔介质、风扇、压力突 变以及VOF/Mixture 多相流 – Body Force Weighted – 用于大体积力的情况，如高瑞利数自然 对流或高旋流
Calculate a solution
Modify solution parameters or grid
– 稳定性
• 设置松弛因子 • 设置 Courant number
Yes
Check for convergence No
– 加速收敛
精度 – 网格无关性 – 自适应网格
Yes
Check for accuracy
• 和 NITA 合用，类似 PISO.
A Pera Global Company © PERA China
初始化
FLUENT 要求所有的求解变 量有初始值 – 更真实的初值能提高收敛稳 定性，加速收敛过程. – 有些情况需要一个好的初值
在特定区域对特定变量单独 赋值 – 自由射流（喷射区高速） – 燃烧问题 （高温激活反应） – 单元标注（自适应）
Pressure-Based (segregated)
Solve U-Momentum Solve V-Momentum Solve W-Momentum Solve Mass Continuity; Update Velocity
Pressure-Based Density-Based (coupled) (coupled)
监测收敛-残差
残差图显示何时收敛达到指定标准
All equations converged. 10-3
10-6
A Pera Global Company © PERA China
监测收敛-力和面上的变量
除了残差外，也可以监测升力、阻力和 力矩系数 边界或其他定义的面上的导出变量或函 数（如面积分）
– 减小残差标准或关闭监测残差的窗口 – 继续迭代直至计算收敛
在Convergence Criterion
窗口选择 None关闭监测
残差的窗口
A Pera Global Company © PERA China
收敛遇到的困难
对一些病态问题，差质量的网格或不合适的求解设置，都 可能出现数值不稳定性
– Quadratic Upwind Interpolation (QUICK) – 适用于四边形/六面体以及 混合网格，对旋转流动有用，在均匀网格上能达到三阶精度
A Pera Global Company © PERA China
插值方法（梯度）
为了得到扩散通量、速度导数，以及高阶离散格式，都需要求解 变量的梯度
Solve Mass & Momentum
Solve Mass, Momentum, Energy, Species
Solve Energy Solve Species Solve Turbulence Equation(s) Solve Other Transport Equations as required
对显式求解器:
– 稳定性约束限制了最大Courant 数
• 不能超过 2 （缺省为1） • 有收敛困难时减少 Courant 数
对隐式求解器:
– Courant 没有稳定性约束限制
• 缺省值为 5.
A Pera Global Company © PERA China