CFD Online Discussion Forums - totalPressure Error

Hi everyone,

I am currently attempting to simulate a shock train phenomenon as shown in the first picture with rhoCentralFoam. In my case, I require a inlet pressure at 47000pa and outlet pressure at 104000pa, with inlet Mach number at 1.61 and temperature at 197K. I am using the boundary conditions below for the simulation, and the pressure at the inlet matches the target pressure magnitude, but not at the outlet. I wonder if I implied the totalpressure wrongly that causes the pressure there basically equals to the total pressure I defined. I have included the pressure contour as the second picture. I would be very grateful if anyone could lend my a helping hand.

0/p

Code:

dimensions      [1 -1 -2 0 0 0 0];



internalField     uniform 47984;



boundaryField

{

    Inlet

    {   

        type           totalPressure;

        p0              uniform 207000;

        gamma      1.4;

        psi              thermo:psi;

        value          uniform 47984;

    } 



    Outlet

    {

        type           totalPressure;

        p0              uniform 451260;

        gamma      1.4;

        psi              thermo:psi;

        value          uniform 104859;

    }



    Wall

    {

        type            zeroGradient;

    }



    Sides

    {

        type            empty;

    }

0/T

Code:

dimensions      [0 0 0 1 0 0 0];



internalField   uniform 197.6;



boundaryField

{

    Inlet

    {

       type            fixedValue;

       value          uniform 197.6;

    }



    Outlet

    {

        type            zeroGradient;

    }



    Wall

    {

        type            zeroGradient;

    }



    Sides

    {

        type            empty;

    }

0/U

Code:

dimensions      [0 1 -1 0 0 0 0];



internalField   uniform (453.65 0 0);



boundaryField

{

    Inlet

    {

        type            pressureInletOutletVelocity;

        value           $internalField;

    }



    Outlet

    {

        type            pressureInletOutletVelocity;

        value           $internalField;

    }



    Wall

    {

         type            noSlip;

    }



    Sides

    {

        type            empty;

    }