[论文解读] Numerical modelling of convection-diffusion-reaction problems with free boundary in 1D

Gabriela Kacurova|ArXiv.org|Sep 2, 2009

Fluid Dynamics and Thin Films参考文献 7被引用 3

一句话总结

本文提出了一种用于一维对流-扩散-反应问题的新型数值方法，其自由边界通过从局部行波解导出的常微分方程（ODE）建模界面演化。通过将问题转换到固定域，并在界面附近采用自适应空间离散化，该方法能够高精度捕捉尖锐前缘，其有效性通过与解析解以及多孔介质和反应-扩散方程的先前数值研究结果的对比得到验证。

ABSTRACT

We discuss a numerical method for convection-diffusion-reaction problems with a free boundary in 1D. The method is based on the numerical modelling of the interface evolution, the transformation to a fixed domain problem and the approximation by an ODE system. The interface evolution is modelled by means of the local shape of the corresponding travelling wave solution. The method can be applied to many free boundary problems with a finite speed of the interface. The presented method can also approximate some problems with an infinite speed of the interface for damped travelling wave type solutions. In the numerical experiments we compare our numerical solution with the analytical ones for some problems.

研究动机与目标

• 开发一种高效且准确的一维对流-扩散-反应自由边界问题数值方法，特别适用于具有有限速度界面的问题。
• 通过从行波解中导出的局部ODE建模界面随时间的演化 $ s(t) $，避免使用昂贵的自适应网格技术。
• 通过 $ y = x/s(t) $ 的变换将移动边界问题转化为固定域问题，从而实现稳定且高精度的空间离散化。
• 将该方法应用于退化扩散、多孔介质和反应-扩散方程，包括非退化和阻尼行波解的情形。
• 通过与解析解及基准问题（如多孔介质方程和湍流流模型）的现有数值结果对比，验证该方法的有效性。

提出的方法

• 通过从行波解的局部形状导出的ODE建模界面演化 $ s(t) $，特别地，对于多孔介质型方程，使用 $ \dot{s}(t) = -\frac{n}{n-1} \partial_x u^{n-1} \big|_{x \nearrow s(t)} $。
• 通过变量替换 $ y = x/s(t) $，将移动域 $ (0, s(t)) $ 转换为固定参考域 $ (0, 1) $，从而实现标准空间离散化。
• 在 $ y = 1 $ 附近（即 $ x = s(t) $）应用专门的空间离散化方法，以高分辨率准确捕捉尖锐前缘。
• 将变换后的PDE离散化为时间方向的ODE系统，包括界面ODE，形成一个刚性ODE系统，可由鲁棒的数值求解器求解。
• 通过扰动非线性项中 $ p \nearrow 1 $ 的方式（如 $ \Psi(s^p) $）近似非退化情形（例如线性吸附），避免界面形成的同时保持精度。
• 通过将数值解和界面与解析解（如Barenblatt-Pattle解）及已发表的数值基准进行对比，验证结果。

实验结果

研究问题

• RQ1 是否可以使用从行波解局部形状导出的局部ODE，准确建模一维对流-扩散-反应问题中的界面演化？
• RQ2 通过 $ y = x/s(t) $ 将移动边界问题转换为固定域，是否能够实现稳定且高精度的数值求解，并在界面附近实现自适应空间分辨率？
• RQ3 该方法是否能准确近似退化（有限速度界面）和非退化（无限速度、阻尼波）情形，包括吸附等温线？
• RQ4 与现有方法相比，该方法在基准问题（如多孔介质方程和湍流流模型）上的精度和效率如何？
• RQ5 是否可以通过 $ p \nearrow 1 $ 的非线性扰动扩展该方法以近似非退化问题，且该近似是否稳定？

主要发现

• 该方法成功以高精度捕捉了解的尖锐前缘，多孔介质情形（$ n=6 $）中数值界面与解析界面近乎完美重合。
• 对于多孔介质方程，数值界面 $ s(t) $ 与解析解 $ s(t) = [\frac{2n(n+1)}{n-1}(t+1)]^{1/(n+1)} $ 的误差可忽略不计，证实了界面ODE模型的有效性。
• 在非退化扩散情形下，$ p=0.95 $ 时，该方法对线性吸附提供了稳定且精确的近似，结果与 [10] 和 [24] 的结果高度一致。
• 对于湍流流模型 $ \partial_t u = \partial_x^2 u^{3/2} - (u^{3/2} - u^{1/2}) $，数值界面与解析界面在视觉上无法区分，且误差随时间保持较低水平。
• 该方法准确模拟了泡沫排液模型 $ \partial_t u = \partial_x^2 u^{3/2} + \partial_x(u^2) $ 和黏性液体推进模型 $ \partial_t u = \partial_x^2 u^4 + \partial_x(u^3) $，重现了对称扩展和前缘传播等预期物理行为。
• 与 [24] 的算子分裂结果对比显示，在 $ D=0.1 $ 和 $ D=0.001 $ 情况下均表现出极佳的一致性，证实了该方法在不同扩散系数下的鲁棒性。