[Tobi]

Hi Likun,

I tryed to change the files you mentioned but the flamelet thermodynamic is not compiling. I think I will reinvestigate time to make it working with 2.3 again.

Maybe you did more changes.
Regards Tobi

[Tobi]

Hi all,

@Likun: you are right. Only the necessary lines has to be added. There are two more functions for the boundary conditions implemented in the newest version. I did not compare the results of the tutorial with measurements but it seems to be the same.

@all: flameletModel-2.3.x is available

git clone https://github.com/shor-ty/flameletModel-2.3.x.git

Enjoy
Regards

PS: BlockMesh is not working at the moment

[Tobi]

Hi all,

blockMesh is "working" now.
Just use blockMesh and then run the added script for changing the BC.

Regards Tobi

[Likun]

Hi Tobi,

Great, thanks for your wonderful work.

I am on holiday these two weeks. I will check the new code when I go back to work.

Best regards,
Likun

[babakflame]

Dear Fellows

"Posing a point that maybe useful for others working on distinct fuels"

I have successfully used runFlameletGeneration.sh utility written by Tobi for methane and hybrid methane/hydrogen fuels. However when I tried to build PDF-Library for other fuels, The resulting PDF-Library was incapable of providing combustion. After testing Alberto 's original Look-up table tools, My PDF-Library was correct.
I have checked Tobi's code written in shell, I think its correct ( and many thanks to him), however this problem has occurred for me.

Best,

Bobi

[Tobi]

Hi Bobi,

Well my bash script is nothing special. It is only doing the same steps as you do in a manual way for different enthalpy defects... I use the same binary as in manual step. Did you change the kinetics in the bash script for your other calculations?

If the manuel steps work, then the bash script should work too (:

@all followers: I started my project for the flamelet constructor. This project is a own and personality project. I dont know if that flamelet generator will be helpful for anybody ... If it is ready and able to create flamelets I keep you posted.

If some one want to be with me and help me, let me know.

Regards
Tobi

[babakflame]

Greetings Tobi

I used the same kinetic (PolimiC1C3) , Just wanted to have combustion for propane/hydrogen as fuel. However, the process of making PDF-Library was suspiciously fast. (Maximum Temperature was 292 for all chi values). With alberto old code, I got combustion for chi values lower than quenching. (Really Don't Know why this happens)

Best

Bobi

[Tobi]

Hi Bobi,

Is it possible to send me that two cases? I want to have a look at these phenomena.

Regards Tobi

[Tobi]

Dear all and Bobi,

I had a look at my script and the problems you told. I found the problem and it is a mistake of input data. You started with an scalar dissipation rate of 0.0001 in with my flamelet script. Therefor no combustion start and all flamelets are "cold". If you set 1e-6 to the first scalar dissipation rate in "Data/Input.inp" then everything is fine. I think in the other case, when you build step by step not using my script, you have a scalar dissipationrate less than 0.0001 at the beginning. Otherwise it should not work.

Regards
Tobi

[babakflame]

Greetings Tobi

I remember that some times ago I could make PDF-Libraries with your bash script code, so that was the problem

You mean that it has been fixed in alberto's binary code that initial chi must be lower that 1e-04?

Best,

Bobi

[Tobi]

Hi,

if you have to high dissipation rate for the beginning then there is no time for reaction (chi = dissipation rate [1/s] = inverse of time) In some literature you find that chi means the time which species have for reaction. Very low chi means a lot of time for chemical reaction, high chi means short time for chemical reaction. In other literature you find a definition of energy transfer out of a cell. That means with low chi you have less energy transport out of the cell and vice versa.

So if you start building flamelets with too high chi, then there is no ignition.

I can not tell you how that code is working and how the flamelets are generated. So that are only things that I can imagine due to the fact that we only use the binarys and are not able to get into the code

Regards Tobi

[sempionato]

Hi friends,

Does anyone here have already tried to insert a lagrangian momentun source term inside a governing equation of this flamelet PDF solver? I have tryed insert a reactingparcelfoam parcel in the flamelet solver.

I'm asking this because there is a SLGThermo error in the execution. Part of the log file of this error is:

Selecting RAS turbulence model kEpsilon
kEpsilonCoeffs
{
Cmu 0.09;
C1 1.47;
C2 1.92;
C3 -0.33;
alphah 1;
alphak 1;
alphaEps 0.76923;
muLimiter on;
Lsgs 0.0002;
sigmak 1;
sigmaEps 1.3;
Prt 1;
}

Creating field dpdt

Creating field kinetic energy K

Reading flameletProperties dictionary

Preparing field Qrad (radiative heat transfer)


Constructing reacting cloud
Constructing particle forces
Selecting particle force sphereDrag
Selecting particle force gravity
Constructing cloud functions
none
Constructing particle injection models
Creating injector: model1
Selecting injection model manualInjection
Constructing 2-D injection
Selecting distribution model uniform
Selecting dispersion model none
Selecting patch interaction model standardWallInteraction
Selecting surface film model none
Selecting U integration scheme Euler
Selecting heat transfer model RanzMarshall
Selecting T integration scheme Euler


--> FOAM FATAL ERROR:
carrier requested, but object is not allocated

From function const Foam::basicMultiComponentMixture& Foam::SLGThermo::carrier() const
in file SLGThermo/SLGThermo.C at line 116.

FOAM aborting

#0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::error::abort() in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 Foam::SLGThermo::carrier() const in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libSLGThermo.so"
#3
in "/home/ghsss/OpenFOAM/ghsss-2.2.2/platforms/linux64GccDPOpt/bin/sprayFlameletPisoFoam"
#4
in "/home/ghsss/OpenFOAM/ghsss-2.2.2/platforms/linux64GccDPOpt/bin/sprayFlameletPisoFoam"
#5
in "/home/ghsss/OpenFOAM/ghsss-2.2.2/platforms/linux64GccDPOpt/bin/sprayFlameletPisoFoam"
#6 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#7
in "/home/ghsss/OpenFOAM/ghsss-2.2.2/platforms/linux64GccDPOpt/bin/sprayFlameletPisoFoam"
Aborted (core dumped)
ghsss@Ares:~/OpenFOAM/ghsss-2.2.2/run/flameletTutorials/sprayFlameletPisoFoam/Sandia_COH2N2$ sprayFlameletPisoFoam > log


--> FOAM FATAL ERROR:
carrier requested, but object is not allocated

From function const Foam::basicMultiComponentMixture& Foam::SLGThermo::carrier() const
in file SLGThermo/SLGThermo.C at line 116.

FOAM aborting

#0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::error::abort() in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 Foam::SLGThermo::carrier() const in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libSLGThermo.so"
#3
in "/home/ghsss/OpenFOAM/ghsss-2.2.2/platforms/linux64GccDPOpt/bin/sprayFlameletPisoFoam"
#4
in "/home/ghsss/OpenFOAM/ghsss-2.2.2/platforms/linux64GccDPOpt/bin/sprayFlameletPisoFoam"
#5
in "/home/ghsss/OpenFOAM/ghsss-2.2.2/platforms/linux64GccDPOpt/bin/sprayFlameletPisoFoam"
#6 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#7
in "/home/ghsss/OpenFOAM/ghsss-2.2.2/platforms/linux64GccDPOpt/bin/sprayFlameletPisoFoam"
Aborted (core dumped)
ghsss@Ares:~/OpenFOAM/ghsss-2.2.2/run/flameletTutorials/sprayFlameletPisoFoam/Sandia_COH2N2$ ^C
ghsss@Ares:~/OpenFOAM/ghsss-2.2.2/run/flameletTutorials/sprayFlameletPisoFoam/Sandia_COH2N2$ ^C
ghsss@Ares:~/OpenFOAM/ghsss-2.2.2/run/flameletTutorials/sprayFlameletPisoFoam/Sandia_COH2N2$ ^C
ghsss@Ares:~/OpenFOAM/ghsss-2.2.2/run/flameletTutorials/sprayFlameletPisoFoam/Sandia_COH2N2$

[Tobi]

Hi,

I am not familiar with parcel solvers and this question does not really belong to that thread but as you see your solver requests a field/scalar/what ever called carrier which does not exist. So please be sure that you implemented your particels as well as in parcel solvers. Maybe you have to make changes in the thermodynamics of the flamelet lib but I dont know due to the fact that I am not sure how that implementation is working

[sempionato]

Thanks Tobias, I will try to review my thermodynamics implementation.

[Tobi]

Dear Bobi,

I read that you are validating LES model in the hagens code. Can you make a comparison between the alberto thermodynamic and hagens thermodynamic code with some example -- maybe Sandy flame?

[babakflame]

Dear Tobi

I will put some comparisons as soon as I got some time (probably next week).


Best,
Bobi

[Tobi]

Normally it should be the same!

[Tobi]

Dear Bobi,

you build a LES model for the old flameletPisoFoam solver I created. Yesterday I had some time to check out some interessting things and build a pimple foam flamelet solver based on the cuoci code. Is is possible to send me your LES case that I can check out and/or modify the pimple solver for LES simulations too?

Thanks in advance
Kind regards
Tobi

[babakflame]

Greetings All

Has anybody ever tried to map the results (no matter LES or URAS) from flameletFoam (Hagen Code) onto a finer grid?

I think (due to construction of Look-up tables on grid points) this is impossible.

Any hint or comment is appreciated.

Best
Bobi

[Tobi]

Hi Bobi, that should be possible. That's not a problem of look up tables or do I misunderstand you?