
[Sponsors] 
April 11, 2013, 14:46 
flameletModel 2.1

#161  
New Member
harikrishnan
Join Date: Sep 2012
Posts: 7
Rep Power: 4 
Quote:
Hi Tobi, Thanks for the update i have started using FlameletModel 2.1.x but i i am not able to create PDF library with bigger chemical kinetics file.if the number of species/reactions are more it is giving an error like Clear folders removed `Flamelet_0/Flamelet_1_Hz.out' removed `Flamelet_0/Flamelet_3_Hz.out' removed `Flamelet_0/Flamelet_2_Hz.out' removed `Flamelet_0/Flamelet_5_Hz.out' removed `Flamelet_0/Flamelet_8_Hz.out' removed directory: `Flamelet_0' removed `Logs/LogLookUpTable0' removed `Logs/LogEnthalpydefect_0' removed directory: `Logs' removed `BzzFile.txt' Calculation of flamelet lib with enthalpy defect 0 Finished... Flamelets generated... LookUpTableEnsemble LookUpTable for Flamelets with enthalpy defect 0 terminate called after throwing an instance of 'std:ut_of_range' what(): basic_string::substr ./runFlameletGeneration.sh: line 86: 95530 Aborted ./OpenSMOKE_LookUpTable.sh input "Input/Input_$i.inp" kinetics "$kinetics" buildlibrary > Logs/LogLookUpTable$i LookUpTable for enthalpy defect 0 generated Remove Fluent libs Finished... I dont know whats happening while calculating the look up table. i was using a bigger chemical kinetics file for DME. It is working well with small CRECK kinetics and GRI Mech. please help Thank you Hari 

April 12, 2013, 01:42 
FlameletModel 2.1.x

#162 
New Member
harikrishnan
Join Date: Sep 2012
Posts: 7
Rep Power: 4 
What is the maximum number of species that can be handled by OpenSMOKE_CHEMKINInterpreter.sh
when i run this, i am getting an error terminate called after throwing an instance of 'std:ut_of_range' what(): basic_string::substr ./Run.sh: line 1: 64711 Aborted Please help Thank you Hari 

April 12, 2013, 08:51 

#163 
Senior Member
Bobi
Join Date: Oct 2012
Posts: 286
Rep Power: 5 
Hi hayhari
I experienced the same problem, when I tried to use C1C16 polimi mechanism, Unfortunately I couldn't fix it.I think that due to lack of access to flamelet generation code there might not be a clear answer for this problem. Bobi 

April 14, 2013, 05:33 

#164  
Senior Member
Tobias Holzmann
Join Date: Oct 2010
Location: Leoben (Austria)
Posts: 1,089
Blog Entries: 6
Rep Power: 19 
Quote:
1. all species must be avaiable in the three folders / files > kinetics > thermodynamics > transportproperties 2. As bobi told you, there is not a clear answer for you due to the binaries. Sorry But if there are some users that want to programm with c++ the flamelet generation new I promise that I will work with you! There are some nice c++ libraries for betapdf calculations and so on! 

April 19, 2013, 18:06 

#165 
Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 8,312
Blog Entries: 34
Rep Power: 84 
Greetings to all!
Just passing by to give some news... looks like the CRECK Modeling team has gone pro! Check it out: I found out from their main website: http://creckmodeling.chem.polimi.it/ Best regards, Bruno
__________________


April 20, 2013, 04:37 

#166  
Senior Member
Tobias Holzmann
Join Date: Oct 2010
Location: Leoben (Austria)
Posts: 1,089
Blog Entries: 6
Rep Power: 19 
Quote:
Thanks bruno! 

April 21, 2013, 04:10 
libOpenSMoke

#167 
New Member
harikrishnan
Join Date: Sep 2012
Posts: 7
Rep Power: 4 
Flamelet n.1
DONE! Flamelet n.2 DONE! Flamelet n.3 DONE! Flamelet n.4 DONE! Flamelet n.5 DONE! Flamelet n.6 DONE! Flamelet n.7 DONE! Flamelet n. 1  Scalar dissipation rate (st.): 1 Hz Number Of Species: 98 Number Of Points: 106 terminate called after throwing an instance of 'std:ut_of_range' what(): basic_string::substr Aborted while excecuting OpenSMOKE_LookUpTable.sh.... what is the possible fix for this error..? Thanks Hari 

April 21, 2013, 05:32 

#168 
Senior Member
Tobias Holzmann
Join Date: Oct 2010
Location: Leoben (Austria)
Posts: 1,089
Blog Entries: 6
Rep Power: 19 
Hi Hari,
do you use the flamlet model from my github? If yes I can not give you a solution for that. The problem (i mentioned many posts befor) is the binarys for the falmetlet generation. Its not possible to have a look into that !!! Therefor you are not able to resolve that problem! But maybe there is a othre tool in the laminar tool?! 

April 23, 2013, 09:41 

#170 
New Member
chengyi
Join Date: Jan 2013
Posts: 3
Rep Power: 4 
thank you. the question is solved. I update the linux, and try again. Then it succeeds.


April 29, 2013, 11:42 

#171 
New Member
Join Date: Feb 2013
Posts: 6
Rep Power: 4 
Hi Tobi,
I've posted here some questions about the flamelet libraries one and half month ago. And during that time, I was working on another model first and then on CANTERA to build my own flamelet librairies. Now, I've got the libraries and my next step is to be able to create an interface between the CFD solver and the libraries. I've one and half month to have interesting results... Is it possible? Thanks Danz 

April 29, 2013, 13:52 

#172  
Senior Member
Tobias Holzmann
Join Date: Oct 2010
Location: Leoben (Austria)
Posts: 1,089
Blog Entries: 6
Rep Power: 19 
Quote:
I am interested in that stuff. 

April 30, 2013, 03:34 

#173 
New Member
Join Date: Feb 2013
Posts: 6
Rep Power: 4 
Hi Tobi,
Unfortunately, my cfd code is not free... Now, I'm trying to create a module to read my library 

May 6, 2013, 02:55 
libOpenSmoke

#174 
New Member
harikrishnan
Join Date: Sep 2012
Posts: 7
Rep Power: 4 
Hi,
I have made a a complex geometry in fluent and exported to OpenFoam using Fluent3dToFoam.it is having 20 patches. Then i tried to do the simulation using Flamelet model then i got the following error.. /**\  =========    \\ / F ield  OpenFOAM: The Open Source CFD Toolbox   \\ / O peration  Version: 2.1.1   \\ / A nd  Web: www.OpenFOAM.org   \\/ M anipulation   \**/ Build : 2.1.1221db2718bbb Exec : rhoSimpleFoamFlamelets_1007 Date : May 06 2013 Time : 12:17:23 Host : "haripc" PID : 4800 Case : /home/harikrishnan/libOpenSMOKE/tutorials/rhoSimpleFoamFlamelets/fluent2 nProcs : 1 sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE). fileModificationChecking : Monitoring runtime modified files using timeStampMaster allowSystemOperations : Disallowing usersupplied system call operations // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time Create mesh for time = 0 Reading g Flamelet thermodynamics... Selecting thermodynamics package hPdfThermo<pureMixture<sutherlandTransport<specieT hermo<hConstThermo<perfectGas>>>>> Patch Fixed Temperature 3 Patch Fixed Temperature 4 Patch Fixed Temperature 5 Patch Fixed Temperature 6 Patch Fixed Temperature 7 Patch Fixed Temperature 8 Patch Fixed Temperature 9 Patch Fixed Temperature 10 Patch Fixed Temperature 11 Patch Fixed Temperature 12 Patch Fixed Temperature 13 Patch Fixed Temperature 14 Patch Fixed Temperature 15 Patch Fixed Temperature 16 Patch Fixed Temperature 17 Patch Fixed Temperature 18 Patch Fixed Temperature 19 Patch Fixed Temperature 20 Adiabatic mode... DeltaDirac Distribution for scalar dissipation rate... OpenSMOKE_PDF_Flamelet Index: 1 File name: PDFLibrary/PDF0/SR_1.bin Mixture fraction points: 327 Variance Mixture fraction points: 32 Nominal strain rate: 4.94066e324 Hz Stoic. scalar dissipation rate: 1e06 Hz Temperature min: 292 K Temperature max: 2386.94 K Density min: 0.157304 kg/m3 Density max: 1.1982 kg/m3 OpenSMOKE_PDF_Flamelet Index: 2 File name: PDFLibrary/PDF0/SR_2.bin Mixture fraction points: 327 Variance Mixture fraction points: 32 Nominal strain rate: 4.94066e324 Hz Stoic. scalar dissipation rate: 0.0001 Hz Temperature min: 292 K Temperature max: 2386.93 K Density min: 0.157304 kg/m3 Density max: 1.1982 kg/m3 OpenSMOKE_PDF_Flamelet Index: 3 File name: PDFLibrary/PDF0/SR_3.bin Mixture fraction points: 327 Variance Mixture fraction points: 32 Nominal strain rate: 4.94066e324 Hz Stoic. scalar dissipation rate: 0.01 Hz Temperature min: 292 K Temperature max: 2386.17 K Density min: 0.157356 kg/m3 Density max: 1.1982 kg/m3 OpenSMOKE_PDF_Flamelet Index: 4 File name: PDFLibrary/PDF0/SR_4.bin Mixture fraction points: 327 Variance Mixture fraction points: 32 Nominal strain rate: 4.94066e324 Hz Stoic. scalar dissipation rate: 0.1 Hz Temperature min: 292 K Temperature max: 2380.62 K Density min: 0.157668 kg/m3 Density max: 1.1982 kg/m3 OpenSMOKE_PDF_Flamelet Index: 5 File name: PDFLibrary/PDF0/SR_5.bin Mixture fraction points: 327 Variance Mixture fraction points: 32 Nominal strain rate: 4.94066e324 Hz Stoic. scalar dissipation rate: 1 Hz Temperature min: 292 K Temperature max: 2337.71 K Density min: 0.159747 kg/m3 Density max: 1.1982 kg/m3 OpenSMOKE_PDF_Flamelet Index: 6 File name: PDFLibrary/PDF0/SR_6.bin Mixture fraction points: 327 Variance Mixture fraction points: 32 Nominal strain rate: 4.94066e324 Hz Stoic. scalar dissipation rate: 5 Hz Temperature min: 292 K Temperature max: 2250.72 K Density min: 0.164261 kg/m3 Density max: 1.1982 kg/m3 OpenSMOKE_PDF_Flamelet Index: 7 File name: PDFLibrary/PDF0/SR_7.bin Mixture fraction points: 327 Variance Mixture fraction points: 32 Nominal strain rate: 4.94066e324 Hz Stoic. scalar dissipation rate: 10 Hz Temperature min: 292 K Temperature max: 2197.33 K Density min: 0.167249 kg/m3 Density max: 1.1982 kg/m3 OpenSMOKE_PDF_Flamelet Index: 8 File name: PDFLibrary/PDF0/SR_8.bin Mixture fraction points: 327 Variance Mixture fraction points: 32 Nominal strain rate: 4.94066e324 Hz Stoic. scalar dissipation rate: 20 Hz Temperature min: 292 K Temperature max: 2134.8 K Density min: 0.170962 kg/m3 Density max: 1.1982 kg/m3 OpenSMOKE_PDF_Flamelet Index: 9 File name: PDFLibrary/PDF0/SR_9.bin Mixture fraction points: 327 Variance Mixture fraction points: 32 Nominal strain rate: 4.94066e324 Hz Stoic. scalar dissipation rate: 30 Hz Temperature min: 292 K Temperature max: 2093.9 K Density min: 0.173518 kg/m3 Density max: 1.1982 kg/m3 OpenSMOKE_PDF_Flamelet Index: 10 File name: PDFLibrary/PDF0/SR_10.bin Mixture fraction points: 327 Variance Mixture fraction points: 32 Nominal strain rate: 4.94066e324 Hz Stoic. scalar dissipation rate: 40 Hz Temperature min: 292 K Temperature max: 2062.88 K Density min: 0.175527 kg/m3 Density max: 1.1982 kg/m3 OpenSMOKE_PDF_Flamelet Index: 11 File name: PDFLibrary/PDF0/SR_11.bin Mixture fraction points: 327 Variance Mixture fraction points: 32 Nominal strain rate: 4.94066e324 Hz Stoic. scalar dissipation rate: 100 Hz Temperature min: 292 K Temperature max: 1952.04 K Density min: 0.183232 kg/m3 Density max: 1.1982 kg/m3 OpenSMOKE_PDF_NonAdiabaticFlamelet_Library File name: PDFLibrary/LookUpTable.out Number of Enthalpy defects: 1 Minimum enthalpy defect: 0 kJ/kg Maximum enthalpy defect: 0 kJ/kg Adiabatic flamelets: 1 Temperature fuel: 300 K Temperature oxidizer: 292 K Density fuel: 0.926635 kg/m3 Density oxidizer: 1.1982 kg/m3 Enthalpy fuel: 3.87397e+06 J/kg Enthalpy oxidizer: 115713 J/kg Initialize basic fields... Initialize enthalpy field (field)... pdfThermo h(const scalarField& T, const label patchi)... end.. pdfThermo h(const scalarField& T, const label patchi)... end.. pdfThermo h(const scalarField& T, const label patchi)... end.. pdfThermo h(const scalarField& T, const label patchi)... end.. pdfThermo h(const scalarField& T, const label patchi)... end.. pdfThermo h(const scalarField& T, const label patchi)... end.. pdfThermo h(const scalarField& T, const label patchi)... end.. pdfThermo h(const scalarField& T, const label patchi)... end.. pdfThermo h(const scalarField& T, const label patchi)... end.. pdfThermo h(const scalarField& T, const label patchi)... end.. pdfThermo h(const scalarField& T, const label patchi)... end.. pdfThermo h(const scalarField& T, const label patchi)... end.. pdfThermo h(const scalarField& T, const label patchi)... end.. pdfThermo h(const scalarField& T, const label patchi)... end.. pdfThermo h(const scalarField& T, const label patchi)... end.. pdfThermo h(const scalarField& T, const label patchi)... end.. pdfThermo h(const scalarField& T, const label patchi)... end.. pdfThermo h(const scalarField& T, const label patchi)... end.. pdfThermo h(const scalarField& T, const label patchi)... end.. pdfThermo h(const scalarField& T, const label patchi)... end.. pdfThermo h(const scalarField& T, const label patchi)... end.. Calculate (first time)... #0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #2 in "/lib/x86_64linuxgnu/libc.so.6" #3 Foam::hPdfThermo<Foam:ureMixture<Foam::sutherlan dTransport<Foam::specieThermo<Foam::hConstThermo<F oam:erfectGas> > > > >::calculate() in "/home/harikrishnan/OpenFOAM/harikrishnan2.1.1/platforms/linux64GccDPOpt/lib/libOpenSMOKE_pdfThermo.so" #4 Foam::hPdfThermo<Foam:ureMixture<Foam::sutherlan dTransport<Foam::specieThermo<Foam::hConstThermo<F oam:erfectGas> > > > >::hPdfThermo(Foam::fvMesh const&) in "/home/harikrishnan/OpenFOAM/harikrishnan2.1.1/platforms/linux64GccDPOpt/lib/libOpenSMOKE_pdfThermo.so" #5 Foam::basicPdfThermo::addfvMeshConstructorToTable< Foam::hPdfThermo<Foam:ureMixture<Foam::sutherlan dTransport<Foam::specieThermo<Foam::hConstThermo<F oam:erfectGas> > > > > >::New(Foam::fvMesh const&) in "/home/harikrishnan/OpenFOAM/harikrishnan2.1.1/platforms/linux64GccDPOpt/lib/libOpenSMOKE_pdfThermo.so" #6 Foam::basicPdfThermo::New(Foam::fvMesh const&) in "/home/harikrishnan/OpenFOAM/harikrishnan2.1.1/platforms/linux64GccDPOpt/lib/libOpenSMOKE_pdfThermo.so" #7 in "/home/harikrishnan/OpenFOAM/harikrishnan2.1.1/platforms/linux64GccDPOpt/bin/rhoSimpleFoamFlamelets_1007" #8 __libc_start_main in "/lib/x86_64linuxgnu/libc.so.6" #9 in "/home/harikrishnan/OpenFOAM/harikrishnan2.1.1/platforms/linux64GccDPOpt/bin/rhoSimpleFoamFlamelets_1007" Floating point exception (core dumped) I am relatively new to openFoam and not able to identify the reason for this error. Also H in 0 directory is which enthalpy? enthalpy of formation?? Please help. Thanks Hari 

May 23, 2013, 07:20 

#175 
Senior Member
Tobias Holzmann
Join Date: Oct 2010
Location: Leoben (Austria)
Posts: 1,089
Blog Entries: 6
Rep Power: 19 
Hi Hari,
that is a Setup mistake. It seems that some BC has no "value uniform x;" Check out your BC!! If you wanna know if that Problem occure due to the BC you can insert a line into the hPdfThermo.C (i think that is the file) > function calculate(); Tobi 

May 25, 2013, 03:44 

#176 
Senior Member
Bobi
Join Date: Oct 2012
Posts: 286
Rep Power: 5 
Hi Tobi
Thanks very much for your quick reply. When I tried simulating DLR flame, the compiler got a error from dimension of the ChiUpdate equation. Although I checked the formula from your modified guide and the peters turbulent combustion book and the same equations were written there. However, from dimension error it seems that the R.H.S. must be divided by rho. If you remember, Once I asked you about a missing rho. I think now it makes sense. Again thanks for your quick and kind reply Regards Bobi 

May 25, 2013, 04:29 

#178 
Senior Member
Tobias Holzmann
Join Date: Oct 2010
Location: Leoben (Austria)
Posts: 1,089
Blog Entries: 6
Rep Power: 19 
[quote="Tobis  Masterthesis"]
The turbulent diffusion term is modelled with the approximation of gradient diffusion: So you see that with RANS Modelling the density is in the mu_t and mu_t is modelled with the turbulence models! (k/eps). In LES (I think) you have a other modellation! So it could be the reason in your calculation (missing rho). But I have no Idea how LES is working  so try out and read litrerature Last edited by Tobi; May 27, 2013 at 06:14. 

May 27, 2013, 10:31 

#180 
Senior Member
Bobi
Join Date: Oct 2012
Posts: 286
Rep Power: 5 
Hi Tobi
Now, I am trying to solve my problem with fireFoam (flamelet/progress variable). Very soon, I will test the LESFlamelet code and will tell you the results. Regards Bobi 

Thread Tools  
Display Modes  

