压力泊松方程植入问题

[strive_DUT]

本人在修改求解器过程中需要植入一个新的压力泊松方程，形式如下所示，方程来自文献。试了多次，但都没有成功，感觉是方程植入形式的问题，求助各位大佬指点，该方程应该怎么植入，谢谢啦

[李东岳]

@strive_dut

    surfaceScalarField phiHbyA = mesh.Sf() & fvc::interpolate(HbyA);
    
    fvc::ddt(rho)
  + psi*correction(fvm::ddt(p))
  + fvc::div(rho, phiHbyA)
  - fvm::laplacian(rhobyAf, p)

[strive_DUT]

@李东岳 感谢李老师百忙之中指点！这么来看correction(fvm::ddt(p))就代表前后时间步压力差值的时间变化率啦