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

BuoyantsimpleFoam error

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

Reply
 
LinkBack Thread Tools Display Modes
Old   April 5, 2012, 03:51
Default BuoyantsimpleFoam error
  #1
New Member
 
Anone
Join Date: Feb 2012
Posts: 16
Rep Power: 4
CFD_user_2012 is on a distinguished road
Dear all

I am running a simulation of the air flow in a room using buoyant simple foam.

I run the program and get the following error message

Not sure whats wrong

Please help

CFD_user


Create time

Create mesh for time = 0


Reading g
Reading thermophysical properties

Selecting thermodynamics package hPsiThermo<pureMixture<constTransport<specieThermo <hConstThermo<perfectGas>>>>>
Reading field U

Reading/calculating face flux field phi

Creating turbulence model

Selecting RAS turbulence model kEpsilon
kEpsilonCoeffs
{
Cmu 0.09;
C1 1.44;
C2 1.92;
C3 -0.33;
sigmak 1;
sigmaEps 1.3;
Prt 1;
}

Calculating field g.h

Reading field p_rgh


SIMPLE: convergence criteria
field p_rgh tolerance 0.01
field U tolerance 0.001
field h tolerance 0.001
field "(k|epsilon|omega)" tolerance 0.001


Starting time loop

Time = 1

DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 0.00488841, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 0.000946431, No Iterations 1
DILUPBiCG: Solving for Uz, Initial residual = 1, Final residual = 0.00347572, No Iterations 1
DILUPBiCG: Solving for h, Initial residual = 1, Final residual = 0.00107994, No Iterations 1
DICPCG: Solving for p_rgh, Initial residual = 0.653293, Final residual = 0.00648429, No Iterations 82
time step continuity errors : sum local = 0.658898, global = -4.7659e-17, cumulative = -4.7659e-17
rho max/min : 1.70542 0.61239
DILUPBiCG: Solving for epsilon, Initial residual = 0.99778, Final residual = 0.0385943, No Iterations 1
DILUPBiCG: Solving for k, Initial residual = 1, Final residual = 0.022947, No Iterations 1
ExecutionTime = 11.81 s ClockTime = 12 s

Time = 2

DILUPBiCG: Solving for Ux, Initial residual = 0.132736, Final residual = 0.00281505, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 0.0335339, Final residual = 0.000215383, No Iterations 1
DILUPBiCG: Solving for Uz, Initial residual = 0.100024, Final residual = 0.00168842, No Iterations 1
DILUPBiCG: Solving for h, Initial residual = 1, Final residual = 0.00449013, No Iterations 1
DICPCG: Solving for p_rgh, Initial residual = 0.999999, Final residual = 0.859615, No Iterations 1001
time step continuity errors : sum local = 8.87527, global = 5.38735e-17, cumulative = 6.21446e-18
rho max/min : 362356 -1.21341e+06
DILUPBiCG: Solving for epsilon, Initial residual = 0.957201, Final residual = 0.00617101, No Iterations 2
#0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam210/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam210/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#2 Uninterpreted:
#3 log in "/lib/i386-linux-gnu/libm.so.6"
#4 Foam::compressible::RASModels::mutUWallFunctionFvP atchScalarField::calcYPlus(Foam::Field<double> const&) const in "/opt/openfoam210/platforms/linuxGccDPOpt/lib/libcompressibleRASModels.so"
#5 Foam::compressible::RASModels::mutUWallFunctionFvP atchScalarField::calcMut() const in "/opt/openfoam210/platforms/linuxGccDPOpt/lib/libcompressibleRASModels.so"
#6 Foam::compressible::RASModels::mutkWallFunctionFvP atchScalarField::updateCoeffs() in "/opt/openfoam210/platforms/linuxGccDPOpt/lib/libcompressibleRASModels.so"
#7 Foam::fvMatrix<double>::fvMatrix(Foam::GeometricFi eld<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::dimensionSet const&) in "/opt/openfoam210/platforms/linuxGccDPOpt/bin/buoyantSimpleFoam"
#8 Foam::tmp<Foam::fvMatrix<double> > Foam::fvm::Sp<double>(Foam:imensionedField<doubl e, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/openfoam210/platforms/linuxGccDPOpt/bin/buoyantSimpleFoam"
#9 Foam::compressible::RASModels::kEpsilon::correct() in "/opt/openfoam210/platforms/linuxGccDPOpt/lib/libcompressibleRASModels.so"
#10
in "/opt/openfoam210/platforms/linuxGccDPOpt/bin/buoyantSimpleFoam"
#11 __libc_start_main in "/lib/i386-linux-gnu/libc.so.6"
#12
in "/opt/openfoam210/platforms/linuxGccDPOpt/bin/buoyantSimpleFoam"
Floating point exception
openfoam@openfoam-VirtualBox:~/OpenFOAM/openfoam-2.1.0/run/tutorials/heatTransfer/buoyantSimpleFoam/hotRoom$
CFD_user_2012 is offline   Reply With Quote

Old   December 5, 2013, 17:15
Default
  #2
New Member
 
Join Date: Dec 2013
Posts: 1
Rep Power: 0
MDeza is on a distinguished road
I've got the same error. Did you ever figure out what was causing it?
MDeza is offline   Reply With Quote

Old   December 6, 2013, 06:27
Default
  #3
Senior Member
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Augsburg (Bavaria, Germany) Leoben (Styria, Austria)
Posts: 845
Rep Power: 17
Tobi will become famous soon enough
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
Hi,

do you have this error on your case or the tutorial case?

As you see the problem depend on the calculation of some equations:
Code:
DICPCG: Solving for p_rgh, Initial residual = 0.999999, Final residual = 0.859615, No Iterations 1001
time step continuity errors : sum local = 8.87527, global = 5.38735e-17, cumulative = 6.21446e-18
rho max/min : 362356 -1.21341e+06
DILUPBiCG: Solving for epsilon, Initial residual = 0.957201, Final residual = 0.00617101, No Iterations 2
In his case I think the pressure BC or volocity BC are not set correct. Therefor he got a blow up on p_rgh equation. Iterations of 1000 indicate a Problem with BC or Initial solution - especially in the first time steps. Additionally you see the residuals and the min/max values of rho.

I think a fluid never will get negativ density or 362356 kg/m - seems that he find a new - till now non found - fluid.

You should be able to solve your problem now.

Additionally:
Code:
Floating point exception
This error message give you a hint that you devide with Zero.
Tobi is offline   Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
GroovyBC the dynamic cousin of funkySetFields that lives on the suburb of the mesh gschaider OpenFOAM 294 July 12, 2013 05:39
Compile problem ivanyao OpenFOAM Running, Solving & CFD 1 October 12, 2012 09:31
c++ libraries and solver compiling vaina74 OpenFOAM Installation 13 February 3, 2012 17:43
Problem with compile the setParabolicInlet ivanyao OpenFOAM Running, Solving & CFD 6 September 5, 2008 20:50
user defined function cfduser CFX 0 April 29, 2006 10:58


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