第七章 FLOW-3D 数值选项PPT幻灯片

1,第七章、FLOW-3D数值选项,FLOW-3Dv9.3,2,数值选项,模拟的主要目标是提供准确的结果。所有数值选项都默认设置。大多数模拟用默认设置就可。大多数求解器使用默认显式的方法。数值选项包括。压力解算器。时间步长控制。隐式求解传热，隐式表面张力等。这章节的主要内容是：什么时候用户需要选择其他数值选项。这些选项在哪里选择。,3,不可压缩流动的控制方程,要求：动量方程（1，2或3的方向）。质量守恒。可选：流体分数（自由表面）能源平流（传热）密度平流（可变密度）湍流颗粒派生数量：应变率，固相分数，非牛顿粘度。,4,动量方程,U=(u,v,w)流速;P压强;G重力和非惯性体加速;t粘性应力张量;KU拖动（多孔挡板，障碍，两相区);RSORU/r大量注射液致零速度加速;F其他力：表面张力，电，质量/动力源，颗粒，用户定义的力.,不可压缩，非牛顿/牛顿流动性：,5,质量守恒定律方程,RSOR质量源/汇RDIF密度扩散,不可压缩流动（用质量源）,变密度或可压缩流动（用质量源和漫射）,6,热能传送方程,fs固相分数（是T和溶质浓度的功能）;L潜热;k导热系数;C(T)比热;h液体/墙体传热系数;Twall壁温;RISOR能量源/汇;RIDIF湍流扩散,7,流体的体积对流方程,运动学上的运动方程：VOF功能运动是根据在流动性里的速度场。,FDIF-流体分数扩散。FSOR-流体体积源/汇。,8,控制体积解决方法,FLOW-3D解决在交错网格有限差分方程。,在控制体积单元里的位置和条件(i,j,k)压力，温度，标量储存在单元中心。强调速度和计算是在单元的表面。,9,计算单元编号条件,第一次真正的单元是i=2,j=2,k=2,第一次真正的单元是i=3,j=3,k=3,单元是自动被加在网格边执行边界条件：1层为大多数的边界条件(symmetry,wall,velocity,outflow,gridoverlay)。压力、周期和跨层块边界需要两层。只有真正的单元（见在白色下面）都显示在图里。,10,FLOW-3D数值选项,时间步距控制,压力解算器,显式或隐式求解,流体界面平流,动量平流,速度场选项,11,时间步长控制,时间步距的主要用途（称为dt和delt）：模拟时间步距向前（每个称为循环步）保持准确和稳定的模拟时间步长受到各种因素的控制：Dt在以前的时间步距。稳定极限(Dt不能大于稳定性)。一些关于前一个周期的压力迭代（当AUTOT=1）。用户定义最大的Dt(DTMAX)FLOW-3D控制整个模拟时间步长：掠过所有的计算单元.在每个单元格相关稳定性准则计算.确定所有单元的最低标准.最大稳定时间步长=Dt=全球稳定的最低标准.,12,控制时间步长,在模拟完成时间上的时间步长控制通过FLOW-3D的基础来选定。时间步长控制-ModelSetup-Numerics标签.初始时间步长-默认=1.0E-6完成时间;可以有效地启动模拟.最小时间步长-默认=1.0E-06x初始的时间步长;防止模拟失控（非常小的Dt）.最大时间步长-默认值是稳定极限。,13,在FLOW-3D时间步长控制,选择的间步长有三个因素影响：稳定精确度效率稳定-最明确的方法需要对时间步长限制所谓稳定性准则。时间步长始终是低于最低标准的稳定性。稳定的标准自动确定在每个周期的开始.精确度-每个周期时间步长增加不能超过5。效率-平衡时间步长vs收敛。较大的时间步长需要更多的压力/粘滞/热量收敛。FLOW-3D可减少时间步长（稳定以下）来实现更快的收敛。,稳定性和收敛,固定大小,14,时间步长的稳定性标准,稳定性标准确定在激活物理模型的每个周期开始。两个字母的代码表明适用于该物理模型稳定性的标准。求解摘要在“长期印刷”间隔列出了所有稳定性标准。求解消息文件在“短期印刷”的间隔列出控制稳定标准.,15,隐式方法加快模拟？,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,典型的信息为了稳定性极限的时间步长.,显式求解的诊断综述搜索“短的印刷时间步长的稳定性极限”什么是物理过程的模拟限制？所指出的最小极限.如果删除最小的限度（隐式）下一个最大的？,16,压力解算器,在FLOW-3D有三个压力求解器可用到.SOR从头到尾连续松驰.相对良好的连接几何.相当统一的网格单元.短小的模拟.SOR收敛应该在1015左右.ADI交替方向隐式.大单元的aspectratios（10:1）在SMP平行面没有良好的比例.GMRES不能适应特殊环境的最小残留的方法.处理复杂，适应计算良好的领域.大单元aspectratios的收敛.内存要求大-当使用SOR和ADI不衔接.在SMP平行面比例非常大.GMRES方法收敛次数应为27.,17,压力迭代收敛,输入小于1的正值,迭代可以收紧。不建议使用大于1的值,输入正值,迭代也可以收紧,18,稳态流动,以稳态流来加速方式可以通过设置一个很大的数字.如输入1.0e6,这个选项只能用在没有传热的ConfinedFlows(onefluid,nosharpinterface),19,显式/隐式求解选项,显式和隐式求解器可用于粘性应力，热传导，表面张力，弹性应力，泡沫压力，GMO/流动偶合,和平流动力。通常默认显式方法。显式方法：精确的.当几何和流动是相对简单时液体快速流动的实现简单数字.需要时间步长限制（稳定性标准）以确保稳定.隐式的方法：允许比显式的方法更大的时间步长.通常需要收敛求解。可能会或可能不会产生精确的解决方案。,20,界面平流选项,通过FLOW-3D模拟类型的基础上选择VOF平流方法。Automatic通过FLOW-3D模拟类型的基础上自动选择的方法。默认的选择是通常良好的工作高级平流方法可供选择。当对流体积时使用误差较大（少数的）。,21,F-Packing,F-packing算法用于在具有重大自由表面崩溃流动来消除小洞或“foaming”的。该算法通过建立在流体单元内部小的不确定的分歧，其中哪个流体分数小于0.99。这是只使用在一个自由表面流动的液体。,默认值为1.0,当你输入0时Packing就没发生,22,IFVOF=46建议具有较强的流动循环,23,VOF平流方法对角流动的比较IFVOF=4vs5,IFVOF=4,IFVOF=5,24,在高压压铸VOF平流方法的比较,填充阶段:填充模拟.填充时间0.048s.高湍流.模拟目标：找到被困气体.VOF4vs5比较.IFVOF=4方法液量保持良好，预计被困气体。IFVOF=5的方法显示大型对流流动错误。,25,VOF平流方法的总结,26,动量平流,4个动量（速度）平流选项:IORDER=11st动量平流.IORDER=22nd为中心差动.IORDER=32nd为逆风单一性保留.IORDER=43rd为逆风单一性保留.大部分模拟默认1st:计算效率预计大部分流动特性建议用高阶动力方法有如下:旋转流动.二次流动特征.,27,高阶动力的应用涡流/旋转流动的模拟,涡流/旋转流动的特点:旋转的主要流程.二次流动的特点往往是重要的.旋涡流动的例子:沉淀池.涡流喷嘴.,28,SedimentCollectionTankSimulationGeneralConcept,淤泥流入,出水,固体缺陷,主要旋流两面中心+二次流动.利用“茶杯效应”.较重的固体聚集在中心.,目标：模拟主旋转流动和二次流动。只有流动的计算-没有固体.,29,在圆槽总共速度的实验结果,上图显示槽附近收集数字代表一小部分的入口速度预计沿中心线强涡旋在涡流1.2倍速度大于入口速度,涡流,30,模拟结果IORDER=1,Initallysolid-bodyrotationNovortexingseennearaxis1storderadvectiondissapatesvortexTryhigherordermomentumadvection,默认的1st平流IORDER=1，计算效率是很多问题，但太多的对角线速度分量扩散而准确预测旋流。,31,模拟结果IORDER=2,速度是不是合理不稳定IORDER=2不建议自由面模拟,不稳定,由于IORDER=2是已知的强循环与自由表面流动的不稳定，但预计不会履行好这个问题。事实上，模拟很快变得不稳定，产生不合理的流速。若要使用IORDER=2针对此问题，用户需要取代对称边界条件的自由表面。然而，并非涡旋预测以及自由表面时被删除，因为不能动用流动在槽中心。这降低了在轴角动量增加。,32,模拟结果IORDER=4,流体运动初始化为solid-body旋转良好的预测轴附近看到涡旋三阶捕获良好的流动以及平流,IORDER=4在流动中的一个圆形槽旋流模拟预测。在流体运动是初始化为坚实的物体旋转，启动是有利益于稳态流。使用IORDER=4在槽附近的出口有强烈的涡旋预测。水坑（沉没在槽中心）是勉强可以解决涡旋约1