📚 Advanced Mathematics: Types and Solution Methods of Partial Differential Equations | 数学进阶:偏微分方程的类型与解法
Partial differential equations (PDEs) are equations that involve partial derivatives of an unknown function of two or more independent variables. They are fundamental in describing physical phenomena such as heat conduction, wave propagation, electrostatics, and fluid dynamics. This article introduces the main types of PDEs and the most common solution techniques used in advanced mathematics.
偏微分方程(PDE)是含有多个自变量未知函数的偏导数的方程。它们在描述热传导、波动传播、静电学、流体动力学等物理现象中至关重要。本文介绍偏微分方程的主要类型及进阶数学中常用的求解方法。
1. What is a Partial Differential Equation? | 什么是偏微分方程?
A partial differential equation is a relation between an unknown function u(x,y,…) and its partial derivatives. The order of a PDE is the highest order of differentiation appearing in it.
偏微分方程是未知函数 u(x,y,…) 与其偏导数之间的关系式。PDE 的阶数是方程中出现的最高导数的阶数。
For example, the first-order equation ∂u/∂x + ∂u/∂y = 0 is a first-order PDE, while ∂²u/∂x² + ∂²u/∂y² = 0 is a second-order PDE.
例如,一阶方程 ∂u/∂x + ∂u/∂y = 0 是一阶偏微分方程,而 ∂²u/∂x² + ∂²u/∂y² = 0 是二阶偏微分方程。
∂²u/∂x² + ∂²u/∂y² = 0
This is the two-dimensional Laplace equation, the prototypical elliptic equation.
这是二维拉普拉斯方程,是典型的椭圆型方程。
2. Linear and Nonlinear PDEs | 线性与非线性偏微分方程
A linear PDE is one in which the unknown function and all its derivatives appear linearly (to the first power) with coefficients that depend only on the independent variables. Otherwise it is nonlinear.
线性偏微分方程是指未知函数及其所有导数都以一次幂形式出现,且系数仅依赖于自变量的方程。否则为非线性。
Because linear equations obey the superposition principle, their solutions can be combined to build new solutions. This property underlies many solution methods.
线性方程满足叠加原理,因此其解可以组合成新的解。这一性质是许多解法的理论基础。
For instance, the heat equation ∂u/∂t = α ∂²u/∂x² is linear and homogeneous. If u₁ and u₂ are solutions, then c₁u₁ + c₂u₂ is also a solution.
例如,热传导方程 ∂u/∂t = α ∂²u/∂x² 是线性齐次的。若 u₁ 和 u₂ 是解,则 c₁u₁ + c₂u₂ 也是解。
3. Classification: Elliptic, Parabolic, Hyperbolic | 分类:椭圆型、抛物型、双曲型
Second-order linear PDEs in two variables can be classified by the discriminant B² – 4AC using the general form:
两个变量的二阶线性偏微分方程可通过判别式 B² – 4AC 进行分类,其一般形式为:
A ∂²u/∂x² + B ∂²u/∂x∂y + C ∂²u/∂y² + D ∂u/∂x + E ∂u/∂y + F u = G
If B² – 4AC < 0, the equation is elliptic (e.g. Laplace's equation). If B² - 4AC = 0, it is parabolic (e.g. heat equation). If B² - 4AC > 0, it is hyperbolic (e.g. wave equation).
若 B² – 4AC < 0,方程为椭圆型(如拉普拉斯方程);若 B² - 4AC = 0,方程为抛物型(如热传导方程);若 B² - 4AC > 0,方程为双曲型(如波动方程)。
Elliptic equations model steady-state phenomena, parabolic equations describe diffusion processes, and hyperbolic equations represent wave propagation.
椭圆型方程模拟稳态现象,抛物型方程描述扩散过程,双曲型方程表达波动传播。
4. Elliptic Equations: Laplace and Poisson | 椭圆型方程:拉普拉斯与泊松
The Laplace equation is ∇²u = 0, and the Poisson equation is ∇²u = f(x,y), where ∇² is the Laplacian operator. In two dimensions, ∇²u = ∂²u/∂x² + ∂²u/∂y².
拉普拉斯方程为 ∇²u = 0,泊松方程为 ∇²u = f(x,y),其中 ∇² 是拉普拉斯算子。在二维中,∇²u = ∂²u/∂x² + ∂²u/∂y²。
Elliptic equations are usually solved with boundary conditions, not initial conditions. The solution is smooth inside the domain and represents equilibrium states.
椭圆型方程通常使用边界条件而不是初始条件来求解。解在区域内光滑,代表平衡状态。
A typical method is to use separation of variables in rectangular, polar, or spherical coordinates combined with Fourier series.
典型解法是在直角坐标、极坐标或球坐标中使用分离变量法,并结合傅里叶级数。
For example, solving Laplace’s equation on a rectangle with boundary values can be expressed as a sum of sinusoidal eigenfunctions.
例如,在矩形区域上求解拉普拉斯方程并给定边界值,解可表示为正弦本征函数的和。
5. Parabolic Equations: Heat Equation | 抛物型方程:热传导方程
The one-dimensional heat equation is ∂u/∂t = α ∂²u/∂x², where α > 0 is the thermal diffusivity. It involves a first derivative in time and a second derivative in space.
一维热传导方程为 ∂u/∂t = α ∂²u/∂x²,其中 α > 0 是热扩散率。该方程含时间一阶导数和空间二阶导数。
Parabolic problems require an initial condition u(x,0) = f(x) and boundary conditions at both ends of the spatial interval.
抛物型问题需要初始条件 u(x,0) = f(x) 以及空间区间两端的边界条件。
Solutions exhibit smoothing: high-frequency oscillations decay rapidly. This reflects the irreversible nature of diffusion.
解具有平滑化效果:高频振荡迅速衰减。这反映了扩散过程的不可逆性。
The method of separation of variables leads to exponential decay in time multiplied by spatial eigenfunctions, e.g. u(x,t) = Σ Aₙ sin(nπx/L) e^(-α(nπ/L)²t).
分离变量法得出时间上的指数衰减乘以空间本征函数,例如 u(x,t) = Σ Aₙ sin(nπx/L) e^(-α(nπ/L)²t)。
6. Hyperbolic Equations: Wave Equation | 双曲型方程:波动方程
The one-dimensional wave equation is ∂²u/∂t² = c² ∂²u/∂x², where c is the wave speed. It has two initial conditions: initial displacement u(x,0) = f(x) and initial velocity ∂u/∂t(x,0) = g(x).
一维波动方程为 ∂²u/∂t² = c² ∂²u/∂x²,其中 c 是波速。它需要两个初始条件:初始位移 u(x,0) = f(x) 和初始速度 ∂u/∂t(x,0) = g(x)。
The general solution is d’Alembert’s formula: u(x,t) = ½[f(x-ct) + f(x+ct)] + (1/(2c)) ∫ g(s) ds from x-ct to x+ct.
一般解为达朗贝尔公式:u(x,t) = ½[f(x-ct) + f(x+ct)] + (1/(2c)) ∫ g(s) ds(从 x-ct 到 x+ct)。
This formula shows that wave signals propagate along characteristic lines x ± ct = constant, preserving the shape of the initial disturbance.
该公式表明波信号沿特征线 x ± ct = 常数传播,并保持初始扰动的形状。
Hyperbolic equations also conserve energy in many ideal settings, in contrast to parabolic equations.
与抛物型方程不同,双曲型方程在许多理想情况下保持能量守恒。
7. Method of Separation of Variables | 分离变量法
Separation of variables assumes a solution of the form u(x,t) = X(x)T(t). Substituting into the PDE splits it into ordinary differential equations linked by a separation constant λ.
分离变量法假设解的形式为 u(x,t) = X(x)T(t)。将其代入PDE后,偏微分方程分解为常微分方程,并由分离常数 λ 联系。
For example, for the heat equation, we obtain X”/X = T’/(αT) = -λ, leading to spatial ODE X” + λX = 0 and temporal ODE T’ + αλT = 0.
例如,对于热传导方程,可得 X”/X = T’/(αT) = -λ,从而得到空间ODE X” + λX = 0 和时间ODE T’ + αλT = 0。
The steps are: (1) separate variables, (2) solve the resulting eigenvalue problems with boundary conditions, (3) sum the product solutions, and (4) use Fourier coefficients from the initial condition.
步骤是:(1) 分离变量;(2) 在边界条件下求解本征值问题;(3) 叠加乘积解;(4) 利用初始条件确定傅里叶系数。
This method works for linear homogeneous equations with simple geometry, such as rectangular plates, circular disks, and spherical shells.
此方法适用于具有简单几何形状的线性齐次方程,如矩形板、圆盘和球壳。
8. Method of Characteristics | 特征线法
For first-order PDEs of the form a(x,y) ∂u/∂x + b(x,y) ∂u/∂y = c(x,y,u), the method of characteristics reduces the PDE to a system of ordinary differential equations along so-called characteristic curves.
对于形如 a(x,y) ∂u/∂x + b(x,y) ∂u/∂y = c(x,y,u) 的一阶PDE,特征线法沿所谓特征曲线将PDE化为常微分方程组。
For the simple advection equation ∂u/∂t + c ∂u/∂x = 0, the characteristic equations are dx/dt = c and du/dt = 0. Hence u is constant along lines x – ct = constant.
对于简单对流方程 ∂u/∂t + c ∂u/∂x = 0,特征方程为 dx/dt = c 和 du/dt = 0。因此 u 沿直线 x – ct = 常数保持不变。
This method also applies to the wave equation by factoring ∂²u/∂t² – c²∂²u/∂x² = (∂/∂t + c∂/∂x)(∂/∂t – c∂/∂x)u = 0.
该方法也适用于波动方程,可通过因式分解 ∂²u/∂t² – c²∂²u/∂x² = (∂/∂t + c∂/∂x)(∂/∂t – c∂/∂x)u = 0。
Characteristics provide insight into the domain of dependence and range of influence of initial data.
特征线揭示了初始数据的依赖区域和影响区域。
9. Integral Transform Methods | 积分变换法
Fourier and Laplace transforms convert PDEs into algebraic equations or ODEs in the transform variable. This simplifies the differentiation operators.
傅里叶变换和拉普拉斯变换将PDE转换为关于变换变量的代数方程或常微分方程,从而简化微分算子。
For the heat equation on an infinite line, the Fourier transform of ∂u/∂t = α ∂²u/∂x² gives dût/dt = -α k² û(k,t), whose solution is a simple exponential.
对于无限直线上的热传导方程,对 ∂u/∂t = α ∂²u/∂x² 取傅里叶变换得 dût/dt = -α k² û(k,t),其解是简单指数函数。
The inverse transform yields the heat kernel: u(x,t) = (1/√(4παt)) ∫ f(s) exp(-(x-s)²/(4αt)) ds.
逆变换得到热核:u(x,t) = (1/√(4παt)) ∫ f(s) exp(-(x-s)²/(4αt)) ds。
Laplace transforms are especially useful for initial value problems with time as the independent variable, as they incorporate the initial conditions automatically.
拉普拉斯变换特别适用于以时间为自变量的初值问题,因为它自动包含初始条件。
10. Green’s Function Method | 格林函数法
The Green’s function G(x;x₀) represents the response of a linear PDE to a point source. Solving a problem with a distributed source f(x) reduces to integrating f(x₀) against G.
格林函数 G(x;x₀) 表示线性PDE对点源的响应。求解具有分布源 f(x) 的问题可归结为 f(x₀) 与 G 的积分。
For the Poisson equation ∇²u = -f(x)/ε₀, the free-space Green’s function in two dimensions is G = -(1/(2π)) ln|x-x₀|, and in three dimensions it is G = 1/(4π|x-x₀|).
对于泊松方程 ∇²u = -f(x)/ε₀,二维无界域格林函数为 G = -(1/(2π)) ln|x-x₀|,三维为 G = 1/(4π|x-x₀|)。
The solution is u(x) = ∫ G(x;x₀) f(x₀) dx₀, plus terms from boundary conditions.
解为 u(x) = ∫ G(x;x₀) f(x₀) dx₀,并加上边界条件的贡献。
This method is powerful for elliptic and parabolic problems, and it provides a direct physical interpretation of superposition.
该方法对椭圆型和抛物型问题非常有效,并为叠加原理提供了直接的物理解释。
11. Numerical Methods for PDEs | 偏微分方程的数值方法
Most real-world PDEs cannot be solved analytically, so numerical approximations are essential. The finite difference method replaces derivatives by differences on a grid.
大多数实际PDE无法解析求解,因此数值近似至关重要。有限差分法在网格上用差分代替导数。
For example, the explicit scheme for the heat equation is u(i,n+1) = u(i,n) + r(u(i-1,n) – 2u(i,n) + u(i+1,n)), where r = α Δt/Δx².
例如,热传导方程的显式格式为 u(i,n+1) = u(i,n) + r(u(i-1,n) – 2u(i,n) + u(i+1,n)),其中 r = α Δt/Δx²。
Stability requires the CFL condition, e.g. r ≤ ½ for the explicit heat scheme, and r ≤ 1 for the wave equation.
稳定性需要满足CFL条件,例如显式热传导格式需要 r ≤ ½,波动方程需要 r ≤ 1。
The finite element and spectral methods are more flexible for complex geometries and higher accuracy. These methods are widely used in engineering and applied science.
有限元法和谱方法在处理复杂几何形状方面更灵活且精度更高。这些方法在工程和应用科学中广泛使用。
12. Summary and Tips | 总结与技巧
To choose a solution method, first identify the type of PDE using the discriminant, then check the boundary and initial conditions.
选择解法时,首先用判别式确定PDE的类型,然后检查边界条件和初始条件。
Separation of variables works for linear, homogeneous equations on simple domains. Characteristic methods suit first-order and wave equations. Integral transforms and Green’s functions handle unbounded regions and sources.
分离变量法适用于简单区域上的线性齐次方程;特征线法适用于一阶方程和波动方程;积分变换法和格林函数法适用于无界区域和源项。
Always verify that the solution satisfies the original PDE and all conditions. Practice by solving canonical problems: Laplace, heat, and wave equations in various coordinate systems.
始终验证解满足原PDE和所有条件。通过求解典型问题来练习:不同坐标系下的拉普拉斯方程、热传导方程和波动方程。
Mastering PDE types and methods builds a strong foundation for advanced physics, engineering, and financial mathematics.
掌握PDE类型和解法为高等物理、工程和金融数学奠定坚实基础。
Published by TutorHao | Mathematics Revision Series | aleveler.com
更多咨询请联系16621398022(同微信)
屏轩国际教育cambridge primary/secondary checkpoint, cat4, ukiset,ukcat,igcse,alevel,PAT,STEP,MAT, ibdp,ap,ssat,sat,sat2课程辅导,国外大学本科硕士研究生博士课程论文辅导