CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Community Contributions

[openSmoke] libOpenSMOKE

Register Blogs Community New Posts Updated Threads Search

Like Tree133Likes

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   February 4, 2015, 15:32
Default
  #441
New Member
 
serkan
Join Date: Dec 2014
Posts: 8
Rep Power: 12
Sermengi is on a distinguished road
Hi Tobi,

thanks for your effort about the flamelet model. We are trying to simulate the Sandia D flame using the flamelet model. Now we are creating the adiabatic flamelet tables using the libOpenSMOKE. We modified the Data.inp according to the properties of the Sandia flame D. Also we use GRI Mech 3.0 for the chemical kinetics. We didn't change the operations.inp file. I am adding content of these files to my massage.

According to our input file, we found that flamelets are quenched when scalar dissipation rate is equal to 572 Hz. On the other hand, flamelets are quenched for scalar dissipation rate is equal to 433 Hz in the tables located in tutorial case of another program. There is a huge difference between the scalar dissipation rates and we can not figure out where the mistake is. We will be so happy if you give us a hint about this situation. Thanks four your help.

Best regards,
Serkan.
Attached Files
File Type: pdf Sandia D flame.pdf (20.4 KB, 86 views)
Sermengi is offline   Reply With Quote

Old   February 4, 2015, 20:59
Default
  #442
Senior Member
 
Freedom
Join Date: May 2014
Posts: 209
Rep Power: 13
wenxu is on a distinguished road
I think the different scalar dissipation rate lead to quench means that the flame quenched at the different scalar dissipation rate values but not means that flame START to quenched!!! For 572 Hz, you may find that 550 may be quenched also.

I think FlameMaster will be a good tool to generate flamelet library. There you will find if 573HZ quenched but 572HZ not quenched, then it will calculate the flamelet equation with scalar dissipation rate between 573--572HZ.

regards,
wenxu
wenxu is offline   Reply With Quote

Old   February 5, 2015, 05:34
Default Bug
  #443
Super Moderator
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Bad Wörishofen
Posts: 2,711
Blog Entries: 6
Rep Power: 52
Tobi has a spectacular aura aboutTobi has a spectacular aura aboutTobi has a spectacular aura about
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
Dear all,

I am working to reorganize the flameletModel (libOpenSMOKE) and implement LES as soon as possible. Also I realized, that no gravity term was included in momentum or pressure equation which could lead to the problem of to high velocity prediction at the centerlines compared to my validations.
__________________
Keep foaming,
Tobias Holzmann
Tobi is offline   Reply With Quote

Old   February 5, 2015, 05:38
Default
  #444
Super Moderator
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Bad Wörishofen
Posts: 2,711
Blog Entries: 6
Rep Power: 52
Tobi has a spectacular aura aboutTobi has a spectacular aura aboutTobi has a spectacular aura about
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
Quote:
Originally Posted by Sermengi View Post
Hi Tobi,

thanks for your effort about the flamelet model. We are trying to simulate the Sandia D flame using the flamelet model. Now we are creating the adiabatic flamelet tables using the libOpenSMOKE. We modified the Data.inp according to the properties of the Sandia flame D. Also we use GRI Mech 3.0 for the chemical kinetics. We didn't change the operations.inp file. I am adding content of these files to my massage.

According to our input file, we found that flamelets are quenched when scalar dissipation rate is equal to 572 Hz. On the other hand, flamelets are quenched for scalar dissipation rate is equal to 433 Hz in the tables located in tutorial case of another program. There is a huge difference between the scalar dissipation rates and we can not figure out where the mistake is. We will be so happy if you give us a hint about this situation. Thanks four your help.

Best regards,
Serkan.

Hi,

this is old topic -> http://sourceforge.net/projects/open...5.pdf/download
What I ment is that this is already done and should be work properly.
But notice that I found a bug (no gravity is included - I am so stupid that I did not found this the last years )
__________________
Keep foaming,
Tobias Holzmann

Last edited by Tobi; February 5, 2015 at 08:27.
Tobi is offline   Reply With Quote

Old   February 5, 2015, 11:26
Default
  #445
New Member
 
serkan
Join Date: Dec 2014
Posts: 8
Rep Power: 12
Sermengi is on a distinguished road
Thanks for your quick replies and suggestions, we will discuss the results again.

Best regards,
Serkan.
Sermengi is offline   Reply With Quote

Old   February 9, 2015, 06:51
Default
  #446
Senior Member
 
Bobby
Join Date: Oct 2012
Location: Michigan
Posts: 454
Rep Power: 16
babakflame is on a distinguished road
Greetings Tobi

Thanks for sending me your comparison file.

That was a very nice and yet brief work describing flamelet solvers procedure.

Just one point: The LES solver which I made based on your hints and helps exits, there is just a requirement for working on PISO algorithm to be developed to PIMPLE.

I can send you my final version. So it can be tested by you and others interested.


Best,
Bobi
babakflame is offline   Reply With Quote

Old   February 9, 2015, 07:00
Default
  #447
Super Moderator
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Bad Wörishofen
Posts: 2,711
Blog Entries: 6
Rep Power: 52
Tobi has a spectacular aura aboutTobi has a spectacular aura aboutTobi has a spectacular aura about
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
Yes please send it to me!
Yesterday I started to develop the libOpenSMOKE tool with FOAM classes in a more readability way.
I think it could be a very nice library which could be included in the standard OpenFOAM (after I finished it).

Looking forward to finish it.
__________________
Keep foaming,
Tobias Holzmann
Tobi is offline   Reply With Quote

Old   February 11, 2015, 05:34
Default
  #448
Senior Member
 
Freedom
Join Date: May 2014
Posts: 209
Rep Power: 13
wenxu is on a distinguished road
Anyone on this thread is familiar with FlameMaster? I can not get converged results if fuel is diluted. The fuel side boundary condition is as follows:

Fuel Side {
dirichlet {
t = 300.0
y->N2 = 0.5825
y->CH4 = 0.0446
y->CO = 0.348
y->C2H2 = 0.0249
}
}


Could anyone give me some hints? Thank you in advance.

regards,
wenxu
wenxu is offline   Reply With Quote

Old   February 11, 2015, 07:00
Default
  #449
Super Moderator
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Bad Wörishofen
Posts: 2,711
Blog Entries: 6
Rep Power: 52
Tobi has a spectacular aura aboutTobi has a spectacular aura aboutTobi has a spectacular aura about
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
Dear Wenxu,

this thread is about libOpenSMOKE not Flame-Master.
Please do not mix it up. The thread should be straight-forward to that topic.

If you have questions with Flame-Master, you should open a new thread.

PS: I used it once but I am not familiar with Flame-Master.
__________________
Keep foaming,
Tobias Holzmann
Tobi is offline   Reply With Quote

Old   February 11, 2015, 07:29
Default
  #450
Senior Member
 
Freedom
Join Date: May 2014
Posts: 209
Rep Power: 13
wenxu is on a distinguished road
sorry for that.
wenxu is offline   Reply With Quote

Old   February 19, 2015, 13:52
Default
  #451
Member
 
Join Date: Feb 2015
Posts: 31
Rep Power: 11
Stefano Puggelli is on a distinguished road
Hi everybody,
I am a PhD student from University of Florence and I am now beginnig to work on combustion modeling in OpenFOAM. My first goal would be to create and validate an FGM model in FOAM for partially premixed combustion.
First of all I made a review of the model already implemented in OF (xiFoam, reactingFoam, fireFoam..) and then I am now trying to understand exactly the available flamelet model as the one built by Tobias or the one realized by Hagen Muller.
I am writing here since I would like to know if someone can give me more informations about the capabilities of these model in order to start my development program. I understand the theory behind for example the flameletSimpleFoam of flameletFoam solver, but I would like to know if you can give more informations on that models (validations, bugs, LES application test cases) considering also the FGM application
I am very interested obviously about the flamelets models but also on fireFoam model. On this second one I don't understand exactly the coupling between a mixture fraction transport equation and FSD model.. it would be great if I clear me these aspects!
Anyway I am completely willing for working together on this items!
Stefano Puggelli is offline   Reply With Quote

Old   February 19, 2015, 14:02
Default
  #452
Senior Member
 
Bobby
Join Date: Oct 2012
Location: Michigan
Posts: 454
Rep Power: 16
babakflame is on a distinguished road
Greetings Stefano

There is an available FGM model in O.F. written by Hans Kruger. refer to:
http://www.cfd-online.com/Forums/ope...?highlight=FGM

About fireFoam:
AFA I know based on O.F. 2.2.x (I have not tested OF 2.3.x), FSD model has some deficiencies featuring:
  • It only includes single step reaction
  • For incorporating fuels other than propane and methane you have to find the coefficients required (not easy to find)
  • If your fuel stream is a hybrid one, the model does not support it
For understanding better the implementation of FSD model refer to:



http://www.tandfonline.com/doi/abs/1...1/009#abstract

Best
babakflame is offline   Reply With Quote

Old   February 19, 2015, 14:14
Default
  #453
Member
 
Join Date: Feb 2015
Posts: 31
Rep Power: 11
Stefano Puggelli is on a distinguished road
Hi,
thank you for the link and for the material. During this evening and in the next days I will have a look to the FGM model.
If I find somenthing interesting I will write here!
Thank you!
Stefano Puggelli is offline   Reply With Quote

Old   February 23, 2015, 05:52
Default
  #454
New Member
 
Sinasi Bahceci
Join Date: Jan 2015
Posts: 13
Rep Power: 11
sino75 is on a distinguished road
Hi Tobi,

I am trying to simulate a Bluff-body flame (Sandia) with the flameletModel (Holzmann) for my Master thesis. First, i wanted to do a mesh study. I proceeded as following:
- cold flow
- Enthalpy defect is 0
- high dissipation rate
When I want to create a pdf library for the solver, I get following error in my log-file:
Maximum number of iteration for calculating temperature from enthalpy...
So can you please tell me how to solve this problem?

Thanks
sino75 is offline   Reply With Quote

Old   February 23, 2015, 18:30
Default
  #455
Super Moderator
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Bad Wörishofen
Posts: 2,711
Blog Entries: 6
Rep Power: 52
Tobi has a spectacular aura aboutTobi has a spectacular aura aboutTobi has a spectacular aura about
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
Quote:
Originally Posted by sino75 View Post
Hi Tobi,

I am trying to simulate a Bluff-body flame (Sandia) with the flameletModel (Holzmann) for my Master thesis. First, i wanted to do a mesh study. I proceeded as following:
- cold flow
- Enthalpy defect is 0
- high dissipation rate
When I want to create a pdf library for the solver, I get following error in my log-file:
Maximum number of iteration for calculating temperature from enthalpy...
So can you please tell me how to solve this problem?

Thanks
You get the error with which solver?
Otherwise if it is from binary flamelet generator you make something wrong
__________________
Keep foaming,
Tobias Holzmann
Tobi is offline   Reply With Quote

Old   February 24, 2015, 06:00
Default
  #456
New Member
 
Sinasi Bahceci
Join Date: Jan 2015
Posts: 13
Rep Power: 11
sino75 is on a distinguished road
Hi Tobi,

there was a mistake with my dissipation, now I solved it.
Tobi likes this.
sino75 is offline   Reply With Quote

Old   March 18, 2015, 10:17
Default
  #457
Member
 
Join Date: Feb 2015
Posts: 31
Rep Power: 11
Stefano Puggelli is on a distinguished road
Hi Tobi,
did you ever try to test the libOpenSmoke code for premixed or partially premixed combustion? I know that the release is for non premixed combustion but maybe someone has modified it also to simulate premixed phenomena.
Stefano Puggelli is offline   Reply With Quote

Old   March 18, 2015, 12:02
Default
  #458
Super Moderator
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Bad Wörishofen
Posts: 2,711
Blog Entries: 6
Rep Power: 52
Tobi has a spectacular aura aboutTobi has a spectacular aura aboutTobi has a spectacular aura about
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
Dear Stefano,

I never coupled that tool with premixed combustion. I think therefore you also have to take into account ignition which is not implemented. Flamelet model in that kind means mixed equal burned.
__________________
Keep foaming,
Tobias Holzmann
Tobi is offline   Reply With Quote

Old   March 18, 2015, 12:23
Default
  #459
Member
 
Join Date: Feb 2015
Posts: 31
Rep Power: 11
Stefano Puggelli is on a distinguished road
hi Tobi,
thanks for your reply! OK, I will see if the libOpensmoke could be a good starting point also for premixed combustion.
Stefano
Stefano Puggelli is offline   Reply With Quote

Old   April 8, 2015, 13:56
Default Stabilized Flame
  #460
New Member
 
Sinasi Bahceci
Join Date: Jan 2015
Posts: 13
Rep Power: 11
sino75 is on a distinguished road
Hi,

I used flameletSimpleFoam at different grids. I'm really confused because I get a relatively stabilized flame with a grid of 6000 cells and 44000 cells, but when I use the same settings with a 16000cell-grid, no chance. I tried a lot of things like relaxation factors, maxCourant number or changed the schemes, but I had no success. Courant number doesn't seem to have influence at all. Residuals are relatively high.
Did you get such a problem and could you please help? I have attached link to my case:

https://www.dropbox.com/s/nxq2xzvmbf...se.tar.gz?dl=0

Thanks in regard.
Attached Images
File Type: jpg T_wrong.jpg (16.4 KB, 34 views)
File Type: jpg U_wrong.jpg (16.1 KB, 24 views)
sino75 is offline   Reply With Quote

Reply


Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Numerical treatment of the source term in combustion equations Tobi Main CFD Forum 37 September 15, 2020 14:42
[openSmoke] flameletSmoke + new ODESolver (by Alberto Cuoci) Tobi OpenFOAM Community Contributions 1 November 21, 2017 19:24
Unsteady solver with Flamelet Model (libOpenSMOKE) francesco_capuano OpenFOAM Running, Solving & CFD 11 November 26, 2013 05:50
LibOpenSmoke, getting the species in ParaFoam Christoph_84 OpenFOAM 1 May 31, 2012 15:42


All times are GMT -4. The time now is 16:54.