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

having trouble using reactingFoam with reactions turned off

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

Like Tree4Likes

Reply
 
LinkBack Thread Tools Display Modes
Old   October 25, 2010, 02:13
Post having trouble using reactingFoam with reactions turned off
  #1
New Member
 
Join Date: Aug 2010
Location: Mumbai
Posts: 4
Rep Power: 6
kaka22 is on a distinguished road
Hi,

I have been simulating the counterFlame2D tutorial file using reactingFoam in version 1.7.1. As I wanted to just solve for the species transport without reaction, I have switched reactions off by specifying in constant/chemistryProperties reactions off and turbulent reactions off. Temperatures at 2 inlets are 573K and 373K. I should expect that temperature will not increase beyond 573K since there are no reactions and I am just solving for the flow. However, the temperature increased beyond 573 K. May I know what exactly that I have done wrong? I have been trying to solve my problem for a long time and tried to think of all means and change any necessary parameters but this problem still persist. Would appreciate any help. Thank you!
kaka22 is offline   Reply With Quote

Old   November 8, 2010, 12:24
Default Same strange behaviour
  #2
New Member
 
Valerio Novaresio
Join Date: Mar 2009
Location: Polonghera, Cuneo, Italy
Posts: 27
Rep Power: 8
novyno is on a distinguished road
Hi all,

I also have the same problem (reactingFoam + no reactions). In a domain with two inlets (both at 1000 K) I start the simulation setting the internalField for T = 1000. During the simulation the temperature range varies between 985 and 1034 K.

I can't explain me this strange behaviour.

Valerio
__________________
...The best way to acquire new knowledge is to share it...
novyno is offline   Reply With Quote

Old   November 9, 2010, 02:01
Default
  #3
Senior Member
 
Nakul
Join Date: Apr 2010
Location: India
Posts: 147
Rep Power: 7
nakul is on a distinguished road
Hi kaka22,

How much is your maximum temperature?

If its only a little more, it might be due to sudden compression of gas in some regions of your domain.

If its considerably more then I don't have a clue. There might be some problem with your case setup.
nakul is offline   Reply With Quote

Old   November 10, 2010, 02:15
Default
  #4
New Member
 
Join Date: Aug 2010
Location: Mumbai
Posts: 4
Rep Power: 6
kaka22 is on a distinguished road
Hi,

The problem I have mentioned earlier have been reported and I got a reply.

"reactingFoam algorithm is an efficient, approximate method for turbulent reacting flow. However, in the laminar limit the approximation breaks down, leading to temperature inconsistencies.

In the laminar case, if you have a composition gradient, the only way to ensure that the diffusion does not give rise to temperature errors is by enforcing a Lewis no = 1 constraint, e.g. update the sensible enthalpy eqn to:

fvScalarMatrix hsEqn
(
fvm::ddt(rho, hs)
+ mvConvection->fvmDiv(phi, hs)
// - fvm::laplacian(turbulence->alphaEff(), hs)
- fvm::laplacian(turbulence->muEff(), hs) // unit lewis no.
==
DpDt
+ chemistrySh
);
"
So basically, I have commented this line // - fvm::laplacian(turbulence->alphaEff(), hs) and run the case again but with inlet velocities as 0.1 m/s, reactions off and turbulent reaction off, initial temperature as 573 K, inlet temperature at 373 K and 573 K. The temperature remains approximately constant at 573K, basically the temperature does not shoot up high as expected since there's no reaction.

However, I have run a similar second case with varying only the inlet velocities. I have changed inlet velocities to be 100 m/s and 289 m/s, still with the same temperature and reactions off. The temperature gradually increases high during the simulation.

I guess the solver works well only for laminar case? since at low velocity we are not facing the problem. Probably we still have to look further into the solver! Hope this helps, cheers
harshad88 likes this.
kaka22 is offline   Reply With Quote

Old   November 10, 2010, 07:17
Default
  #5
Senior Member
 
Nakul
Join Date: Apr 2010
Location: India
Posts: 147
Rep Power: 7
nakul is on a distinguished road
Hi,
Do you mean that at higher velocities when you did the simulation with unit Lewis No., you again got temperature inconsistencies?

Are you getting agreeable results with combustion(reaction) on?
nakul is offline   Reply With Quote

Old   November 10, 2010, 23:13
Default
  #6
New Member
 
Join Date: Aug 2010
Location: Mumbai
Posts: 4
Rep Power: 6
kaka22 is on a distinguished road
Hi Nakul,

Yes that's right I got temperature inconsistencies still when I simulated with unit Lewis No. for flow at high velocities.

For my case when I simulated it with reactions on, the temperature profile came out fine but I am not getting mass fraction of product as desired. So I can't really tell if I get agreeable results with reactions on since there are some things that are still inconsistent.

Have you run any case with reactingFoam and got the result as per expected?
kaka22 is offline   Reply With Quote

Old   November 11, 2010, 01:30
Default
  #7
Senior Member
 
Nakul
Join Date: Apr 2010
Location: India
Posts: 147
Rep Power: 7
nakul is on a distinguished road
Hi,

I am also trying to understand reactingFOAM and I am in process of making my test case and I am also facing some problems and that's why I thought it might be useful to discuss it with someone whose also using it.

From what I have heard reactingFOAM is a good solver. There is some publication also out there which uses reactingFOAM in ramjet engine combustor.

Have you taken a look at the code in Y-equ.H ? Its also a little bit different from what I expected it to be. There also I think variation of Lewis No. is playing a role (in terms of "mueff()").

According to me the species equation should be something like the one in the attachment. I am still confused about its present formulation.

Do you have any ideas about this?
Attached Images
File Type: png Screenshot.png (6.8 KB, 246 views)
nakul is offline   Reply With Quote

Old   November 12, 2010, 06:55
Default Thanks for your help and something about y.eqn
  #8
New Member
 
Valerio Novaresio
Join Date: Mar 2009
Location: Polonghera, Cuneo, Italy
Posts: 27
Rep Power: 8
novyno is on a distinguished road
Hi all!

First of all thanks to kaka22 for his exaustive reply! It was very useful for me!

Then speaking about y.eqn its form is different than the usual one becouse also there in my opinion there is the hypotesis of Sc = 1 for all species (so we can sobstitute the rho*D term with mu).

This is done (always in my opinion) becouse OpenFOAM doesn't (officially) have a complete library for multispecies diffusion model (multispecies Fick model, Maxwell-Stefan model, etc...). I'm working on it (multispecies diffusion model on OpenFOAM ) and I hope to release soon this library. Actually I'm trying to combine my library of diffusion models with reactingFoam in order to leave the rho*D term in y.eqn and h.eqn without make any hypotesis on Le or Sc numbers.

As soon as I will obtain some relevant results I'll update you!

Regards,

Valerio
Thamali likes this.
__________________
...The best way to acquire new knowledge is to share it...

Last edited by novyno; November 15, 2010 at 04:19.
novyno is offline   Reply With Quote

Old   November 12, 2010, 07:36
Default
  #9
Senior Member
 
Nakul
Join Date: Apr 2010
Location: India
Posts: 147
Rep Power: 7
nakul is on a distinguished road
Hi Valerio,

Correct me if I am wrong but I think that if assumption of Le=1 is there then we should have alpha instead of mu in the species equation.

In hs-equ. also for Le=1 they have used alphaEff().

According to me muEff() is there to account for turbulent diffusion. Please clarify me if I am wrong!!!
nakul is offline   Reply With Quote

Old   November 25, 2010, 01:17
Default
  #10
Senior Member
 
Nakul
Join Date: Apr 2010
Location: India
Posts: 147
Rep Power: 7
nakul is on a distinguished road
Hi,

I am facing another problem with rectingFOAM. Hope that somebody here might be able to help!!!

I am trying to run a case using reactingFOAM. But I am getting the following error:

FOAM FATAL ERROR : attempt to use janafThermo<equationofstate> out of temperature range 200 -> 5000; T = 195.559.

From function janafThermo<equationofstate>::checkT(const scalar T) const
in file /home/dm2/henry/OpenFOAM/OpenFOAM-1.2/src/thermophysicalModels/specie/lnInclude/ janafThermoI.H at line 73

The temperatures that I have supplied as BC are 1000K and 600K for O2 and H2 respectively. (Its H2-O2 combustion.)

My Courant No. is 0.2 and my max. cell skewness = 0.66838.
My "checkMesh" results are all OK.

Can anybody please tell me where am I going wrong?

You may look at the posts in this discussion also :

Problem With reactingFOAM

-Nakul
nakul is offline   Reply With Quote

Old   November 25, 2010, 06:55
Default Which version?
  #11
New Member
 
Valerio Novaresio
Join Date: Mar 2009
Location: Polonghera, Cuneo, Italy
Posts: 27
Rep Power: 8
novyno is on a distinguished road
Hi, nakul

Do you use the standard version of reactingFoam (1.7.* version without any modifications in the hs equation)?
__________________
...The best way to acquire new knowledge is to share it...
novyno is offline   Reply With Quote

Old   November 26, 2010, 01:41
Default
  #12
Senior Member
 
Nakul
Join Date: Apr 2010
Location: India
Posts: 147
Rep Power: 7
nakul is on a distinguished road
No I don't use the standard version of reactingFOAM.

I have modified it for supersonic combustion and hence the energy equation has been modified to solve for supersonic flows.

Actually I have modified the energy equation of rhoCentralFOAM to include turbulence and reaction. However the energy equation still solves for hs only.
-Nakul
nakul is offline   Reply With Quote

Old   November 26, 2010, 05:07
Default
  #13
New Member
 
Valerio Novaresio
Join Date: Mar 2009
Location: Polonghera, Cuneo, Italy
Posts: 27
Rep Power: 8
novyno is on a distinguished road
Can you post your energy equation?
__________________
...The best way to acquire new knowledge is to share it...
novyno is offline   Reply With Quote

Old   November 26, 2010, 05:11
Default
  #14
Senior Member
 
Nakul
Join Date: Apr 2010
Location: India
Posts: 147
Rep Power: 7
nakul is on a distinguished road
Hi
I am attaching my energy equation. If you have any doubts, feel free to ask right now.

-Nakul
nakul is offline   Reply With Quote

Old   November 26, 2010, 05:13
Default
  #15
Senior Member
 
Nakul
Join Date: Apr 2010
Location: India
Posts: 147
Rep Power: 7
nakul is on a distinguished road
Sorry the energy equation is attached here.
Attached Files
File Type: txt energy equation.txt (1.8 KB, 231 views)
nakul is offline   Reply With Quote

Old   November 26, 2010, 05:16
Default
  #16
Senior Member
 
Nakul
Join Date: Apr 2010
Location: India
Posts: 147
Rep Power: 7
nakul is on a distinguished road
My solver is running for around 500 time steps beginning with 0.

It goes to a time of around 2e-6 and then I get the above error that T has exceeded 5000 K limit to become ~5200K.

I have also tried different BC like waveTransmissive, inletOutlet at outlet for p and U respectively and also I have tried zeroGradient.

Any help would be appreciated.
nakul is offline   Reply With Quote

Old   November 26, 2010, 05:38
Default
  #17
New Member
 
Valerio Novaresio
Join Date: Mar 2009
Location: Polonghera, Cuneo, Italy
Posts: 27
Rep Power: 8
novyno is on a distinguished road
Only one question more.

If you switch off all the reactions and if you set all the inlets and the internal field to the same temperature, what do you obtain?

In line of principle you have to obtain good results (uniform temperature or temperature variations only due to the compressibility). Probably if you have in your domain an uniform mixture you will obtain this. But if you have some gradients of the species you will probably obtain some strange results.

I am not an expert in supersonic flows, but looking your energy equations it seems that miss some term related to the gradient of the composition. And this, in my opinion, is the cause of your strange behavior of the temperature (to hight if I have well understood).

Try this test case and update me.
__________________
...The best way to acquire new knowledge is to share it...
novyno is offline   Reply With Quote

Old   May 26, 2011, 04:03
Default non-unit lewis number
  #18
Member
 
Yashar Afarin
Join Date: May 2010
Location: Toronto- Canada
Posts: 40
Rep Power: 7
yashar.afarin is on a distinguished road
Send a message via Skype™ to yashar.afarin
Quote:
Originally Posted by novyno View Post
Hi all!

First of all thanks to kaka22 for his exaustive reply! It was very useful for me!

Then speaking about y.eqn its form is different than the usual one becouse also there in my opinion there is the hypotesis of Sc = 1 for all species (so we can sobstitute the rho*D term with mu).

This is done (always in my opinion) becouse OpenFOAM doesn't (officially) have a complete library for multispecies diffusion model (multispecies Fick model, Maxwell-Stefan model, etc...). I'm working on it (multispecies diffusion model on OpenFOAM ) and I hope to release soon this library. Actually I'm trying to combine my library of diffusion models with reactingFoam in order to leave the rho*D term in y.eqn and h.eqn without make any hypotesis on Le or Sc numbers.

As soon as I will obtain some relevant results I'll update you!

Regards,

Valerio
Hi valerio

I am using reactingFoam and I am very interesting to using multi-species diffusion model and compare its results with my result when lewis was considered as one. have you had any progress to implementing multi-species diffusion model in reactingFoam? I would appreciate it if you could kindly help me.

best regards
yashar
nichkoohi likes this.
yashar.afarin is offline   Reply With Quote

Old   June 8, 2011, 05:08
Default
  #19
Senior Member
 
Illya Shevchuk
Join Date: Aug 2009
Location: Darmstadt, Germany
Posts: 174
Rep Power: 7
linch is on a distinguished road
One question:

I don't get how the term
Code:
fvm::laplacian(alphaEff(),hs)
in the enthalpie equation can work in multicomponent case? Different components have different standard enthalpies of formation, so the term exists even without temperature gradient. pre.cjk { font-family: "DejaVu Sans",monospace; }p { margin-bottom: 0.21cm; }
Code:
fvc::laplacian(k, T)
Seems to make sense, but It can't be treated implicitly, so internal iterations with temperature update needed.

What is the error in my reasoning?

Last edited by linch; June 16, 2011 at 07:12.
linch is offline   Reply With Quote

Old   June 10, 2011, 07:24
Default
  #20
Senior Member
 
Illya Shevchuk
Join Date: Aug 2009
Location: Darmstadt, Germany
Posts: 174
Rep Power: 7
linch is on a distinguished road
I'm surprised, because I thought it must be something obvious and everyone who uses reactingFoam / rhoReactingFoam can help me.

Best,
linch is offline   Reply With Quote

Reply

Tags
hseqn, reactingfoam

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
ReactingFoam with surface reactions robert_mornhinweg OpenFOAM Running, Solving & CFD 8 April 21, 2015 06:20
No reactions in reactingFoam 2.1 OMN OpenFOAM Running, Solving & CFD 16 April 7, 2015 12:14
PaSR + infinite reaction rate in reactingFoam --> no reactions occurring tatu OpenFOAM Running, Solving & CFD 1 November 21, 2012 06:02
reactingFoam wedge handling wrong U dhondupant OpenFOAM Bugs 1 December 9, 2010 08:34
ReactingFoam without reactions lasb OpenFOAM Running, Solving & CFD 5 June 10, 2008 08:50


All times are GMT -4. The time now is 18:09.