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

Compressible WALE model

Register Blogs Community New Posts Updated Threads Search

Like Tree14Likes

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   February 3, 2012, 13:39
Default Compressible WALE model
  #1
Senior Member
 
n/a
Join Date: Sep 2009
Posts: 199
Rep Power: 16
deji is on a distinguished road
Hello to all here at the forum. I am trying to implement the compressible version of the WALE model into OpenFOAM, but thus far haven't found any papers on it. Does anyone know where I can find literature?

Cheers,
Deji
deji is offline   Reply With Quote

Old   August 1, 2012, 14:09
Default
  #2
Senior Member
 
Awais Ali
Join Date: Feb 2010
Location: Germany
Posts: 128
Rep Power: 17
owayz is on a distinguished road
Send a message via MSN to owayz
Hallo Deji,
Were you successful in implementing the WALE model for compressible flow in openFOAM I am also interested in the same topic would be interesting if you could give me some pointers about where did you start from and how much have you done till now?
Regards,
Awais
owayz is offline   Reply With Quote

Old   August 1, 2012, 14:12
Default
  #3
Senior Member
 
Awais Ali
Join Date: Feb 2010
Location: Germany
Posts: 128
Rep Power: 17
owayz is on a distinguished road
Send a message via MSN to owayz
And as far as the paper is concerned I think this one is the paper which describes the details. I just matched the author names and the year published and found this paper via google.
http://www.math.univ-montp2.fr/~nicoud/PDF/wale.pdf
owayz is offline   Reply With Quote

Old   August 1, 2012, 14:13
Default
  #4
Senior Member
 
n/a
Join Date: Sep 2009
Posts: 199
Rep Power: 16
deji is on a distinguished road
Hello. There isn't a compressible version of the WALE model, so I just coded the incompressible version into my compressible LES solver. It does work okay, I can send you the SGS solver if you would like.

Cheers,
Deji
mksingh, owayz and jiujiumin like this.
deji is offline   Reply With Quote

Old   August 1, 2012, 14:14
Default
  #5
Senior Member
 
n/a
Join Date: Sep 2009
Posts: 199
Rep Power: 16
deji is on a distinguished road
I also coded the locally dynamic Smagorinsky compressible SGS model and I prefer it over the WALE model.
owayz and jiujiumin like this.
deji is offline   Reply With Quote

Old   August 1, 2012, 15:27
Default
  #6
Senior Member
 
Awais Ali
Join Date: Feb 2010
Location: Germany
Posts: 128
Rep Power: 17
owayz is on a distinguished road
Send a message via MSN to owayz
I would be very happy if you could send me both the models.
I would like to see how they have been implemented and if I run them I will share the results with you.
I am going to send you my email via message.
Or you can post them here, this way others will also benefit from your implementation of both the models.
Regards,
Awais
owayz is offline   Reply With Quote

Old   August 1, 2012, 15:45
Default
  #7
Senior Member
 
n/a
Join Date: Sep 2009
Posts: 199
Rep Power: 16
deji is on a distinguished road
Sure thing. I will send you both SGS solvers via email and you can post it on the forum if you like.

Cheers,
Deji
deji is offline   Reply With Quote

Old   August 1, 2012, 16:49
Default
  #8
New Member
 
M K Singh
Join Date: Sep 2009
Posts: 19
Rep Power: 16
mksingh is on a distinguished road
Quote:
Originally Posted by deji View Post
Sure thing. I will send you both SGS solvers via email and you can post it on the forum if you like.

Cheers,
Deji
Hi,
Could you please also send me both SGS solvers to me?
Thanks in advance.
M K Singh
mksingh is offline   Reply With Quote

Old   August 2, 2012, 09:54
Default
  #9
Senior Member
 
n/a
Join Date: Sep 2009
Posts: 199
Rep Power: 16
deji is on a distinguished road
Hey guys. My apologies for not sending the solvers yesterday, I am trying to clean up the code with comments to make it more understandable. Hence, I will attach it here on the forum this evening at the very latest. My apologies once again and I look forward to you all giving me feedback if there are any problems or bugs in the codes.

Cheers,
Deji
owayz likes this.
deji is offline   Reply With Quote

Old   August 2, 2012, 21:32
Default
  #10
Senior Member
 
n/a
Join Date: Sep 2009
Posts: 199
Rep Power: 16
deji is on a distinguished road
Hey guys, I have attached the WALE model for the compressible SGS models in OpenFOAM. I have also attached the LESProperties file which gives the needed input for the LES model.

As for the locally dynamic Smagorinsky compressible LES model, I will upload it shortly. I need to make a quick fix to the code.

Please let me know if there are any problems.

Cheers,
Deji
Attached Files
File Type: c oneKEqEddyWale.C (4.9 KB, 330 views)
File Type: h oneKEqEddyWale.H (3.7 KB, 220 views)
File Type: txt LESProperties.txt (2.2 KB, 279 views)
vonboett likes this.
deji is offline   Reply With Quote

Old   August 2, 2012, 21:44
Default
  #11
Senior Member
 
n/a
Join Date: Sep 2009
Posts: 199
Rep Power: 16
deji is on a distinguished road
Here is the locally dynamic Smagorinksy compressible LES SGS model. The LESProperties that I attached previously also has the necessary coefficients to run this SGS model. I ought to say that I have had problems running dynamic models in OpenFOAM at t=0. Therefore, run any static LES model such as Smagorinksy or the WALE model for one timestep, then you can proceed to utilize the dynamic model thereafter.

Cheers,
Deji
Attached Files
File Type: c dynLocalAverageSmagorinsky.C (6.0 KB, 210 views)
File Type: h dynLocalAverageSmagorinsky.H (4.1 KB, 180 views)
rkc.cfd likes this.
deji is offline   Reply With Quote

Old   August 3, 2012, 12:33
Default
  #12
Senior Member
 
Awais Ali
Join Date: Feb 2010
Location: Germany
Posts: 128
Rep Power: 17
owayz is on a distinguished road
Send a message via MSN to owayz
I need sometime I will surely get back to you on this topic.
Thanks for sharing them.
Regards,
Awais
owayz is offline   Reply With Quote

Old   August 6, 2012, 11:47
Default
  #13
New Member
 
M K Singh
Join Date: Sep 2009
Posts: 19
Rep Power: 16
mksingh is on a distinguished road
Dear Deji,
Thanks for your solvers. I will test it and update you in this regard.
M K Singh
mksingh is offline   Reply With Quote

Old   August 6, 2012, 19:50
Default
  #14
Senior Member
 
Awais Ali
Join Date: Feb 2010
Location: Germany
Posts: 128
Rep Power: 17
owayz is on a distinguished road
Send a message via MSN to owayz
Hi deji,
I have never compiled a turbulence model before in openfoam. Can you tell me how should I compile your your code?
Regards,
Awais
owayz is offline   Reply With Quote

Old   August 7, 2012, 13:20
Default
  #15
Senior Member
 
n/a
Join Date: Sep 2009
Posts: 199
Rep Power: 16
deji is on a distinguished road
Hi deji,
I have never compiled a turbulence model before in openfoam. Can you tell me how should I compile your your code?
Regards,
Awais


The LES models I uploaded are for compressible LES solvers in OpenFOAM. You need to copy the code into the directory where the other SGS models are, go through the files in the MAKE directory and add the lines necessary for the compiler to know that it should compile the codes too. Then just type wmake libso, it is quite simple, but you should spend some time going through the directory that has the SGS models.

Cheers,
Deji
deji is offline   Reply With Quote

Old   August 20, 2012, 04:47
Default
  #16
Member
 
cosimo bianchini
Join Date: Mar 2009
Location: Florence, Tuscany, Italy
Posts: 88
Rep Power: 17
cosimobianchini is on a distinguished road
Send a message via Skype™ to cosimobianchini
I also worked in the past a bit on the compressible WALE sgs model.
My solution was indeed closer to original paper notation and slightly more compact, see it below and check if the two versions are equivalent (too much tensorial algebra for me ).

volSymmTensorField S = symm(gradU);
volTensorField Omega = skew(gradU);

volTensorField Sd = (S & S) + (Omega & Omega) - 1.0/3.0*I*( (S && S) - (Omega && Omega));

muSgs_ = rho()*sqr(C_omega_*delta())*pow((Sd && Sd),3.0/2.0)/(pow((S && S),5.0/2.0)+pow((Sd && Sd),5.0/4.0));

Concerning your version I do not understand why you need the computation of a transport equation for the subgrid turbulent kinetic energy. The model is a simple algebraic one so it is not clear what is the purpose of the extra transported quantity.
Regards,
Cosimo
PonchO likes this.
__________________
Cosimo Bianchini

Ergon Research s.r.l.
Via Panciatichi, 92
50127 Florence - ITALY
Tel: +39 055 0763716
Mob: +39 320 9460153
e-mail: cosimo.bianchini@ergonresearch.it
URL: www.ergonresearch.it
cosimobianchini is offline   Reply With Quote

Old   October 1, 2012, 15:06
Default
  #17
Senior Member
 
n/a
Join Date: Sep 2009
Posts: 199
Rep Power: 16
deji is on a distinguished road
Concerning your version I do not understand why you need the computation of a transport equation for the subgrid turbulent kinetic energy. The model is a simple algebraic one so it is not clear what is the purpose of the extra transported quantity.

Hey Cosimo. The only reason why the k transport equation is there is for turbulent combustion problems; such problems tend to require a turbulent k field to solve for models such as the eddy dissipation combustion model.

Sorry for the late response.

Cheers,
Deji
deji is offline   Reply With Quote

Old   January 7, 2013, 11:26
Default
  #18
Senior Member
 
Albrecht vBoetticher
Join Date: Aug 2010
Location: Zürich, Swizerland
Posts: 237
Rep Power: 16
vonboett is on a distinguished road
Thanks Dejio for your work, you are really filling the gap.
Preferring the local dyn. Smagorinsky I think you might be interested in the WALE-Similarity Mixed model, if you are not tired of it jet: http://link.springer.com/chapter/10....LI=true#page-1
vonboett is offline   Reply With Quote

Old   January 16, 2014, 14:39
Default
  #19
New Member
 
Sagnik
Join Date: Oct 2012
Posts: 28
Rep Power: 13
sagnikmazumdar is on a distinguished road
Quote:
Originally Posted by deji View Post
I also coded the locally dynamic Smagorinsky compressible SGS model and I prefer it over the WALE model.
Hi Deji, I was going through your CFD online posts and got to know that you have a compressible version of 'dynLocalAverageSmagorinsky' model. Would it be possible to share it with us ! This model does not exist in the new version of OpenFOAM but we found that its incompressible version does a better job than other existing model for some cases. We would be happy to acknowledge you as one of the contributors in our work.

I can be reached at 'sagnikmazumdar@gmail.com'.

Thanks for all the help.

Sagnik
jiujiumin likes this.
sagnikmazumdar is offline   Reply With Quote

Old   January 17, 2014, 12:26
Default
  #20
New Member
 
Sagnik
Join Date: Oct 2012
Posts: 28
Rep Power: 13
sagnikmazumdar is on a distinguished road
Hi Deji, I got the files. I did not see it was already attached to another post. Thanks for the help.

Sagnik
sagnikmazumdar is offline   Reply With Quote

Reply


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
LES Compressible Smagorinsky Model iyer_arvind OpenFOAM Running, Solving & CFD 26 September 9, 2014 07:22
Low-Re turbulence model for compressible flow volker OpenFOAM Programming & Development 0 March 15, 2010 10:20
LES and combustion model Margherita Cadorin CFX 0 October 29, 2008 05:24
Problem calculate Y in compressible turbulence model luca OpenFOAM Running, Solving & CFD 5 June 1, 2006 05:53
simulation of RSM in compressible mixing layer lars Main CFD Forum 0 February 19, 2004 22:40


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