常系数传输方程

常系数传输方程（Constant Coefficient Transport Equation）。这是一个最基本、最重要的一类偏微分方程，它清晰地展示了波动和传输现象的本质。

定义与方程形式

常系数传输方程的标准形式是：

$$ \frac{\partial u}{\partial t} + c\frac{\partial u}{\partial x} = 0 $$

其中：$u(x,t)$是未知函数（例如流体的密度、浓度、温度等）。$c$ 是一个常数（实数），代表传输速度或波速。$\frac{\partial u}{\partial t}$是时间变化率，$\frac{\partial u}{\partial x}$是空间变化率

方程的含义：物理量$u$ 随时间的变化，完全是由其以恒定速度$c$ 在空间中的“流动”或“传输”引起的。

有时也会遇到更一般形式的源项方程：

$$ \frac{\partial u}{\partial t} + c\frac{\partial u}{\partial x} = f(x,t)$$。但核心和特性主要研究齐次方程($f=0$)

物理意义

这个方程描述了一种无扩散、无耗散、形状保持的平移运动。几个典型例子：

河流中的污染团：如果一股污染物以恒定流速 $c$随河水流向下游，且不考虑扩散，那么其浓度分布 $u (x,t)$就满足此方程。
交通流模型（简化）：在一段均匀道路上，车流以恒定速度 $c$行驶，车流密度 $u(x,t)$的演化。
波动现象（单方向）：比如一根非常长的弦，只向一个方向传播一个波形。
颜色传输：：在传送带上，某种颜色的染料图案以恒定速度向前移动。

核心思想：初始时刻$u$的“形状”或“剖面” $u(x,0)$，在运动过程中不发生任何形变，只是单纯地以速度 $c$向某个方向“平移”。

求解方法

特征线法

1.特征方程：将PDE转化为沿着一族特殊曲线的ODE。

$$ \frac{dx}{dt} = c$$解为：

$$ x-ct=\xi (常数)$$

这族直线$x=ct+ξ$ 称为特征线。参数 $ξ$ 标志了不同的特征线。

2.沿特征线化简：沿着任何一条特征线$x−ct=\xi$，原PDE变为：

$$ \frac{du}{dt} = \frac{\partial u}{\partial t} + \frac{\partial u}{\partial x} \cdot \frac{d x}{d t} = 0$$

这是因为,考虑函数 $u(x,t)$ 沿着其中一条特征线 $x(t)$（满足 $dx/dt=c$）的变化率。由于$x$本身又是 $t$的函数 $x=x(t)$，那么$u$最终可以看作只是$t$的复合函数：$u(x(t),t)$。根据多元函数的链式法则，这个复合函数对 $t$的导数就是全导数：

$$ \frac{du}{dt} = \frac{\partial u}{\partial t}\cdot \frac{dt}{dt} + \frac{\partial u}{\partial x}\cdot \frac{dx}{dt} = \frac{\partial u}{\partial t} + c \frac{\partial u}{\partial x} \implies \frac{du}{dt} = u_t + cu_x = 0 $$

$du/dt=0$ 意味着函数 $u$ 关于变量 $t$的变化率为零。一个量如果随时间的变化率为零，那么它不随时间改变，即它是一个常数。

3.得到通解:由于在特征线上$u$为常数，且该常数由该特征线($x=ct+\xi$)的起点（如$t=0$）决定，因此解可以写为：

$$ u(x,t) = u(\xi, 0) = u(x-ct,0)$$

这就是著名的行波解。

解的形式与性质

通解

$$ u(x,t) = f(x-ct)$$ 这里 $f$是一个任意函数，其具体形式由初始条件决定。

给定初始条件: 若已知初始时刻的分布$u(x,0)=\phi (x)$，则对比可得$f(x)=\phi(x)$。因此，初值问题（柯西问题）的解为：

$$ u(x,t) = \phi(x-ct)$$

$$ $$

在QkCalc中，常系数传输方程的输入模板为：

$$ \LARGE \begin{cases} \frac{\partial^2\color{red}u}{\partial\color{#08c365}t} + c \frac{\partial^2\color{red}u}{\partial\color{#08c365}x} = 0 \\ {\color{red}u}({\color{red}x,t}) = \Box \end{cases}$$

其中$u$是未知函数的符号。$t,x$是自变量。如果$u(x,t)=$右侧表达式为空，则表明求解的结果使用通解。例如上述表达式的通解是：

$$ u(x,t) = F(-tc + x) $$

有初值的情况：

非齐次常系数传输方程

非齐次常系数传输方程引入了源项（或强迫项），使得物理模型更加丰富（如考虑外部热源、化学反应、外力等）。我们来看如何用特征线法求解。方程形式：

$$\frac{\partial u}{\partial t} + c \frac{\partial u}{\partial x} = f(x,t)， \quad x \in \mathbb{R}, t > 0 $$

其中 $f(x,t)$ 是已知的源函数， $c$ 是常数

求解步骤（特征线法）

第一步：写出特征方程与齐次情况相同：

$$ \frac{dx}{dt} = c $$

第二步：求解特征线

$$ x = ct+\xi$$其中 $\xi$ 是常数，标记不同的特征线。

第三步：沿特征线将PDE转化为ODE沿特征线 $x=ct+\xi$，计算 $u$的全导数：

$$ \frac{du}{dt} = \frac{\partial u}{\partial t} + \frac{\partial u}{\partial x} \cdot \frac{dx}{dt} = u_t + cu_x = f(x,t) $$

关键区别：这里沿特征线的ODE不是$du/dt=0$，而是等于源项$f$。

第四步：将ODE中的 $x,t$用特征线参数表示

在特征线 $x=ct+\xi$ 上， $x$和 $t$不是独立的。我们可以将源项 $f(x,t)$表示为沿这条线的函数：

$$ x=ct+\xi$$，所以：

$$ f(x,t) = f(ct+\xi,t)$$

记：$\tilde{f}(\xi, t) = f(ct+\xi,t)$，则ODE变为：

$$ \frac{du}{dt} = \tilde{f}(\xi, t)$$

第五步：积分ODE

沿固定 $\xi$ 的特征线，将上述ODE从初始时刻 $t=0$ 积分到当前时刻 $t$：

$$ u(t) - u(0) = \int_{0}^t \tilde{f}(\xi, s)ds $$

$u(t)$ 表示在特征线上，对应于时间 $s=t$ 时的 $u$ 值，即 $u(x,t)$，其中 $x=ct+\xi$。
$u(0)$ 表示在特征线上，对应于时间 $s=0$ 时的 $u$ 值，即 $u(\xi,0)=\phi(\xi)$。

所以：

$$ u(x,t) = \phi(\xi) + \int_{0}^t f(cs + \xi, s)ds$$

第六步：将$\xi$用 $x,t$ 表示

由特征线方程$\xi=x−ct$，代入上式：

$$ u(x,t)=\phi(x−ct)+\int_0^t f(x−c(t−s),s)ds. $$

这是最终的求解公式。

在QkCalc中，非齐次的计算模板：

$$ \LARGE \begin{cases} \frac{\partial^2\color{red}u}{\partial\color{#08c365}t} + c \frac{\partial^2\color{red}u}{\partial\color{#08c365}x} = e^{-t} \\ {\color{red}u}({\color{red}x,t}) = 0\end{cases}$$

上面的方程通解是：

$$ u(x,t) = \phi(-tc+x)+1-e^{-t}$$

其中$\phi$表示初始波形

常系数非齐次，且有初始函数：

在QkCalc中它显示的解是：

$$ u(x,t)=\sin(-2t+x)+1-e^{-t}$$