🔢 mathematics

A volume penalization method for solving conjugate scalar transport with interfacial jump conditions

本文提出了一种新颖的体积惩罚方法，该方法能够有效处理复杂几何形状下共轭标量输运问题中的界面跳跃条件，与体拟合网格模拟相比，其精度高，偏差小于3%。

原作者： Ming Liu, Yosuke Hasegawa

发布于 2026-08-11

📖 1 分钟阅读🧠 深度阅读

原作者： Ming Liu, Yosuke Hasegawa

原始论文采用 CC BY 4.0 许可（http://creativecommons.org/licenses/by/4.0/）。 ✨ 这是对下方论文的AI生成解释。它不是由作者撰写或认可的。如需技术准确性，请参阅原始论文。阅读完整免责声明

想象一个热量和化学物质不再只是平滑地流过管道或房间，而是必须跳过隐形障碍的世界。在现实世界中，当热气体接触到冰冷的固体壁面，或者当化学反应恰好发生在液相与气相的边界时，情况会变得非常复杂。温度可能会突然飙升，或者化学物质的浓度可能会发生瞬时变化。科学家们称之为“具有界面跳跃条件的共轭标量输运”。这是一种高级的说法，意在表达：“当我们追踪热量或污染物时，如何处理它们在墙壁两侧表现出的差异，以及它们有时甚至如何跨越墙壁的问题？”

为了解决这个问题，研究人员通常使用计算机模拟。传统的方法是构建一个能够完美契合物体每一个曲线的数字地图（网格），就像用黏土包裹雕像一样。这种方法很精确，但如果物体的形状很奇特，比如一张揉皱的纸或一块粗糙的岩石，制作起来就会是一场噩梦。一种更新颖、更酷的方法被称为“体积惩罚法”（Volume Penalization Method, VPM）。把 VPM 想象成一面“幽灵墙”。与其用黏土去包裹物体，你只需在一个简单的网格上将物体“画”出来。计算机将固体物体视为一个超厚、具有粘性的海绵，阻止物质穿过它。这种方法快速且简便，即使面对极其复杂的形状也游刃有余。但直到目前为止，这个“幽灵墙”技巧一直存在一个缺陷：它无法处理边界处热量或化学物质的突然跳跃，否则会在固体内部产生虚假且荒谬的结果。

本文介绍了一个针对该缺陷的巧妙修复方案。作者刘明（Ming Liu）和长谷川洋介（Yosuke Hasegawa）开发了一种使用体积惩罚法的新方法，该方法允许在不破坏模拟的前提下处理这些“跳跃条件”。他们并没有仅仅强制墙壁阻挡物质，而是添加了一个特殊的“源项”——可以将其想象成位于边界处的一个微小的、隐形的泵——它能推动适量的热量或化学物质跨越界线，即使其量发生了突变。他们通过几种场景测试了这种新方法，包括一个简单的单维扩散问题和更复杂的流固相互作用。结果显示，他们的新型“幽灵墙”具有极高的准确性，与传统的、制作困难的“黏土包裹法”相比，平均误差小于 3.0%。至关重要的是，与旧版本的方法不同，他们的新方法不会在固体内部产生虚假数据，这使其成为模拟从微通道到工业反应器等各种复杂热力学和化学过程的可靠工具。

技术摘要：一种用于求解具有界面跳跃条件的共轭标量输运的体积惩罚法

问题陈述
涉及界面跳跃条件的共轭标量输运（例如热量或质量传递）在热学和化学过程中广泛存在，如微通道传热、界面处的化学反应以及相变系统。在这些场景中，标量本身（例如温度或浓度跳跃）以及跨流固界面的通量往往存在不连续性。虽然体网格法（BFM）模拟可以通过显式耦合独立的相域来处理这些条件，但其在处理复杂几何形状时需要繁琐的网格生成。浸没边界法（IBM），特别是体积惩罚法（VPM），提供了一种将任意几何形状嵌入笛卡尔网格的方法，无需复杂的网格划分。然而，现有的 VPM 方案在处理 Neumann 边界条件和一般的界面跳跃条件方面表现不佳。传统的处理方法会在固体区域内产生非物理的标量场，并且在没有额外的、复杂的耦合程序的情况下，无法精确强制执行通量跳跃或标量不连续。

方法论
作者提出了一种在 VPM 框架内的创新处理方法，用于求解具有一般界面跳跃条件的共轭标量输运。该方法分为三个主要阶段：

水平集表示： 任意流固界面由水平集函数 ( $\phi_0$ ) 表示，并转换为相指示函数 ( $\phi$ )，该函数在厚度为 $2\delta$ 的薄界面区域内从 0（流体）平滑过渡到 1（固体）。
Neumann 边界条件处理： 对于具有 Neumann 条件的单相平流扩散问题，作者推导出了一个散度形式的源项 $\nabla \cdot (q_w \phi \mathbf{n})$ ，将其添加到控制方程中。与以往通过在整个固体体积内施加惩罚项来强制导数的方法不同，这个新项仅在界面区域内非零。这确保了所规定的标量通量 $q_w$ 被精确施加，同时保持固体内部的分子扩散率为零，从而防止了向固体相的非物理扩散。
推广至具有跳跃条件的共轭输运： 为了处理标量跳跃（ $c_f = \alpha c_s$ ）和通量跳跃，作者引入了一个“等效标量” $h = \alpha c$ 。这种变换消除了标量不连续性，从而允许对两种相使用统一的控制方程。通量跳跃条件被重新表述为等效通量跳跃 $q^*_w$ ，它结合了预设的通量跳跃和标量跳跃的影响。等效标量 $h$ 的统一控制方程为：
$\frac{\partial h}{\partial t} + \mathbf{u} \cdot \nabla h = \nabla \cdot (D \nabla h) + \nabla \cdot (q^*_w \phi \mathbf{n})$
其中 $D$ 是相相关的扩散率。

该方案在开源有限体积计算代码 OpenFOAM（版本 8）中实现，空间项采用二阶线性离散，瞬态项采用 Euler 隐式格式。

核心贡献

创新的源项公式： 本文引入了一种散度形式的源项用于 VPM，能够精确强制执行界面通量跳跃，而不会在固体区域内产生非物理解。
统一的控制方程： 通过定义等效标量和通量，该方法统一了流体相和固体相的控制方程，消除了对显式界面耦合或独立域求解的需求。
推广至一般跳跃条件： 该方法成功地将 VPM 扩展到能够同时处理标量值和通量的跳跃，这是此前基于有限体积的浸没边界方案所缺乏的能力。
误差分析： 作者提供了解析误差分析，证明即使在存在界面曲率的情况下，该方案相对于网格尺寸仍能达到一阶精度，这与传统的 VPM 一致。

结果
所提方法通过三个不同的测试用例进行了验证：

一维扩散： 解决了一个带有 Neumann 边界条件的瞬态扩散问题。结果显示与解析解吻合良好，平均相对误差约为 1.25%。至关重要的是，与现有的 VPM 方案不同，本方法在固体内部保持了零标量场，避免了非物理扩散。误差在不同的扩散率值下保持稳定，而现有方案在扩散率较高时表现出误差增加的现象。
流固耦合扩散： 在正方形区域内对一个带有标量和通量跳跃的圆柱体进行了模拟。VPM 结果与使用 chtMultiRegionFoam 的体网格法（BFM）模拟进行了对比。在四种具有不同扩散率比和跳跃强度的测试用例中，VPM 与 BFM 结果之间的相对偏差在 3.0% 以内（平均为 1.46%）。
流固耦合平流-扩散： 模拟了带有界面跳跃的圆柱体周围稳态层流（$Re=40$）过程。VPM 结果再次显示出与 BFM 基准的高度一致性，相对偏差为 1.59%。该方法准确地重现了界面附近的标量分布和通量跳跃。

意义与主张
本文声称，所提出的 VPM 为涉及复杂界面几何形状和跳跃条件的共轭标量输运问题提供了一种鲁棒且高效的替代方案。其主要意义在于能够使用单一的统一控制方程解决耦合流固问题，同时精确强制执行局部跳跃条件。作者指出，这种方法避免了传统 VPM 实现中常见的针对 Neumann 条件的非物理伪影。

作者谦虚地将这项工作定位为基础性的一步，指出该方法处理复杂几何形状的能力使其适用于未来的粗糙壁面和复杂换热表面的应用。他们还强调了将此 VPM 与伴随方法相结合，用于辐射、湍流及共轭换热问题中的形状和拓扑优化的潜力，尽管这些具体应用被确定为未来研究的方向，而非当前研究的直接结果。

类似论文