# Running dieselFoam error

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

 November 6, 2008, 05:29 I have include temperature fie #81 emilianyassenov Guest   Posts: n/a I have include temperature field in IcoFoam and I have run it after a while it is stopping. Diameter of pipe is 0.001m and length is 1m, velocity is 0.212 m/s. it gives me this message Time = 0.085 Courant Number mean: 5.92365e+55 max: 2.22812e+59 DILUPBiCG: Solving for Ux, Initial residual = 0.999827, Final residual = 28.78, No Iterations 1001 DILUPBiCG: Solving for Uy, Initial residual = 0.999891, Final residual = 13.7978, No Iterations 1001 DILUPBiCG: Solving for Uz, Initial residual = 0.993418, Final residual = 5.55113, No Iterations 1001 DICPCG: Solving for p, Initial residual = 1, Final residual = 0.728859, No Iterations 1001 time step continuity errors : sum local = 1.93343e+63, global = -4.09326e+57, cumulative = -4.09326e+57 DICPCG: Solving for p, Initial residual = 0.978532, Final residual = 18.7317, No Iterations 1001 time step continuity errors : sum local = 1.56696e+66, global = 1.27684e+62, cumulative = 1.2768e+62 #0 Foam::error::printStack(Foam:stream&) in "/home/rkahraman/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so" #1 Foam::sigFpe::sigFpeHandler(int) in "/home/rkahraman/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so" #2 ?? in "/lib64/libc.so.6" #3 Foam::PBiCG::solve(Foam::Field&, Foam::Field const&, unsigned char) const in "/home/rkahraman/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so" #4 Foam::fvMatrix::solve(Foam::Istream&) in "/home/rkahraman/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libfiniteVolume.so" #5 main in "/home/rkahraman/OpenFOAM/OpenFOAM-1.5/applications/bin/linux64GccDPOpt/my_icoFo am" #6 __libc_start_main in "/lib64/libc.so.6" #7 Foam::regIOobject::readIfModified() in "/home/rkahraman/OpenFOAM/OpenFOAM-1.5/applications/bin/linux64GccDPOpt/my_icoFo am" Gleitkomma-Ausnahme thanks for suggestion Emo

 November 6, 2008, 06:00 Hello Sebastian, it gives t #82 emilianyassenov Guest   Posts: n/a Hello Sebastian, it gives the same error like previous one with solver rhoPimpleFoam. it is stopping..What is wrong can you help me? thanks in advance Emo

 November 6, 2008, 08:38 Hello Sebastian again, I ha #83 emilianyassenov Guest   Posts: n/a Hello Sebastian again, I have run my case with solver rhoPimpleFoam. it gives me the following problem Starting time loop Courant Number mean: 0 max: -0 Time = 1 #0 Foam::error::printStack(Foam:stream&) in "/home/rkahraman/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so" #1 Foam::sigFpe::sigFpeHandler(int) in "/home/rkahraman/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so" #2 ?? in "/lib64/libc.so.6" #3 void Foam::fvc::surfaceIntegrate >(Foam::Field >&, Foam::GeometricField, Foam::fvsPatchField, Foam::surfaceMesh> const&) in "/home/rkahraman/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libfiniteVolume.so" #4 Foam::tmp, Foam::fvPatchField, Foam::volMesh> > Foam::fvc::surfaceIntegrate >(Foam::GeometricField, Foam::fvsPatchField, Foam::surfaceMesh> const&) in "/home/rkahraman/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libfiniteVolume.so" #5 Foam::tmp, Foam::fvPatchField, Foam::volMesh> > Foam::fvc::div >(Foam::GeometricField, Foam::fvsPatchField, Foam::surfaceMesh> const&) in "/home/rkahraman/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libfiniteVolume.so" #6 Foam::tmp, Foam::fvPatchField, Foam::volMesh> > Foam::fvc::div >(Foam::tmp, Foam::fvsPatchField, Foam::surfaceMesh> > const&) in "/home/rkahraman/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libfiniteVolume.so" #7 Foam::fv::gaussLaplacianScheme, double>::fvmLaplacian(Foam::GeometricField const&, Foam::GeometricField, Foam::fvPatchField, Foam::volMesh>&) in "/home/rkahraman/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libfiniteVolume.so" #8 Foam::fv::laplacianScheme, double>::fvmLaplacian(Foam::GeometricField const&, Foam::GeometricField, Foam::fvPatchField, Foam::volMesh>&) in "/home/rkahraman/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libfiniteVolume.so" #9 Foam::tmp > > Foam::fvm::laplacian, double>(Foam::GeometricField const&, Foam::GeometricField, Foam::fvPatchField, Foam::volMesh>&, Foam::word const&) in "/home/rkahraman/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libcompressibleRASMod els.so" #10 Foam::tmp > > Foam::fvm::laplacian, double>(Foam::GeometricField const&, Foam::GeometricField, Foam::fvPatchField, Foam::volMesh>&) in "/home/rkahraman/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libcompressibleRASMod els.so" #11 Foam::compressible::RASModels::kEpsilon::divDevRho Reff(Foam::GeometricField, Foam::fvPatchField, Foam::volMesh>&) const in "/home/rkahraman/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libcompressibleRASMod els.so" #12 main in "/home/rkahraman/OpenFOAM/rkahraman-1.5/applications/bin/linux64GccDPOpt/rhoPimp leFoam" #13 __libc_start_main in "/lib64/libc.so.6" #14 Foam::regIOobject::readIfModified() in "/home/rkahraman/OpenFOAM/rkahraman-1.5/applications/bin/linux64GccDPOpt/rhoPimp leFoam" Gleitkomma-Ausnahme

 November 6, 2008, 10:26 can someone help me? best reg #84 emilianyassenov Guest   Posts: n/a can someone help me? best regards Emo

 November 6, 2008, 10:40 Hi Emilian about your 03:29 #85 Senior Member   Cedric DUPRAT Join Date: Mar 2009 Location: Nantes, France Posts: 195 Rep Power: 16 Hi Emilian about your 03:29 post, your mistake occure before, because, this time step show a time step continuity error of 1.56696e+66 which is ..... wrong. about your last post now, I don't know the solution but, keep in mind that if in the message there is something with Foam::sigFpe it just mean, usualy that you divide by 0. So, ..... check your boundary conditions (initial values) and you will solve your problem by yourself hope it helps, Cedric vivek05 likes this.

 November 6, 2008, 11:13 juuppiii...I have done...thank #86 emilianyassenov Guest   Posts: n/a juuppiii...I have done...thanks very much Cedric

 November 6, 2008, 11:30 be carefull, It's the begining #87 Senior Member   Cedric DUPRAT Join Date: Mar 2009 Location: Nantes, France Posts: 195 Rep Power: 16 be carefull, It's the begining of a new world now !! :-)

 November 6, 2008, 12:14 Hello Emilian, in addition #88 Member   Sebastian Vogl Join Date: Mar 2009 Location: Munich, Germany Posts: 62 Rep Power: 16 Hello Emilian, in addition to Cedric's post I would also recommend that you take care of the Courant number which you specify in the /system/condrolDict dictionary. In your case (of your post on 3:29) it is 2.22812e+59. It shoud, however be between 0 and 1: 0/system/fvSolution. There is some good information about it in the User Guide chapter 4.4/4.5. Referring to your post of 6:38: A possible mistake (I'm not sure) is, that rhoPimpleFoam is a solver for turbulent flows. So if your flow is laminar you will have to go into the /constant/RASPropertis file and switch of turbulence. As turbulence model you take laminar. This is a dummy-turbulence model. Maybe that could help you. Best regards, Sebastian vivek05 likes this.

 November 7, 2008, 02:30 Thanks Sebastian it is working #89 emilianyassenov Guest   Posts: n/a Thanks Sebastian it is working well... best regards EMO

 February 3, 2009, 10:12 Dear all I am quite new in #90 New Member   Lara Aleluia Reis Join Date: Mar 2009 Posts: 9 Rep Power: 16 Dear all I am quite new in CFD and OpenFoam. I am trying to simulate de dispersion of gas pollutants in the the atmosphere using dieselFoam, with chemistry off. I used the ammonia case that Niklas Nordin post on the forum as a build up point. I have changed the mesh to a much simpler case in witch the inlet is the whole left wall. It seems as though I have been having some problem with the injectorProperties. What I wanted to know to be able to continue is what is really the function of the injectorProperties? Do these properties control the entrance of the pollutants at the inlet? or can I simply define a velocity boundary condition and a concentration at the inlet, and not use the injectorProperties? Thank you in advance Best regards Lara

 February 3, 2009, 10:16 Dear all I am quite new in #91 New Member   Lara Aleluia Reis Join Date: Mar 2009 Posts: 9 Rep Power: 16 Dear all I am quite new in CFD and OpenFoam. I am trying to simulate de dispersion of gas pollutants in the the atmosphere using dieselFoam, with chemistry off. I used the ammonia case that Niklas Nordin post on the forum as a build up point. I have changed the mesh to a much simpler case in witch the inlet is the whole left wall. It seems as though I have been having some problem with the injectorProperties. What I wanted to know to be able to continue is what is really the function of the injectorProperties? Do these properties control the entrance of the pollutants at the inlet? or can I simply define a velocity boundary condition and a concentration at the inlet, and not use the injectorProperties? Thank you in advance Best regards Lara

 February 4, 2009, 03:47 hello again I am running th #92 New Member   Lara Aleluia Reis Join Date: Mar 2009 Posts: 9 Rep Power: 16 hello again I am running the ammonia case from Niklas Nordin and I have changed only the mesh and the controlDict. At time 2.466568e+04 stops here: DILUPBiCG: Solving for h, Initial residual = 0.000156685, Final residual = 5.38316e-10, No Iterations 2 and it keeps giving me this error message attempt to use janafThermo out of temperature range 200 -> 5000; T = 199.974#0 Foam::error::printStack(Foam:stream&) in "/home/aleluia/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libOpenFOAM.so" #1 Foam::error::abort() in "/home/aleluia/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libOpenFOAM.so" #2 Foam::specieThermo >::H(double) const in "/home/aleluia/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libcombustionThermophysic alModels.so" #3 Foam::hMixtureThermo::calcu late() in "/home/aleluia/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libcombustionThermophysic alModels.so" #4 Foam::hMixtureThermo::corre ct() in "/home/aleluia/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libcombustionThermophysic alModels.so" #5 main in "/home/aleluia/OpenFOAM/OpenFOAM-1.5/applications/bin/linuxGccDPOpt/dieselFoam" #6 __libc_start_main in "/lib/tls/i686/cmov/libc.so.6" #7 Foam::regIOobject::writeObject(Foam::IOstream::str eamFormat, Foam::IOstream::versionNumber, Foam::IOstream::compressionType) const in "/home/aleluia/OpenFOAM/OpenFOAM-1.5/applications/bin/linuxGccDPOpt/dieselFoam" From function janafThermo::checkT(const scalar T) const in file /home/dm2/henry/OpenFOAM/OpenFOAM-1.5/src/thermophysicalModels/specie/lnInclude/ janafThermoI.H at line 70. FOAM aborting Aborted can some one help me? regards Lara

 February 4, 2009, 04:30 Thank you sebastian I have #94 New Member   Lara Aleluia Reis Join Date: Mar 2009 Posts: 9 Rep Power: 16 Thank you sebastian I have specified the injection time for a time after my endTime. But I am now getting the error that I have posted before. I'll use your hint to look for the post about distribution of CH4. Thank you very much Lara

 February 4, 2009, 04:36 Hi Lara, referring to your #95 Member   Sebastian Vogl Join Date: Mar 2009 Location: Munich, Germany Posts: 62 Rep Power: 16 Hi Lara, referring to your last post. As you can see in the line: attempt to use janafThermo out of temperature range 200 -> 5000; T = 199.974#0 its a little bit cold within your domain. The properties of your gas species is represented by functions which are based on the coefficients, which are written within your therm.dat file. When you open this file, you will see the collection of numbers for each species. You can see from your error message that the functions for one ore more species are only valid within the given temperature range: 200 -> 5000K So as for an colder environment the coefficients for this functions are not valid, your simulation crashed. You should first restart the simulation from the last written time step and find out whether this error occurs again, then do the post processing and find out why it is so cold in your domain and change the setting so that the temperature doesn't decrease so much. Best regards, Sebastian

 February 4, 2009, 08:42 Hi Sebastian Thank you very #96 New Member   Lara Aleluia Reis Join Date: Mar 2009 Posts: 9 Rep Power: 16 Hi Sebastian Thank you very much for the answer. I have noticed that Temperature was falling out of the given range .. My question is then: Should OpenFoam be using the therm.dat information if chemistry is switch off? Because in my understating, if species are set to not react... how can temperature fall to ~199K? Thank you very much regards Lara

 February 4, 2009, 08:47 Hi Sebastian Thank you very #97 New Member   Lara Aleluia Reis Join Date: Mar 2009 Posts: 9 Rep Power: 16 Hi Sebastian Thank you very much for the answer. I have noticed that Temperature was falling out of the given range .. My question is then: Should OpenFoam be using the therm.dat information if chemistry is switch off? Because in my understating, if species are set to not react... how can temperature fall to ~199K? Thank you very much regards Lara

 February 5, 2009, 05:19 Thank you once again sebastian #99 New Member   Lara Aleluia Reis Join Date: Mar 2009 Posts: 9 Rep Power: 16 Thank you once again sebastian. I'll follow your advise ;) regards Lara

 February 5, 2009, 06:59 Dear all, Please let me kn #100 Member   Hamed Aghajani Join Date: Mar 2009 Location: London, UK Posts: 77 Rep Power: 16 Dear all, Please let me know, How does dieselFoam work? Is it possible to study a LIQUID JET (Liquid Hydrogen) coming out from a simple orifice (d = 1 cm); which its physics involve atomization/spray seconds after spill and vaporization very quickly in atmosphere (flash vaporization) due to its very low boiling point. (droplets will form in atomsphere after coming out from orifice) I run the aachenbomb case with following settings in "injectorProperties" to investigate how the C7H16 Liquid comes out, injectorType unitInjector; diameter 0.005; cd 1; mass 6e-4; temperature 320; nParcels 10000; massFlowRateProfile ( (0 0.12) (0.005 0.12) ); Results in paraFoam do not show any LIQUID JET or liquid particle in system? (How can i track/visualize them?) In case this is not the proper solver, which other ones do you suggest? Best, Hamed