[CFD] The Spalart-Allmaras Turbulence Model

Ғылым және технология

A brief introduction to the Spalart-Allmaras turbulence model. The following topics are covered:
1) 3:04 Why was the Spalart-Allmaras Turbulence Model Proposed?
2) 10:03 What do each of the terms in the model mean?
3) 18:49 What boundary conditions should be used with the model?
#spalartAllmaras #fluidmechanics101 #turbulence
==================================
Some useful references:
1) G. Kalitzin, G. Medic, G. Iaccarino and P. Durbin
Near-wall behaviour of RANS turbulene models and implications for wall functions
Journal of Computational Physics 204 (2005) 265-291
2) Spalart, P. R. and Allmaras, S. R.
A One-Equation Turbulence Model for Aerodynamic Flows
AIAA Paper 92-0439 (1992)
www.researchgate.net/publicat...
============================================
Want to learn more?
============================================
Grab a copy of my CFD Fundamentals Course (for beginners):
www.udemy.com/course/computat...
Learn how to write your own CFD code in MATLAB and python (for intermediates):
dr-aidan-wimshurst-s-school.t...
Learn how I draw my figures and diagrams in Inkscape (for everyone):
dr-aidan-wimshurst-s-school.t...
============================================
Did you like the video?
============================================
Download the lecture slides from my website:
www.fluidmechanics101.com/pag...
Buy me a coffee to say thanks:
www.buymeacoffee.com/NKEZrhvg0
Support the channel on Patreon (and get useful extras for your CFD studies):
/ fluidmechanics101
============================================
Donations:
============================================
1) PayPal
www.paypal.me/fluidmechanics101
==================================
Disclaimer
==================================
The methods, algorithms, equations, formulae, diagrams and explanations in this talk are for educational and demonstrative purposes only. They should never be used to analyse, design, accredit or validate real scientific / engineering / mathematical structures and flow systems. For such applications, appropriate trained, qualified and accredited (SQEP) engineers / scientists should be consulted along with the appropriate documentation, procedures and engineering standards. Furthermore, the information contained within this talk has not been verified, peer reviewed or checked in any way and is likely to contain several errors. It is therefore not appropriate to use this talk itself (or any of the algorithms, equations, formulae, diagrams and explanations contained within this talk) as an academic or technical reference. The reader should consult the original references and follow the verification and validation processes adopted by your company / institution when carrying out engineering calculations and analyses. Fluid Mechanics 101 and Dr. Aidan Wimshurst are not accountable or liable in any form for the use or misuse of the information contained in this talk beyond the specific educational and demonstrative purposes for which it was intended.

Пікірлер: 69

  • @SyedHussain-wl6tc
    @SyedHussain-wl6tc2 жыл бұрын

    I cannot express how much of a service this is to the world, I'm doing a PhD and these videos are literally saving my ass. I come back and watch these videos over again every so often, every time I learn and understand a little more, every time I feel a sense of relief when the penny drops. Thank you so much.

  • @morr1177
    @morr11774 жыл бұрын

    I have been digging through papers for the last week or so and this explanation was very well laid out and could have saved me a lot of time on coming up to speed. Really really really well done!!! You are asset to the new folk in the CFD community, dont stop!!!

  • @fluidmechanics101

    @fluidmechanics101

    4 жыл бұрын

    Thanks Bryan 😊 much appreciated!

  • @shubhamsangodkar6992
    @shubhamsangodkar69924 жыл бұрын

    All, I can say is "I FOUND YOU, NuTILDA". This video has answered most of the questions with which i've been struggling far too long . Probably the most concise video on SA turbulence model on KZread. Couple of questions: 1) How would they have measured Turbulent Kinematic Viscosity experimentally , to obtain that profile? 2) When we speak about diffusion terms in the Momentum equation, What is getting diffused ? (Velocity?) How is diffusion related to viscosity? and how do we relate viscosity to concentration gradient intuitively ? 3) Can you differentiate how the SA Curvature Correction model is different from the Standard SA model? From , my experience the SA curvature correction model provides a much better correlation at higher AoA (i.e separation dominated flows), I would like to understand , why? Finally , thank you soo much for this one. P.S (Watching it on replay)

  • @noahzhang892
    @noahzhang8924 жыл бұрын

    I used to be super confused with the value nuTilda in OpenFOAM. Thank you so much for the explanation!

  • @mateuselias6018
    @mateuselias60184 жыл бұрын

    I have no words to describe how happy I am to have watched this video. Thank you a lot man. You could describe whats I have looking in just simple words that make me understand and increment my wish to keep going to the CFD world. Thank you a lot.

  • @professionalprocrastinator8103
    @professionalprocrastinator81034 жыл бұрын

    As usual, concise yet precise, crystal clear presentation! Great job

  • @roshniiyer3351
    @roshniiyer33514 жыл бұрын

    All that you have explained is quite lucid and informative. It gives a good overview to understand what actually happens when we use this model!

  • @jeremyyi7619
    @jeremyyi76193 жыл бұрын

    A FANTASTIC VIDEO! I finally start to understand some of the concepts. Thank you so much!

  • @itssachink
    @itssachink4 жыл бұрын

    thank you so much! you're explainations are the best. it is helping in understanding CFD a lot. never stop making videos. 👍👌😊

  • @AshokAshok-dt6lu
    @AshokAshok-dt6lu2 жыл бұрын

    Splendid as usual. Excellent explanation

  • @EclecticVibe
    @EclecticVibe4 жыл бұрын

    I was waiting for this! Thank you.

  • @fluidmechanics101

    @fluidmechanics101

    4 жыл бұрын

    Fantastic! I knew a few of you were looking for this video. It was a tricky one but well worth it 😊

  • @juanpablovaca-lago5659
    @juanpablovaca-lago56592 жыл бұрын

    Awsome video! Very clear and informative.

  • @lizarettflavour
    @lizarettflavour4 жыл бұрын

    Thank you! Very clear and informative =)

  • @jairoandres92
    @jairoandres924 жыл бұрын

    Thank you for the concise explanation. I am doing some work that requires me to understand this model and reading the original paper is sometimes. At least now I can read it with a clearer idea.

  • @saikrishna-kz4gs
    @saikrishna-kz4gs4 жыл бұрын

    Please do a video on LES/DES models

  • @emmanuelmaduakookorie8808
    @emmanuelmaduakookorie88083 жыл бұрын

    Wao! Great Job. Thanks for sharing.

  • @whalewolve
    @whalewolve5 ай бұрын

    Thank you for the very comprehensible explanation. One remark: comparing with my old cfd lecture manuscript, isnt there a c_mu missing in front of the right-hand expression in (25) ?

  • @shuangtu8067
    @shuangtu80676 ай бұрын

    I think 'S' in the production term in the S-A model is the magnitude of vorticity which is related to the rotation tensor instead of the shear-rate tensor. Correct me if I am wrong. Thanks

  • @florianschittl8079
    @florianschittl80794 жыл бұрын

    Fantastic! Can‘t wait to see more about turbulence modeling...maybe v2-f or reynolds-stress models? 😉👌

  • @fluidmechanics101

    @fluidmechanics101

    4 жыл бұрын

    🤔 hmmmm which one to do next ...

  • @florianschittl8079

    @florianschittl8079

    4 жыл бұрын

    Maybe the v2-f-model to complete the linear Eddy viscosity models 😉

  • @abdbot2509
    @abdbot2509 Жыл бұрын

    Good job bro congratulations

  • @Liontarix93
    @Liontarix932 жыл бұрын

    Thank you very much, that was very helpful, keep it up!

  • @rpapa
    @rpapa3 жыл бұрын

    These videos are great !!!!

  • @christopheryang6416
    @christopheryang64164 жыл бұрын

    Good lecture.

  • @DanielMolina-st9sp
    @DanielMolina-st9sp4 жыл бұрын

    just fantastic!!!!! easy way to learn

  • @user-oy3dk5jd5e
    @user-oy3dk5jd5e Жыл бұрын

    Thx

  • @javierlazaroortega8326
    @javierlazaroortega83264 ай бұрын

    For practical use, I understand SA is only applicable in exterior aerodynamics, and even there is less accurate that k-w SST. The only advantage it has is that it is computationally cheaper?

  • @5020alexi
    @5020alexi4 жыл бұрын

    Hey buddy...big fan of your videos...would you mind doing a video on the LES and DES models as it is commonly used in CFD simulations...would really appreciate it..:)

  • @fluidmechanics101

    @fluidmechanics101

    4 жыл бұрын

    I will start looking at LES soon. It is quite difficult, so i want to make sure i get it right 🙃

  • @nikhilnick8886
    @nikhilnick88864 жыл бұрын

    Thanks Dr.Aidan

  • @shangguicai211
    @shangguicai2114 жыл бұрын

    Nice episode. Could you please give a little more explanation about why nut is proportional to (y+)^4 in the viscous sublayer ?

  • @fluidmechanics101

    @fluidmechanics101

    4 жыл бұрын

    Im not too sure on this. I think it is from experimental measurements (DNS) of the viscous sub-layer. Let me get back to you on this one!

  • @macornman
    @macornman4 жыл бұрын

    Great talk, thanks so much. Can you make a video about the Rotational corrections applied to the SA turb model (SA-RC) to capture streamline curvature in RANS simulations?

  • @MrSrira
    @MrSrira4 жыл бұрын

    Hey...could u make a lecture video on Moving Reference Frame method!

  • @fluidmechanics101

    @fluidmechanics101

    4 жыл бұрын

    Oooh yes, thats a great idea! Will add it to my list!

  • @punkisinthedetails1470
    @punkisinthedetails14704 жыл бұрын

    Isn't it for external flow in the main aerofoils etc

  • @thomashamori3144
    @thomashamori3144 Жыл бұрын

    Great video, thanks! A question I had is: Why do you say that u_t = \tilde u f_{v1} has "quartic behavior"? It appears as (linear)*(cubic)/(cubic+C), which I am not sure how this is quartic?

  • @biaohuang7957
    @biaohuang79574 жыл бұрын

    Great video! By the way, is it possible for you to make a lecture video on dealing with the pressure boundaries (Prgh and/or P), in both incompressible and compressible solvers? It should be very helpful to understand how OpenFOAM handles the term before choosing the Boundary/Initial conditions provided by OpenFOAM, like prgh, p, prghTotalHydrostaticPressure, prghTotalPressure, prghPressure, and many more out there. I know there are many discussions on the subject on the cfd online or other forums, but a comprehensive introduction is still absent. Anyone interested in this, please let me know. Thanks in advance.

  • @fluidmechanics101

    @fluidmechanics101

    4 жыл бұрын

    Yep, i agree. Sounds like a good introductory video is needed!

  • @NITESHKSAHU
    @NITESHKSAHU Жыл бұрын

    Hi Aidan, Just wondering where are the turbulent kinetic energy and other 2/3 term in the Reynolds averaged equation? Since they were present at the begging of k-epsilon model video. Waiting for your reply.

  • @fluidmechanics101

    @fluidmechanics101

    Жыл бұрын

    As quite an old model, I think they did not bother with the additional terms (the effect of turbulence is only accounted through adding nut to the Laplacian terms). You are of course correct, and the full divergence of the Reynolds Stress tensor is actually needed. I am looking into this at the moment because it seems to be a problem with a lot of turbulence model implementations

  • @akbarravan5604

    @akbarravan5604

    20 күн бұрын

    @@fluidmechanics101 Hello Dr Wimshurst and thank you as always for this awesome content you created. I think @NITESHKSAHU meant the difference between Equation 1 in the present video and Equation 40 in kzread.info/dash/bejne/hYqNurCCc8i0h6Q.htmlsi=cg8wpmDK9IQqh5Qc In the Background slides there is no -2/3*rho*k*I anymore.

  • @biaohuang7957
    @biaohuang79574 жыл бұрын

    I was wondering that there might be a typo in the momentum equation. abla \dot \tau = abla \dot ((\mu+\mu_t)(( ablaU)+( ablaU)^T)) + abla \dot (-2/3\mu( abla \dot U)I) = abla \dot ((\mu+\mu_t)(( ablaU)+( ablaU)^T)) + abla(-2/3\mu( abla \dot U)). Correct me if I am wrong. Thanks.

  • @Ankit-df8rw
    @Ankit-df8rw4 жыл бұрын

    Ausome job

  • @pawankumar-ir1fd
    @pawankumar-ir1fd Жыл бұрын

    best

  • @jialiangzhou1428
    @jialiangzhou14283 жыл бұрын

    hi Dr.Aidan Thnak you for sharing this wonderful vedio. Could you please answer a stupid question from me. What happened if we apply Multiple reference Frame.?Which velocity are we going to use in Spalart-Allmaras transport equation? My opinion is to use a relative velocity in the LHS and all absoulte Velocity in the Rhs( production term,dissipate term and diffusion term).

  • @fluidmechanics101

    @fluidmechanics101

    3 жыл бұрын

    Good question! I think in the MRF approach we use absolute velocity for everything, so I would go with that here. Maybe have another watch of my 'MRF video' ??

  • @jialiangzhou1428

    @jialiangzhou1428

    3 жыл бұрын

    @@fluidmechanics101 would be any extra source term?

  • @varunrepo6913
    @varunrepo69134 жыл бұрын

    Thanks for the awesome videos! Are the slides that you use free available?

  • @fluidmechanics101

    @fluidmechanics101

    4 жыл бұрын

    I sell them on my website for super cheap (£10 for 10 sets of slides). Its just a little bit of money that helps me pay the webpage upkeep and buy research papers to help with the videos. Or you can get them all on Patreon if you prefer

  • @varunrepo6913

    @varunrepo6913

    4 жыл бұрын

    @@fluidmechanics101 Thanks

  • @marcorosatti4023
    @marcorosatti40234 жыл бұрын

    Hello! thank you for your explanation. :D I've tried to run SA model in an airfoil with C-domain (structured mesh) I got y^+ below 1 but I couldn't get good results in Cd coefficient. I've tried to solve a flow of Re = 7e5 (Transitional) in OpenFOAM... I think that this Reynolds number isn't compatible with OpenFOAM's SA model formulation. What I mean is the SA model hasn't a transitional model incorporated. In numbers I've could get results with a difference around 5% and 10 % in Cl and Cm coefficients in relation to experiment, but a difference of 60% in Cd. Could you give me your opinion of what happens whit this? Best regards!

  • @marcorosatti4023

    @marcorosatti4023

    4 жыл бұрын

    I found the answer. Effectively the OpenFOAM formulation hasn't the ft2 term, it has the same formulation that you explain in the video. I could code the standard version of SA model in OpenFOAM having good results. The formulation that you explain is only compatible with high Re cases!.

  • @stefano.lovato
    @stefano.lovato4 жыл бұрын

    Stupid question: In the momentum equations there should also be a term related to the Boussinesq approximation, namely -2/3 k delta_ij, right? If so, how can we find the turbulent kinetic energy k from nu_tilde? Could we also find the dissipation epsilon? Thank you for your help Aidan:)

  • @georgeelhaber3758

    @georgeelhaber3758

    3 жыл бұрын

    The turbulent kinetic energy in Boussinesq approximation of shear stresses is usually neglected when using the SA model ( the term is simply droped out). However, some variants of the model exist were the KE is approximated (rarely used). You can find a good summary of all the models and the contribution of each variant to the original model in the below link: turbmodels.larc.nasa.gov/spalart.html. Hope this helps :)

  • @fluidmechanics101

    @fluidmechanics101

    3 жыл бұрын

    Awesome! Thanks George

  • @gauravchhaperwal7519
    @gauravchhaperwal75194 жыл бұрын

    Hi, nice videos, solved many doubts about basics of CFD. I have one request/suggestion. You should also make videos on complete simulation of few standard problems. For example, I am working on vortex shedding about circular cylinder in a rectangular duct( 300 mm wide and 25 mm depth). I can found videos on the basics or the shallow simulations. But no videos are there that explain simulation of shedding phenomena using correct fundamentals of CFD and Fluid Mechanics. Please make a video on vortex shedding , Reynolds number 3000-15000 considering all basic details correctly. It will help a lot in practice/implementation of fundamentals.

  • @fluidmechanics101

    @fluidmechanics101

    4 жыл бұрын

    Yes, i am going to start doing some example problems on this channel. The difficulty is getting hold of some good meshing software, as CFD products and licenses tend to be very expensive! Also vortex shedding is one of the most difficult flow phenomena to simulate accurately using CFD, so dont worry if you are having difficulties. It is very hard!

  • @vivekchaurasia9807
    @vivekchaurasia98073 жыл бұрын

    Hi Aidan sir, I think equation 8 is not correct. should we not use del . (-2/3 meu.....)????? in equation 8

  • @fluidmechanics101

    @fluidmechanics101

    3 жыл бұрын

    Yep, I think there might be a typo here!

  • @punkisinthedetails1470
    @punkisinthedetails14704 жыл бұрын

    Is S-A model suitable for rotor stator interaction under MRF and under Sliding mesh.

  • @Ankit-df8rw

    @Ankit-df8rw

    4 жыл бұрын

    Maybe

  • @thortt515

    @thortt515

    4 жыл бұрын

    Yes, it is. At least it is so for most centrifual and axial compressors.

  • @punkisinthedetails1470

    @punkisinthedetails1470

    4 жыл бұрын

    @@thortt515 thanks a bunch man S-A runs fast too right? Like faster than Ke. What are the y+ needs? 5》30》200

  • @fluidmechanics101

    @fluidmechanics101

    4 жыл бұрын

    Yep, if you are going to run SA, try and get y+ = 1. Remember that aerofoils have adverse pressure gradients and the SA model was derived from flat plate boundary layer results with no pressure gradients. ... besides having y+ = 1 always gets reviewers/approvers/professors to leave you alone and stop being picky about the mesh ... 😏

  • @punkisinthedetails1470

    @punkisinthedetails1470

    4 жыл бұрын

    @@fluidmechanics101 Thanks !!!!!

Келесі