Compressible WALE model
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 |
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 |
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 |
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 |
I also coded the locally dynamic Smagorinsky compressible SGS model and I prefer it over the WALE model.:D
|
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 |
Sure thing. I will send you both SGS solvers via email and you can post it on the forum if you like.
Cheers, Deji |
Quote:
Could you please also send me both SGS solvers to me? Thanks in advance. M K Singh |
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 |
3 Attachment(s)
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 |
2 Attachment(s)
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 |
I need sometime I will surely get back to you on this topic. :)
Thanks for sharing them. Regards, Awais |
Dear Deji,
Thanks for your solvers. I will test it and update you in this regard. M K Singh |
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 |
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 |
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 |
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 |
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 |
Quote:
I can be reached at 'sagnikmazumdar@gmail.com'. Thanks for all the help. Sagnik |
Hi Deji, I got the files. I did not see it was already attached to another post. Thanks for the help.
Sagnik |
All times are GMT -4. The time now is 12:24. |