CFD中文网

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

@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];

@丑居 $M_j$ 是摩尔质量，单位：$\rm kg/mol$。

@尚善若水 看了一下，开发版和of2012差别很大，开发版里不能用specieThermo_了。依据要不要单位，组分的生成焓在开发版里可以这么写：

dimensionedScalar hi = chemistryModel_.thermo().hfi(speciei); //- Enthalpy of formation [J/kg]
或者
scalar hi = chemistryModel_.thermo().hfiValue(speciei); //- Enthalpy of formation [J/kg]

@尚善若水 没太看懂const scalar h[speciei] = specieThermos_h[speciei].hf(); //如何访问？这是啥意思。
如果需要得到每个组分的化学焓（生成焓），可以这么写：

const scalar hi = specieThermo_[speciei].Hc(); // 化学焓 J/kg

OpenFOAM版本不同，可能有些许差异，我用的是OpenFOAM-v2012。

@李东岳 哈哈，李老师，我不是大佬，我们组里没人用PINN研究化学反应。看论坛里讨论得这么火热，我也想学习学习PINN的，就是还不熟悉。搞这个要向李老师多学习呀！

@研究牲 粒子的相对速度与碰撞横截面积的乘积，单位是$\rm m^3/s$。源码里应该有说明的。

@研究牲 是的。