1.1 Laplace's equation

Laplace's equation is a fundamental concept in potential theory, describing harmonic functions in various fields. It's a second-order partial differential equation that models equilibrium states in physics and engineering, from electrostatics to fluid dynamics.

The equation involves the Laplace operator, which measures local curvature of scalar fields. Solutions, called harmonic functions, have important properties like infinite differentiability. Understanding Laplace's equation is key to solving many boundary value problems in mathematical physics.

Definition of Laplace's equation

• Laplace's equation is a second-order partial differential equation (PDE) that describes the behavior of harmonic functions in various fields, such as electrostatics, fluid dynamics, and heat conduction
• It is named after the French mathematician Pierre-Simon Laplace, who made significant contributions to the study of PDEs and their applications in mathematical physics

Laplace operator

Top images from around the web for Laplace operator

• 

  Laplace distribution - Wikipedia View original

  Is this image relevant?

• 

  Elliptic PDEs - Essential skills for reproducible research computing View original

  Is this image relevant?

• 

  CFD Intro - Essential skills for reproducible research computing View original

  Is this image relevant?

• 

  Laplace distribution - Wikipedia View original

  Is this image relevant?

• 

  Elliptic PDEs - Essential skills for reproducible research computing View original

  Is this image relevant?

1 of 3

Top images from around the web for Laplace operator

• 

  Laplace distribution - Wikipedia View original

  Is this image relevant?

• 

  Elliptic PDEs - Essential skills for reproducible research computing View original

  Is this image relevant?

• 

  CFD Intro - Essential skills for reproducible research computing View original

  Is this image relevant?

• 

  Laplace distribution - Wikipedia View original

  Is this image relevant?

• 

  Elliptic PDEs - Essential skills for reproducible research computing View original

  Is this image relevant?

1 of 3

• The Laplace operator, denoted as $\nabla^2$ or $\Delta$, is a differential operator that acts on a scalar field $\phi(x, y, z)$
• In Cartesian coordinates, the Laplace operator is defined as: $\nabla^2 \phi = \frac{\partial^2 \phi}{\partial x^2} + \frac{\partial^2 \phi}{\partial y^2} + \frac{\partial^2 \phi}{\partial z^2}$
• The Laplace operator measures the local curvature or divergence of a scalar field

Harmonic functions

• A function $\phi(x, y, z)$ is called harmonic if it satisfies Laplace's equation: $\nabla^2 \phi = 0$
• Harmonic functions have several important properties, such as being infinitely differentiable and satisfying the mean value property
• Examples of harmonic functions include the potential functions in electrostatics and the velocity potential in irrotational fluid flow

Cartesian, cylindrical, and spherical coordinates

• Laplace's equation can be expressed in different coordinate systems depending on the geometry of the problem
• In Cartesian coordinates $(x, y, z)$, Laplace's equation takes the form: $\frac{\partial^2 \phi}{\partial x^2} + \frac{\partial^2 \phi}{\partial y^2} + \frac{\partial^2 \phi}{\partial z^2} = 0$
• In cylindrical coordinates $(r, \theta, z)$, Laplace's equation is: $\frac{1}{r} \frac{\partial}{\partial r} \left(r \frac{\partial \phi}{\partial r}\right) + \frac{1}{r^2} \frac{\partial^2 \phi}{\partial \theta^2} + \frac{\partial^2 \phi}{\partial z^2} = 0$
• In spherical coordinates $(r, \theta, \phi)$, Laplace's equation becomes: $\frac{1}{r^2} \frac{\partial}{\partial r} \left(r^2 \frac{\partial \phi}{\partial r}\right) + \frac{1}{r^2 \sin \theta} \frac{\partial}{\partial \theta} \left(\sin \theta \frac{\partial \phi}{\partial \theta}\right) + \frac{1}{r^2 \sin^2 \theta} \frac{\partial^2 \phi}{\partial \phi^2} = 0$

Derivation of Laplace's equation

• Laplace's equation can be derived from the fundamental concepts of vector calculus, such as the gradient and divergence operators
• The derivation involves applying these operators to a scalar field and setting the result equal to zero

Gradient of a scalar field

• The gradient of a scalar field $\phi(x, y, z)$ is a vector field that points in the direction of the greatest rate of increase of $\phi$ and has a magnitude equal to the rate of change
• In Cartesian coordinates, the gradient is given by: $\nabla \phi = \left(\frac{\partial \phi}{\partial x}, \frac{\partial \phi}{\partial y}, \frac{\partial \phi}{\partial z}\right)$
• The gradient is used to calculate the electric field from an electric potential or the velocity field from a velocity potential

Divergence of a vector field

• The divergence of a vector field $\mathbf{F}(x, y, z) = (F_x, F_y, F_z)$ measures the net outward flux of the field per unit volume
• In Cartesian coordinates, the divergence is defined as: $\nabla \cdot \mathbf{F} = \frac{\partial F_x}{\partial x} + \frac{\partial F_y}{\partial y} + \frac{\partial F_z}{\partial z}$
• The divergence appears in physical laws such as Gauss's law in electrostatics and the continuity equation in fluid dynamics

Laplacian operator

• The Laplacian operator $\nabla^2$ is obtained by taking the divergence of the gradient of a scalar field: $\nabla^2 \phi = \nabla \cdot (\nabla \phi)$
• Expanding this expression using the definitions of gradient and divergence leads to the Laplacian in Cartesian coordinates: $\nabla^2 \phi = \frac{\partial^2 \phi}{\partial x^2} + \frac{\partial^2 \phi}{\partial y^2} + \frac{\partial^2 \phi}{\partial z^2}$
• Setting the Laplacian equal to zero yields Laplace's equation, which describes the behavior of harmonic functions

Fundamental solutions

• Fundamental solutions are specific solutions to Laplace's equation that can be used to construct more general solutions using the principle of superposition
• These solutions often have a physical interpretation and serve as building blocks for solving boundary value problems

Green's functions

• Green's functions are fundamental solutions that represent the response of a system to a point source or a delta function input
• For Laplace's equation in 3D, the Green's function is given by: $G(\mathbf{r}, \mathbf{r}') = -\frac{1}{4\pi|\mathbf{r} - \mathbf{r}'|}$, where $\mathbf{r}$ and $\mathbf{r}'$ are position vectors
• Green's functions can be used to solve inhomogeneous PDEs, such as Poisson's equation, by convolving the Green's function with the source term

Method of images

• The method of images is a technique for solving Laplace's equation in the presence of boundaries by introducing fictitious sources or sinks that satisfy the boundary conditions
• The solution is then obtained by superimposing the fields generated by the real and imaginary sources
• Examples of the method of images include the electric field of a point charge near a grounded conducting plane and the potential flow around a circular cylinder

Poisson's equation

• Poisson's equation is a generalization of Laplace's equation that includes a source term: $\nabla^2 \phi = -\rho$, where $\rho$ is the source density
• In electrostatics, Poisson's equation relates the electric potential $\phi$ to the charge density $\rho$: $\nabla^2 \phi = -\frac{\rho}{\epsilon_0}$, where $\epsilon_0$ is the permittivity of free space
• Poisson's equation can be solved using Green's functions or by first finding a particular solution and then adding a homogeneous solution that satisfies the boundary conditions

Boundary value problems

• Boundary value problems (BVPs) involve solving Laplace's equation subject to specified conditions on the boundaries of the domain
• The solution to a BVP is unique if the boundary conditions are well-posed and the domain is sufficiently regular

Dirichlet boundary conditions

• Dirichlet boundary conditions specify the value of the function $\phi$ on the boundary of the domain
• For example, in electrostatics, Dirichlet conditions correspond to fixed potentials on conducting surfaces
• Dirichlet problems can be solved using techniques such as separation of variables, Green's functions, or numerical methods

Neumann boundary conditions

• Neumann boundary conditions specify the normal derivative of the function $\phi$ on the boundary of the domain
• In physical problems, Neumann conditions often represent the flux of a quantity across the boundary, such as the electric field or the heat flux
• Neumann problems can be solved using similar techniques as Dirichlet problems, but the solutions may not be unique up to an additive constant

Mixed boundary conditions

• Mixed boundary conditions involve a combination of Dirichlet and Neumann conditions on different parts of the boundary
• For example, a conducting surface may have a fixed potential on one part and a specified charge density on another part
• Mixed problems can be more challenging to solve and may require the use of specialized techniques, such as conformal mapping or the Wiener-Hopf method

Uniqueness of solutions

• The uniqueness of solutions to Laplace's equation with given boundary conditions is an important property that ensures the well-posedness of the problem
• For Dirichlet problems, the solution is unique if the boundary values are continuous and the domain is connected
• For Neumann problems, the solution is unique up to an additive constant if the boundary values are consistent with the divergence theorem and the domain is connected

Separation of variables

• Separation of variables is a powerful technique for solving Laplace's equation in separable coordinate systems, such as rectangular, cylindrical, and spherical coordinates
• The method involves assuming that the solution can be written as a product of functions, each depending on only one variable, and then solving the resulting ordinary differential equations (ODEs)

Rectangular coordinates

• In rectangular coordinates $(x, y, z)$, the separated solution takes the form: $\phi(x, y, z) = X(x)Y(y)Z(z)$
• Substituting this ansatz into Laplace's equation leads to three ODEs: $X''(x) = -\lambda_x X(x)$, $Y''(y) = -\lambda_y Y(y)$, and $Z''(z) = -\lambda_z Z(z)$, where $\lambda_x$, $\lambda_y$, and $\lambda_z$ are separation constants
• The solutions to these ODEs are trigonometric or exponential functions, depending on the signs of the separation constants

Cylindrical coordinates

• In cylindrical coordinates $(r, \theta, z)$, the separated solution is: $\phi(r, \theta, z) = R(r)\Theta(\theta)Z(z)$
• The resulting ODEs are: $r^2 R''(r) + rR'(r) - (\lambda_r + m^2)R(r) = 0$, $\Theta''(\theta) = -m^2 \Theta(\theta)$, and $Z''(z) = -\lambda_z Z(z)$, where $m$ is an integer and $\lambda_r$ and $\lambda_z$ are separation constants
• The solutions involve Bessel functions for $R(r)$, trigonometric functions for $\Theta(\theta)$, and exponential functions for $Z(z)$

Spherical coordinates

• In spherical coordinates $(r, \theta, \phi)$, the separated solution is: $\phi(r, \theta, \phi) = R(r)\Theta(\theta)\Phi(\phi)$
• The ODEs are: $r^2 R''(r) + 2rR'(r) - [l(l+1)]R(r) = 0$, $\sin \theta \frac{d}{d\theta} \left(\sin \theta \frac{d\Theta}{d\theta}\right) + [l(l+1)\sin^2 \theta - m^2]\Theta = 0$, and $\Phi''(\phi) = -m^2 \Phi(\phi)$, where $l$ and $m$ are integers
• The solutions involve spherical Bessel functions for $R(r)$, associated Legendre functions for $\Theta(\theta)$, and trigonometric functions for $\Phi(\phi)$

Eigenfunction expansions

• The separated solutions form a complete set of eigenfunctions for the Laplacian operator in the given coordinate system
• The general solution to Laplace's equation can be expressed as an infinite series of these eigenfunctions, with coefficients determined by the boundary conditions
• Eigenfunction expansions provide a systematic way to solve boundary value problems and analyze the properties of the solutions, such as symmetry and asymptotic behavior

Laplace transforms

• Laplace transforms are an integral transform technique that converts a function of a real variable $t$ to a function of a complex variable $s$
• Laplace transforms can be used to solve initial value problems for linear ODEs and PDEs, including Laplace's equation

Definition and properties

• The Laplace transform of a function $f(t)$ is defined as: $F(s) = \mathcal{L}\{f(t)\} = \int_0^\infty f(t)e^{-st} dt$, where $s$ is a complex number
• Laplace transforms have several important properties, such as linearity, scaling, shifting, and differentiation
• The convolution theorem states that the Laplace transform of the convolution of two functions is the product of their Laplace transforms

Application to Laplace's equation

• Laplace transforms can be applied to Laplace's equation by taking the transform with respect to one of the variables, usually time in time-dependent problems
• The transformed equation is an ODE in the remaining spatial variables, which can be solved using standard techniques
• The solution in the transformed domain is then inverted to obtain the solution in the original domain

Inversion of Laplace transforms

• Inverting a Laplace transform involves finding the original function $f(t)$ given its transform $F(s)$
• The inverse Laplace transform is defined as: $f(t) = \mathcal{L}^{-1}\{F(s)\} = \frac{1}{2\pi i} \int_{\gamma-i\infty}^{\gamma+i\infty} F(s)e^{st} ds$, where $\gamma$ is a real number greater than the real part of all singularities of $F(s)$
• Inversion can be performed using techniques such as partial fraction decomposition, complex contour integration, or numerical methods

Numerical methods

• Numerical methods are computational techniques for solving Laplace's equation and related problems when analytical solutions are not available or practical
• These methods discretize the domain into a finite number of elements or points and approximate the solution using algebraic equations or iterative algorithms

Finite difference methods

• Finite difference methods approximate the derivatives in Laplace's equation using finite differences based on the values of the function at neighboring grid points
• The resulting system of linear equations is solved using methods such as Jacobi iteration, Gauss-Seidel iteration, or successive over-relaxation (SOR)
• Finite difference methods are simple to implement but may suffer from issues such as numerical dispersion and difficulty in handling complex geometries

Finite element methods

• Finite element methods (FEM) divide the domain into a mesh of elements, such as triangles or tetrahedra, and approximate the solution using piecewise polynomial functions on each element
• The weak form of Laplace's equation is obtained by multiplying the equation by a test function and integrating by parts, which leads to a system of linear equations for the nodal values of the solution
• FEM is more flexible than finite difference methods in handling irregular domains and boundary conditions but requires more computational resources

Boundary element methods

• Boundary element methods (BEM) reformulate Laplace's equation as an integral equation on the boundary of the domain using Green's identities
• The boundary is discretized into elements, and the solution is represented using a linear combination of fundamental solutions or boundary integrals
• BEM reduces the dimensionality of the problem by one and is particularly efficient for problems with small surface-to-volume ratios, but it leads to dense matrices that may be computationally expensive to solve

Applications of Laplace's equation

• Laplace's equation arises in numerous physical and engineering applications, where it describes the behavior of potential fields, steady-state diffusion, and other equilibrium phenomena

Electrostatics and magnetostatics

• In electrostatics, Laplace's equation governs the electric potential $\phi$ in charge-free regions: $\nabla^2 \phi = 0$
• The electric field is obtained from the potential as: $\mathbf{E} = -\nabla \phi$
• In magnetostatics, Laplace's equation describes the magnetic scalar potential $\psi$ in current-free regions: $\nabla^2 \psi = 0$, and the magnetic field is given by: $\mathbf{H} = -\nabla \psi$

Heat conduction and diffusion

• In steady-state heat conduction, Laplace's equation relates the temperature $T$ to the thermal conductivity $k$: $\nabla \cdot (k\nabla T) = 0$
• For constant thermal conductivity, this reduces to: $\nabla^2 T = 0$
• In diffusion problems, such as the concentration of a solute in a solvent, Laplace's equation describes the equilibrium distribution of the concentration $c$: $\nabla^2 c = 0$

Fluid dynamics and potential flow

• In irrot