[karoltomek]

Hello

I use buoyantSimpleFoam solver. I have a box 1x1x1.4m with inlet, outlet and two heating boxes inside.
My BC

U

dimensions      [0 1 -1 0 0 0 0];
internalField   uniform (0 0 0);
boundaryField
{
    wall
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }
    inlet
    {
        type            fixedValue;
        value           uniform (0 0 1);
    }
    outlet
    {
        type            inletOutlet;
        inletValue      uniform (0 0 0);
        value           uniform (0 0 0);
    }
    grzalka_1
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }
    grzalka_2
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }
}

T

dimensions      [ 0 0 0 1 0 0 0 ];
internalField   uniform 293;
boundaryField
{
    wall
    {
        type            zeroGradient;
    }
    inlet
    {
        type            fixedValue;
        value           uniform 293;
    }
    outlet
    {
        type            zeroGradient;
    }
    grzalka_1
    {
        type            compressible::turbulentHeatFluxTemperature;
        heatSource      power;
        q               uniform 0;
        kappa           fluidThermo;
        Qr              none;
        kappaName       none;
        value           uniform 293;
    }
    grzalka_2
    {
        type            compressible::turbulentHeatFluxTemperature;
        heatSource      power;
        q               uniform 55;
        kappa           fluidThermo;
        Qr              none;
        kappaName       none;
        value           uniform 293;
     }
}
 
// ************************************************************************* //

p_rgh

 
dimensions      [ 1 -1 -2 0 0 0 0 ];
internalField   uniform 101325;
boundaryField
{
    wall
    {
        type            zeroGradient;
    }
    inlet
    {
        type            zeroGradient;
    }
    outlet
    {
        type            fixedValue;
        value           uniform 101325;
    }
    grzalka_1
    {
        type            zeroGradient;
    }
    grzalka_2
    {
        type            zeroGradient;
    }
}
 
// ************************************************************************* //

When buoyantSimpleFoam solves this case, there is a additional massFlow (kg/h) from heatting box (grzalka_2) >>> see my attachments.

But if I set g=0 [m/s^2] additional massflow disappear (!)

I have the same problem with buoyantBoussinesqSimpleFoam (addtional mass flow and very high temperatures).

My simulation crashes after 1800 steps. I have this problem with different meshes.

It is a bug? Anyone have similar problem?

[karoltomek]

Hello

I know what was wrong.

The correct set for p_rgh for all walls (grzalka_1, grzalka_2 and wall) is fixedFluxPressure, not zeroGradient.