[danz]

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...

[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]

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

[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

[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:

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  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

[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]

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

[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]

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:

     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

 -I$(LIB_SRC)/thermophysicalModels/basic/lnInclude \

and

    -lfluidThermophysicalModels

where added.

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

    -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]

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

[tomf]

Hi Tobi,

Yes the solver did work and the result looked reasonable.

Cheers,
happy Easter,
Tom

[lcy]

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

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

[Tobi]

Originally Posted by lcy

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 ?

[lcy]

Openfoam 1.7.0

Centos 6.3 x86_64

I donot know how to solve the problem, Thanks

Originally Posted by Tobi

Hi,

what OF version do you use ?

[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.

[hayhari]

Originally Posted by Tobi

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:

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

[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]

Originally Posted by 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

 -I$(LIB_SRC)/thermophysicalModels/basic/lnInclude \

and

    -lfluidThermophysicalModels

where added.

    -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

[babakflame]

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

[Tobi]

Originally Posted by babakflame

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:

//- 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.