openfoam化学反应模拟如何导出组分的反应源项

[江山不如凉竹]

我用openfoam进行化学反应模拟，如何在计算结束之后，导出组分的化学反应源项，比如RR_C7H16。如果不行，有什么办法可以在求解器中定义出组分的化学反应源项，运行再输出？

[尚善若水]

计算时应该可以直接输出，后处理获取没试过，可以使用writeNow计算一步输出。

[wangfei9088]

@江山不如凉竹 我就补充下求解器中保存组分的化学反应源项吧。
不知道你用的什么chemistryModel。假设你用的是StandardChemistryModel，那就将StandardChemistryModel.C文件中下面代码里的NO_WRITE改成AUTO_WRITE，编译完成后算例再运行一步保存就可以了。

     // Create the fields for the chemistry sources
     forAll(RR_, fieldi)
     {
         RR_.set
         (
             fieldi,
             new volScalarField::Internal
             (
                 IOobject
                 (
                     "RR." + Y_[fieldi].name(),
                     this->mesh().time().timeName(),
                     this->mesh(),
                     IOobject::NO_READ,
                     //IOobject::NO_WRITE
                     IOobject::AUTO_WRITE
                 ),
                 this->mesh(),
                 dimensionedScalar(dimMass/dimVolume/dimTime, Zero)
             )
         );
     }

注：要拷贝到自己的工程路径下修改编译，不建议直接在原始代码里改。