[论文解读] Fractional Laplacian, Levy stable distribution, and time-space models for linear and nonlinear frequency-dependent lossy media

W. Chen, Sverre Holm|ArXiv.org|Dec 30, 2002

Fractional Differential Equations Solutions参考文献 36被引用 13

一句话总结

本文提出了一种基于分数阶拉普拉斯算子的波动模型，用于描述损耗介质中频率相关的衰减，其中衰减遵循指数在0到2之间的幂律。通过引入一种具有内在正则化和边界处理能力的新分数阶拉普拉斯定义，该模型能够准确捕捉线性和非线性波动方程中的异常（非整数幂律）衰减，并将其与Lévy稳定分布和分形微观结构联系起来。

ABSTRACT

The frequency-dependent attenuation typically obeys an empirical power law with an exponent ranging from 0 to 2. The standard time-domain partial differential equation models can describe merely two extreme cases of frequency independent and frequency-squared dependent attenuations. The otherwise non-zeor and non-square frequency dependency occuring in many cases of practical interest is thus often called the anomalous attenuation. In this study, we developed a linear integro-differential equation wave model for the anomalous attenuation by using the space fractional Laplacian operation, and the strategy is then extended to the nonlinear Burgers, KZK, and Westervelt equations. A new definition of the fractional Laplacian is also introduced which naturally includes the boundary conditions and has inherent regularization to ease the hyper-singularity in the conventional fractional Laplacian. Under the Szabo's smallness approximation where attenuation is assumed to be much smaller than the wave number, our linear model is found consistent with arbitrary frequency dependencies. According to the fact that the physical attenuation can be understood a statistic process, the empirical range [0,2] of the power law exponent is explained via the Levy stable distribution theory. It is noted that the power law attentuation underlies fractal microstructures of anomalously attenuative media.

研究动机与目标

• 解决标准波动模型仅能描述频率无关或与频率平方成正比衰减的局限性。
• 利用分数阶微积分对实际介质中常见但经典偏微分方程无法捕捉的异常衰减进行建模。
• 为频率相关损耗介质中的线性和非线性波动传播提供物理解释一致的时间-空间模型。
• 通过Lévy稳定分布理论解释经验观测到的幂律指数范围[0,2]的物理成因。
• 建立异常衰减与介质中分形微观结构之间的联系。

提出的方法

• 提出一种新的分数阶拉普拉斯定义，其天然地包含边界条件并能对超奇点进行正则化。
• 通过空间分数阶拉普拉斯算子构建线性积分微分波动方程，以实现对任意频率相关衰减的建模。
• 通过分数阶算子将线性模型推广至非线性方程，包括Burgers方程、KZK方程和Westervelt方程。
• 应用Szabo的小量近似方法，确保时间域中对任意幂律衰减的一致性。
• 利用Lévy稳定分布理论解释幂律指数范围[0,2]的物理起源。
• 使用分数阶微积分推导时间-空间域中的波动模型，从而实现对衰减机制的直接物理解释。

实验结果

研究问题

• RQ1 如何利用分数阶微积分对介于标准频率无关和频率平方衰减之间的频率相关衰减进行建模？
• RQ2 在损耗介质中经验观测到的幂律衰减指数范围[0,2]的物理基础是什么？
• RQ3 如何重新定义分数阶拉普拉斯算子，使其自然地包含边界条件并减少数学奇点？
• RQ4 线性和非线性波动方程（Burgers、KZK、Westervelt）在多大程度上可以通过分数阶算子进行推广？
• RQ5 异常衰减与介质的分形微观结构之间存在何种关系？

主要发现

• 所提出的分数阶拉普拉斯模型成功捕捉了[0,2]范围内任意幂律衰减，突破了经典偏微分方程模型的局限性。
• 新的分数阶拉普拉斯定义消除了对外部边界处理的需求，并减少了超奇点，从而提高了数值稳定性。
• 该模型与Szabo的小量近似一致，确保在衰减水平远小于波数时仍保持有效性。
• 经验观测的幂律指数范围[0,2]通过Lévy稳定分布理论得到物理解释，将统计力学与波动传播联系起来。
• 研究表明，异常衰减源于介质中的分形微观结构，为观测到的幂律行为提供了物理解释。
• 该模型向非线性方程（Burgers、KZK、Westervelt）的推广表明其在损耗介质中真实非线性波动现象中的广泛应用潜力。