[join4us]

Hi. there.

I want to solve a hotpot simulation.

The pot has water and air in half on pot.

Base solver is interFoam.

I modified twoPhaseMixture.C and twoPhaseMixture.H for thermal expansion and thermal diffusive.
and add
TEqn.H
fvScalarMatrix TEqn
(
fvm::ddt(T)
+ fvm::div(phi, T)
- fvm::laplacian(DTEff, T)
);

UEqn.H
fvc::reconstruct
(
fvc::interpolate(interface.sigmaK())*fvc::snGrad(a lpha1)
// - ghf*fvc::snGrad(rho)
+ ghf*fvc::interpolate(rho)*twoPhaseProperties.TEf() *fvc::snGrad(T)
- fvc::snGrad(p_rgh)
) * mesh.magSf()

pEqn.H
surfaceScalarField phig
(
(
fvc::interpolate(interface.sigmaK())*fvc::snGrad(a lpha1)
// - ghf*fvc::snGrad(rho)
+ ghf*fvc::interpolate(rho)*twoPhaseProperties.TEf() *fvc::snGrad(T)
)*rAUf*mesh.magSf()
);

The simulation gives good result except pressure field.
p=p_rgh + rho*gh

I attached flow field pictures.

Why pressure field look like step?

Could you help me to solve the problem well?

[Bernhard]

Because rho is discontinuous.