[agustinvo]

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

boundaryField
{
    ".*Wall"
    {
        type    noSlip;
    }
    
    ".*Periodic"
    {
        type    $thirdDim; // cyclic for 3D simulation
    }
    
    ".*Cyclic"
    {
        type    cyclic;
    }
 }

p

boundaryField
{
    ".*Wall"
    {
        type    zeroGradient;
    }

    ".*Periodic"
    {
        type    $thirdDim;
    }

    ".*Cyclic"
    {
        type    cyclic;
    }
}

T

boundaryField
{

    ".*Wall"
    {
        type     fixedValue;
        value    uniform $Tw;  // equal to 1
    }

    ".*Periodic"
    {
        type    $thirdDim;
    }
    
    ".*Cyclic"
    {
        type    cyclic;
    }   
}

rho

boundaryField
{
    ".*Wall"
    {
        type    zeroGradient;
    }

    ".*Periodic"
    {
        type    $thirdDim;
    }
    
    ".*Cyclic"
    {
        type    cyclic;
    }    
}

Am I missing something?


Thanks in advance!

[dlahaye]

I will be happy to try the same case by replacing the cyclic boundary conditions with inlet and outlet boundary conditions as in tutorials of rhoSimpleFoam and rhoPimpleFoam.

What are dimension of the channel and velocity at the inlet?

Best, Domenico.

[agustinvo]

You can check the reference I'm following:
A numerical study of turbulent supersonic isothermal-wall channel flow

[Ale22pic]

Hi Agustin, I know it's been a while since you wrote it, but I'm trying to replicate Coleman results with my own solver. Is it possible to have the case or at least give me some info in order to understand where start? Thank you in advance.

Kind regard,
Alessio