CFD Online Discussion Forums - Supersonic compressible channel flow

Hello,

I am trying to reproduce the results of Coleman on a compressible channel flow with Mach number M=1.5. To do so, I am using rhoPimpleFoam on a channel flow configuration, using the same domain size than the authors. I triggered the turbulent state thanks to the utilities of wyldcat with that purpose.

In order to set the Mach number, I followed the instructions of the forward facing step from the OF tutorials. I am fixing the velocity using the meanVelocityForce fvOption. I am searchg for any option to fix the mass flux (as the authors) rather than the velocity, but for the moment it suits for me.

However, comparing with the results of Coleman (check attached Figure 1) I observe that my temperature and density fields are inversed: temperature should increase at the centreline, while the density should do the opposite.

Have you ever faced this channel flow before? My BC are:

U

Code:

boundaryField

{

    ".*Wall"

    {

        type    noSlip;

    }

    

    ".*Periodic"

    {

        type    $thirdDim; // cyclic for 3D simulation

    }

    

    ".*Cyclic"

    {

        type    cyclic;

    }

 }

Code:

boundaryField

{

    ".*Wall"

    {

        type    zeroGradient;

    }



    ".*Periodic"

    {

        type    $thirdDim;

    }



    ".*Cyclic"

    {

        type    cyclic;

    }

}

Code:

boundaryField

{



    ".*Wall"

    {

        type     fixedValue;

        value    uniform $Tw;  // equal to 1

    }



    ".*Periodic"

    {

        type    $thirdDim;

    }

    

    ".*Cyclic"

    {

        type    cyclic;

    }   

}

rho

Code:

boundaryField

{

    ".*Wall"

    {

        type    zeroGradient;

    }



    ".*Periodic"

    {

        type    $thirdDim;

    }

    

    ".*Cyclic"

    {

        type    cyclic;

    }    

}

Am I missing something?

Thanks in advance!