[MarkXI]

Hi everyone,


I am simulating a ideal compressible gas using dbnsTurbFoam (foam-extend 4.1) setting total pressure at the inlet. The pressure is not constant but it rises in time from the ambient value to a certain value therefore I have used timeVaryingTotalPressure as boundary condition. The solver evaluates the density at the inlet boundary using the ideal gas law





but it takes as p the value written in "p0" instead of using the actual pressure value, therefore I think that the density is not correctly calculated.


inlet
{
type timeVaryingTotalPressure;
rho none;
psi none;
gamma 1.1643;
p0 5e+06;
fileName "$FOAM_CASE/constant/p0vsTime";
outOfBounds clamp;
value uniform 5e+06;
}

I would very much appreciate it if anyone knows how to resolve this issue and is willing to share a solution. Thank you