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

[openSmoke] libOpenSMOKE

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

Like Tree133Likes

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   February 4, 2013, 07:55
Default
  #121
Super Moderator
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Tussenhausen
Posts: 2,708
Blog Entries: 6
Rep Power: 51
Tobi has a spectacular aura aboutTobi has a spectacular aura aboutTobi has a spectacular aura about
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
Quote:
Originally Posted by babakflame View Post
Dear Tobi
Hi
Would u plz hint me that in the Sandia CH4H2N2 2D case, How does the values for k & epsilon be calculated for the inlet air and fuel? (Forexample the following lines

inletfuel
{
type turbulentIntensityKineticEnergyInlet;
intensity 0.2;
U U;
phi phi;
value uniform 239.65;
}

For instance, If I change the fuel from CH4 to C3H8 without any changes in velocity magnitude, Does these values change?
Regards
Bobi
Hi Bobi,

this value is calculated with the cfd-online calculation tool: http://www.cfd-online.com/Tools/turbulence.php

For that you know the free-stream velocity and the turbulent length scale (d_hyd * 0.07). At least you have to set the intensity. 20% for high turbulent inlets.

Then you get the values for k epsilon and omega.
Tobi is offline   Reply With Quote

Old   February 12, 2013, 13:28
Default
  #122
Senior Member
 
Bobby
Join Date: Oct 2012
Location: Michigan
Posts: 454
Rep Power: 15
babakflame is on a distinguished road
Dear Tobi
Hi
I am trying to understand the mixture fraction and it's variance in the code:
Mixture Fraction:
(
fvm::div(phi, csi)
- fvm::Sp(fvc::div(phi), csi)
- fvm::laplacian(turbulence->muEff()/sigmat, csi)
)
The equation has different forms in different references, lets start with your equation:

Thanks buddy , I got a problem with internet, I mean exactly what you have written, would u plz tell me that how exactly each term in the code is related to the equation, I have read the openFoam user guide , However I can not derive in openFoam format exactly as it is written in the code,
I have problem with sigmat (How does it appear) and also fvm::SP()

Last edited by babakflame; February 12, 2013 at 13:44.
babakflame is offline   Reply With Quote

Old   February 12, 2013, 13:32
Default
  #123
Super Moderator
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Tussenhausen
Posts: 2,708
Blog Entries: 6
Rep Power: 51
Tobi has a spectacular aura aboutTobi has a spectacular aura aboutTobi has a spectacular aura about
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
Hi,

is your sentence complete ?
That equation is a normal conserved scalar euqation?

What do you need an where are the differences?
But that Code is like that one:


\rho \frac{\partial \tilde Z}{\partial t} + \tilde  \rho \tilde v \cdot \nabla \tilde Z = \nabla (\tilde \rho D \nabla \tilde Z - \tilde \rho \widetilde{ v'' Z''})

The term on the RHS is the diffusion term and is modelled in RANS with the viscosity turbulence model. So the turbulent diffusion term is:

\widetilde{v'' Z''} = - D_t \nabla \tilde Z

now it follows that the diffusion term on the RHS is:

\nabla (\tilde \rho D \nabla \tilde Z  + \tilde \rho \tilde D_t \nabla \tilde Z) = \nabla (\tilde \rho D \nabla \tilde Z  + \tilde \rho \frac{\mu_t}{\sigma_t}\nabla \tilde Z)

Most references neglect the molecular diffusivity couse it has no influence: it follows:

\nabla (\tilde \rho D \nabla \tilde Z  + \tilde \rho  \frac{\mu_t}{\sigma_t}\nabla \tilde Z) = \nabla ( \tilde \rho  \frac{\mu_t}{\sigma_t}\nabla \tilde Z)

But I implemented it with the molecular diffusion so you get the following line:

\nabla (\tilde \rho D \nabla \tilde Z  +  \frac{\mu_t}{\sigma_t}\nabla \tilde Z) = \nabla (  \frac{\mu_{eff}}{\sigma_t}\nabla \tilde Z)

Then your equation is (with no time derivation -steady state):

\tilde  \rho \tilde v \cdot  \nabla \tilde Z = \nabla (\frac{\mu_{eff}}{\sigma_t} \nabla \tilde Z )

and that equation is that one I (and Alberto) implemented in the flamelet solver. Alberto use the mu_t for that. I am using the mu_eff couse I was developing an laminar flamelet solver that uses the molecular diffusion. Thats the reason for it. That flamelet solver is build for turbulent flows so mu_lam <<<< mu_turbulent and therefor mu_eff =~ mu_turbulent
TrII4d and babakflame like this.

Last edited by Tobi; February 16, 2013 at 09:56. Reason: forgot rho
Tobi is offline   Reply With Quote

Old   February 12, 2013, 14:24
Default
  #124
Senior Member
 
Bobby
Join Date: Oct 2012
Location: Michigan
Posts: 454
Rep Power: 15
babakflame is on a distinguished road
Many thanks Tobi , for ur detailed explanations, I got exactly what happens to the R.H.S. , However can I ask to explain me how does
ρ (V . grad Z) is divided in to
fvm::div(phi, csi) - fvm::Sp(fvc::div(phi), csi)

I think we can use this tensor expression: ∂_i (ρ〖V)〗_i=ρ ∂_i V_i+ V_i ∂_i ρ

However, I am still trying to write it in OpenFoam format


Thanks
Bobi

Last edited by babakflame; February 12, 2013 at 15:21.
babakflame is offline   Reply With Quote

Old   February 12, 2013, 15:16
Default
  #125
Super Moderator
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Tussenhausen
Posts: 2,708
Blog Entries: 6
Rep Power: 51
Tobi has a spectacular aura aboutTobi has a spectacular aura aboutTobi has a spectacular aura about
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
Quote:
Originally Posted by babakflame View Post
Many thanks Tobi , for ur detailed explanations, I got exactly what happens to the R.H.S. , However can I ask to explain me how does
ρ (V . grad Z) is divided in to
fvm::div(phi, csi) + fvm::Sp(fvc::div(phi), csi)
Thanks
Bobi
Hi Bobi,

no problem - its my topic and want to help everyone. Step by step:

\tilde \rho \tilde v \nabla \tilde Z

is

Code:
fvm::div(phi, csi)      // its the convection term
fvm is standing in the programmers guide:
Code:
Implicit / Explicit    ::::   fvm:: / fvc::
I think thats clear. So the convection term is treated implicit.

The diffusion term:

\nabla (\frac{\mu_{eff}}{\sigma_t} \nabla \tilde Z)

is

Code:
fvm::laplacian(turbulence->muEff()/sigmat, csi)
Thats it. The other term:
Code:
fvm::Sp(fvc::div(phi), csi)
is a numerical trick for stabilisation.
Lets consider the following equation:

\frac{dU}{dt} = \underline{A} + \underline{B}U

The underlined variables are matrixes. Okay if you use the eulerian procedure to solve that problem you get (explicit):

U^n - U^{n-1} = dt \cdot \underline{A} + dt \cdot \underline{B} U^{n-1}

You want to know U^(n), then you get:

U^n = dt \cdot \underline{A} + dt \cdot [\underline{B} + 1 ] U^{n-1}

You calculate U^n with the known variables. That method (explicit) is in some cases very instable so that you use an implicit method for the matrix B (implicit) :

U^n - U^{n-1} = dt \cdot \underline{A} + dt \cdot \underline{B}U^{n}

after transform that equation you get:

[1 - dt \underline{B}] U^n = dt \cdot \underline{A} + U^{n-1}

Now its possible to make a stabilisation for explicit source terms (here the matrix A). We go back to the first equation:

\frac{dU}{dt} = \underline{A} + \underline{B}U = \underline{A} \cdot 1 + \underline{B} U

You can see that the matrix A is multiplyed by 1. So now the trick:


\frac{dU}{dt} = \underline{A} + \underline{B}U = \underline{A} \cdot 1 + \underline{B} U = \underline{A} \cdot \frac{U}{U} + \underline{B} U

Go back to the implicit eulerian procedure:

U^n - U^{n-1} = dt \cdot \underline{A} \cdot 1 + dt \cdot \underline{B}U^{n} = dt \cdot \underline{A} \cdot \frac{U^n}{U^{n-1}} + dt \cdot \underline{B}U^{n}

Now you see that the matrix A is not mulitplyed by 1 anymore but with a complete different value. That gives you the possibility to linearize explicit source terms (and thats the trick).

If you reach steady state conditions U^n/U^n-1 = 1 and everything is fine again

And thats what the following line is doing:
Code:
fvm::Sp(fvc::div(phi), csi)
Tobi

PS: sigmat = turbulent prandtl number

Last edited by Tobi; February 12, 2013 at 15:35. Reason: Adding desciripton of sigmat
Tobi is offline   Reply With Quote

Old   February 16, 2013, 07:53
Default
  #126
Senior Member
 
Bobby
Join Date: Oct 2012
Location: Michigan
Posts: 454
Rep Power: 15
babakflame is on a distinguished road
Dear My friend Tobi
Hi

Many Thanks for all your earlier helps. I have started changing the code to LES .
In the csiEqn.H file I need to rewrite the stoichiometric scalar dissipation rate formula for LES.
In LES we have: chi= 2(Dm+Dt) (GradZ)^2
Cause I am still novice with openFOAM format , would you plz check out that I have written the mentioned formula right?

// Updating stoichiometric scalar dissipation rate

chi_st = 2*(turbulence->muEff()/(sigmat*rho))*magSqr(fvc::grad(csi));



Regards
Bobi
babakflame is offline   Reply With Quote

Old   February 16, 2013, 09:59
Default
  #127
Super Moderator
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Tussenhausen
Posts: 2,708
Blog Entries: 6
Rep Power: 51
Tobi has a spectacular aura aboutTobi has a spectacular aura aboutTobi has a spectacular aura about
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
Hi,

I had a look at my first equations.
I think I missed the rho in the turbulent diffusion (i corrected it now but I am not 100% sure).

I ll have a look at my litrature.

To you.
Why do you use (sigmat * rho) ?

And magSqr is wrong... you need POW()
Tobi is offline   Reply With Quote

Old   February 17, 2013, 06:00
Default
  #128
Senior Member
 
Bobby
Join Date: Oct 2012
Location: Michigan
Posts: 454
Rep Power: 15
babakflame is on a distinguished road
Dear Tobi
Hi
Thanks For POW()
About Rho: (would u plz check out my relation)
As we have
(Dm+Dt) = (Deff) = (mu-eff)/(rho* Schmidt) = (if we cosider Le =1, then we have)

= (mu-eff) / (rho * Prandtl)
Tobi likes this.
babakflame is offline   Reply With Quote

Old   February 18, 2013, 07:28
Default
  #129
Super Moderator
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Tussenhausen
Posts: 2,708
Blog Entries: 6
Rep Power: 51
Tobi has a spectacular aura aboutTobi has a spectacular aura aboutTobi has a spectacular aura about
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
Hi Bobi,

I had a look into the reference: B. Cuenot: The Flamelet Model for Non-Premixed Combustion, Springer Science+Business Media B.V (2011)

and into my masterthesis (references Ferziger and Peric).

Due to the RANS turbulence modelling you get the following for the navier stokes eqn (Ferziger):

\frac {\partial (\rho \bar u_i)}{\partial t} + \frac{\partial}{\partial x_i}(\rho \bar u_i \bar u_j + \rho \overline{u_i'u_j'}) = -\frac{\partial \bar p}{\partial x_i} + \frac{\partial \bar \tau_{ij}}{\partial x_j}


Now you can see the same "molecular and turbulent diffusion" there. The diffusion term is modeled by:

-\rho\overline{u_i'u_j'} = \mu_t (\frac{\partial \bar u_i}{\partial x_j} + \frac{\partial \bar u_j}{\partial x_i}) - \frac{2}{3} \rho \delta_{ij} k

where k is the turbulent kinematic energy.

Now you have to know that mu_t is calculated by:

\mu_t = C_\mu \rho q L = C_\mu \rho \sqrt{2k} L

As you can see in the variable "mu_t" there is rho.
So if you are using mu_t you have already rho. If you neglegt the molecluar diffusion:

\nabla (\bar \rho D \nabla \tilde Z - \bar \rho \widetilde{v''Z''}) = \nabla (\frac{\mu_t}{\sigma_t} \nabla \tilde Z)

So if you are using mu_eff I think the rho is invoken in the calculation of mu_eff. You should have a look into the file where mu_eff is calculated.

As far as I know the following equation should be right:


\nabla (\bar \rho D \nabla \tilde Z - \bar \rho \widetilde{v''Z''}) = \nabla (\frac{\mu_{eff}}{\sigma_{eff}} \nabla \tilde Z)
babakflame likes this.

Last edited by Tobi; February 18, 2013 at 17:10.
Tobi is offline   Reply With Quote

Old   February 18, 2013, 07:44
Default
  #130
Senior Member
 
Bobby
Join Date: Oct 2012
Location: Michigan
Posts: 454
Rep Power: 15
babakflame is on a distinguished road
Dear Tobi

I will look into the references You have mentioned, very soon. Thanks For all buddy.
babakflame is offline   Reply With Quote

Old   February 18, 2013, 07:51
Default
  #131
Super Moderator
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Tussenhausen
Posts: 2,708
Blog Entries: 6
Rep Power: 51
Tobi has a spectacular aura aboutTobi has a spectacular aura aboutTobi has a spectacular aura about
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
Quote:
Originally Posted by babakflame View Post
Dear Tobi

I will look into the references You have mentioned, very soon. Thanks For all buddy.

Ferziger Page: 346 - 349
The Paper of Cuenot:

Varianz Z'' of the mixture fraction:

\bar \rho \frac{\partial \tilde Z''}{\partial t} + \bar \rho \tilde v \nabla \tilde Z'' = \nabla \cdot (\bar \rho D_t \nabla \tilde Z'') + 2 \bar \rho D_t (\nabla \tilde Z)^2 - \rho \tilde \chi

Compared with the c++ code in alberto cuocis flamelet model you can see that the term:

\bar \rho D_t

is modelled with

\bar \rho D_t = \frac {\mu_t}{\sigma_t}

Page 49. In that section you can read the following sentence: "The molecular diffusivity is often neglected as it is much smaller than the turbulent diffusivity D_t" (like I told befor). Hope it helps you.
Tobi
babakflame likes this.
Tobi is offline   Reply With Quote

Old   February 22, 2013, 12:04
Default
  #132
Senior Member
 
Bobby
Join Date: Oct 2012
Location: Michigan
Posts: 454
Rep Power: 15
babakflame is on a distinguished road
Dear Tobi
I have made the midifications needed for your code to be changed into LES. However , when I run the code in serial, I got this strange error:

Quote:
Creating turbulence model

Selecting turbulence model type LESModel
Selecting LES turbulence model oneEqEddy
oneEqEddyCoeffs
{
ce 1.048;
Prt 1;
ck 0.094;
}

Creating field dpdt

Creating field kinetic energy K

Reading flamelet dictionary

Preparing field Qrad (radiative heat transfer)

Courant Number mean: 3.01006e-05 max: 0.0379248

PIMPLE: Operating solver in PISO mode


Starting time loop

Reading/calculating field UMean

Reading/calculating field pMean

Reading/calculating field UPrime2Mean

Reading/calculating field pPrime2Mean

fieldAverage: starting averaging at time 0

Courant Number mean: 3.01006e-05 max: 0.0379248
deltaT = 1.2e-06
Time = 1.2e-06

#0 Foam::error:rintStack(Foam::Ostream&) in "/home/babak/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/home/babak/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/home/babak/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#4 Foam:perator/(Foam::UList<double> const&, Foam::UList<double> const&) in "/home/babak/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#5 Foam::diagonalSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/home/babak/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#6 Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/home/babak/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#7 Foam::fvMatrix<double>::solve() in "/home/babak/OpenFOAM/babak-2.1.x/platforms/linux64GccDPOpt/bin/turbulentFlameletRhoPimpleFoam"
#8
in "/home/babak/OpenFOAM/babak-2.1.x/platforms/linux64GccDPOpt/bin/turbulentFlameletRhoPimpleFoam"
#9 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#10
in "/home/babak/OpenFOAM/babak-2.1.x/platforms/linux64GccDPOpt/bin/turbulentFlameletRhoPimpleFoam"
Floating point exception
rm: cannot remove `pro*': No such file or directory


I should note that , I am trying to run your SANDIA_3D case in LES.
Have you ever experienced this error?

Regards
Bobi
babakflame is offline   Reply With Quote

Old   February 22, 2013, 12:16
Default
  #133
Super Moderator
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Tussenhausen
Posts: 2,708
Blog Entries: 6
Rep Power: 51
Tobi has a spectacular aura aboutTobi has a spectacular aura aboutTobi has a spectacular aura about
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
Hi Bobi,

hmmmm ... I got a similar error for my transient flamelet solver. Till now I am not sure how to solve that but I am not working on that solver anymore.

Did you implement your time-derivation in all equations?

\frac{\partial Z'}{\partial t}
.
.
.

Further more you have to refine the mesh (but I think you know that)

Error #3 seems to be a division with zero so you get the floating point exeption (but I am not sure)

I think gschaider or bruno knows more about the numbered outputs.
Tobi is offline   Reply With Quote

Old   February 22, 2013, 14:13
Default
  #134
Senior Member
 
Bobby
Join Date: Oct 2012
Location: Michigan
Posts: 454
Rep Power: 15
babakflame is on a distinguished road
Hi Tobi
I will contact Bernhard and Bruno soon. However, I have a conceptual problem with the issue. Why should I need the time derivatives of mixture fraction and it's variance?
I have seen lots of papers, like Prof. Kempf's ones that he has combined steady flamelet (Equations without time-derivatives) with LES (that needs transient solvers like Pimple).

If we import time-derivatives, then it will be unsteady flamelet. Although with unsteady flamelet we have more accuracy due to widespread libraries, But I think that we can still combine steady flamelet with LES solvers.

The problem is in the fvSolution file Tobi, Cause I made some changes in that file and the error changed. I think I should coincide solvers for pimple algorithm. I will try my best to solve it.

Regards
Bobi

Last edited by babakflame; February 22, 2013 at 15:23.
babakflame is offline   Reply With Quote

Old   February 25, 2013, 04:25
Default
  #135
Super Moderator
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Tussenhausen
Posts: 2,708
Blog Entries: 6
Rep Power: 51
Tobi has a spectacular aura aboutTobi has a spectacular aura aboutTobi has a spectacular aura about
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
Quote:
Originally Posted by babakflame View Post
Hi Tobi
I will contact Bernhard and Bruno soon. However, I have a conceptual problem with the issue. Why should I need the time derivatives of mixture fraction and it's variance?
I have seen lots of papers, like Prof. Kempf's ones that he has combined steady flamelet (Equations without time-derivatives) with LES (that needs transient solvers like Pimple).

If we import time-derivatives, then it will be unsteady flamelet. Although with unsteady flamelet we have more accuracy due to widespread libraries, But I think that we can still combine steady flamelet with LES solvers.

The problem is in the fvSolution file Tobi, Cause I made some changes in that file and the error changed. I think I should coincide solvers for pimple algorithm. I will try my best to solve it.

Regards
Bobi
I think that there are differences between | calculating flamelets steady state and solving the flow structur.

If you wanna make steady state flamelets (like Albertos tool does) you have to neglect the time derivations. BUT in the flow structure you should use the time derivation (in my opinion).
Tobi is offline   Reply With Quote

Old   March 6, 2013, 13:30
Default
  #136
New Member
 
Join Date: Feb 2013
Posts: 7
Rep Power: 13
danz is on a distinguished road
Hi Tobi,

I'm doing my master thesis in the field of flamelet models. The problem is that I already read so many papers and I do not know how to generate the flamelet libraries (computational technic). Can you help me by giving to me some basic ideas or plan to do that?

Best
danz is offline   Reply With Quote

Old   March 6, 2013, 15:31
Default
  #137
Super Moderator
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Tussenhausen
Posts: 2,708
Blog Entries: 6
Rep Power: 51
Tobi has a spectacular aura aboutTobi has a spectacular aura aboutTobi has a spectacular aura about
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
Hi Danz,

you mean the calculation of the flamelet libs ?

Do you want to build your flamelets by your own or programming your own flamelet calculator?

To generate the flamelets you first have to solve the flamelet quations!
After that you have all physical and chemical variable as a function of Z

\phi = f(Z)

You can use that libraries only for diffusion flames couse the flamelet model is build on the diffusion flames. After that and if you want to use the flamelet model for turbulent RANS eqn. you have to expand your flamelets to the variance.

For that you have to use propability distribution functions (PDF) depend on the variables. For chi you have to use the log-normal function, for the enthalpy defect you should use the dirac delta function and for the rest you have to use the beta-pdf.

Then you have every value phi depend on the following:

\tilde \phi = f (\tilde Z, \tilde Z'', \tilde \chi, \tilde \phi_H)


In my master thesis I wanted to calculate flamelets by my own but there was no time to build a programm

The reason was due to alberto cuocis flamelet calculator. There are just two binarys to execute and I had no ability to get through the code

But I hope it helps you a bit.
IF you build a own flamelet - calculater it would be very nice to get involved there.


If you are searching for a Flamelet calculator you can use Pitsch´s calculator. There you can have a look into the whole c++ code but the code is annoying


Tobi
Tobi is offline   Reply With Quote

Old   March 7, 2013, 04:40
Default
  #138
New Member
 
Join Date: Feb 2013
Posts: 7
Rep Power: 13
danz is on a distinguished road
Hi,

Thank you for your reply. Actually I just start my master thesis but the goal is to generate my own flamelet libraries
So that why I am asking to you how did you do that ...
danz is offline   Reply With Quote

Old   March 7, 2013, 05:13
Default
  #139
Super Moderator
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Tussenhausen
Posts: 2,708
Blog Entries: 6
Rep Power: 51
Tobi has a spectacular aura aboutTobi has a spectacular aura aboutTobi has a spectacular aura about
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
Hi,

i build my flamelets with the binarys included in the libOpenSmoke libraries
Tobi is offline   Reply With Quote

Old   March 12, 2013, 14:23
Default
  #140
Super Moderator
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Tussenhausen
Posts: 2,708
Blog Entries: 6
Rep Power: 51
Tobi has a spectacular aura aboutTobi has a spectacular aura aboutTobi has a spectacular aura about
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
Hi danz,

are you working on a creator for flamelets now ?

Would you like to use c++ ? I know a few good libraries for that.

If you wanna do that I want to be involved if its okay ?

After that it would be a greate expansion to the flamelet tool by using a compilable flamelet generator which everyone can modify.
Tobi 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
Numerical treatment of the source term in combustion equations Tobi Main CFD Forum 37 September 15, 2020 14:42
[openSmoke] flameletSmoke + new ODESolver (by Alberto Cuoci) Tobi OpenFOAM Community Contributions 1 November 21, 2017 19:24
Unsteady solver with Flamelet Model (libOpenSMOKE) francesco_capuano OpenFOAM Running, Solving & CFD 11 November 26, 2013 05:50
LibOpenSmoke, getting the species in ParaFoam Christoph_84 OpenFOAM 1 May 31, 2012 15:42


All times are GMT -4. The time now is 18:29.