[Toorop]

Hi,

I would like to specify atmosphere boundary condition in the hootRoom tutorials. The default case uses walls on the sides and at the top, so the only modification is the change in boundary conditions.

alphat
{
        type            calculated;
        value           uniform 0;
}

epsilon
{
        type            inletOutlet;
        inletValue      uniform 0.01;
        value           uniform 0.01;
}

k
{
        type            inletOutlet;
        inletValue      uniform 0.1;
        value           uniform 0.1;
}

kappat
{
        type            calculated;
        value           uniform 0;
}

nut
{
        type            calculated;
        value           uniform 0;
}

p
{
        type            calculated;
        value           0;
}

p_rgh
{
        type           totalPressure;
        p0              uniform 0;
        U               U;
        phi             phi;
        rho             rhok;
        psi             none;
        gamma      1;
        value         uniform 0;
}

T
{
        type            inletOutlet;
        inletValue      uniform 300;
        value           uniform 300;
}

U
{
        type            pressureInletOutletVelocity;
        value           uniform (0 0 0);
}

Unfortunately, the BC that I set is not giving me the results I'm looking for. The air won't exit the domain in the upward direction as expected, it makes an U shape path (picture attached). The solvers were buoyantBoussinesqPimpleFoam and buoyantBoussinesqSimpleFoam.

It would be great if someone could share his setup to achieve correct behaviour at the boundary. Thanks!