Conformal Mapping

Why conformal mapping matters

• 2D fluid mechanics. Incompressible irrotational flow in any simply connected 2D domain is solvable via conformal mapping to the disc or upper half plane, where the answer is known. Joukowski airfoils, Karman-Trefftz airfoils, and Theodorsen's iteration build on this. Lift, drag, and circulation around realistic airfoil shapes drop out as boundary integrals on the unit circle.
• Electrostatics in 2D. The electric potential φ in a charge-free region satisfies Laplace's equation. Conformal maps transport φ between domains: the field around a square electrode is computed by mapping square → disc, where the answer is the Poisson integral. Capacitor edge effects, microstrip impedances, and parallel-plate field non-uniformity are all conformal-map exercises.
• Mesh parametrization in graphics. Unfolding a 3D mesh to a 2D plane while preserving shape requires a discrete conformal map. Algorithms — Least-Squares Conformal Maps, Discrete Natural Conformal Parametrization, Boundary First Flattening — minimize an energy that vanishes on conformal maps. Texture mapping, normal mapping, and UV unwrapping in production pipelines (Maya, Blender, Houdini) all run on these.
• Complex dynamics — Mandelbrot and Julia sets. The boundary of the Mandelbrot set is conformally equivalent to the unit circle (Douady-Hubbard). External rays — pre-images of radial lines under the uniformizing conformal map — give a combinatorial parametrization of the boundary that powers most of the modern theory of complex dynamics.
• Brain mapping. The cortical surface of the brain is conformally flattened to a disc or sphere for atlas registration. FreeSurfer (Harvard/Mass General) and BrainSuite use conformal maps to compare anatomy across subjects; the maps minimize angular distortion, preserving local shape signals important for fMRI region matching.
• Heat conduction and diffusion. Steady-state heat transport in 2D thin plates with complicated boundary geometries is solved by conformal map to the disc, applying Poisson formula, mapping back. Industrial heat-sink design, thermal bridging in buildings, and chip-cooling layouts use conformal-map software (PLTMG, Driscoll's SC Toolbox).
• String theory worldsheets. The worldsheet of a string is a 2D Riemann surface; conformal invariance of the action is the foundational symmetry of bosonic and superstring theories. Vertex operator algebras, conformal field theories, and modular invariance of partition functions all live in conformal-mapping land.

Common misconceptions

• "Preserves area." No. Only angles are preserved; area scales by |f'(z)|² locally — generically non-uniform. A conformal map can dramatically distort areas while keeping angles intact. Map the upper half plane to the disc: bounded regions become finite-area in the disc but infinite-area in the half plane.
• "Is linear." No. The simplest conformal maps after similarities are the inversions (z ↦ 1/z) and Möbius transformations (z ↦ (az + b)/(cz + d)) — non-linear. The whole point of conformal mapping is bending non-trivial domains while keeping angles; linear maps would be useless.
• "Any holomorphic = conformal." Holomorphic gives angle-preservation only where f'(z) ≠ 0. At zeros of f', the local map is z ↦ z^k (k = order of vanishing), and angles are multiplied by k. Conformality fails at ramification points; an analytic map that has zeros of derivative is conformal except at those points.
• "Conformal in 3D too." In ℝⁿ for n ≥ 3, Liouville's theorem (different from the complex one) says conformal maps are restricted to similarities, inversions, and their compositions — a 10-parameter family in 3D, vs. infinite-dimensional in 2D. The richness of 2D conformal mapping is unique.
• "Maps boundaries pointwise." The Riemann mapping extends continuously to the boundary only under regularity conditions on the boundary (Carathéodory's theorem: the extension is a homeomorphism iff the boundary is a Jordan curve). For domains with cusps, slits, or pathological boundaries, the extension may be merely measurable.
• "Numerically easy." Computing the Riemann map for a given domain is a non-trivial numerical task. The Schwarz-Christoffel formula gives an integral with parameters to be solved (n − 3 conditions for a polygon with n sides); badly-conditioned domains require sophisticated numerics. Driscoll, Trefethen, and Wegmann's algorithms are state of the art.

Definition and basic properties

Let U ⊂ ℂ be open. A function f: U → ℂ is conformal at z₀ ∈ U if it preserves angles at z₀: for any two smooth curves γ₁, γ₂ passing through z₀ at angle θ, the image curves f∘γ₁, f∘γ₂ pass through f(z₀) at the same angle θ.

Theorem: f is conformal at every point of U iff f is holomorphic on U with f'(z) ≠ 0 throughout. Proof in two parts:

• Holomorphic + f'(z) ≠ 0 ⇒ conformal: the linear approximation f(z₀ + h) ≈ f(z₀) + f'(z₀) · h is multiplication by f'(z₀), a scaled rotation. Vectors in any direction are rotated by the same arg f'(z₀); angles are preserved.
• Conformal ⇒ holomorphic with f' ≠ 0: a real-differentiable map is conformal iff its Jacobian at every point is a scaled rotation matrix [[a, −b], [b, a]], iff the Cauchy-Riemann equations hold and the determinant a² + b² is non-zero, iff f is holomorphic with |f'|² ≠ 0.

If f' has zeros, conformality fails at those points; the local behavior is z ↦ z^k for k ≥ 2 (k = order of f' zero plus 1), and angles are multiplied by k. Such points are called "ramification points" or "critical points" of f.

The Riemann Mapping Theorem

Statement: let U ⊂ ℂ be a simply connected open set, U ≠ ℂ. Then there exists a holomorphic bijection f: U → 𝔻 (the open unit disc) with holomorphic inverse — a "conformal equivalence" between U and 𝔻. After fixing f(z₀) = 0 and arg f'(z₀) = 0 for some z₀ ∈ U, the map is unique.

This is one of the most striking theorems in mathematics. It says simply connected open sets in ℂ — apart from ℂ itself — are conformally indistinguishable. A square, a triangle, a snowflake-shaped region, the inside of a Jordan curve, the slit plane — all conformally equivalent to the disc. The disc is "the universal model" up to conformality.

Why the exception "U ≠ ℂ"? Because there is no holomorphic bijection from ℂ to 𝔻: if f: ℂ → 𝔻 were holomorphic, |f| would be bounded (≤ 1), so by Liouville's theorem f would be constant — contradicting bijectivity. ℂ and 𝔻 are conformally distinct.

Sketch of proof (Carathéodory's): consider the family ℱ of holomorphic injections from U into 𝔻 mapping z₀ to 0 with positive derivative. Show ℱ is non-empty and the supremum of f'(z₀) over ℱ is achieved; the maximizer is the Riemann map. Compactness of bounded holomorphic functions (Montel's theorem) and the Schwarz lemma (any holomorphic self-map of 𝔻 fixing 0 has |f'(0)| ≤ 1, with equality iff a rotation) drive the argument.

Möbius transformations

The Möbius (or fractional linear) transformations are

f(z) = (az + b) / (cz + d),    ad − bc ≠ 0.

They form a group under composition: PSL(2, ℂ), the projective linear group. Properties:

• Conformal automorphisms of the Riemann sphere ℂ ∪ {∞} (with the conventions f(∞) = a/c, f(−d/c) = ∞).
• Send circles and lines to circles and lines (treating lines as "circles through ∞"). Three points determine a unique Möbius transformation; preservation of cross-ratio of four points is the defining property.
• Decompose into elementary moves: translation z ↦ z + b, dilation z ↦ az, inversion z ↦ 1/z. Every Möbius is a composition of these.
• Subgroups of importance: SO(3) ≅ rotations of the sphere (Möbius transformations preserving the round metric), PSL(2, ℝ) = automorphisms of the upper half plane, automorphisms of the disc (3-real-parameter family fixing |z| = 1).

Möbius transformations are the most-used building blocks: in geometric group theory (Kleinian groups), in complex dynamics (parabolic, elliptic, hyperbolic classification), in hyperbolic geometry (isometries of ℍ² are PSL(2, ℝ) Möbius transformations).

Schwarz-Christoffel formula

The Schwarz-Christoffel formula constructs explicit conformal maps from the upper half plane (or the disc) to a polygonal region. Suppose the target polygon has interior angles α₁π, α₂π, …, αₙπ at vertices w₁, …, wₙ (with Σ αⱼ = n − 2 for a closed polygon). Choose pre-images x₁ < x₂ < … < xₙ on the real axis. Then

f(z) = A + C · ∫_{z₀}^z (w − x₁)^{α₁ − 1} (w − x₂)^{α₂ − 1} … (w − xₙ)^{αₙ − 1} dw

maps the upper half plane H = {z : Im(z) > 0} to the polygon's interior, sending xⱼ to wⱼ. The constants A, C and the parameters xⱼ are determined by:

• 3 of the xⱼ can be chosen freely (Möbius freedom of the upper half plane).
• n − 3 parameters (xⱼ) and 2 (A, C) are solved numerically to match the polygon's side lengths and orientation.

The factors (w − xⱼ)^{αⱼ − 1} introduce the right corner behavior near each xⱼ via branch-point structure. The integral is single-valued on H provided we use a consistent branch of each factor.

Driscoll's MATLAB toolbox (the SC Toolbox) automates parameter solving and integral evaluation for arbitrary polygons; Trefethen's chebfun ecosystem now includes a Python port. Schwarz-Christoffel is the workhorse for explicit conformal mapping in engineering applications.

Worked examples

Example 1: Map the upper half plane to the unit disc. The Cayley transform

f(z) = (z − i) / (z + i)

is a Möbius transformation. Check: f(i) = 0 (image center), f(0) = (−i)/i = −1 (boundary), f(∞) = 1 (boundary), f(−i) = (−2i)/0 = ∞. It maps H bijectively to 𝔻 conformally.

Example 2: Joukowski transform z ↦ z + 1/z. On the exterior of the unit circle |z| > 1, this is a 2-to-1 map (combine z and 1/z giving the same image), but restricted to one branch (say |z| > 1), it maps to the exterior of the segment [−2, 2]. Slightly offset the original circle (to a circle through z = 1 with center near z = −0.1 + 0.1i), and the image is an airfoil-shaped curve. Flow past a cylinder (textbook problem) maps to flow past the airfoil; lift and circulation drop out by Kutta-Joukowski.

Example 3: Conformal map of the unit disc to itself, fixing 0. By Schwarz lemma, any such map is a rotation z ↦ e^{iθ} z. More generally, conformal automorphisms of 𝔻 are Möbius transformations of the form

f(z) = e^{iθ} (z − a) / (1 − ā z),    a ∈ 𝔻, θ ∈ ℝ.

This 3-real-parameter family is exactly Aut(𝔻) ≅ PSU(1, 1). Same as PSL(2, ℝ) acting on the upper half plane via the Cayley transform.

Example 4: Schwarz-Christoffel for a rectangle. Map the upper half plane to the rectangle with vertices at ±K + iK', ±K − iK' (where K = K(k), K' = K(√(1−k²)) are complete elliptic integrals). Pre-images: x₁ = −1/k, x₂ = −1, x₃ = 1, x₄ = 1/k. Each interior angle is π/2 = (1/2)π, so each αⱼ = 1/2 and αⱼ − 1 = −1/2. The integral is

f(z) = ∫₀^z dw / √[(1 − w²)(1 − k² w²)]

which is the elliptic integral of the first kind. The Riemann map of the upper half plane to a rectangle is exactly this elliptic integral; its inverse is the Jacobi sn function.

Conformal maps and harmonic functions

Conformal maps preserve harmonic functions: if u is harmonic on the target domain V and f: U → V is conformal, then u ∘ f is harmonic on U. Reason: the Laplacian transforms as Δ(u ∘ f)(z) = |f'(z)|² · (Δu)(f(z)); if Δu = 0, then Δ(u ∘ f) = 0.

This is the engineering use case. Solve Laplace's equation on a complicated domain U: conformal-map U to the disc 𝔻 (or half plane), where Laplace's equation is solved by the Poisson integral, then pull back the answer via f. Boundary data on ∂U pulls back to boundary data on ∂𝔻; the Poisson integral solves the Dirichlet problem on the disc; composition with f gives the answer on U.

Concrete: the steady-state temperature in a square plate with prescribed boundary temperatures becomes the Dirichlet problem on a square. Conformally map square → disc via a Schwarz-Christoffel derivative inversion; apply Poisson formula on the disc; map back. The numerical work is in the conformal map; the rest is one integral.

Where conformal mapping shows up

• Aerodynamics — airfoil design. Joukowski (1910) and Karman-Trefftz transforms construct airfoil-shaped conformal images of circles. Flow past the airfoil = flow past the circle in transformed coordinates. Computational tools like XFOIL still implement variants of this.
• Electromagnetics — microstrip and stripline impedance. The characteristic impedance of a microstrip transmission line on a substrate is computed by conformal mapping the cross-section to a parallel-plate geometry where impedance is trivial. Wheeler's 1965 paper and modern tools like Sonnet derive impedance formulas this way.
• Computer graphics — UV mapping. Production tools (Maya, ZBrush, Blender) flatten 3D meshes to 2D for texturing using least-squares conformal maps. The energy ∫ |∂̄ f|² is minimized over piecewise-linear maps; the result is a near-conformal flattening that preserves local shape.
• Brain imaging — surface registration. FreeSurfer (Massachusetts General Hospital) flattens the cortical surface to a sphere via conformal map, then aligns subjects to an atlas. The atlas-comparison step depends on the angle-preserving property to keep functional regions correctly aligned across subjects.
• Geodynamics — mantle flow models. 2D approximations of mantle convection in cross-sections are solved by conformally mapping the spherical shell sector to a rectangle, where finite-difference solvers efficiently compute viscous flow. The map is part of the standard preprocessing pipeline.
• Complex dynamics — Mandelbrot set boundary. Douady and Hubbard proved the boundary of the Mandelbrot set is conformally equivalent to the unit circle's complement. The "external rays" — radial lines pulled back via this map — carry combinatorial information that classifies hyperbolic components and external angles. Thurston's lamination theory builds on this.
• String theory and CFT. The worldsheet of a string is a 2D Riemann surface; conformal invariance of the action is the foundational symmetry. The Virasoro algebra, modular invariance, and partition functions on tori all live in conformal-mapping language. Polyakov's path integral over conformal classes of metrics is the rigorous version.

Moduli of multiply connected domains

The Riemann mapping theorem applies only to simply connected domains. For multiply connected domains (with holes), conformal classification has finite-dimensional moduli.

Annulus: two annuli {r < |z| < R} and {r' < |z| < R'} are conformally equivalent iff R/r = R'/r'. The single invariant μ = (1/2π) log(R/r) — the "modulus" of the annulus — completely classifies up to conformality. Any doubly connected domain is conformally equivalent to a unique annulus.

For domains with n holes (n ≥ 2), 3n − 6 real parameters classify (up to Möbius equivalence). This is the Teichmüller space of n-holed spheres. Riemann surfaces of genus g have Teichmüller space of complex dimension 3g − 3 (for g ≥ 2). The classification of conformal structures on surfaces is a deep area of mathematics with connections to hyperbolic geometry, mathematical physics (string theory), and the Mumford-Mehrotra-Mumford theorems on the moduli space of curves.

Quasi-conformal maps

Conformal maps preserve angles exactly. Quasi-conformal maps preserve angles up to bounded distortion: at every point, the maximum length-to-width ratio of an infinitesimal ellipse mapped from a circle is bounded by a constant K (the dilatation). For K = 1 these are exactly conformal maps.

Quasi-conformal maps form a much larger class with more flexibility — they need not be holomorphic, allowing non-trivial homeomorphisms of multiply connected domains. The Teichmüller theorem says any homeomorphism between Riemann surfaces of the same topology is homotopic to a unique quasi-conformal map of minimal dilatation (Teichmüller geodesic).

Applications: in cartography (Mercator and conformal projections vs. their quasi-conformal generalizations); in cosmology (CMB sky maps); in computer graphics (mesh deformation respecting structure); in dynamical systems (Sullivan's no-wandering-domains theorem uses qc maps to perturb dynamical systems). Bers' simultaneous uniformization theorem combines two Riemann surfaces into a quasi-Fuchsian deformation, foundational for hyperbolic 3-manifolds.