数学物理方法 06 定解问题

$\color{blue}{第六章 定解问题}$

$\color{blue}{\S 6.1 引言}$

$\color{blue}{6.1.1 数理方程简介}$

$1.数学物理方程概念$
$数学物理方程是指从物理、工程问题中,导出的反映客观物理量在各个地点、时刻之间\\ 相互制约关系的一些偏微分方程.$
$数学物理方程 \left \{ \begin{array}{l}线性方程 \\ 非线性方程 \end{array} \right.$

$2.数理方程的产生和发展$
$(1) 18世纪初期: taylor:u_{tt} = a^2 u_{xx} + f$
$(2) 19世纪中期: 三类数学物理方程:$
$\color{blue}{波动方程:u_{tt} = a^2 \Delta u + f}$
$u - 波动, a-波速, f-与源有关的函数$

$\color{blue}{输运方程: u_{t} = D \Delta u + f}$
$u - 浓度, D - 系数, f - 与源有关的已知量$

$\color{blue}{泊松方程: \Delta u = - h}$
$h - 与源有关的已知量, u - 表示稳定物理量$

$(3) 19世纪末到20世纪初$
$\color{blue}{高阶方程(梁的横震动): u_{tt} = a^2 u_{xxxx} + f(x, t)}$
$\color{blue}{非线性方程: KdV: u_t + \sigma u u_x + u_{xxx} = 0}$
$\color{blue}{schr \ddot o -dinger: i\hbar \dfrac{\partial \psi}{\partial t} = - \dfrac{\hbar ^2}{2 u} \Delta \psi + U(r) \psi}$

$\color{blue}{6.1.2 用数理方法研究问题的步骤}$

$1.写出定解问题$
$\left \{ \begin{array}{l} 泛定方程:数理方程(一般规律) \\ 定解条件: 初始、边界、衔接条件(个性) \end{array} \right.$
$如:y^{\prime \prime}(t) - 4y = 0 -- 泛定方程$
$y = C_1e^{2t} +C_2e^{-2t} -- 通解$
$\left \{ \begin{array}{l} y^{\prime \prime} - 4y = 0 \\ \left. \begin{array}{l} y(0) = 0 \\ y^{\prime}(0) = 4 \end{array} \right \} -- 定解条件 \end{array} \right.$

$2求解$
$求解方法:\\ 行波法、分离变量法、积分变换法、格林函数法、保角变换法、复变函数法、变分法$

$3.分析解答:\left \{ \begin{array}{l} 物理意义 \\ 适定性 \left \{ \begin{array}{l} 存在 \\ 唯一 \\ 稳定 \end{array} \right. \end{array} \right.$

$\color{blue}{\S 6.2 三类数理方程的导出}$

$\color{blue}{6.2.1 弦的横震动}$

$1.物理模型:$
$细长而柔软的弦线,紧绷与A、B两点之间,作振幅极微小的横震动,求其运动规律.$

$2.分析:$
$(1) 研究的问题:u(x, t) -- 弦的位移$
$(2) 已知:a. 密度\rho (x, t) = \rho(t),重力 p = 0；$
$b. 无抗弯力$
$c.张力T沿切向;$
$d.u_x是小量, u_x^2 \approx 0$
$(3)研究方法:微积分思想、任意性.$

$3.建立方程:$
$(1)考虑任意段, \Delta x 受力:$
$x: \left \{ \begin{array}{l} -T_1 \cos \alpha_1 \\ T_2 \cos \alpha_2 \end{array} \right.$
$y: \left \{ \begin{array}{l} -T_1 \sin \alpha_1 \\ T_2 \sin \alpha_2 \\ F(x + \eta_1 \Delta x, t) \cdot \Delta x \quad (0 \leq \eta_1 \leq 1) (F是单位长度所受外力) \end{array} \right.$
$(2) 按牛顿运动定律写出方程$
$T_2 \cos \alpha_2 - T_1 \cos \alpha_1 = 0 \quad (1)$
$T_2 \sin \alpha_2 - T_1 \sin \alpha_1 + F(x + \eta_1 \Delta x, t) \Delta x = u_{tt}(x + \eta_2 \Delta x, t) \rho \Delta x \quad (2)$
$(3)简化整理$
$u_{tt} = a^2 u_{xx} + f -- 弦的横震动$
$其中 a^2 = \dfrac{T}{\rho}, 量纲: \dfrac{g \cdot cm/s^2}{g/cm} = (\dfrac{cm}{s})^2$
$f = \dfrac{F}{\rho} -- 单位质量所受力(即力密度)$
$注意:$
$(1) f = 0称为齐次方程: u_{tt} = a^2 u_{xx} + f \to u_{tt} = a^2u_{xx}$
$(2) 三维波动方程: u_{tt} = a^2 \Delta u + f$
$(3)建立方程的步骤:A.从内部划出一小块, B.由物理规律写出算式, C.化简整理得方程$

$\color{blue}{6.2.2 热传导方程}$

$0.热量的几个概念:$
$设:Q - 热量, S - 面积, V - 体积, t - 时间, \rho - 密度, T - 温度, 则:$
$(1) 比热C:单位物质,温度升高一度所需热量. C = \dfrac{Q}{(\rho V) T}$
$(2) 热流密度q: 单位时间流过单位面积的热量. q = \dfrac{Q}{tS}$
$(3)傅里叶实验定理: q = -k \dfrac{\partial T}{\partial n} \quad k - 导热率$
$热流密度与温度的下降率成正比.$
$(4)热源强度F:单位时间,单位体积发出热量. F = \dfrac{Q}{tV}$

$1.物理模型:$
$截面积为A的均匀细杆, 侧面绝热, 沿杆长方向有温差, 求热量的流动.$

$2.分析:$
$(1)研究的问题: u(x, t) -- 温度$
$(2)已知:c、\rho、k是常数; u = u(x, t)是一维问题$
$(3)研究方法:微积分思想、任意性.$

$3建立方程:$
$(1) 考虑任一\Delta x 段在 \Delta t 时间热量情况:$
$流入x面: Q_1 = -k \dfrac{\partial u}{\partial x} | _ {x} \cdot A \Delta t$
$流出x + \Delta x面: Q_2 = -k \dfrac{\partial u}{\partial x} | _ {x + \Delta x} \cdot A \Delta t$
$热源产生: Q_3 = F \cdot \Delta t (A \Delta x) \quad (设其热源强度为F)$
$升温所需热量: Q = C \cdot(\rho A \Delta x) [u(x, t + \Delta t) - u(x, t)]$
$(2) 根据热量守恒定律: Q = Q_1 - Q_2 + Q_3$

$(3)简化整理得: u_t = D u_{xx} + f$
$其中 D = \dfrac{k}{c \rho }, f = \dfrac{F}{c \rho}$
$此即为一维的热传导方程, 中子扩散,高频电流分布皆属于此类方程.$

$\color{blue}{6.2.3 泊松公式}$

$1.物理模型:$
$设在充满了介电常数\varepsilon的区域中, 有体电荷密度为\rho(x, y, z)的电荷,求静电场.$

$2.分析:$
$(1) 研究的问题: \because \vec E = - \nabla V, V -- 标量势$
$(2) 已知:稳定场$
$(3)方法:与上面方法相同$

$3.建立方程:$
$(1)考虑封闭曲面S中的情况$
$(2)由电学中奥-高定理,有: \oint _s \bar E \cdot d \bar s = 4 \pi \cdot \dfrac{1}{ 4 \pi \varepsilon} \int _ {\tau} \rho d \tau$
$(通过一封闭面的静余电通量,等于该平面内所有电荷的代数和)$
$(3)简化整理得:\Delta V = - \dfrac{\rho}{\varepsilon} \to Poisson方程$
$若\rho = 0, 则 \Delta V = 0 \to Laplace方程$
$注意:$
$在稳定温度场中 u_{t} = 0$
$u_{t} = D \Delta u + f \to \Delta u = -\dfrac{f}{D}$

$\color{blue}{\S 6.3 定解条件}$

$引入定解条件的必要性:$
$a.从物理的角度看:数理方程仅能表示一般性$
$b.从数学的角度看:微分方程的解的任意性也需附加条件来确定$
$定解条件 \left \{ \begin{array}{l} 初始条件 \\ 边界条件 \\ 其它条件 \end{array} \right.$

$\color{blue}{6.3.1 初始条件}$

$1.定义:$
$物理过程初始状况的数学表达式为初始条件.$
$弦振动: \left \lbrace \begin{array}{l} u|_{t=0} = \varphi(x) \\ u_t | _{t = 0} = \psi(x) \end{array} \right.$

$2.注意$
$(1)整个系统的初始状况:$
$弦的横震动,当t = 0时$
$\left \lbrace \begin{array}{l} u | _{t=0} = \left \lbrace \begin{array}{l} (2h/l)x \quad (0 \leq x \leq l / 2) \\ (2h/l)(l-x) \quad (l/2 \leq x \leq l) \end{array} \right. \\ u_t | _{t=0} = 0 \end{array} \right.$

$杆的纵震动,当t = 0时$
$u | _{t=0} = \dfrac{b}{l} x$
$(2) 时间t的n阶方程需要n个条件:$
$如:u_{tt} = a^2 u_{xx} + f \to \left \lbrace \begin{array}{l} u| _{t=0} = \varphi(x) \\ u_t| _ {t=0} = \psi(x) \end{array} \right.$
$u_t = D \Delta u + f \to u | _{t = 0 } = \varphi(x)$
$\Delta u = - \dfrac{\rho}{\varepsilon}$

$\color{blue}{6.3.2 边界条件}$

$1.定义:$
$物理过程边界状况的数学表达式为边界条件$

$2.三类边界条件:$
$(1)第一类边界条件(Dirichlet条件)$
$u | _{边} = f(M, t) \quad (已知函数)$
$杆的导热问题: u | _{x = l} = T_0e^{-t}$
$两端固定弦横震动: u | _{x=0}, u | _{x = l} = 0$

$(2)第二类边界条件(Neuman条件)$
$u_n | _{边} = f(M, t) \quad(已知函数)$
$杆的纵震动问题:一端固定,另一端单位面积受力为F(t)$
$u_x | _{x = l} = F / E,$

$杆的导热问题:$
$a) x = l 端有热量流出,热流密度为\psi(t): \dfrac{\partial u}{\partial x} | _{x = l} = - \dfrac{\psi(t)}{k}$
$b)x = l端有热量流入,热流密度为\psi(t): \dfrac{\partial u}{\partial x} | _{x=l} = \dfrac{\psi(t)}{k}$
$c)x=0端有热量流入,热流密度为\psi(t): \dfrac{\partial u}{\partial x} | _{x=0} = - \dfrac{\psi(t)}{k}$

$(3)第三类边界条件(混合边界条件):$
$(u + hu_n) | _{边} = f(M, t) \quad (已知函数)$
$杆的导热问题:一端自由冷却(即牛顿冷却问题)$
$\because [u + hu_x]_{x=l} = u_0, h = \dfrac{k}{H}$
$杆的纵震动问题:若一端固定,一端与弹簧相连$
$u|_{x=0} = 0, (u + \dfrac{Es}{k}u_x)|_{x=l} = 0$

$3.注意：$
$(1)区别边界条件和外源$
$例:长为l的均匀杆,一端固定于x = 0，在t = 0时, 一个沿着杆长方向的力F\\ (单位面积上)加在杆的另一端上,求t > 0时杆长各点的位移。$
$u_x | _ {x = l} = \dfrac{F}{E}$

$(2)一个边界只有一个边界条件$
$例:长为l的均匀杆一端固定于以匀速v前进的撤壁上, 另一端自由, 突然静止,\\ 写出杆做纵震动的定解条件$
$u_x | _{x=l} = 0, u|_{x=0} = 0$

$(3)当f = 0时, 分别称为第1、2、3类齐次边界条件$

$\color{blue}{6.3.3 其它条件}$

$1.衔接条件$
$若所研究的问题由不同部分组成,相接处有衔接条件.$

$杆的纵震动问题:若由两段不同材料组成,则其衔接处有$
$\left \lbrace \begin{array}{l} u_1 | _{x = x_0} = u_2 | _{x=x_0} \\ E_1 \dfrac{\partial u_1}{\partial x} | _ {x = x_0} = E_2 \dfrac{\partial u_2}{\partial x} | _{x = x_0} \end{array} \right.$

$静电场问题:若由两种介质组成,则其衔接处有$
$电势连续: u_1 |_{\sigma} = u_2 | _{\sigma}$
$电位移矢量连续: \varepsilon_1 \dfrac{\partial u_1}{\partial n} | _{\sigma} = \varepsilon_2 \dfrac{\partial u_2}{\partial n} | _ {\sigma}$

$2.自然边界条件:$
$物理问题的有限性、单值性、周期所决定的条件.$
$例:(Euler)方程通解为:$
$y = A x^{l} + B x^{-(l+1)}$
$[0, a]: y = A x_l, \quad [a, \infty): y = B x ^{-(l+1)}$

$\color{blue}{6.3.4 三类定解问题}$

$1.初值问题:$
$只有初始条件而无边界条件的定解问题:$
$例:\left \lbrace \begin{array}{l}u_{tt} = a^2 u_{xx}, -\infty < x < \infty \\ u | _{t=0} = \varphi(x), u_t | _{t=0} = \psi(x) \end{array} \right.$

$2.边值问题:$
$只有边界条件而无初始条件的定解问题:$
$例\left \lbrace \begin{array}{l} \Delta u = 0 \\ u | _{\sigma} = f(M) \end{array} \right.$

$3.混合问题:$
$既有边界条件又有初始条件的定解问题:$
$\left \lbrace \begin{array}{l} u_{tt} = a^2 u_{xx} \\ u| _{x=0} = u| _{x = l} = 0 \\ u | _{t=0} = \varphi(x), u_t | _ {t=0} = \psi(x) \end{array} \right.$