求助fluent计算两相流鼓泡塔,连续性方程残差收敛困难?

采用2D轴对称模型,模拟鼓泡塔。使用压力基求解器,瞬态计算模式。模型选用欧拉-欧拉两相流模型(初始相为水相,第二相为气相)。湍流模型选用k-e RNG 模型,标准壁面函数,相湍流选择dispersed。相间作用力考虑:虚拟质量力、曳力(Tomiyama)、湍流作用选择Sato模型,表现张力设置0.07.边界条件,入口为velocity-inlet,气相体积分数为1,气速0.025m/s。出口边界设置pressure-out,气相回流体积分数为1.离散格式均采用二阶,体积分数采用QUICK格式。采用Coupled计算模式(采用phase-coupled simple总是在压力出口边界产生回流),时间步长为0.01,最大迭代50.

模型一开始计算的时候残差曲线是收敛的,计算100步时间步长后连续性方程的残差怎么都不收敛(降低courant数,减小松弛因子,减小迭代步长,调整网格密度和质量都试过了)。表现为连续性残差在一个水平线附近震荡,一下10^-3一下10^-2就是无法收敛。

模型是一个0.1*2m的矩形,初始化后patch左边0.1*1部分是水,其他部分为气体。

残差震荡后分析residual contour显示在塔高1m气液交界面处的mass imbalance最大。

请问有没有相关方面的大牛知道为什么该问题的残差曲线难以收敛?还有就是采用phase-coupled simple时总是产生回流,而使用coupled计算时 ,却没有回流的问题?


求助fluent计算两相流鼓泡塔,连续性方程残差收敛困难?的图1求助fluent计算两相流鼓泡塔,连续性方程残差收敛困难?的图2


邀请回答 举报

共1个回答

默认排序
  • 随波逐流
    2019年3月31日

    难以收敛可以采用更小的时间步长试试;
    耦合算法和分离算法的具体区别可以看帮助;简要而言,耦合算法是连续性方程和动量方程一起来算的,消耗内存更多,但也容易收敛;分离算法是连续性方程得到速度,再由该速度算动量方程,并不断重复修正的过程,消耗内存少,但收敛相对迭代更久;

    0评论 采纳答案 举报

  • 求助fluent计算两相流鼓泡塔,连续性方程残差收敛困难?的相关案例教程

    基于Fluent风沙两相流路基周围流场数值模拟分析 本课程基于Fluent风沙两相流路基周围流场数值模拟分析,从建模到网格划分,以及相关设置进行了讲解,并运用tecplot软件进行了后处理;其中许多点作者花费了大量时间进行理解,例如风速轮廓线的设置以及求解设置中相关参数的取值,读者可在本课程基础进行拓展延申,节省前期大量时间,使得可以让许多时间花费在研究上。 课程内容: 章节1、建立3d路基模型.
    1. Continuity残差不收敛:这和Simple有关,Simple根据连续性方程推导出压力修正方法求解压力,流场耦合被过渡简化,使压力修正方程不能准确反映流场变化; 解决办法:1. 试验Simplec方法; 2. 加密网格; 3. 对于结构网格,建议使用高阶格式(二阶迎风格式);对于非结构网格,除pressure格式(压力离散格式)保持standard格式不变,其他格式改用高阶; 4. 分离
    0. 伪瞬态作用 为什么要使用伪瞬态的算法?伪瞬态的作用实际上是增加收敛性的,当你的稳态计算收敛性不好时,可以将稳态计算更改为伪瞬态计算,收敛性会增强。 当然还可以通过前面所说的降低松弛因子的方式来增强收敛性。 但是,伪瞬态并不是真正的瞬态,它虽然会出现时间步长这种概念,但是在每个时间步长并不收敛,而只是最终的计算结果收敛,因此当计算只考虑稳态结果时可以使用伪瞬态算法,而如果考虑某时刻的结果,则必
    最近尝试了Fluent与Rocky DEM单向耦合,记录下使用的心得。Rocky DEM被Ansys收购之后,它们之间的耦合几乎无需复杂的配置了,RockyDEM安装上即可实现和Fluent的耦合。这里有一个双层筛的清选段模型,出口装有风机,设计风量为12000m³/h。右上入口落入原粮、清杂(包含长短不一的水稻茎秆),其中原粮占混合物总重的97.5%;右侧通道中有一料滚以52rpm的速度顺时针旋
    01 前言 最常用的热水器形式为真空玻璃管热水器,这种热水器由真空集热管和贮热水箱组成,真空集热管由内外两根同轴硼硅玻璃管组成,外形犹如一根拉伸过后的热水壶内胆,一端开口连接于上部的热水箱,另一端密封成半圆球头。 在内玻璃管上附有选择性吸收涂层,其辐射吸收特性会随着入射辐射波长的改变而改变;为了防止对流和热传导带来的热损失,内玻璃管和外玻璃管之间被抽成真空。为防止不必要的热损,在贮热水箱外包有保温