[yanchong13]

I'm A PHD student major in combustion. Recently, I want to simulate a supersonic combustion case using OpenFoam. I thought rhoReactingFoam maybe capable for the case. But the position of the shock wave is ayways incorrect. I thought that may caused by th boundry condition of U or P or something else,but i cant figure is out. can anyon give me some suggestions?


p

/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  dev
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    location    "0";
    object      p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField   uniform 101325;

boundaryField
{
    inletfuel
    {
        type            fixedValue;
        value           uniform 101325;
    }
    inletair
    {
        type            fixedValue;
        value           uniform 101325;
    }
    outlet
    {
        type            zeroGradient;
    }

    frontAndBackPlanes
    {
        type            empty;
    }
    WALL_TOP
    {
        type            zeroGradient;
    }
    WALL_BOT
    {
        type            zeroGradient;
    }
    WALL_INSIDE
    {
        type            zeroGradient;
    }
}


// ************************************************************************* //

/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  dev
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volVectorField;
    location    "0";
    object      U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField   uniform (730 0 0);

boundaryField
{
    inletfuel
    {
        type            fixedValue;
        value           uniform (1200 0 0);
    }
    inletair
    {
        type            fixedValue;
        value           uniform (730 0 0);
    }
    outlet
    {
        type            inletOutlet;
        inletValue      uniform (730 0 0);
        value           uniform (730 0 0);
    }
    frontAndBackPlanes
    {
        type            empty;
    }
    WALL_TOP
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }
    WALL_BOT
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }
    WALL_INSIDE
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }
}


// ************************************************************************* //

/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  dev
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    location    "0";
    object      T;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 0 0 1 0 0 0];

internalField   uniform 1200;

boundaryField
{
    inletfuel
    {
        type            fixedValue;
        value           uniform 250;
    }
    inletair
    {
        type            fixedValue;
        value           uniform 1200;
    }
    outlet
    {
        type            inletOutlet;
        inletValue      uniform 340;
        value           uniform 340;
    }
    frontAndBackPlanes
    {
        type            empty;
    }
    WALL_TOP
    {
        type            zeroGradient;
    }
    WALL_BOT
    {
        type            zeroGradient;
    }
    WALL_INSIDE
    {
        type            zeroGradient;
    }
}


// ************************************************************************* //

fvSchemes

/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  dev
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

ddtSchemes
{
    default         Euler;
}

gradSchemes
{
    default         Gauss linear;
}

divSchemes
{
    default             none;

    div(phi,U)      Gauss linearUpwindV grad(U);
    div(phi,Yi)         Gauss linearUpwind;
    div(phi,h)          Gauss linearUpwind;
    div(phi,K)          Gauss linearUpwind grad(U);
    div(phid,p)         Gauss linearUpwind;
    div(phi,omega)      Gauss upwind;
    div(phi,epsilon)    Gauss upwind;
    div(phi,Yi_h)       Gauss upwind;
    div(phi,k)          Gauss upwind;
    div(((rho*nuEff)*dev2(T(grad(U)))))     Gauss linear;
    div(phiv,p)         Gauss linearUpwind grad(U);
}

laplacianSchemes
{
    default         Gauss linear corrected;
}    

interpolationSchemes
{
    default         linear;
    //interpolate(rho)              upwind phi;
    //interpolate((rho*(1|A(U))))   upwind phi;
    //flux(HbyA)                    upwind phi;
    //dotInterpolate(S,(rho_0*U_0)) linear;
}

snGradSchemes
{
    default         corrected;
}

wallDist
{
    method meshWave;
}

// ************************************************************************* //

fvSolution

/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  dev
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

solvers
{
    "rho.*"
    {
        solver          diagonal;
    }

    p
    {
        solver           PBiCGStab;
        preconditioner   DIC;
        tolerance        1e-6;
        relTol           1e-6;
    }

    pFinal
    {
        $p;
        relTol           0;
    }

    "(U|e|k|epsilon|omega)"
    {
        solver          PBiCGStab;
        preconditioner  DILU;
        tolerance       1e-6;
        relTol          1e-6;
    }

    "(U|e|k|epsilon|omega)Final"
    {
        $U;
    }

    "Yi.*"
    {
        solver          PBiCGStab;
        preconditioner  DILU;
        tolerance       1e-8;
        relTol          1e-6;
    }
}

PIMPLE
{
    momentumPredictor   yes;
    nOuterCorrectors    1;
    nCorrectors         2;
    nNonOrthogonalCorrectors 1;

    maxDeltaT           1e-7;
    maxCo               0.25;
    alphaTemp           0.05;
    alphaY              0.05;
    Yref
    {
        O2                  0.1;
        H2                  0.1;
    }
    rDeltaTSmoothingCoeff   0.05;
    rDeltaTDampingCoeff     1;
}

relaxationFactors
{
    equations
    {
        ".*"    1;
    }
}

// ************************************************************************* //

[as020002]

Hi,

I am doing combustion research on this geometry at the moment.
You can try the "waveTransmissive" boundary on the p.outlet

Regards,
Rick

[ozyz21]

hello. im simulating supersonic air hydrogen mixing in DLR Scramjet combustion too if you solve this problem or not can yol share case file? thank you.