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
|
http://www.cfd-online.com/Open
|
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? |
All times are GMT -4. The time now is 08:48. |