CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Ignition modeling in OpenFOAM

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   September 5, 2011, 11:13
Question Ignition modeling in OpenFOAM
  #1
New Member
 
M. Timaji
Join Date: Sep 2010
Posts: 6
Rep Power: 15
timaji is on a distinguished road
Send a message via Skype™ to timaji
Dear Foamers,
What is the method in ignition process of solver XiFoam and engineFoam?
I found that the final results are so sensitive to ignition properties.
in OF ignition is characterized by "strength","duration","radius" in combustionProperties file, what are the true values for them for some ignitor? where can we find these values?

Timaji

Last edited by timaji; September 10, 2011 at 12:31.
timaji is offline   Reply With Quote

Old   September 11, 2011, 06:12
Default
  #2
New Member
 
Jacob Sardoei
Join Date: Sep 2011
Posts: 1
Rep Power: 0
Jacob is on a distinguished road
The method can be found in SRC/engine/ignition/ignite.H
here is the formulation:
bEqn.diag()[ignSite.cells()[icelli]] +=
(
ignSite.strength()*ignSite.cellVolumes()[icelli]
*rhou[ignSite.cells()[icelli]]/ignSite.duration()
) /(b[ignSite.cells()[icelli]] + 0.001);


Jacob
Jacob is offline   Reply With Quote

Old   September 11, 2011, 06:54
Default
  #3
New Member
 
M. Timaji
Join Date: Sep 2010
Posts: 6
Rep Power: 15
timaji is on a distinguished road
Send a message via Skype™ to timaji
Quote:
Originally Posted by Jacob View Post
The method can be found in SRC/engine/ignition/ignite.H
here is the formulation:
bEqn.diag()[ignSite.cells()[icelli]] +=
(
ignSite.strength()*ignSite.cellVolumes()[icelli]
*rhou[ignSite.cells()[icelli]]/ignSite.duration()
) /(b[ignSite.cells()[icelli]] + 0.001);


Jacob
Hi Jacob

I already saw this equation but the problem is that i dont know how to determine the ignition properties like strength,duration,radious?

according to "Solution adaptive CFD simulation of premixed flame propagation over various solid obstructions" by P.Naamansen, the ignition is carried out by changing the regress variable uniformly in space from one to zero in a sharply defined region.the region has the shape of semicircle in 2D simulation or semi-sphere in 3D simulation, and during the ramping which is a linear function of time no other combustion is allowed to occur.the time it takes to consume the reactants during that period is calculated from the laminar flame speed of gas.
as u see the radius and duration of ignition is proportional to laminar flame speed but there is nothing about the strength! I think the strength is probably the slope of ramp function.
but as i mentioned in previous post the final results are so sensitive to these quantities.

In the article there is no equation for relating the diameter and duration to laminar flame speed!

for example in XiFoam tutorial-moriyoshiHomogeneous- for propane with phi=1 and Su=0.434 how can we evaluate the diameter=0.003 and duration=0.003??


Regards

timaji
timaji is offline   Reply With Quote

Old   June 5, 2013, 17:18
Default
  #4
Member
 
Camille
Join Date: Oct 2012
Posts: 54
Rep Power: 13
camille131 is on a distinguished road
Hello Timaji,

Have you finally found the solution to your problem?
because I ask myself the same quesitons today ..

Could you please send me the article by P. Naamansen you mention in your previous post? (camilletienne75@gmail.com)

Thank you in advance,

Cam
camille131 is offline   Reply With Quote

Old   April 24, 2016, 22:33
Default
  #5
New Member
 
Pierluigi Olmati
Join Date: Apr 2016
Location: Tokyo
Posts: 21
Rep Power: 10
Pier84 is on a distinguished road
Hi All

It is a shame that there is not a minimum of documentation on that and no answer is given. Only bug reports are welcome for the convenience of the openFOAM(c) distributor.

PO
Pier84 is offline   Reply With Quote

Old   September 8, 2016, 06:57
Default
  #6
Senior Member
 
sheaker's Avatar
 
Oskar
Join Date: Nov 2015
Location: Poland
Posts: 184
Rep Power: 10
sheaker is on a distinguished road
Hello.
Because there is lack of information about ignition settings I would like to share with You my little compare that I've done during my studies.

Engine:
conRodLenght
0.150
bore
0.086
stroke
0.086
clearance
0.001
rmp
2000
deltaZ = 0.086
Piston speed = 5.73333 m/s
Vmax = 0.000547039, Vmin = 4.75896e-05
Vmax/Vmin = 11.4949

Plots and ignition settings:
See attachments.

I hope it helps a little for the beginners.
sheaker
Attached Images
File Type: png legend.png (2.9 KB, 116 views)
File Type: png plot.png (54.0 KB, 123 views)
File Type: png plot2.png (31.3 KB, 89 views)
sheaker is offline   Reply With Quote

Reply

Tags
enginefoam, ignition, openfoam, xifoam


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
Cross-compiling OpenFOAM 1.6 on Linux for Windows 32 and 64bits with Mingw-w64 wyldckat OpenFOAM Announcements from Other Sources 7 January 19, 2010 15:39
OpenFOAM Install problem masb OpenFOAM 3 May 25, 2009 11:32
The OpenFOAM extensions project mbeaudoin OpenFOAM 16 October 9, 2007 09:33
OpenFOAM Training and Workshop Hrvoje Jasak Main CFD Forum 0 October 7, 2005 07:14
Ignition modeling Julie Polyakh Siemens 0 May 20, 2004 08:05


All times are GMT -4. The time now is 00:24.