版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、continuity收敛问题fluent求解器设置fluent求解器设置主要包括:1、压力-速度耦合方程格式选择2、对流插值3、 梯度插值4、压力插值下面对这几种设置做详细说明。一、压力-速度耦合方程求解算法fluent 中主要有四种算法:simple, stmplec, ptso, fsm(1) simple (semi-implicit method for pressure-linked equations)半隐式连 接压力方程方法,是fluent的默认格式。(2) slmplec(slmple-consistent) o对于简单的问题收敛非常快速,不对压力进 行修正,所以压力松弛因子可
2、以设置为1(3) pressure-implicit with splitting of operators (piso)。对非定常流动 问题或者包含比平均网格倾斜度更高的网格适用(4) fractional step method (fsm)对非定常流的分步方法。用于nita格式,与 p1so具有相同的特性。二、对流插值(动量方程)fluent有五种方法:一阶迎风格式、幕率格式、二阶迎风格式、musl三阶格式、quick格式(1) fluent默认采用一阶格式。容易收敛,但精度较差,主要用于初值计算。(2) power lar.幕率格式,当雷诺数低于5时,计算精度比一阶格式要高。(3) 二阶
3、迎风格式。二阶迎风格式相对于一阶格式來说,使用更小的截断误差, 适用于三角形、四面体网格或流动与网格不在同一直线上;二阶格式收敛可能比 较慢。(4) musl(monotone upstream-centered schemes for conservation 1aws) 当地3阶离散格式。主要用于非结构网格,在预测二次流,漩涡,力等时更精确。(5) quick (quadratic upwind interpolation)格式。此格式用于四边形/六 面体时具有三阶精度,用于朵交网格或三角形/四面体时只具有二阶精度。三、梯度插值梯度插值主要是针对扩散项。fluent 有三种梯度插值方案:g
4、reen-gauss cell-based, green-gauss node-based, least-quares cell based.(1) 格林-高斯基于单元体。求解方法可能会出现伪扩散。(2) 格林-高斯基于节点。求解更精确,最小化伪扩散,推荐用于三角形网榕 上(3) 基于单元体的最小二乘法插值。推荐用于多面体网格,与基于节点的格林- 高斯格式具有相同的精度和格式。四、压力插值压力基分离求解器主要有五种压力插值算法。(1) 标准格式(standard) o为fluent缺省格式,对大表妹边界层附近的曲线 发现圧力梯度流动求解精度会降低(但不能用丁流动中压力急剧变化的地方 此时应该使
5、用presto!格式代替)(2) presto!主要用于高旋流,压力急剧变化流(如多孔介质、风扇模型等), 或剧烈弯曲的区域。(3) linear (线性格式)。当其他选项导致收敛困难或出现非物理解时使用此格 式。(4) second order (一阶格式)。用于可压缩流动,不能用于多孔介质、阶跃、 风扇、vof/m1xture多相流。(5) body force weighted体积力。当体积力很犬时,如高雷诺数自然对流或 高凹旋流动屮采用此格式。(1) 连续性方程不收敛是怎么回事?在计算过程中其它指数都收敛了 ,就continuity不收敛是怎么回事。这和fluent程序的求解方法sim
6、ple有关。simple根据连续方程推导出压 力修正方法求解压力。曲于连续方程中流场耦合项被过渡简化,使得压力修正方 程不能准确反映流场的变化,从而导致该方程收敛缓慢。你可以试验simplec方法,应该会收敛快些。在计算模拟屮,continuity总不收敛,除了加密网格,还有别的办法吗? 别的条件都己经收敛了,就差它口己了,述有收敛的标准是什么?是不是到了一 定的尺度就能收敛了,比如10-e5具体的数量级就收敛了continuity是质量残差,具体是表示本次计算结果与上次计算结果的差别, 如果别的条件收敛了,就差它。可以点report,打开里面flux选项,算出进口 与岀口的质量流量差,看它是
7、否小于0.5%如果小于,可以判断它收敛.(2) fluent残差曲线图中continuity是什么含义?是质量守恒方程的反映,也就是连续性的残差。这个收敛的快并不能说明你的计 算就一定正确,还要看动量方程的迭代计算。表示某次迭代与上一次迭代在所有 cells积分的差值,continuty表示连续性方程的残差(3) 要加速continuity收敛该设置那些参数?感觉需要调整courant numberfluent中courant number是在耦合求解的时候才出现的。正确的调整,可 以更好地加速收敛和解的增强稳定性。courant number实际上是指时间步长和空间步长的相对关系,系统自动减
8、 小courant数,这种情况一般出现在存在尖锐外形的计算域,当局部的流速过 大或者压差过人时出错,把局部的网格加密再试一 ho在fluent中,用courant number来调节计算的稳定性与收敛性。一般来 说,随着courantnumber的从小到大的变化,收敛速度逐渐加快,但是稳定性 逐渐降低。所以具体的问题,在计算的过程屮,最好是把ourant number从小 开始设置,看看迭代残差的收敛情况,如杲收敛速度较慢而且比较稳雄的话,可 以适当的增加courant number的大小,根据自己具体的问题,找出一个比较合 适的courant number,让收敛速度能够足够的快,而且能够保
9、持它的稳定性。 个人认为这应该和你采用的算法有关simple算法是根据连续方程推导出压力修正方法求解压力。由于连续方程屮流场耦合项被过渡简化,使得压力修正方程不能准确反映流场的 变化,从而导致该方程收敛缓慢。试着用stmplec算法看看。(4) 正在学习fluent,模拟圆管内的流动,速度入口,出口。utflow运行后xy 的速度很快就到le-06了,但是continuity老是降不下去,维持在le-00和le-0 3之间,减小松弛因子好像也没什么变化大家有什么建议吗?你查看了流量是否平衡吗?在rcport->f lux里面操作,mass f 1 ow rate, 把所有进出口都选上,c
10、ompute 下,看看nut flux是什么水平,如果它的值 小于总进口流量的1%,并且其他检测量在继续迭代之后不会发生波动,也可以认 为你的解是收敛的。造成连续方程高残差不收敛的原因主要有以下几点:1 网格质量,主要可能是相邻单元的尺寸大小相養较大,它们的尺寸z比最好控 制在1.2以内,不能超过1.4.2. 离散格式及压力速度耦合方法,如果是结构网格,建议使用高阶榕式,如2阶 迎风格式等,如果是非结构网格,除pressure保持standard格式不变外,其他 格式改用高阶格式;压力速度耦合关系,如果使用simple, s1mplec, p1s0等 segcratcd solver对联系方程
11、收敛没有提高的话,可以尝试使用coupled solveto另外,对于梯度的计算,不论使用结构或非结构网格,都可以改用 node-based来提高计算精度。一些情况:1. 监测流场某个变量来判断收敛更合理一些.2. 网格质量.3. velocity inlet boundary conditions are not appropriate for compressible flowproblems.(4) 要加速continuity收敛该设置那些参数?感觉需要调整courant numberfluent屮courant number是在耦合求解的时候才出现的。正确的调整,可 以更好地加速收敛和
12、解的增强稳定性。courant number实际上是指时间步长和空间步氏的相对关系,系统口动减 小courant数,这种情况一般出现在存在尖锐外形的计算域,当局部的流速过 大或者压差过大时出错,把局部的网格加密再试一下。在fluent中,用courant number来调节计算的稳定性与收敛性。-般来 说,随着courantnumber的从小到大的变化,收敛速度逐渐加快,但是稳定性 逐渐降低。所以具体的问题,在计算的过程中,最好是把ourant number从小 开始设置,看看迭代残差的收敛情况,如果收敛速度较慢而且比较稳定的话,可以适当的增加courant number的大小,根据自己具体的
13、问题,找出一个比较合 适的courant number,让收敛速度能够足够的快,而冃能够保持它的稳定性。 个人认为这应该和你采用的算法有关simple算法是根据连续方程推导岀压力修正方法求解压力。曲于连续方程中流场耦合项被过渡简化,使得压力修正方程不能准确反映流场的 变化,从而导致该方程收敛缓慢。试着用stmplec算法看看。ansys仿真计算代做:模态分析,瞬态动力学,谐响应分析和谱分析、械结构的疲劳、损伤,cfd流体;结构的强度评估和优化;企鹅:690294845利用fluent不收敛通常怎么解决? 、一般首先是改变初值,尝试不同的初始化,事实上好像初始化很关键,对于收敛。 、fluent
14、的收敛最基础的是网格的质量,计算的时候看怎样选择cfl数,这个靠经验 、首先查找网格问题,如果问题复杂比如多相流问题,与模型、边界、初始条件都有关系。 、有时初始条件和边界条件严重影响收敛性,曾经作过一个计算反反复复,通过修改网格,重新 肚义初始条件,包括具体的选择的模型,还有老师经常用的方法就是看看哪个因索不收敛,然后寻 找和它有关的条件,改变相应参数。就收敛了 、a.检查是否哪里设眾有误:比方用mni的unit建构的mesh,忘了scale;比方给足的边界条件不 合理。b从算至发散前几步,看presure分布,看不岀來的话,再算几步,看看问题大概岀在那个区域。 c网格,配合第二点作修正,就
15、重建个更漂亮的,或是更粗略的来处理。d再找不出来的话,换个solvere 、解决的办法是设几个监测点,比如出流或参数变化较大的地方,若这些地方的参数变化很小, 就可以认为是收敛了,尽管此时残值曲线还没有降下來。 、调节松弛因子也能影响收敛,不过代价是收敛速度。亚松弛因子对收敛的影响所谓亚松驰就是将本层次计算结果与上一层次结果的差值作适当缩减,以避免由于差值过人而引起 非线性迭代过程的发散。用通用变量来写出时,为松驰因子(relaxation factors)/数值传热学-214fluent中的亚松驰:曲于fluent所解方程组的非线性,我们有必要控制变化。一般用亚松驰方法來 实现控制,该方法在
16、每一部迭代中减少了变化量。亚松驰最简单的形式为:单元内变量等于原来的值 加上亚松驰因子a与变化的积:分离解算器使用亚松驰来控制每一步迭代屮的计算变量的更新。这就意味着使用分离解算器解的方 程,包括耦合解算器所解的非耦合方程(湍流和其他标量)都会有一个和关的亚松驰因子。在fluent中,所有变量的默认亚松驰因子都是对大多数问题的最优值。这个值适合于很多问题,但是 对于一些特殊的非线性问题(如:某些湍流或若高rayleigh数自然对流问题),在计算开始时要慎重 减小亚松驰因子。使用默认的亚松驰因子开始计算是很好的习惯。如果经过4到5步的迭代残差仍然增长,你就需耍减 小亚松驰因子。有时候,如果发现残
17、差开始增加,你可以改变亚松驰因子重新计算。在亚松驰因子过大时通常会出现 这种情况。最为安全的方法就是在对亚松驰因子做任何修改z前先保存数据文件,并对解的算法做几 步迭代以调节到新的参数。最典型的悄况是,亚松驰因子的增加会使残差有少量的増加,但是随着解 的进行残差的増加又消失了。如果残差变化有几个量级你就需要考虑停止计算并回到最后保存的较好 的数据文件。注意:粘性和密度的亚松驰是在每一次迭代z间的。而且,如果直接解焰方程而不是温度方程(iu: 对pdf计算),基于焰的温度的更新是要进行亚松驰的。要查看默认的亚松弛因了的值,你可以在解控 制面板点击默认按钮。对于大多数流动,不需要修改默认亚松弛因子
18、。但是,如果出现不稳定或者发散你就需要减小默认 的亚松弛因子了,其中压力、动量、k和e的亚松弛因子默认值分别为02, 0. 5, 0. 5和0. 5。对j: simplec 格式一般不需要减小压力的亚松弛因了。在密度和温度强烈耦合的问题中,如相当高的rayleigh数的 自然或混合对流流动,应该对温度和/或密度(所用的亚松弛因子小于1.0)进行亚松弛。相反,当温 度和动量方程没有耦合或者耦合较弱时,流动密度是常数,温度的亚松弛因子可以设为1.0。对于其它的标量方程,如漩涡,组分,pdf变量,对于某些问题默认的亚更松弛可能过人,尤其是对 于初始计算。你可以将松弛因子设为0.8以使得收敛容易。 看
19、了流量是否平衡在report->flux里面操作,mass flow rate,把所有进出口都选上,compute一下,看看nut flux 是什么水平,如果它的值小于总进口流量的1%,并且其他检测量在继续迭代z后不会发生波动,也可 以认为你的解是收敛的。造成连续方程高残差不收敛的原因主要有以下几点:1网格质量,主要可能是相邻单元的尺寸大小相差较大,它们的尺寸之比最好控制在12以内,不 能超过1. 42.离散格式及压力速度耦合方法,如果是结构网格,建议使用高阶格式,如2阶迎风格式等,如果 是非结构网格,除pressure保持standard格式不变外,其他格式改用高阶格式;压力速度耦合关
20、系, 如果使用simple, s1mplec, piso等scgcratcd solver对联系方程收敛没有提髙的话,可以尝试使用 coupled solvero另外,对于梯度的计算,不论使用结构或非结构网格,都可以改用node-based来提 高计算精度。在fluent的中文帮助里,对收敛有比较详细地描述,建议去看看。收敛性可能会受到很多因素影响。大量的计算单元,过于保守的亚松驰因了和复杂的流动物理性质 常常是主要的原因。有时候很难确定你是否得到了收敛阶。没有判断收敛性的普遍准则。残差定义对于一类问题是有用的,但是有时候对其它类型问题会造成 谋导。因此,最好的方法就是不仅用残差来判断收敛性而
21、且还要监视诸如阻力、热传导系数等相关的 积分量。在fluent的帮助文件里给出了下面几种典型的情况。1、如果你对流场的初始猜测很好,初始的连续性残差会很小从而导致连续性方程的标度残差很人。 在这种情况下,检查未标度的残差并与适当的标度如入口的质量流速相比较是很有用的。什么是标度 残茅?就是选作用来标准化的残魁值,一般是取第五步吧,所以,一开是残弟就很小,那么,后面的 残差和它一比,值也很难收敛到很小数。2、对丁某些方程,如湍流量,较差的初始猜测可能会造成较高的标度因子。在这种情况下,标度 的残差最开始会很小,随后会呈非线性增长,最后减小。因此,最好是从残差变化的行为来判断收敛 性而不仅仅是残丼
22、的木身值来判断收敛性。你应该确认在几步迭代(比如说50步)之后残劳继续减小 或者仍然保持较低值,才能得出收敛的结论。另一个判断收敛性的流行方法就是要求未标度的残差减小到三阶量级。为了实现这一方法,fluent 提供了残差标准化,有关残差标准化的信息请参阅分离求解器残差定义和耦合求解器残差定义两节。 在这种方法中,要求标准化的未标度残差降到10=3。但是这种要求在很多情况卜fj能是不合适的。1、如果你提供了较好的初始猜测,残差可能不会降到三阶量级。比方说,在等温流动中,如果 温度的初始猜测非常接近最终值,那么能量残差根本就不会降到三阶量级。2、如果控制方程中包括的非线性源项在计算开始时是零,但是
23、在计算过程中缓慢增加,残差是 不会降到三阶量级的。例如,在封闭区域内部的口然对流问题,出于初始的均一温度猜测不会产生浮 力,所以初始的动量残差可能非常接近零。在这种情况下,初始的接近零的残差就不适合作为残差的 较好的标度。3、如果所感兴趣的变量在所有的地方都接近零,残差不会降到三阶量级。例如,在完全发展的 管流中,截面上的速度为零。如果这些速度初始化为零,那么初始的和最终的残差都接近零,因此也 就不能期待降三阶量级。在这种情况下,最好监视诸如阻力、总热传导系数等积分量来判断解的收敛。检查非标准化未标度 的残差来确定这个残差和适当的标度和比是不是很小也是很有用的。相反,如果初始猜测很差,初始的残差过大以至于残差下降三阶量级也不能保证收敛。这种情况对 于初始猜测很难的k和e方程尤其常见。在这里,检查你所感兴趣的所有积分量就很有用了。如果解是 不收敛的,你可以减少收敛公差fluent运行过程屮,出现残差曲线震荡原因?如何解决残差震荡的问题?残羞震荡对计算收敛性 和计算结果有什么影响?一.残差波动的主
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- Spark实时日志挖掘框架开发教程课程设计
- 教学设计:信用工具和外汇
- 教学媒体理论与实践教案
- 2026年oyo招聘测试题答案
- 2026年工厂保全测试题及答案
- 2026年特斯拉labview测试题及答案
- 2026年结业鉴定测试题及答案
- 2026年思吴江歌测试题及答案
- 2026年名校分班测试题及答案
- 成考毕业论文答辩
- 2026届山东青岛市平度市高三下学期模拟检测(一)历史试卷
- 2026年大理供电局项目制用工招聘(60人)笔试备考试题及答案解析
- 2026重庆市航空应急救援总队公开招聘34人笔试模拟试题及答案解析
- GB/T 47355-2026外包指南
- 中国鼻咽癌诊治指南(2026版)
- 【2025年】组织行为学试题及答案
- 2026年国开电大建筑制图基础形考考前冲刺测试卷及完整答案详解(历年真题)
- 市场监督管理局全流程市场监管工作手册(标准版)
- UPS电源售后服务承诺书范本
- 统编人教版五年级语文下册《田忌赛马》示范教学课件
- 急性气管支气管炎课件教学
评论
0/150
提交评论