
[Sponsors] 
March 12, 2013, 15:31 

#141 
New Member
Join Date: Feb 2013
Posts: 6
Rep Power: 6 
Hi Tobi,
Thank you for your message! In my case, I will use cantera and I will start from the EDC. From now, I just begun my internship so I study the literature and then I will let you know how does it will go... 

March 18, 2013, 03:22 
openfoam 2.2.x

#142 
Super Moderator
Tobias Holzmann
Join Date: Oct 2010
Location: Leoben (Austria)
Posts: 1,766
Blog Entries: 6
Rep Power: 31 
Hi all,
I tryed to use the flameletmodel in openfoam 2.2.x but there are some problems with compiling the thermodynamics. I ll try To rebuild the libraries for 2.2.x and Bild the instationary solver "turbulentFlameletRhoPimpelFoam". Tobi 

March 18, 2013, 13:00 

#143 
Senior Member
Bobi
Join Date: Oct 2012
Location: Chicago, Illinois
Posts: 403
Rep Power: 8 
Dear Tobi
Hi I hope u will succeed in this route, however I am working on the code too. I have also worked with firefoam solver in 2.1.x . This solver uses flame surface density model which is a flamelet model that instead of integrating using PDF, employs some kind of sigma formulation (according to Tarek echekki book Turbulent Combustion ModelingAdvances and new trends) At the moment I do not have enough time to check the new firefoam solver (I will check it later). Nevertheless, I have read in openfoam site that in the new solver, multistep combustion has added to the firefoam. in the prior version, it only had single step. Also the cfd main solver is RhoPimple in this solver. Regards Bobi 

March 23, 2013, 12:15 

#145 
Super Moderator
Tobias Holzmann
Join Date: Oct 2010
Location: Leoben (Austria)
Posts: 1,766
Blog Entries: 6
Rep Power: 31 
Hi all,
I am working the last day on rebuilding the thermomodel for 2.2.x. I added two features for output the flamelet properties (two new variables in the flameletProperties) and cleared the thermocode (named flameletThermo) Here is an output at the beginning: Code:
/**\  =========    \\ / F ield  OpenFOAM: The Open Source CFD Toolbox   \\ / O peration  Version: 2.2.x   \\ / A nd  Web: www.OpenFOAM.org   \\/ M anipulation   \**/ Build : 2.2.x020332035ede Exec : flamaletSimpleFoam Date : Mar 23 2013 Time : 17:12:28 Host : "cfd" PID : 14013 Case : /home/shorty/OpenFOAM/shorty2.2.x/run/SANDIA_CH4H2N2 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 SIMPLE: convergence criteria field p tolerance 0.01 field U tolerance 0.0001 field h tolerance 0.001 field "(kepsilonomega)" tolerance 0.001 Reading thermophysical properties Selecting thermodynamics package { type pdfFlameletThermo; mixture pureMixture; transport sutherland; thermo hConst; equationOfState perfectGas; specie specie; energy sensibleEnthalpy; } /**\  Flamelet Thermo initialisation     Rebuild by Tobias Holzmann M.Eng.   www.Holzmanncfd.de  \**/ Enthalpy: + inletfuel <fixedValue> + inletair <fixedValue> + outlet + axis + sidewall + burnerwall + front + back Temperature: + inletfuel <fixedValue> + inletair <fixedValue> + outlet + axis + sidewall + burnerwall + front + back Mixture fraction: + inletfuel <fixedValue> + inletair <fixedValue> + outlet + axis + sidewall + burnerwall + front + back Flamelet thermo mode is <nonadiabatic> Flamelet thermo uses the <deltadirac distribution> for the scalar dissipation rate Flamelet thermo shows all single flamelet properties Flamelet thermo shows all flamelet library properties Flamelet property: + Index: 1 + File name: PDFLibrary/PDF0/SR_1.bin + Mixture fraction points: 44 + Variance Mixture fraction points: 30 + Nominal strain rate: 6.91521e310 Hz + Stoic. scalar dissipation rate: 1e06 Hz + Temperature min: 292.15 K + Temperature max: 2124.88 K + Density min: 0.152907 kg/m3 + Density max: 1.20255 kg/m3 Flamelet property: + Index: 2 + File name: PDFLibrary/PDF0/SR_2.bin + Mixture fraction points: 44 + Variance Mixture fraction points: 30 + Nominal strain rate: 6.91521e310 Hz + Stoic. scalar dissipation rate: 1 Hz + Temperature min: 292.15 K + Temperature max: 2033.91 K + Density min: 0.157558 kg/m3 + Density max: 1.20255 kg/m3 Flamelet property: + Index: 3 + File name: PDFLibrary/PDF0/SR_3.bin + Mixture fraction points: 44 + Variance Mixture fraction points: 30 + Nominal strain rate: 6.91521e310 Hz + Stoic. scalar dissipation rate: 5 Hz + Temperature min: 292.15 K + Temperature max: 1989.67 K + Density min: 0.159646 kg/m3 + Density max: 1.20255 kg/m3 Flamelet property: + Index: 4 + File name: PDFLibrary/PDF0/SR_4.bin + Mixture fraction points: 44 + Variance Mixture fraction points: 30 + Nominal strain rate: 6.91521e310 Hz + Stoic. scalar dissipation rate: 30 Hz + Temperature min: 292.15 K + Temperature max: 1923.46 K + Density min: 0.163351 kg/m3 + Density max: 1.20255 kg/m3 Flamelet property: + Index: 5 + File name: PDFLibrary/PDF0/SR_5.bin + Mixture fraction points: 44 + Variance Mixture fraction points: 30 + Nominal strain rate: 6.91521e310 Hz + Stoic. scalar dissipation rate: 75 Hz + Temperature min: 292.15 K + Temperature max: 1868.15 K + Density min: 0.167195 kg/m3 + Density max: 1.20255 kg/m3 Flamelet property: + Index: 6 + File name: PDFLibrary/PDF0/SR_6.bin + Mixture fraction points: 44 + Variance Mixture fraction points: 30 + Nominal strain rate: 6.91521e310 Hz + Stoic. scalar dissipation rate: 200 Hz + Temperature min: 292.15 K + Temperature max: 1760.13 K + Density min: 0.176681 kg/m3 + Density max: 1.20255 kg/m3 Flamelet property: + Index: 7 + File name: PDFLibrary/PDF0/SR_7.bin + Mixture fraction points: 44 + Variance Mixture fraction points: 30 + Nominal strain rate: 6.91521e310 Hz + Stoic. scalar dissipation rate: 800 Hz + Temperature min: 292.15 K + Temperature max: 292.15 K + Density min: 0.699245 kg/m3 + Density max: 1.20255 kg/m3 + Cold flame Flamelet library property: + File name: PDFLibrary/PDF0/LookUpTable.out + Number of Flamelets: 6 + Minimum stoic. scalar diss. rate: 1e06 Hz + Maximum stoic. scalar diss. rate: 200 Hz + Temperature min: 292.15 K + Temperature max: 2124.88 K + Density min: 0.152907 kg/m3 + Density max: 1.20255 kg/m3 + Temperature fuel: 292.15 K + Temperature oxidizer: 292.15 K + Density fuel: 0.699245 kg/m3 + Density oxidizer: 1.20255 kg/m3 + Enthalpy fuel: 998143 J/kg + Enthalpy oxidizer: 15974.1 J/kg /**\ > Finished loading LookUpTable number: 1 out of 5 < \**/ . . . . . . Flamelet library property: + File name: PDFLibrary/PDF500/LookUpTable.out + Number of Flamelets: 6 + Minimum stoic. scalar diss. rate: 1e06 Hz + Maximum stoic. scalar diss. rate: 200 Hz + Temperature min: 292.15 K + Temperature max: 1816.43 K + Density min: 0.17916 kg/m3 + Density max: 1.20255 kg/m3 + Temperature fuel: 292.15 K + Temperature oxidizer: 292.15 K + Density fuel: 0.699245 kg/m3 + Density oxidizer: 1.20255 kg/m3 + Enthalpy fuel: 998143 J/kg + Enthalpy oxidizer: 15974.1 J/kg /**\ > Finished loading LookUpTable number: 5 out of 5 < \**/ Flamelet thermo loaded all libraries Summarize the whole flamelet library: + File name: PDFLibrary/LookUpTable.out + Number of Enthalpy defects: 5 + Minimum enthalpy defect: 0 kJ/kg + Maximum enthalpy defect: 500 kJ/kg + Adiabatic flamelets: 1 + Temperature fuel: 292.15 K + Temperature oxidizer: 292.15 K + Density fuel: 0.699245 kg/m3 + Density oxidizer: 1.20255 kg/m3 + Enthalpy fuel: 998143 J/kg + Enthalpy oxidizer: 15974.1 J/kg 

March 24, 2013, 09:39 

#146 
Super Moderator
Tobias Holzmann
Join Date: Oct 2010
Location: Leoben (Austria)
Posts: 1,766
Blog Entries: 6
Rep Power: 31 
Hi all,
the new thermodynamic for the flamelet model is build in 2.2.x I finished code developement at the moment and it seems that everything is working. Tomorrow I ll validate the solver and if everything is working fine I ll share it with you. 

March 24, 2013, 12:12 

#147 
Senior Member
Bobi
Join Date: Oct 2012
Location: Chicago, Illinois
Posts: 403
Rep Power: 8 
Dear Tobi
Hi Love u man. I wish I were in deutschland and could buy u a beer. Did u have any success on pimple algorithm? Bobi 

March 25, 2013, 06:21 

#148 
Super Moderator
Tobias Holzmann
Join Date: Oct 2010
Location: Leoben (Austria)
Posts: 1,766
Blog Entries: 6
Rep Power: 31 
Hi all,
the new flamelet model is avaiable at my homepage: http://www.holzmanncfd.de/index.php...openfoam22x (not finished yet) or directly on my github: https://github.com/shorty/flameletModel2.2.x + I added a very very brief description (have no time to build a very good one  its all my spare time ) + I added a validation doc (I ll add this in the validation forum too) + Validation: http://www.cfdonline.com/Forums/ope...tml#post416171 @bobi: yesterday I tryed to build the flameletPimpleFoam solver but without access. I think the problem is the pEqn.H file. I realized that I have to use the old pEqn.H file from 1.7.x for getting a correct result with the flamelet thermophysical model. Last edited by Tobi; March 25, 2013 at 07:59. 

March 28, 2013, 11:38 

#149 
Senior Member

Hi Tobi,
Thanks a lot for all the work. I had some problems with compiling, this may be related to the relatively large amounts of updates on the 2.2.x repository. Anyhow for some reason the "basic" folder did not compile at first, I had to change the file Make/option into: Code:
EXE_INC = \ I$(LIB_SRC)/finiteVolume/lnInclude \ I$(LIB_SRC)/thermophysicalModels/specie/lnInclude \ I$(LIB_SRC)/thermophysicalModels/basic/lnInclude \ I$(LIB_SRC)/meshTools/lnInclude \ I$(realpath ../flameletExtraction/flamelets/turbulent/lnInclude) \ I$(realpath ../flameletExtraction/common/lnInclude) LIB_LIBS = \ lfiniteVolume \ lfluidThermophysicalModels Code:
I$(LIB_SRC)/thermophysicalModels/basic/lnInclude \ Code:
lfluidThermophysicalModels And for the solver I changed the last line in the EXE_INC part into: Code:
I$(realpath ../../../../thermophysicalModels/basic/flameletThermo/) At least everything compiled and the tutorial did run, but I am not entirely sure if everything is now setup correctly. Once again thanks for your effort, much appreciated. 

March 29, 2013, 06:31 

#150 
Super Moderator
Tobias Holzmann
Join Date: Oct 2010
Location: Leoben (Austria)
Posts: 1,766
Blog Entries: 6
Rep Power: 31 
Hi Tom,
I ll have a look at that till monday. I ll compile OF2.2.x on my laptop and will try to use the flamelet model for 2.2.x on that computer. Maybe there are some wrong or missing entrys in the Make folder like you told. If your solver is working everything should work Tobi 

March 29, 2013, 06:50 

#151 
Senior Member

Hi Tobi,
Yes the solver did work and the result looked reasonable. Cheers, happy Easter, Tom 

March 29, 2013, 09:05 

#152  
New Member
chengyi
Join Date: Jan 2013
Posts: 3
Rep Power: 6 
Hi, all
I am learning to the libopensmoke. but when I want to try to compile the libs and solver , some problem happen and not succeed to compile the solver, Do you know what happen? Thank you Quote:
Last edited by lcy; March 29, 2013 at 09:21. Reason: when I want to try to compile the solver, some problem happen and not succeed to compile the solver, Do you know what happen? 

March 29, 2013, 11:31 

#154 
New Member
chengyi
Join Date: Jan 2013
Posts: 3
Rep Power: 6 

March 29, 2013, 11:33 

#155 
Super Moderator
Tobias Holzmann
Join Date: Oct 2010
Location: Leoben (Austria)
Posts: 1,766
Blog Entries: 6
Rep Power: 31 
Your solver should work and may be the soot solver is not working but I never used used the 1.7.0 Version and this version has some bugs in it. You should use the latest one!
But its your decision. 

April 1, 2013, 03:12 

#156  
New Member
harikrishnan
Join Date: Sep 2012
Posts: 7
Rep Power: 6 
Quote:
Dear tobi i am facing this problem now. Did u find the solution for this? when i updated GSL_INST_DIR i an getting a new error rhoSimpleFoamFlamelets_1007: error while loading shared libraries: libOpenSMOKE_pdfThermo.so: cannot open shared object file: No such file or directory Please help me solve this issue. Thanks Hari 

April 1, 2013, 05:49 

#157 
Super Moderator
Tobias Holzmann
Join Date: Oct 2010
Location: Leoben (Austria)
Posts: 1,766
Blog Entries: 6
Rep Power: 31 
Hi Harry,
1. you do not have to compile/install GSL 2. use the FlameletModel 2.1.x on my git repository. Its newer and some userunfriendly setups are made there. 3. if you have the possibility use OpenFOAM 2.2.x couse I made a few output changes in the FlameletModel 2.2.x https://github.com/shorty?tab=repositories 

April 1, 2013, 07:23 

#158  
Super Moderator
Tobias Holzmann
Join Date: Oct 2010
Location: Leoben (Austria)
Posts: 1,766
Blog Entries: 6
Rep Power: 31 
Quote:
I compiled the libs and solver on my laptop now. There were some problems like you told. But the additional libFluid... i did not have to add Well I updated the flamelet libraries on github and now it should work Thanks for your feedback 

April 3, 2013, 04:44 

#159 
Senior Member
Bobi
Join Date: Oct 2012
Location: Chicago, Illinois
Posts: 403
Rep Power: 8 
Hi Tobi
I have a question about the code results. As u know, the code gives the mass fraction of any specie that we want. However, for computing emission index I need the mole fraction of CO and CO2 at a certain position. I know that in literature we have a correlation for commuting mass to mole fraction. But, as u know in that correlation in addition to the molecular weight of the considering specie we need to have the molecular weight of the mixture which is now the products molecular weight. I think since the code computes species mass fractions, therefore it calculates the molecular weight of products as a whole mixture. So, would you plz hint me how can I exploit CO and CO2 mole fractions from the results or maybe change the code to get mole fractions data? Regards Bobi Last edited by babakflame; April 3, 2013 at 05:11. 

April 3, 2013, 06:14 

#160  
Super Moderator
Tobias Holzmann
Join Date: Oct 2010
Location: Leoben (Austria)
Posts: 1,766
Blog Entries: 6
Rep Power: 31 
Quote:
you are right. The code gives you the mass fraction. In my thesis I used the important species and the molecular weight for that to calculate the mole fraction in paraview. The problem is, that the flameletsolver and thermodynamic is not calculating the mass fraction. The thermodynamics just look for the mass fraction into the LookUpTables and there for it was not easy for me to build a new function that is calculating that. It 's possible to do that in the thermodynamics, so that you build new functions or stuff like that. The problem is here, that you have do many things: 1. the thermodynamics just realize the species you want to output: Code:
// Initialise the wanted mass fraction fields omega_.setSize(flamelets_library.number_of_species()+1); for (int j=0;j<flamelets_library.number_of_species()+1;j++) { if(j < flamelets_library.number_of_species()) { std::string name_of_species = "omega_" + flamelets_library.species()[j+1]; omega_.set ( j, new volScalarField ( IOobject ( name_of_species, mesh.time().timeName(), mesh, IOobject::NO_READ, IOobject::AUTO_WRITE), mesh, dimensionedScalar("zero", dimensionSet(0,0,0,0,0,0,0), 0.0 ) ) ); } } 3. its annoying that we have no chance to have a look into the flamelet generatin code So what you have to do is: 1. set up a function that get all involved species for you case 2. set up a list of constants (molecular weights of ALL species) 3. building a new function calculate_molefraction() (the simplest thing) Well maybe the mole fraction is also avaiable in the lookuptable but I dont know how to have a look into that. If there are the mole fractions included, it would be the easiest way to get the values. But as I told befor... its annoying that we dont have the chance to have a look into the flamelet generation! I would build a own free one but I have no time for that. 

Thread Tools  
Display Modes  

