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

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

基本信息

  • 论文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, Colombia)
  • 分类: math.NA, cs.NA, math.AP
  • 发表时间: October 14, 2025 (arXiv预印本)
  • 论文链接: https://arxiv.org/abs/2510.10829

摘要

本文研究了一个具有变系数的非线性色散模型的适定性问题，该模型描述了表面波在具有不规则底部地形的一维有限长度浅水通道中的传播演化。为补充理论分析，作者利用其在先前工作中开发的数值求解器（有限元方法）来近似求解有界空间区间上的模型解，考虑了Boussinesq系统中变系数的各种参数值和函数形式。此外，本文还提出了初步的数值实验来解决一个反问题：从最终时刻的测量数据重构初始波高和流体速度。这通过构建一个优化问题来实现，其中初始条件被估计为最小化目标泛函的解，该泛函量化了观测的最终状态与从试验初始状态演化的数值解之间的差异。

研究背景与动机

  1. 要解决的问题：
    • 研究具有变系数的一维Boussinesq系统的数学适定性（存在性、唯一性、连续依赖性）
    • 开发有效的数值方法求解该变系数系统
    • 探索基于最终时刻观测数据重构初始条件的反问题
  2. 问题重要性：
    • Boussinesq系统是描述浅水波传播的重要数学模型，在海洋工程、海岸防护等领域有重要应用
    • 实际水域往往具有不均匀的底部地形，因此变系数模型更贴近实际情况
    • 反问题在波浪预测、海洋监测等方面具有重要实用价值
  3. 现有方法局限性：
    • 大多数研究集中在常系数Boussinesq系统
    • 对于变系数情况的理论分析相对缺乏
    • 针对此类色散系统的反问题研究尚属空白
  4. 研究动机：
    • 填补变系数Boussinesq系统理论分析的空白
    • 为实际海洋工程应用提供更精确的数学模型和数值工具
    • 开拓色散波动系统反问题研究的新方向

核心贡献

  1. 建立了变系数一维Boussinesq系统的适定性理论：证明了解的存在性、唯一性和对初始数据的连续依赖性
  2. 发展了基于Green函数技术的理论分析方法：将问题转化为非线性算子的不动点问题
  3. 提供了能量守恒定律：证明了系统的能量泛函在时间上保持守恒
  4. 实现了高效的有限元数值求解器：结合FEniCS库开发了稳定的数值方案
  5. 首次研究了此类色散系统的反问题：提出了基于变分优化的初始条件重构方法

方法详解

任务定义

研究如下变系数Boussinesq系统： {Nt+[(1+αc2(ξ)N)V]ξ−β6Nξξt=0Vt+(c(ξ)N)ξ+12α(c2V2)ξ−β6Vξξt=0\begin{cases} 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}

其中 (ξ,t)∈[0,L]×[0,T](\xi,t) \in [0,L] \times [0,T]，配以初始条件： N(ξ,0)=N0(ξ),V(ξ,0)=V0(ξ)N(\xi,0) = N_0(\xi), \quad V(\xi,0) = V_0(\xi)

和Dirichlet边界条件： N(0,t)=N(L,t)=0,V(0,t)=V(L,t)=0N(0,t) = N(L,t) = 0, \quad V(0,t) = V(L,t) = 0

理论分析方法

1. Green函数技术 将系统重写为： (I−β6∂ξ2)Nt=−∂ξ[(1+αc2(ξ)N)V]\left(I - \frac{\beta}{6}\partial_\xi^2\right)N_t = -\partial_\xi[(1 + \alpha c^2(\xi)N)V](I−β6∂ξ2)Vt=−∂ξ[c(ξ)N+12αc2V2]\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−β6∂ξ2P = I - \frac{\beta}{6}\partial_\xi^2 的Green函数： G(ξ,s)=12β/6cosh⁡(L−∣ξ−s∣β/6)−cosh⁡(L−(ξ+s)β/6)sinh⁡(Lβ/6)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. 不动点定理 将问题转化为积分方程组，定义非线性算子 AA，通过Banach不动点定理证明解的存在唯一性。

数值方法

有限元离散 采用Galerkin方法，弱形式为： ∂t(⟨N,v1⟩L2+β6⟨Nξ,v1ξ⟩L2)=⟨[(1+αc2N)V],v1ξ⟩L2\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-方法（θ=1/2\theta = 1/2）进行时间积分，确保数值稳定性。

反问题方法

定义目标泛函： J(N,V):=12∫0L(∣N(ξ,T;N0,V0)−NT(ξ)∣2+∣V(ξ,T;N0,V0)−VT(ξ)∣2)dξ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][-20, 40]
  • 网格点数：3000个空间格点
  • 空间分辨率：Δξ≈0.02\Delta\xi \approx 0.02
  • 时间步长：Δt=8/3000≈2.7×10−3\Delta t = 8/3000 \approx 2.7 \times 10^{-3}
  • 模型参数：α=β=0.1\alpha = \beta = 0.1

初始条件

采用局部化高斯轮廓： N(ξ,0)=V(ξ,0)=e−(ξ−18)2N(\xi,0) = V(\xi,0) = e^{-(\xi-18)^2}

变系数函数

情况1：光滑振荡-高斯扰动 c(ξ)=1+0.3sin⁡(π5ξ)+0.6e−(ξ−8)2c(\xi) = 1 + 0.3\sin\left(\frac{\pi}{5}\xi\right) + 0.6e^{-(\xi-8)^2}

情况2：分段常数函数（模拟分层介质） c(ξ)={0.8,ξ<201.5,20<ξ≤222.0,22<ξ≤27⋮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（局部适定性）： 存在 T0>0T_0 > 0，使得系统在 LT02×LT02L^2_{T_0} \times L^2_{T_0} 中有唯一解 (N,V)(N,V)，且满足估计： ∥(N,V)∥LT02×LT02≤∥(N0,V0)∥H01×H01T01/2exp⁡(L1/2(β/6)1/2DT0)\|(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∥L22+α∥c∣N∣1/2V∥L22+∥∣c∣1/2N∥L22=常数\|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. Boussinesq系统理论：
    • Bona & Chen (1998, 2002, 2004)：建立了常系数情况的基础理论
    • 本文扩展到变系数情况，填补了理论空白
  2. 数值方法：
    • 有限元方法在色散波动方程中的应用
    • 本文结合FEniCS库实现了高效求解器
  3. 反问题：
    • 波动方程反问题的传统方法
    • 本文首次将变分优化应用于Boussinesq系统反问题

结论与讨论

主要结论

  1. 成功建立了变系数一维Boussinesq系统的数学理论基础
  2. 发展了有效的数值求解方法
  3. 开创了此类色散系统反问题研究的新方向
  4. 为实际海洋工程应用提供了理论和计算工具

局限性

  1. 理论方面：
    • 仅证明了局部适定性，全局存在性有待研究
    • 反问题的唯一性和稳定性理论分析缺失
  2. 数值方面：
    • 反问题实验仅考虑了理想的无噪声情况
    • 缺乏对观测噪声敏感性的分析
  3. 应用方面：
    • 限于一维情况，实际应用需要扩展到二维
    • 变系数函数的选择相对简单

未来方向

  1. 理论发展：
    • 全局适定性理论
    • 反问题的唯一性和稳定性分析
    • 正则化理论
  2. 数值改进：
    • 噪声数据的鲁棒重构方法
    • 高维情况的扩展
    • 机器学习技术的融合
  3. 实际应用：
    • 真实海洋数据的验证
    • 与海洋观测系统的结合

深度评价

优点

  1. 理论贡献显著：首次系统研究变系数Boussinesq系统的适定性
  2. 方法创新：巧妙运用Green函数技术处理变系数问题
  3. 数值实现完整：提供了完整的数值求解框架
  4. 问题新颖：开创了色散系统反问题研究新领域
  5. 写作清晰：数学推导严谨，表述清晰

不足

  1. 理论深度有限：仅建立了局部理论，缺乏全局分析
  2. 实验验证不足：反问题实验相对简单，缺乏复杂场景验证
  3. 实用性待验证：需要更多实际数据的验证
  4. 计算复杂度分析缺失：未分析算法的计算复杂度

影响力

  1. 学术价值：为变系数色散系统理论发展奠定基础
  2. 应用前景：在海洋工程、波浪预测等领域有潜在应用价值
  3. 方法论贡献：Green函数方法和变分优化的结合为相关问题提供了新思路

适用场景

  1. 浅水波传播建模
  2. 海岸工程中的波浪分析
  3. 海洋环境监测中的反演问题
  4. 数值海洋学中的数据同化

参考文献

论文引用了20篇重要文献，涵盖了Boussinesq系统理论、数值方法、函数分析等相关领域的经典工作，为研究提供了坚实的理论基础。