[论文解读] The Riccati Differential Equation and a Diffusion-Type Equation

Erwin Suazo, Sergeĭ K. Suslov+1 more|ArXiv.org|Jul 28, 2008

Differential Equations and Numerical Methods被引用 19

一句话总结

本文为一类一维扩散型方程在时变系数下的柯西初值问题提供了显式解析解。通过将问题约化为里卡蒂方程及其关联的二阶线性特征方程，作者推导出格林函数（热核）的闭式表达式，该表达式由初等函数和积分构成，从而可借助杜哈梅尔原理对齐次与非齐次情形实现精确求解。

ABSTRACT

We construct an explicit solution of the Cauchy initial value problem for certain diffusion-type equations with variable coefficients on the entire real line. The corresponding Green function (heat kernel) is given in terms of elementary functions and certain integrals involving a characteristic function, which should be found as an analytic or numerical solution of the second order linear differential equation with time-dependent coefficients. Some special and limiting cases are outlined. Solution of the corresponding non-homogeneous equation is also found.

研究动机与目标

• 为实直线上具有变系数的一般扩散型方程的柯西初值问题构造显式解。
• 确定格林函数（热核）可表示为初等函数与积分的条件。
• 建立扩散方程解与由特征函数导出的时变系数里卡蒂方程解之间的联系。
• 通过杜哈梅尔原理将该方法推广至非齐次方程。
• 提供特殊情形的显式解，包括具有双曲与三角函数形式时变系数的方程。

提出的方法

• 通过假设形式 $ u = A(t) e^{S(x,y,t)} $ 构造解，其中 $ A(t) = 1 / \sqrt{2\pi\mu(t)} $，且 $ S $ 是关于 $ x $ 和 $ y $ 的二次型，其系数为时变函数。
• 将该假设代入PDE后，通过求解所得常微分方程组确定 $ S $ 中的系数，从而导出 $ \alpha(t) $ 满足的里卡蒂方程。
• 将里卡蒂方程转化为关于特征函数 $ \mu(t) $ 的二阶线性常微分方程，形式为 $ \mu'' - \tau(t)\mu' - 4\sigma(t)\mu = 0 $，其中 $ \tau $ 与 $ \sigma $ 由时变系数定义。
• 格林函数显式给出为 $ K(x,y,t) = \frac{1}{\sqrt{2\pi\mu(t)}} \exp(S(x,y,t)) $，其中 $ S $ 依赖于该常微分方程组的解。
• 该方法依赖于求解特征方程 $ \mu'' - \tau(t)\mu' - 4\sigma(t)\mu = 0 $，初始条件为 $ \mu(0) = 0 $，$ \mu'(0) = 2a(0) \neq 0 $，以确保正确的初始行为。
• 对于非齐次方程，应用杜哈梅尔原理，将解表示为包含齐次解算子与源项的时间积分。

实验结果

研究问题

• RQ1 在何种条件下，具有时变系数的扩散型PDE的基本解可使用初等函数与积分以闭式表达？
• RQ2 扩散方程的解与具有时变系数的里卡蒂微分方程的解之间有何关系？
• RQ3 特征函数 $ \mu(t) $ 的作用是什么？它如何决定格林函数的结构？
• RQ4 如何应用杜哈梅尔原理推导具有变系数的非齐次扩散方程的解？
• RQ5 对于特定的时变系数函数（如双曲与三角函数形式），热核的显式形式是什么？

主要发现

• 扩散方程的格林函数被显式构造为 $ K(x,y,t) = \frac{1}{\sqrt{2\pi\mu(t)}} \exp(\alpha x^2 + \beta xy + \gamma y^2 + \delta x + \varepsilon y + \kappa) $，其中系数由求解常微分方程组确定。
• 系数 $ \alpha(t) $ 满足里卡蒂方程：$ \frac{d\alpha}{dt} + b(t) - 2c(t)\alpha - 4a(t)\alpha^2 = 0 $，该方程是解法的核心。
• 特征函数 $ \mu(t) $ 满足二阶线性常微分方程 $ \mu'' - \tau(t)\mu' - 4\sigma(t)\mu = 0 $，其中 $ \tau(t) = \frac{a'}{a} + 2c - 4d $，$ \sigma(t) = ab + cd - d^2 + \frac{d}{2}\left(\frac{a'}{a} - \frac{d'}{d}\right) $。
• 当 $ a(t) = \sinh^2 t $，$ b(t) = \cosh^2 t $，$ c(t) = \frac{1}{2}\sinh 2t $ 时，特征方程变为 $ \mu'' - 2\tanh t\, \mu' + 2\mu = 0 $，其解为 $ \mu_1(t) = \cos t \sinh t + \sin t \cosh t $，在 $ 0 < t < T_1 \approx 0.93755 $ 范围内有效。
• 该情形下的热核显式给出，其收敛性条件与 $ \gamma(t) $ 的符号相关，确保对合适初始数据的可积性。
• 对于方程 $ a(t) = \cos^2 t $，$ b(t) = \sin^2 t $，$ c(t) = -\frac{1}{2}\sin 2t $，也导出了类似解，热核形式相同但 $ x $ 与 $ y $ 互换，有效范围为 $ 0 < t < T_2 \approx 2.347 $。