反应速率，摩尔分数问题

[CarmeloSun]

请问大家，oopenfoam里怎么输出多步反应的每一步的反应速率，还有怎么输出摩尔分数，谢谢啦

[李东岳]

最近在搞反应ODE看了一下代码，做个记录。

在chemistryModel.C中有定义：

forAll(RR_, fieldI)
    {
        RR_.set
        (
            fieldI,
            new DimensionedField<scalar, volMesh>
            (
                IOobject
                (
                    "RR." + Y_[fieldI].name(),
                    mesh.time().timeName(),
                    mesh,
                    IOobject::NO_READ,
                    IOobject::NO_WRITE
                ),
                mesh,
                dimensionedScalar("zero", dimMass/dimVolume/dimTime, 0.0)
            )
        );
    }

RR并不是体场。所以在求解器层面，可以定义一个volScalarField，然后把RR的值复制给这个volScalarField然后进行输出，并在paraview查看。

也可以吧上面的NO_WRITE改成AUTO_WRITE，就是不能再paraview查看。