[Tegdlemat]

Hi guys,

I am currently working on a case simulating a quiescent flow using rhoCentralFoam, with the velocity of the flow equals to 0 at the start (act like a detonation in an actual experiment), across a compression corner. The case is running fine, however the flow is separated at the leading edge from the beginning of the calculation.

I hypothetically think that the flow should be same as a normal case run as the first picture, yet the flow with quiescent condition is show in the second picture. (See attached files)

I suspect the reason is the 0 velocity of air flow plus the no-slip condition of the wall, which made the flow separated at the leading edge. However, I am not sure if it makes sense. Therefore I would like to seek help from u guys, since the case is a bit specific and I can find nothing to take reference from.

Below are the 0/ codes:

0/U

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

internalField   uniform (0 0 0);

boundaryField
{
    inlet
    {
        type            fixedValue;
        value           uniform (539.5 0 0);
    }

    outlet
    {
        type            zeroGradient;	
    }

    top
    {
        type            supersonicFreestream;
        pInf            922.4;
        TInf            143.1;
        UInf           (539.5 0 0);
        gamma           1.4;
        value           uniform (539.5 0 0);
    }

    wall
    {
         type            noSlip;
    }

    sides
    {
        type            empty;
    }
}

0/p

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

internalField   uniform 922.4;

boundaryField
{
    inlet
    {
        type           fixedValue;
        value          uniform 922.4;
    } 

    outlet
    {
        type            waveTransmissive;
        field           p;
        psi             thermo:psi;
        gamma           1.4;
        fieldInf       922.4;
        lInf            1;
        value          uniform 922.4;
    }

    top
    {
        type           zeroGradient;
    }

    wall
    {
        type            zeroGradient;
    }

    sides
    {
        type            empty;
    }

}

0/T

dimensions      [0 0 0 1 0 0 0];

internalField   uniform 143.1;

boundaryField
{
    inlet
    {
        type            fixedValue;
        value           $internalField;
    }

    outlet
    {
        type            inletOutlet;
        value           $internalField;
        inletValue      $internalField;
    }

    top
    {
        type            inletOutlet;
        value           $internalField;
        inletValue      $internalField;
    }

    wall
    {
        type            fixedValue;
        value           uniform 300;
    }

    sides
    {
        type            empty;
    }
}