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 Tree127Likes

Reply
 
LinkBack Thread Tools Display Modes
Old   July 6, 2013, 08:49
Default
  #201
Senior Member
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Leoben (Austria)
Posts: 1,089
Blog Entries: 6
Rep Power: 19
Tobi will become famous soon enough
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
Hi bobi and all,


I compiled OF-2.2.x on my laptop today and tried the flamelet model 2.2.x:
-------------------------------------------------------------------------------------------------------------

1. git clone ... is working fine (no errors)
2. flamelet libs compiled well
3. thermodynamics do not compile

That problem is solved by removing a deb file

4. flameletSimpleFoam is not compiling

That problem is going to solved. The files are "crashed" there are some nonsence lines in the files. I remove them now and clean the files

5. flameletPisoFoam is compiling

6. tutorial SIMPLE is not working due to changes in the PDF library

That problem is going to be solved

7. tutorial PISO is not working due to changes in the PDF library

That problem is going to be solved



Tobi
babakflame likes this.
Tobi is offline   Reply With Quote

Old   July 6, 2013, 09:18
Default flameletModel 2.2.x
  #202
Senior Member
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Leoben (Austria)
Posts: 1,089
Blog Entries: 6
Rep Power: 19
Tobi will become famous soon enough
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
Flamelet Model 2.2.x is now cleared and working without any problems.


Please test it!
Tobi is offline   Reply With Quote

Old   July 6, 2013, 11:09
Default
  #203
Senior Member
 
Bobi
Join Date: Oct 2012
Posts: 286
Rep Power: 5
babakflame is on a distinguished road
Many Thanks to my friend Tobi
babakflame is offline   Reply With Quote

Old   July 24, 2013, 01:31
Default
  #204
Senior Member
 
Bobi
Join Date: Oct 2012
Posts: 286
Rep Power: 5
babakflame is on a distinguished road
Hi Tobi

What is the convergence criterion in your simple and piso codes and how can we increase or decrease the convergence accuracy?
I mean where should we implement the modifications?

Regards]
Bobi
babakflame is offline   Reply With Quote

Old   July 29, 2013, 06:31
Default convergence
  #205
Senior Member
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Leoben (Austria)
Posts: 1,089
Blog Entries: 6
Rep Power: 19
Tobi will become famous soon enough
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
Hi Bobi,

as you can see in the sandia-flame the residuals never go below 1e-3 - 1e-4. In that case I wait till i reach a "steady-state" residual state; like a horizontal line.

But that depend on your case. In my masterthesis it was possible to get a residual to 1e-5 but for that I decreased the relax factor strong. Further more I always save the last time steps (like every 100 iterations) and compare the differences between the slices/Clips/values to be sure that I reached a steady state solution.

Problem in most cases I had is the fact, that the combustion Chambers (in my case) never be a steady. I always realize transient movements in my values on certain Points in the geomtry.



For the transient (piso) algo I use the time, residuals and visualisation.

In the sandia-flame I get a better convergence in my Simulation but that was all I simulated with my transient solver.
babakflame likes this.
Tobi is offline   Reply With Quote

Old   August 13, 2013, 04:33
Default About simpleFoamflamelet solver
  #206
Member
 
sonu
Join Date: Jul 2013
Location: delhi
Posts: 79
Rep Power: 4
yash.aesi is on a distinguished road
helo ,
After successfully running the tutorial of Sandia COH2N2 flame in flamelet model under the simplefoamflamelet solver now i want to try it for my geometry and case . i had converted my mesh file into OF , so now can anybody tell me which file i have to copy from the tutorial and kept in my working folder to solve my case ?
yash.aesi is offline   Reply With Quote

Old   August 13, 2013, 11:34
Default
  #207
Senior Member
 
Bobi
Join Date: Oct 2012
Posts: 286
Rep Power: 5
babakflame is on a distinguished road
Hi

You should make your own PDF file, 0 (initial conditions) and constant folder
for required coefficients. Just the same as Tobi tutorial.

Bobi
babakflame is offline   Reply With Quote

Old   August 14, 2013, 03:48
Default
  #208
Member
 
sonu
Join Date: Jul 2013
Location: delhi
Posts: 79
Rep Power: 4
yash.aesi is on a distinguished road
thanks alot bobi .....i will try as u said ....
yash.aesi is offline   Reply With Quote

Old   August 21, 2013, 09:04
Default Mixing of two gases with turbulentFlameletRhoSimplecFoam
  #209
TBO
Member
 
Join Date: May 2013
Location: Netherlands
Posts: 30
Rep Power: 4
TBO is on a distinguished road
The turbulentFlameletRhoSimplecFoam creates a lot of opportunities for me for steady state combustion modelling in OpenFOAM. In the end, I am very interested in simulations with chemical reactions and combustion involved and that is what this solver is doing. Also thank you for the good documentation about the chemistry involved in this solver, that is what is missing for some other chemistry solvers on CFD online.

However, I am also interested in steady state modelling of the mixing between two different gases (fuel and air). So a cold flow case without the actual chemical reaction which shows the mass or volume fractions of the different species. Is it possible to do this with turbulentFlameletRhoSimplecFoam and what do I need to change for this (e.g. in the library). I am using OpenFOAM 2.1.1.
TBO is offline   Reply With Quote

Old   August 21, 2013, 11:23
Default
  #210
Senior Member
 
Bobi
Join Date: Oct 2012
Posts: 286
Rep Power: 5
babakflame is on a distinguished road
Hi TBO

My suggestion for you for simulating cold flows is using RhoPimpleFoam and RhoSimpleFoam solvers of OpenFoam, cause their flow solvers are exactly the same as libOpenSmoke modified versions. LibOpenSmoke and its modified bversions are steady flamelet solvers which automatically spread combustion according to mixture fraction value in the domain.

Good Luck

Bobi
babakflame is offline   Reply With Quote

Old   August 22, 2013, 02:36
Default
  #211
TBO
Member
 
Join Date: May 2013
Location: Netherlands
Posts: 30
Rep Power: 4
TBO is on a distinguished road
Quote:
Originally Posted by babakflame View Post
Hi TBO

My suggestion for you for simulating cold flows is using RhoPimpleFoam and RhoSimpleFoam solvers of OpenFoam, cause their flow solvers are exactly the same as libOpenSmoke modified versions. LibOpenSmoke and its modified bversions are steady flamelet solvers which automatically spread combustion according to mixture fraction value in the domain.

Good Luck

Bobi

Hi Bobi,

Thank you for your reply, please correct me if I'm wrong, since I'm not very experienced with OpenFOAM so far. I am already using rhoSimplecFoam for cold flow modelling, however what I would like to do is calculate the mixture fractions of different species in the domain. So for example the mixing of methane jets in an air flow. As far as I know rhoSimplecFoam can only cope with one species (if anyone can tell me how to implement several species in rhoSimplecFoam, I am very satisfied).

I already looked into other solvers, e.g. reactingFoam and there you have the possibility to switch off chemical reactions (however this solver is transient, resulting in extreme calculation times). Is there any possibility to switch off the chemical reactions (or modify the library) for the libOpenSmoke?
immortality likes this.
TBO is offline   Reply With Quote

Old   August 22, 2013, 05:29
Default
  #212
Senior Member
 
Bobi
Join Date: Oct 2012
Posts: 286
Rep Power: 5
babakflame is on a distinguished road
Hi TBO
HTML Code:
if anyone can tell me how to implement several species in rhoSimplecFoam, I am very satisfied).
I did not get what you mean by this sentence. However, If you mean various inlet constituents with not having combustion , you might use fireFoam cause this solver supports this option.

HTML Code:
Is there any possibility to switch off the chemical reactions (or modify the library) for the libOpenSmoke?
You can insert this option for libOpenSmoke if you have sufficient knowledge of OpenFOAM and combustion, However my proposal for you is using fireFOAM in O.F. 2.2.x , cause in this version this option (the same as you want ) is added to the current combustion models. In addition, it calculates the mixture fraction of the domain without having combustion though.

I hope you find this helpful

Bobi
babakflame is offline   Reply With Quote

Old   August 22, 2013, 05:46
Default
  #213
TBO
Member
 
Join Date: May 2013
Location: Netherlands
Posts: 30
Rep Power: 4
TBO is on a distinguished road
Hi Bobi,

Thank you for your answer. Indeed I mean various inlet constituents without having combustion.

I already looked into fireFoam (as well as reactingFoam) before, however the main problem with these solvers is that they are transient and I would like to model steady state (since this saves a lot of computational time).

I already looked into several other options (alternateSteadyReactingFoam, edcSimpleFoam). However results for these solvers are not satisfying jet and documentation/CFD online discussions are very limited.

Regards,
TBO is offline   Reply With Quote

Old   August 22, 2013, 11:15
Default
  #214
Senior Member
 
Bobi
Join Date: Oct 2012
Posts: 286
Rep Power: 5
babakflame is on a distinguished road
Hi TBO

HTML Code:
I already looked into fireFoam (as well as reactingFoam)
I should mention that these solvers are totally different and if you use them you can find out that why I proposed you firefoam with no combustion model (so much faster than reactingfoam with off chemistry). Although you are right. Both of them use piso algorithm which is a transient flow solver.

HTML Code:
I would like to model steady state
I am still sure that you can use RhoSimpleFoam for your aim

HTML Code:
documentation/CFD online discussions are very limited
Ask your question in forum (not this one), You will get your answer

Bobi
babakflame is offline   Reply With Quote

Old   September 16, 2013, 05:41
Default
  #215
New Member
 
parsa
Join Date: Sep 2013
Posts: 3
Rep Power: 3
parsa_abbasi is on a distinguished road
Dear All,

I'm doing also the same job with OpenFoam, but until now, i didn't get the semialer result as I got in Transient mode. I made several steady-State Solvers, but seems that No Combustion happens! I'm using the 'psi Combustion Model' which is used in reactingFoam Solver as well. If any of you could run the Transient Solver in Steady-State mode ( with delta T=1 for example , and not with deltaT= 10^-6) or already created a new steady-state solver , I would be so grateful if help me this way,

looking forward to your replies,

Sincerely,

Parsa
parsa_abbasi is offline   Reply With Quote

Old   September 16, 2013, 12:03
Default
  #216
Senior Member
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Leoben (Austria)
Posts: 1,089
Blog Entries: 6
Rep Power: 19
Tobi will become famous soon enough
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
Dear all,

this thread is for libOpenSMOKE and not for other combustion solver.

To have a clear thread, please make new threads with your problems

Thanks
Tobi
babakflame and TBO like this.
Tobi is offline   Reply With Quote

Old   September 16, 2013, 14:28
Default
  #217
Senior Member
 
Bobi
Join Date: Oct 2012
Posts: 286
Rep Power: 5
babakflame is on a distinguished road
Hi Tobi

I have an issue with k parameter in simulating bluff-body flames. what should I take for this value? For instance; In your Sandia_CO/H/N2 you have mentioned as:

Code:
  outlet
    {
        type            zeroGradient;
    }
    axis
    {
        type            empty;
    }
    sidewall
    {
        type            compressible::kqRWallFunction;
        value           uniform 1;
    }
    burnerwall
    {
        type            compressible::kqRWallFunction;
        value           uniform 1;
    }
for the boundary walls. Also in another tutorial for a bluff-body flame I found:

Code:
bluffbody
    {
        type            compressible::kqRWallFunction;
        value           uniform 17.34;
    }
    back
    {
        type            wedge;
    }

    outlet
    {
        type            inletOutlet;
        inletValue      uniform 1;
        value           uniform 1;
    }
Would you PLZ hint me about it that how these parameters are taken?

My specific issue is its value for the bluff-body boundary as a kind of wall.


It is the initial value, So is it important choosing the exact right value. I suggest maybe just estimating a value and starting my simulation.

Am I right Buddy?

Regards

Bobi

Last edited by babakflame; September 17, 2013 at 04:22.
babakflame is offline   Reply With Quote

Old   September 17, 2013, 06:10
Default
  #218
Senior Member
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Leoben (Austria)
Posts: 1,089
Blog Entries: 6
Rep Power: 19
Tobi will become famous soon enough
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
Quote:
Originally Posted by parsa_abbasi View Post
Dear All,

I'm doing also the same job with OpenFoam, but until now, i didn't get the semialer result as I got in Transient mode. I made several steady-State Solvers, but seems that No Combustion happens! I'm using the 'psi Combustion Model' which is used in reactingFoam Solver as well. If any of you could run the Transient Solver in Steady-State mode ( with delta T=1 for example , and not with deltaT= 10^-6) or already created a new steady-state solver , I would be so grateful if help me this way,

looking forward to your replies,

Sincerely,

Parsa

Hi Parsa,

just one question. Are you using the libOpenSMOKE tool and the flamelet solver or a other solver?

It s confusing cause you mentioned you are using the ''psi Combustion Model" ?

Tobi
Tobi is offline   Reply With Quote

Old   September 17, 2013, 06:16
Default
  #219
Senior Member
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Leoben (Austria)
Posts: 1,089
Blog Entries: 6
Rep Power: 19
Tobi 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

My specific issue is its value for the bluff-body boundary as a kind of wall.


It is the initial value, So is it important choosing the exact right value. I suggest maybe just estimating a value and starting my simulation.

Am I right Buddy?

Regards

Bobi
Hi Bobi,

you are right. The values you choose are just initial values. After the first iteration the value change. But if you choose your value very good your simulation get to a steady state solution faster.

Additionally the first timesteps could be hard for openfoam to calculate due to bounding things and stuff like that. So its good to have a accurate value.

In my tutorial case its just a standard value if I do not know a correct one
I always begin with 0.01 to 1 ... due to the velocity and geometry.


The inlet value can be calculated (not exact but good for initialising).
For the wall I could use 10 or 100 too

So just start with 0.1 or 1 and check the first 10 or 50 timesteps / iterations ... if you have no problems with the wall value then everything is okay.

I read in an other thread that People set the k and Epsilon values to 1e5 or sth. like that

Well in my opinion you should set a correct dimension.

So if you know that the inlet has k=0.0032 then I choose the same Dimension for the walls like 0.001 ... but not 10


Hope this answer will help you for your Problem.

Tobi
immortality and babakflame like this.
Tobi is offline   Reply With Quote

Old   September 18, 2013, 14:11
Default
  #220
Member
 
sonu
Join Date: Jul 2013
Location: delhi
Posts: 79
Rep Power: 4
yash.aesi is on a distinguished road
greetings oll ,

i have a small doubt that in 0 folder the value of each species (like CH4) in boundary fields is in mole fraction . rite ?? not the mass fraction

correct me if i am wrong




Thank you ,
Regards
sonu
yash.aesi 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 03:39.