|
[Sponsors] |
PaSR + infinite reaction rate in reactingFoam --> no reactions occurring |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
November 15, 2012, 08:31 |
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 |
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 |
|
November 21, 2012, 05:02 |
|
#2 | |
New Member
Tatu Pinomaa
Join Date: Oct 2012
Location: Finland
Posts: 16
Rep Power: 13 |
Quote:
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 |
||
November 24, 2016, 18:34 |
|
#3 |
Member
Anonymouse
Join Date: Dec 2015
Posts: 98
Rep Power: 10 |
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. |
|
June 2, 2024, 10:04 |
|
#4 |
New Member
Join Date: Mar 2022
Posts: 7
Rep Power: 4 |
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. |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
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 |