CFD中文网

@李东岳 好的老师

@李东岳 李老师好，我跟我老师说了一下，他对于指导和合作论文什么的是没啥问题的，而且能指导话对我自己来说也是个很好的学习机会。我是建筑环境与能源应用工程专业的，现在准备做室内流场快速计算和这个颗粒物矩方法这块工作，原来我们老师跟一个做室内快速流场计算和颗粒物快速计算的老师有合作，那个老师做的时候就是用openfoam做的，所以openfoam植入这块和方向上我们应该是没问题。

@bestucan 好的，感谢老师，我操作学习一下

@李东岳 好的谢谢李老师，我看一下这篇论文。 这个算是博士课题，主要是学术上的，openfoam是想要弄，不过用的不熟，还在学习中，本来是想在Fluent上整明白再去整openfoam的。

@bestucan 是的老师，我这个是对数分布的矩方法不是TEMOM，我是看了文献[1]他用的这个经典矩方法和TEMOM处理这个颗粒物沉积计算的，右边的表达式是在文献[2]中看见的。然后，老师这个无量纲化我算三维的也是需要先把变量都无量纲化，然后解三个三维的无量纲标量方程么。
感谢老师的解答，我看看这个文献具体操作
[1]Guangping Xu, Jiasong Wang,CFD modeling of particle dispersion and deposition coupled with particle dynamical models in a ventilated room. DOI:https://doi.org/10.1016/j.atmosenv.2017.07.027.
[2]S.H. Park, K.W. Lee,Analytical solution to change in size distribution of polydisperse particles in closed chamber due to diffusion and sedimentation. DOI:https://doi.org/10.1016/S1352-2310(02)00673-8.

@李东岳 是的李老师，我计算0维没问题，但是到三维计算就有问题，不过Fluent里面的格式我都试了一下都不好使，植入格式这个我还不太了解，我再看看

各位老师、各位前辈好：
最近在通过CFD和MOM计算颗粒物浓度的时候遇到了一些问题和困惑，希望有前辈、老师能够指点一下。
首先这个模型仅考虑了颗粒物的沉积损失，不考虑成核等过程。


然后通过对数正态分布矩方法把这个GDE方程写出矩的形式，如下：（下面分别是i=2和i=2.6的情况下）


为了确保准确，我先照着文献数值计算了一下相关内容：（包括沉降率、粒径分布变化等等）

然后这个仅考虑了前三阶矩，矩的初始值由下面的式子进行计算，分别是用平均体积和平均半径计算，这就出现了第一个问题，这两个计算公式算出来的矩初始值不一样。我再确保跟文献的结果没太大差别后，进行了下一步工作，就是通过Fluent中的UDF将这个颗粒物沉降编写为源项，将三阶矩分别编写成三个UDS方程来进行求解。但是在求解时候发现老是发散，然后基于平均体积、i=2和i=2.6计算的时候，发现M1和M2总是变成了负数，或者M1也是负数。


下面是计算的一些图片，


下面是M0、M1和M2的残差变化，这个一看就有问题

因此，想跟各位老师请教一下这个问题，首先是为什么这两个公式算出来的矩初值不一样，然后就是发散的问题，我想请教下这么用UDF进行计算是正确的么，流场和MOM共同求解的话正确解法是什么样的。

@李东岳 都设置沉积边界条件的话，颗粒物算的就大很多了。它文献中是1微米的颗粒物稳态情况下充满整个房间，但我算的10微米，也是充满整个房间，跟文献中的数据差别很大。

@李东岳 李老师好，下面这个图是现在算的速度场

各位大佬前辈好，我在用UDF照着Chen F, Yu SCM, Lai ACK (2006). Modeling particle distribution and deposition in indoor environments with a new drift–flux model. . 这篇文献做颗粒物模拟时候遇到些问题。我用沉积通量做边界条件或者沉积源项时，模拟结果跟文献结果相差很大。按照文献Effect of supply/exhaust diffuser configurations on the contaminant distribution in ultra clean environments:Eulerian and Lagrangian approaches.中定义的仅在地板上有沉积通量边界条件的话能跟结果对上，但是感觉很有问题，因为它把其他边界处的颗粒物浓度设定为0，并且我算其他粒径颗粒物的话，跟数据有差别。下面的图依次是沉积源项、仅地板有沉积边界条件和都设置有沉积通量边界条件，最后两张图是其他文献计算的云图和我用沉积源项计算的对比。请教下有没有做过相关工作的前辈，有没有好的建议。

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

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

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

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

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

各位大佬前辈好：
最近照下面这篇文献的气流场进行模拟时候，我的模拟结果与文章给出的实验数据对不上(比较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.