我在无痛苦ns方程笔记里面增加了对这一行代码的解释:
fvScalarMatrix alpha1Eqn ( fvm::ddt(alpha1) - fvc::ddt(alpha1) // WHY??? - fvm::laplacian(alpha1alpha2f()*pPrimeByA_(), alpha1) ); alpha1Eqn.relax(); alpha1Eqn.solve();尤其是其中的
fvm::ddt(alpha1) - fvc::ddt(alpha1) // WHY???OpenFOAM交流区