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

RhoSimpleFoam

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

Reply
 
LinkBack Thread Tools Display Modes
Old   February 9, 2009, 16:11
Default Hi. I'm just wondering if anyo
  #1
New Member
 
John Paul Handrigan
Join Date: Mar 2009
Location: Ottawa, Ontario, Canada
Posts: 16
Rep Power: 8
jphandrigan is on a distinguished road
Hi. I'm just wondering if anyone can tell me what this error means. Thanks in advance.

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:stream&) 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
jphandrigan is offline   Reply With Quote

Old   February 10, 2009, 03:49
Default Hi John, With a floating po
  #2
Senior Member
 
Mark Olesen
Join Date: Mar 2009
Location: http://olesenm.github.io/
Posts: 777
Rep Power: 18
olesen will become famous soon enough
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.
olesen is offline   Reply With Quote

Old   February 10, 2009, 05:55
Default Also check that (absolute) pre
  #3
Member
 
Ville Tossavainen
Join Date: Mar 2009
Location: Helsinki, Finland
Posts: 60
Rep Power: 8
villet is on a distinguished road
Also check that (absolute) pressure isn't zero anywhere in the field. Same thing with temperature (in Kelvins).
villet is offline   Reply With Quote

Old   February 10, 2009, 07:25
Default Thanks. I'll give it a try thi
  #4
New Member
 
John Paul Handrigan
Join Date: Mar 2009
Location: Ottawa, Ontario, Canada
Posts: 16
Rep Power: 8
jphandrigan is on a distinguished road
Thanks. I'll give it a try this morning, and let you know how it turns up.
jphandrigan is offline   Reply With Quote

Old   February 18, 2009, 05:14
Default Hi, maybe anyone knows this e
  #5
New Member
 
Steven Darkes
Join Date: Mar 2009
Posts: 11
Rep Power: 8
stevendarkes is on a distinguished road
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:stream&) 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
stevendarkes is offline   Reply With Quote

Old   February 18, 2009, 10:22
Default It means, that you've got nega
  #6
Senior Member
 
mkraposhin's Avatar
 
Matvej Kraposhin
Join Date: Mar 2009
Location: Moscow, Russian Federation
Posts: 197
Rep Power: 9
mkraposhin is on a distinguished road
It means, that you've got negative temperature
mkraposhin is offline   Reply With Quote

Old   February 18, 2009, 11:02
Default That's it! Thank's a lot!!
  #7
New Member
 
Steven Darkes
Join Date: Mar 2009
Posts: 11
Rep Power: 8
stevendarkes is on a distinguished road
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?
stevendarkes is offline   Reply With Quote

Old   February 19, 2009, 04:16
Default I think, you must look for in
  #8
Senior Member
 
mkraposhin's Avatar
 
Matvej Kraposhin
Join Date: Mar 2009
Location: Moscow, Russian Federation
Posts: 197
Rep Power: 9
mkraposhin is on a distinguished road
I think, you must look for in special literature
mkraposhin is offline   Reply With Quote

Old   February 19, 2009, 06:21
Default Now I got positive temperature
  #9
New Member
 
Steven Darkes
Join Date: Mar 2009
Posts: 11
Rep Power: 8
stevendarkes is on a distinguished road
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:stream&) 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)
stevendarkes is offline   Reply With Quote

Old   February 19, 2009, 08:06
Default The error in converting h-valu
  #10
Senior Member
 
mkraposhin's Avatar
 
Matvej Kraposhin
Join Date: Mar 2009
Location: Moscow, Russian Federation
Posts: 197
Rep Power: 9
mkraposhin is on a distinguished road
The error in converting h-values to T. Can you post description of your case?
mkraposhin is offline   Reply With Quote

Old   February 19, 2009, 10:01
Default I sent it to you by mail!
  #11
New Member
 
Steven Darkes
Join Date: Mar 2009
Posts: 11
Rep Power: 8
stevendarkes is on a distinguished road
I sent it to you by mail!
stevendarkes is offline   Reply With Quote

Old   March 3, 2009, 12:14
Default I don't know why, but I cannot
  #12
New Member
 
Steven Darkes
Join Date: Mar 2009
Posts: 11
Rep Power: 8
stevendarkes is on a distinguished road
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?


stevendarkes is offline   Reply With Quote

Old   March 3, 2009, 12:20
Default http://www.cfd-online.com/Open
  #13
New Member
 
Steven Darkes
Join Date: Mar 2009
Posts: 11
Rep Power: 8
stevendarkes is on a distinguished road

stevendarkes is offline   Reply With Quote

Old   March 3, 2009, 12:25
Default http://www.cfd-online.com/Open
  #14
New Member
 
Steven Darkes
Join Date: Mar 2009
Posts: 11
Rep Power: 8
stevendarkes is on a distinguished road

stevendarkes is offline   Reply With Quote

Old   March 3, 2009, 13:56
Default http://www.cfd-online.com/Open
  #15
New Member
 
Steven Darkes
Join Date: Mar 2009
Posts: 11
Rep Power: 8
stevendarkes is on a distinguished road

stevendarkes is offline   Reply With Quote

Old   March 3, 2009, 14:19
Default hi, Steven! 1) i can"t see
  #16
Senior Member
 
mkraposhin's Avatar
 
Matvej Kraposhin
Join Date: Mar 2009
Location: Moscow, Russian Federation
Posts: 197
Rep Power: 9
mkraposhin is on a distinguished road
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.
mkraposhin is offline   Reply With Quote

Old   March 4, 2009, 12:04
Default Hi Matvej, thank's a lot fo
  #17
New Member
 
Steven Darkes
Join Date: Mar 2009
Posts: 11
Rep Power: 8
stevendarkes is on a distinguished road
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!
stevendarkes is offline   Reply With Quote

Old   March 27, 2009, 10:44
Default
  #18
Member
 
ehsan
Join Date: Mar 2009
Posts: 98
Rep Power: 8
ehsan is on a distinguished road
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
ehsan is offline   Reply With Quote

Old   March 29, 2009, 11:05
Default
  #19
Senior Member
 
mkraposhin's Avatar
 
Matvej Kraposhin
Join Date: Mar 2009
Location: Moscow, Russian Federation
Posts: 197
Rep Power: 9
mkraposhin is on a distinguished road
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?
mkraposhin 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
BC for rhoSimpleFoam uli_u OpenFOAM Running, Solving & CFD 77 April 4, 2013 11:28
RhoSimpleFoam maritozzo OpenFOAM Running, Solving & CFD 1 April 30, 2010 18:18
RhoPorousSimpleFoam RhoSimpleFoam spv24 OpenFOAM Running, Solving & CFD 4 November 12, 2008 10:19
Problem with rhoSimpleFoam matteo_gautero OpenFOAM Running, Solving & CFD 0 February 28, 2008 07:51
RhoSimpleFoam luca OpenFOAM Running, Solving & CFD 0 August 22, 2005 12:26


All times are GMT -4. The time now is 03:24.