pEqn.H

Go to the documentation of this file.

 
{
    volScalarField rAU(permeability/mu); //K/mu
    surfaceScalarField rhorAUf("rhorAUf", fvc::interpolate(rho*rAU));//rho/A on surface
    volVectorField HbyA(U*0);
 
    surfaceScalarField phig("phig",(fvc::interpolate(rho)*rhorAUf * g) & mesh.Sf());
    
    surfaceScalarField phiHbyA
    (
        "phiHbyA",
        phig
    );
    // Update the pressure BCs to ensure flux consistency
    constrainPressure(p, rho, U, phiHbyA, rhorAUf);
 
    fvScalarMatrix p_rghDDtEqn
    (
        porosity*rho*betaT*fvm::ddt(p)
        -porosity*rho*alphaP*fvc::ddt(T)
        +fvc::div(phiHbyA)
    );
    while (pimple.correctNonOrthogonal())
    {
        fvScalarMatrix pEqn
        (
            p_rghDDtEqn - fvm::laplacian(rhorAUf, p)
        );
        pEqn.solve();
        if (pimple.finalNonOrthogonalIter())
        {
            // option 2: using flux reconstruct velocity, magic function of reconstruct
            phi = phiHbyA + pEqn.flux();
            U = HbyA + rAU*fvc::reconstruct((phig + pEqn.flux())/rhorAUf);
            U.correctBoundaryConditions();
        }
    }
}