[Butcher]

Below are the boundary condition I have given to a buoyantBoussinesqSimpleFoam problem. I am getting the following error -
Continuity error cannot be removed by adjusting the outflow.
Please check the velocity boundary conditions and/or run potentialFoam to initialise the outflow.
Total flux : 3.52534
Specified mass inflow : 0.03
Specified mass outflow : 0
Adjustable mass outflow : 0

Please check and see where the error is?

FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object alphat; //turbulent thermal diffusivity
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField uniform 0;

boundaryField
{
Inlet
{
type zeroGradient;
}


outlet
{
type zeroGradient;
}


cylwall
{
type alphatJayatillekeWallFunction;
Prt 0.85;
value uniform 0;
}
sqwall
{
type alphatJayatillekeWallFunction;
Prt 0.85;
value uniform 0;
}
}


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

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

internalField uniform 0.3;

boundaryField
{
Inlet
{
type turbulentMixingLengthDissipationRateInlet;
mixingLength 0.0035; //a good estimation is 0.07*characteristic_length
value uniform 0.3;
}

outlet
{
type zeroGradient;
}

cylwall
{
type epsilonWallFunction;
value $internalField;
}
sqwall
{
type epsilonWallFunction;
value $internalField;
}

}

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

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

internalField uniform 0.034;

boundaryField
{
Inlet
{
type turbulentIntensityKineticEnergyInlet;
intensity 0.0505;//turbulence intensity
value uniform 0.034;
}

outlet
{
type zeroGradient;
}

cylwall
{
type kqRWallFunction;
value $internalField;
}

sqwall
{
type kqRWallFunction;
value $internalField;
}

}

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

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

internalField uniform 0;

boundaryField
{
Inlet
{
type calculated;
value uniform 0;
}

outlet
{
type calculated;
value uniform 0;
}

cylwall
{
type nutkWallFunction;
value uniform 0;
}
sqwall
{
type nutkWallFunction;
value uniform 0;
}

}

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

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

internalField uniform 1e5;

boundaryField
{
Inlet
{
type fixedFluxPressure;
value uniform 1e5;
}

outlet
{
type fixedValue;
value uniform 1e5;
}

cylwall
{
type fixedFluxPressure;
value uniform 1e5;
}

sqwall
{
type fixedFluxPressure;
value uniform 1e5;
}
}

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

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

internalField uniform 0;

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

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

cylwall
{
type fixedFluxPressure;
rho rhok;
value uniform 0;
}
sqwall
{
type fixedFluxPressure;
rho rhok;
value uniform 0;
}
}

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

dimensions [0 0 0 1 0 0 0];

internalField uniform 300;

boundaryField
{
Inlet
{
type zeroGradient;
}
outlet
{
type zeroGradient;
}
sqwall
{
type zeroGradient;
}

cylwall
{
type fixedValue;
value uniform 500;
}
}

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

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

internalField uniform (0 0 0);

boundaryField
{
Inlet
{
type fixedValue;
value uniform (0 3.0 0);
}

outlet
{
type zeroGradient;
//value uniform (0 0 0);
}

cylwall
{
type noSlip;
}

sqwall
{
type noSlip;
}
}

[Yann]

Hello,

I'm not sure this is the only problem in your setup, but with solvers using p_rgh variable, you need to define pressure boundary conditions on p_rgh and use the calculated boundary condition on p. Have a look at the tutorials to have examples.

Also your message is very hard to read, please use [CODE] blocks to post files contents.

Regards,
Yann

[Butcher]

Originally Posted by Yann

Hello,

I'm not sure this is the only problem in your setup, but with solvers using p_rgh variable, you need to define pressure boundary conditions on p_rgh and use the calculated boundary condition on p. Have a look at the tutorials to have examples.

Also you're message is very hard to read, please use [CODE] blocks to post files contents.

Regards,
Yann

Thanks! I will try this out.