CFD中文网

@bestucan 我看了目前用到这个 -I. 的地方，包括reactingParcelFoam跟sprayFoam，都有子文件夹，比如sprayFoam下的egineFoam，应该跟compile子求解器有关

@bestucan
我做的是droplet，属于颗粒的。组分说的是这个droplet是一种溶液，只是我不清楚这里面的溶液到底算不算组分，溶液里面的水是可以蒸发的，溶质又无法蒸发，溶质跟溶液到底算不算俩种phase，还是按照同一种phase的不同组分来处理

@bestucan 这个混乱是可以解决的，在引用的库文件前面加入路径就行了

-L$(FOAM_USER_LIBBIN) \ -llagrangianIntermediate \

而且我在其他论坛上也看到，lagrangian的代码修改时必须把名字改成一致的，不然会出现很多warning，

气体多孔介质的扩散系数是一个与孔隙率和曲曲率有关的量，注意关注这两个量的设置。

@东岳 李博，我准备试试，之前写过fluent多面体的网格转换程序，不过这个发不了论文:mihu:

Really? 那得谢谢王博士！

@Izumi 这个没添加，我回去试试

@东岳 这个边界条件已经整合到exteranlHeatFluxTemperature里面了，算的也是温度梯度，也就是$q=-\kappa \frac{\partial T}{\partial n}$，这种边界按照FDS技术手册里面应该是DNS的计算方法。
LES计算时给定的方法是$q=h \Delta T$,$h$是对流换热系数。openfoam里面的边界条件要么是固定值要么是梯度，这种差值不知道在of里面能不能实现

在模拟co2、水蒸气的对流扩散的时候，准备用[Fluent theory guide]{http://www.afs.enea.it/project/neptunius/docs/fluent/html/th/node128.htm}的species transport方程，其中对流项和扩散项都没有什么问题，主要是重力作用的加入，因为co2水蒸气都比空气重，在实际中重力作用很明显的
对于重力项之前别人推荐参考buoyantFoam，我也根据[东岳老师]{http://www.dyfluid.com/buoyantPimpleFoam.html}的这篇讲解，在OpenFOAM里面加入speices transport引起的重力变化。主要还是加入了密度场，通过不同密度对速度方程中的浮力项的作用，来实现重力的影响，并没有直接对co2对流扩散有作用。所以主要还是靠的速度变化，引起co2向下流动的现象。这种变化感觉非常不明显，比如人呼出的co2才占整体重量的4%，实际模拟中呼出co2基本不太向下聚集，Fluent在species transport中是不是也加入了源项？0_1519956843102_scalarTransportBuoyant.zip

@李东岳 最后还是用了后处理的方法进行计算，又要做实验又要做模拟有点赶，用的是老板之前研究的方法几何平均粒径、几何标准差与浓度的关系，算了浓度再来得到颗粒数关系。不过布朗运动要加进去我感觉得用群体平衡方法。

@一二 在 LES介绍的文章的一个公式 中说：

嗯，下面是我推的
\begin{equation}
\overline{\mathbf{S}}:\mathbf{\tau}+C_{e}\frac{k_{sgs}^{1.5}}{\Delta}=0 \
left = \overline{\mathbf{S}}:[\frac{2}{3}k_{sgs}\mathbf{I}-2\mu_{sgs}\mathbf{dev}(\overline{\mathbf{S}})]+k_{sgs}^{1.5}\frac{C_{\epsilon}}{\Delta}
\end{equation}
因为
\begin{equation}
\mu_{sgs}=C_{k}\Delta\sqrt{k_{sgs}} \
= \overline{\mathbf{S}}:[\frac{2}{3}k_{sgs}\mathbf{I}-2C_{k}\Delta\sqrt{k_{sgs}}]+k_{sgs}^{1.5}\frac{C_{\epsilon}}{\Delta} \
= \sqrt{k_{sgs}}(\frac{C_{\epsilon}}{\Delta}k_{sgs}+\frac{2}{3}\mathbf{tr}(\overline{\mathbf{S}})\sqrt{k_{sgs}}-2C_{k}(\mathbf{dev}(\overline{\mathbf{S}}):\overline{\mathbf{S}})) \
= ak_{sgs}+b\sqrt{k_{sgs}}-c\
= right = 0
\end{equation}
其中
\begin{equation}
a=\frac{C_{\epsilon}}{\Delta} \
b=\frac{2}{3}tr(\overline{\mathbf{S}}) \
c=2 \Delta C_{k}(\mathbf{dev}(\overline{\mathbf{S}}):\overline{\mathbf{S}} \
\sqrt{k_{sgs}}=\frac{-b+\sqrt{b^2+4ac}}{2a}
\end{equation}
当为不可压缩流体时$tr{\overline{\mathbf{S}}}=0$,那么$b=0$、$c=2 \Delta C_{k} (\mathbf{dev}(\overline{\mathbf{S}}):\overline{\mathbf{S}})$,就可以得到$k_{sgs}=\frac{c}{a}=\frac{2C_{k}||\overline{\mathbf{S}}||^{2}{}\Delta^2}{C_{\epsilon}}$

不是很清楚你说的那种涡脱落现象。不过用户指南中我记得有一章曾经表示，大体意思为：“如果用稳态求解器模拟比较强的非稳态问题会出现残差的震荡”。

能否附稳态求解涡街的文献？谢谢。

@李东岳 哈哈，其实是从nodebb的帮助中抄来的:cheeky:

@micro 网格数量越多，网格尺度越小，计算肯定更加不稳定，您可以查一下CFL数