2025-11-18T01:22:13.423549

Existence and numerical approximation of a one-dimensional Boussinesq system with variable coefficients on a finite interval

Grajales, Pizo
In this paper, we investigate the well-posedness of a nonlinear dispersive model with variable coefficients that describes the evolution of surface waves propagating through a one-dimensional shallow water channel of finite length with irregular bottom topography. To complement the theoretical analysis, we utilize the numerical solver developed by the authors in \cite{PizoMunoz} to approximate solutions of the model on a finite spatial interval, considering various parameter values and forms of the variable coefficients in the Boussinesq system under study. Additionally, we present preliminary numerical experiments addressing an inverse problem: the reconstruction of the initial wave elevation and fluid velocity from measurements taken at a final time. This is achieved by formulating an optimization problem in which the initial conditions are estimated as minimizers of a functional that quantifies the discrepancy between the observed final state and the numerical solution evolved from a trial initial state.
academic

Существование и численная аппроксимация одномерной системы Буссинеска с переменными коэффициентами на конечном интервале

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

  • ID статьи: 2510.10829
  • Название: Existence and numerical approximation of a one-dimensional Boussinesq system with variable coefficients on a finite interval
  • Авторы: Deissy Marcela Pizo, Juan Carlos Muñoz Grajales (Universidad del Valle, Колумбия)
  • Классификация: math.NA, cs.NA, math.AP
  • Дата публикации: 14 октября 2025 г. (препринт arXiv)
  • Ссылка на статью: https://arxiv.org/abs/2510.10829

Аннотация

В данной работе исследуется корректность постановки нелинейной дисперсионной модели с переменными коэффициентами, описывающей эволюцию поверхностных волн в одномерном канале конечной длины с неправильной топографией дна. Для дополнения теоретического анализа авторы используют численный решатель, разработанный в предыдущих работах (метод конечных элементов), для аппроксимации решений модели на ограниченном пространственном интервале с учётом различных значений параметров и функциональных форм переменных коэффициентов в системе Буссинеска. Кроме того, в работе предложены предварительные численные эксперименты для решения обратной задачи: восстановление начальной высоты волны и скорости жидкости по измерениям данных в конечный момент времени. Это достигается путём построения оптимизационной задачи, в которой начальные условия оцениваются как решение, минимизирующее целевой функционал, количественно определяющий разницу между наблюдаемым конечным состоянием и численным решением, эволюционирующим из пробных начальных условий.

Научный контекст и мотивация

  1. Решаемые проблемы:
    • Исследование математической корректности одномерной системы Буссинеска с переменными коэффициентами (существование, единственность, непрерывная зависимость)
    • Разработка эффективных численных методов решения системы с переменными коэффициентами
    • Исследование обратной задачи восстановления начальных условий по наблюдениям в конечный момент времени
  2. Значимость проблемы:
    • Система Буссинеска является важной математической моделью для описания распространения волн на мелкой воде с приложениями в морской инженерии и береговой защите
    • Реальные водоёмы часто имеют неоднородную топографию дна, поэтому модели с переменными коэффициентами более адекватны реальным ситуациям
    • Обратные задачи имеют важное практическое значение для прогнозирования волн и мониторинга морской среды
  3. Ограничения существующих методов:
    • Большинство исследований сосредоточены на системах Буссинеска с постоянными коэффициентами
    • Теоретический анализ случая переменных коэффициентов относительно недостаточен
    • Исследования обратных задач для таких дисперсионных систем практически отсутствуют
  4. Научная мотивация:
    • Заполнить пробел в теоретическом анализе систем Буссинеска с переменными коэффициентами
    • Предоставить более точные математические модели и численные инструменты для практических приложений в морской инженерии
    • Открыть новые направления в исследовании обратных задач для дисперсионных волновых систем

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

  1. Установлена теория корректности для одномерной системы Буссинеска с переменными коэффициентами: доказаны существование, единственность решения и непрерывная зависимость от начальных данных
  2. Разработан метод теоретического анализа на основе техники функций Грина: задача преобразована в задачу о неподвижной точке нелинейного оператора
  3. Получены законы сохранения энергии: доказано, что энергетический функционал системы сохраняется во времени
  4. Реализован эффективный численный решатель на основе метода конечных элементов: разработана стабильная численная схема с использованием библиотеки FEniCS
  5. Впервые исследована обратная задача для таких дисперсионных систем: предложен метод восстановления начальных условий на основе вариационной оптимизации

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

Постановка задачи

Исследуется следующая система Буссинеска с переменными коэффициентами:

N_t + [(1 + \alpha c^2(\xi)N)V]_\xi - \frac{\beta}{6}N_{\xi\xi t} = 0 \\ V_t + (c(\xi)N)_\xi + \frac{1}{2}\alpha(c^2V^2)_\xi - \frac{\beta}{6}V_{\xi\xi t} = 0 \end{cases}$$ где $(\xi,t) \in [0,L] \times [0,T]$, с начальными условиями: $$N(\xi,0) = N_0(\xi), \quad V(\xi,0) = V_0(\xi)$$ и граничными условиями Дирихле: $$N(0,t) = N(L,t) = 0, \quad V(0,t) = V(L,t) = 0$$ ### Методы теоретического анализа **1. Техника функций Грина** Система переписывается в виде: $$\left(I - \frac{\beta}{6}\partial_\xi^2\right)N_t = -\partial_\xi[(1 + \alpha c^2(\xi)N)V]$$ $$\left(I - \frac{\beta}{6}\partial_\xi^2\right)V_t = -\partial_\xi\left[c(\xi)N + \frac{1}{2}\alpha c^2V^2\right]$$ Используется функция Грина оператора $P = I - \frac{\beta}{6}\partial_\xi^2$: $$G(\xi,s) = \frac{1}{2\sqrt{\beta/6}} \frac{\cosh\left(\frac{L-|\xi-s|}{\sqrt{\beta/6}}\right) - \cosh\left(\frac{L-(\xi+s)}{\sqrt{\beta/6}}\right)}{\sinh\left(\frac{L}{\sqrt{\beta/6}}\right)}$$ **2. Теорема о неподвижной точке** Задача преобразуется в систему интегральных уравнений, определяется нелинейный оператор $A$, и через теорему Банаха о неподвижной точке доказывается существование и единственность решения. ### Численные методы **Дискретизация методом конечных элементов** Применяется метод Галёркина, слабая формулировка: $$\partial_t\left(\langle N,v_1\rangle_{L^2} + \frac{\beta}{6}\langle N_\xi, v_{1\xi}\rangle_{L^2}\right) = \langle[(1+\alpha c^2N)V], v_{1\xi}\rangle_{L^2}$$ **Дискретизация по времени** Используется $\theta$-метод ($\theta = 1/2$) для интегрирования по времени, обеспечивающий численную стабильность. ### Методология обратной задачи Определяется целевой функционал: $$J(N,V) := \frac{1}{2}\int_0^L \left(|N(\xi,T;N_0,V_0) - N_T(\xi)|^2 + |V(\xi,T;N_0,V_0) - V_T(\xi)|^2\right)d\xi$$ Для минимизации функционала используется алгоритм L-BFGS-B в сочетании с библиотекой Dolfin-Adjoint для автоматического дифференцирования. ## Экспериментальная установка ### Численные параметры - **Пространственная область**: $[-20, 40]$ - **Количество узлов сетки**: 3000 пространственных узлов - **Пространственное разрешение**: $\Delta\xi \approx 0.02$ - **Временной шаг**: $\Delta t = 8/3000 \approx 2.7 \times 10^{-3}$ - **Параметры модели**: $\alpha = \beta = 0.1$ ### Начальные условия Используется локализованный гауссов профиль: $$N(\xi,0) = V(\xi,0) = e^{-(\xi-18)^2}$$ ### Функции переменных коэффициентов **Случай 1**: Гладкие осциллирующие гауссовы возмущения $$c(\xi) = 1 + 0.3\sin\left(\frac{\pi}{5}\xi\right) + 0.6e^{-(\xi-8)^2}$$ **Случай 2**: Кусочно-постоянная функция (моделирование слоистой среды) $$c(\xi) = \begin{cases} 0.8, & \xi < 20 \\ 1.5, & 20 < \xi \leq 22 \\ 2.0, & 22 < \xi \leq 27 \\ \vdots \end{cases}$$ ## Результаты экспериментов ### Теоретические результаты **Теорема 3.1 (локальная корректность)**: Существует $T_0 > 0$ такое, что система имеет единственное решение $(N,V)$ в $L^2_{T_0} \times L^2_{T_0}$, удовлетворяющее оценке: $$\|(N,V)\|_{L^2_{T_0} \times L^2_{T_0}} \leq \|(N_0,V_0)\|_{H^1_0 \times H^1_0} T_0^{1/2} \exp\left(\frac{L^{1/2}}{(\beta/6)^{1/2}}DT_0\right)$$ **Теорема 3.2 (сохранение энергии)**: $$\|V\|_{L^2}^2 + \alpha\|c|N|^{1/2}V\|_{L^2}^2 + \||c|^{1/2}N\|_{L^2}^2 = \text{константа}$$ ### Результаты численных экспериментов 1. **Решение прямой задачи**: - Успешно смоделирована эволюция волн в среде с переменными коэффициентами - Наблюдаются изменения скорости волны, модуляция амплитуды и дисперсионные эффекты, вызванные изменениями топографии дна - Численный метод демонстрирует хорошую стабильность и точность 2. **Восстановление в обратной задаче**: - Алгоритм L-BFGS-B значительно улучшает восстановление начальных условий за 4 итерации - Восстановленные высота волны и скорость жидкости постепенно приближаются к целевым профилям - Доказана эффективность стратегии вариационной ассимиляции ### Численная верификация - Закон сохранения энергии хорошо сохраняется в численных расчётах - Согласованность результатов при различных разрешениях сетки подтверждает сходимость численного метода ## Связанные работы 1. **Теория систем Буссинеска**: - Bona & Chen (1998, 2002, 2004): основополагающие теории для случая постоянных коэффициентов - Данная работа расширяет результаты на случай переменных коэффициентов, заполняя теоретический пробел 2. **Численные методы**: - Применение метода конечных элементов к дисперсионным волновым уравнениям - В данной работе реализован эффективный решатель с использованием библиотеки FEniCS 3. **Обратные задачи**: - Традиционные методы для обратных задач волновых уравнений - Впервые применена вариационная оптимизация к обратной задаче для системы Буссинеска ## Заключение и обсуждение ### Основные выводы 1. Успешно установлена математическая теория одномерной системы Буссинеска с переменными коэффициентами 2. Разработаны эффективные численные методы решения 3. Открыто новое направление в исследовании обратных задач для дисперсионных систем 4. Предоставлены теоретические и вычислительные инструменты для практических приложений в морской инженерии ### Ограничения 1. **Теоретический аспект**: - Доказана только локальная корректность; глобальное существование требует дальнейшего исследования - Отсутствует теоретический анализ единственности и устойчивости обратной задачи 2. **Численный аспект**: - Эксперименты по обратной задаче рассматривают только идеальный случай без шума - Отсутствует анализ чувствительности к шуму в наблюдениях 3. **Практический аспект**: - Ограничение одномерным случаем; практические приложения требуют расширения на двумерный случай - Выбор функций переменных коэффициентов относительно прост ### Направления будущих исследований 1. **Теоретическое развитие**: - Теория глобальной корректности - Анализ единственности и устойчивости обратной задачи - Теория регуляризации 2. **Численные улучшения**: - Робастные методы восстановления при наличии шума в данных - Расширение на многомерные случаи - Интеграция методов машинного обучения 3. **Практические приложения**: - Верификация на реальных морских данных - Интеграция с системами морских наблюдений ## Глубокая оценка ### Преимущества 1. **Значительный теоретический вклад**: первое систематическое исследование корректности систем Буссинеска с переменными коэффициентами 2. **Методологическая инновация**: умелое применение техники функций Грина для работы с переменными коэффициентами 3. **Полная численная реализация**: предоставлена полная схема численного решения 4. **Новизна проблемы**: открыто новое направление в исследовании обратных задач для дисперсионных систем 5. **Ясное изложение**: строгие математические выводы и чёткое представление материала ### Недостатки 1. **Ограниченная теоретическая глубина**: установлена только локальная теория; отсутствует глобальный анализ 2. **Недостаточная экспериментальная верификация**: эксперименты по обратной задаче относительно простые; отсутствует верификация на сложных сценариях 3. **Неподтверждённая практическая применимость**: требуется верификация на реальных данных 4. **Отсутствие анализа вычислительной сложности**: не проанализирована сложность алгоритмов ### Влияние 1. **Академическая ценность**: закладывает основы для развития теории дисперсионных систем с переменными коэффициентами 2. **Перспективы применения**: потенциальное применение в морской инженерии, прогнозировании волн и других областях 3. **Методологический вклад**: сочетание техники функций Грина и вариационной оптимизации предлагает новые подходы для решения аналогичных задач ### Области применения 1. Моделирование распространения волн на мелкой воде 2. Анализ волн в прибрежной инженерии 3. Обратные задачи в мониторинге морской среды 4. Ассимиляция данных в численной океанографии ## Библиография В работе цитируется 20 важных источников, охватывающих теорию систем Буссинеска, численные методы, функциональный анализ и другие смежные области, обеспечивая прочную теоретическую базу для исследования.