|
[Sponsors] |
February 9, 2009, 16:11 |
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: 17 |
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 |
|
February 10, 2009, 03:49 |
Hi John,
With a floating po
|
#2 |
Senior Member
Mark Olesen
Join Date: Mar 2009
Location: https://olesenm.github.io/
Posts: 1,714
Rep Power: 40 |
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. |
|
February 10, 2009, 05:55 |
Also check that (absolute) pre
|
#3 |
Member
Ville Tossavainen
Join Date: Mar 2009
Posts: 60
Rep Power: 17 |
Also check that (absolute) pressure isn't zero anywhere in the field. Same thing with temperature (in Kelvins).
|
|
February 10, 2009, 07:25 |
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: 17 |
Thanks. I'll give it a try this morning, and let you know how it turns up.
|
|
February 18, 2009, 05:14 |
Hi,
maybe anyone knows this e
|
#5 |
New Member
Steven Darkes
Join Date: Mar 2009
Posts: 11
Rep Power: 17 |
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 |
|
February 18, 2009, 10:22 |
It means, that you've got nega
|
#6 |
Senior Member
Matvey Kraposhin
Join Date: Mar 2009
Location: Moscow, Russian Federation
Posts: 355
Rep Power: 21 |
It means, that you've got negative temperature
__________________
MDPI Fluids (Q2) special issue for OSS software: https://www.mdpi.com/journal/fluids/..._modelling_OSS GitHub: https://github.com/unicfdlab Linkedin: https://linkedin.com/in/matvey-kraposhin-413869163 RG: https://www.researchgate.net/profile/Matvey_Kraposhin |
|
February 18, 2009, 11:02 |
That's it! Thank's a lot!!
|
#7 |
New Member
Steven Darkes
Join Date: Mar 2009
Posts: 11
Rep Power: 17 |
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? |
|
February 19, 2009, 04:16 |
I think, you must look for in
|
#8 |
Senior Member
Matvey Kraposhin
Join Date: Mar 2009
Location: Moscow, Russian Federation
Posts: 355
Rep Power: 21 |
I think, you must look for in special literature
__________________
MDPI Fluids (Q2) special issue for OSS software: https://www.mdpi.com/journal/fluids/..._modelling_OSS GitHub: https://github.com/unicfdlab Linkedin: https://linkedin.com/in/matvey-kraposhin-413869163 RG: https://www.researchgate.net/profile/Matvey_Kraposhin |
|
February 19, 2009, 06:21 |
Now I got positive temperature
|
#9 |
New Member
Steven Darkes
Join Date: Mar 2009
Posts: 11
Rep Power: 17 |
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) |
|
February 19, 2009, 08:06 |
The error in converting h-valu
|
#10 |
Senior Member
Matvey Kraposhin
Join Date: Mar 2009
Location: Moscow, Russian Federation
Posts: 355
Rep Power: 21 |
The error in converting h-values to T. Can you post description of your case?
__________________
MDPI Fluids (Q2) special issue for OSS software: https://www.mdpi.com/journal/fluids/..._modelling_OSS GitHub: https://github.com/unicfdlab Linkedin: https://linkedin.com/in/matvey-kraposhin-413869163 RG: https://www.researchgate.net/profile/Matvey_Kraposhin |
|
February 19, 2009, 10:01 |
I sent it to you by mail!
|
#11 |
New Member
Steven Darkes
Join Date: Mar 2009
Posts: 11
Rep Power: 17 |
I sent it to you by mail!
|
|
March 3, 2009, 12:14 |
I don't know why, but I cannot
|
#12 |
New Member
Steven Darkes
Join Date: Mar 2009
Posts: 11
Rep Power: 17 |
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? |
|
March 3, 2009, 12:20 |
http://www.cfd-online.com/Open
|
#13 |
New Member
Steven Darkes
Join Date: Mar 2009
Posts: 11
Rep Power: 17 |
|
|
March 3, 2009, 12:25 |
http://www.cfd-online.com/Open
|
#14 |
New Member
Steven Darkes
Join Date: Mar 2009
Posts: 11
Rep Power: 17 |
|
|
March 3, 2009, 13:56 |
http://www.cfd-online.com/Open
|
#15 |
New Member
Steven Darkes
Join Date: Mar 2009
Posts: 11
Rep Power: 17 |
|
|
March 3, 2009, 14:19 |
hi, Steven!
1) i can"t see
|
#16 |
Senior Member
Matvey Kraposhin
Join Date: Mar 2009
Location: Moscow, Russian Federation
Posts: 355
Rep Power: 21 |
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.
__________________
MDPI Fluids (Q2) special issue for OSS software: https://www.mdpi.com/journal/fluids/..._modelling_OSS GitHub: https://github.com/unicfdlab Linkedin: https://linkedin.com/in/matvey-kraposhin-413869163 RG: https://www.researchgate.net/profile/Matvey_Kraposhin |
|
March 4, 2009, 12:04 |
Hi Matvej,
thank's a lot fo
|
#17 |
New Member
Steven Darkes
Join Date: Mar 2009
Posts: 11
Rep Power: 17 |
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! |
|
March 27, 2009, 10:44 |
|
#18 |
Senior Member
Ehsan
Join Date: Mar 2009
Posts: 112
Rep Power: 17 |
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 |
|
March 29, 2009, 12:05 |
|
#19 |
Senior Member
Matvey Kraposhin
Join Date: Mar 2009
Location: Moscow, Russian Federation
Posts: 355
Rep Power: 21 |
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? |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
BC for rhoSimpleFoam | uli_u | OpenFOAM Running, Solving & CFD | 77 | April 4, 2013 12:28 |
RhoSimpleFoam | maritozzo | OpenFOAM Running, Solving & CFD | 1 | April 30, 2010 19: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 13:26 |