divq中的correction是怎么来的？

[BlookCFD]

查看OpenFOAM10的代码可以发现，在src/ThermophysicalTransportModels/turbulence/eddyDiffusivity/eddyDiffusivity.C中的divq的实现为：

-correction(fvm::laplacian(this->alpha()*this->alphaEff(), he))
-fvc::laplacian(this->alpha()*this->kappaEff(), this->thermo().T());

而按我的理解，这里应该只是求解 $-\mathrm{\nabla }(k_{eff}\mathrm{\nabla }T)$，即：

-fvc::laplacian(this->alpha()*this->kappaEff(), this->thermo().T())

那这里多出来的correction应该怎么理解？

[李东岳]

《无痛笔记》5.5.13 fvMatrix 的correction() 函数

http://www.dyfluid.com/theory.pdf

[BlookCFD]

@李东岳 看得懵懵懂懂，还是理不出来相应的数学表达式。感觉是一种线化处理策略，如果没有correction，就变成了纯显式处理

[李东岳]

就是同一个时间步下的${\mathrm{\nabla }}^2{q}^{n+1}-{\mathrm{\nabla }}^2{q}^n$，在收敛的时候为0，具体为什么会有这一项，还需要研究下

[BlookCFD]

@李东岳 多谢大佬点拨

[09116333]

@李东岳 《笔记》里没有5.5.13这一节啊？