{"id":4867,"date":"2013-04-12T18:29:44","date_gmt":"2013-04-12T07:29:44","guid":{"rendered":"https:\/\/www.computationalfluiddynamics.com.au\/?p=932"},"modified":"2025-02-07T11:01:16","modified_gmt":"2025-02-07T00:01:16","slug":"turbulence-part-3-selection-of-wall-functions-and-y-to-best-capture-the-turbulent-boundary-layer","status":"publish","type":"post","link":"https:\/\/www.leapaust.com.au\/blog\/cfd\/turbulence-part-3-selection-of-wall-functions-and-y-to-best-capture-the-turbulent-boundary-layer\/","title":{"rendered":"Turbulence Part 3 &#8211; Selection of wall functions and Y+ to best capture the Turbulent Boundary Layer"},"content":{"rendered":"<div id=\"bsf_rt_marker\"><\/div>\t\t<div data-elementor-type=\"wp-post\" data-elementor-id=\"4867\" class=\"elementor elementor-4867\" data-elementor-post-type=\"post\">\n\t\t\t\t<div class=\"elementor-element elementor-element-36cee4b0 e-flex e-con-boxed e-con e-parent\" data-id=\"36cee4b0\" data-element_type=\"container\" data-e-type=\"container\">\n\t\t\t\t\t<div class=\"e-con-inner\">\n\t\t\t\t<div class=\"elementor-element elementor-element-845c367 elementor-widget elementor-widget-text-editor\" data-id=\"845c367\" data-element_type=\"widget\" data-e-type=\"widget\" data-widget_type=\"text-editor.default\">\n\t\t\t\t<div class=\"elementor-widget-container\">\n\t\t\t\t\t\t\t\t\t<p>Note: this is an old post.\u00a0 The updated post series from 2020 is LEAP&#8217;s 3-Part Series on &#8220;What y+ should I use in my simulations?&#8221; which is available here:<\/p><ul><li><span style=\"text-decoration: underline;\"><a href=\"https:\/\/www.leapaust.com.au\/blog\/cfd\/y-plus_part1_understanding-the-physics-of-boundary-layers\/\" target=\"_blank\" rel=\"noopener noreferrer\"><strong>Part 1 \u2013 Understanding the physics of boundary layer<\/strong><strong>s<\/strong><\/a><\/span><\/li><li><span style=\"text-decoration: underline;\"><a href=\"https:\/\/www.leapaust.com.au\/blog\/cfd\/y-plus_part2_resolving-each-region-of-the-boundary-layer\/\" target=\"_blank\" rel=\"noopener noreferrer\"><strong>Part 2 \u2013 Resolving each region of the boundary layer<\/strong><\/a><\/span><\/li><li><span style=\"text-decoration: underline;\"><a href=\"https:\/\/www.leapaust.com.au\/blog\/cfd\/y-plus_part3_understanding-impact-of-y-and-number-of-prism-layers-on-flow-resolution\/\" target=\"_blank\" rel=\"noopener noreferrer\"><strong>Part 3 \u2013 Understanding impact of Y+ and number of prism layers on flow resolution<\/strong><\/a><\/span><\/li><\/ul><p>Old Post continues:<\/p><p>In recent posts in our <a href=\"https:\/\/www.computationalfluiddynamics.com.au\/tag\/turbulence-modelling\/\">series of Turbulence Modelling<\/a> posts, we have covered boundary layer theory and touched on some useful meshing and post-processing guidelines\u00a0to check you are appropriately resolving the boundary layer profile.\u00a0 Today we will consider three critical questions that are often asked by CFD engineers when developing or refining a CFD simulation:<\/p><p><em>&#8211; Am I using the correct turbulence model for the type of results I am looking for?<\/em><\/p><p><em>&#8211; Do I have an appropriate Y+ value and a sufficient number of inflation layers?<\/em><\/p><p><em>&#8211; Am I using the right wall function for my problem?<\/em><\/p><p>This topic is so important because we know that in turbulent flows the\u00a0velocity fluctuations within\u00a0the turbulent boundary layer can be a significant percentage of the mean flow velocity, so it is\u00a0critical that we capture these effects with accuracy. A Reynolds averaging approach using turbulence models will provides us with an estimate of the increased levels of stress within the boundary layer, termed the Reynolds stresses.\u00a0In order to appreciate the use of wall functions and the\u00a0influence of walls on\u00a0the turbulent flowfield, we should first gain familiarity with the composite regions of the turbulent boundary layer:<\/p>\t\t\t\t\t\t\t\t<\/div>\n\t\t\t\t<\/div>\n\t\t\t\t<div class=\"elementor-element elementor-element-409cc350 elementor-widget elementor-widget-text-editor\" data-id=\"409cc350\" data-element_type=\"widget\" data-e-type=\"widget\" data-widget_type=\"text-editor.default\">\n\t\t\t\t<div class=\"elementor-widget-container\">\n\t\t\t\t\t\t\t\t\t<p>\u00a0<\/p><figure id=\"attachment_1009\" aria-describedby=\"caption-attachment-1009\" style=\"width: 300px\" class=\"wp-caption alignleft\"><a href=\"https:\/\/www.leapaust.com.au\/blog\/wp-content\/uploads\/2013\/04\/wallfunction.png\"><img fetchpriority=\"high\" decoding=\"async\" class=\"size-medium wp-image-1009\" title=\"wallfunction\" src=\"https:\/\/www.leapaust.com.au\/blog\/wp-content\/uploads\/2013\/04\/wallfunction-300x195.png\" alt=\"\" width=\"300\" height=\"195\" \/><\/a><figcaption id=\"caption-attachment-1009\" class=\"wp-caption-text\">Composite regions of the turbulent boundary layer<\/figcaption><\/figure><p>\u00a0<\/p><p>\u00a0<\/p><p>\u00a0<\/p><p>\u00a0<\/p><p>\u00a0<\/p><p>\u00a0<\/p><p>\u00a0<\/p><p>\u00a0<\/p><p>\u00a0<\/p><p>In the laminar sub-layer region (Y+ &lt; 5) inertial forces are less domineering and the flow\u00a0exhibits laminar characteristics, which is why\u00a0this is known as the low-Re region. Low-Re turbulent models (e.g. the SST model) aim to resolve this area\u00a0and therefore require an appropriate mesh resolution to do this with accuracy. This is most critical for flows with a changing pressure gradient where we expect to see separation, as observed below.<\/p><figure id=\"attachment_937\" aria-describedby=\"caption-attachment-937\" style=\"width: 300px\" class=\"wp-caption alignleft\"><a href=\"https:\/\/www.leapaust.com.au\/blog\/wp-content\/uploads\/2013\/03\/separation.png\"><img decoding=\"async\" class=\"size-medium wp-image-937\" title=\"separation\" src=\"https:\/\/www.leapaust.com.au\/blog\/wp-content\/uploads\/2013\/03\/separation-300x275.png\" alt=\"\" width=\"300\" height=\"275\" \/><\/a><figcaption id=\"caption-attachment-937\" class=\"wp-caption-text\">Predicting separation in a diffuser-type geometry<\/figcaption><\/figure><p>\u00a0<\/p><p>\u00a0<\/p><p>\u00a0<\/p><p>\u00a0<\/p><p>\u00a0<\/p><p>\u00a0<\/p><p>\u00a0<\/p><p>\u00a0<\/p><p>\u00a0<\/p><p>\u00a0<\/p><p>\u00a0<\/p><p>In the law of the wall region, inertial forces strongly dominate over viscous forces and we have a high presence of turbulent stresses (this is known as the high-Re composite region).\u00a0If using a low-Re model, the whole turbulent boundary layer will be resolved including the log-law region. However, it possible to use semi-empirical expressions known as <strong>wall functions<\/strong> to bridge the viscosity-affected region between the wall and the fully-turbulent region.<\/p><figure id=\"attachment_941\" aria-describedby=\"caption-attachment-941\" style=\"width: 300px\" class=\"wp-caption alignleft\"><a href=\"https:\/\/www.leapaust.com.au\/blog\/wp-content\/uploads\/2013\/03\/ssteddyvisc.png\"><img decoding=\"async\" class=\"size-medium wp-image-941\" title=\"ssteddyvisc\" src=\"https:\/\/www.leapaust.com.au\/blog\/wp-content\/uploads\/2013\/03\/ssteddyvisc-300x213.png\" alt=\"\" width=\"300\" height=\"213\" srcset=\"https:\/\/www.leapaust.com.au\/blog\/wp-content\/uploads\/2013\/03\/ssteddyvisc-300x213.png 300w, https:\/\/www.leapaust.com.au\/blog\/wp-content\/uploads\/2013\/03\/ssteddyvisc-1024x728.png 1024w, https:\/\/www.leapaust.com.au\/blog\/wp-content\/uploads\/2013\/03\/ssteddyvisc-768x546.png 768w, https:\/\/www.leapaust.com.au\/blog\/wp-content\/uploads\/2013\/03\/ssteddyvisc.png 1268w\" sizes=\"(max-width: 300px) 100vw, 300px\" \/><\/a><figcaption id=\"caption-attachment-941\" class=\"wp-caption-text\">Contours of the eddy viscosity ratio on a low-Re grid illustrating high turbulent viscosity in the log-law region as opposed to the laminar sub-layer<\/figcaption><\/figure><p>\u00a0<\/p><p>\u00a0<\/p><p>\u00a0<\/p><p>\u00a0<\/p><p>\u00a0<\/p><p>\u00a0<\/p><p>\u00a0<\/p><p>\u00a0<\/p><p>\u00a0<\/p><p>\u00a0<\/p><p>The main benefit of\u00a0this wall function approach lies in the significant reduction in mesh resolution and thus reduction in simulation time. However, the shortcoming lies in numerical results deteriorating under subsequent refinement of the grid in wall normal direction (thus reducing the Y+ value\u00a0into the buffer layer zone). Continued reduction of Y+ to below 15\u00a0can gradually result in unbounded errors in wall shear stress and wall heat transfer (due to the damping functions inherent within the wall function approach).<\/p><p>Bearing all of the above in mind, and keeping our eye on finding the right balance between accuracy, stability and speed, we can tackle a wide variety of CFD problems using the following guidelines:<\/p><p><strong>What results am I interested in<\/strong>\u00a0<strong>and am I using the right turbulence model?<\/strong><\/p><p>If our aim is to accurately predict the boundary layer velocity or thermal profile, or if the developing boundary layer will tend to separate (due to a changing pressure gradient \u2013 and not because of sharp edges or discontinuities in the geometry), then we recommend the use of a low-Re model. Low-Re models are also required for accurate pressure-drop or drag calculations. We highly recommend the use of the Shear Stress Transport (SST) model, but all <em>\u03c9\u00ad<\/em>-based models or <em>\u03b5<\/em>-based models with enhanced wall treatment may be used. For high speed external aerodynamic flows, the one-equation Spalart-Allmaras model (with Y+ &lt; 2)\u00a0may also be considered to reduce the computational time. Alternatively, for flows where wall-bounded effects are not a priority, or if separation is expected to occur only due to sharp changes in the geometry, an <em>\u03b5<\/em>-based wall function approach is more than sufficient. In ANSYS CFD, all <em>\u03c9<\/em>-based models and the SST model are capable of resorting to a wall function formulation (automatic wall treatment) in the presence of coarse mesh resolutions near the wall without any further user input. Wall function models are also useful for calibrating our CFD models, due to the decreased simulation time.<\/p><figure id=\"attachment_944\" aria-describedby=\"caption-attachment-944\" style=\"width: 300px\" class=\"wp-caption alignleft\"><a href=\"https:\/\/www.leapaust.com.au\/blog\/wp-content\/uploads\/2013\/03\/image2.png\"><img loading=\"lazy\" decoding=\"async\" class=\"size-medium wp-image-944\" title=\"image2\" src=\"https:\/\/www.leapaust.com.au\/blog\/wp-content\/uploads\/2013\/03\/image2-300x248.png\" alt=\"\" width=\"300\" height=\"248\" \/><\/a><figcaption id=\"caption-attachment-944\" class=\"wp-caption-text\">Flow pattern with separation and reattachment on a rotor blade<\/figcaption><\/figure><p>\u00a0<\/p><p>\u00a0<\/p><p>\u00a0<\/p><p>\u00a0<\/p><p>\u00a0<\/p><p>\u00a0<\/p><p>\u00a0<\/p><p>\u00a0<\/p><p>\u00a0<\/p><p>\u00a0<\/p><p><strong>What is my Y+ value and do I have a sufficient number of prism layers?<\/strong><\/p><p>When using low-Re models or any models with enhanced wall treatment, the average Y+ value should be on the order of ~1 to ensure we are capturing the laminar sub-layer. When using wall function models, the Y+ value should ideally be above 15 to avoid erroneous modelling in the buffer layer and the laminar sub-layer. High quality numerical results for the boundary layer will only be obtained if the overall resolution of the boundary layer is sufficient. This requirement in some cases is more important than achieving certain Y+ values. The minimum number of cells to cover a boundary layer accurately is around 10, but values of 20 are desirable. The total thickness of the prisms should be implemented such that around 15 or more nodes are actually covering the boundary layer. Our <a title=\"Turbulence Part 4 \u2013 Reviewing how well you have resolved the Boundary Layer\" href=\"https:\/\/www.computationalfluiddynamics.com.au\/tips-tricks-turbulence-part-4-reviewing-how-well-you-have-resolved-the-boundary-layer\/\">next post<\/a> in this series on the turbulent boundary layer will cover a very useful and practical technique to post-process the resolution of the boundary layer, and offer insight into modifications required to improve accuracy.<\/p><figure id=\"attachment_945\" aria-describedby=\"caption-attachment-945\" style=\"width: 300px\" class=\"wp-caption alignleft\"><a href=\"https:\/\/www.leapaust.com.au\/blog\/wp-content\/uploads\/2013\/03\/image3.png\"><img loading=\"lazy\" decoding=\"async\" class=\"size-medium wp-image-945\" title=\"image3\" src=\"https:\/\/www.leapaust.com.au\/blog\/wp-content\/uploads\/2013\/03\/image3-300x189.png\" alt=\"\" width=\"300\" height=\"189\" \/><\/a><figcaption id=\"caption-attachment-945\" class=\"wp-caption-text\">Boundary layer velocity profile modeled with standard k-e for three different mesh densities using Enhanced Wall Treatment<\/figcaption><\/figure><p>\u00a0<\/p><p>\u00a0<\/p><p>\u00a0<\/p><p>\u00a0<\/p><p>\u00a0<\/p><p>\u00a0<\/p><p>\u00a0<\/p><p>\u00a0<\/p><p>\u00a0<\/p><p><strong>Am I using the right wall function?<\/strong><\/p><p>In ANSYS CFD, all turbulence models are y-plus independent. However selecting the <em>most appropriate<\/em> wall function is dependent on level of refinement of our wall adjacent mesh, or the relative scales in our flow. Use of the standard wall function (<em>\u03b5<\/em>-based models) implies that our boundary layer mesh lies entirely within the log-law region of the boundary layer. For industrial applications, this\u00a0in fact\u00a0might be difficult to achieve due to varying geometrical and velocity scales associated with our model \u2013 and therefore grids inherently designed with arbitrary refinement. We highly recommend the use of the scalable wall function, which offers an elegant solution to this ambiguity often encountered. This wall function virtually displaces the mesh to a Y+ ~ 11.225 (transition to the log-law composite layer) irrespective of the level of refinement, thereby avoiding the erroneous modelling of the laminar sub-layer and buffer region. It is also important to note that for grids designed with a Y+ &gt; 11.225, the scalable wall function will provide identical results to the standard wall function. Enhanced wall treatment may further be selected for <em>\u03b5<\/em>-based models on refined low-Re grids, and is also formulated such that it can perform well for meshes of intermediate resolution. However the use of enhanced (or non-equilibrium) wall treatment for low-Re modelling of the turbulent boundary layer is generally not recommended and more confidence in our solution can be obtained by selecting a suitable <em>\u03c9\u00ad<\/em>-based formulation, such as the SST model.[\/vc_column_text][\/vc_column][\/vc_row]<\/p>\t\t\t\t\t\t\t\t<\/div>\n\t\t\t\t<\/div>\n\t\t\t\t\t<\/div>\n\t\t\t\t<\/div>\n\t\t\t\t<\/div>\n\t\t","protected":false},"excerpt":{"rendered":"<p>Note: this is an old post.\u00a0 The updated post series from 2020 is LEAP&#8217;s 3-Part Series on &#8220;What y+ should I use in my simulations?&#8221; which is available here: Part 1 \u2013 Understanding the physics of boundary layers Part 2 \u2013 Resolving each region of the boundary layer Part 3 \u2013 Understanding impact of Y+&hellip;&nbsp;<a href=\"https:\/\/www.leapaust.com.au\/blog\/cfd\/turbulence-part-3-selection-of-wall-functions-and-y-to-best-capture-the-turbulent-boundary-layer\/\" rel=\"bookmark\">Read More &raquo;<span class=\"screen-reader-text\">Turbulence Part 3 &#8211; Selection of wall functions and Y+ to best capture the Turbulent Boundary Layer<\/span><\/a><\/p>\n","protected":false},"author":3,"featured_media":4149,"comment_status":"open","ping_status":"closed","sticky":false,"template":"","format":"standard","meta":{"neve_meta_sidebar":"","neve_meta_container":"","neve_meta_enable_content_width":"","neve_meta_content_width":0,"neve_meta_title_alignment":"","neve_meta_author_avatar":"","neve_post_elements_order":"","neve_meta_disable_header":"","neve_meta_disable_footer":"","neve_meta_disable_title":"","neve_meta_reading_time":"","footnotes":""},"categories":[323],"tags":[169,393,395,396,446,273,487,508,515],"class_list":["post-4867","post","type-post","status-publish","format-standard","has-post-thumbnail","hentry","category-cfd","tag-aerodynamics","tag-cfd-meshing-methods","tag-cfd-modelling-of-turbulent-flows","tag-cfd-simulation-accuracy","tag-inflation-layer-meshing","tag-prediction-of-aerodynamic-flows","tag-shear-stress-transport-sst-turbulence-model","tag-turbulence-modelling","tag-wall-functions"],"_links":{"self":[{"href":"https:\/\/www.leapaust.com.au\/blog\/wp-json\/wp\/v2\/posts\/4867","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/www.leapaust.com.au\/blog\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/www.leapaust.com.au\/blog\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/www.leapaust.com.au\/blog\/wp-json\/wp\/v2\/users\/3"}],"replies":[{"embeddable":true,"href":"https:\/\/www.leapaust.com.au\/blog\/wp-json\/wp\/v2\/comments?post=4867"}],"version-history":[{"count":3,"href":"https:\/\/www.leapaust.com.au\/blog\/wp-json\/wp\/v2\/posts\/4867\/revisions"}],"predecessor-version":[{"id":5047,"href":"https:\/\/www.leapaust.com.au\/blog\/wp-json\/wp\/v2\/posts\/4867\/revisions\/5047"}],"wp:featuredmedia":[{"embeddable":true,"href":"https:\/\/www.leapaust.com.au\/blog\/wp-json\/wp\/v2\/media\/4149"}],"wp:attachment":[{"href":"https:\/\/www.leapaust.com.au\/blog\/wp-json\/wp\/v2\/media?parent=4867"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/www.leapaust.com.au\/blog\/wp-json\/wp\/v2\/categories?post=4867"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/www.leapaust.com.au\/blog\/wp-json\/wp\/v2\/tags?post=4867"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}