{"id":4911,"date":"2020-05-11T15:33:25","date_gmt":"2020-05-11T04:33:25","guid":{"rendered":"https:\/\/www.computationalfluiddynamics.com.au\/?p=3214"},"modified":"2025-02-05T10:17:46","modified_gmt":"2025-02-04T23:17:46","slug":"y-plus_part1_understanding-the-physics-of-boundary-layers","status":"publish","type":"post","link":"https:\/\/www.leapaust.com.au\/blog\/cfd\/y-plus_part1_understanding-the-physics-of-boundary-layers\/","title":{"rendered":"What y+ should I use? Part 1 &#8211; Understanding the physics of boundary layers"},"content":{"rendered":"<div id=\"bsf_rt_marker\"><\/div>\t\t<div data-elementor-type=\"wp-post\" data-elementor-id=\"4911\" class=\"elementor elementor-4911\" data-elementor-post-type=\"post\">\n\t\t\t\t<div class=\"elementor-element elementor-element-57715724 e-flex e-con-boxed e-con e-parent\" data-id=\"57715724\" 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-992c130 elementor-widget elementor-widget-text-editor\" data-id=\"992c130\" 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>LEAP Australia has been working with university departments and engineering teams providing Ansys simulation software and technical support. The largest student engineering competition in Australia is Formula Student, where students learn how to design, build and race Formula 1 (F1) like cars around Winton Racetrack in Victoria. We actively sponsor university teams around Australia and New Zealand, providing licenses and technical support, giving us a unique insight into common questions posed by students.<\/p><p>Formula Student teams designing aerodynamics kits should understand the underlying physics, mathematical models and Computational Fluid Dynamics (CFD) concepts, which is difficult due to limited available public research. This blog post focuses on a common question: What y+ should I use in my simulations? To answer this question, the blog is broken into 3 parts in this series:<\/p><ul><li><strong>Part 1 \u2013 Understanding the physics of boundary layer<\/strong><strong>s<\/strong><\/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><h1>Part 1 \u2013 The underlying physics<\/h1><h3>What is the boundary layer?<\/h3><p>The boundary layer refers to the fluid region near a wall where the viscous forces are significant, and for Formula Student, understanding how it behaves is necessary to determine how much force can be produced by a wing. An important aspect of the boundary layer is how resolving it helps predict flow separation. Flow separation on a front wing changes the amount of force it creates, and the flow field behind it, affecting cooling, the diffuser and rear wing. By understanding and controlling the boundary layer, engineers can maximise the effect of the front wing.<\/p><p>With a better understanding of the boundary layer, F1 engineers can control the flow from the front wing. Check out F1 engineer and Formula Student Aerodynamics judge\u00a0<span style=\"text-decoration: underline;\"><a href=\"https:\/\/www.linkedin.com\/in\/willemtoet1\/\" target=\"_blank\" rel=\"noopener noreferrer\">Willem Toet<\/a><\/span> as he explains how the <a href=\"https:\/\/www.youtube.com\/watch?v=xjsUwFHQbWo\" target=\"_blank\" rel=\"noopener noreferrer\">front wing manages airflow around an F1 car in this video<\/a>.<\/p><p>Now that we understand the importance of the boundary layer, let us look at its structure in detail. In the boundary layer, the velocity profile changes rapidly and can be broken up into three regions: viscous sublayer, log layer and outer layer, as shown in Figure 1 below.<\/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-496b20b elementor-widget elementor-widget-image\" data-id=\"496b20b\" data-element_type=\"widget\" data-e-type=\"widget\" data-widget_type=\"image.default\">\n\t\t\t\t<div class=\"elementor-widget-container\">\n\t\t\t\t\t\t\t\t\t\t\t\t\t\t\t<img fetchpriority=\"high\" decoding=\"async\" width=\"778\" height=\"456\" src=\"https:\/\/www.leapaust.com.au\/blog\/wp-content\/uploads\/2020\/05\/y-image-1-1.png\" class=\"attachment-full size-full wp-image-5000\" alt=\"\" srcset=\"https:\/\/www.leapaust.com.au\/blog\/wp-content\/uploads\/2020\/05\/y-image-1-1.png 778w, https:\/\/www.leapaust.com.au\/blog\/wp-content\/uploads\/2020\/05\/y-image-1-1-300x176.png 300w, https:\/\/www.leapaust.com.au\/blog\/wp-content\/uploads\/2020\/05\/y-image-1-1-768x450.png 768w\" sizes=\"(max-width: 778px) 100vw, 778px\" \/>\t\t\t\t\t\t\t\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-375c28b elementor-widget elementor-widget-text-editor\" data-id=\"375c28b\" 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><em>Figure 1 \u2013 Schematic of the velocity boundary layer (Source: Ansys Introduction training material)<\/em><\/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-c95edbb elementor-widget elementor-widget-text-editor\" data-id=\"c95edbb\" 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>Resolving the boundary layer close to the wall provides an accurate representation of its profile, leading to accurate predictions of wall shear stress, surface pressure, the effect of adverse pressure gradients and forces.<\/p><p>In FSAE, due to the shape of the wings, there are large pressure gradients that must be resolved to accurately predict downforce and drag.<\/p><p>For more information on boundary layers check out these videos:\u00a0<span style=\"text-decoration: underline;\"><a href=\"https:\/\/www.youtube.com\/playlist?list=PLwjkNXfc1xQYVdk4G0iYIhGVrVWQOnKPY\" data-wplink-edit=\"true\">Click here<\/a><\/span><\/p><h3>How does CFD resolve the boundary layer?<\/h3><p>By using a log-scale and empirical equations to non-dimensionalise the velocity (u<sup>+<\/sup>) and the distance from the wall (y<sup>+<\/sup>), the velocity profile of the boundary layers takes on a predictable form, as shown below by the red line in Figure 2. There are two other equations, shown via blue lines, that have been derived capturing the profile of the viscous sublayer and the log-region. For more information about the mathematical equations check out here: <span style=\"text-decoration: underline;\"><a href=\"http:\/\/brennen.caltech.edu\/fluidbook\/basicfluiddynamics\/turbulence\/lawofthewall.pdf\" target=\"_blank\" rel=\"noopener noreferrer\">Click Me <\/a><\/span>\u00a0(in this article u* and y* can be treated the same as u<sup>+<\/sup> and y<sup>+<\/sup>).<\/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-cbc856f elementor-widget elementor-widget-image\" data-id=\"cbc856f\" data-element_type=\"widget\" data-e-type=\"widget\" data-widget_type=\"image.default\">\n\t\t\t\t<div class=\"elementor-widget-container\">\n\t\t\t\t\t\t\t\t\t\t\t\t\t\t\t<img decoding=\"async\" width=\"831\" height=\"783\" src=\"https:\/\/www.leapaust.com.au\/blog\/wp-content\/uploads\/2020\/05\/y-image-2-1.png\" class=\"attachment-full size-full wp-image-5001\" alt=\"\" srcset=\"https:\/\/www.leapaust.com.au\/blog\/wp-content\/uploads\/2020\/05\/y-image-2-1.png 831w, https:\/\/www.leapaust.com.au\/blog\/wp-content\/uploads\/2020\/05\/y-image-2-1-300x283.png 300w, https:\/\/www.leapaust.com.au\/blog\/wp-content\/uploads\/2020\/05\/y-image-2-1-768x724.png 768w\" sizes=\"(max-width: 831px) 100vw, 831px\" \/>\t\t\t\t\t\t\t\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-60fe6fd elementor-widget elementor-widget-text-editor\" data-id=\"60fe6fd\" 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><em>Figure 2 \u2013 Non-dimensional velocity versus distance from the wall (Source: <a href=\"https:\/\/en.wikipedia.org\/wiki\/Law_of_the_wall\" target=\"_blank\" rel=\"noopener noreferrer\"><span style=\"text-decoration: underline;\">https:\/\/en.wikipedia.org\/wiki\/Law_of_the_wall<\/span><\/a>)<\/em><\/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-e500593 elementor-widget elementor-widget-text-editor\" data-id=\"e500593\" 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>In CFD software the boundary layer is modelled by extruding the surface mesh elements at the wall. The surface mesh can be made up of different element types, but for the purposes of this discussion, we will consider triangular surface elements which are extruded into prism cells. The boundary layer resolution depends on the first cell height and the number of prism layers. The solver calculates a value, y<sup>+<\/sup>, which is computed from the first cell height, local flow speed and fluid properties. For the same fluid and flow conditions, we can say that reducing the first prism height by a factor of 2 will reduce y<sup>+<\/sup>\u00a0by a factor of 2. Obtaining an initial value for y<sup>+ <\/sup>and then adjusting the prism heights in the meshing tool or via adaption is a common practice in CFD when tuning the mesh for best results.<\/p><p>Different turbulence models are used to capture the boundary layer in different ways. The k-\u03c9 SST model has been widely researched, tested and is accepted as a good choice for many industrial problems, and has the robustness to resolve the viscous sublayer or the log layer. If the first cell height is in the viscous layer, it will resolve the boundary layer by evaluating the velocity at each cell. If the first cell is in the log layer, the software will make an empirical assumption about the viscous sublayer, referred to as using a <em>wall function<\/em>.<\/p><p>A careful selection of how you would like to model the boundary layer is essential.<\/p><p>When resolving the viscous sublayer, a value for y<sup>+<\/sup>~1 is required, which will require significantly more mesh elements and effort, Figure 3, left-hand image. This will lead to higher meshing and solution times and the requirement for more computer system resources. If you determine, that is not important to resolve the viscous sublayer, i.e. to use a wall function, a value for y<sup>+<\/sup> &gt; 30 is valid; this will result in less meshing effort and computational time, Figure 3, right-hand image. Note, however, that using wall functions can result in modelling errors, the significance of which will vary depending on flow conditions.<\/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-7c418f8 elementor-widget elementor-widget-image\" data-id=\"7c418f8\" data-element_type=\"widget\" data-e-type=\"widget\" data-widget_type=\"image.default\">\n\t\t\t\t<div class=\"elementor-widget-container\">\n\t\t\t\t\t\t\t\t\t\t\t\t\t\t\t<img decoding=\"async\" width=\"867\" height=\"351\" src=\"https:\/\/www.leapaust.com.au\/blog\/wp-content\/uploads\/2020\/05\/y-image-3-1.png\" class=\"attachment-full size-full wp-image-5002\" alt=\"\" srcset=\"https:\/\/www.leapaust.com.au\/blog\/wp-content\/uploads\/2020\/05\/y-image-3-1.png 867w, https:\/\/www.leapaust.com.au\/blog\/wp-content\/uploads\/2020\/05\/y-image-3-1-300x121.png 300w, https:\/\/www.leapaust.com.au\/blog\/wp-content\/uploads\/2020\/05\/y-image-3-1-768x311.png 768w\" sizes=\"(max-width: 867px) 100vw, 867px\" \/>\t\t\t\t\t\t\t\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-03c3b86 elementor-widget elementor-widget-text-editor\" data-id=\"03c3b86\" 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><em>Figure 3 \u2013 Left \u2013 Near wall mesh when resolving the viscous sublayer, Right \u2013 Near wall mesh when using wall functions<\/em><br \/><em>(Source: Ansys Introduction training material)<\/em><\/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-34d39761 elementor-widget elementor-widget-text-editor\" data-id=\"34d39761\" 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>It is recommended to either resolve the viscous sublayer by creating a mesh with y<sup>+<\/sup>~1\u00a0or to use the wall function approach by having the first element in the log region (y<sup>+ <\/sup>&gt;30) and avoid having a y<sup>+ <\/sup>between those values. The models in Ansys CFD can manage the entire regime; however, it is good practice to make the intention clear and use one approach or the other.<\/p><p>The truth is most Formula Student teams do not have the computational resources to run a y<sup>+<\/sup>~1\u00a0for every 3D simulation of every design across the different yaw, pitch and roll conditions they desire.<\/p><p>With a better understanding of the boundary layer, how it is modelled in CFD and what y<sup>+<\/sup>\u00a0is, we can now address the question \u201cWhich one should we pick, and how many inflation layers are needed?\u201d<\/p><p>Stay tuned\u00a0for the next in this 3-part blog series:<\/p><ul><li><strong>Part 2 \u2013 <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\">Resolving each region of the boundary layer<\/a><\/span><\/strong><\/li><li><strong>Part 3 \u2013 <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\">Understanding impact of Y+ and number of prism layers on flow resolution<\/a><\/span><\/strong><\/li><\/ul>\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>This blog series focuses on a common question: What y+ should I use in my CFD simulations? This first post is designed to help you understand the physics of boundary layers in relation to CFD meshes and Y+ values.<\/p>\n","protected":false},"author":3,"featured_media":4518,"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,174,181,395,396,446,449,273,508,515],"class_list":["post-4911","post","type-post","status-publish","format-standard","has-post-thumbnail","hentry","category-cfd","tag-aerodynamics","tag-ansys-cfd","tag-ansys-meshing","tag-cfd-modelling-of-turbulent-flows","tag-cfd-simulation-accuracy","tag-inflation-layer-meshing","tag-mesh-control-settings","tag-prediction-of-aerodynamic-flows","tag-turbulence-modelling","tag-wall-functions"],"_links":{"self":[{"href":"https:\/\/www.leapaust.com.au\/blog\/wp-json\/wp\/v2\/posts\/4911","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=4911"}],"version-history":[{"count":6,"href":"https:\/\/www.leapaust.com.au\/blog\/wp-json\/wp\/v2\/posts\/4911\/revisions"}],"predecessor-version":[{"id":5008,"href":"https:\/\/www.leapaust.com.au\/blog\/wp-json\/wp\/v2\/posts\/4911\/revisions\/5008"}],"wp:featuredmedia":[{"embeddable":true,"href":"https:\/\/www.leapaust.com.au\/blog\/wp-json\/wp\/v2\/media\/4518"}],"wp:attachment":[{"href":"https:\/\/www.leapaust.com.au\/blog\/wp-json\/wp\/v2\/media?parent=4911"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/www.leapaust.com.au\/blog\/wp-json\/wp\/v2\/categories?post=4911"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/www.leapaust.com.au\/blog\/wp-json\/wp\/v2\/tags?post=4911"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}