{"id":1323,"date":"2020-06-26T23:30:25","date_gmt":"2020-06-27T03:30:25","guid":{"rendered":"http:\/\/underthehood.blogwyrm.com\/?p=1323"},"modified":"2022-09-21T08:48:55","modified_gmt":"2022-09-21T12:48:55","slug":"elementary-viscous-flow-part-1","status":"publish","type":"post","link":"https:\/\/underthehood.blogwyrm.com\/?p=1323","title":{"rendered":"Elementary Viscous Flow – Part 1"},"content":{"rendered":"

This month\u2019s column deals with viscous flow in a variety of situations.\u00a0 Three main physical effects can drive elementary fluid flow: the pull of gravity coaxing a fluid downhill; a pressure differential, such as in a pipe, causing the fluid to flow from inlet to outlet; or shearing motion, which, in the usual scenario, involves a moving boundary dragging the fluid.\u00a0 We\u2019ll be looking at four examples covering different aspects of these three effects.\u00a0 The first two will illustrate the effect of gravitational body forces and pressure differentials to drive the fluid and will be covered in this column.\u00a0 The next two examples deal with steady and unsteady shear-induced flows and will be featured next month.<\/p>\n

But, before we deal with specific examples, we will need to settle on a model of the fluid.\u00a0 Usually, when seeking analytic solutions describing the motion of non-ideal Newtonian fluids (already a simplification compared to the amazingly rich and diverse types of fluids that nature offers), one either chooses to assume the flow is viscous and incompressible or that it is inviscid and compressible.\u00a0 For these examples, we will choose the former.<\/p>\n

The starting point will be the incompressible Navier-Stokes equations, of which the flux-conservative form is<\/p>\n

\\[ \\frac{\\partial \\rho {\\vec V}}{\\partial t} + \\nabla \\cdot (\\rho {\\vec V} {\\vec V}) = – \\nabla P + \\mu \\nabla^2 {\\vec V} + \\frac{1}{3} \\mu \\nabla \\left( \\nabla \\cdot {\\vec V} \\right) + \\rho {\\vec g} \\; .\\]<\/p>\n

Applying the incompressibility assumption ($$\\nabla \\cdot {\\vec V} = 0$$) is easy to do using the Lagrangian form of the equations, but it is more instructive to play with the flux-conservative form.\u00a0 Most of the action takes place on the left-hand side so let\u2019s focus on that by expanding the terms to get<\/p>\n

\\[ \\rho \\frac{\\partial {\\vec V}}{\\partial t} + {\\vec V} \\frac{\\partial \\rho}{\\partial t} + {\\vec V} ( {\\vec V} \\cdot \\nabla \\rho ) + \\rho ({\\vec V} \\cdot \\nabla) {\\vec V} + {\\vec V} \\rho (\\nabla \\cdot {\\vec V}) \\; .\\]<\/p>\n

The\u00a0last\u00a0term vanishes under the incompressibility assumption.\u00a0 The\u00a0second and third\u00a0terms, combined, are<\/p>\n

\\[ {\\vec V} \\left( \\frac{\\partial \\rho}{\\partial t} + ( {\\vec V} \\cdot \\nabla ) \\rho \\right) \\; .\\]<\/p>\n

The quantity in the parentheses vanishes when applying incompressibility to the continuity equation, as follows:<\/p>\n

\\[ \\frac{\\partial \\rho}{\\partial t} + ({\\vec V} \\cdot \\nabla ) \\rho + \\rho \\nabla \\cdot {\\vec V} = \\frac{\\partial \\rho}{\\partial t} + ({\\vec V} \\cdot \\nabla ) \\rho = 0 \\; ,\\]<\/p>\n

where the first expression is the full continuity equation and the second is what is left after the incompressibility assumption is applied.\u00a0 The left-hand side is now simply<\/p>\n

\\[\\rho \\left( \\frac{\\partial {\\vec V}}{\\partial t} + ({\\vec V} \\cdot \\nabla) {\\vec V} \\right) \\;. \\]<\/p>\n

The third term on the right-hand side vanishes due to the incompressibility assumption.\u00a0 Dividing by the density gives the general form of the Navier-Stokes equations for incompressible fluid flow:<\/p>\n

\\[ \\frac{\\partial {\\vec V}}{\\partial t} + ({\\vec V} \\cdot \\nabla) {\\vec V} = – \\frac{1}{\\rho} \\nabla P + \\nu \\nabla^2 {\\vec V} + {\\vec g} \\; , \\]<\/p>\n

where $$\\nu = \\mu\/\\rho$$ is the kinematic viscosity<\/a>.<\/p>\n

This equation is what we will use in the following examples.<\/p>\n

Gravity-driven Flow<\/h1>\n

This example derives from the very nice YouTube lecture by Rick Sellens<\/a>.\u00a0 Here we imagine flow of a fluid, such as oil, down an inclined plane (flow in the $$x$$ direction), as shown in the figure below.\u00a0 The height of the oil is larger at the top of the ramp but it quickly thins as gravity accelerates the oil.\u00a0 Eventually, the viscous force equals the component of the gravitational force parallel to the plane and, thereafter, the height of the oil is constant (the flow is said to be fully developed) and the velocity profile depends only on $$y$$.\u00a0 Once the flow establishes itself, this height profile remains constant in time and the flow is steady, allowing us to set all time derivatives to zero.\u00a0 We will assume that the incline is wide enough in the $$z$$-direction (not pictured) that we can approximate the flow as being two-dimensional.\u00a0 Next, we will confine our description to the fully developed portion of the flow ($$x > 0$$).\u00a0 Finally, we will assume a no-slip boundary condition at the incline plane and, since the flow is open to the air on top and air is not particularly viscous, we will assume the shear stress on the upper surface to be zero.<\/p>\n

\"\"<\/a><\/p>\n

Applying the steady, fully-developed, and 2-D flow assumptions leads to the continuity equation<\/p>\n

\\[ \\partial_x V_x = 0 \\; \\]<\/p>\n

and the $$x$$- and $$y$$-component Navier-Stokes equations of<\/p>\n

\\[ 0 = g_x + \\nu \\frac{\\partial^2}{\\partial_y^2} V_X \\; \\]<\/p>\n

and<\/p>\n

\\[ 0 = -\\partial_y P + \\rho g_y \\; ,\\]<\/p>\n

where $$g_x = g \\sin(\\theta) $$ and $$g_y = g \\cos(\\theta)$$, are the components of the gravitational acceleration parallel and perpendicular to the surface of the incline, respectively.<\/p>\n

The continuity equation merely confirms what we already assumed, $$V_x = V_x(y)$$.\u00a0 The $$y$$ component of Navier-Stokes says that the pressure at the surface of the incline is higher than at the top due to the weight of the oil above.\u00a0 All the action is in the $$x$$-component equation.<\/p>\n

Integrating the equation once yields<\/p>\n

\\[ \\frac{d}{dy} V_x = – \\frac{g_x}{\\nu} y + C_1 \\; . \\]<\/p>\n

The constant of integration is determined by the shear-stress-free boundary condition at the free surface and gives $$C_1 = g_x h\/\\nu$$.<\/p>\n

Integrating the equation a second time yields<\/p>\n

\\[ V_x = \\frac{g_x}{\\nu} \\left(h y – \\frac{1}{2} y^2 \\right) + C_2 \\; .\\]<\/p>\n

The constant of integration is identically zero after the application of the no-slip boundary condition at the incline ($$V_x(y=0)=0$$).<\/p>\n

Pressure-driven Flow<\/h1>\n

This example draws from the Faculty of Kahn video on Poiseuille\u2019s law<\/a>.\u00a0 In this scenario, we imagine a horizontal pipe with an inlet pressure, $$P_{in}$$, and an outlet pressure, $$P_{out}$$, with $$P_{in} > P_{out}$$ and that the flow has had time to become steady.\u00a0 Cylindrical geometry clearly suggests itself and we expect that the only flow is in the $$z$$-direction ($$V_r = V_{\\theta} = 0$$) and its only variation is radial ($$V_z = V_z(r)$$).<\/p>\n

\"\"<\/a><\/p>\n

As discussed in an earlier post, usually one must take care to differentiate the basis vectors as well as the components when setting up fluid equations in curvilinear coordinates.\u00a0 However, the above assumptions are sufficiently limiting that only the $$z$$-component equation survives in the form<\/p>\n

\\[ \\frac{d P}{d z} = \\nu \\frac{1}{r} \\frac{d}{dr} \\left( r \\frac{d V_z}{dr} \\right) \\; . \\]<\/p>\n

The general solution to this is<\/p>\n

\\[ V_z(r) = \\frac{P_{in} – P_{out}}{4 \\mu L} r^2 + C_1 \\ln r + C_2 \\; .\\]<\/p>\n

The constants of integration follow from a boundary condition that requires $$V_z(r=0)$$ to be finite and the usual no-slip boundary condition at the walls.\u00a0 The final form is then<\/p>\n

\\[ V_z(r) = \\frac{P_{in} – P_{out}}{4 \\mu L}(r^2 – R^2) \\; .\\]<\/p>\n

Integrating this form over the surface area of the pipe gives the famous Poisueille flow rate relationship that says, for a given flow rate, the pressure differential is inversely proportional to the radius to the fourth power.\u00a0 This relationship has profound consequences in the field of medicine.<\/p>\n","protected":false},"excerpt":{"rendered":"

This month\u2019s column deals with viscous flow in a variety of situations.\u00a0 Three main physical effects can drive elementary fluid flow: the pull of gravity coaxing a fluid downhill; a… Read more ><\/a><\/p>\n","protected":false},"author":2,"featured_media":0,"comment_status":"closed","ping_status":"closed","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[1],"tags":[],"class_list":["post-1323","post","type-post","status-publish","format-standard","hentry","category-uncategorized"],"_links":{"self":[{"href":"https:\/\/underthehood.blogwyrm.com\/index.php?rest_route=\/wp\/v2\/posts\/1323","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/underthehood.blogwyrm.com\/index.php?rest_route=\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/underthehood.blogwyrm.com\/index.php?rest_route=\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/underthehood.blogwyrm.com\/index.php?rest_route=\/wp\/v2\/users\/2"}],"replies":[{"embeddable":true,"href":"https:\/\/underthehood.blogwyrm.com\/index.php?rest_route=%2Fwp%2Fv2%2Fcomments&post=1323"}],"version-history":[{"count":16,"href":"https:\/\/underthehood.blogwyrm.com\/index.php?rest_route=\/wp\/v2\/posts\/1323\/revisions"}],"predecessor-version":[{"id":1813,"href":"https:\/\/underthehood.blogwyrm.com\/index.php?rest_route=\/wp\/v2\/posts\/1323\/revisions\/1813"}],"wp:attachment":[{"href":"https:\/\/underthehood.blogwyrm.com\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=1323"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/underthehood.blogwyrm.com\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=1323"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/underthehood.blogwyrm.com\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=1323"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}