第七章flow-3d数值选项.ppt

第七章、FLOW-3D数值选项 FLOW-3D® v9.3 数值选项 模拟的主要目标是提供准确的结果。 所有数值选项都默认设置。 —大多数模拟用默认设置就可。 —大多数求解器使用默认显式的方法。 数值选项包括。 压力解算器。 时间步长控制。 隐式求解传热，隐式表面张力等。 这章节的主要内容是： 什么时候用户需要选择其他数值选项。 这些选项在哪里选择。 不可压缩流动的控制方程 要求： —动量方程（1，2或3的方向）。 —质量守恒。 可选： 流体分数（自由表面） 能源平流（传热） 密度平流（可变密度） 湍流 颗粒 派生数量： —应变率，固相分数，非牛顿粘度。 动量方程 U = (u, v, w) —流速; P —压强; G —重力和非惯性体加速; t —粘性应力张量; K U —拖动（多孔挡板，障碍，两相区); RSOR U/r —大量注射液致零速度加速; F —其他力：表面张力，电，质量/动力 源，颗粒，用户定义的力. 不可压缩，非牛顿/牛顿流动性： 质量守恒定律方程 RSOR —质量源/汇 RDIF —密度扩散 不可压缩流动（用质量源） 变密度或可压缩流动（用质量源和漫射） 热能传送方程 fs —固相分数（是T和溶质浓度的功能）; L —潜热; k —导热系数; C(T) —比热; h —液体/墙体传热系数; Twall —壁温; RISOR —能量源/汇; RIDIF —湍流扩散 流体的体积对流方程 运动学上的运动方程：VOF功能运动是根据在流动性里的速度场。 FDIF -流体分数扩散。 FSOR -流体体积源/汇。 控制体积解决方法 FLOW-3D解决在交错网格有限差分方程。 在控制体积单元里的位置和条件(i, j, k) 压力，温度，标量储存在单元中心。 强调速度和计算是在单元的表面。 计算单元编号条件 第一次真正的单元是i =2, j =2, k =2 第一次真正的单元是i =3, j =3, k =3 单元是自动被加在网格边执行边界条件： 1层为大多数的边界条件(symmetry, wall, velocity, outflow, grid overlay)。 压力、周期和跨层块边界需要两层。 只有真正的单元（见在白色下面）都显示在图里。 FLOW-3D数值选项 时间步距 控制 压力解算器 显式或隐式求解 流体界面平流 动量平流 速度场选项 时间步长控制 时间步距的主要用途（称为dt和delt）： 模拟时间步距向前（每个称为循环步） 保持准确和稳定的模拟 时间步长受到各种因素的控制： Dt 在以前的时间步距。 稳定极限(Dt不能大于稳定性)。 一些关于前一个周期的压力迭代（当AUTOT = 1）。 用户定义最大的Dt (DTMAX) FLOW-3D®控制整个模拟时间步长： 掠过所有的计算单元. 在每个单元格相关稳定性准则计算. 确定所有单元的最低标准. 最大稳定时间步长= Dt =全球稳定的最低标准. 控制时间步长 在模拟完成时间上的时间步长控制通过FLOW-3D 的基础来选定。 时间步长控制-Model Setup-Numerics标签. 初始时间步长----默认=1.0E-6×完成时间;可以有效地启动模拟. 最小时间步长----默认=1.0E-06 x初始的时间步长;防止模拟失控（非常小的Dt ）. 最大时间步长----默认值是稳定极限。 在FLOW-3D®时间步长控制 选择的间步长有三个因素影响： 稳定 精确度 效率 稳定----最明确的方法需要对时间步长限制所谓稳定性准则。时间步长始终是低于最低标准的稳定性。稳定的标准自动确定在每个周期的开始. 精确度----每个周期时间步长增加不能超过5％。 效率-平衡时间步长vs收敛。 较大的时间步长需要更多的压力/粘滞/热量收敛。FLOW-3D可减少时间步长（稳定以下）来实现更快的收敛。 稳定性和收敛 固定大小 时间步长的稳定性标准 稳定性标准确定在激活物理模型的每个周期开始。 两个字母的代码表明适用于该物理模型稳定性的标准。 求解摘要在“长期印刷”间隔列出了所有稳定性标准。 求解消息文件在“短期印刷”的间隔列出控制稳定标准. 隐式方法加快模拟？ x- explicit convective time-step limit = 1.51646E-06 at cell ( 10, 57, 6) in mesh block 2 y- explicit convective time-step l