[论文解读] Analytic solutions for the Burgers equation with source terms

Gino I. Montecinos|arXiv (Cornell University)|Mar 29, 2015

Nonlinear Waves and Solitons参考文献 3被引用 11

一句话总结

本文提出了一种基于特征线的解析方法，用于求解带有源项的Burgers方程，包括刚性及非线性形式。通过将偏微分方程（PDE）转化为两个可解的常微分方程（ODE）——等价ODE（控制沿特征线的演化）和特征线ODE（定义时空中的曲线），当源项的倒数具有原函数时，该方法可获得精确解。关键贡献在于提出了一套系统化程序，用于生成解析解以验证数值格式，即使在求解特征线曲线时需通过二分法进行数值求根亦可适用。

ABSTRACT

Analytic solutions for Burgers equations with source terms, possibly stiff, represent an important element to assess numerical schemes. Here we present a procedure, based on the characteristic technique to obtain analytic solutions for these equations with smooth initial conditions.

研究动机与目标

• 开发一种系统化方法，推导带有源项（包括刚性和非线性形式）的Burgers方程的精确解析解。
• 通过提供能捕捉非线性和刚性的精确解，实现对数值格式的严格验证。
• 建立基于特征线的ODE系统可解的条件，特别关注源项及其倒数是否存在原函数。
• 通过引入二分法求解特征线曲线所导致的代数方程，解决解对初值的隐式依赖问题。

提出的方法

• 该方法采用特征线法，将PDE约化为两个常微分方程（ODE）组成的系统。
• 沿特征线由PDE导出的等价ODE为 $ \frac{d}{dt}\hat{q}(t) = s(\hat{q}(t)) $，初值条件为 $ \hat{q}(0) = h_0(y) $。
• 若 $ s(q)^{-1} $ 具有原函数 $ G(q) $，则等价ODE的可解性得以保证，解可表示为 $ \mathcal{E}(t, h_0(y)) = G^{-1}(t + G(h_0(y))) $。
• 特征线ODE $ \frac{d}{dt}x(t) = \mathcal{E}(t, h_0(y)) $ 通过积分 $ \mathcal{E} $ 求解，得到 $ x = y + \mathcal{F}(t, h_0(y)) $，其中 $ \mathcal{F} $ 是 $ \mathcal{E} $ 的原函数。
• 完整解通过隐式关系定义：$ q(x,t) = \mathcal{E}(t, h_0(y)) $，其中 $ y $ 满足 $ x = y + \mathcal{F}(t, h_0(y)) $，在非显式情况下需借助数值根求解。
• 当 $ y $ 无法显式求解时，采用二分法求解由特征线曲线定义所导出的代数方程。

实验结果

研究问题

• RQ1 能否基于特征线方法，为带有刚性或非线性源项的Burgers方程推导出解析解？
• RQ2 在何种源项条件下，等价ODE是可解的？
• RQ3 当特征线方程为隐式时，解如何用初值表示？
• RQ4 在特征线公式中出现隐式代数方程时，哪种数值方法能有效求解？

主要发现

• 若源项 $ s(q)^{-1} $ 具有原函数 $ G(q) $，则等价ODE可解，解为 $ \mathcal{E}(t, h_0(y)) = G^{-1}(t + G(h_0(y))) $。
• 若解 $ \mathcal{E}(t, h_0(y)) $ 具有原函数 $ \mathcal{F}(t, h_0(y)) $，则特征线ODE可解，得到曲线 $ x = y + \mathcal{F}(t, h_0(y)) $。
• 对于线性源项 $ s(q) = \beta q $，解为 $ \hat{q}(t) = \frac{h_0(y)}{1 - \beta t h_0(y)} $，特征线为 $ x = y - \frac{\ln(1 - \beta t h_0(y))}{\beta} $。
• 对于二次源项 $ s(q) = \beta q^2 $，解为 $ \hat{q}(t) = \frac{h_0(y)}{1 - \beta t h_0(y)} $，且特征线可显式求解。
• 对于非线性源项 $ s(q) = e^{\beta q} $，解为 $ \hat{q}(t) = -\frac{1}{\beta} \ln(e^{-\beta h_0(y)} - \beta t) $，特征线通过一个复杂的原函数表达。
• 当特征线方程对 $ y $ 无法显式求解时，采用二分法数值求解 $ y(x,t) $，从而可在任意 $ (x,t) $ 处计算解的值。