buoyantSimpleFoam divergence
 

Hi,

I'm trying to model a buoyant thermal plume from a chimney in a cross flow of air using buoyantSimpleFoam.

My rho min/max values quickly degenerate into nonsense and I'm sure it must be a boundary condition issue. I've trawled the links and tried a lot of different boundary conditions without success.

The relevant b.c.'s etc are below.

I'd be very grateful for any advice/suggestions.

Thanks,

Tom

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

dimensions [0 0 0 1 0 0 0];

internalField uniform 293;

boundaryField
{
pressureOutlet
{
type inletOutlet;
inletValue uniform 293;
value uniform 293;
}
velocityInlet
{
type fixedValue;
value uniform 293;
}
atmosphere
{
type zeroGradient;
}
sea
{
type zeroGradient;
}
funnel_patch0
{
type fixedValue;
value uniform 476;
}
ship_patch0
{
type zeroGradient;
}
}

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

FoamFile
{
version 2.0;
format ascii;
class volVectorField;
object U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField uniform (0 15.88 0);

boundaryField
{

velocityInlet
{
type fixedValue;
value uniform (0 15.88 0);
}

pressureOutlet
{
type inletOutlet;
inletValue 0;
value 0;
}

sea
{
type fixedValue;
value uniform (0 0 0);
}

funnel_patch0
{
type fixedValue;
value uniform (-2.4248 0 17.3017);
}

ship_patch0
{
type fixedValue;
value uniform (0 0 0);
}

atmosphere
{
type zeroGradient;
}

}

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

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

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

internalField uniform 101325.;

boundaryField
{
pressureOutlet
{
type calculated;
value uniform 101325;
}
velocityInlet
{
// type zeroGradient;

type calculated;
value uniform 101325;
}
atmosphere
{
type calculated;
value uniform 101325;
}
sea
{
type calculated;
value uniform 101325;
}
funnel_patch0
{
type calculated;
value uniform 101325;
}
ship_patch0
{
type calculated;
value uniform 101325;
}
}

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

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

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

internalField uniform 0;

boundaryField
{

pressureOutlet
{
type buoyantPressure;
gradient uniform 0;
value uniform 101325;
}
velocityInlet
{
type buoyantPressure;
gradient uniform 0;
value uniform 101325;
}
atmosphere
{
type buoyantPressure;
gradient uniform 0;
value uniform 101325;
}
sea
{
type wallBuoyantPressure;
gradient uniform 0;
value uniform 101325;
}
funnel_patch0
{
type buoyantPressure;
gradient uniform 0;
value uniform 101325;
}
ship_patch0
{
type wallBuoyantPressure;
gradient uniform 0;
value uniform 101325;
}
}

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

FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "constant";
object thermophysicalProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

thermoType hPsiThermo<pureMixture<constTransport<specieThermo <hConstThermo<perfectGas>>>>>;

pRef 101325;

mixture
{
specie
{
nMoles 1;
molWeight 28.96;
}
thermodynamics
{
Cp 1004.4;
Hf 0;
}
transport
{
mu 1.831e-05;
Pr 0.705;
}
}


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