[OpenFoamBegins]

Hello everyone,

I'm trying to simulate the interaction of a jet (secondary flow) in a turbulent mainstream flow using the buoyantSimpleFoam solver. The hot turbulent freestream (air at T=540K) is flowing over a flat plate (wall) and the cold secondary flow (air at T=300K) is injected verticaly from a plenium then a pipe (both walls) inclined at 35° compared to the plate (picture of the domain given in the file attached).

I've been able to run the simpleFoam solver on this domain, however i would like now to simulate the adiabatic temperature over the plate (with buoyantSimpleFoam), but I don't know how to set up the pressure/p_rgh initial conditions (for the others 0 files I took it from my simpleFoam simulation), and my simulation collapses after 6 or 7 iterations (it crashes at the p_rgh calculation step).

From what I've read from the OpenFoam tutos, so far I've tried to set the pressure with the "calculated" BC for all boundaries, and for my p_rgh initial BC i've put "fixedFluxPressure" (with zero gradient and 1atm value) for the walls and the two inlets, and a fixed value at the outlet (1atm as well).

My p and p_rgh 0 files are given below ($pressure is 101325):

FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    location    "0";
    object      p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

#include        "include/initialConditions"

internalField   uniform $pressure;

boundaryField
{

    Blade
    {
        type            calculated;
        value		 uniform $pressure;
    }

    Pipe
    {
        type            calculated;
        value		 uniform $pressure;
    }

    CoolantInlet
    {
        type            calculated;
        value		 uniform $pressure;
    }
    
    Plenium
    {
        type            calculated;
        value		 uniform $pressure;
    }
    
    Inlet
    {
        type            calculated;
        value		 uniform $pressure;
    }
    
    Outlet
    {
        type            fixedValue;
        value		 uniform $pressure;
    }
    
    TopBoundary
    {
        type            slip;
    }
    
    FrontAndBack
    {
        type            symmetry;
    }

}

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

FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    location    "0";
    object      p_rgh;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

#include        "include/initialConditions"

internalField   uniform $pressure;

boundaryField
{

    Blade
    {
        type            fixedFluxPressure;
        gradient        uniform 0;
        value           uniform $pressure;
    }

    Pipe
    {
        type            fixedFluxPressure;
        gradient        uniform 0;
        value           uniform $pressure;
    }

    CoolantInlet
    {
        type            fixedFluxPressure;
        gradient        uniform 0;
        value           uniform $pressure;
    }
    
    Plenium
    {
        type            fixedFluxPressure;
        gradient        uniform 0;
        value           uniform $pressure;
    }
    
    Inlet
    {
        type            fixedFluxPressure;
        gradient        uniform 0;
        value           uniform $pressure;
    }
    
    Outlet
    {
        type            fixedValue;
        value           uniform $pressure;
    }
    
    TopBoundary
    {
        type            slip;
    }
    
    FrontAndBack
    {
        type            symmetry;
    }

}

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

A picture of the domain and of the p_rgh field at the latest time before the simulation crashes are attached. If someone could help me on this I would be really grateful !

Thank you