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

species mass fraction in reacting flow

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

Reply
 
LinkBack Thread Tools Display Modes
Old   June 16, 2004, 17:01
Default species mass fraction in reacting flow
  #1
ma
Guest
 
Posts: n/a
In reacting flow simulation, the species mass fraction may become less than 0 or larger than 1 (i.e., Z<0, Z>1), which is unphysical.

What's the general way to overcome it?

Thanks.

-ma
  Reply With Quote

Old   June 17, 2004, 02:25
Default Re: species mass fraction in reacting flow
  #2
Rami
Guest
 
Posts: n/a
For negative mass fraction - correct by resetting to 0. For excessively high values, actually the condition is even more strict than you suggest: the sum of mass fractions should be 1. A corrective measure I use is to re-normalize the mass fractions of all species such that the sum is 1. These corrections should be applied within the iterative solution.

I hope this helps.
  Reply With Quote

Old   June 17, 2004, 04:17
Default Re: species mass fraction in reacting flow
  #3
jdc
Guest
 
Posts: n/a
Hi Ma,

It is also usefull to identify why you are observing this unphysical behavior. It could be a usefull exercise to check whether there is any coding problem. the known problems that results in this kind of behavior are usually:

1) unbounded numerical schemes;

2) diffusion fluxes that can not lead to \sum Yk = 1 (for example, if you are using unity Lewis number).

3) implicit segragated solver.

The solution are (1) to use more appropriate scheme/approximation/solver or (2) to do Rami's tricks.

Note: there is an additional trick that is as follow (a) solve for NS-1 species, (b) force all solved mass fractions to be between 0. and 1., (c) compute the unsolved species as Y_{NS} = 1. - \sum_{NS-1} Yk. This trick works resonnably if species NS is not reacting (such as helium, nitrogen....)

Hoppe this help. Julien

  Reply With Quote

Old   June 17, 2004, 14:53
Default Re: species mass fraction in reacting flow
  #4
ma
Guest
 
Posts: n/a
Hi, Rami and Julien,

Thanks a lot for your helpful suggestions.

Since my code/solver has already been set up, I'll have to use the post-remedy as you and Rami mentioned. I later checked Oran and Poris's book "numerical simulation of reactive flow", which also referred to this kind of remedy.

A consequent question is that: When we post-remedy the mass fractions, shall we also modify the conserved varable (Q4=rho*et, in 2D)?

I tested my reacting Euler solver for a non-reacting flow problem (a species contact surface moving from a tube to a unconfined region), and found that: (1) without post-remedy of mass fraction, the temperature field is correct but the mass fraction may reach 1.2 somewhere even with small CFL. (2) with post-remedy of mass fraction, the mass fraction field becomes correct but the temperature reaches unphysically high somewhere. (3) with post-remedy of mass fraction and an ad-hoc remedy of conserved variable Q4, both temperature and mass fraction field are correct.

I would like to know if someone has encountered the similar problem or situation.

Thanks a lot.

- Ma
  Reply With Quote

Old   June 19, 2004, 00:51
Default Re: species mass fraction in reacting flow
  #5
versi
Guest
 
Posts: n/a
I met with problems of Y around 0 or 1, say -0.001 or 1.002 , but not so large as 1.2. In that case, I would first to check if there is unconvervative thing that artificailly produces mass. It may be BC, dissipation scheme, a diffusive flux that did not met conservation. etc.
  Reply With Quote

Old   June 21, 2004, 07:00
Default Re: species mass fraction in reacting flow
  #6
barbaro
Guest
 
Posts: n/a
A trick is suggested by Patankar in his book "Numerical Heat Tranfer and Fluid Flow" on pp. 144-145.
  Reply With Quote

Old   June 23, 2004, 12:12
Default Re: species mass fraction in reacting flow
  #7
Dean
Guest
 
Posts: n/a
The "tricks" applied to Y=1.2 can lead to serious loss of accuracy, up to and including unphysical behavior. This is a case where the best cure is to use decent numerical methods to begin with. For example, the species equations in strongly conservative form is essential, as is a molecular diffusion method in which the fluxes sum to exactly zero without some kind of "fudge" or "trick". There can still be problems with chemical reactions overshooting the 0 and 1 limits, but good numerical methods and sufficient temporal resolution will keep that problem under control. While the tricks that were mentioned have all been widely used, they are best avoided when possible, in my experience, as they introduce uncontrolled approximations into the numerical results. It is worth the small extra effort to use appropriate numerical methods from the start rather than patch poor ones after the fact.
  Reply With Quote

Old   June 24, 2004, 03:01
Default Re: species mass fraction in reacting flow
  #8
Rami
Guest
 
Posts: n/a
Dean,

I quite agree with you that a good scheme should be used in the first place. However, as you also mentioned, "There can still be problems with chemical reactions overshooting the 0 and 1 limits, but good numerical methods and sufficient temporal resolution will keep that problem under control". My tricks are simple correction to these deviations.

If you have better ways to control these - especially when dealing with "stiff" chemical kinetics, I would be glad to learn them from you (or any other poster).

Rami
  Reply With Quote

Old   February 20, 2013, 17:06
Default
  #9
New Member
 
Join Date: Jan 2013
Posts: 6
Rep Power: 3
Joshua_X is on a distinguished road
Hi guys?
I really like this kind of discussion. Is anybody interested in the negative mass fraction of product at the cold boundary?

Negative mass Fraction on cold boundary

anybody could help? Because to me the species equation:

d2Y/dx2-dY/dx=-w

always let the mass fraction of reactant increase and make the mass fraction of product decrease (no matter what kind of ignition concept used)
Joshua_X is offline   Reply With Quote

Reply

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
Mass fraction and volume fraction eric weddle CFX 0 September 26, 2011 05:02
Difference between Volume Fraction and Mass Fraction in CEL foo7 CFX 0 September 21, 2011 08:45
PB: Saving mass fraction of species on two points with an udf biget.a FLUENT 2 July 23, 2009 03:06
species mass fraction unbalanced ashishpande FLUENT 0 March 31, 2009 17:32
How set volume fraction and species mass fraction Reginaldo Cotto FLUENT 0 August 31, 2006 04:58


All times are GMT -4. The time now is 14:07.