第七章-FLOW-3D-数值选项

第七章、FLOW-3D

数值选项FLOW-3D®v9.3数值选项模拟的主要目标是提供准确的结果。所有数值选项都默认设置。—大多数模拟用默认设置就可。—大多数求解器使用默认显式的方法。数值选项包括。压力解算器。时间步长控制。隐式求解传热，隐式表面张力等。这章节的主要内容是：什么时候用户需要选择其他数值选项。这些选项在哪里选择。不可压缩流动的控制方程要求：—动量方程（1，2或3的方向）。—质量守恒。可选：流体分数（自由表面）能源平流（传热）密度平流（可变密度）湍流颗粒派生数量：—应变率，固相分数，非牛顿粘度。动量方程U=(u,v,w)—流速;P—压强;G

—重力和非惯性体加速;t

—粘性应力张量;KU

—拖动（多孔挡板，障碍，两相区);RSOR

U/r—大量注射液致零速度加速;F

—其他力：表面张力，电，质量/动力源，颗粒，用户定义的力.不可压缩，非牛顿/牛顿流动性：zzzFwRSORKwρGzPzwwywvxwutw---D-+-=þýüîíì+++rt¶¶r¶¶¶¶¶¶¶¶11质量守恒定律方程RSOR—质量源/汇RDIF

—密度扩散不可压缩流动（用质量源）变密度或可压缩流动（用质量源和漫射）热能传送方程fs —固相分数（是T和溶质浓度的功能）;L —潜热;k —导热系数;C(T) —比热;h —液体/墙体传热系数;Twall —壁温;RISOR —能量源/汇;RIDIF

—湍流扩散LfsdTTCIT

)1(

)(

-+=ò流体的体积对流方程运动学上的运动方程：VOF功能运动是根据在流动性里的速度场。F=1.0-fluidF=0.0-voidUFDIF-流体分数扩散。FSOR-流体体积源/汇。控制体积解决方法FLOW-3D解决在交错网格有限差分方程。在控制体积单元里的位置和条件(i,j,k)压力，温度，标量储存在单元中心。强调速度和计算是在单元的表面。计算单元编号条件第一次真正的单元是i=2,j=2,k=2第一次真正的单元是i=3,j=3,k=3单元是自动被加在网格边执行边界条件：1层为大多数的边界条件(symmetry,wall,velocity,outflow,gridoverlay)。压力、周期和跨层块边界需要两层。只有真正的单元（见在白色下面）都显示在图里。FLOW-3D数值选项时间步距控制压力解算器显式或隐式求解流体界面平流动量平流速度场选项时间步长控制时间步距的主要用途（称为dt和delt）：模拟时间步距向前（每个称为循环步）保持准确和稳定的模拟时间步长受到各种因素的控制：Dt在以前的时间步距。稳定极限(Dt不能大于稳定性)。一些关于前一个周期的压力迭代（当AUTOT=1）。用户定义最大的Dt(DTMAX)FLOW-3D®控制整个模拟时间步长：掠过所有的计算单元.在每个单元格相关稳定性准则计算.确定所有单元的最低标准.最大稳定时间步长=

Dt=全球稳定的最低标准.控制时间步长在模拟完成时间上的时间步长控制通过FLOW-3D的基础来选定。时间步长控制-'ModelSetup-Numerics标签'.初始时间步长----默认=1.0E-6×完成时间;可以有效地启动模拟.最小时间步长----默认=1.0E-06x初始的时间步长;防止模拟失控（非常小的Dt）.最大时间步长----默认值是稳定极限。

在FLOW-3D®时间步长控制选择的间步长有三个因素影响：稳定精确度效率稳定----最明确的方法需要对时间步长限制所谓'稳定性准则'。时间步长始终是低于最低标准的稳定性。稳定的标准自动确定在每个周期的开始.精确度----每个周期时间步长增加不能超过5％。效率-平衡时间步长vs收敛。较大的时间步长需要更多的压力/粘滞/热量收敛。FLOW-3D可减少时间步长（稳定以下）来实现更快的收敛。稳定性和收敛固定大小时间步长的稳定性标准稳定性标准确定在激活物理模型的每个周期开始。两个字母的代码表明适用于该物理模型稳定性的标准。求解摘要在“长期印刷”间隔列出了所有稳定性标准。求解消息文件在“短期印刷”的间隔列出控制稳定标准.隐式方法加快模拟？x-explicitconvectivetime-steplimit=1.51646E-06atcell(10,57,6)inmeshblock2y-explicitconvectivetime-steplimit=2.46870E-06atcell(11,57,7)inmeshblock2z-explicitconvectivetime-steplimit=1.15658E-06atcell(10,57,6)inmeshblock2x-implicitconvectivetime-steplimit=6.06859E-06atcell(26,56,14)inmeshblock2y-implicitconvectivetime-steplimit=6.87478E-06atcell(38,21,44)inmeshblock1z-implicitconvectivetime-steplimit=7.64813E-06atcell(38,30,43)inmeshblock1freesurfacetime-steplimit= 1.37319E-03atcell(16,6,7)inmeshblock2thermalconductiontime-steplimit=2.15091E-04atcell(18,17,8)inmeshblock2wallheattransfertime-steplimit=1.55914E-04atcell(7,56,9)inmeshblock2典型的信息为了稳定性极限的时间步长.显式求解的诊断综述搜索“短的印刷时间步长的稳定性极限”什么是物理过程的模拟限制？

—所指出的最小极限.如果删除最小的限度（隐式）下一个最大的？压力解算器在FLOW-3D有三个压力求解器可用到.SOR–从头到尾连续松驰.

–相对良好的连接几何.

–相当统一的网格单元.

–短小的模拟.

–SOR收敛应该在10〜15左右.ADI–交替方向隐式.

–大单元的aspectratios（〜10:1）

–在SMP平行面没有良好的比例.GMRES–不能适应特殊环境的最小残留的方法.

–处理复杂，适应计算良好的领域.

–大单元aspectratios的收敛.

–内存要求大-当使用SOR和ADI不衔接.

–在SMP平行面比例非常大.

–GMRES方法收敛次数应为2〜7.压力迭代收敛输入小于1的正值,迭代可以收紧。不建议使用大于1的值输入正值,迭代也可以收紧稳态流动以稳态流来加速方式可以通过设置一个很大的数字.如输入1.0e6这个选项只能用在没有传热的ConfinedFlows(onefluid,nosharpinterface)显式/隐式求解选项显式和隐式求解器可用于粘性应力，热传导，表面张力，弹性应力，泡沫压力，GMO/流动偶合,和平流动力。通常默认显式方法。显式方法：

–精确的.–当几何和流动是相对简单时液体快速流动的

–实现简单数字.

–需要时间步长限制（稳定性标准）以确保稳定.隐式的方法：–允许比显式的方法更大的时间步长.–通常需要收敛求解。–可能会或可能不会产生精确的解决方案。界面平流选项通过FLOW-3D模拟类型的基础上选择VOF平流方法。Automatic–通过FLOW-3D模拟类型的基础上自动选择的方法。默认的选择是通常良好的工作高级平流方法可供选择。当对流体积时使用误差较大（少数的％）。F-PackingF-packing算法用于在具有重大自由表面崩溃流动来消除小洞或“foaming”的。该算法通过建立在流体单元内部小的不确定的分歧，其中哪个流体分数小于0.99。这是只使用在一个自由表面流动的液体。默认值为1.0,当你输入0时Packing就没发生IFVOF=4＆6建议具有较强的流动循环VOF平流方法对角流动的比较IFVOF=4vs5IFVOF=4IFVOF=5在高压压铸VOF平流方法的比较填充阶段:–填充模拟.–填充时间〜0.048s.–高湍流.–模拟目标：找到被困气体.VOF4vs5比较.IFVOF=4方法液量保持良好，预计被困气体。IFVOF=5的方法显示大型对流流动错误。VOF平流方法的总结方法优点缺点IFVOF=4简单地实现数字.包含的所有信息：F,NF,AF和VF。良好的单向准确性.有许多特殊情况的逻辑，例如在出口边界.在2D和3D流动量守恒和界面跟踪减少准确性。IFVOF=5在跟踪复杂流动急剧界面里的良好性。改进表面张力行为.没有圆柱坐标的支持。在FAVOR是大约的部分块单元准确性减少。在旋涡流动有大量的错误结果(流动损失）。IFVOF=6在跟踪复杂流动急剧界面里的良好性。较好的均衡表面张力行为.圆柱坐标支持9.3版本进度。在部分块单元IFVOF=4。非保守的制定，体积误差不代表大量的损失/增加.动量平流4个动量（速度）平流选项:–IORDER=11st动量平流.–

IORDER=22nd为中心差动.–IORDER=32nd为逆风单一性保留.

–IORDER=43rd为逆风单一性保留.大部分模拟默认1st:–计算效率–预计大部分流动特性建议用高阶动力方法有如下:旋转流动.二次流动特征.高阶动力的应用涡流/旋转流动的模拟涡流/旋转流动的特点:旋转的主要流程.二次流动的特点往往是重要的.旋涡流动的例子:沉淀池.涡流喷嘴.SedimentCollectionTankSimulation

GeneralConcept淤泥流入出水固体缺陷主要旋流两面中心+二次流动.利用“茶杯效应”.较重的固体聚集在中心.目标：模拟主旋转流动和二次流动。只有流动的计算-没有固体.在圆槽总共速度的实验结果上图显示槽附近收集数字代表一小部分的入口速度预计沿中心线强涡旋在涡流1.2倍速度大于入口速度涡流模拟结果IORDER=1Initallysolid-bodyrotationNovortexingseennearaxis1storderadvectiondissapatesvortexTryhigherordermomentumadvection默认的1st平流IORDER=1，计算效率是很多问题，但太多的对角线速度分量扩散而准确预测旋流。模拟结果IORDER=2速度是不是合理－不稳定．IORDER=2不建议自由面模拟．不稳定由于IORDER=2是已知的强循环与自由表