[J_here]

Hi All,

I am trying to simulate a vertical heated plate on a column surrounded by atmospheric boundary conditions (see attached image) and these unexplained velocities develop either at the bottom left corner or at the intersection of the column and the bottom on the heated surfaces side.

My 0 files are as follows:

alphat

 dimensions      [0 2 -1 0 0 0 0];  //m^2/s

internalField   uniform 0;

boundaryField
{
    heatedSurface
    {
        type            alphatJayatillekeWallFunction;
        Prt             1;
        value           uniform 0;
    }

    fixedWalls
    {
        type            alphatJayatillekeWallFunction;
        Prt             1;
        value           uniform 0;
    }

    openAir
    {
        type            calculated;
	value		uniform 0;
    }
    
    frontAndBack
    {
        type            empty;
    }

epsilon

dimensions      [0 2 -3 0 0 0 0];

internalField   uniform 0.01;

boundaryField
{
    heatedSurface
    {
        type            epsilonWallFunction;
        value           uniform 0.01;
    }

    fixedWalls
    {
        type            epsilonWallFunction;
        value           uniform 0.01;
    }

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

    frontAndBack
    {
        type            empty;
    }
}

k

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

internalField   uniform 0.1;

boundaryField
{
    heatedSurface
    {
        type            kqRWallFunction;
        value           uniform 0.1;
    }

    fixedWalls
    {
        type            kqRWallFunction;
        value           uniform 0.1;
    }

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

    frontAndBack
    {
        type            empty;
    }
}

nut

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

internalField   uniform 0;

boundaryField
{
    heatedSurface
    {
        type            nutkWallFunction;
        value           uniform 0;
    }

    fixedWalls
    {
        type            nutkWallFunction;
        value           uniform 0;
    }

    openAir
    {
        type            calculated;
        value           uniform 0;
    }


    frontAndBack
    {
        type            empty;
    }
}

p

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

internalField   uniform 0;

boundaryField
{
    heatedSurface
    {
        type            calculated;
        value           $internalField;
    }

    fixedWalls
    {
        type            calculated;
        value           $internalField;
    }

    openAir
    {
        type            calculated;
        value           $internalField;
    }

    frontAndBack
    {
        type            empty;
    }
}

p_rgh

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

internalField   uniform 0;

boundaryField
{
    heatedSurface
    {
        type            fixedFluxPressure;
        rho             rhok;
        value           uniform 0;
    }

    fixedWalls
    {
        type            fixedFluxPressure;
        rho             rhok;
        value           uniform 0;
    }

    openAir
    {
        type            fixedFluxPressure;
        rho             rhok;
        value           uniform 0;
    }


    frontAndBack
    {
        type            empty;
    }
}

T

dimensions      [0 0 0 1 0 0 0];

internalField   uniform 300;

boundaryField
{
    heatedSurface
    {
        type            fixedValue;
        value           uniform 305;
    }

    fixedWalls
    {
        type           zeroGradient;
    }

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

    frontAndBack
    {
        type            empty;
    }
}

and finally U

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

internalField   uniform (0 0 0);

boundaryField
{
    heatedSurface
    {
        type            noSlip;
    }

    openAir
    {
        type            pressureInletOutletVelocity;
        outletValue     uniform (0 0 0);
        value           uniform (0 0 0);
    }


    fixedWalls
    {
	  type           noSlip;

    }

    frontAndBack
    {
        type            empty;
    }
}

I have looked elsewhere on the forum but nothing specifically deals with this. Or if it is similar like 'atmospheric boundary conditions for hotRoom' it is with an old version of openFoam.

Thanks in advance for any help or guidance you can spare