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

Janaf coefficient for gaz mixture + nMoles ?

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

Like Tree3Likes
  • 1 Post By Uyan
  • 2 Post By Uyan

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   August 8, 2016, 09:59
Default Janaf coefficient for gaz mixture + nMoles ?
  #1
Member
 
AdOo
Join Date: Mar 2016
Location: Bordeaux
Posts: 91
Rep Power: 7
adrieno is on a distinguished road
Hi all,

I would like to change the gas in the solver XiFoam, (running methane instead of propane). The thing is, I don't know how to calculate the Janaf coefficient for gas mixture. Does anybody know ?
Another question more specific about this solver I think, WHY is there nMoles = 1 for products, whereas it's set to 24.8095 for reactants (in the combustionProperties file). Where does this 24.8095 comes from ?

Thanks a lot, hopping one of you hold the answer,
Adrien
adrieno is offline   Reply With Quote

Old   August 10, 2016, 08:11
Default
  #2
Member
 
AdOo
Join Date: Mar 2016
Location: Bordeaux
Posts: 91
Rep Power: 7
adrieno is on a distinguished road
Ok, about the gas mixture,I've found my answer, using the Dalton's law in order to determine the janaf coeff of a gas mixture.

https://en.wikipedia.org/wiki/Dalton%27s_law

Leading to : Cp = xi . Cpi where xi is the mole fraction.

Unfortunately, I can't handle why there's nmoles = 24.8095 for reactants and nmoles = 1 for products...(According to me, they all should be at 1...) If someone have the answer, let me know.

Have a nice ay,
Adrien
adrieno is offline   Reply With Quote

Old   August 13, 2016, 15:55
Default nmoles ?
  #3
Member
 
Join Date: Feb 2014
Posts: 61
Rep Power: 9
Uyan is on a distinguished road
Hi adrieno,

To change the cp coefficients you can use the adiabaticFlameT utility.

cd $FOAM_UTILITIES/thermophysical/adiabaticFlameT
adiabaticFlameT controlDict

you ll get values for the fuel specified in controlDict.


However I have the same problem about 'nmoles' could not find an answer.
If you find an answer please share.
Uyan is offline   Reply With Quote

Old   August 16, 2016, 05:01
Default
  #4
Member
 
AdOo
Join Date: Mar 2016
Location: Bordeaux
Posts: 91
Rep Power: 7
adrieno is on a distinguished road
Hi Uyan,

I don't understand how to use the controlDict you've told me about :


Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.2.2                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    note        "settings for calculating the adiabatic flame temperature";
    class       dictionary;
    object      controlDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
P       1e5;
T0      300.0;
H2
{
    fuel H2;
    n    0;
    m    2;
}
CH4
{
    fuel CH4___ANHARMONIC;
    n    1;
    m    4;
}
ETHYLENE
{
    fuel C2H4;
    n    2;
    m    4;
}
PROPANE
{
    fuel C3H8;
    n    3;
    m    8;
}
OCTANE
{
    fuel C8H18(L)_isooctane;
    n    8;
    m    18;
}
N-HEPTANE
{
    fuel C7H16_n-heptane;
    n    7;
    m    16;
}
$H2;
// $CH4;
// $PROPANE;
// ************************************************************************* //
Could you be more specific please, sorry for my ignorance.

Another question : Do you know how to adapt the coeff in the combustionPorperties file ? for another gas like methane in my case (I just handle the laminar speed with Su but I'm a bit lost with all the other coefficients)

Thanks a lot,
Adrien
adrieno is offline   Reply With Quote

Old   August 16, 2016, 21:31
Default
  #5
Member
 
Join Date: Feb 2014
Posts: 61
Rep Power: 9
Uyan is on a distinguished road
Hi Adrian,

Type these two commands on an OpenFOAM terminal,

Code:
cd $FOAM_UTILITIES/thermophysical/adiabaticFlameT 
adiabaticFlameT controlDict

this will print out cp coefficients you need for the fuel specified in the controlDict. In the controlDict you have shown here it is H2.

Well for other values in combustionProperties dictionary you ll have to look for values from literature and try them out.
arvindpj likes this.
Uyan is offline   Reply With Quote

Old   August 17, 2016, 05:26
Default
  #6
Member
 
AdOo
Join Date: Mar 2016
Location: Bordeaux
Posts: 91
Rep Power: 7
adrieno is on a distinguished road
Thank you very much for your time Uyan,

Does it gives the Cp value for a stoechiometric mixture ? (This is what I'm looking for).
And just to know, is it possible to have it for another stoechiometric value ?

For the other coefficients, this is what I was trying to do... but it's hard to find especially as the Weller combustion model doesn't seem very common. Even the laminar flame speed seems to differ from a revue to another ! Aaaa combustion is such a complex science for neophyte !

Thank you again Uyan.
adrieno is offline   Reply With Quote

Old   August 24, 2016, 09:14
Default nMoles
  #7
Member
 
Join Date: Feb 2014
Posts: 61
Rep Power: 9
Uyan is on a distinguished road
Hi Adrian,

For you and others who wondered why the "nMoles" is not equal to 1 in the reactants side but 24.8. I finally found the answer how to calculate. So I post it here for anyone who is stuck there.

In the XiFoam tutorial file (moriyoshiHomogeneous), thermophysicalProperties dictionary contains properties for reactants and products as following.

Code:
reactants
{
    specie
    {
        nMoles          24.8095;
        molWeight       29.4649;
    }

   .....

Code:
products
{
    specie
    {
        nMoles          1;
        molWeight       28.3233;
    }
This case is for a stoichoemetric mixture of Propane and air. So Propane and Air mixture is considered as reactants.

Reactants ---------------> Products

C_3H_8 + 5\bigg( O_2 + 3.76 N2 \bigg) -----> 3CO_2 + 4H_2O + 18.8 N2

For reactants
nMoles = 1 + 5 + 18.8 = 24.8 \sim 24.8095.

molWeight = \frac{44 + 5\times (32  + (3.76 \times 28) ) }{1 + 5\times(1+3.76)} \sim 29.4

For products it is considered as one species, so nMoles = 1
Uyan is offline   Reply With Quote

Old   August 24, 2016, 09:23
Default
  #8
Member
 
AdOo
Join Date: Mar 2016
Location: Bordeaux
Posts: 91
Rep Power: 7
adrieno is on a distinguished road
Thank you Uyan, that's what I've done finally in my simulation, but that was just a guess finding the way the weller combustion model was working...
adrieno is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
compile of new combination of property in thermophisicalpropery milad653279 OpenFOAM Running, Solving & CFD 4 April 6, 2017 04:35
Janaf (Old ?) coefficient adrieno OpenFOAM Pre-Processing 1 April 18, 2016 07:07
Ansys CFX problem: unexpected very high temperatures in premix laminar combustion faizan_habib7 CFX 4 February 1, 2016 17:00
Question about heat transfer coefficient setting for CFX Anna Tian CFX 1 June 16, 2013 06:28
Automotive test case vinz OpenFOAM Running, Solving & CFD 98 October 27, 2008 08:43


All times are GMT -4. The time now is 17:06.