[Tegdlemat]

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

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

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

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;
    }

[Sakun]

Hi,


Did you manage to sort this out ?