ALE法在流体-结构相互作用中的应用
理论与物理
ALE法的基本概念
老师，ALE法是“Arbitrary Lagrangian-Eulerian”的缩写对吧？它结合了拉格朗日描述和欧拉描述的哪些优点呢？
欧拉描述中网格固定，流体流过网格。它擅长处理大变形但不擅长界面追踪。拉格朗日描述中网格随物质一起移动，因此界面清晰，但在大变形下网格会畸变。ALE法是两者的折中，可以任意（Arbitrary）设置网格速度。在界面处让网格随物质一起移动（拉格朗日式），在内部则适度移动网格以保持网格质量（欧拉式）。
ALE描述下的Navier-Stokes方程动量守恒定律如下。
这里 $\mathbf{c} = \mathbf{u} - \hat{\mathbf{u}}$ 是流体速度 $\mathbf{u}$ 与网格速度 $\hat{\mathbf{u}}$ 的相对速度（对流速度）。
也就是说，如果网格速度 $\hat{\mathbf{u}}$ 为零就是欧拉法，$\hat{\mathbf{u}} = \mathbf{u}$ 就是拉格朗日法，对吧？
没错。在FSI问题中，在流体-结构界面上令 $\hat{\mathbf{u}} = \dot{\mathbf{d}}_s$（结构速度），内部的网格速度则通过平滑算法确定。
GCL（几何守恒定律）
ALE法有什么需要特别注意的地方吗？
最重要的是满足GCL（几何守恒定律）。当网格移动时，如果守恒定律的离散化不一致，就会对均匀流产生虚假的源项。
这个公式要求单元体积的时间变化与网格速度的通量一致。Fluent或STAR-CCM+会自动满足GCL，但在OpenFOAM的自定义求解器中需要特别注意。
如果GCL不满足会发生什么？
即使输入均匀流，也会产生非物理的压力或能量波动。特别是时间精度为一阶时，GCL误差会很明显。使用二阶精度时间积分，并在计算网格速度时正确使用上一步的网格位置非常重要。
网格平滑方法
用于确定网格内部速度的平滑算法有哪些？
我们来比较一下代表性的方法。
| 方法 | 原理 | 大变形耐受性 | 成本 |
|---|---|---|---|
| 拉普拉斯平滑 | 移动到相邻节点的平均位置 | 低 | 非常低 |
| 弹簧类比法 | 通过弹簧网络传播变形 | 中 | 低 |
| 基于扩散的方法 | 求解扩散方程传播位移 | 高 | 中 |
| RBF（径向基函数） | 从界面位移进行RBF插值 | 非常高 | 高 |
| 基于弹性的方法 | 通过弹性体方程传播位移 | 高 | 高 |
Fluent的基于扩散的平滑方法中，如果将扩散率设为与壁面距离的倒数（Boundary Distance选项），则可以抑制壁面附近的网格变形，保持棱柱层质量。
拉格朗日与欧拉——经典力学的两大视角在FSI中的融合
拉格朗日描述是“跟随粒子一起运动”的视角，欧拉描述是“在空间固定点观察流动”的视角。结构分析擅长拉格朗日法，流体分析则以欧拉法为主流。ALE法诞生于“适度结合两者”的构想。19世纪数学家提出的描述方法，成为了解决21世纪复杂FSI问题的关键，可见基础理论的生命力出乎意料地长久。
各项的物理意义
- 时间项 $\partial(\rho\phi)/\partial t$：想象一下拧开水龙头的瞬间。最初水流会不稳定地喷溅，过一会儿才会变成稳定的水流，对吧？描述这个“变化过程中”的就是时间项。心脏搏动导致血流脉动，发动机阀门每次开闭引起流动变化，这些都是非定常现象。那么定常分析是什么？就是只看“经过足够时间流动稳定之后”——也就是将此项设为零。计算成本会大幅下降，因此先用定常分析求解是CFD的基本策略。
- 对流项 $\nabla \cdot (\rho \mathbf{u} \phi)$：把落叶扔进河里会怎样？会被水流带着往下游漂，对吧？这就是“对流”——流体运动搬运物质的效果。暖风的暖气能到达房间的另一端，也是因为空气这个“搬运工”通过对流输送了热量。这里有趣的是——这项包含了“速度×速度”，因此是非线性的。也就是说，流速加快时此项会急剧增强，变得难以控制。这就是湍流的根本原因。常见的误解：“对流和传导差不多”→ 完全不一样！对流是流动搬运，传导是分子传递。效率有天壤之别。
- 扩散项 $\nabla \cdot (\Gamma \nabla \phi)$：有过在咖啡里倒入牛奶后放置不管的经历吗？即使不搅拌，过一会儿也会自然混合，对吧？那就是分子扩散。那么下一个问题——蜂蜜和水，哪个更容易流动？当然是水，对吧？因为蜂蜜的粘度（$\mu$）高，所以不易流动。粘度越大，扩散项越强，流体的运动就越“粘稠”。雷诺数小的流动（缓慢、粘稠）中扩散占主导。相反，Re数大的流动中对流占压倒性优势，扩散则成为配角。
- 压力项 $-\nabla p$：推注射器的活塞，液体就会从针头有力地射出，对吧？为什么呢？因为活塞侧压力高，针头侧压力低——这个压力差产生了推动流体的力。水坝放水也是同样的原理。天气图上等压线密集的地方会怎样？没错，会刮强风。“有压力差的地方就会产生流动”——这就是纳维-斯托克斯方程压力项的物理意义。这里的误解点：CFD中的“压力”多为表压而非绝对压力。切换到可压缩分析时结果突然出错，原因可能就是混淆了绝对压力/表压。
- 源项 $S_\phi$：被加热的空气会上升——为什么呢？因为比周围空气轻（密度低），被浮力推上去了。这个浮力作为源项添加到方程中。此外，燃气灶火焰产生化学反应热、工厂电磁泵对金属熔液施加洛伦兹力……这些都是“从外部向流体注入能量或力”的作用，都用源项表示。忘记源项会怎样？自然对流分析中忘记加入浮力，流体就完全不动——就像冬天房间里开了暖气但暖空气不上升一样，会得到物理上不可能的结果。
假设条件与适用范围
- 连续介质假设：克努森数 Kn < 0.01（分子平均自由程 ≪ 特征长度）时成立
- 牛顿流体假设：剪切应力与应变速率呈线性关系（非牛顿流体需要粘度模型）
- 不可压缩假设（Ma < 0.3时）：将密度视为常数。马赫数0.3以上需考虑可压缩性效应
- Boussinesq近似（自然对流）：仅在浮力项中考虑密度变化，其他项使用恒定密度
- 不适用的情形：稀薄气体（Kn > 0.1）、超音速/高超音速流动（需要激波捕捉）、自由表面流动（需要VOF/Level Set等）
量纲分析与单位制
| 变量 | SI单位 | 注意事项·换算备忘 |
|---|---|---|
| 速度 $u$ | m/s | 从入口条件的体积流量换算时，注意截面积单位 |
| 压力 $p$ | Pa | 区分表压与绝对压力。可压缩分析中使用绝对压力 |
| 密度 $\rho$ | kg/m³ | 空气：约1.225 kg/m³@20°C，水：约998 kg/m³@20°C |
| 粘性系数 $\mu$ | Pa·s | 注意与运动粘性系数 $\nu = \mu/\rho$ [m²/s] 混淆 |
| 雷诺数 $Re$ | 无量纲 | $Re = \rho u L / \mu$。层流/湍流转换的判断指标 |
| CFL数 | 无量纲 | $CFL = u \Delta t / \Delta x$。直接关系到时间步长的稳定性 |
数值解法与实现
弱耦合与强耦合
流体与结构耦合算法的选择，如何影响分析的稳定性？
耦合的强弱会极大地影响稳定性和精度。
| 耦合方式 | 各步骤流程 | 稳定性 | 精度 | 成本 |
|---|---|---|---|---|
| 弱耦合（显式） | F→S→F（仅一次） | 低 | 低〜中 | 低 |
| 强耦合（隐式） | F⇆S（迭代至收敛） | 高 | 高 | 高 |
| 整体式 | F+S同时求解 | 最高 | 最高 | 非常高 |
弱耦合存在一个时间步的滞后（time lag）。当密度比 $\rho_s/\rho_f$ 较大时（如金属结构+空气），弱耦合也能稳定；但当 $\rho_s/\rho_f \approx 1$（如血管+血液、橡胶+水）时，会发生附加质量不稳定性，无法使用弱耦合。
飞机机翼颤振分析应该用哪种？
对于空气中的结构（$\rho_s/\rho_f \gg 1$），弱耦合在多数情况下也能稳定。但如果重视精度，则应使用强耦合。对于水中的结构（$\rho_s/\rho_f \sim 1-10$），则必须使用强耦合。
Aitken松弛法
有没有方法可以加速强耦合迭代的收敛？
Aitken $\Delta^2$ 加速法是最简单且最有效的。它能动态优化耦合迭代的松弛系数。
这里 $\mathbf{r}^k$ 是第 $k$ 次迭代的界面残差。使用固定的松弛系数（例如 $\omega = 0.5$）可能需要10-20次迭代，而使用Aitken加速通常3-5次就能收敛。
哪些软件可以使用Aitken加速？
preCICE已内置。Ansys System Coupling最近也添加了此功能。在OpenFOAM中，solidDisplacementFoam和pimpleFoam的耦合需要自定义实现，但如果使用preCICE适配器则很简单。
ALE法的局限：重划网格
如果ALE法中网格无论如何都会失效，该怎么办？
当结构变形量达到与原始网格尺寸相当时，仅靠平滑无法应对。此时需要自动重划网格。
Remeshing选项。当单元质量低于阈值时，会自动局部重新生成四面体网格dynamicRefineFvMesh进行局部细化/粗化，使用tetDecomposition进行重划网格但是，重划网格会产生解的插值误差。特别是边界层棱柱网格的重划，其质量管控很困难。当变形量非常大时，应考虑切换到重叠网格法或浸入边界法。
ALE法中“网格速度”这个奇妙的量
在ALE法中，除了“流体速度”、“结构速度”，还出现了“网格速度”这第三个速度。在实际工作中，工程师第一次看到ALE设置界面时困惑地问“这是什么？”并不少见。网格速度终究是为了计算方便而产生的虚拟量，没有物理意义。但是，通过巧妙地控制这个“虚拟的速度”，即使是不合理的变形也能让网格不畸变而继续计算，可见数值分析真是深奥。
迎风格式（Upwind）
一阶迎风：数值扩散大但稳定。二阶迎风：精度提高但有振荡风险。高雷诺数流动中必须使用。
中心差分（Central Differencing）
二阶精度，但当Pe数 > 2时会产生数值振荡。适用于低雷诺数的扩散主导流动。
TVD格式（MUSCL、QUICK等）
通过限制器函数抑制数值振荡，同时保持高精度。对捕捉激波或陡峭梯度有效。
有限体积法 vs 有限元法
FVM: 自然满足守恒定律。CFD的主流。FEM: 对复杂形状、多物理场有利。SPH等无网格法也在发展中。
CFL条件（库朗数）
显式方法: CFL ≤ 1为
相关主题
なった
詳しく
報告