[论文解读] Going with the Flow: a Lagrangian approach to self-similar dynamics and its consequences

D. G. Aronson, S. Betelú+1 more|ArXiv.org|Nov 23, 2001

Advanced Mathematical Modeling in Engineering参考文献 14被引用 14

一句话总结

本文提出了一种拉格朗日动态钉扎方法——称为'MN-动力学'——通过基于模板的变换，系统地从非线性偏微分方程中提取自相似解。该方法将原始演化方程转化为一个在自相似参考系中随时间演化的非局部方程，从而能够同时计算第一类（Barenblatt）和第二类（Graveleau）解，并自动获得自相似指数。

ABSTRACT

We present a systematic computational approach to the study of self-similar dynamics. The approach, through the use of what can be thought of as a ``dynamic pinning condition" factors out self-similarity, and yields a transformed, non-local evolution equation. The approach, which is capable of treating both first and second kind self-similar solutions, yields as a byproduct the self-similarity exponents of the original dynamics. We illustrate the approach through the porous medium equation, showing how both the Barenblatt (first kind) and the Graveleau (second kind) self-similar solutions arise in this framework. We also discuss certain implications of the dynamics of the transformed equation (which we will name "MN-dynamics"); in particular we discuss the discrete-time implementation of the approach, and connections with time-stepper based methods for the "coarse" integration/bifurcation analysis of microscopic simulators.

研究动机与目标

• 开发一种系统化的计算方法，用于识别非线性演化方程中的自相似解。
• 将先前用于行波的基于模板的动态钉扎方法，扩展至自相似动力学。
• 在连续时间内动态地消除自相似性，得到一个揭示趋近自相似态过程的变换演化方程。
• 在统一框架内实现第一类（如Barenblatt）和第二类（如Graveleau）自相似解的计算。
• 通过避免显式推导宏观方程，支持基于时间步长的微观模拟器的粗粒度分岔分析。

提出的方法

• 应用基于模板的动态钉扎条件以消除自相似性，将原始PDE转化为自相似参考系中的非局部演化方程。
• 使用缩放变换定义一个新变量，以消除自相似依赖性，从而得到'MN-动力学'方程。
• 通过数值求解变换后的方程，提取自相似解及其关联的指数。
• 在连续时间和离散时间两种框架下实现该方法，使其适用于基于时间步长的微观动力学模拟。
• 利用变换后的动力学，应用动力系统工具（如分岔理论、延续法）于自相似解。
• 利用离散时间版本分析微观模拟器（如蒙特卡洛、格子玻尔兹曼）中的自相似性，而无需推导粗粒度方程。

实验结果

研究问题

• RQ1 如何在连续时间内动态地消除自相似性，类似于用于行波的动态钉扎方法？
• RQ2 该基于模板的方法能否从平移不变性推广至自相似性，从而得到一个非局部演化方程？
• RQ3 第一类和第二类自相似解如何自然地从变换后的动力学中涌现？
• RQ4 将自相似解视为MN-动力学的稳态时，对数值分岔分析有何影响？
• RQ5 该方法能否应用于微观模拟器（如MD、MC）以检测自相似行为，而无需闭式宏观方程？

主要发现

• MN-动力学框架成功计算了多孔介质方程的第一类（Barenblatt）和第二类（Graveleau）自相似解。
• 该方法在求解过程中自动获得自相似指数作为副产品。
• 变换后的方程使得高级动力系统技术（如延续法、分岔理论）可应用于自相似解。
• 离散时间实现使得无需显式构建粗粒度MN方程，即可对微观时间步长模拟器进行有效的分岔分析。
• 该方法为在宏观方程不可得或过于复杂而难以推导的系统中，实现自相似动力学的计算机辅助分析开辟了道路。
• 该方法与多尺度积分器（如投影积分器）兼容，可用于为微观模拟定义物理信息引导的系综。