[wayne14]

Hi everyone,

I already have a good understanding of buoyantSimpleFoam, but in the pEqn.H of buoyantPimpleFoam, the following code makes it difficult to understand the logic of the solver:

// Thermodynamic density needs to be updated by psi*d(p) after the
    // pressure solution - done in 2 parts. Part 1:
    thermo.rho() -= psi*p_rgh;

    volScalarField rAU(1.0/UEqn.A());
    surfaceScalarField rAUf("rAUf", fvc::interpolate(rho*rAU));

    volVectorField HbyA("HbyA", U);
    HbyA = rAU*UEqn.H();

    surfaceScalarField phig(-rAUf*ghf*fvc::snGrad(rho)*mesh.magSf());

Why we need thermo.rho() -= psi*p_rgh ?

Many thanks,
Yan Wang

[wayne14]

Originally Posted by wayne14

Hi everyone,

I already have a good understanding of buoyantSimpleFoam, but in the pEqn.H of buoyantPimpleFoam, the following code makes it difficult to understand the logic of the solver:

// Thermodynamic density needs to be updated by psi*d(p) after the
    // pressure solution - done in 2 parts. Part 1:
    thermo.rho() -= psi*p_rgh;

    volScalarField rAU(1.0/UEqn.A());
    surfaceScalarField rAUf("rAUf", fvc::interpolate(rho*rAU));

    volVectorField HbyA("HbyA", U);
    HbyA = rAU*UEqn.H();

    surfaceScalarField phig(-rAUf*ghf*fvc::snGrad(rho)*mesh.magSf());

Why we need thermo.rho() -= psi*p_rgh ?

Many thanks,
Yan Wang

It seems that I have solved the question.