[raghav.venky]

I'm trying to simulate a forced convection problem using the buoyantBoussinesqSimpleFoam solver.

When I run the wallHeatFlux utility, I get this error:

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0

Time = 0
Selecting thermodynamics package 
{
    type            heRhoThermo;
    mixture         pureMixture;
    transport       const;
    thermo          hConst;
    equationOfState perfectGas;
    specie          specie;
    energy          sensibleEnthalpy;
}

Reading/calculating face flux field phi

Selecting turbulence model type RASModel
Selecting RAS turbulence model kEpsilon
#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2   in "/lib/x86_64-linux-gnu/libc.so.6"
#3  Foam::compressible::mutkWallFunctionFvPatchScalarField::calcMut() const at ??:?
#4  Foam::compressible::mutWallFunctionFvPatchScalarField::updateCoeffs() at ??:?
#5  Foam::fvPatchField<double>::evaluate(Foam::UPstream::commsTypes) at ??:?
#6  Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::GeometricBoundaryField::evaluate() at ??:?
#7  Foam::compressible::RASModels::kEpsilon::kEpsilon(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::fluidThermo const&, Foam::word const&, Foam::word const&) at ??:?
#8  Foam::compressible::RASModel::adddictionaryConstructorToTable<Foam::compressible::RASModels::kEpsilon>::New(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::fluidThermo const&, Foam::word const&) at ??:?
#9  Foam::compressible::RASModel::New(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::fluidThermo const&, Foam::word const&) at ??:?
#10  Foam::compressible::turbulenceModel::addturbulenceModelConstructorToTable<Foam::compressible::RASModel>::NewturbulenceModel(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::fluidThermo const&, Foam::word const&) at ??:?
#11  Foam::compressible::turbulenceModel::New(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::fluidThermo const&, Foam::word const&) at ??:?
#12  
 at ??:?
#13  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#14  
 at ??:?
Floating point exception (core dumped)

the simulationType in the turbulenceProperties file is RASModel. If i change it to laminar, the utility works and calculates the fluxes on the walls. What could be the problem?
contents of the constant folder attached.

Thanks

[jherb]

I am confused: You use buoyantBoussinesqSimpleFoam but with compressible turbulence model/wall function? Have a look at the tutorials for this solver. If I remember correctly it only uses an expansion factor to calculate the density change by temperature but not a full thermophysical model.

[Muhammad Omer Mughal]

Dear Raghav

I am having same error. Can you please advise how you solved it. I shall be grateful.

Thanks

[BSengupta]

Hi Can any one please help me to overcome it.. Thank you very much.



Overall domain bounding box (0 0 0) (1 0.4 0.1)





Error:
Time = 10


GAMG: Solving for Ux, Initial residual = 0.424802, Final residual = 0.017996, No Iterations 1
GAMG: Solving for Uy, Initial residual = 0.332272, Final residual = 0.0140346, No Iterations 1
GAMG: Solving for Uz, Initial residual = 0.559406, Final residual = 0.0285128, No Iterations 1
GAMG: Solving for e, Initial residual = 0.574963, Final residual = 0.0496275, No Iterations 1
#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::hePsiThermo<Foam:siThermo, Foam:ureMixture<Foam::sutherlandTransport<Foam:: species::thermo<Foam::hConstThermo<Foam:erfectGa s<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::calculate() at ??:?
#4 Foam::hePsiThermo<Foam:siThermo, Foam:ureMixture<Foam::sutherlandTransport<Foam:: species::thermo<Foam::hConstThermo<Foam:erfectGa s<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::correct() at ??:?
#5 ? in "/opt/openfoam6/platforms/linux64GccDPInt32Opt/bin/rhoSimpleFoam"
#6 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#7 ? in "/opt/openfoam6/platforms/linux64GccDPInt32Opt/bin/rhoSimpleFoam"
Floating point exception (core dumped)


/*---------------------------------------------------------------------------------*/
Model: rhoSimpleFoam
/*---------------------------------------------------------------------------------*/
Boundary Condition:
1. Pressure
frontandback
{
type symmetry;


}
inlet
{
type zeroGradient;
//type mixed;
//refValue uniform 100000;
//refGradient uniform 0;
//valueFraction uniform 0.3;
}
outlet
{
type fixedValue;
value uniform 100000;
}
wall
{
type zeroGradient;
}


2. Velocity
frontandback
{
type symmetry;


}
inlet
{
type fixedValue;
value uniform (5 0 0);
rhoInlet 1.5;
}
outlet
{
type zeroGradient;
}
wall
{
type noSlip;
}


3. Temperature
frontandback
{
type symmetry;


}
inlet
{
type fixedValue;
value uniform 300;
}
outlet
{
type zeroGradient;
}
wall
{
type zeroGradient;
}


4. k
frontandback
{
type symmetry;


}
inlet
{
type turbulentMixingLengthDissipationRateInlet;
mixingLength 0.14;
value uniform 0.1;
}
outlet
{
type inletOutlet;
inletValue uniform 1;
value uniform 1;
}
wall
{
type kqRWallFunction;
value uniform 1;
}


5. epsilon
frontandback
{
type symmetry;


}
inlet
{
type turbulentMixingLengthDissipationRateInlet;
mixingLength 0.14;
value uniform 0.006;
}
outlet
{
type inletOutlet;
inletValue uniform 0.006;
value uniform 0.006;
}
wall
{
type epsilonWallFunction;
Cmu 0.09;
kappa 0.41;
E 9.8;
value uniform 0.006;
}


6. nut


frontandback
{
type symmetry;


}
inlet
{
type calculated;
value uniform 0;
}
outlet
{
type calculated;
value uniform 0;
}
wall
{
type nutkWallFunction;
Cmu 0.09;
kappa 0.41;
E 9.8;
value uniform 0;
}


7. alphat
frontandback
{
type symmetry;


}
inlet
{
type calculated;
value uniform 0;
}
outlet
{
type calculated;
value uniform 0;
}
wall
{
type compressible::alphatWallFunction;
Prt 0.85;
value uniform 0;
}


/*------------------------------------------------------------------------------*/
fvSolution
ddtSchemes
{
default steadyState;
}


gradSchemes
{
default Gauss linear;
}


divSchemes
{
default none;


div(phi,U) bounded Gauss upwind;
div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear;
div(phi,e) bounded Gauss upwind;
div(phi,epsilon) bounded Gauss upwind;
div(phi,k) bounded Gauss upwind;


div(phid,p) Gauss upwind;
div(phi,Ekp) bounded Gauss upwind;
div((phi|interpolate(rho)),p) Gauss upwind;
}


laplacianSchemes
{
default Gauss linear corrected;
}


interpolationSchemes
{
default linear;
}


snGradSchemes
{
default corrected;
}


/*----------------------------------------------------------------------*/
fvSolution
solvers
{
p
{
solver GAMG;
tolerance 1e-08;
relTol 0.1;
smoother GaussSeidel;
nCellsInCoarsestLevel 20;
}


"(U|e|k|epsilon)"
{
solver GAMG;
tolerance 1e-08;
relTol 0.1;
smoother GaussSeidel;
nCellsInCoarsestLevel 20;
}
}


SIMPLE
{
nNonOrthogonalCorrectors 0;
rhoMin 0.1;
rhoMax 1.5;
//pMaxFactor 2;
//pMinFactor 0.1;
transonic yes;
consistent yes;


residualControl
{
p 1e-3;
U 1e-4;
e 1e-3;


// possibly check turbulence fields
"(k|epsilon|omega)" 1e-3;
}
}


relaxationFactors
{
fields
{
p 0.3;
rho 0.7;
}
equations
{
p 0.7;
U 0.7;
e 0.7;
k 0.7;
epsilon 0.7;
}
}