wolfindark |
March 11, 2013 00:42 |
buoyantBoussinesqPimpleFoam & turbulentHeatFluxTemperature
Dear Foamers
I try to simulate natural convection and temporal variation of temperature distribution inside a water pool.
with:
OpenFOAM version: 2.1.1.
Solver : buoyantBoussinesqPimpleFoam
RAS : kEps turbulent
I use turbulentHeatFluxTemperature for a constant heat input from a surface with following block in 0/T file:
Code:
PIPEOUT
{
type turbulentHeatFluxTemperature;
heatSource flux; // power [W]; flux [W/m2]
q uniform 10; // heat power or flux
alphaEff kappat; // alphaEff field name;
// alphaEff in [kg/m/s]
Cp Cp; // Cp field name; Cp in [J/kg/K]
value uniform 300; // initial temperature value
}
Code runs but it gives following error in the second time step. I would appreciate if you could give me an idea where the error is originated?
Error:
Code:
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.1.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 2.1.1-221db2718bbb
Exec : buoyantBoussinesqPimpleFoam
Date : Mar 11 2013
Time : 14:37:12
Host : "ubun"
PID : 8648
Case : /home/neo/OpenFOAM/neo-2.1.1/run/exp_tr
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time
Create mesh for time = 0
Reading g
Reading thermophysical properties
Reading field T
Reading field p_rgh
Reading field U
Reading/calculating face flux field phi
Selecting incompressible transport model Newtonian
Creating turbulence model
Selecting RAS turbulence model kEpsilon
kEpsilonCoeffs
{
Cmu 0.09;
C1 1.44;
C2 1.92;
sigmaEps 1.3;
}
Reading field kappat
Calculating field g.h
Courant Number mean: 0 max: 0
PIMPLE: Operating solver in PISO mode
Starting time loop
Time = 1
Courant Number mean: 0 max: 0
DILUPBiCG: Solving for T, Initial residual = 1, Final residual = 9.88958e-07, No Iterations 208
DICPCG: Solving for p_rgh, Initial residual = 1, Final residual = 0.0093319, No Iterations 160
time step continuity errors : sum local = 1.97326e-09, global = 1.3291e-21, cumulative = 1.3291e-21
DICPCG: Solving for p_rgh, Initial residual = 0.00686217, Final residual = 9.33739e-09, No Iterations 212
time step continuity errors : sum local = 4.3591e-14, global = 3.96046e-14, cumulative = 3.96046e-14
DILUPBiCG: Solving for epsilon, Initial residual = 0.0714385, Final residual = 8.4943e-07, No Iterations 50
DILUPBiCG: Solving for k, Initial residual = 1, Final residual = 9.70453e-07, No Iterations 100
ExecutionTime = 5.43 s ClockTime = 5 s
Time = 2
Courant Number mean: 5.30969e-08 max: 3.99768e-07
#0 Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#4 Foam::operator/(Foam::UList<double> const&, Foam::tmp<Foam::Field<double> > const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#5 Foam::incompressible::turbulentHeatFluxTemperatureFvPatchScalarField::updateCoeffs() in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
#6 at gaussLaplacianSchemes.C:0
#7 Foam::fv::gaussLaplacianScheme<double, double>::fvmLaplacianUncorrected(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#8 Foam::fv::gaussLaplacianScheme<double, double>::fvmLaplacian(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#9 Foam::fv::laplacianScheme<double, double>::fvmLaplacian(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#10
in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/buoyantBoussinesqPimpleFoam"
#11 Foam::tmp<Foam::fvMatrix<double> > Foam::fvm::laplacian<double, double>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/buoyantBoussinesqPimpleFoam"
#12
in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/buoyantBoussinesqPimpleFoam"
#13 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#14
in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/buoyantBoussinesqPimpleFoam"
Floating point exception (core dumped)
|