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{\mathrm{\partial }T}{\mathrm{\partial }n}$，这种边界按照FDS技术手册里面应该是DNS的计算方法。
LES计算时给定的方法是$q=h\mathrm{\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{array}{}\text{(1)}& \overline{\mathbf{S}}:\tau +C_e\frac{k_{sgs}^{1.5}}{\mathrm{\Delta }}=0left=\overline{\mathbf{S}}:[\frac{2}{3}{k}_{sgs}\mathbf{I}-2{\mu }_{sgs}\mathbf{d}\mathbf{e}\mathbf{v}(\overline{\mathbf{S}})]+k_{sgs}^{1.5}\frac{C_{ϵ}}{\mathrm{\Delta }}\end{array}$$
因为
$$\begin{array}{}\text{(2)}& {\mu }_{sgs}=C_k\mathrm{\Delta }\sqrt{k_{sgs}}=\overline{\mathbf{S}}:[\frac{2}{3}{k}_{sgs}\mathbf{I}-2C_k\mathrm{\Delta }\sqrt{k_{sgs}}]+k_{sgs}^{1.5}\frac{C_{ϵ}}{\mathrm{\Delta }}=\sqrt{k_{sgs}}(\frac{C_{ϵ}}{\mathrm{\Delta }}{k}_{sgs}+\frac{2}{3}\mathbf{t}\mathbf{r}(\overline{\mathbf{S}})\sqrt{k_{sgs}}-2C_k(\mathbf{d}\mathbf{e}\mathbf{v}(\overline{\mathbf{S}}):\overline{\mathbf{S}}))=a{k}_{sgs}+b\sqrt{k_{sgs}}-c=right=0\end{array}$$
其中
$$\begin{array}{}\text{(3)}& a=\frac{C_{ϵ}}{\mathrm{\Delta }}b=\frac{2}{3}tr(\overline{\mathbf{S}})c=2\mathrm{\Delta }{C}_k(\mathbf{d}\mathbf{e}\mathbf{v}(\overline{\mathbf{S}}):\overline{\mathbf{S}}\sqrt{k_{sgs}}=\frac{-b+\sqrt{b^2+4ac}}{2a}\end{array}$$
当为不可压缩流体时$tr\overline{\mathbf{S}}=0$,那么$b=0$、$c=2\mathrm{\Delta }{C}_k(\mathbf{d}\mathbf{e}\mathbf{v}(\overline{\mathbf{S}}):\overline{\mathbf{S}})$,就可以得到$k_{sgs}=\frac{c}{a}=\frac{2C_k||\overline{\mathbf{S}}|{|}^2{\mathrm{\Delta }}^2}{C_{ϵ}}$

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

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

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

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