2025-11-16T20:16:12.924612

On linear waves with memory in a Bessel-like medium

Giusti, Colombaro, Mentrelli
We discuss the propagation of harmonic and transient waves for systems governed by a wave equation with memory whose integral kernel involves ratios of modified Bessel functions of the first kind in the Laplace domain. In particular, the investigation of transient waves is carried out by means of a fully numerical approach based on the Talbot method for the numerical inversion of Laplace transforms.
academic

О линейных волнах с памятью в среде типа Бесселя

Основная информация

  • ID статьи: 2510.11493
  • Название: On linear waves with memory in a Bessel-like medium
  • Авторы: A. Giusti, I. Colombaro, A. Mentrelli
  • Классификация: math-ph math.MP
  • Дата публикации: 14 октября 2025
  • Ссылка на статью: https://arxiv.org/abs/2510.11493

Аннотация

В данной работе рассматривается распространение гармонических и переходных волн в системах, описываемых волновым уравнением с памятью, интегральное ядро которого в области преобразования Лапласа включает отношение модифицированных функций Бесселя первого рода. В частности, исследование переходных волн проводится с использованием полностью численного метода, основанного на методе Талбота для численного обращения преобразования Лапласа.

Исследовательский контекст и мотивация

Описание проблемы

Основная проблема, которую решает данное исследование, касается поведения распространения волн в среде типа Бесселя с эффектами памяти. Конкретно, объектом исследования является волновое уравнение с памятью вида:

∂2Y∂t2−c2[1−Φν(t)∗]∂2Y∂x2=0\frac{\partial^2 Y}{\partial t^2} - c^2\left[1 - \Phi_\nu(t) * \right]\frac{\partial^2 Y}{\partial x^2} = 0

где ядро памяти в области преобразования Лапласа представляется как: Φ~ν(s)=2(ν+1)sτIν+1(sτ)Iν(sτ)\tilde{\Phi}_\nu(s) = \frac{2(\nu+1)}{\sqrt{s\tau}}\frac{I_{\nu+1}(\sqrt{s\tau})}{I_\nu(\sqrt{s\tau})}

Значимость и приложения

Уравнения такого типа имеют важные приложения в нескольких областях физики:

  1. Гемодинамика: описание распространения волн в заполненных жидкостью упругих трубках
  2. Механика вязкоупругих материалов: моделирование эффектов памяти материалов
  3. Электрические системы: анализ электрических цепей с характеристиками памяти

Ограничения существующих методов

Предыдущие исследования в основном использовали полуаналитический алгоритм Бухена-Майнарди, однако этот метод надёжен только вблизи волнового фронта и не может обеспечить полную картину эволюции в пространстве-времени.

Исследовательская мотивация

Целью данной работы является преодоление ограничений существующих методов посредством численных подходов, обеспечивающих полный анализ поведения распространения волн во всей пространственно-временной области.

Основные вклады

  1. Установление полной структуры дисперсионного анализа: для случая ν=0 проведён детальный анализ дисперсионного соотношения гармонических волн с пространственным затуханием и временной периодичностью (SATP)
  2. Вывод точных дисперсионных законов и коэффициентов затухания: получены явные выражения через функции Кельвина
  3. Реализация полностью численного анализа переходных волн: применён метод Талбота для численного вычисления обратного преобразования Лапласа
  4. Выявление аномальной дисперсии: доказано, что система проявляет аномальную дисперсию, при которой фазовая скорость меньше групповой скорости

Подробное описание методов

Определение задачи

Исследование характеристик распространения волн в среде типа Бесселя с памятью включает:

  • Анализ дисперсии и затухания гармонических волн
  • Расчёт переходной характеристики волн при ступенчатом возбуждении
  • Определение фазовой и групповой скоростей

Архитектура модели

1. Основное волновое уравнение

Рассматривается частный случай при ν=0: ∂2Y∂t2−c2[1−Φ(t)∗]∂2Y∂x2=0\frac{\partial^2 Y}{\partial t^2} - c^2\left[1 - \Phi(t) * \right]\frac{\partial^2 Y}{\partial x^2} = 0

где: Φ~(s)=2sτI1(sτ)I0(sτ)\tilde{\Phi}(s) = \frac{2}{\sqrt{s\tau}}\frac{I_1(\sqrt{s\tau})}{I_0(\sqrt{s\tau})}

2. Метод дисперсионного анализа

Используется стандартная форма решения в виде гармонической волны: Y(t,x)=Aei(ωt−kx)Y(t,x) = A e^{i(\omega t - kx)}

Получается дисперсионное соотношение: (iω)2+k2c2[1−Φ^(ω)]=0(i\omega)^2 + k^2 c^2[1 - \hat{\Phi}(\omega)] = 0

3. Анализ волн SATP

Для волн с пространственным затуханием и временной периодичностью устанавливается:

  • ω∈R\omega \in \mathbb{R} (действительная частота)
  • k=κ−iδattk = \kappa - i\delta_{att} (комплексное волновое число)

Путём решения системы нелинейных уравнений:

\kappa^2 - \delta_{att}^2 = A(\omega) \\ -2\kappa\delta_{att} = B(\omega) \end{cases}$$ где A(ω) и B(ω) выражаются через функции Кельвина. ### Технические инновации #### 1. Представление через функции Кельвина Преобразование отношения комплексных функций Бесселя в комбинацию действительных функций Кельвина, что избегает численной неустойчивости при комплексных вычислениях. #### 2. Применение метода Талбота Использование метода Талбота для обратного преобразования Лапласа: $$Y(t,x) = \frac{1}{2\pi i}\int_C e^{st}\tilde{Y}(s,x)ds$$ Деформация контура интегрирования преобразует интеграл Бромвича в интеграл по действительной оси, повышая численную устойчивость. #### 3. Численное решение во всей области Преодоление ограничений традиционных методов, действительных только вблизи волнового фронта, и достижение точных вычислений во всей пространственно-временной области. ## Экспериментальная установка ### Начально-граничная задача Рассматривается задача о переходной характеристике на полубесконечной области: $$\begin{cases} Y_{tt} - c^2[1 - \Phi(t)*]Y_{xx} = 0, & x > 0, t > 0 \\ Y(t,0) = H(t), & t > 0 \\ Y(0,x) = 0, Y_t(0,x) = 0, & x > 0 \end{cases}$$ ### Численная реализация - Использование Wolfram Mathematica® для символических вычислений - Метод конечных разностей для численного дифференцирования групповой скорости - Высокоточное численное интегрирование методом Талбота ### Параметры Использование безразмерных параметров: - Параметр частоты: $\omega\tau$ - Пространственный параметр: $x/(c\tau)$ - Временной параметр: $(ct-x)/(c\tau)$ ## Результаты исследования ### Основные результаты #### 1. Дисперсионный закон и коэффициент затухания На рисунке 1 показаны безразмерные дисперсионный закон $c\tau\kappa$ и коэффициент затухания $c\tau\delta_{att}$ в зависимости от $\omega\tau$: - Дисперсионный закон демонстрирует нелинейный рост - Коэффициент затухания достигает максимума на средних частотах, затем убывает #### 2. Фазовая и групповая скорости На рисунке 2 представлено сравнение фазовой и групповой скоростей: - Фазовая скорость: $v_p(\omega) = \frac{\omega}{\kappa(\omega)}$ - Групповая скорость: $v_g(\omega) = \left[\frac{d\kappa(\omega)}{d\omega}\right]^{-1}$ - Численные результаты показывают: $v_p(\omega) < v_g(\omega) < c$, что подтверждает аномальную дисперсию #### 3. Распространение переходных волн На рисунке 3 показаны переходные характеристики в различных точках: - Волновой фронт распространяется со скоростью c - Волновая форма претерпевает значительную дисперсию при распространении - Области вдали от волнового фронта демонстрируют сложное колебательное затухание ### Важные открытия #### 1. Подтверждение аномальной дисперсии Система проявляет типичные характеристики аномальной дисперсии: - $v_p < v_g < c$ - Обе скорости стремятся к c в пределе высоких частот #### 2. Влияние эффектов памяти Член памяти приводит к: - Затуханию и дисперсии волн - Сложному переходному поведению - Неэкспоненциальной эволюции волновой формы ## Связанные работы ### Историческое развитие 1. **Mainardi и др. (1982-1983)**: первое введение подобных моделей в гемодинамику 2. **Giusti & Mainardi (2016)**: установление аналогии с динамической вязкоупругостью 3. **Colombaro и др. (2017)**: модели линейной вязкоупругости на основе функций Бесселя ### Сравнение методов - **Алгоритм Бухена-Майнарди**: полуаналитический, действителен только вблизи волнового фронта - **Метод данной работы**: полностью численный, действителен во всей области - **Асимптотические методы**: высокочастотные приближения, сильные ограничения ## Заключение и обсуждение ### Основные выводы 1. **Дисперсионные характеристики**: система проявляет аномальную дисперсию, при которой фазовая скорость всегда меньше групповой скорости, и обе не превышают скорость волнового фронта c 2. **Эффективность численного метода**: метод Талбота успешно преодолевает ограничения традиционных методов 3. **Эффекты памяти**: ядро памяти типа Бесселя приводит к сложному поведению распространения волн ### Ограничения 1. **Вычислительная сложность**: переходный этап требует большого количества узлов, что приводит к высоким вычислительным затратам 2. **Ограничения параметров**: рассмотрен только частный случай ν=0 3. **Недостаточность теоретического анализа**: строгие свойства групповой скорости трудно доказать аналитически ### Направления будущих исследований 1. **Гибридные методы**: комбинирование полуаналитических методов для обработки волнового фронта с численными методами для дальнего поля 2. **Обобщение**: расширение на случай произвольных значений ν 3. **Расширение приложений**: проверка применимости в конкретных физических системах ## Глубокая оценка ### Преимущества 1. **Методологическая инновация**: первое применение метода Талбота к задаче распространения волн в среде с памятью типа Бесселя 2. **Теоретическая полнота**: предоставление полной структуры от дисперсионного анализа до расчёта переходных процессов 3. **Численная точность**: преодоление ограничений существующих методов, достижение высокоточных вычислений во всей области 4. **Физические инсайты**: выявление важных физических явлений, таких как аномальная дисперсия ### Недостатки 1. **Вычислительная эффективность**: высокая сложность расчётов переходных процессов ограничивает практические приложения 2. **Диапазон параметров**: ограничение конкретными формами ядра памяти 3. **Экспериментальная верификация**: отсутствие сравнения с результатами реальных физических экспериментов ### Влияние 1. **Теоретический вклад**: предоставление новых инструментов численного анализа для распространения волн в средах с памятью 2. **Прикладная ценность**: потенциальные приложения в биомедицине, материаловедении и других областях 3. **Методологическое значение**: применение метода Талбота служит ориентиром для решения аналогичных задач ### Области применения 1. **Академические исследования**: теория распространения волн, моделирование материалов с памятью 2. **Инженерные приложения**: анализ потока в трубопроводах, управление вибрациями 3. **Численные вычисления**: высокоточные расчёты обратного преобразования Лапласа ## Библиография Статья ссылается на 23 важные работы, охватывающие: - Теорию функций Бесселя (Abramowitz & Stegun, 1972) - Теорию распространения волн (Mainardi, 2022; Whitham, 2011) - Численные методы (Talbot, 1979; Garrappa и др., 2021) - Прикладной контекст (литература по гемодинамике, вязкоупругости и связанным темам) --- Данная работа вносит значительный вклад как в теоретический анализ, так и в численные методы, предоставляя новые инструменты и физические инсайты для исследования распространения волн в средах с памятью. Несмотря на определённые ограничения в вычислительной эффективности и диапазоне параметров, её методологическая ценность и физические открытия делают её важным прогрессом в данной области.