Potential flows [PDF]

(Notice that for 2-D incompressible irrotational flows, both velocity potential, φ, and stream- ... where m and C are i

7 downloads 4 Views 393KB Size

Recommend Stories


PDF Book Financial Market Rates and Flows
Everything in the universe is within you. Ask all from yourself. Rumi

PDF Financial Market Rates and Flows
You're not going to master the rest of your life in one day. Just relax. Master the day. Than just keep

Material Flows
Respond to every call that excites your spirit. Rumi

Call Flows
Those who bring sunshine to the lives of others cannot keep it from themselves. J. M. Barrie

Cash Flows
What you seek is seeking you. Rumi

Inviscid Flows
Knock, And He'll open the door. Vanish, And He'll make you shine like the sun. Fall, And He'll raise

Network Flows
You have to expect things of yourself before you can do them. Michael Jordan

Stocks =====> Flows
You can never cross the ocean unless you have the courage to lose sight of the shore. Andrè Gide

Maximum Flows
Life isn't about getting and having, it's about giving and being. Kevin Kruse

Potential
Your task is not to seek for love, but merely to seek and find all the barriers within yourself that

Idea Transcript


Chapter 4

Potential flows Contents

4.1

4.1

Velocity potential . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27

4.2

Kinematic boundary conditions . . . . . . . . . . . . . . . . . . . . . 28

4.3

Elementary potential flows . . . . . . . . . . . . . . . . . . . . . . . . 28

4.4

Properties of Laplace’s equation . . . . . . . . . . . . . . . . . . . . 31

4.5

Flow past an obstacle . . . . . . . . . . . . . . . . . . . . . . . . . . . 32

4.6

Method of images . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35

4.7

Method of separation of variables . . . . . . . . . . . . . . . . . . . . 35

Velocity potential

We shall now consider the special case of irrotational flows, i.e. flows with no vorticity, such that ω = ∇ × u = 0. (4.1) If a velocity field u is irrotational, that is if ∇ × u = 0, then there exists a velocity potential φ(x, t) defined by u = ∇φ. (4.2) This is a result from vector calculus; the converse is trivially true since ∀φ, ∇ × ∇φ ≡ 0. If in addition the flow is incompressible, the velocity potential φ satisfies Laplace’s equation ∇2 φ = 0.

(4.3)

Indeed, for incompressible irrotational flows one has ∇ · u = ∇ · ∇φ = ∇2 φ = 0. Hence, incompressible irrotational flows can be computed by solving Laplace’s equation (4.3) and imposing appropriate boundary conditions (or conditions at infinity) on the solution. (Notice that for 2-D incompressible irrotational flows, both velocity potential, φ, and streamfunction, ψ, are solutions to Laplace’s equation, ∇2 ψ = −ω = 0 and ∇2 φ = 0; boundary conditions on ψ and φ are different however.) 27

28

4.2

4.2 Kinematic boundary conditions

Kinematic boundary conditions n

U

Consider a flow past a solid body moving at velocity U. If n is the unit vector normal to the surface of the solid, then, locally, the surface advances (i.e. moves in the direction of n) at the velocity (U·n) n.

Since the fluid cannot penetrate into the solid body, its velocity normal the surface, (u·n) n, must locally equal that of the solid, u · n = U · n.

So, since u = ∇φ,

∂φ = U · n; (4.4) ∂n the velocity potential satisfies Neumann boundary conditions at the solid body surface. n · ∇φ =

4.3 4.3.1

Elementary potential flows Source and sink of fluid

Line source/sink Consider an axisymmetric potential φ ≡ φ(r). From Laplace’s equation in plane polar coordinates,   1 d m dφ dφ 2 ∇ φ= = , r =0⇔ r dr dr dr r one finds φ(r) = m ln r + C,

(4.5)

where m and C are integration constants. This potential produces the planar radial velocity u = ∇φ =

m ˆ er r

corresponding to a source (m > 0) or sink (m < 0) of fluid of strength m. Notice that the constant C in φ is arbitrary and does not affect u. By convention the constant m = Q/2π where Q is the flow rate. This flow could be produced approximately using a perforated hose.

Source (m > 0)

Sink (m < 0)

29

Chapter 4 – Potential flows Point source/sink

Consider a spherically symmetric potential φ ≡ φ(r). From Laplace’s equation in spherical polar coordinates,   1 d m dφ 2 2 dφ ∇ φ= 2 = 2, r =0⇔ r dr dr dr r one finds φ(r) = −

m + C, r

(4.6)

where m and C are integration constants. This potential produces the three-dimensional radial velocity m u = ∇φ = 2 ˆ er . r corresponding to a source (m > 0) or sink (m < 0) of fluid of strength m. By convention the constant m = Q/4π where Q is the flow rate or volume flux.

Source (m > 0)

4.3.2

Sink (m < 0)

Line vortex

For the potential φ(θ) = kθ, solution to Laplace’s equation in plane polar coordinates, one has ur =

∂φ =0 ∂r

and uθ =

1 ∂φ k = , r ∂θ r

where the strength of the flow k is a constant. By convention k = Γ/2π if Γ is the circulation of the flow. This represents a rotating fluid (bath-plug vortex) around a line vortex at r = 0; it has zero vorticity but is singular at the origin.

4.3.3

Uniform stream

For a uniform flow along the z-axis, u = (0, 0, U ), the velocity potential φ(z) = U z. (The integration constant is set to zero.)

4.3.4

Dipole (doublet flow)

Since Laplace’s equation is linear we can add two solutions together to form a new one. A dipole is the superposition of a sink and an source of equal but opposite strength next to each other.

30

4.3 Elementary potential flows

Three-dimensional flow Consider a point sink of strength −m at the origin and a point source of strength m at the position (0, 0, δ).

Sink

Source z

O δ

The velocity potential of the flow is formed by adding the potentials of the source and sink, m m φ= p −p , 2 2 2 2 2 x + y + (z − δ)2 x +y +z m m −√ = , 2 r r − 2zδ + δ 2

where m > 0 is constant and r = (x2 + y 2 + z 2 )1/2 . Expanding the potential to first order in δ,  m m z φ= − 1 + 2 δ + O(δ 2 ) , r r r and taking the limit δ → 0, leads to the potential of a dipole z φ = −mδ 3 , r where m → ∞ as δ → 0 so that the strength of the dipole µ = mδ remains finite. Thus,   1 µ·r , φ=− 3 =µ·∇ r r

where µ = µˆ ez and r ≡ x is the vector position. The three components of the fluid velocity, u = ∇φ, are for a dipole of strength µ, xz , r5 yz uy = 3µ 5 , r  z2 µ uz = − 3 1 − 3 2 . r r

ux = 3µ

Planar flow Similarly, combining a line sink at the origin with a line source of equal but opposite strength at (δ, 0) gives    m (x − δ)2 + y 2 m 2 2 2 2 ln(x + y ) − ln (x − δ) + y = ln φ=− , m > 0. 2 2 x2 + y 2

As in the three-dimensional case, we consider the limit δ → 0, with µ = mδ fixed. The expression of the potential for a two-dimensional dipole of strength µ then becomes µ·r φ = − 2 = −µ · ∇ ln r, r where µ = µˆ ex and r ≡ x is the vector position.

31

Chapter 4 – Potential flows

4.4

Properties of Laplace’s equation

4.4.1

Identity from vector calculus

Let f (x) be a function defined in a simply connected domain V with boundary S. From vector calculus, ∇ · (f ∇f ) = f ∇2 f + |∇f |2 Z Z Z ⇒ ∇ · (f ∇f ) dV = f ∇2 f dV + |∇f |2 dV.

S n

V

V

V

So, using the divergence theorem Z Z Z 2 f (∇f ) · n dS = f ∇ f dV + |∇f |2 dV. S

4.4.2

V

V

(4.7)

V

Uniqueness of solutions of Laplace’s equation

Given the value of the normal component of the fluid velocity, u · n, on the surface S (i.e. the boundary condition), there exists a unique flow satisfying both ∇ · u = 0 and ∇ × u = 0 (i.e. incompressible and irrotational).

Proof. Suppose there exists two distinct solutions to the boundary value problem, u1 = ∇φ1 and u2 = ∇φ2 . Let f = φ1 − φ2 , then ∇2 f = ∇2 φ 1 − ∇ 2 φ 2 = 0 in the domain V and (∇f ) · n = (∇φ1 ) · n − (∇φ2 ) · n = u1 · n − u2 · n = 0 on the boundary S. Hence, from identity (4.7),

Z

V

|∇φ1 − ∇φ2 |2 dV = 0. However, since

|∇φ1 − ∇φ2 |2 ≥ 0 one must have ∇φ1 = ∇φ2 everywhere. Therefore u1 = u2 and the solution to the boundary value problem is unique. 

4.4.3

Uniqueness for an infinite domain S′

The proof above holds for flows in a finite domain. What about flows in an infinite domain — e.g. flow around an obstacle? The above argument holds by considering the volume V as shown and letting S ′ → ∞. (See, e.g. Patterson p. 211.)

V

11 00 00 11 00 11 00 11 S

n

32

4.4.4

4.5 Flow past an obstacle

Kelvin’s minimum energy theorem

Of all possible fluid motions satisfying the boundary condition for u · n on the surface S and ∇ · u = 0 in domain V , the potential flow is the flow with the smallest kinetic energy, Z 1 K= ρ |u|2 dV. 2 V Proof. Let u′ be another incompressible but non vorticity-free flow such that u · n = u′ · n on S and ∇ · u′ = 0 in V but with ∇ × u′ 6= 0. The fluid flow u is potential, so let u = ∇φ such that Z Z Z 2 2 ρ |u| dV = ρ|∇φ| dV = ρ |∇φ|2 dV, V VZ V = ρ φ u · n dS (by identity (4.7) with f = φ), ZS = ρ φ u′ · n dS (boundary condition), ZS =ρ ∇ · (φ u′ ) dV (divergence theorem), ZV =ρ u′ · ∇φ dV (∇ · u′ = 0), ZV =ρ u′ · u dV. (4.8) V

So, Z

′ 2

V

ρ (u − u ) dV = =

Z

ZV V

 ρ |u|2 − 2ρu · u′ + ρ |u′ |2 dV,

 ρ |u′ |2 − ρ |u|2 dV

(from (4.8)).

Therefore, since (u − u′ )2 ≥ 0, Z Z Z Z ρ |u′ |2 dV = ρ |u|2 dV + ρ (u − u′ )2 dV ≥ ρ |u|2 dV. V

V

V

V



4.5

Flow past an obstacle

Since the solution to Laplace’s equation for given boundary conditions is unique, if we find a solution, we have found the solution. (This is only true if the domain is simply-connected; if the domain is multiply connected, multiple solutions become possible.) One technique to calculate non elementary potential flows involves adding together simple known solutions to Laplace’s equation to get the solution that satisfies the boundary conditions.

4.5.1

Flow around a sphere

We seek an axisymmetric flow of the form u = ur ˆ er + u z ˆ ez in cylindrical polar coordinates (r, θ, z).

33

Chapter 4 – Potential flows n r

U

r

a O

Z

z

At large distances from the sphere of radius a the flow is asymptotic to a uniform stream, ur = 0, uz = U , and at the sphere’s surface, r = a, the fluid velocity must satisfy u · n = 0 since the solid body forms a non-penetrable boundary.

The unit vector normal to surface of the sphere is n = nr ˆ er + n z ˆ ez

with

nr =

r a

and

nz =

z . a

So, the boundary condition u · n = 0 implies that z r ur + uz = 0 ⇔ rur + zuz = 0 a a at the spherical surface of equation r2 + z 2 = a2 . At large distances, the flow is essentially uniform along the z-axis, φ ≃ U z,

for krk ≫ a.

Now, add to the uniform stream a dipole velocity field of strength µ = µˆ ez at the origin, µz , φ(r, z) = U z − 2 (r + z 2 )3/2 so that   ∂φ ∂φ 3z 2 3µrz µ ur = and uz = = 2 =U+ 2 −1 . ∂r ∂z (r + z 2 )5/2 (r + z 2 )3/2 r2 + z 2 Thus, at the sphere’s surface,   z z 3µ(r2 + z 2 ) µ r U+ 2 − , u · n = ur + uz = a a a (r + z 2 )5/2 (r2 + z 2 )3/2     z 2µ 2µ z = U+ 2 U+ 3 , = a a a (r + z 2 )3/2

since r2 + z 2 = a2 . Hence the boundary condition u · n = 0 at the sphere’s surface determines the strength of the dipole, U a3 . µ=− 2 The velocity potential for a uniform flow past a stationary sphere is therefore given by   a3 . (4.9) φ(r, z) = U z 1 + 2(r2 + z 2 )3/2

The corresponding Stokes streamfunction is given by   U r2 a3 Ψ(r, z) = 1− 2 . (4.10) 2 (r + z 2 )3/2

Outside the sphere Ψ > 0, but we also obtain a solution inside the sphere with Ψ < 0. This flow is not real; it is a “virtual flow” that allows for fluid velocity to be consistent with the boundary condition on a solid sphere.

34

4.5.2

4.5 Flow past an obstacle

Rankine half-body

Suppose that, in the velocity potential of a flow past a sphere, we replace the dipole with a point source (m > 0), so that ! m mz mr φ(r, z) = U z − . ⇒ u = ∇φ = ,U + (r2 + z 2 )1/2 (r2 + z 2 )3/2 (r2 + z 2 )3/2 p This flow has a single stagnation point ur = uz = 0 at r = 0 and z = − m/U . To find the streamlines of the flow we calculate the Stokes streamfunction using ur = −

1 ∂Ψ r ∂z

and

uz =

1 ∂Ψ . r ∂r

Thus, ∂Ψ mrz mz U r2 = Ur + − ⇒ Ψ = + α(z), 3/2 2 2 2 ∂r 2 (r + z ) (r + z 2 )1/2 1 ∂Ψ mz 2 α′ (z) m ⇒ + + =− , r ∂z r r (r2 + z 2 )1/2 r (r2 + z 2 )3/2  mr α′ (z) α′ (z) m r2 + z 2 − z 2 = − , + + =− r (r2 + z 2 )3/2 r r (r2 + z 2 )3/2 mr = −ur = − . 2 (r + z 2 )3/2 So, since α′ (z) = 0, α is a constant (set to zero). The Stokes streamfunction is therefore Ψ(r, z) =

U r2 mz . − 2 (r2 + z 2 )1/2

p At the stagnation point (r = 0, z = − m/U ), Ψ = m. Hence, the equation of the streamline, or streamtube, passing through this stagnation point is ! U r2 z Ψ(r, z) = m ⇔ =m 1+ . 2 (r2 + z 2 )1/2 Notice that the straight line r = 0 with z < 0 satisfies the equation of the streamline Ψ = m. For large positive z, the equation of the streamtube Ψ = m becomes r m U r2 ≃ 2m ⇒ r ≃ 2 . 2 U Thus, the velocity potential and the Stokes streamfunction ! ! 2z U a2 a and Ψ(r, z) = φ(r, z) = U z − r2 − 2 4 (r2 + z 2 )1/2 2 (r2 + z 2 )1/2 p provide a model for a long slender body of radius a = 2 m/U .

35

Chapter 4 – Potential flows

4.6

Method of images

In previous examples we introduced flow singularities (e.g. sources and dipoles) outside of the domain of fluid flow in order to satisfy boundary conditions at a solid surface. This technique can also be used to calculate the flow produced by a singularity near a boundary; it is then called method of images. Example 4.1 (Point source near a wall) Consider a point-source of fluid placed at the position (d, 0, 0) (Cartesian coordinates) near a solid wall at x = 0. In free space (no wall), the potential of the source is m

, (x − d)2 + y 2 + z 2 ∂φ∞ m(x − d) = . = ∂x [(x − d)2 + y 2 + z 2 ]3/2

φ∞ = − p

⇒ u∞

y

So that, at x = 0, u∞ = −

O

x

d

md

6= 0, (d2 + y 2 + z 2 )3/2 which is inconsistent with the boundary condition u · n = u · ˆ ex = u = 0 at the wall.

To rectify this problem, (i.e. for the flow to satisfy the boundary condition at the wall), we add a source of equal strength m outside the domain, at (−d, 0, 0). By symmetry, this source will produce an equal but opposite velocity field at x = 0, so that the boundary condition for the combined flow can be satisfied. The velocity potential for both sources becomes m m −p , φ = −p (x − d)2 + y 2 + z 2 (x + d)2 + y 2 + z 2

and the velocity field along the x-axis, u=

m(x − d) ∂φ m(x + d) = + . 3/2 2 2 2 ∂x [(x − d) + y + z ] [(x + d)2 + y 2 + z 2 ]3/2

Clearly, at x = 0, now u = 0 as required. The fluid can slip along the wall however as, for x = 0,

S

2my v= 2 , (d + y 2 + z 2 )3/2 2mz . w= 2 (d + y 2 + z 2 )3/2

−d

4.7

0

d

x

Method of separation of variables

This is a standard method for solving linear partial differential equations with compatible boundary conditions. We shall seek separable solutions to Laplace’s equations, of the form φ(x, y) = f (x)g(y) in Cartesian coordinates or φ(r, θ) = f (r)g(θ) in polar coordinates.

36

4.7 Method of separation of variables

Plane polar coordinates. We substitute a potential of the form φ(r, θ) = f (r)g(θ) in Laplace’s equation expressed in plane polar coordinates,   1 ∂ ∂φ 1 ∂2φ 2 ∇ φ= r + 2 2 = 0, r ∂r ∂r r ∂θ   g d df f d2 g ⇒ r + 2 2 = 0, r dr dr r dθ   df 1 d2 g r d = 0, (division by f (r)g(θ)/r2 ) r + ⇒ f dr dr g dθ2   r d df 1 d2 g ⇒ . r =− f dr dr g dθ2 Since the terms on the left and right sides of the equation are functions of independent variables, r and θ respectively, they must take a constant value, k 2 say. Thus we have transformed a partial differential equation for φ into two ordinary differential equations for f and g,     df d df r d 2 r =k ⇒r r − k 2 f = 0, f dr dr dr dr 1 d2 g d2 g 2 = −k ⇒ + k 2 g = 0. g dθ2 dθ2 Thus, g(θ) = A cos(kθ) + B sin(kθ). For a 2π-periodic function g, such that g(θ) = g(θ + 2π), k must be integer. So g(θ) = A cos(nθ) + B sin(nθ), n ∈ Z, and f is solution to df d2 f +r − n2 f = 0. 2 dr dr Substituting nontrivial functions of the form f = arα gives,   α(α − 1) + α − n2 arα = 0 ⇔ α2 = n2 . r2

The two independent solutions have α = ±n; the general separable solution to Laplace’s equation in plane polar coordinates is therefore   φ(r, θ) = Arn + Br−n cos(nθ) + Crn + Dr−n sin(nθ), n ∈ Z, (4.11) where A, B, C and D are constants to be determined by the boundary conditions.

Separable solutions to Laplace’s equation in spherical polar coordinates can be obtained in a similar manner, but involves Legendre polynomials Pl (cos(θ)). Example 4.2 (Cylinder in an extensional flow) Consider the velocity potential  φ(r, θ) = Ar2 + Br−2 cos(2θ)

corresponding to a particular solution to Laplace’s equation of the form (4.11), with n = 2. The radial velocity of this flow is   B ∂φ = 2r A − 4 cos(2θ). ur = ∂r r

37

Chapter 4 – Potential flows

It vanishes at the surface of a solid cylinder of radius a placed at the origin if B = a4 A. Therefore the velocity field     a4 a4 ur = 2Ar 1 − 4 cos(2θ) and uθ = −2Ar 1 + 4 sin(2θ) r r produced by the potential φ(r, θ) = Ar

2



a4 1+ 4 r



cos(2θ)

represents a fluid flow past a solid cylinder placed in an extensional flow. Notice that at large distances, i.e. if r ≫ a, the fluid velocity is that of an extensional flow ur ≃ 2Ar cos(2θ)

and

uθ ≃ −2Ar sin(2θ),

in polar coordinates, or equivalently u ≃ 2Ax

and

in Cartesian coordinates.

v ≃ −2Ay,

38

4.7 Method of separation of variables

Smile Life

When life gives you a hundred reasons to cry, show life that you have a thousand reasons to smile

Get in touch

© Copyright 2015 - 2024 PDFFOX.COM - All rights reserved.