有关通风室气流模拟的问题

[WYBoo]

各位大佬前辈好：
最近照下面这篇文献的气流场进行模拟时候，我的模拟结果与文章给出的实验数据对不上(比较x方向速度)，特别是地板附近和靠近出风口上半空间的速度。几何模型很简单，进风口设为速度入口(0.225m/s)，出发口为outflow，尺寸一样都是0.04m×0.04m，风口中心距离最近的水平壁面均为0.04m。结构化网格和非结构化网格，各种壁面函数都用了一下，但还是相差很大。而且壁面附近的网格粗化后，Y+还是在10一下。想请教下各位前辈有没有好的建议。
Chen F, Yu SCM, Lai ACK (2006). Modeling particle distribution and deposition 22 in indoor environments with a new drift–flux model. Atmospheric Environment, 23 40:357-367.

[李东岳]

这是单相流还是多相流？稳态瞬态？RANS？什么湍流模型

[WYBoo]

@李东岳 李老师好，这个是单相流，稳态模拟，用的RNG湍流模型。

[李东岳]

simpleFoam？

但是我看文章里面是两相流。你的实验也是单相流数据么

[WYBoo]

@李东岳 文章是气-固两相流模型，但是做的时候是先算气流场，再算的固相。我这个是按照它的模型先做的气流场模拟，但是跟文章中的速度实验数据对不上。

[李东岳]

那就不太好弄了。有时候确实对不上。不过你这个大体上没问题，就是下面有些逆压梯度没模拟出来。只能建议用别的RANS算算。或者用LES算，然后算调取UMean对比下。其他也没啥好办法。

[WYBoo]

@李东岳 好的好的，谢谢李老师，那我再试试把

[ibelief]

是否可以考虑先进行网格无关性验证，或者 可以修改一下RANS中的湍流模型

[WYBoo]

@ibelief 谢谢大佬，我做了一下网格无关性验证，现在气流场算的和文献数据差距不是特别大了。不过后面还需要算颗粒物场，用不同的湍流模型，得到湍流粘度不一样，然后颗粒物场算的也不一样。壁面函数我用了standard wall function和scalable wall function，standard wall function算出的y+小于1，scalable wall function算出的y+最大在7。这地方有些不明白，wall function是按照什么原则确定的。

[ibelief]

@wyboo 这部分一言难尽，可以参考一下《数值传热学》第九章相关内容。
颗粒流动中，湍流脉动是一个非常重要的因素，在大量的管道通风、沉积等研究中广泛讨论，为了降低工作量，你可以先查阅一下相关文献。如果考虑壁面沉积的话，建议看一下各向异性湍流模型和各向同性湍流模型之间的差异。
你引用的LAI 课题组就有相关的文章，还有和他们关系紧密的香港几个课题组，他们对这方面讨论比较充分。

[WYBoo]

@ibelief 好的，谢谢大佬。我看书研究研究