Skip to content

方法Method

拟线性方程的特征线法

Method of characteristics · Noncharacteristic Cauchy problem · 特征线法

把一阶拟线性方程沿曲线变成扩展ODE,以非特征初始面和位置投影的局部可逆性构造唯一经典解,并区分切向数据与投影退化。

如果温度随流体移动,我们可以追着一个流体粒子看温度如何变化。这样得到的是一条时间曲线,而不是整个空间的函数。特征线法先解这些曲线,再检查它们能否拼成一张单值的函数图像;最后这一步正是初始面不能沿着流向摆放的原因。

形式陈述 ​

方程与非特征数据 ​

考虑标量一阶拟线性偏微分方程

(1)a(x,u(x))⋅Du(x)=b(x,u(x)),x∈Rn.

a 取值于 Rn,b 取实值,二者在所讨论点的开邻域内为 C1。初始面用 C1 参数化 x=σ(ξ) 给出,ξ∈Rn−1,Dσ 秩为 n−1;指定 C1 数据

u(σ(ξ))=g(ξ).

在 ξ0 处要求下列行列式非零

(2)det⁡[Dσ(ξ0), a(σ(ξ0),g(ξ0))]≠0.

矩阵前 n−1 列是初始面的切方向,最后一列是方程给出的流向。条件(2)称数据在该点非特征:最后一列不在前面各列张成的切平面中。流向可以依赖 u,所以这里不只检查初始面的形状,还要代入初值 g。

在这些条件下,σ(ξ0) 附近存在唯一的局部 C1 解,满足(1)和初始数据。局部唯一指任意两解在足够小的共同邻域内相同;不保证一张初始面能生成全空间、全时间的经典解。

要解的扩展ODE ​

对每个初始标签 ξ 解

(3)X˙=a(X,U),U˙=b(X,U),X(0,ξ)=σ(ξ),U(0,ξ)=g(ξ).

参数 s 表示沿曲线走了多远。X∈Rn 记录位置,U∈R 记录该位置携带的函数值。这是一套 n+1 维常微分方程,未知量暂时不是空间函数 u。

直觉

为什么沿曲线只剩一个导数 ​

若已经有经典解,沿满足 X˙=a(X,u(X)) 的位置曲线,由链式法则得到

ddsu(X(s))=Du(X(s))⋅a(X(s),u(X(s)))=b(X(s),u(X(s))).

于是 (X,u(X)) 必须满足(3)。这个推导告诉我们该解什么ODE,但还没有证明它的解能还原为一张函数图像。

从曲线恢复函数需要一个逆映射 ​

C1 系数在小紧邻域内局部Lipschitz,Picard–Lindelöf定理给出(3)的局部唯一解。ODE初值可微依赖又保证 X,U 对 (ξ,s) 为 C1。

令位置投影

Φ(ξ,s)=X(s,ξ).

在 (ξ0,0),它的导数恰是(2)中的矩阵。反函数定理因而给出局部坐标 (ξ(x),s(x))。定义

u(x)=U(s(x),ξ(x)).

它是 C1,在初始面给出 g,沿每条参数曲线又满足 Du⋅a=b。这完成存在性。

若 v 是另一局部经典解,从同一初始位置出发沿 X˙=a(X,v(X)) 行走,链式法则使 (X,v(X)) 满足同一套(3)。ODE唯一性使它与刚才的曲线及所携函数值一致。位置投影覆盖一个开邻域,故该邻域内 v=u,局部唯一性也得到证明。

例子与边界

变速输运还带有源项 ​

解

(4)ut+xux+u=x,u(0,x)=x2.

把 (t,x) 作为两个位置坐标,初始曲线为 (0,ξ);流向 (1,x) 不与它相切。取特征参数就是 t,得到

X′=X,X(0)=ξ,U′+U=X,U(0)=ξ2.

故 X=etξ。第二式乘积分因子 et 后积分,得到

U=e−tξ2+ξ2(et−e−t).

代入 ξ=e−tx,最后的空间函数为

(5)u(t,x)=e−3tx2+x2(1−e−2t).

直接计算 ut+xux+u:平方项系数为 −3+2+1=0,一次项合为 x,初值也等于 x2。这两项回代验证同时检查了沿流衰减与外部源的作用。

特征初始面可能不相容,也可能不唯一 ​

对 ux=0,在水平直线 σ(ξ)=(ξ,0) 上给初值。流向 (1,0) 正好沿着初始面,(2)失败。

若 g(ξ)=ξ,沿该线初值要求 ux=1,与方程矛盾,没有经典解。若 g 恒为常数 c,则任意 u(x,y)=c+ψ(y)、ψ∈C1 且 ψ(0)=0 都满足数据,解不唯一。非特征条件的失败不是一律“无解”;它使上述存在唯一性保证失效,实际情况取决于数据。

扩展曲线唯一,位置图像仍会折叠 ​

对 ut+uux=0、u(0,x)=−x,有

X(t,ξ)=(1−t)ξ,U(t,ξ)=−ξ,u(t,x)=−x1−t(t<1).

在 t=1,所有位置都投到 x=0,不同标签却仍携带不同的 U。扩展空间中的ODE曲线并没有相交;失去可逆性的是从曲线族到位置的投影。经典图像因此不能继续穿过这次坍缩。这里与已有Burgers守恒律相接,后续弱解和熵选择属于另一套接口。

推论与应用

一般输运的可复算公式 ​

对 ut+v(t,x)⋅Du+c(t,x)u=f(t,x),先解 X′=v(t,X)。沿轨道定义

C(t)=∫0tc(s,X(s))ds.

则

u(t,X(t))=e−C(t)[g(X(0))+∫0teC(s)f(s,X(s))ds].

公式先要求所需轨道确实存在,并能追到初始面;若问题有侧边界,还要判断轨道最先从哪一部分边界进入,不能给所有边界同时任意指定数据。

当方程对梯度非线性时,仅追踪位置和函数值通常不够。Hamilton–Jacobi特征构造再追踪梯度 P,用作用量恒等式证明它确实是还原函数的空间导数;位置投影是否可逆仍是关键检查。

参考资料
关系图谱14 个相邻概念 · 2 类关系

拖动节点调整位置。

显示关系

显示:依赖

  1. 前置三跳
  2. 前置二跳
  3. 前置一跳
  4. 当前条目
  5. 后续一跳
  6. 后续二跳
  7. 后续三跳
文字版关系按与当前条目的最短距离分组
类型化关系