Two-dimensional flow

Two-dimensional flowMain

Community hub

7 pages, 0 posts

0 subscribers

Recent from talks

Be the first to start a discussion here.

Recent from talks

Be the first to start a discussion here.

Contribute something

About hubMembersContent overviewUpdatesRules

Main reference articles

Two-dimensional flow

View on Wikipedia

from Wikipedia

Not found

Revisions and contributors Edit on Wikipedia Read on Wikipedia

View on Grokipedia

from Grokipedia

Two-dimensional flow in fluid dynamics describes the motion of a fluid where the velocity at every point is parallel to a fixed plane and remains uniform along directions perpendicular to that plane, effectively varying only in two spatial dimensions.^[1] This assumption simplifies analysis for scenarios where one dimension (such as width) is significantly larger than the others, allowing flow parameters like velocity, pressure, and density to be treated as independent of the third coordinate.^[2] Key assumptions in two-dimensional flow models often include incompressibility (constant density), irrotationality (zero vorticity, enabling potential flow theory), and inviscidity (negligible viscous effects), which facilitate mathematical tractability.^[1] The velocity field can be expressed using a stream function

\psi(x, y, t)

, where the components are

v_x = \frac{\partial \psi}{\partial y}

and

v_y = -\frac{\partial \psi}{\partial x}

, ensuring the flow is divergence-free for incompressible cases; for irrotational flows,

\psi

satisfies Laplace's equation

\nabla^2 \psi = 0

.^[3] In complex analysis, two-dimensional irrotational and incompressible flows are represented by analytic functions via the complex potential

\Phi(z) = \phi + i\psi

, where

z = x + iy

, linking velocity to

\Phi'(z)

.^[4] Applications of two-dimensional flow theory are prominent in aerospace engineering, such as modeling airflow over airfoils or wings with high aspect ratios, and in hydrodynamics for analyzing flow around ship hulls or in channels.^[2] It also underpins computational simulations and experimental studies of boundary layers, vortices, and uniform flows combined with sources or sinks, providing foundational insights into more complex three-dimensional phenomena.^[4]

Fundamentals

Definition and assumptions

Two-dimensional flow in fluid dynamics describes a type of fluid motion in which the velocity components vary only in two spatial directions, typically within a plane, while assuming uniformity or infinite extent in the third perpendicular direction, as seen in plane flow configurations. This simplification models scenarios where the flow is effectively confined between parallel planes or extends indefinitely without variation along the spanwise axis, allowing the velocity at any point to be identical along lines normal to the flow plane.^[1]^[2]^[5] The analysis of two-dimensional flow relies on several key assumptions to reduce the complexity of the governing equations. These include treating the fluid as incompressible, meaning density remains constant; inviscid, or ideal, where viscosity and frictional effects are neglected; and often steady-state, with time-independent flow unless otherwise specified. Additionally, three-dimensional effects, such as spanwise variations, are disregarded to focus solely on planar dynamics. These assumptions stem from the Euler equations for ideal fluids and enable tractable mathematical models for many engineering applications.^[6]^[4]^[7] The theoretical foundations of two-dimensional flow emerged in the 19th century within hydrodynamics and aerodynamics, building on Leonhard Euler's 18th-century equations for inviscid fluid motion. Pioneering work by figures such as Hermann von Helmholtz and Gustav Kirchhoff advanced the use of potential theory for such flows, laying the groundwork for analytical treatments in subsequent decades. This development was particularly influential in early airfoil design and hydrodynamic problems.^[8]^[9] Unlike three-dimensional flows, which involve variations in all spatial coordinates and often require numerical solutions due to their inherent complexity, two-dimensional flow benefits from reduced dimensionality that permits exact analytical solutions, frequently employing complex variable techniques to represent the velocity field. This distinction facilitates deeper insights into fundamental flow behaviors without the full intricacy of volumetric effects.^[4]^[6]

Velocity field in two dimensions

In two-dimensional flow, the velocity field is represented by a velocity vector

\vec{v}

that varies with position in a plane, typically the $xy$ -plane, and possibly time. In Cartesian coordinates, this is expressed as $\vec{v} = u(x,y,t) \hat{i} + v(x,y,t) \hat{j}$

=u(x,y,t)i^+v(x,y,t)j^, where $u$ and $v$ are the horizontal and vertical velocity components, respectively.^[10] These components describe the local fluid motion at any point $(x,y)$ within the flow domain.^[11] For flows exhibiting radial symmetry or involving circular geometries, such as those around cylinders, the velocity field is conveniently expressed in cylindrical (polar) coordinates as $\vec{v} = v_r(r,\theta,t) \hat{e}_r + v_\theta(r,\theta,t) \hat{e}_\theta$

=vr(r,θ,t)e^r+vθ(r,θ,t)e^θ, where $v_r$ is the radial component and $v_\theta$ is the azimuthal (tangential) component, with no axial velocity in the $z$ -direction for strictly two-dimensional cases.^[10] This representation facilitates analysis in problems where the flow depends on radial distance $r$ and angular position $\theta$ .^[11] To visualize the flow kinematics, pathlines, streamlines, and streaklines are used, each providing distinct insights into particle trajectories. A pathline is the actual path traced by an individual fluid particle over time, governed by the ordinary differential equations $\frac{dx}{dt} = u(x,y,t)$ and $\frac{dy}{dt} = v(x,y,t)$ in Cartesian coordinates.^[12] A streamline is an instantaneous curve tangent to the velocity vector at every point, defined by the differential equation $\frac{dy}{dx} = \frac{v}{u}$ (for $u \neq 0$ ), or equivalently in polar coordinates $\frac{r d\theta}{dr} = \frac{v_\theta}{v_r}$ (for $v_r \neq 0$ ).^[10] A streakline connects the positions of all particles that have passed through a fixed point in the flow at different times, revealing injection patterns but coinciding with streamlines only in steady flows.^[12] In two-dimensional flows, these lines lie within the plane and simplify visualization compared to three-dimensional cases.^[13] The magnitude of the velocity, or speed $q$ , quantifies the flow intensity and is given by $q = \sqrt{u^2 + v^2}$

in Cartesian coordinates or $q = \sqrt{v_r^2 + v_\theta^2}$

in cylindrical coordinates.^[11] The direction of the velocity vector is specified by the angle $\alpha = \tan^{-1}(v/u)$ relative to the $x$ -axis in Cartesian coordinates, indicating the local flow orientation.^[10] These kinematic descriptors assume the flow adheres to the standard two-dimensional approximations, such as uniformity in the third dimension.^[10]

Governing principles

Continuity and momentum equations

In two-dimensional incompressible flow, the continuity equation enforces mass conservation by requiring that the divergence of the velocity field vanishes. For a fluid with constant density

\rho

, this simplifies to the condition that the sum of the partial derivatives of the velocity components with respect to their respective spatial coordinates is zero:

\frac{\partial u}{\partial x} + \frac{\partial v}{\partial y} = 0,

where

u

and

v

are the velocity components in the

x

and

y

directions, respectively.^[14] This equation ensures that the net mass flux through any closed curve in the plane is zero, reflecting the incompressibility assumption that volume is preserved under deformation.^[14] The momentum conservation in two-dimensional inviscid flow is governed by Euler's equations, which describe the acceleration of fluid particles under pressure gradients. In the

x

-direction, the equation is

\frac{\partial u}{\partial t} + u \frac{\partial u}{\partial x} + v \frac{\partial u}{\partial y} = -\frac{1}{\rho} \frac{\partial p}{\partial x},

and in the

y

-direction, it takes the analogous form

\frac{\partial v}{\partial t} + u \frac{\partial v}{\partial x} + v \frac{\partial v}{\partial y} = -\frac{1}{\rho} \frac{\partial p}{\partial y},

where

p

is the pressure and

t

is time.^[14] These nonlinear partial differential equations capture the convective and local acceleration terms balanced against the pressure force per unit mass, assuming no viscous effects or body forces other than pressure.^[14] For steady, inviscid flow along a streamline, integrating Euler's equations yields Bernoulli's equation, which expresses conservation of mechanical energy per unit mass. The equation states that

\frac{p}{\rho} + \frac{1}{2} q^2 + g z = \constant,

where

q = \sqrt{u^2 + v^2}

q=u2+v2

is the speed, $g$ is gravitational acceleration, and $z$ is the elevation (often negligible in horizontal flows).^[15] This relation holds under the assumptions of constant density, no friction, and steady conditions, allowing pressure, kinetic energy, and potential energy to interchange without loss along the flow path.^[15] In irrotational two-dimensional flows, where the vorticity vanishes, a velocity potential $\phi$ exists such that the velocity components are $u = \partial \phi / \partial x$ and $v = \partial \phi / \partial y$ . Substituting into the continuity equation for incompressible flow leads to Laplace's equation for the potential: $\nabla^2 \phi = \frac{\partial^2 \phi}{\partial x^2} + \frac{\partial^2 \phi}{\partial y^2} = 0.$ This elliptic partial differential equation simplifies the analysis of potential flows by reducing the vector velocity field to a scalar potential that satisfies harmonic properties.^[11]

Vorticity in two-dimensional flows

In two-dimensional flows, vorticity is defined as the out-of-plane component of the curl of the velocity field, representing a scalar measure of the local rotation of fluid elements. For a velocity field with components

u

in the

x

-direction and

v

in the

y

-direction, the vorticity

\omega

(specifically

\omega_z

) is given by

\omega = \frac{\partial v}{\partial x} - \frac{\partial u}{\partial y}.

This expression quantifies the tendency of fluid parcels to rotate about the

z

-axis perpendicular to the flow plane.^[8] The sign of vorticity indicates the direction of rotation: positive

\omega

corresponds to counterclockwise rotation, while negative

\omega

corresponds to clockwise rotation, assuming the standard right-hand rule convention for the coordinate system. Zero vorticity implies irrotational flow, where the velocity field has no local spinning motion.^[16]^[17] In cylindrical coordinates

(r, \theta, z)

for two-dimensional flow in the

r

\theta

plane, with radial velocity

v_r

and azimuthal velocity

v_\theta

, the vorticity component

\omega_z

takes the form

\omega = \frac{1}{r} \frac{\partial (r v_\theta)}{\partial r} - \frac{1}{r} \frac{\partial v_r}{\partial \theta}.

This formulation is particularly useful for analyzing flows with radial symmetry, such as those involving vortices or swirling motions.^[17] For incompressible two-dimensional flows, the evolution of vorticity is governed by the vorticity transport equation, which describes how vorticity is advected by the flow and diffused by viscosity:

\frac{D \omega}{Dt} = \nu \nabla^2 \omega,

where

\frac{D}{Dt}

is the material derivative,

\nu

is the kinematic viscosity, and

\nabla^2

is the Laplacian operator. This simplified equation neglects baroclinic torque terms, which arise from density gradients and are absent in constant-density flows; it highlights that vorticity is conserved along fluid particle paths in inviscid conditions but spreads through viscous diffusion. The incompressibility assumption relates to the continuity equation, ensuring the velocity field is divergence-free.^[18]

Potential flow formulation

Velocity potential and stream function

In two-dimensional irrotational flows, the velocity field can be represented by a scalar velocity potential

\phi

, defined such that the velocity vector

\vec{v}

is the gradient of $\phi$ , i.e., $\vec{v} = \nabla \phi$

=∇ϕ.^[11] This representation implies that the curl of the velocity is zero, $\nabla \times \vec{v} = 0$

=0, which is the mathematical condition for irrotationality.^[19] For incompressible flows, where the divergence of the velocity is zero, $\nabla \cdot \vec{v} = 0$

=0, the velocity potential satisfies Laplace's equation, $\nabla^2 \phi = 0$ , or in Cartesian coordinates, $\frac{\partial^2 \phi}{\partial x^2} + \frac{\partial^2 \phi}{\partial y^2} = 0$ .^[11] Solutions to this equation enable analytical determination of the velocity components $u = \frac{\partial \phi}{\partial x}$ and $v = \frac{\partial \phi}{\partial y}$ .^[8] For two-dimensional incompressible flows, the velocity field can alternatively be described using a stream function $\psi$ , a scalar field such that the velocity components are given by $u = -\frac{\partial \psi}{\partial y}$ and $v = \frac{\partial \psi}{\partial x}$ .^[11] These relations automatically satisfy the continuity equation for incompressible flow, $\frac{\partial u}{\partial x} + \frac{\partial v}{\partial y} = 0$ , as they represent the curl of a vector normal to the flow plane.^[19] When the flow is also irrotational, the stream function likewise satisfies Laplace's equation, $\nabla^2 \psi = 0$ , or $\frac{\partial^2 \psi}{\partial x^2} + \frac{\partial^2 \psi}{\partial y^2} = 0$ .^[11] The lines of constant $\psi$ , known as streamlines, trace the direction of the velocity field at every point.^[20] In such flows, the equipotential lines ( $\phi = \constant$ ) are everywhere perpendicular to the streamlines ( $\psi = \constant$ ), forming an orthogonal network that aligns with the flow direction.^[8] This orthogonality arises from the dot product of the gradients being zero, $\nabla \phi \cdot \nabla \psi = 0$ .^[20] The velocity potential $\phi$ and stream function $\psi$ are mathematically linked through the Cauchy-Riemann conditions, $\frac{\partial \phi}{\partial x} = -\frac{\partial \psi}{\partial y}$ and $\frac{\partial \phi}{\partial y} = \frac{\partial \psi}{\partial x}$ , which ensure both functions are harmonic conjugates (up to sign) in the complex plane.^[11] These conditions stem directly from the shared velocity definitions and confirm that the pair $(\phi, \psi)$ represents a valid potential flow solution.^[19]

Complex potential representation

In two-dimensional potential flow, the complex potential

w(z)

provides a unified representation of the flow field using complex analysis, where

z = x + i y

is the complex position variable,

x

and

y

are the Cartesian coordinates, and

w(z) = \phi(x, y) - i \psi(x, y)

. Here,

\phi

is the velocity potential and

\psi

is the stream function, serving as the real and negative imaginary parts, respectively.^[11]^[8] The velocity components are obtained from the derivative of the complex potential as

\frac{dw}{dz} = u - i v

, where

u = \frac{\partial \phi}{\partial x} = -\frac{\partial \psi}{\partial y}

and

v = \frac{\partial \phi}{\partial y} = \frac{\partial \psi}{\partial x}

are the Cartesian velocity components. For

w(z)

to represent a valid potential flow, it must be an analytic function of

z

, meaning it is differentiable in the complex sense everywhere in the flow domain except possibly at singularities. This analyticity implies that

\phi

and

\psi

satisfy the Cauchy-Riemann equations,

\frac{\partial \phi}{\partial x} = -\frac{\partial \psi}{\partial y}

and

\frac{\partial \phi}{\partial y} = \frac{\partial \psi}{\partial x}

, which in turn guarantee that the flow is irrotational (

\nabla \times \mathbf{V} = 0

) and incompressible (

\nabla \cdot \mathbf{V} = 0

), with both potentials satisfying Laplace's equation

\nabla^2 \phi = 0

and

\nabla^2 \psi = 0

.^[21]^[11]^[8] In polar coordinates, the transformation

z = r e^{i \theta}

is particularly useful for flows involving circular or cylindrical geometries, allowing

w(z)

to be expressed in terms of

r

and

\theta

. The corresponding Cauchy-Riemann equations become

\frac{\partial \phi}{\partial r} = -\frac{1}{r} \frac{\partial \psi}{\partial \theta}

and

\frac{1}{r} \frac{\partial \phi}{\partial \theta} = \frac{\partial \psi}{\partial r}

, facilitating analysis of radial and angular variations in the flow.^[11] A key advantage of the complex potential formulation is its exploitation of conformal mapping, where analytic functions preserve local angles and orientations between curves, enabling the transformation of complex flow domains into simpler ones while maintaining the physical flow properties. This simplifies the solution of boundary value problems, such as mapping a circular boundary to an airfoil shape, as demonstrated in classical applications.^[8]^[21]

Basic flow elements

Uniform flow

Uniform flow represents the simplest case of two-dimensional potential flow, characterized by a constant velocity vector

\vec{v} = U \hat{i}

=Ui^, where $U$ is the constant speed in the positive x-direction.^[22] This flow is steady, irrotational, and incompressible, with no spatial variation in velocity magnitude or direction.^[11] The uniformity implies zero divergence and zero vorticity, ensuring the flow aligns with the assumptions of potential theory.^[23] The velocity potential for uniform flow is given by $\phi = U x$ , from which the velocity components are derived as $u = \frac{\partial \phi}{\partial x} = U$ and $v = \frac{\partial \phi}{\partial y} = 0$ .^[22] The corresponding stream function is $\psi = U y$ , satisfying the Cauchy-Riemann conditions with the potential and yielding the same velocity components via $u = \frac{\partial \psi}{\partial y} = U$ and $v = -\frac{\partial \psi}{\partial x} = 0$ .^[11] These functions satisfy Laplace's equation $\nabla^2 \phi = 0$ and $\nabla^2 \psi = 0$ , confirming the irrotational nature.^[23] In complex variable formulation, the complex potential is $w(z) = U z$ , where $z = x + i y$ , combining the real part $\phi = U x$ and imaginary part $\psi = U y$ .^[22] The streamlines, defined by lines of constant $\psi$ , are parallel straight lines perpendicular to the equipotential lines of constant $\phi$ , resulting in horizontal flow paths along the x-direction.^[11]

Source and sink flows

In two-dimensional potential flow, a source represents an idealized point singularity from which fluid emanates radially outward in the plane. The flow is irrotational and incompressible, with fluid particles moving away from the origin along straight radial paths. The strength of the source, denoted by

m

, is defined as the volume flow rate per unit depth perpendicular to the flow plane.^[24] The radial velocity component is given by

v_r = \frac{m}{2\pi r}

, while the tangential component is

v_\theta = 0

, resulting in a purely radial velocity field that diminishes inversely with distance

r

from the source.^[24] The velocity potential for a source is

\phi = \frac{m}{2\pi} \ln r

, satisfying Laplace's equation

\nabla^2 \phi = 0

everywhere except at the origin. The corresponding stream function is

\psi = \frac{m}{2\pi} \theta

, where

\theta

is the polar angle. These functions ensure that the velocity components derive correctly as

v_r = \frac{\partial \phi}{\partial r} = \frac{1}{r} \frac{\partial \psi}{\partial \theta}

and

v_\theta = \frac{1}{r} \frac{\partial \phi}{\partial \theta} = -\frac{\partial \psi}{\partial r}

. In complex variable notation, the complex potential is

w(z) = \frac{m}{2\pi} \ln z

, where

z = x + iy

is the complex position, with

\phi = \operatorname{Re}(w)

and

\psi = \operatorname{Im}(w)

. The streamlines, defined by

\psi = \constant

, are radial lines emanating from the origin at constant

\theta

.^[24] A sink flow is analogous to the source but directed inward toward the singularity, representing fluid absorption. It is modeled by assigning a negative strength

-m

to the source expressions, yielding

v_r = -\frac{m}{2\pi r} < 0

\phi = -\frac{m}{2\pi} \ln r

\psi = -\frac{m}{2\pi} \theta

, and

w(z) = -\frac{m}{2\pi} \ln z

. The streamlines remain radial lines, but fluid motion converges toward the origin. This configuration approximates the two-dimensional limit of a line source or sink extending uniformly in the third dimension.^[24]

Vortex and multipole flows

Irrotational vortex

The irrotational vortex, also known as a potential or free vortex, represents a fundamental element of two-dimensional incompressible potential flow characterized by circulatory motion around a central singularity with constant circulation

\Gamma

. In this flow, the velocity field consists solely of a tangential component, with no radial motion, given by

v_\theta = \frac{\Gamma}{2\pi r}

and

v_r = 0

, where

r

is the radial distance from the origin.^[11] This velocity distribution ensures that the circulation around any closed contour encircling the origin equals

\Gamma

, independent of the contour's size, while the flow remains purely azimuthal and decays inversely with radius.^[11] Despite its circulatory nature, the irrotational vortex satisfies the condition of zero vorticity everywhere except at the origin, where the vorticity manifests as a Dirac delta function singularity to account for the concentrated rotation.^[25] This irrotationality outside the core allows the flow to be described using a velocity potential

\phi

and stream function

\psi

, which are harmonic conjugates in two dimensions. The velocity potential is

\phi = -\frac{\Gamma}{2\pi} \theta

, where

\theta

is the polar angle, and the stream function is

\psi = -\frac{\Gamma}{2\pi} \ln r

.^[11] These functions yield the velocity components through their gradients:

v_r = \frac{\partial \phi}{\partial r} = \frac{\partial \psi}{\partial \theta} = 0

and

v_\theta = \frac{1}{r} \frac{\partial \phi}{\partial \theta} = -\frac{\partial \psi}{\partial r} = \frac{\Gamma}{2\pi r}

.^[11] In complex variable formulation, the irrotational vortex is elegantly represented by the complex potential

w(z) = -\frac{i \Gamma}{2\pi} \ln z

, where

z = x + i y

is the complex position.^[11] Here, the real part of

w(z)

corresponds to

\phi

and the imaginary part to

\psi

, confirming the orthogonal nature of equipotential lines (rays from the origin) and streamlines. The streamlines, defined by

\psi =

constant, form concentric circles centered at the origin, illustrating the purely rotational, non-divergent character of the flow.^[11] This representation underscores the vortex's role as a building block for more complex potential flows via superposition, while highlighting the infinite velocity and singular behavior at

r = 0

.^[11]

Doublet

A two-dimensional doublet arises as the limiting case of a source-sink pair in potential flow, where a source and sink each of strength

m

are separated by a small distance

d

along the x-axis, and the limit is taken as

d \to 0

while keeping the product

\mu = m d

constant; here,

\mu

represents the strength of the doublet.^[8] For an x-directed doublet at the origin, the velocity potential

\phi

and stream function

\psi

are

\phi = -\frac{\mu}{2\pi} \frac{x}{x^2 + y^2},

\psi = \frac{\mu}{2\pi} \frac{y}{x^2 + y^2}.

These satisfy Laplace's equation and the Cauchy-Riemann conditions, ensuring the flow is irrotational and incompressible.^[26] The velocity components in Cartesian coordinates are

u = \frac{\mu}{2\pi} \frac{x^2 - y^2}{(x^2 + y^2)^2}, \quad v = \frac{\mu}{2\pi} \frac{2 x y}{(x^2 + y^2)^2}.

These expressions describe the dipole-like velocity distribution, with singularities at the origin.^[27] In complex notation, the complex potential is

w(z) = -\frac{\mu}{2\pi z},

where

z = x + i y

. The streamlines, given by

\psi = \ constant

, consist of a family of circles tangent to the x-axis at the origin, reflecting the symmetric dipole structure of the flow.^[11]

Applications and combinations

Superposition of basic flows

In two-dimensional potential flow, the superposition principle allows the combination of fundamental flow solutions to construct more complex velocity fields, as the governing Laplace's equation for the velocity potential

\phi

or stream function

\psi

is linear.^[11] This linearity implies that if

\phi_1

and

\phi_2

are solutions satisfying

\nabla^2 \phi = 0

, then the total potential

\phi = \phi_1 + \phi_2

also satisfies the equation, enabling the algebraic addition of individual flows without altering their irrotational and incompressible properties.^[8] Similarly, in complex potential representation, the total

w(z) = w_1(z) + w_2(z)

preserves analyticity, facilitating analytical solutions for composite flows.^[28] A key example is the superposition of uniform flow and a source flow, which models the Rankine half-body—a streamlined shape open at the downstream end, useful for approximating blunt-nosed bodies in aerodynamics.^[29] The uniform flow has potential

\phi_U = U x

, where

U

is the free-stream velocity, while the two-dimensional source at the origin has

\phi_s = \frac{m}{2\pi} \ln r

with strength

m

. The combined potential is

\phi = U x + \frac{m}{2\pi} \ln r,

where

r = \sqrt{x^2 + y^2}

r=x2+y2

.^[11] The corresponding stream function $\psi = U y + \frac{m}{2\pi} \theta$ reveals a stagnation streamline $\psi = 0$ that forms the half-body boundary, with the stagnation point located at $x = -m/(2\pi U)$ on the upstream axis; far downstream, the body half-width asymptotes to $m/(2U)$ .^[29] Another basic superposition involves a source and a sink of equal strength $m$ separated by a small distance, approximating a dipole (or doublet) flow when the separation approaches zero while maintaining constant product of strength and separation.^[11] The source potential $\phi_s = \frac{m}{2\pi} \ln r_1$ and sink $\phi_k = -\frac{m}{2\pi} \ln r_2$ , where $r_1$ and $r_2$ are distances from the points, yield a combined field that, in the limit, produces $\phi_d = \frac{\mu \cos \theta}{2\pi r}$ with dipole moment $\mu$ , representing a flow with radial symmetry useful as a building block for more advanced configurations.^[8] The method of images extends superposition to enforce boundary conditions, such as no normal flow through a wall, by placing mirror-image singularities outside the domain of interest.^[28] For a straight wall at $y=0$ and a source of strength $m$ at $(0, a)$ , an image source of strength $m$ is placed at $(0, -a)$ ; the total potential $\phi = \frac{m}{2\pi} \ln (r_1 r_2)$ ensures $\partial \phi / \partial y = 0$ on the wall by symmetry, where $r_1$ and $r_2$ are distances to the source and image, respectively.^[28] This technique simplifies problems involving solid boundaries, such as flow in channels or near flat plates, by transforming them into unbounded flows with added singularities.^[28]

Flow around airfoils and cylinders

In two-dimensional potential flow, the flow past a circular cylinder of radius

a

in a uniform stream of speed

U

at infinity is obtained by superposing a uniform flow and a doublet. The complex potential for this configuration is given by

w(z) = U \left( z + \frac{a^2}{z} \right),

where

z

is the complex coordinate in the plane perpendicular to the cylinder axis.^[30]^[31] This representation yields streamlines that form closed loops around the cylinder, with no penetration through its surface, and the flow remains symmetric fore and aft, resulting in zero lift and zero drag in the inviscid approximation.^[32]^[33] To generate lift, an irrotational vortex of circulation

\Gamma

is added to the uniform flow plus doublet potential, modeling the Magnus effect observed when a cylinder rotates in a fluid stream. The modified complex potential becomes

w(z) = U \left( z + \frac{a^2}{z} \right) - \frac{i \Gamma}{2\pi} \log z,

which introduces asymmetry in the velocity distribution around the cylinder.^[34]^[35] This circulation shifts higher velocities to one side and lower to the other, producing a pressure difference that yields a lift force per unit length of

L = \rho U \Gamma

, where

\rho

is the fluid density; the force acts perpendicular to the free-stream direction.^[36]^[37] The Joukowski airfoil extends this approach by using conformal mapping to transform the flow around a circular cylinder into flow around an airfoil-shaped body. The mapping is defined as

z = \zeta + \frac{b}{\zeta},

where

\zeta

is the complex coordinate in the circle plane,

z

is in the airfoil plane, and

b

controls the asymmetry, producing a cusped trailing edge when the circle is offset appropriately.^[38]^[39] Applying the cylinder flow potential in the

\zeta

-plane and transforming yields the potential flow over the airfoil, with circulation added to satisfy the Kutta condition of smooth flow off the trailing edge.^[40] The lift on such airfoils, and more generally on any two-dimensional body in potential flow with circulation

\Gamma

, is given by the Kutta-Joukowski theorem:

L = \rho U \Gamma

per unit span, directed perpendicular to the oncoming flow.^[41]^[42] This theorem derives from the Blasius formula for the force on a body, which expresses the complex force as a contour integral of the complex velocity

dw/dz

around the body:

F_x - i F_y = \frac{i \rho}{2} \oint \left( \frac{dw}{dz} \right)^2 dz,

where evaluation for circulatory flows simplifies to the linear relation with

\Gamma

.^[43]^[44] The theorem underscores that lift in inviscid theory arises solely from circulation, independent of body shape details beyond enforcing the Kutta condition.^[45]

Info Pages

Talk Pages

Special Pages

Two-dimensional flow

Recent from talks

Recent from talks

Contribute something

Contribute something

Media Pages

Timelines

Articles

Notes collections

Notes

Notes

Days in Chronicle

Two-dimensional flow