运行结束后输出密度rho

[尚善若水]

这里改一下就好了，不然详细机理太占空间了

[欧阳]

@尚善若水 写个functionObject，后处理，这个我没有试过，我之前一直都是在调用模型的原代码上修改的。我去学学，
“按照cantera的处理，输出每步反应里面各个组分的放热和每步反应的总放热”
这个是我现在在实现的，每个化学反应步在每个网格的Qdot，然后一些敏感性ODCantera那边的已经实现了，但是这边一直编译的或多或多出现一些问题

[欧阳]

@尚善若水 大佬有输出过tc吗？我按```
tc[celli] = nReaction*cSum/tc[celli];

StandardChemistryModel.C中的修改出现的结果，tc的值特别的大。反应几乎不发生时，特征时间趋于无穷大。但是这些极大值没有实际意义，反而会造成数值问题。

[尚善若水]

@欧阳 先请教下，组分总放热有什么分析场景吗？我看到较多的是分步的。tc的话你是用的什么燃烧模型？PaSR吗？我以前输出过，我找下

[欧阳]

@尚善若水 喷雾燃烧，当关注一些燃烧的重要时间点，低温反应着火与高温反应着火时间点，还有燃烧衰退时，组分的总放热蛮重要的，尤其当分析出，这个物质的对流，扩散与化学反应哪个占主导后。是的，PaSR

[尚善若水]

@欧阳 我以前是这样输出的，10版本之后好像不行了，我没测试

[欧阳]

@尚善若水

                    IOobject
                    (
                        "Qdot." + Y_[i].name(),
                        this->mesh_.time().timeName(),
                        this->mesh_,
                        IOobject::NO_READ,
                        IOobject::NO_WRITE,  
                        IOobject::REGISTER    
                    ),

这样修改后，我使用

writeObjects
(
    "Qdot.OH" 
);

后，没有输出Qdot.OH

[尚善若水]

@欧阳 先存储一下试试？

 
//必须得在functions里面定义启用才会计算，然后才能缓存
cacheTemporaryObjects
(
    //RR:methaneReaction:CH4 //speciekReactioni
    //RR:methaneReaction //reactioni
 
    //Qdot:methaneReaction:CH4
    //Qdot:methaneReaction
);

[欧阳]

@尚善若水 我现在用的是v2306，测试出来的tc值，如果没有反应就特别大

[尚善若水]

@欧阳 我以前 OpenFOAM-8 用的上面发的代码，好像没问题的话