关于sonicLiquidFoam中pEqn.H的一个疑问
-
U = rAU*UEqn.H(); surfaceScalarField phid ( "phid", psi *( (fvc::interpolate(U) & mesh.Sf()) + rhorAUf*fvc::ddtCorr(rho, U, phi)/fvc::interpolate(rho) ) ); phi = (rhoO/psi)*phid; fvScalarMatrix pEqn//for solution of pEqn, construction with respect to p0 and rho0 ( fvm::ddt(psi, p) + fvc::div(phi)////?????? + fvm::div(phid, p)//=0, if psi=0 - fvm::laplacian(rhorAUf, p) );//
与rhoPimpleFoam相比,为什么多出了第二项 fvc::div(phi)////??????
-
是否与readThermodynamicProperties.H里面的
dimensionedScalar rhoO("rhoO", rho0 - psi*p0);有关? -
咦?看看这个:sonicLiquidFoam解析
-
@李东岳 :kiss:
论坛登录问题反馈可联系 li.dy@dyfluid.com