[hayhari]

Originally Posted by Tobi

Hi Harry,

1. you do not have to compile/install GSL
2. use the Flamelet-Model 2.1.x on my git repository. Its newer and some user-unfriendly setups are made there.

3. if you have the possibility use OpenFOAM 2.2.x couse I made a few output changes in the Flamelet-Model 2.2.x

https://github.com/shor-ty?tab=repositories

Hi Tobi,

Thanks for the update i have started using Flamelet-Model 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/Log-LookUpTable0'
removed `Logs/Log-Enthalpydefect_0'
removed directory: `Logs'
removed `BzzFile.txt'
Calculation of flamelet lib with enthalpy defect 0
Finished...
Flamelets generated...
LookUpTable-Ensemble
Look-Up-Table 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" -build-library > Logs/Log-LookUpTable$i
Look-Up-Table 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

[hayhari]

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

[babakflame]

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

[Tobi]

Originally Posted by hayhari

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

Normally everything should work:

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 beta-pdf- calculations and so on!

[wyldckat]

Greetings to all!

Just passing by to give some news... looks like the CRECK Modeling team has gone pro! Check it out:

• http://www.opensmoke.polimi.it/
• http://www.catalyticfoam.polimi.it/

I found out from their main website: http://creckmodeling.chem.polimi.it/

Best regards,
Bruno

[Tobi]

Originally Posted by wyldckat

Greetings to all!

Just passing by to give some news... looks like the CRECK Modeling team has gone pro! Check it out:

• http://www.opensmoke.polimi.it/
• http://www.catalyticfoam.polimi.it/

I found out from their main website: http://creckmodeling.chem.polimi.it/

Best regards,
Bruno

Sounds and looks good
Thanks bruno!

[hayhari]

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

[Tobi]

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?!

[Tobi]

Well the laminarSmoke uses a binary execution file too

[lcy]

thank you. the question is solved. I update the linux, and try again. Then it succeeds.

[danz]

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

[Tobi]

Originally Posted by danz

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

Yes the implementation of your Flamelets into the OF Solver is very easy. Is your code free ?

I am interested in that stuff.

[danz]

Hi Tobi,

Unfortunately, my cfd code is not free...
Now, I'm trying to create a module to read my library

[hayhari]

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.1-221db2718bbb
Exec : rhoSimpleFoamFlamelets_1007
Date : May 06 2013
Time : 12:17:23
Host : "hari-pc"
PID : 4800
Case : /home/harikrishnan/libOpenSMOKE/tutorials/rhoSimpleFoamFlamelets/fluent2
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied 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...
Delta-Dirac Distribution for scalar dissipation rate...

OpenSMOKE_PDF_Flamelet
Index: 1
File name: PDF-Library/PDF-0/SR_1.bin
Mixture fraction points: 327
Variance Mixture fraction points: 32
Nominal strain rate: 4.94066e-324 Hz
Stoic. scalar dissipation rate: 1e-06 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: PDF-Library/PDF-0/SR_2.bin
Mixture fraction points: 327
Variance Mixture fraction points: 32
Nominal strain rate: 4.94066e-324 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: PDF-Library/PDF-0/SR_3.bin
Mixture fraction points: 327
Variance Mixture fraction points: 32
Nominal strain rate: 4.94066e-324 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: PDF-Library/PDF-0/SR_4.bin
Mixture fraction points: 327
Variance Mixture fraction points: 32
Nominal strain rate: 4.94066e-324 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: PDF-Library/PDF-0/SR_5.bin
Mixture fraction points: 327
Variance Mixture fraction points: 32
Nominal strain rate: 4.94066e-324 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: PDF-Library/PDF-0/SR_6.bin
Mixture fraction points: 327
Variance Mixture fraction points: 32
Nominal strain rate: 4.94066e-324 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: PDF-Library/PDF-0/SR_7.bin
Mixture fraction points: 327
Variance Mixture fraction points: 32
Nominal strain rate: 4.94066e-324 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: PDF-Library/PDF-0/SR_8.bin
Mixture fraction points: 327
Variance Mixture fraction points: 32
Nominal strain rate: 4.94066e-324 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: PDF-Library/PDF-0/SR_9.bin
Mixture fraction points: 327
Variance Mixture fraction points: 32
Nominal strain rate: 4.94066e-324 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: PDF-Library/PDF-0/SR_10.bin
Mixture fraction points: 327
Variance Mixture fraction points: 32
Nominal strain rate: 4.94066e-324 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: PDF-Library/PDF-0/SR_11.bin
Mixture fraction points: 327
Variance Mixture fraction points: 32
Nominal strain rate: 4.94066e-324 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: PDF-Library/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_64-linux-gnu/libc.so.6"
#3 Foam::hPdfThermo<Foam:ureMixture<Foam::sutherlan dTransport<Foam::specieThermo<Foam::hConstThermo<F oam:erfectGas> > > > >::calculate() in "/home/harikrishnan/OpenFOAM/harikrishnan-2.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/harikrishnan-2.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/harikrishnan-2.1.1/platforms/linux64GccDPOpt/lib/libOpenSMOKE_pdfThermo.so"
#6 Foam::basicPdfThermo::New(Foam::fvMesh const&) in "/home/harikrishnan/OpenFOAM/harikrishnan-2.1.1/platforms/linux64GccDPOpt/lib/libOpenSMOKE_pdfThermo.so"
#7
in "/home/harikrishnan/OpenFOAM/harikrishnan-2.1.1/platforms/linux64GccDPOpt/bin/rhoSimpleFoamFlamelets_1007"
#8 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#9
in "/home/harikrishnan/OpenFOAM/harikrishnan-2.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

[Tobi]

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

[babakflame]

Hi Tobi

Thanks very much for your quick reply. When I tried simulating DLR flame, the compiler got a error from dimension of the Chi-Update 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

[Tobi]

...


Removed due to the reason in the edit line!

[Tobi]

[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

[Tobi]

Is it working now?

How is your eqn. for Z ? Especially the Diffusion Terms!

[babakflame]

Hi Tobi

Now, I am trying to solve my problem with fireFoam (flamelet/progress variable).
Very soon, I will test the LES-Flamelet code and will tell you the results.

Regards
Bobi