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

libOpenSMOKE

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

Like Tree128Likes

Reply
 
LinkBack Thread Tools Display Modes
Old   March 12, 2013, 15:31
Default
  #141
New Member
 
Join Date: Feb 2013
Posts: 6
Rep Power: 5
danz is on a distinguished road
Hi Tobi,

Thank you for your message! In my case, I will use cantera and I will start from the EDC. From now, I just begun my internship so I study the literature and then I will let you know how does it will go...
danz is offline   Reply With Quote

Old   March 18, 2013, 03:22
Default openfoam 2.2.x
  #142
Senior Member
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Leoben (Austria)
Posts: 1,554
Blog Entries: 6
Rep Power: 27
Tobi will become famous soon enoughTobi will become famous soon enough
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
Hi all,

I tryed to use the flameletmodel in openfoam 2.2.x but there are some problems with compiling the thermodynamics. I ll try To rebuild the libraries for 2.2.x and Bild the instationary solver "turbulentFlameletRhoPimpelFoam".

Tobi
babakflame likes this.
Tobi is offline   Reply With Quote

Old   March 18, 2013, 13:00
Default
  #143
Senior Member
 
Bobi
Join Date: Oct 2012
Location: Chicago, Illinois
Posts: 372
Rep Power: 7
babakflame is on a distinguished road
Dear Tobi
Hi
I hope u will succeed in this route, however I am working on the code too. I have also worked with firefoam solver in 2.1.x . This solver uses flame surface density model which is a flamelet model that instead of integrating using PDF, employs some kind of sigma formulation (according to Tarek echekki book- Turbulent Combustion Modeling-Advances and new trends)
At the moment I do not have enough time to check the new firefoam solver (I will check it later). Nevertheless, I have read in openfoam site that in the new solver, multistep combustion has added to the firefoam. in the prior version, it only had single step. Also the cfd main solver is RhoPimple in this solver.
Regards
Bobi
babakflame is offline   Reply With Quote

Old   March 22, 2013, 09:19
Default
  #144
Senior Member
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Leoben (Austria)
Posts: 1,554
Blog Entries: 6
Rep Power: 27
Tobi will become famous soon enoughTobi will become famous soon enough
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
Hi all,

i think next week I ll get the libs for 2.2.x.

@bruno: I am sorry but I think I dont need your help
wyldckat likes this.
Tobi is offline   Reply With Quote

Old   March 23, 2013, 12:15
Default
  #145
Senior Member
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Leoben (Austria)
Posts: 1,554
Blog Entries: 6
Rep Power: 27
Tobi will become famous soon enoughTobi will become famous soon enough
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
Hi all,

I am working the last day on rebuilding the thermo-model for 2.2.x.

I added two features for output the flamelet properties (two new variables in the flameletProperties) and cleared the thermo-code (named flameletThermo)

Here is an output at the beginning:
Code:
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.2.x                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 2.2.x-020332035ede
Exec   : flamaletSimpleFoam
Date   : Mar 23 2013
Time   : 17:12:28
Host   : "cfd"
PID    : 14013
Case   : /home/shorty/OpenFOAM/shorty-2.2.x/run/SANDIA_CH4H2N2
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0


SIMPLE: convergence criteria
    field p     tolerance 0.01
    field U     tolerance 0.0001
    field h     tolerance 0.001
    field "(k|epsilon|omega)"     tolerance 0.001

Reading thermophysical properties

Selecting thermodynamics package 
{
    type            pdfFlameletThermo;
    mixture         pureMixture;
    transport       sutherland;
    thermo          hConst;
    equationOfState perfectGas;
    specie          specie;
    energy          sensibleEnthalpy;
}

/*-------------------------------------*\
|    Flamelet Thermo initialisation     |
|                                       |
|   Rebuild by Tobias Holzmann M.Eng.   |
|          www.Holzmann-cfd.de          |
\*-------------------------------------*/

Enthalpy:
     + inletfuel <fixedValue>
     + inletair <fixedValue>
     + outlet
     + axis
     + sidewall
     + burnerwall
     + front
     + back

Temperature: 
     + inletfuel <fixedValue>
     + inletair <fixedValue>
     + outlet
     + axis
     + sidewall
     + burnerwall
     + front
     + back

Mixture fraction:
     + inletfuel <fixedValue>
     + inletair <fixedValue>
     + outlet
     + axis
     + sidewall
     + burnerwall
     + front
     + back

Flamelet thermo mode is <non-adiabatic>
Flamelet thermo uses the <delta-dirac distribution> for the scalar dissipation rate
Flamelet thermo shows all single flamelet properties
Flamelet thermo shows all flamelet library properties

Flamelet property:

     + Index:                             1
     + File name:                         PDF-Library/PDF0/SR_1.bin
     + Mixture fraction points:           44
     + Variance Mixture fraction points:  30
     + Nominal strain rate:               6.91521e-310 Hz
     + Stoic. scalar dissipation rate:    1e-06 Hz
     + Temperature min:                   292.15 K
     + Temperature max:                   2124.88 K
     + Density min:                       0.152907 kg/m3
     + Density max:                       1.20255 kg/m3


Flamelet property:

     + Index:                             2
     + File name:                         PDF-Library/PDF0/SR_2.bin
     + Mixture fraction points:           44
     + Variance Mixture fraction points:  30
     + Nominal strain rate:               6.91521e-310 Hz
     + Stoic. scalar dissipation rate:    1 Hz
     + Temperature min:                   292.15 K
     + Temperature max:                   2033.91 K
     + Density min:                       0.157558 kg/m3
     + Density max:                       1.20255 kg/m3


Flamelet property:

     + Index:                             3
     + File name:                         PDF-Library/PDF0/SR_3.bin
     + Mixture fraction points:           44
     + Variance Mixture fraction points:  30
     + Nominal strain rate:               6.91521e-310 Hz
     + Stoic. scalar dissipation rate:    5 Hz
     + Temperature min:                   292.15 K
     + Temperature max:                   1989.67 K
     + Density min:                       0.159646 kg/m3
     + Density max:                       1.20255 kg/m3


Flamelet property:

     + Index:                             4
     + File name:                         PDF-Library/PDF0/SR_4.bin
     + Mixture fraction points:           44
     + Variance Mixture fraction points:  30
     + Nominal strain rate:               6.91521e-310 Hz
     + Stoic. scalar dissipation rate:    30 Hz
     + Temperature min:                   292.15 K
     + Temperature max:                   1923.46 K
     + Density min:                       0.163351 kg/m3
     + Density max:                       1.20255 kg/m3


Flamelet property:

     + Index:                             5
     + File name:                         PDF-Library/PDF0/SR_5.bin
     + Mixture fraction points:           44
     + Variance Mixture fraction points:  30
     + Nominal strain rate:               6.91521e-310 Hz
     + Stoic. scalar dissipation rate:    75 Hz
     + Temperature min:                   292.15 K
     + Temperature max:                   1868.15 K
     + Density min:                       0.167195 kg/m3
     + Density max:                       1.20255 kg/m3


Flamelet property:

     + Index:                             6
     + File name:                         PDF-Library/PDF0/SR_6.bin
     + Mixture fraction points:           44
     + Variance Mixture fraction points:  30
     + Nominal strain rate:               6.91521e-310 Hz
     + Stoic. scalar dissipation rate:    200 Hz
     + Temperature min:                   292.15 K
     + Temperature max:                   1760.13 K
     + Density min:                       0.176681 kg/m3
     + Density max:                       1.20255 kg/m3


Flamelet property:

     + Index:                             7
     + File name:                         PDF-Library/PDF0/SR_7.bin
     + Mixture fraction points:           44
     + Variance Mixture fraction points:  30
     + Nominal strain rate:               6.91521e-310 Hz
     + Stoic. scalar dissipation rate:    800 Hz
     + Temperature min:                   292.15 K
     + Temperature max:                   292.15 K
     + Density min:                       0.699245 kg/m3
     + Density max:                       1.20255 kg/m3

     + Cold flame


Flamelet library property:

     + File name:                         PDF-Library/PDF0/LookUpTable.out
     + Number of Flamelets:               6
     + Minimum stoic. scalar diss. rate:  1e-06 Hz
     + Maximum stoic. scalar diss. rate:  200 Hz
     + Temperature min:                   292.15 K
     + Temperature max:                   2124.88 K
     + Density min:                       0.152907 kg/m3
     + Density max:                       1.20255 kg/m3
     + Temperature fuel:                  292.15 K
     + Temperature oxidizer:              292.15 K
     + Density fuel:                      0.699245 kg/m3
     + Density oxidizer:                  1.20255 kg/m3
     + Enthalpy fuel:                     -998143 J/kg
     + Enthalpy oxidizer:                 -15974.1 J/kg

/*------------------------------------------------------------------------------*\
|--------->     Finished loading Look-Up-Table number: 1 out of 5      <---------|
\*------------------------------------------------------------------------------*/

.
.
.
.
.
.

Flamelet library property:

     + File name:                         PDF-Library/PDF-500/LookUpTable.out
     + Number of Flamelets:               6
     + Minimum stoic. scalar diss. rate:  1e-06 Hz
     + Maximum stoic. scalar diss. rate:  200 Hz
     + Temperature min:                   292.15 K
     + Temperature max:                   1816.43 K
     + Density min:                       0.17916 kg/m3
     + Density max:                       1.20255 kg/m3
     + Temperature fuel:                  292.15 K
     + Temperature oxidizer:              292.15 K
     + Density fuel:                      0.699245 kg/m3
     + Density oxidizer:                  1.20255 kg/m3
     + Enthalpy fuel:                     -998143 J/kg
     + Enthalpy oxidizer:                 -15974.1 J/kg

/*------------------------------------------------------------------------------*\
|--------->     Finished loading Look-Up-Table number: 5 out of 5      <---------|
\*------------------------------------------------------------------------------*/


Flamelet thermo loaded all libraries
Summarize the whole flamelet library: 

     + File name:                         PDF-Library/LookUpTable.out
     + Number of Enthalpy defects:        5
     + Minimum enthalpy defect:           0 kJ/kg
     + Maximum enthalpy defect:           -500 kJ/kg
     + Adiabatic flamelets:               1
     + Temperature fuel:                  292.15 K
     + Temperature oxidizer:              292.15 K
     + Density fuel:                      0.699245 kg/m3
     + Density oxidizer:                  1.20255 kg/m3
     + Enthalpy fuel:                     -998143 J/kg
     + Enthalpy oxidizer:                 -15974.1 J/kg
MPJ and babakflame like this.
Tobi is offline   Reply With Quote

Old   March 24, 2013, 09:39
Default
  #146
Senior Member
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Leoben (Austria)
Posts: 1,554
Blog Entries: 6
Rep Power: 27
Tobi will become famous soon enoughTobi will become famous soon enough
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
Hi all,

the new thermodynamic for the flamelet model is build in 2.2.x

I finished code developement at the moment and it seems that everything is working. Tomorrow I ll validate the solver and if everything is working fine I ll share it with you.
babakflame likes this.
Tobi is offline   Reply With Quote

Old   March 24, 2013, 12:12
Default
  #147
Senior Member
 
Bobi
Join Date: Oct 2012
Location: Chicago, Illinois
Posts: 372
Rep Power: 7
babakflame is on a distinguished road
Dear Tobi
Hi

Love u man.
I wish I were in deutschland and could buy u a beer.

Did u have any success on pimple algorithm?

Bobi
babakflame is offline   Reply With Quote

Old   March 25, 2013, 06:21
Default
  #148
Senior Member
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Leoben (Austria)
Posts: 1,554
Blog Entries: 6
Rep Power: 27
Tobi will become famous soon enoughTobi will become famous soon enough
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
Hi all,

the new flamelet model is avaiable at my homepage: http://www.holzmann-cfd.de/index.php...openfoam-2-2-x (not finished yet)

or directly on my github:

https://github.com/shor-ty/flameletModel-2.2.x

+ I added a very very brief description (have no time to build a very good one - its all my spare time )
+ I added a validation doc (I ll add this in the validation forum too)
+ Validation: http://www.cfd-online.com/Forums/ope...tml#post416171


@bobi:

yesterday I tryed to build the flameletPimpleFoam solver but without access. I think the problem is the pEqn.H file.

I realized that I have to use the old pEqn.H file from 1.7.x for getting a correct result with the flamelet thermophysical model.
tomf, wyldckat, MPJ and 1 others like this.

Last edited by Tobi; March 25, 2013 at 07:59.
Tobi is offline   Reply With Quote

Old   March 28, 2013, 11:38
Default
  #149
Senior Member
 
Tom Fahner
Join Date: Mar 2009
Location: Delft, Netherlands
Posts: 328
Rep Power: 12
tomf is on a distinguished road
Send a message via MSN to tomf Send a message via Skype™ to tomf
Hi Tobi,

Thanks a lot for all the work. I had some problems with compiling, this may be related to the relatively large amounts of updates on the 2.2.x repository.

Anyhow for some reason the "basic" folder did not compile at first, I had to change the file

Make/option into:
Code:
     EXE_INC = \
    -I$(LIB_SRC)/finiteVolume/lnInclude \
    -I$(LIB_SRC)/thermophysicalModels/specie/lnInclude \
    -I$(LIB_SRC)/thermophysicalModels/basic/lnInclude \
    -I$(LIB_SRC)/meshTools/lnInclude \
    -I$(realpath ../flameletExtraction/flamelets/turbulent/lnInclude) \
    -I$(realpath ../flameletExtraction/common/lnInclude)


LIB_LIBS = \
    -lfiniteVolume \
    -lfluidThermophysicalModels
The lines

Code:
 -I$(LIB_SRC)/thermophysicalModels/basic/lnInclude \
and
Code:
    -lfluidThermophysicalModels
where added.

And for the solver I changed the last line in the EXE_INC part into:

Code:
    -I$(realpath ../../../../thermophysicalModels/basic/flameletThermo/)
(I did not compile in the $FOAM_RUN folder, so that may be related)

At least everything compiled and the tutorial did run, but I am not entirely sure if everything is now setup correctly.

Once again thanks for your effort, much appreciated.
Tobi likes this.
tomf is offline   Reply With Quote

Old   March 29, 2013, 06:31
Default
  #150
Senior Member
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Leoben (Austria)
Posts: 1,554
Blog Entries: 6
Rep Power: 27
Tobi will become famous soon enoughTobi will become famous soon enough
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
Hi Tom,

I ll have a look at that till monday. I ll compile OF2.2.x on my laptop and will try to use the flamelet model for 2.2.x on that computer.

Maybe there are some wrong or missing entrys in the Make folder like you told.

If your solver is working everything should work

Tobi
Tobi is offline   Reply With Quote

Old   March 29, 2013, 06:50
Default
  #151
Senior Member
 
Tom Fahner
Join Date: Mar 2009
Location: Delft, Netherlands
Posts: 328
Rep Power: 12
tomf is on a distinguished road
Send a message via MSN to tomf Send a message via Skype™ to tomf
Hi Tobi,

Yes the solver did work and the result looked reasonable.

Cheers,
happy Easter,
Tom
tomf is offline   Reply With Quote

Old   March 29, 2013, 09:05
Default
  #152
lcy
New Member
 
chengyi
Join Date: Jan 2013
Posts: 3
Rep Power: 5
lcy is on a distinguished road
Hi, all
I am learning to the libopensmoke. but when I want to try to compile the libs and solver , some problem happen and not succeed to compile the solver, Do you know what happen? Thank you

Quote:
wmake libso common
+ wmake libso flamelets
+ wmake libso twoequations
+ wmake libso qmom
+ wmake libso pdfThermo
+ wmake rhoSimpleFoamFlamelets_1007
/root/OpenFOAM/OpenFOAM-1.7.0/src/finiteVolume/lnInclude/readSIMPLEControls.H: In function ‘int main(int, char**)’:
/root/OpenFOAM/OpenFOAM-1.7.0/src/finiteVolume/lnInclude/readSIMPLEControls.H:6: warning: unused variable ‘momentumPredictor’
/root/OpenFOAM/OpenFOAM-1.7.0/src/finiteVolume/lnInclude/readSIMPLEControls.H:9: warning: unused variable ‘transonic’
+ wmake rhoSimpleFoamPostProcessorSoot2E_1007
/root/OpenFOAM/OpenFOAM-1.7.0/src/finiteVolume/lnInclude/readSIMPLEControls.H: In function ‘int main(int, char**)’:
/root/OpenFOAM/OpenFOAM-1.7.0/src/finiteVolume/lnInclude/readSIMPLEControls.H:3: warning: unused variable ‘nNonOrthCorr’
/root/OpenFOAM/OpenFOAM-1.7.0/src/finiteVolume/lnInclude/readSIMPLEControls.H:6: warning: unused variable ‘momentumPredictor’
/root/OpenFOAM/OpenFOAM-1.7.0/src/finiteVolume/lnInclude/readSIMPLEControls.H:9: warning: unused variable ‘transonic’
/root/OpenFOAM/OpenFOAM-1.7.0/src/finiteVolume/lnInclude/initContinuityErrs.H:37: warning: unused variable ‘cumulativeContErr’
sootTwoEquationsOptions.H:5: warning: ‘soot_2e’ may be used uninitialized in this function
sootTwoEquationsOptions.H:5: note: ‘soot_2e’ was declared here
sootTwoEquationsOptions.H:6: warning: ‘soot_2e_uncorrelated’ may be used uninitialized in this function
sootTwoEquationsOptions.H:6: note: ‘soot_2e_uncorrelated’ was declared here

Last edited by lcy; March 29, 2013 at 09:21. Reason: when I want to try to compile the solver, some problem happen and not succeed to compile the solver, Do you know what happen?
lcy is offline   Reply With Quote

Old   March 29, 2013, 11:29
Default
  #153
Senior Member
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Leoben (Austria)
Posts: 1,554
Blog Entries: 6
Rep Power: 27
Tobi will become famous soon enoughTobi will become famous soon enough
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
Quote:
Originally Posted by lcy View Post
Hi, all
I am learning to the libopensmoke. but when I want to try to compile the libs and solver , some problem happen and not succeed to compile the solver, Do you know what happen? Thank you
Hi,

what OF version do you use ?
Tobi is offline   Reply With Quote

Old   March 29, 2013, 11:31
Default
  #154
lcy
New Member
 
chengyi
Join Date: Jan 2013
Posts: 3
Rep Power: 5
lcy is on a distinguished road
Openfoam 1.7.0

Centos 6.3 x86_64

I donot know how to solve the problem, Thanks

Quote:
Originally Posted by Tobi View Post
Hi,

what OF version do you use ?
lcy is offline   Reply With Quote

Old   March 29, 2013, 11:33
Default
  #155
Senior Member
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Leoben (Austria)
Posts: 1,554
Blog Entries: 6
Rep Power: 27
Tobi will become famous soon enoughTobi will become famous soon enough
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
Your solver should work and may be the soot solver is not working but I never used used the 1.7.0 Version and this version has some bugs in it. You should use the latest one!


But its your decision.
Tobi is offline   Reply With Quote

Old   April 1, 2013, 03:12
Unhappy
  #156
New Member
 
harikrishnan
Join Date: Sep 2012
Posts: 7
Rep Power: 6
hayhari is on a distinguished road
Quote:
Originally Posted by Tobi View Post
Tom, do you use openSMOKE with 2.1.x ?

With 1.7.x everythings work fine. But with 2.1.x I get the following error:

Code:
flameletRhoSimpleFoam: error while loading shared libraries: libgsl.so.0: cannot open shared object file: No such file or directory
maybe I did something wrong while I copied everything to other folders, but I can compile everything without errors.

hmm...


Bruno I cleared the tutorial files, and updated the header to 2.1.x
you can update your replace the old one with those but oyu have to put in the "PDF's" - its too big with that files.






Dear tobi
i am facing this problem now. Did u find the solution for this? when i updated GSL_INST_DIR i an getting a new error
rhoSimpleFoamFlamelets_1007: error while loading shared libraries: libOpenSMOKE_pdfThermo.so: cannot open shared object file: No such file or directory

Please help me solve this issue.

Thanks
Hari
hayhari is offline   Reply With Quote

Old   April 1, 2013, 05:49
Default
  #157
Senior Member
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Leoben (Austria)
Posts: 1,554
Blog Entries: 6
Rep Power: 27
Tobi will become famous soon enoughTobi will become famous soon enough
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
Hi Harry,

1. you do not have to compile/install GSL
2. use the Flamelet-Model 2.1.x on my git repository. Its newer and some user-unfriendly setups are made there.

3. if you have the possibility use OpenFOAM 2.2.x couse I made a few output changes in the Flamelet-Model 2.2.x

https://github.com/shor-ty?tab=repositories
Tobi is offline   Reply With Quote

Old   April 1, 2013, 07:23
Default
  #158
Senior Member
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Leoben (Austria)
Posts: 1,554
Blog Entries: 6
Rep Power: 27
Tobi will become famous soon enoughTobi will become famous soon enough
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
Quote:
Originally Posted by tomf View Post
Hi Tobi,

Thanks a lot for all the work. I had some problems with compiling, this may be related to the relatively large amounts of updates on the 2.2.x repository.

Anyhow for some reason the "basic" folder did not compile at first, I had to change the file
Code:
 -I$(LIB_SRC)/thermophysicalModels/basic/lnInclude \
and
Code:
    -lfluidThermophysicalModels
where added.

Code:
    -I$(realpath ../../../../thermophysicalModels/basic/flameletThermo/)
(I did not compile in the $FOAM_RUN folder, so that may be related)

At least everything compiled and the tutorial did run, but I am not entirely sure if everything is now setup correctly.

Once again thanks for your effort, much appreciated.
Hi Tom,

I compiled the libs and solver on my laptop now.
There were some problems like you told.

But the additional libFluid... i did not have to add

Well I updated the flamelet libraries on github and now it should work

Thanks for your feedback
Tobi is offline   Reply With Quote

Old   April 3, 2013, 04:44
Default
  #159
Senior Member
 
Bobi
Join Date: Oct 2012
Location: Chicago, Illinois
Posts: 372
Rep Power: 7
babakflame is on a distinguished road
Hi Tobi

I have a question about the code results. As u know, the code gives the mass fraction of any specie that we want. However, for computing emission index I need the mole fraction of CO and CO2 at a certain position. I know that in literature we have a correlation for commuting mass to mole fraction. But, as u know in that correlation in addition to the molecular weight of the considering specie we need to have the molecular weight of the mixture which is now the products molecular weight. I think since the code computes species mass fractions, therefore it calculates the molecular weight of products as a whole mixture. So, would you plz hint me how can I exploit CO and CO2 mole fractions from the results or maybe change the code to get mole fractions data?
Regards
Bobi

Last edited by babakflame; April 3, 2013 at 05:11.
babakflame is offline   Reply With Quote

Old   April 3, 2013, 06:14
Default
  #160
Senior Member
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Leoben (Austria)
Posts: 1,554
Blog Entries: 6
Rep Power: 27
Tobi will become famous soon enoughTobi will become famous soon enough
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
Quote:
Originally Posted by babakflame View Post
Hi Tobi

I have a question about the code results. As u know, the code gives the mass fraction of any specie that we want. However, for computing emission index I need the mole fraction of CO and CO2 at a certain position. I know that in literature we have a correlation for commuting mass to mole fraction. But, as u know in that correlation in addition to the molecular weight of the considering specie we need to have the molecular weight of the mixture which is now the products molecular weight. I think since the code computes species mass fractions, therefore it calculates the molecular weight of products as a whole mixture. So, would you plz hint me how can I exploit CO and CO2 mole fractions from the results or maybe change the code to get mole fractions data?
Regards
Bobi
Hi Bobi,

you are right. The code gives you the mass fraction. In my thesis I used the important species and the molecular weight for that to calculate the mole fraction in paraview. The problem is, that the flamelet-solver and thermodynamic is not calculating the mass fraction. The thermodynamics just look for the mass fraction into the Look-Up-Tables and there for it was not easy for me to build a new function that is calculating that.

It 's possible to do that in the thermodynamics, so that you build new functions or stuff like that.

The problem is here, that you have do many things:

1. the thermodynamics just realize the species you want to output:
Code:
//- Initialise the wanted mass fraction fields
        omega_.setSize(flamelets_library.number_of_species()+1);

        for (int j=0;j<flamelets_library.number_of_species()+1;j++)
        {
            if(j < flamelets_library.number_of_species())
            {
                std::string name_of_species = "omega_" + flamelets_library.species()[j+1];

                omega_.set
                (
                    j,
                    new volScalarField
                    (
                        IOobject
                        (
                            name_of_species,
                            mesh.time().timeName(),
                            mesh,
                            IOobject::NO_READ,
                            IOobject::AUTO_WRITE),
                            mesh,
                            dimensionedScalar("zero", dimensionSet(0,0,0,0,0,0,0),
                            0.0
                        )
                    )
                );
            }
        }
2. you have no molecular weights in the thermodynamics
3. its annoying that we have no chance to have a look into the flamelet generatin code


So what you have to do is:

1. set up a function that get all involved species for you case
2. set up a list of constants (molecular weights of ALL species)
3. building a new function calculate_molefraction() (the simplest thing)



Well maybe the mole fraction is also avaiable in the look-up-table but I dont know how to have a look into that. If there are the mole fractions included, it would be the easiest way to get the values.
But as I told befor... its annoying that we dont have the chance to have a look into the flamelet generation!

I would build a own free one but I have no time for that.
babakflame likes this.
Tobi 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 17:20.