请教关于射流模拟的一个问题

ThermalSpray · 发表于 2003-12-2 13:16:02

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？注册

x

[这个贴子最后由ThermalSpray在 2003/12/02 01:22pm 第 2 次编辑]

我在模拟拉瓦尔喷嘴射流的时候，fluent计算不能得到收敛的结果，请大虾多多给与指点。
拉瓦尔喷管的入口直径是26mm，出口是24mm，喉部直径是4.2mm。用gambit将拉瓦尔喷嘴进行网格划分，枪外射流计算域去5D×12D（D为出口直径），其网格划分如下：
Gambit中边界条件设置如下：边界1：PRESSURE_INET；边界2：AXIS；
边界3：PRESSURE_OUTLET；边界4、5： PRESSURE_FAR_FIELD
将网格存为nozzle.msh后导入Flent中后，Fluent的设置如下：
Solver：Coupled－Implicit
Energy：打开
Viscious：Spalart－Allmaras，参数设置为默认值
Materials：Ideal gas
Operating Pressure：0
Boundary Conditions：
      Pressure Inlet 设置如下：
               Gauge Total Pressure[pascal]：650000
               Super/Initial Gauge Pressure[pascal]：615000
               Total Temperature[k]：293
               Direction Specification Method: Normal to Boundary
               Turbulence Specification Method: Turbulent Viscosity Ratio
               Turbulent Viscosity Ratio: 1,
      Pressure Outlet:
               Gauge Pressure[pascal]：101325
               Backflow Total Temperature[k]：293
               Direction Specification Method: Normal to Boundary
               Turbulence Specification Method: Turbulent Viscosity Ratio
               Turbulent Viscosity Ratio: 10,
      Pressure Far-field:
            Gauge Pressure[pascal]:101325
            Mach Number:0
            Temperature[k]:293
            Turbulence Specification Method: Turbulent Viscosity Ratio
            Turbulent Viscosity Ratio:10
      Solution Initialization:
            Compute from: pressure inlet
这样设置下进行计算，结果不能收敛。实际上计算一开始的时候就提示divergence detected。我是一位新手，不知道我上面的设置哪里出了问题。请大侠不吝赐教呀！
Many thanks!

cjjj2000 · 发表于 2003-12-2 18:22:38

1,首先Turbulent Viscosity Ratio:10，太大了
以上不是主要问题，主要是因为你的网格问题。
在喉部附近还要再密一些，其余部位的网格是有喉部开始逐渐变疏得网格（和喉部一样密也行）
再则不要追求结构化网格，全场非结构化也挺好。重要的是得到结果，对吧？

ThermalSpray · 发表于 2003-12-2 19:41:00

非常感谢你的指点，
我再把网格修改一下，看能不能得到收敛的结果。
thanks again!

ThermalSpray · 发表于 2003-12-3 12:22:24

哈，计算终于收敛了。非常感谢jeffchenjw大虾（不知他上不上这个论坛了，但在这里还是要感谢一下了：P）的悉心指点，百忙之中还抽出时间帮我算了一下

还要感谢cjjj2000的指点。结合jeffchenjw的意见，我原来存在的计算在网格划分、边界条件以及参数设置上都有些问题。在这里，我把各个计算参数再列一些，以供作射流分析的新手们作个简单的参考。
Gambit中边界条件设置如下：边界1：PRESSURE_INET；边界2：AXIS；
边界3、4、5：PRESSURE_OUTLET
将网格存为nozzle.msh后导入Flent中后，Fluent的设置如下：
Solver：Coupled－Explicit
Energy：打开
Viscious：k-epsilon(2 eqn), k-epsilon Model: RNG
Materials：Ideal gas, Viscosity:sutherland
Operating Pressure：0
Boundary Conditions：
      Pressure Inlet 设置如下：
            Gauge Total Pressure[pascal]：650000
            Super/Initial Gauge Pressure[pascal]：615000
            Total Temperature[k]：293
            Direction Specification Method: Normal to Boundary
            Turbulence Specification Method: Intensity and Viscosity Ratio
            Turbulence Intensity: 5
            Turbulent Viscosity Ratio: 2,
      Pressure Outlet:
            Gauge Pressure[pascal]：101325
            Backflow Total Temperature[k]：293
            Direction Specification Method: Normal to Boundary
            Turbulence Specification Method: Intensity and Viscosity Ratio
            Turbulence Intensity: 2
            Turbulent Viscosity Ratio: 2,
      Solution Controls: Courant Number:0.5
                     Flow: First Order Upwind
                     Multigrid Levels:5
      Solution Initialization:
            Gauge Pressure[pascal]:101325
  这样设置下进行计算，结果能很好的收敛。

松鸟菜菜子 · 发表于 2003-12-3 14:35:55

计算结果能贴上来看看吗？

jeffchenjw3 · 发表于 2003-12-3 17:28:04

不要那么客气。我帮你重新作了一套网格，二阶精度的收敛很好，计算效果自然更好，1400步左右即可收敛，不知你是否需要。
这个站点我有时会上来看看。

ThermalSpray · 发表于 2003-12-4 10:17:03

下图为一阶精度下的速度分布图

ThermalSpray · 发表于 2003-12-4 10:20:43

上图为一阶精度下的Much Number图

zhangshaowei · 发表于 2003-12-4 12:12:36

还需要看波系分布

fengbeibin · 发表于 2005-10-13 10:28:56

[这个贴子最后由fengbeibin在 2005/12/22 02:45pm 第 2 次编辑]

[DISABLELBCODE]

richard001 · 发表于 2005-10-13 12:35:00

什么配置的电脑，算了多长时间啊？网格数多少啊？

眼镜蛇 · 发表于 2005-10-17 10:54:59

有收获

porsche75339 · 发表于 2005-10-17 16:41:42

不錯,果然高手如雲啊!

yandongjie01 · 发表于 2005-11-14 14:17:47

有帮助

blck · 发表于 2005-12-4 10:10:59

如果要计算一个速度不高的水射流但是加上vof如何计算呢？

		自动登录	找回密码
密码			注册

请教关于射流模拟的一个问题

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。