CFD中文网

@先先贝 那瞧不出来啥问题，得慢慢debug了。

@先先贝 看上去是#4 sqrt 开方函数有问题，是不是出现负数了？

@李东岳 CFD中文网陪伴了我整个读博生涯，感谢给我答疑解惑的大佬们。感谢李老师！

@ann ECN网站现在可以进去了。

@钟余 如果 $\nu$ 和 $\varepsilon$ 边界条件都没有设置为 0，那就要试试一个核跑去慢慢 debug 了。

@钟余 边界条件里湍动能 $k$ 给一个很小的值，不要设置为0。EDC模型里需要把 $k$ 作为分母。

@李东岳 哈哈，不是大佬。向东岳老师学习！

@wangfei9088 追踪了cpc的文章，看到名字，第一想到了你

@尚善若水 哈哈

@wangfei9088 再多问大佬一句，代码是open source的吗？我看到了CPC Library link to program files。我们用的求解器是基于OpenFOAM 12的内部求解器和v2312俄罗斯的一个求解器，都没考虑混合平均多组分扩散。感觉这对算缓燃还是挺重要的。

@尚善若水 是open source的，校对完正式发布就可以看到完整代码，提供了三个版本，可以满足大多数用户的需求。

@尚善若水 哈哈，是的，刚online的。消息这么灵通啊

@BlancheDH 上面的文件有好几个打不开。可以仿照这个帖子的#15楼附件代码写，用的就是OpenFOAM-10。

https://cfd-china.com/topic/6744/wray-agarwal湍流模型/15?_=1732935616393

@zju_klh cSpace和cTime需要用户给定，对应文献 (https://doi.org/10.1016/j.compfluid.2015.01.018) 公式里的$C_1$和$C_2$。

$$
\begin{equation}
P_{ij}
=\left(\frac {r_i + r_j} {\max (r_i+r_j, \mu_0)} \right)^{C_1}
\exp \left( -C_2 \frac{|\alpha_0 - \beta_0|}{\Delta t} \right)
\end{equation}
$$
其中，模型常数 $C_1 = 1$, $ C_2 = 0-0.3$。

一般就取1和0.3。

trajectoryCoeffs
{
    cSpace       1;
    cTime        0.3;
    coalescence  on;
}

@李东岳 没有下载权限，看不了

@yjc95 在 OpenFoam不可压+瞬态 对流项隐性离散推导看不懂 中说：

为啥是二重积分？

分别对体积 $V$ 和时间 $t$ 积分。

@yjc95 在 OpenFoam不可压+瞬态 对流项隐性离散推导看不懂 中说：

等式左边是如何转化到右边的呢？

体积分转为面积分，用的高数里的高斯定理。李老师的《笔记》里有推导的。

@HIT_lly 试试这个：

mpirun --use-hwthread-cpus -np 32 compressibleInterIsoFoam -parallel

@tiankai U 没问题，试试将fields后面的大括号{}改成小括号()，最后面再加一个分号;，如下。

    fields
    //{
    (
        U
        {
            mean       on;      
            prime2Mean off;    
            base       time;     
        }
    //}
    );

@13615617682 我不清楚fluent里具体用的单位，可能得查手册。但我见过的文献指前因子A的单位都是 $\rm m^3/(mol \cdot s)$。

@13615617682 啊，我不知道fluent中的指前因子是什么样的，它的单位是什么？能看一下不？

@13615617682 看 $A_i$ 的单位 $\rm kg/(m^2sPa)$，反应速率 $R_i$（单位：kg/s）应该是这样的：
$$
R_i =A_i \exp \left (-\frac{E_i}{RT} \right) S_i P
$$
其中，$S_i$是面积（$\rm m^2$），$P$是压力（$\rm Pa$）。

@hehaohao2013 报错提示，湍流模型只能选择kEpsilon或kOmegaSST，没有kineticTheory这个模型。

@尚善若水 应该可以吧。RR[speciei][celli]是celli网格里speciei的反应速率，那该网格里该组分的Qdot应该可以这么写：

scalar Qdot-speciei-celli = hi * RR_[speciei][celli];