有关通风室气流模拟的问题
-
各位大佬前辈好:
最近照下面这篇文献的气流场进行模拟时候,我的模拟结果与文章给出的实验数据对不上(比较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?什么湍流模型
-
@李东岳 李老师好,这个是单相流,稳态模拟,用的RNG湍流模型。
-
simpleFoam?
但是我看文章里面是两相流。你的实验也是单相流数据么
-
@李东岳 文章是气-固两相流模型,但是做的时候是先算气流场,再算的固相。我这个是按照它的模型先做的气流场模拟,但是跟文章中的速度实验数据对不上。
-
那就不太好弄了。有时候确实对不上。不过你这个大体上没问题,就是下面有些逆压梯度没模拟出来。只能建议用别的RANS算算。或者用LES算,然后算调取UMean对比下。其他也没啥好办法。
-
@李东岳 好的好的,谢谢李老师,那我再试试把
-
是否可以考虑先进行网格无关性验证,或者 可以修改一下RANS中的湍流模型
-
@ibelief 谢谢大佬,我做了一下网格无关性验证,现在气流场算的和文献数据差距不是特别大了。不过后面还需要算颗粒物场,用不同的湍流模型,得到湍流粘度不一样,然后颗粒物场算的也不一样。壁面函数我用了standard wall function和scalable wall function,standard wall function算出的y+小于1,scalable wall function算出的y+最大在7。这地方有些不明白,wall function是按照什么原则确定的。
-
@wyboo 这部分一言难尽,可以参考一下《数值传热学》第九章相关内容。
颗粒流动中,湍流脉动是一个非常重要的因素,在大量的管道通风、沉积等研究中广泛讨论,为了降低工作量,你可以先查阅一下相关文献。如果考虑壁面沉积的话,建议看一下各向异性湍流模型和各向同性湍流模型之间的差异。
你引用的LAI 课题组就有相关的文章,还有和他们关系紧密的香港几个课题组,他们对这方面讨论比较充分。 -
@ibelief 好的,谢谢大佬。我看书研究研究