CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (http://www.cfd-online.com/Forums/openfoam-solving/)
-   -   RhoSimpleFoam (http://www.cfd-online.com/Forums/openfoam-solving/57860-rhosimplefoam.html)

jphandrigan February 9, 2009 15:11

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

olesen February 10, 2009 02:49

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.

villet February 10, 2009 04:55

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

jphandrigan February 10, 2009 06:25

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

stevendarkes February 18, 2009 04:14

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

mkraposhin February 18, 2009 09:22

It means, that you've got nega
 
It means, that you've got negative temperature

stevendarkes February 18, 2009 10:02

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?

mkraposhin February 19, 2009 03:16

I think, you must look for in
 
I think, you must look for in special literature

stevendarkes February 19, 2009 05:21

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)

mkraposhin February 19, 2009 07:06

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

stevendarkes February 19, 2009 09:01

I sent it to you by mail!
 
I sent it to you by mail!

stevendarkes March 3, 2009 11:14

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

stevendarkes March 3, 2009 11:20

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

stevendarkes March 3, 2009 11:25

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

stevendarkes March 3, 2009 12:56

http://www.cfd-online.com/Open
 
http://www.cfd-online.com/OpenFOAM_D...es/1/11491.jpg

mkraposhin March 3, 2009 13:19

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.

stevendarkes March 4, 2009 11:04

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!

ehsan March 27, 2009 10:44

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

mkraposhin March 29, 2009 11:05

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?


All times are GMT -4. The time now is 18:59.