CFD中文网

想计算每个基元反应的反应速率， 需要调用chemistryModelI.H里定义的reactions()函数，

template<class CompType, class ThermoType>
inline const Foam::PtrList<Foam::Reaction<ThermoType> >&
Foam::chemistryModel<CompType, ThermoType>::reactions() const
{
    return reactions_;
}

怎么做才能调用到这个函数

@Wayne 我最近想把reactingFoam中的能量方程中加上质量扩散引起的焓输运项， Ji 是组分i 的扩散流量， hi 是每个组分的显焓，hi 在openfoam中可以得到吗？下面是reactingFoam的thermophsicalproperties中的代码：

thermoType
{
    type            hePsiThermo;
    mixture         reactingMixture;
    transport       sutherland;
    thermo          janaf;
    energy          sensibleEnthalpy;
    equationOfState perfectGas;
    specie          specie;
}

我在basicSpecieMixture.H文件中找到了计算单个组分的显焓的函数

 //- Sensible enthalpy [J/kg]
            virtual scalar Hs
            (
                const label speciei,
                const scalar p,
                const scalar T
            ) const = 0;

问题是const label speciei,是什么，我是要把组分的名字写上吗

这个我没有研究过，但是我有些链接是讲这个的，希望对你有用
http://xiaopingqiu.github.io/2016/03/20/fvOptions1/
http://www.sourceflux.de/blog/software-design-fvoptions/

我最近在用Pimple算法，是simple和piso的结合，simple的截止条件应该是在fvsolution中设置收敛残差，而piso不用设置收敛残差，直接给迭代次数，一般2-3次就可以了，本人也是刚接触of, 如果有错误请指出。

solve括号里的内容怎么理解，和fvSolution里的参数设置是否有关系？