关于置换通风室内气流模拟不收敛问题

[深蓝]

我想要模拟室内置换通风的温度速度场的分布情况，几何模型很简单，如图1所示。我使用的模拟工具是fluent和airpak。湍流模型选用的标准k-e模型、RNGk-e模型、室内零方程湍流模型（都有试过），空气密度设为不可压理想气体、Boussinesq模型（都有试过），墙面边界条件为固定热流密度，人员和电脑边界条件为固定热流密度，风口入口边界条件：速度=0.2m/s，温度=20度，湍动能=4.25%，水力直径=0.58m，回风口边界条件试过自由出流、压力出口。使用SIMPLE、SIMPLEC算过，压力采用PRESTO！，其他都是二阶迎风精度（一阶的也试过）。网格采用ICEM的结构网格，数量20万、50W、80W，后来也考虑过Y+值。

我的问题是，我现在不管怎么处理都收敛不了，尝试过调温度、密度的松弛因子都不收敛，而且收敛残差值都类似，都是先下降，然后上升，再不断振荡下去，我使用airpak建立相同模型，残差也是类似结果。如图2。

这个问题卡了我快两个月，找了各种原因都不行，请大家帮帮找找原因，谢谢大家了。
我的模型传到网盘上了：链接: https://pan.baidu.com/s/1-glgZNYdPPbJ-xnvyuH_GQ 提取码: tmun

[李东岳]

说一下你灯、回风口、电脑人员桌子的边界条件，很有可能是边界条件的问题

[a602657563]

建议列表表示下边界条件设置，实在不行可以试试CFX，它有个Openning 出口边界，对不可压缩稳态问题收敛性较好。

[仲阳]

您好！请问您的问题有解决了嘛？我在模拟置换通风的过程中也遇到了该问题。

[李东岳]

这个要是能提供网格，我可以用openfoam试试

[仲阳]

@李东岳 链接：李老师，您好！这是我的msh以及case文件
https://pan.baidu.com/s/1NR7js6Lmw1W5_0VoLwr0Jg?pwd=11kt
提取码：11kt
--来自百度网盘超级会员V4的分享

[李东岳]

@仲阳 好像你不是楼主。fluent的设置我不会看。你那个要设置的话，恐怕要说一下边界条件风速啥的。

[仲阳]

@李东岳 李老师，下面是我的边界条件的设置
进风口：速度入口0.5m/s，温度26℃，湍流强度4%，水利直径0.266m
回风口：压力出口（标准大气压），温度26℃
壁面：常壁温/固定热流密度都有试过
使用SIMPLE，压力采用PRESTO！，其余项是二阶迎风精度。
网格采用ICEM的结构网格

[李东岳]

@仲阳 我试试

[李东岳]

@仲阳 这个网格用openfoam转化出错了

[仲阳]

@李东岳 好的好的 谢谢李老师 我再研究一下

[深蓝]

@仲阳 迟复为歉，我的方法是修改壁面普朗特数为0.005，参考张腾飞论文的论文：ZHANG T T, ZHOU H, WANG S. An adjustment to the standard temperature wall function for CFD modeling of indoor convective heat transfer[J]. Building and environment, 2013, 68: 159–169. 如果有问题可以进一步与我联系，乐意回答，我已经给你发送邮件。

[仲阳]

@深蓝 感谢感谢！