[论文解读] Numerical Analysis of Time-Dependent Galbrun Equation in an Infinite Duct

Kamel Berriri, Anne-Sophie Bonnet-Ben Dhia+1 more|arXiv (Cornell University)|Mar 22, 2006

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

一句话总结

本文提出了一种针对具有均匀流的刚性管道中时间依赖的Galbrun方程的稳定有限元方法，采用受麦克斯韦方程启发的正则化技术，以实现稳定的Lagrange有限元近似。关键贡献在于引入人工涡量正则化并精心构造吸收边界条件的变分格式，从而防止数值不稳定，并实现对流动流体中声波传播的精确模拟。

ABSTRACT

In this paper we are interested in the mathematical and numerical analysis of the time-dependent Galbrun equa- tion in a rigid duct. This equation models the acoustic propagation in presence of flow [1]. We propose a regu- larized variational formulation of the problem, in the sub- sonic case, suitable for an approximation by Lagrange finite elements, and corresponding absorbing boundary conditions.

研究动机与目标

• 解决标准节点有限元在求解时间依赖Galbrun方程以模拟流动流体中声传播时的数值不稳定性问题。
• 克服时间瞬态情况下Galbrun方程变分公式缺乏自然函数框架的问题。
• 为亚音速、均匀流状态下的标准Lagrange有限元开发一种合适的正则化变分公式。
• 构建稳定且变分兼容的吸收边界条件，以最小化无限管道模拟中的虚假反射。
• 实现对包含流动对流效应的瞬态声波传播的精确数值模拟。

提出的方法

• 在Galbrun方程中引入正则化项 $ s\,\text{curl}(\text{curl}\,\boldsymbol{\xi} - \psi) $，其中 $ \psi = \text{curl}\,\boldsymbol{\xi} $，以稳定变分格式。
• 在均匀流条件下（$ M' = 0 $）解耦涡量方程，使 $ \psi $ 可被视为已知场。
• 在希尔伯特空间框架下使用 $ \mathbf{H}_0(\Omega) \times L^2(\Omega)^2 $ 构造正则化问题，通过Hille-Yosida理论确保适定性。
• 推导涉及双线性形式 $ a(\cdot,\cdot) $、$ b(\cdot,\cdot) $、$ c^{\Gamma^{\pm}}(\cdot,\cdot) $ 和 $ d^{\Gamma^{\pm}}(\cdot,\cdot) $ 的变分格式，包含材料导数和边界项。
• 利用 $ c^{\Gamma^{\pm}} $ 和 $ d^{\Gamma^{\pm}} $ 项构造与变分结构兼容的吸收边界条件，以减少波的反射。
• 采用连续Lagrange有限元对系统进行空间离散化，并使用中心二阶时间有限差分格式，得到半离散的常微分方程系统。

实验结果

研究问题

• RQ1 能否在缺乏自然变分框架的情况下，为具有均匀流的管道中时间依赖的Galbrun方程开发一种稳定且收敛的有限元方法？
• RQ2 在瞬态声学问题的Galbrun公式中，如何克服标准节点有限元的数值不稳定性？
• RQ3 正则化在使用标准有限元实现Galbrun方程一致变分公式中起什么作用？
• RQ4 如何设计吸收边界条件以最小化虚假反射，同时保持问题的变分结构？
• RQ5 正则化参数 $ s $ 对数值解的稳定性和准确性有何影响？

主要发现

• 当不使用正则化（$ s = 0 $）时，标准有限元方法导致数值不稳定，数值模拟中波振幅无界增长。
• 当使用正则化（$ s = 1 $）时，方法产生稳定且物理解释一致的结果，能清晰区分无旋（向外传播）与有旋（半径不变）波前。
• 位移的无旋部分形成随时间扩展的圆形波前，而有旋部分保持局部化并随流场对流，证实了预期的物理行为。
• 吸收边界条件显著减少了虚假反射，第二组数值实验中使用高斯时间信号时，波反射部分被有效抑制。
• 带有正则化项的变分格式确保了能量估计成立，能量泛函满足 $ \frac{d}{dt} \left( \|\boldsymbol{\xi}_t\|^2 + \|\boldsymbol{\nabla}\boldsymbol{\xi}\|^2 - M^2 \|\partial_x \boldsymbol{\xi}\|^2 \right) \leq 0 $，表明能量衰减与稳定性。
• 该方法成功以弱形式处理材料导数 $ \frac{D}{Dt} = \partial_t + M(y)\partial_x $，实现了对流动介质中对流波传播的精确模拟。