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

PaSR + infinite reaction rate in reactingFoam --> no reactions occurring

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

Like Tree1Likes
  • 1 Post By KarenRei

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   November 15, 2012, 08:31
Default PaSR + infinite reaction rate in reactingFoam --> no reactions occurring
  #1
New Member
 
Tatu Pinomaa
Join Date: Oct 2012
Location: Finland
Posts: 16
Rep Power: 13
tatu is on a distinguished road
Hello!

Background:
First of all I want to say I'm a noobie with OpenFOAM, and I'm using OF 2.1.1. My goal is to model hydrogen combustion in a high temperature as a one-step reaction, so that the combustion reaction will not be a pure 2H2 + O2 => 2H202, but will include dissociated molecules (e.g. O, OH). Therefore Arrhenius data, i.e. reaction rate coefficient A or activation energy Ea, will not be available.
Also, the flow is compressible, so it is important to make the chemistry model as computationally light as possible. Therefore, I will be using (at least for now) eddy-dissipation model with infinite chemical reaction rates, i.e. I assume that turbulent mixing is the bottleneck for the combustion.

Problem:
I am studying the case counterFlowFlame2D under reactingFoam solver, where the reaction rate is controlled by both the turbulent mixing and Arrhenius chemical kinetics (i.e. PaSR model). Now my goal is to simulate the combustion in this test case simply with the eddy-dissipation model, so that I switch off Arrhenius chemistry.

Q1: To use the simpler eddy-dissipation model without Arrhenius chemistry, should I merely change the reaction type from irreversibleArrheniusReaction to irreversibleinfiniteReaction in file constant/reactions? Are there other options I should consider for this case?

Q2: Why are there practically no reactions occurring when I use irreversibleinfiniteReaction? Shouldn't the infinite chemical reaction rate model increase the reaction rate, and not decrease it?

As you can see from the two figures (infinite reaction rate and Arrhenius chemistry) for time 0.3, the number of reaction products are practically non-existing for the infinite reaction rate. Also, for the case of infinite reaction rate, the maximum temperature doesn't exceed 800K, which is the inlet temperature of the oxygen and CH4.

Thanks in advance,

Tatu
Attached Images
File Type: png arrhenius.png (16.0 KB, 69 views)
File Type: png infinite.png (23.2 KB, 62 views)
tatu is offline   Reply With Quote

Old   November 21, 2012, 05:02
Default
  #2
New Member
 
Tatu Pinomaa
Join Date: Oct 2012
Location: Finland
Posts: 16
Rep Power: 13
tatu is on a distinguished road
Quote:
Originally Posted by tatu View Post

Q1: To use the simpler eddy-dissipation model without Arrhenius chemistry, should I merely change the reaction type from irreversibleArrheniusReaction to irreversibleinfiniteReaction in file constant/reactions? Are there other options I should consider for this case?

Tatu
I tried to study the reactingFoam solver, and here's where I am now.

I figured out that by setting turbulentReaction to false, the reactive volume fraction (kappa in the PaSR model), that scales the fuel consumption rate RR(specieI) for each species i, will be automatically set to 1. By doing this, the reaction rate increases significantly, as it should.

I cant find any other dependencies on Arrhenius data in PaSR.C, than on lines 82-89 (useReactionRate == true), and lines 95-132 (turbulentReaction == true).
However, setting useReactionRate and turbulentReaction as false in constant/combustionProperties and then fiddling with the Arrhenius data, affects the reaction rate (increasing the activation temperature Ta in file constant/reactions decreases the reaction rate significantly).

Am I right that all the chemistry in the reactingFoam solver can be seen from PaSR.C?

Where exactly does the Arrhenius data come into play in the chemistry of reactingFoam solver (PaSR)?

And still, by setting reaction type to irreversibleinfiniteReaction from infiniteArrheniusReaction, the reaction rate decreases to almost zero. Why doesn't the reaction rate increase instead?

Cheers,

Tatu
tatu is offline   Reply With Quote

Old   November 24, 2016, 18:34
Default
  #3
Member
 
Anonymouse
Join Date: Dec 2015
Posts: 98
Rep Power: 10
KarenRei is on a distinguished road
This thread is years old, but I've spent about 4 days scratching my head on this one and I can confirm it: "irreversibleinfiniteReaction" reaction is actually incredibly slow for some reason.

Don't use it.
lumpor likes this.
KarenRei is offline   Reply With Quote

Old   June 2, 2024, 10:04
Default
  #4
New Member
 
Join Date: Mar 2022
Posts: 7
Rep Power: 4
ngbidi is on a distinguished road
This thread is old, but apparently the problem is still there as I am experiencing this in OF-v2012. I think it is a bug in infiniteReactionRateI.H.
Have you found a solution to this?

Last edited by ngbidi; June 2, 2024 at 11:12.
ngbidi is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

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


Similar Threads
Thread Thread Starter Forum Replies Last Post
Implement code to calculate the reaction rate in reactingFoam matari OpenFOAM Programming & Development 17 November 16, 2018 09:43
Segmentation fault in running alternateSteadyReactingFoam,why? NewKid OpenFOAM 18 January 20, 2011 16:55
UDF for catalyst volumetric reaction, temperature overshoot though rate is small elmcmaster FLUENT 0 October 11, 2010 09:37
UDF changing the rate exponent of a reaction Stefan H Fluent UDF and Scheme Programming 0 September 16, 2009 13:20
chemical reaction - decompostition La S. Hyuck CFX 1 May 23, 2001 00:07


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