CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)

- OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)

- - RhoSimpleFoam (https://www.cfd-online.com/Forums/openfoam-solving/57860-rhosimplefoam.html)

Hi. I'm just wondering if anyo

Hi. I'm just wondering if anyone can tell me what this error means. Thanks in advance. http://www.cfd-online.com/OpenFOAM_D...part/happy.gif

Here's the print-out I get when I run the solver:

[actmodel@r4-bcc-a010590 rhoSimpleFoam]$ rhoSimpleFoam
/*---------------------------------------------------------------------------*\
| ========= | |
| \ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \ / O peration | Version: 1.5 |
| \ / A nd | Web: http://www.OpenFOAM.org |
| \/ M anipulation | |
\*---------------------------------------------------------------------------*/
Exec : rhoSimpleFoam
Date : Feb 09 2009
Time : 15:09:20
PID : 20392
Case : /mnt/hypatia/linux_home/rhel4_home/actmodel/OpenFOAM/junk-1.5/run/rhoSimpleFoam
nProcs : 1

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

Create mesh for time = 0

Reading thermophysical properties

Selecting thermodynamics package hThermo<puremixture<sutherlandtransport<speciether mo<hconstthermo<perfectgas>>>> >
Reading field U

Reading/calculating face flux field phi

Creating turbulence model

Selecting RAS turbulence model kEpsilon
#0 Foam::error::printStack(Foam:http://www.cfd-online.com/OpenFOAM_D...part/proud.gifstream&) in "/mnt/hypatia/linux_home/rhel4_home/actmodel/OpenFOAM/OpenFOAM-1.5/lib/linux64Gc cDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/mnt/hypatia/linux_home/rhel4_home/actmodel/OpenFOAM/OpenFOAM-1.5/lib/linux64Gc cDPOpt/libOpenFOAM.so"
#2 __restore_rt in "/lib64/tls/libc.so.6"
#3 Foam::compressible::RASModels::kEpsilon::kEpsilon( Foam::GeometricField<double,> const&, Foam::GeometricField<foam::vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double,> const&, Foam::basicThermo&) in "/mnt/hypatia/linux_home/rhel4_home/actmodel/OpenFOAM/OpenFOAM-1.5/lib/linux64Gc cDPOpt/libcompressibleRASModels.so"
#4 Foam::compressible::RASModel::adddictionaryConstru ctorToTable<foam::compressible ::rasmodels::kepsilon>::New(Foam::GeometricField<d ouble,> const&, Foam::GeometricField<foam::vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double,> const&, Foam::basicThermo&) in "/mnt/hypatia/linux_home/rhel4_home/actmodel/OpenFOAM/OpenFOAM-1.5/lib/linux64Gc cDPOpt/libcompressibleRASModels.so"
#5 Foam::compressible::RASModel::New(Foam::GeometricF ield<double,> const&, Foam::GeometricField<foam::vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double,> const&, Foam::basicThermo&) in "/mnt/hypatia/linux_home/rhel4_home/actmodel/OpenFOAM/OpenFOAM-1.5/lib/linux64Gc cDPOpt/libcompressibleRASModels.so"
#6 main in "/mnt/hypatia/linux_home/rhel4_home/actmodel/OpenFOAM/OpenFOAM-1.5/applications/ bin/linux64GccDPOpt/rhoSimpleFoam"
#7 __libc_start_main in "/lib64/tls/libc.so.6"
#8 Foam::regIOobject::readIfModified() in "/mnt/hypatia/linux_home/rhel4_home/actmodel/OpenFOAM/OpenFOAM-1.5/applications/ bin/linux64GccDPOpt/rhoSimpleFoam"
Floating exception

Hi John, With a floating po

Hi John,

With a floating point exception at the very beginning, check that your fields are initialized to something sensible. Eg, don't use zero for the tke.

Also check that (absolute) pre

Also check that (absolute) pressure isn't zero anywhere in the field. Same thing with temperature (in Kelvins).

Thanks. I'll give it a try thi

Thanks. I'll give it a try this morning, and let you know how it turns up. http://www.cfd-online.com/OpenFOAM_D...part/happy.gif

Hi, maybe anyone knows this e

Hi,
maybe anyone knows this error as well. Thanks in advance!

Here's the print-out I get when I run the rhoSimpleFoam solver:

...

Time = 6

smoothSolver: Solving for Ux, Initial residual = 0.300191, Final residual = 0.0129647, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 0.0539535, Final residual = 0.00221057, No Iterations 4
smoothSolver: Solving for Uz, Initial residual = 0.996496, Final residual = 0.0279006, No Iterations 2
DILUPBiCG: Solving for h, Initial residual = 0.0954665, Final residual = 0.00199884, No Iterations 1

Maximum number of iterations exceeded#0 Foam::error::printStack(Foam:http://www.cfd-online.com/OpenFOAM_D...part/proud.gifstream&) in "/home/sda2sh/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libOpenFOAM.so"
#1 Foam::error::abort() in "/home/sda2sh/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libOpenFOAM.so"
#2 Foam::hThermo<foam::puremixture<foam::consttranspo rt<foam::speciethermo<foam::hc onstthermo<foam::perfectgas> > > > >::calculate() in "/home/sda2sh/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libbasicThermophysicalMode ls.so"
#3 Foam::hThermo<foam::puremixture<foam::consttranspo rt<foam::speciethermo<foam::hc onstthermo<foam::perfectgas> > > > >::correct() in "/home/sda2sh/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libbasicThermophysicalMode ls.so"
#4 main in "/home/sda2sh/OpenFOAM/OpenFOAM-1.5/applications/bin/linuxGccDPOpt/rhoSimpleFoam "
#5 __libc_start_main in "/lib/libc.so.6"
#6 Foam::regIOobject::readIfModified() in "/home/sda2sh/OpenFOAM/OpenFOAM-1.5/applications/bin/linuxGccDPOpt/rhoSimpleFoam "

From function specieThermo<thermo>::T(scalar f, scalar T0, scalar (specieThermo<thermo>::*F)(const scalar) const, scalar (specieThermo<thermo>::*dFdT)(const scalar) const) const
in file /home/dm2/henry/OpenFOAM/OpenFOAM-1.5/src/thermophysicalModels/specie/lnInclude/ specieThermoI.H at line 78.

FOAM aborting

It means, that you've got nega

It means, that you've got negative temperature

That's it! Thank's a lot!!

That's it! Thank's a lot!!

Maybe I should change my viscosity. E.g. define it as a function of the Temperature.
What I didn't understand up to now in the thermoPhysicalProperties is the heat fusion Hf. How can I get this value?

I think, you must look for in

I think, you must look for in special literature

Now I got positive temperature

Now I got positive temperatures, but anyway the solver aborts after a couple of iterations, probably because of division by zero!

Maybe you know what the reason could be?

Here the error message and below my thermoPhysicalProperties, because I think the problem is coming from there.

-----------------------------------------

Time = 7

smoothSolver: Solving for Ux, Initial residual = 0.999964, Final residual = 0.0151332, No Iterations 4
smoothSolver: Solving for Uy, Initial residual = 0.999954, Final residual = 0.03379, No Iterations 4
smoothSolver: Solving for Uz, Initial residual = 0.655577, Final residual = 0.0122258, No Iterations 2
DILUPBiCG: Solving for h, Initial residual = 0.598926, Final residual = 0.026901, No Iterations 1
#0 Foam::error::printStack(Foam:http://www.cfd-online.com/OpenFOAM_D...part/proud.gifstream&) in "/home/sda2sh/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/home/sda2sh/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libOpenFOAM.so"
#2 Uninterpreted: [0xffffe420]
#3 Foam::hThermo<foam::puremixture<foam::sutherlandtr ansport<foam::speciethermo<foa m::hconstthermo<foam::perfectgas> > > > >::calculate() in "/home/sda2sh/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libbasicThermophysicalMode ls.so"
#4 Foam::hThermo<foam::puremixture<foam::sutherlandtr ansport<foam::speciethermo<foa m::hconstthermo<foam::perfectgas> > > > >::correct() in "/home/sda2sh/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libbasicThermophysicalMode ls.so"
#5 main in "/home/sda2sh/OpenFOAM/OpenFOAM-1.5/applications/bin/linuxGccDPOpt/rhoSimpleFoam "
#6 __libc_start_main in "/lib/libc.so.6"
#7 Foam::regIOobject::readIfModified() in "/home/sda2sh/OpenFOAM/OpenFOAM-1.5/applications/bin/linuxGccDPOpt/rhoSimpleFoam "
Gleitkomma-Ausnahme

-----------------------------------------

thermoType hThermo<puremixture<sutherlandtransport<speciether mo<hconstthermo<perfectgas>>>> >;

mixture
argon 1 39.948 // specieThermo
520 0 // hConstTherm
1.96566e-06 144.4 // transportCoeffs/sutherlandTransport (Star-CD)

The error in converting h-valu

The error in converting h-values to T. Can you post description of your case?

I sent it to you by mail!

I sent it to you by mail!

I don't know why, but I cannot

I don't know why, but I cannot get a converging result for this geometry using the left edge as symmetry axis and consequently the front and back surface as wedge boundary!
Do you think this problem could be resolved by a full 360 degree geometry?

http://www.cfd-online.com/OpenFOAM_D...your_image.gif

http://www.cfd-online.com/Open

hi, Steven! 1) i can"t see

hi, Steven!

1) i can"t see your images

2) why you want to use 'wedge' patch? maybe you must try symmetryPlane for front and back surfaces.

Hi Matvej, thank's a lot fo

Hi Matvej,

thank's a lot for your manipulation of my files. they worked fine, but as soon as I put wedge boundaries it doesn't converge any more.
Even symmetryPlanes cause the same problem!

Hi Steven Darkes,

I am facing the same your problem hthermo before. So did you have a solution for this? If so please let me know. Thanks

Hi, Hung Nguyen

after some experience with rhoSimpleFoam (or rhoTurbFoam), i found, that all errors, that arises when converting h-values to T-values means an improper case setup (pressure boundaries, turbulence initial setup, physical BC, etc). Can you tell more about your case?