CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM (https://www.cfd-online.com/Forums/openfoam/)
-   -   buoyantSimpleFoam (https://www.cfd-online.com/Forums/openfoam/118777-buoyantsimplefoam.html)

samiam1000 June 4, 2013 05:45

buoyantSimpleFoam
 
Dear All,

I am trying a 2D simulation, using the buoyantSimpleFoam solver.

I have a problem with the very first iteration, since I get this error:
Code:

lab@lab-laptop:~/Documenti/Ethics/FRISBEE/CFD/testVerticalTN_PCM/NOPCM/steady$ buoyantSimpleFoam
/*---------------------------------------------------------------------------*\
| =========                |                                                |
| \\      /  F ield        | OpenFOAM: The Open Source CFD Toolbox          |
|  \\    /  O peration    | Version:  2.1.0                                |
|  \\  /    A nd          | Web:      www.OpenFOAM.org                      |
|    \\/    M anipulation  |                                                |
\*---------------------------------------------------------------------------*/
Build  : 2.1.0-0bc225064152
Exec  : buoyantSimpleFoam
Date  : Jun 04 2013
Time  : 11:34:52
Host  : "lab-laptop"
PID    : 4155
Case  : /home/lab/Documenti/Ethics/FRISBEE/CFD/testVerticalTN_PCM/NOPCM/steady
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
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 laminar
Calculating field g.h

Reading field p_rgh


SIMPLE: convergence criteria
    field p_rgh    tolerance 0.001
    field U    tolerance 0.0001
    field T    tolerance 0.001


Starting time loop

Time = 1

DILUPBiCG:  Solving for Uy, Initial residual = 1, Final residual = 0.00427575, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 1, Final residual = 0.0103943, No Iterations 1
DILUPBiCG:  Solving for h, Initial residual = 1, Final residual = 0.0995407, No Iterations 7
DICPCG:  Solving for p_rgh, Initial residual = 0, Final residual = 0, No Iterations 0
#0  Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1  Foam::sigFpe::sigHandler(int) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2  in "/lib/libc.so.6"
#3  Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#4  void Foam::divide<Foam::fvsPatchField, Foam::surfaceMesh>(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&) in "/opt/openfoam210/platforms/linux64GccDPOpt/bin/buoyantSimpleFoam"
#5  Foam::tmp<Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> > Foam::operator/<Foam::fvsPatchField, Foam::surfaceMesh>(Foam::tmp<Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> > const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&) in "/opt/openfoam210/platforms/linux64GccDPOpt/bin/buoyantSimpleFoam"
#6 
 in "/opt/openfoam210/platforms/linux64GccDPOpt/bin/buoyantSimpleFoam"
#7  __libc_start_main in "/lib/libc.so.6"
#8 
 in "/opt/openfoam210/platforms/linux64GccDPOpt/bin/buoyantSimpleFoam"
Floating point exception
lab@lab-laptop:~/Documenti/Ethics/FRISBEE/CFD/testVerticalTN_PCM/NOPCM/steady$

Does anyone have an idea about how to solve this?

Thanks a lot,
Samuele


All times are GMT -4. The time now is 22:00.