ConceptComplete

Harmonic Functions

Harmonic functions are solutions to Laplace's equation and arise naturally as the real and imaginary parts of holomorphic functions. They play a central role in potential theory, physics, and the study of conformal mappings.


Definition and Basic Properties

Definition8.1Harmonic function

A C2C^2 function u:Dβ†’Ru: D \to \mathbb{R} on a domain DβŠ†R2D \subseteq \mathbb{R}^2 is harmonic if it satisfies Laplace's equation:

Ξ”u=βˆ‚2uβˆ‚x2+βˆ‚2uβˆ‚y2=0.\Delta u = \frac{\partial^2 u}{\partial x^2} + \frac{\partial^2 u}{\partial y^2} = 0.

Theorem8.1Real and imaginary parts of holomorphic functions

If f=u+ivf = u + iv is holomorphic on a domain DD, then uu and vv are harmonic on DD. Conversely, if uu is harmonic on a simply connected domain DD, there exists a harmonic function vv (unique up to a constant) such that f=u+ivf = u + iv is holomorphic. The function vv is called the harmonic conjugate of uu.

ExampleExamples of harmonic functions
  • u(x,y)=x2βˆ’y2u(x,y) = x^2 - y^2 is harmonic (it is Re(z2)\text{Re}(z^2)), with conjugate v=2xyv = 2xy.
  • u(x,y)=ln⁑x2+y2=ln⁑∣z∣u(x,y) = \ln\sqrt{x^2+y^2} = \ln|z| is harmonic on R2βˆ–{0}\mathbb{R}^2 \setminus \{0\} (it is Re(log⁑z)\text{Re}(\log z)).
  • u(x,y)=excos⁑yu(x,y) = e^x\cos y is harmonic (it is Re(ez)\text{Re}(e^z)), with conjugate v=exsin⁑yv = e^x\sin y.
  • u(x,y)=x/(x2+y2)u(x,y) = x/(x^2+y^2) is harmonic on R2βˆ–{0}\mathbb{R}^2 \setminus \{0\}.

The Cauchy-Riemann Equations

RemarkCauchy-Riemann and harmonicity

The Cauchy-Riemann equations ux=vyu_x = v_y, uy=βˆ’vxu_y = -v_x link harmonic conjugates. Differentiating:

uxx=vyx,uyy=βˆ’vxy.u_{xx} = v_{yx}, \quad u_{yy} = -v_{xy}.

If vv is C2C^2, then vxy=vyxv_{xy} = v_{yx}, so uxx+uyy=0u_{xx} + u_{yy} = 0. The Cauchy-Riemann equations thus imply harmonicity and provide a method to find the conjugate.

ExampleFinding a harmonic conjugate

Given u(x,y)=x3βˆ’3xy2u(x,y) = x^3 - 3xy^2, find vv:

vy=ux=3x2βˆ’3y2β€…β€ŠβŸΉβ€…β€Šv=3x2yβˆ’y3+Ο•(x)v_y = u_x = 3x^2 - 3y^2 \implies v = 3x^2y - y^3 + \phi(x).

vx=6xy+Ο•β€²(x)=βˆ’uy=6xyβ€…β€ŠβŸΉβ€…β€ŠΟ•β€²(x)=0β€…β€ŠβŸΉβ€…β€ŠΟ•=Cv_x = 6xy + \phi'(x) = -u_y = 6xy \implies \phi'(x) = 0 \implies \phi = C.

So v=3x2yβˆ’y3+Cv = 3x^2y - y^3 + C and f=(x3βˆ’3xy2)+i(3x2yβˆ’y3)=z3f = (x^3 - 3xy^2) + i(3x^2y - y^3) = z^3.


Mean Value Property and Maximum Principle

Theorem8.2Mean value property

A continuous function uu on DD is harmonic if and only if it satisfies the mean value property: for every disk Dβ€Ύ(z0,r)βŠ‚D\overline{D}(z_0, r) \subset D,

u(z0)=12Ο€βˆ«02Ο€u(z0+reiΞΈ) dΞΈ.u(z_0) = \frac{1}{2\pi}\int_0^{2\pi} u(z_0 + re^{i\theta})\,d\theta.

Theorem8.3Maximum principle for harmonic functions

If uu is harmonic and non-constant on a domain DD, then uu has no local maximum or minimum in DD. If DD is bounded and uu is continuous on Dβ€Ύ\overline{D}, then

minβ‘βˆ‚Du≀u(z)≀maxβ‘βˆ‚DuforΒ allΒ z∈D.\min_{\partial D} u \leq u(z) \leq \max_{\partial D} u \quad \text{for all } z \in D.

ExampleTemperature distribution

In steady-state heat conduction, the temperature uu is harmonic. The maximum principle says the hottest and coldest points must occur on the boundary β€” the temperature inside is always an average. This has immediate physical content: heat flows from hot to cold until equilibrium is reached.


Poisson Integral Formula

Theorem8.4Poisson integral formula

If uu is harmonic on ∣z∣<R|z| < R and continuous on ∣zβˆ£β‰€R|z| \leq R, then for ∣z∣<R|z| < R:

u(z)=12Ο€βˆ«02Ο€R2βˆ’βˆ£z∣2∣ReiΞΈβˆ’z∣2 u(ReiΞΈ) dΞΈ.u(z) = \frac{1}{2\pi}\int_0^{2\pi} \frac{R^2 - |z|^2}{|Re^{i\theta} - z|^2}\,u(Re^{i\theta})\,d\theta.

The kernel PR(z,ΞΈ)=R2βˆ’βˆ£z∣2∣ReiΞΈβˆ’z∣2P_R(z, \theta) = \frac{R^2 - |z|^2}{|Re^{i\theta} - z|^2} is the Poisson kernel.

RemarkSolving the Dirichlet problem on the disk

The Poisson formula solves the Dirichlet problem: given continuous boundary data gg on ∣z∣=R|z| = R, the function u(z)=12Ο€βˆ«02Ο€PR(z,ΞΈ)g(ReiΞΈ) dΞΈu(z) = \frac{1}{2\pi}\int_0^{2\pi} P_R(z,\theta)g(Re^{i\theta})\,d\theta is harmonic in ∣z∣<R|z| < R and continuous on ∣zβˆ£β‰€R|z| \leq R with u=gu = g on the boundary.