CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Programming & Development

New solver: ddtFoam

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

Like Tree8Likes

Reply
 
LinkBack Thread Tools Display Modes
Old   August 25, 2013, 07:30
Default New solver: ddtFoam
  #1
Member
 
Florian Ettner
Join Date: Mar 2009
Location: Munich, Germany
Posts: 40
Rep Power: 8
dohnie is on a distinguished road
Dear Foamers,
I developed an OpenFOAM solver for the simulation of the deflagration-to-detonation transition in gases.
Some features:
  • accurate shock capturing with a density based solver (largely based on Oliver Borm's development densityBasedTurbo)
  • deflagration modelling via a reaction progress variable (based on XiFoam)
  • auto-ignition sub-grid model separating each computational cell into a "shocked" and an "unshocked" part
  • auto-ignition from tabulated ignition delay times, gained from a detailed chemical mechanism
  • low computational cost due to applicability to coarse grids and usage of pre-tabulated chemistry
Even if you're not interested in detonations/combustion, it might be a nice alternative to existing OpenFOAM solvers when it comes to shock resolution in transonic flows.

The code incl. a short documentation and tutorials can be downloaded from
http://sourceforge.net/projects/ddtfoam/
Attached Images
File Type: png HLLC_limited_small.png (64.0 KB, 160 views)
File Type: jpg DDT_mit_Legende_small.jpg (72.8 KB, 182 views)
dohnie is offline   Reply With Quote

Old   August 27, 2013, 13:56
Default
  #2
Member
 
Chris
Join Date: Aug 2012
Location: Calgary, Alberta, Canada
Posts: 73
Rep Power: 5
ChrisA is on a distinguished road
Great work, I haven't had much time to look at it but from first glances the solver looks pretty good. The guys at my lab working on deflagration to detonation will certainly find it interesting. Is there going to be a paper to accompany the solver? I'd be very interested in reading that and seeing some results from the solver.
ChrisA is offline   Reply With Quote

Old   August 27, 2013, 14:11
Default
  #3
Member
 
Florian Ettner
Join Date: Mar 2009
Location: Munich, Germany
Posts: 40
Rep Power: 8
dohnie is on a distinguished road
Chris,
I have a paper in preparation, but at the moment I can only refer to my thesis (in German, unfortunately):
Effiziente numerische Simulation des Deflagrations-Detonations-Übergangs.pdf
[8.6 MB]
Tobi likes this.
dohnie is offline   Reply With Quote

Old   August 27, 2013, 14:14
Default
  #4
Member
 
Chris
Join Date: Aug 2012
Location: Calgary, Alberta, Canada
Posts: 73
Rep Power: 5
ChrisA is on a distinguished road
Ah, yes, I found the thesis shortly after posting that, unfortunately my German is terrible. However, I do look forward to the other paper you mentioned
ChrisA is offline   Reply With Quote

Old   October 30, 2013, 06:46
Default
  #5
New Member
 
Knight Huang
Join Date: Oct 2012
Posts: 9
Rep Power: 4
uwknight is on a distinguished road
Hi dohnie,

Thanks for your sharing this good solver. I am working on this ddt solvers.

Do you HAVE some ideads about how to change the reaction model? Can I use this code to simulate the problem about oblique detonation wave? Could you give me some comments? You comments should be very useful to me. Thanks a lot.

Btw, shall I have your email? You can reply me using priviate message. Thanks.

Best regards,

Yue
uwknight is offline   Reply With Quote

Old   November 11, 2013, 15:38
Default
  #6
Member
 
Florian Ettner
Join Date: Mar 2009
Location: Munich, Germany
Posts: 40
Rep Power: 8
dohnie is on a distinguished road
Yue,
sorry for the late reply, I don't work on this topic any longer.
You can use this solver to simulate oblique detonations as well, and of course you can also modify the reaction model.
From your short description I cannot see what kind of problem you really want to simulate, but it seems to me that you are not interested in simulating the deflagration-to-detonation transition. Thus, if you want to investigate only detonations (probably even looking for steady solutions?), you might be better off using simpler solvers. For example you could use a compressible solver (like rhoCentralFoam) and add 1-step Arrhenius kinetics - the publications of Elaine Oran et al. might give you some orientation.

I hope that helps!
Florian
dohnie is offline   Reply With Quote

Old   November 21, 2013, 04:39
Default
  #7
New Member
 
Knight Huang
Join Date: Oct 2012
Posts: 9
Rep Power: 4
uwknight is on a distinguished road
[QUOTE=dohnie;461587]Yue,
sorry for the late reply, I don't work on this topic any longer.
You can use this solver to simulate oblique detonations as well, and of course you can also modify the reaction model.
From your short description I cannot see what kind of problem you really want to simulate, but it seems to me that you are not interested in simulating the deflagration-to-detonation transition. Thus, if you want to investigate only detonations (probably even looking for steady solutions?), you might be better off using simpler solvers. For example you could use a compressible solver (like rhoCentralFoam) and add 1-step Arrhenius kinetics - the publications of Elaine Oran et al. might give you some orientation.

I hope that helps!
Florian[/QUOTE

Hi Florian,

Thanks for your reply. I would do some works on your mentioned.

Yue
uwknight is offline   Reply With Quote

Old   November 28, 2013, 03:47
Default
  #8
New Member
 
Zhong Wei
Join Date: Nov 2013
Posts: 1
Rep Power: 0
lengshui222 is on a distinguished road
[QUOTE=uwknight;462854]
Quote:
Originally Posted by dohnie View Post
Yue,
sorry for the late reply, I don't work on this topic any longer.
You can use this solver to simulate oblique detonations as well, and of course you can also modify the reaction model.
From your short description I cannot see what kind of problem you really want to simulate, but it seems to me that you are not interested in simulating the deflagration-to-detonation transition. Thus, if you want to investigate only detonations (probably even looking for steady solutions?), you might be better off using simpler solvers. For example you could use a compressible solver (like rhoCentralFoam) and add 1-step Arrhenius kinetics - the publications of Elaine Oran et al. might give you some orientation.

I hope that helps!
Florian[/QUOTE

Hi Florian,

Thanks for your reply. I would do some works on your mentioned.

Yue
Hi Yue,

I want to work on this ddt solvers, but when I am installing it, it shows some mistake like "expected unqualified-id before '.' token". Have you met the similar mistake when you were installing the ddt solver? Can you show me the detailed steps of installing this solver, if you remember the installing steps? Thank you very much.

Best regards,

Zhong
lengshui222 is offline   Reply With Quote

Old   February 4, 2014, 07:00
Default Hi Dohnie
  #9
New Member
 
pratap
Join Date: Nov 2013
Posts: 21
Rep Power: 3
cfdsat is on a distinguished road
You have used two different solver one based on densityBasedTurbo solver package and other based on sonicFoam. What is your experience regarding which one is best for solving compressible and DDT and detonation problems.

I have got your PhD thesis, since it is written in German, I am not able to understand it. Could you please provide me a link to your papers if published ?
cfdsat is offline   Reply With Quote

Old   February 11, 2014, 16:18
Default
  #10
Member
 
Florian Ettner
Join Date: Mar 2009
Location: Munich, Germany
Posts: 40
Rep Power: 8
dohnie is on a distinguished road
The density based solver is definitely better for simulating DDT and detonation.
I use the pressure based solver only to get an initial solution for unsteady problems where I start in quiescent or very low-speed (Mach number < 0.1) flow, because the density based solver tends to get unstable at very low Mach numbers.

I have prepared a journal paper which is currently in the review process. I'll let you know once it has been accepted.

Regards,
Florian
dohnie is offline   Reply With Quote

Old   February 12, 2014, 06:03
Default
  #11
New Member
 
pratap
Join Date: Nov 2013
Posts: 21
Rep Power: 3
cfdsat is on a distinguished road
Hi Dohnie

Thanks for the comments. Which openfoam version did you run this solver ?

Can you already send me the pdf file of the publication you have submitted ?
cfdsat is offline   Reply With Quote

Old   June 2, 2014, 05:56
Post Paper published
  #12
Member
 
Florian Ettner
Join Date: Mar 2009
Location: Munich, Germany
Posts: 40
Rep Power: 8
dohnie is on a distinguished road
Dear all,
a paper describing the features of ddtFoam has now been published:

http://www.hindawi.com/journals/jc/2014/686347/

Florian Ettner, Klaus G. Vollmer, and Thomas Sattelmayer, “Numerical Simulation of the Deflagration-to-Detonation Transition in Inhomogeneous Mixtures,” Journal of Combustion, vol. 2014, Article ID 686347, 14 pages, 2014. doi:10.1155/2014/686347
dohnie is offline   Reply With Quote

Old   June 5, 2014, 10:08
Default
  #13
New Member
 
Anh
Join Date: Jan 2012
Posts: 1
Rep Power: 0
anh065 is on a distinguished road
Great. Thanks, Florian, for sharing this. I would want to apply this code to problems with non-confined domains and wonder what types of boundary condition should be applied to the open (non-reflective) boundaries?
anh065 is offline   Reply With Quote

Old   June 7, 2014, 15:42
Default
  #14
Member
 
Florian Ettner
Join Date: Mar 2009
Location: Munich, Germany
Posts: 40
Rep Power: 8
dohnie is on a distinguished road
Dear Anh,
you should get the best result by using the waveTransmissive BC.
However, this will require some extra programming work. The standard wave transmissive BC in OpenFOAM is based on the assumption of a constant ratio of specific heats (gamma) and a flux field which is called phi in most OpenFOAM solvers. In ddtFoam, where the fluxes are computed using a Riemann solver, this flux corresponds to Riemann.rhoUFlux() divided by Riemann.rhoFlux().
I hope that helps!
Florian
dohnie is offline   Reply With Quote

Old   June 11, 2014, 03:21
Smile
  #15
New Member
 
Yoshiaki SENDA
Join Date: Jun 2014
Location: Kyoto, Japan
Posts: 2
Rep Power: 0
panda1100 is on a distinguished road
Dear dohnie,

Thank you for sharing this great work

Is there any information about how to change H2+air to Diesel mixture?
And If its possible, Can you share your utility codes?
for example,
  • loolupIgn_06.m
  • automaticSetFields_tfc_zu_Grad3s_03.c etc
I want to try this solver with another type of gas.
(I tried your tutorial and that works pretty good.)

Any information you can provide me would be greatly appreciated.
I wish I could read German...
panda1100
panda1100 is offline   Reply With Quote

Old   June 18, 2014, 04:05
Default
  #16
Member
 
Florian Ettner
Join Date: Mar 2009
Location: Munich, Germany
Posts: 40
Rep Power: 8
dohnie is on a distinguished road
Dear Yoshiaki,
in principle you just have to replace the tabulated ignition delay times in the respective case. However, you should be aware that the heat release of long-chained hydrocarbons is not as instantaneous as the heat release of hydrogen.
See, for example:
O. Colin, A. Pires da Cruz and S. Jay: "Detailed Chemistry-Based Auto-Ignition Model Including Low Temperature Phenomena Applied to 3-D Engine Calculations", Proceedings of the Combustion Institute, Vol.30, pp.2649-2656 (2005)

Send me your e-mail address in a PM and I can send you the files you requested.

Regards,
Florian
dohnie is offline   Reply With Quote

Old   June 19, 2014, 00:12
Default
  #17
New Member
 
Yoshiaki SENDA
Join Date: Jun 2014
Location: Kyoto, Japan
Posts: 2
Rep Power: 0
panda1100 is on a distinguished road
Dear Florian,

Thank you so much for your advice and your help.
I sent you PM with my e-mail address .

Thank you again for sparing your time for me.

Best regards,
Yoshiaki
panda1100 is offline   Reply With Quote

Old   June 19, 2014, 08:41
Default Great!
  #18
New Member
 
Bruno Blais
Join Date: Sep 2013
Location: Canada
Posts: 24
Rep Power: 3
blais.bruno is on a distinguished road
Looks both very good and very interesting!

Thanks for posting it online this way...
blais.bruno is offline   Reply With Quote

Old   September 1, 2014, 11:26
Default
  #19
New Member
 
Join Date: May 2011
Posts: 28
Rep Power: 6
dowlee is on a distinguished road
Dear Florian Ettner
Thanks for your efforts on developing solver DDTFoam. I am very interested to this solver. And at present, I am working on a project concerning detonation combustion. In my case, the combustible mixture close to the end wall of the resonators is heated up to a temperature that creates ignition automatically. And after a short while, the flame front reaches from behind the reflected shock wave and combines with it forming a plannar detonation wave. My question is whether this solver is suitable for automatic ignition and detonation. If not, what should be modified in this solver. I am looking forward to your reply.
Thanks very much
Dowlee
dowlee is offline   Reply With Quote

Old   September 8, 2014, 12:57
Default
  #20
Member
 
Florian Ettner
Join Date: Mar 2009
Location: Munich, Germany
Posts: 40
Rep Power: 8
dohnie is on a distinguished road
Dear Dowlee,
I don't understand from your description HOW the mixture close to the end wall is heated up, but I assume it is done by a shock, as you mention a reflected shock wave later on.
I think the solver should work fine for your problem!
dohnie is offline   Reply With Quote

Reply

Tags
combustion, compressible, density based, detonation, transonic

Thread Tools
Display Modes

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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Fluent Adjoint Solver? ex10148 FLUENT 15 July 21, 2015 01:13
thobois class engineTopoChangerMesh error Peter_600 OpenFOAM 4 August 2, 2014 09:52
Divergence problem Smaras FLUENT 13 February 21, 2013 06:03
3d vof Smaras FLUENT 2 February 19, 2013 07:58
why the solver reject it? Anyone with experience? bearcat CFX 6 April 28, 2008 14:08


All times are GMT -4. The time now is 02:10.