# How to get partial pressure of oxygen?

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

 March 2, 2010, 08:49 How to get partial pressure of oxygen? #1 Senior Member     Sebastian Gatzka Join Date: Mar 2009 Location: Frankfurt, Germany Posts: 729 Rep Power: 12 Hello World. Can anybody tell me how I can access the partial pressure of a species / componente? In my case I need the partial pressure of oxgen! Cheers. S. __________________ Schrödingers wife: "What did you do to the cat? It's half dead!"

 March 2, 2010, 18:52 #2 New Member   Vinay Join Date: Nov 2009 Location: Nagpur Posts: 21 Rep Power: 9 i too need the same!

March 3, 2010, 02:51
#3
Senior Member

Ji Junjie
Join Date: Feb 2010
Location: Shanghai, China
Posts: 109
Rep Power: 8
Quote:
 Originally Posted by sega Hello World. Can anybody tell me how I can access the partial pressure of a species / componente? In my case I need the partial pressure of oxgen! Cheers. S.
You can use the following equation
pi=C_P(c,t)*gas_molar_fraction
where, pi is the partial pressure, C_P(c,t) gets the mixed gas pressure.
gas_molar_fraction is the mole fraction of a species in the mixture, which can be calculated by
- first summing up [mass fraction / molecular weight] of all species and
- then dividing the same quotient for the individual species by the calculated sum.
This is done automatically by the following function call (including variables and preparation):
Material *mix_mat = mixture_material(Get_Domain(1));
Material *spe_mat = NULL;
real all_mass_fracts[MAX_SPE_EQNS];
real all_mole_fracts[MAX_SPE_EQNS];
int i = -1;
mixture_species_loop(mix_mat, spe_mat, i)
{
all_mass_fracts[i] = C_YI(c,t,i);
}
Mole_Fraction(mix_mat, all_mass_fracts, all_mole_fracts)
After this, the mole fractions of all species can be found in the array "all_mole_fracts".

 March 4, 2010, 11:31 #4 Senior Member     Sebastian Gatzka Join Date: Mar 2009 Location: Frankfurt, Germany Posts: 729 Rep Power: 12 Great, Thank you. This leads me to the question to determine which position in the array corresponds to which species? Any idea how I can know this? __________________ Schrödingers wife: "What did you do to the cat? It's half dead!"

March 4, 2010, 11:59
#5
Senior Member

Ji Junjie
Join Date: Feb 2010
Location: Shanghai, China
Posts: 109
Rep Power: 8
Quote:
 Originally Posted by sega Great, Thank you. This leads me to the question to determine which position in the array corresponds to which species? Any idea how I can know this?
The index of the first species in the list of the Fluent interface is 0.
If the list in the Fluent is : O2, CO, CO2......, then
all_mass_fracts[0] means the mass fraction of O2, all_mass_fracts[1] means that of CO , and so on.

March 4, 2010, 12:07
#6
Senior Member

Sebastian Gatzka
Join Date: Mar 2009
Location: Frankfurt, Germany
Posts: 729
Rep Power: 12
Quote:
 Originally Posted by gearboy The index of the first species in the list of the Fluent interface is 0. If the list in the Fluent is : O2, CO, CO2......, then all_mass_fracts[0] means the mass fraction of O2, all_mass_fracts[1] means that of CO , and so on.
Ok. I got this.
And how do I know the actual list?
__________________
Schrödingers wife: "What did you do to the cat? It's half dead!"

 March 4, 2010, 12:29 #7 Senior Member     Sebastian Gatzka Join Date: Mar 2009 Location: Frankfurt, Germany Posts: 729 Rep Power: 12 Ok. I put this so far into my code Code: Material *mixtureMaterial = mixture_material(Get_Domain(1)); Material *speciesMaterial = NULL; real massFractions[MAX_SPE_EQNS]; real moleFractions[MAX_SPE_EQNS]; real partialPressures[MAX_SPE_EQNS]; int i = -1; mixture_species_loop(mixtureMaterial, speciesMaterial, i) { massFractions[i] = C_YI(cell,thread,i); } Mole_Fraction(mixtureMaterial, massFractions, moleFractions); int i = -1; mixture_species_loop(mixtureMaterial, speciesMaterial, i) { partialPressures[i] = C_P(cell,thread) * moleFractions[i]; } Now the code is killing me. It want's a ; bevor the second int i = -1. But there is one - right after Mole_Fraction( ... ). So what may be wrong? __________________ Schrödingers wife: "What did you do to the cat? It's half dead!"

 March 7, 2010, 05:55 How to get partial pressure of oxygen? #8 New Member   Jamshid Montazer Join Date: Mar 2010 Location: Norway Posts: 1 Rep Power: 0 Hi, Well, as new member in OpenFoam (days of experience in OpenFoam) I would say using dsmcFoam (Direct simulation Monte Carlo (DSMC) solver for 3D, transient, multi- species flows) might help Back in my Fluent experience, you have to use the Multiphase model and depends to the nature of component (mixture, Eulerian or VOF) will works fine. Looks the same approach is working in OpenFoam: -compressibleInterFoam: Solver for 2 compressible, isothermal immiscible fluids using a VOF (volume of fluid) phase-fraction based interface capturing approach, ) More could be find at Standard solvers: http://www.openfoam.com/features/standard-solvers.php __________________ Dr.Montazer Norway

March 7, 2010, 06:02
#9
Senior Member

Sebastian Gatzka
Join Date: Mar 2009
Location: Frankfurt, Germany
Posts: 729
Rep Power: 12
Quote:
 Originally Posted by Dr.Montazer Hi, Well, as new member in OpenFoam (days of experience in OpenFoam) I would say using dsmcFoam (Direct simulation Monte Carlo (DSMC) solver for 3D, transient, multi- species flows) might help Back in my Fluent experience, you have to use the Multiphase model and depends to the nature of component (mixture, Eulerian or VOF) will works fine. Looks the same approach is working in OpenFoam: -compressibleInterFoam: Solver for 2 compressible, isothermal immiscible fluids using a VOF (volume of fluid) phase-fraction based interface capturing approach, ) More could be find at Standard solvers: http://www.openfoam.com/features/standard-solvers.php
I'm not sure this is related to the topic.
__________________
Schrödingers wife: "What did you do to the cat? It's half dead!"

August 8, 2011, 04:12
#10
New Member

Join Date: Nov 2010
Posts: 19
Rep Power: 8
Quote:
 Originally Posted by gearboy You can use the following equation pi=C_P(c,t)*gas_molar_fraction where, pi is the partial pressure, C_P(c,t) gets the mixed gas pressure. gas_molar_fraction is the mole fraction of a species in the mixture, which can be calculated by - first summing up [mass fraction / molecular weight] of all species and - then dividing the same quotient for the individual species by the calculated sum. This is done automatically by the following function call (including variables and preparation): Material *mix_mat = mixture_material(Get_Domain(1)); Material *spe_mat = NULL; real all_mass_fracts[MAX_SPE_EQNS]; real all_mole_fracts[MAX_SPE_EQNS]; int i = -1; mixture_species_loop(mix_mat, spe_mat, i) { all_mass_fracts[i] = C_YI(c,t,i); } Mole_Fraction(mix_mat, all_mass_fracts, all_mole_fracts) After this, the mole fractions of all species can be found in the array "all_mole_fracts".
Hi every body

I want to have my species' partial pressure but the problem is that I have back flow in my simulating.(I simulate a 2D bubble column). so the static pressure that I recive from (c_P()) is negetive in some cells.
do you think using the absolute of static pressure is ok? my codeis becom:
Po2= abs(C_P(c,t))*mol_fraction_o2

 Tags oxygen, partial pressure

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post Mihail CFX 7 September 7, 2014 06:27 Gearb0x OpenFOAM Running, Solving & CFD 2 February 28, 2011 22:18 jai CFX 7 December 12, 2008 04:11 Chris CD-adapco 2 June 21, 2005 08:20 DS & HB Main CFD Forum 0 January 8, 2000 16:00

All times are GMT -4. The time now is 11:01.