数学物理方程总览

数学物理方程总览

数学物理方程相关介绍

数学物理方程,通常指从物理学及其他各门自然科学、技术科学中所产生的偏微分方程,有时也包括与此有关的积分方 程、微分积分方程和常微分方程。

数学物理方程的种类: { 常 微 分 方 程 ( O D E ) 偏 微 分 方 程 ( P D E ) 积 分 方 程 \left\{\begin{array}{l} 常微分方程(ODE)\\偏微分方程(PDE)\\ 积分方程\end{array}\right. (ODE)(PDE)

数学物理方程中主要研究二阶线性偏微分方程

二阶线性偏微分方程包括 { ∇ 2 u = f 椭 圆 形 、 泊 松 ( p o i s s o n ) 方 程 u t − D ∇ 2 u = f 抛 物 型 、 扩 散 方 程 u t t − a 2 ∇ 2 u = f 双 曲 型 、 波 动 方 程 ∇ 2 u = 0 拉 普 拉 斯 ( L a p l a c e ) 方 程 ∇ 2 u + k 2 u = f 亥 姆 霍 兹 ( H e l m h o l t z ) 方 程 \left\{\begin{array}{l} \nabla^2 u=f\quad椭圆形、泊松(poisson)方程\\u_t-D\nabla^2u=f\quad 抛物型、扩散方程\\u_{tt}-a^2\nabla^2u=f\quad 双曲型、波动方程\\\nabla^2 u=0 \quad 拉普拉斯(Laplace)方程 \\ \nabla^2u+k^2u=f\quad 亥姆霍兹(Helmholtz)方程 \end{array}\right. 2u=f(poisson)utD2u=futta22u=f2u=0(Laplace)2u+k2u=f(Helmholtz)

相关算符:

梯度算符,Del算符( ∇ \nabla ), ∇ = ( ∂ ∂ x , ∂ ∂ y , ∂ ∂ z ) \nabla=\left(\frac{\partial}{\partial x}, \frac{\partial}{\partial y}, \frac{\partial}{\partial z}\right) =(x,y,z) ,在物理学中被称为哈密顿算子。

其中 ∇ 2 u = ∇ ⋅ ( ∇ u ) , ∇ 2 ≡ ∂ 2 ∂ x 2 + ∂ 2 ∂ y 2 + ∂ 2 ∂ z 2 = ∇ ⋅ ∇ \nabla^{2} u=\nabla \cdot(\nabla u), \quad \nabla^{2} \equiv \frac{\partial^{2}}{\partial x^{2}}+\frac{\partial^{2}}{\partial y^{2}}+\frac{\partial^{2}}{\partial z^{2}}=\nabla \cdot \nabla 2u=(u),2x22+y22+z22=称为拉普拉斯算符/算子,Laplace operator。 Δ 3 \Delta_3 Δ3 表示三维空间的拉普拉斯算符。

物理意义泊松方程 表示静电势或引力势,即电势与电荷密度间的关系

波动方程 表示波的传播,方程 u t t − a 2 ∇ 2 u = f u_{tt}-a^2\nabla^2u=f utta22u=f a a a表示波传播的速度

扩散方程 表示热传导、扩散问题,方程 u t − D ∇ 2 u = f u_{t}-D\nabla^2u=f utD2u=f

二阶线性PDE(偏微分方程)之间的联系

扩散过程达到稳定时, u t = 0 u_t=0 ut=0,温度分布变为Poisson 方程;更进一步,若没有热源,即 f = 0 f=0 f=0,又变为Laplace 方程。

波动方程将时间部分分离则变为亥姆霍兹(Helmholtz)方程【 T ′ ′ ( t ) a 2 T ( t ) = ∇ 2 U ( m ) U ( m ) = − λ \frac{T^{\prime \prime}(t)}{a^{2} T(t)}=\frac{\nabla^{2} U(m)}{U(m)}=-\lambda a2T(t)T(t)=U(m)2U(m)=λ ∇ 2 U ( m ) + λ U ( m ) = 0 \nabla^{2} U(m)+\lambda U(m)=0 2U(m)+λU(m)=0

如果齐次波动方程
∂ 2 u ∂ t 2 − a 2 ∇ 2 u = 0 \frac{\partial^{2} u}{\partial t^{2}}-a^{2} \nabla^{2} u=0 t22ua22u=0
中, u ( x , y , z , t ) u(x, y, z, t) u(x,y,z,t) 随时间周期地变化,频率为 ω \omega ω,
u ( x , y , z , t ) = v ( x , y , z ) e − i ω t u(x, y, z, t)=v(x, y, z) \mathrm{e}^{-\mathrm{i} \omega t} u(x,y,z,t)=v(x,y,z)eiωt
v ( x , y , z ) v(x, y, z) v(x,y,z) 满足 Helmholtz 方程
∇ 2 v ( x , y , z ) + k 2 v ( x , y , z ) = 0 \nabla^{2} v(x, y, z)+k^{2} v(x, y, z)=0 2v(x,y,z)+k2v(x,y,z)=0
其中 k = ω / a k=\omega / a k=ω/a 称为波数。

连续性方程,其中 ρ \rho ρ 是介质的体密度
∂ ρ ∂ t + ∇ ⋅ ρ U = 0 \frac{\partial\rho }{\partial t}+\nabla \cdot \boldsymbol{\rho U}=0 tρ+ρU=0

PDE求解

如常微分方程一样,对于偏微分方程来说,如果我们能求得二阶线性偏微分方程的通解,它应含有两个任意函数。例如,偏微分方程
∂ 2 u ( x , y ) ∂ x 2 = 0 \frac{\partial^{2} u(x, y)}{\partial x^{2}}=0 x22u(x,y)=0
通解就是
u ( x , y ) = c 1 ( y ) + x c 2 ( y ) , u(x, y)=c_{1}(y)+x c_{2}(y), u(x,y)=c1(y)+xc2(y),
其中 c 1 ( y ) c_{1}(y) c1(y) c 2 ( y ) c_{2}(y) c2(y) y y y 的任意函数。

但上面只对介质内部进行了求解,如果介质通过边界面与外部相互作用,或需要描述变化的过程,则除了偏微分方程外,还必须有定解条件,包括初始条件和边界条件。

定解问题 { 范 定 方 程 定 解 条 件 { 初 始 条 件 边 界 条 件 \left\{\begin{array}{l} 范定方程\\定解条件\left\{\begin{array}{l} 初始条件\\边界条件 \\ \end{array}\right. \\ \end{array}\right. {

边界条件

其中 ∂ ∂ n \frac{\partial}{\partial n} n 称为法向微商,它是物理量梯度矢量在外法线方向上的投影,
∂ ∂ n = n ⋅ ∇ ,  即  ∂ u ∂ n = n ⋅ ∇ u .  \frac{\partial}{\partial n}=\boldsymbol{n} \cdot \nabla, \quad \text { 即 } \quad \frac{\partial u}{\partial n}=\boldsymbol{n} \cdot \nabla u \text {. } n=n,  nu=nu
第一类边界条件–狄利克雷(Dirichlet)条件: u ( 0 , t ) = μ ( t ) u(0,t)=\mu(t) u(0,t)=μ(t) 给出边界上各点的函数值;
第二类边界条件–诺依曼边界条件: u x ( 0 , t ) = ν ( t ) u_x(0,t)=\nu(t) ux(0,t)=ν(t) 给出边界上各点函数的法向微商值;
第三类边界条件–罗宾边界条件: u x ( 0 , t ) − h u ( 0 , t ) = θ ( t ) u_x(0,t)-hu(0,t)=\theta(t) ux(0,t)hu(0,t)=θ(t) 给出边界上各点的函数值与法向微商值之间的线性关系。

其他边界条件:内部界面上的连接条件、周期性条件、自然边界条件 R ( r ) = c l r l + d l r l R(r)=c_lr^l+d_lr_l R(r)=clrl+dlrl、无穷远条件、有界条件

初始条件个数应等于方程对时间导数的阶数,n阶导数的通解中需要n个待定参数。

定解问题分类 { 柯 西 问 题 只 有 初 条 , 无 需 边 条 狄 氏 问 题 只 有 边 条 , 无 需 初 条 混 合 形 定 解 问 题 既 有 边 条 , 又 有 初 条 \left\{\begin{array}{l} 柯西问题 \quad 只有初条,无需边条\\狄氏问题 \quad 只有边条,无需初条\\ 混合形定解问题 \quad 既有边条,又有初条\end{array}\right. 西

解的适定性 { 存 在 性 唯 一 性 稳 定 性 \left\{\begin{array}{l} 存在性\\唯一性\\ 稳定性\end{array}\right.

实际计算过程中二阶偏微分方程的解法有:

二阶偏微分方程的解法 { 行 波 方 程 − 行 波 法 − 波 动 方 程 的 达 朗 贝 尔 解 分 离 变 量 法 − 波 动 、 热 传 导 、 直 角 柱 球 坐 标 中 拉 普 拉 斯 方 程 积 分 变 换 法 − 傅 里 叶 变 换 、 拉 普 拉 斯 积 分 变 换 格 林 函 数 − 场 对 点 源 的 响 应 变 分 法 数 值 方 法 − 插 分 法 \left\{\begin{array}{l} 行波方程-行波法-波动方程的达朗贝尔解\\分离变量法-波动、热传导、直角柱球坐标中拉普拉斯方程\\ 积分变换法-傅里叶变换、拉普拉斯积分变换\\格林函数-场对点源的响应 \\变分法 \\数值方法-插分法 \\ \end{array}\right.

接下来我们将对这些解法进行更详细的介绍,点击传送门阅读更多内容。

下一章节:线性偏微分方程的通解

想要直接看分离变量法的同学:分离变量法

数学物理方法专栏阅读指南

posted @ 2022-05-18 21:58  Moster2469  阅读(839)  评论(1编辑  收藏  举报