关于间断对流方程（DGFEM）源项推导及弱解唯一性的技术疑问

咱们先理清楚你碰到的核心矛盾——你混淆了分布意义下的精确解和满足物理约束的弱解的区别，还有间断系数对流方程弱解的特殊性，一步步拆解：

1. 习题要求的“exact solution”到底是什么？ #

题目说的“exact solution”指的是在分布意义下严格满足PDE的解（允许包含δ函数型的源项），而不是你以为的“$g=0$时的弱解”。

你可以直接把给定的$u(x,t)=\sin(\pi(x-a(x)t))$代入PDE计算$g(x,t)$：

• 在$|x|≠0.5$的光滑区域，$a(x)$是常数，求导后显然$\partial_t u +a(x)\partial_x u=0$；
• 在$x=±0.5$的间断处，$a(x)$和$u(x,t)$都有跳变，这时候用分布导数计算会得到δ函数型的项——这部分就是$g(x,t)$的核心组成。

具体来说，利用分布导数的性质：对于在$x=\xi$处有跳变的函数，其导数会包含跳变幅度乘以δ函数。计算后你会发现，$g(x,t)$是由两个δ函数项组成的，对应$x=±0.5$处的跳变贡献，光滑区域则为0。

2. 为什么特征线法得到的解和题目给的$u$不一样？ #

你这里的关键误解是：题目给的$u$根本不是$g=0$时的弱解，它不满足间断对流方程弱解必须遵守的Rankine-Hugoniot条件（跳变条件）。

对于静止间断面$x=±0.5$，$g=0$的对流方程弱解需要满足：$a^+u+ = a^-u-$（$a^+/a-$是间断面两侧的$a(x)$值，$u^+/u-$是$u$的两侧极限）。但题目给的$u$在$x=0.5$处：
$a^-u-=1.5\sin(\pi(0.5-1.5t))$，$a^+u+=1\cdot\sin(\pi(0.5-t))$，显然两者不相等，所以这个$u$不可能是$g=0$的弱解——你之前觉得它满足弱解，是因为拆分积分时忽略了间断面处的分部积分边界项，那些项不为零，必须由$g(x,t)$来抵消。

特征线法得到的解是满足熵条件的物理弱解，它严格遵循特征线的传播规律，在间断面处满足Rankine-Hugoniot条件，和题目给的$u$对应不同的PDE（一个是$g=0$，一个是带δ源项的$g(x,t)$），自然不一样。

3. 弱解的唯一性问题 #

对于带间断系数的对流方程，无附加条件的弱解确实不唯一——你可以构造多个满足弱形式的解，但只有满足熵条件（或类似物理约束）的解才是有实际意义的（比如特征线法得到的解）。

题目给的$u$是另一个弱解，但它对应带源项的PDE，并不是$g=0$的解，这也解释了为什么你用DGFEM求解$g=0$的PDE时，得到的是特征线型的解，而非题目给的$u$。

总结 #

习题本身没问题，只是你误解了“exact solution”的含义——它不是$g=0$时的物理弱解，而是需要你推导对应的$g(x,t)$（包含δ函数项），使得$u$在分布意义下满足PDE。而弱解的唯一性需要附加熵条件才能保证，这也是特征线解和题目给的$u$不同的核心原因。

备注：内容来源于stack exchange，提问作者Victor Hansen