Fluid flow past a cylinder: Difference between revisions

From Citizendium
Jump to navigation Jump to search
imported>Brian Fiedler
mNo edit summary
 
(19 intermediate revisions by 4 users not shown)
Line 1: Line 1:
{{subpages}}
'''Fluid flow past a cylinder''' is classical mathematical solution  for the  flow of an inviscid, incompressible fluid  around a cylinder that is transverse to the flow.  Far from the cylinder, the flow is unidirectional and uniform.
"The '''flow of an incompressible fluid past a cylinder''' is one of the first mathematical models that a student of fluid dynamics encounters. This flow is an excellent vehicle for the study of concepts that will be encountered numerous times in mathematical physics, such as vector fields, coordinate transformations, and most important, the physical interpretation of mathematical results." <ref>http://library.wolfram.com/infocenter/Articles/2731/</ref>
"The '''flow of an incompressible fluid past a cylinder''' is one of the first mathematical models that a student of fluid dynamics encounters. This flow is an excellent vehicle for the study of concepts that will be encountered numerous times in mathematical physics, such as vector fields, coordinate transformations, and most important, the physical interpretation of mathematical results." <ref>http://library.wolfram.com/infocenter/Articles/2731/</ref>


Line 4: Line 7:
{{Image|Flow past a cylinder. Pressure field and velocity vectors.png|right|350px|Colors: pressure field. Red is high and blue is low. Velocity vectors.}}
{{Image|Flow past a cylinder. Pressure field and velocity vectors.png|right|350px|Colors: pressure field. Red is high and blue is low. Velocity vectors.}}
{{Image|Flow past a cylinder. Pressure field and velocity vectors. One quadrant.png|right|350px|Close-up view of one quadrant of the flow. Colors: pressure field. Red is high and blue is low. Velocity vectors.}}
{{Image|Flow past a cylinder. Pressure field and velocity vectors. One quadrant.png|right|350px|Close-up view of one quadrant of the flow. Colors: pressure field. Red is high and blue is low. Velocity vectors.}}
{{Image|Flow past a cylinder. Pressure streamfunction potential.png|right|350px|Pressure field (colors), streamfunction (black), velocity potential (white).}}
{{Image|Flow past a cylinder. Pressure streamfunction potential.png|right|350px|Pressure field (colors), streamfunction (black) with contour interval 0f <math>0.2Ur</math> from bottom to top, velocity potential (white) with contour interval <math>0.2Ur</math> from left to right.}}
A cylinder (or disk) of radius <math>R</math> is placed in two-dimensional, incompressible, inviscid flow.
A cylinder (or disk) of radius <math>R</math> is placed in two-dimensional, incompressible, inviscid flow.
The goal is to find the steady velocity vector <math>\vec{V}</math>  and pressure <math>p</math> in a plane, subject to the condition that
The goal is to find the steady velocity vector <math>\vec{V}</math>  and pressure <math>p</math> in a plane, subject to the condition that
Line 20: Line 23:
The solution for <math>\phi</math> is obtained most easily in polar coordinates <math>r</math> and <math>\theta</math>, related to conventional Cartesian coordinates by <math>x=r\cos\theta</math> and <math>y=r\sin\theta</math>.  In polar coordinates, Laplace's equation is:
The solution for <math>\phi</math> is obtained most easily in polar coordinates <math>r</math> and <math>\theta</math>, related to conventional Cartesian coordinates by <math>x=r\cos\theta</math> and <math>y=r\sin\theta</math>.  In polar coordinates, Laplace's equation is:
:<math>\frac{\partial^2\phi}{\partial r^2} + \frac{1}{r}\frac{\partial \phi}{\partial r} + \frac{1}{r^2} \frac{\partial^2\phi}{\partial \theta^2}=0</math>
:<math>\frac{\partial^2\phi}{\partial r^2} + \frac{1}{r}\frac{\partial \phi}{\partial r} + \frac{1}{r^2} \frac{\partial^2\phi}{\partial \theta^2}=0</math>
The solution that satisfies the boundary conditions is
The solution that satisfies the boundary conditions is <ref> William S. Janna, ''Introduction to Fluid Mechanics'', PWS Publishing Company, Boston (1993)</ref>
:<math>\phi(r,\theta)=U\left(r+\frac{R^2}{r}\right)\cos\theta</math>
:<math>\phi(r,\theta)=U\left(r+\frac{R^2}{r}\right)\cos\theta</math>
The velocity components in polar coordinates are obtained from the components of <math>\nabla\phi</math> in polar coordinates:
The velocity components in polar coordinates are obtained from the components of <math>\nabla\phi</math> in polar coordinates:
Line 28: Line 31:
Being invisicid and irrotational, Bernoulli's equation allows the solution for pressure field to be obtained directly form the velocity field:
Being invisicid and irrotational, Bernoulli's equation allows the solution for pressure field to be obtained directly form the velocity field:
:<math> p=\frac{1}{2}\rho\left(U^2-V^2\right) + p_\infty</math>  
:<math> p=\frac{1}{2}\rho\left(U^2-V^2\right) + p_\infty</math>  
where the constants <math>U</math> and <math>p_\infty</math> appear to that <math>p\rightarrow p_\infty </math> far from the cylinder, where <math>V=U</math>.
where the constants <math>U</math> and <math>p_\infty</math> appear so that <math>p\rightarrow p_\infty </math> far from the cylinder, where <math>V=U</math>.
Using <math>V^2=V_r^2+V_\theta^2</math>,
Using  
:<math>V^2=V_r^2+V_\theta^2</math>,
:<math> p=\frac{1}{2}\rho U^2\left(2\frac{R^2}{r^2}\cos(2\theta)-\frac{R^4}{r^4}\right) + p_\infty </math>
:<math> p=\frac{1}{2}\rho U^2\left(2\frac{R^2}{r^2}\cos(2\theta)-\frac{R^4}{r^4}\right) + p_\infty </math>


Line 35: Line 39:
:<math> 2 \frac{p - p_\infty}{\rho U^2}  =2\frac{R^2}{r^2}\cos(2\theta)-\frac{R^4}{r^4}</math>
:<math> 2 \frac{p - p_\infty}{\rho U^2}  =2\frac{R^2}{r^2}\cos(2\theta)-\frac{R^4}{r^4}</math>
On the surface of the cylinder, or <math>r=R</math>, pressure varies from a maximum of 1 (red color) at the stagnation points at <math>\theta=0</math> and
On the surface of the cylinder, or <math>r=R</math>, pressure varies from a maximum of 1 (red color) at the stagnation points at <math>\theta=0</math> and
<math>\theta=\pi</math> to a minimum of -3 (purple) in the limb of the cylinder at <math>\theta=\frac{1}{2}\pi</math> and <math>\theta=\frac{3}{2}\pi</math>. Likewise, <math>V</math> varies from V=0 at the stagnation points to <math>V=2U</math> on the sides, in the low pressure.
<math>\theta=\pi</math> to a minimum of -3 (purple) on the sides of the cylinder, at <math>\theta=\tfrac{1}{2}\pi</math> and <math>\theta=\tfrac{3}{2}\pi</math>. Likewise, <math>V</math> varies from V=0 at the stagnation points to <math>V=2U</math> on the sides, in the low pressure.


=== Stream function ===
=== Stream function ===
The flow being incompressible, a [http://en.wikipedia.org/wiki/Stream_function  stream function] can be found such that
The flow being incompressible, a [[Stream_function  | stream function]] can be found such that
:<math>\vec{V}=\nabla\psi \times \widehat{k}</math>
:<math>\vec{V}=\nabla\psi \times \widehat{k}</math>
It follows from this definition, using [http://en.wikipedia.org/wiki/List_of_vector_identities vector identities],  
It follows from this definition, using [[Vector_identities |vector identities]],  
:<math>\vec{V}\cdot\nabla{\psi}=0</math>
:<math>\vec{V}\cdot\nabla{\psi}=0</math>
Therefore a contour of a constant value of <math>\psi</math> will also be a [http://en.wikipedia.org/wiki/Streamlines,_streaklines,_and_pathlines stream line], a line tangent to <math>\vec{V}</math>.  For the flow past a cylinder, we find:
Therefore a contour of a constant value of <math>\psi</math> will also be a  stream line, a line tangent to <math>\vec{V}</math>.  For the flow past a cylinder, we find:
:<math> \psi= U \left( r - \frac{R^2}{r} \right) \cos\theta </math>
:<math> \psi= U \left( r - \frac{R^2}{r} \right) \sin\theta </math>


== Physical interpretation ==
== Physical interpretation ==
Laplace's equation  is linear, and is one of the most elementary partial differential equations.  This simple equation yields the entire solution for both <math>\vec{V}</math> and <math>p</math> because of the constraint of irrotation and incompressibility.  Having obtained the solution for <math>\vec{V}</math> and <math>p</math>, the consistency of the pressure gradient with the accelerations can be noted.
The dynamic pressure at the upstream stagnation point has value of <math>\rho U^2/2</math>, a value needed to decelerate the free stream flow of speed U.
This same value appears at the downstream stagnation point, this high pressure is again need to decelerate the flow to zero speed.  This symmetry arises only because the flow is completely frictionless.
The low pressure on sides on the cylinder is need to provide the centripetal acceleration of the flow.
:<math>\frac{\partial p}{\partial r}=\frac{\rho V^2}{L} </math>
where <math>L</math> is the radius of curvature of the flow.  But <math>L \approx R</math>, and <math>V \approx U</math>.  The integral of the equation for centripetal acceleration, which will over a distance <math>\Delta r \approx R</math> will thus yield
:<math>p-p_\infty \approx - \rho U^2</math>
The exact solution has, for the lowest pressure,
:<math>p-p_\infty = - \frac{3}{2} \rho U^2</math>
The low pressure, which must be present to provide the centripetal acceleration, will also increase the flow speed as the fluid travels from higher to lower values of pressure.  Thus we find the maximum speed in the flow, <math>V=2U</math>, in the low pressure on the sides of the cylinder.
A value of <math>V>U</math> is consistent with conservation of the volume of fluid. With the cylinder blocking some of the flow, V must be greater than U somewhere in the plane through the center of the cylinder and transverse to the flow.


== Comparison with flow of a real fluid past a cylinder ==
== Comparison with flow of a real fluid past a cylinder ==
This symmetry of this ideal solution has the peculiar property of having zero net drag on the cylinder, a property known as D'Alembert's paradox.
Unlike an ideal inviscid fluid, a viscous flow past a cylinder, no matter how small the viscosity, will acquire vorticity in a thin boundary layer adjacent to the cylinder.  Boundary layer separation can occur, and a trailing wake will occur behind the cylinder.  The pressure will be lower on the wake side of
the cylinder, than on the upstream side, resulting in a drag force in the downstream direction.


==References==
==References==
<references />
{{reflist}}[[Category:Suggestion Bot Tag]]

Latest revision as of 12:04, 17 August 2024

This article is developing and not approved.
Main Article
Discussion
Related Articles  [?]
Bibliography  [?]
External Links  [?]
Citable Version  [?]
 
This editable Main Article is under development and subject to a disclaimer.

Fluid flow past a cylinder is classical mathematical solution for the flow of an inviscid, incompressible fluid around a cylinder that is transverse to the flow. Far from the cylinder, the flow is unidirectional and uniform.

"The flow of an incompressible fluid past a cylinder is one of the first mathematical models that a student of fluid dynamics encounters. This flow is an excellent vehicle for the study of concepts that will be encountered numerous times in mathematical physics, such as vector fields, coordinate transformations, and most important, the physical interpretation of mathematical results." [1]

Mathematical solution

PD Image
Colors: pressure field. Red is high and blue is low. Velocity vectors.
PD Image
Close-up view of one quadrant of the flow. Colors: pressure field. Red is high and blue is low. Velocity vectors.
Pressure field (colors), streamfunction (black) with contour interval 0f from bottom to top, velocity potential (white) with contour interval from left to right.

A cylinder (or disk) of radius is placed in two-dimensional, incompressible, inviscid flow. The goal is to find the steady velocity vector and pressure in a plane, subject to the condition that far from the cylinder the velocity vector is

where is a constant, and at the boundary of the cylinder

where is vector normal to the cylinder surface. The upstream flow is uniform and has no vorticity. The flow is inviscid, incompressible and has constant mass density . The flow therefore remains without vorticity, or is said to be irrotational, with everywhere. Being irrotational, there must exist a velocity potential :

Being incompressible, , so must satisify Laplace's equation:

The solution for is obtained most easily in polar coordinates and , related to conventional Cartesian coordinates by and . In polar coordinates, Laplace's equation is:

The solution that satisfies the boundary conditions is [2]

The velocity components in polar coordinates are obtained from the components of in polar coordinates:

and

Being invisicid and irrotational, Bernoulli's equation allows the solution for pressure field to be obtained directly form the velocity field:

where the constants and appear so that far from the cylinder, where . Using

,

In the figures, the colorized field referred to as "pressure" is a plot of

On the surface of the cylinder, or , pressure varies from a maximum of 1 (red color) at the stagnation points at and to a minimum of -3 (purple) on the sides of the cylinder, at and . Likewise, varies from V=0 at the stagnation points to on the sides, in the low pressure.

Stream function

The flow being incompressible, a stream function can be found such that

It follows from this definition, using vector identities,

Therefore a contour of a constant value of will also be a stream line, a line tangent to . For the flow past a cylinder, we find:

Physical interpretation

Laplace's equation is linear, and is one of the most elementary partial differential equations. This simple equation yields the entire solution for both and because of the constraint of irrotation and incompressibility. Having obtained the solution for and , the consistency of the pressure gradient with the accelerations can be noted.

The dynamic pressure at the upstream stagnation point has value of , a value needed to decelerate the free stream flow of speed U. This same value appears at the downstream stagnation point, this high pressure is again need to decelerate the flow to zero speed. This symmetry arises only because the flow is completely frictionless.

The low pressure on sides on the cylinder is need to provide the centripetal acceleration of the flow.

where is the radius of curvature of the flow. But , and . The integral of the equation for centripetal acceleration, which will over a distance will thus yield

The exact solution has, for the lowest pressure,

The low pressure, which must be present to provide the centripetal acceleration, will also increase the flow speed as the fluid travels from higher to lower values of pressure. Thus we find the maximum speed in the flow, , in the low pressure on the sides of the cylinder.

A value of is consistent with conservation of the volume of fluid. With the cylinder blocking some of the flow, V must be greater than U somewhere in the plane through the center of the cylinder and transverse to the flow.

Comparison with flow of a real fluid past a cylinder

This symmetry of this ideal solution has the peculiar property of having zero net drag on the cylinder, a property known as D'Alembert's paradox. Unlike an ideal inviscid fluid, a viscous flow past a cylinder, no matter how small the viscosity, will acquire vorticity in a thin boundary layer adjacent to the cylinder. Boundary layer separation can occur, and a trailing wake will occur behind the cylinder. The pressure will be lower on the wake side of the cylinder, than on the upstream side, resulting in a drag force in the downstream direction.

References

  1. http://library.wolfram.com/infocenter/Articles/2731/
  2. William S. Janna, Introduction to Fluid Mechanics, PWS Publishing Company, Boston (1993)