CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Post-Processing

wallHeatFlux error

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   November 23, 2014, 19:38
Default wallHeatFlux error
  #1
New Member
 
Raghav
Join Date: Oct 2014
Posts: 12
Rep Power: 7
raghav.venky is on a distinguished road
I'm trying to simulate a forced convection problem using the buoyantBoussinesqSimpleFoam solver.

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

Code:
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
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
Attached Files
File Type: gz constant.tar.gz (938 Bytes, 4 views)
raghav.venky is offline   Reply With Quote

Old   November 24, 2014, 06:27
Default
  #2
Senior Member
 
Joachim Herb
Join Date: Sep 2010
Posts: 594
Rep Power: 17
jherb is on a distinguished road
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.
jherb is offline   Reply With Quote

Old   December 22, 2014, 07:49
Default
  #3
New Member
 
Muhammad Omer Mughal
Join Date: Jul 2010
Location: Singapore
Posts: 22
Rep Power: 12
Muhammad Omer Mughal is on a distinguished road
Dear Raghav

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

Thanks
Muhammad Omer Mughal is offline   Reply With Quote

Old   December 31, 2018, 03:17
Default
  #4
Member
 
Bidesh Sengupta
Join Date: Sep 2018
Posts: 68
Rep Power: 3
BSengupta is on a distinguished road
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;
}
}
BSengupta is offline   Reply With Quote

Reply

Tags
rasmodel, turbulence, wallheatflux

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
[OpenFOAM.org] compile error in dynamicMesh and thermophysicalModels libraries NickG OpenFOAM Installation 3 December 30, 2019 00:21
[blockMesh] blockMesh with double grading. spwater OpenFOAM Meshing & Mesh Conversion 92 January 12, 2019 09:00
[swak4Foam] GroovyBC the dynamic cousin of funkySetFields that lives on the suburb of the mesh gschaider OpenFOAM Community Contributions 300 October 29, 2014 18:00
OpenFOAM without MPI kokizzu OpenFOAM Installation 4 May 26, 2014 09:17
DecomposePar links against liblamso0 with OpenMPI jens_klostermann OpenFOAM Bugs 11 June 28, 2007 17:51


All times are GMT -4. The time now is 16:38.