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

How to understand "Y[inertIndex] = scalar(1) - Yt" in YEqn.H

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

Reply
 
LinkBack Thread Tools Display Modes
Old   September 21, 2013, 09:59
Default How to understand "Y[inertIndex] = scalar(1) - Yt" in YEqn.H
  #1
Member
 
赵庆良
Join Date: Aug 2013
Posts: 56
Rep Power: 3
zqlhzx is on a distinguished road
Hi Foamers:
As we can see the following codes in YEqn.H:

{
label inertIndex = -1;
volScalarField Yt(0.0*Y[0]);
forAll(Y, i)
{
if (Y[i].name() != inertSpecie)
{
volScalarField& Yi = Y[i];
solve
(
fvm::ddt(rho, Yi)
+ mvConvection->fvmDiv(phi, Yi)
- fvm::laplacian(turbulence->muEff(), Yi)
==
kappa*chemistry.RR(i),
mesh.solver("Yi")
);
Yi.max(0.0);
Yt += Yi;
}
else
{
inertIndex = i;
}
}
Y[inertIndex] = scalar(1) - Yt;
Y[inertIndex].max(0.0);
}
What does the "Y[inertIndex]"meas?How can I understand it?and what does the "Y[inertIndex]"use?
Thanks in advance!
zqlhzx is offline   Reply With Quote

Old   September 22, 2013, 04:55
Default
  #2
New Member
 
M. Sabouri
Join Date: Nov 2011
Posts: 22
Rep Power: 5
Moslem is on a distinguished road
I guess:
It is mass fraction of the inert species.
For example the N2 mass fraction in the following:
CH4+2(O2+3.76N2) --> CO2+2H2O+7.52N2
Moslem is offline   Reply With Quote

Old   September 22, 2013, 06:40
Default
  #3
Member
 
赵庆良
Join Date: Aug 2013
Posts: 56
Rep Power: 3
zqlhzx is on a distinguished road
Thanks Moslem,I think you are right!
zqlhzx 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



All times are GMT -4. The time now is 02:35.