🔬 physics

Monolithic kinetic algorithm for heterogeneous porous media systems using a continuous one-domain approach

本文提出了一种单体格子玻尔兹曼模型，该模型利用连续单域方法，通过恢复具有正确孔隙率定标、伽利略不变性和增强数值稳定性的广义体积平均方程，来精确模拟非均匀多孔介质中的弱可压缩等温流，旨在应用于器件级设计。

原作者： Nikita Gusev, Ilya Karlin

发布于 2026-06-30

📖 1 分钟阅读☕ 轻松阅读

原作者： Nikita Gusev, Ilya Karlin

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

想象一下，你正试图预测水如何在复杂的景观中流动。有时，水在宽阔、开阔的河流中自由流动（自由流体）。有时，它必须挤过茂密的森林或厚厚的海绵（多孔介质）。在现实世界中，这两个世界经常交汇：一条河流可能会流入沼泽，或者一个过滤器可能会安装在管道内。

长期以来，计算机科学家必须使用两套不同的规则手册来模拟这一过程。一套规则手册用于开阔的河流，而另一套完全不同且复杂的规则则用于海绵。当水从河流进入海绵时，计算机必须停止运行，切换规则手册，并尝试将两个模拟过程在边界处粘合在一起。这种做法既混乱，又容易出错，而且计算成本高昂。

核心理念：一套通用的规则手册

这项研究中的研究人员 Nikita Gusev 和 Ilya Karlin 开发了一种全新的“单体式”（一体化）计算机算法。你可以把它想象成一个通用的翻译官，能够同时理解开阔河流和海绵的语言。

与其将河流和海绵视为两个独立的地方，他们的模型将整个计算机屏幕视为一个连续的景观。在这个景观中，“海绵度”（孔隙率）可以从 0%（固体岩石）到 100%（空隙）平滑变化。计算机不需要知道河流在哪里结束，也不需要知道海绵从哪里开始；它只需根据局部区域的“海绵程度”来计算流动。

它是如何工作的：“粒子人群”类比

为了实现这一点，他们使用了一种称为格子玻尔兹曼方法 (Lattice Boltzmann Method, LBM) 的方法。想象一下，流体并不是连续的河流，而是由大量微小的、看不见的粒子组成的，这些粒子在网格方格之间跳跃。

在开阔区域，粒子自由跳跃。
在多孔区域，粒子会撞到看不见的“障碍物”（海绵的固体部分），这会减慢它们的速度并改变它们的运动方向。

该论文的创新之处在于，他们设计了一种新的方式来编程这些粒子，使它们在移动到自由流动区域或多孔区域时，能够自动感知并调整行为，而无需计算机在两者之间画出一条明确的界线。

“交通警察”问题

以往模拟中的一个主要难题是，当“海绵度”发生变化时，数学计算无法正确对齐。

旧的方法： 想象一名试图管理车辆的交通警察。在开阔的高速公路上，警察按正常方式计数。在施工区域（多孔介质）中，警察试图计数，但数学计算变得非常混乱。压力和速度的计算会根据施工障碍物的密度而产生偏差。这意味着如果施工区域不是完全均匀的，模拟就会崩溃。
新的方法： 作者们修正了这位“交通警察”的数学逻辑。他们创建了一套新的规则，确保即使施工障碍物的密度在不同位置发生变化，压力和速度的计算仍能保持准确。他们还修正了“摩擦力”（黏度）的计算方式，确保无论流体移动快还是慢，无论是在海绵区域还是在开阔管道中，其表现都符合物理规律。

他们测试了什么

为了证明这种新的“通用规则手册”行之有效，他们进行了几项测试：

流动通道： 他们模拟了流体通过一个“海绵度”从一侧到另一侧平滑变化的管道。他们的模型与已知的物理特性完美契合。
多孔圆柱体： 他们模拟了风吹过一个由多孔材料（如网状球）制成的圆柱体。他们展示了如果球体的孔隙率很高，风会直接穿过它，从而平滑掉其后方的湍流；如果孔隙率较低，它就像一块固体岩石，会在后方产生旋涡状的气流。他们的模型捕捉到了所有这些细微的变化。
空腔： 他们模拟了一个带有移动顶盖（类似于滑动门）的盒子，并在其中放置了一个多孔障碍物。他们展示了盒子内部的旋转涡流如何根据障碍物的“海绵程度”而改变形状。

为什么这很重要

该论文声称，这种方法是早期设计阶段的强大工具。由于该模型具有极高的灵活性，并且能很好地处理变化的孔隙率，工程师可以用它来设计如下设备：

过滤器和分离器： 设计那些通过改变过滤器密度来优化流量的材料。
燃料电池： 理解气体如何在电池或燃料电池内部复杂的孔隙结构中移动。
热交换器： 设计流体通过填充床以高效传递热量的结构。

简而言之，他们构建了一个单一且稳健的计算机引擎，可以模拟流体通过从空管到致密海绵（以及两者之间任何状态）的一切物体中的流动，而无需切换模式或绘制人工界限。这使得设计依赖复杂流体流动的设备变得更加快速且准确。

技术摘要：用于非均匀多孔介质系统的单体动力学算法

问题陈述
在非均匀多孔介质中模拟流体流动，特别是耦合了自由流区域与多孔结构的系统，面临着显著的挑战。传统方法通常依赖于“双域”表述，即为自由流区域和多孔区域分别求解不同的控制方程，并需要显式的界面边界条件（例如 Beavers-Joseph 条件）来进行耦合。或者，“单域”方法虽然使用一套方程，但在处理非均匀介质时，历史上一直难以保持一致性。

现有的多孔介质格子玻尔兹曼方法（LBM）表述通常存在孔隙率缩放不一致的问题，即无法正确恢复宏观方程中的物理量。具体而言，许多模型在孔隙率随空间变化时，往往无法正确地对压力梯度和粘性应力项进行基于孔隙率的缩放。此外，许多表述是将 Brinkman 粘性应力表示为固有速度（intrinsic velocity）的梯度，而非表观速度（superficial velocity/Darcy velocity），这导致了在非均匀区域中的不准确性。此外，标准的 LBM 模型通常缺乏独立调节有效体粘度或声速的灵活性，从而限制了计算效率和稳定性。

方法论
作者提出了一种基于连续单域方法的单体、单域格子玻尔兹曼模型。该方法建立在以下支柱之上：

宏观模型： 目标流体力学通过体积平均法导出，得到一个广义的 Darcy-Forchheimer-Brinkman 模型。该模型适用于整个计算域，将自由流与多孔区域之间的过渡视为物理性质（孔隙率 $\epsilon$ ）的连续变化。控制方程包括：
- 体积平均密度（ $\epsilon\rho$ ）的连续性方程。
- 包含 Darcy 阻力和 Forchheimer 惯性阻力以及有效粘性应力张量（ $\sigma_{eff}$ ）的动量方程。
- 至关重要的是，有效粘性应力是根据表观速度（ $v$ ）而非固有速度（ $u = v/\epsilon$ ）的梯度来表达的，这确保了在孔隙率变化时与 Brinkman 表述的一致性。
动力学求解器： 该模型在标准晶格（如 D2Q9, D3Q19）上使用单一的广义动力学方程。
- 平衡态分布： 采用乘积形式的离散平衡态，以确保正确的矩和伽利略不变性。这包括特定的修正项（ $\Phi_{ext}$ 和 $\Phi_{\eta}$ ），用于恢复粘性应力张量的各向同性，并允许独立调节体粘度。
- 力项方案： 使用一致的广义力项方案（松弛至力法，Relax-to-Force method）来引入 Darcy-Forchheimer 阻力、体力和必要的修正项。
- 孔隙率缩放： 局部守恒场被定义为体积平均密度和动量，将孔隙率缩放自然地纳入动力学矩中。
- 速度评估： 提出了一个迭代算法，用于求解宏观速度与力项之间的隐式关系，从而处理非线性的 Darcy-Forchheimer 阻力和粘性修正。
边界条件： 作者将标准的湿节点边界条件（平衡态方案、非平衡外推和 Tamm-Mott-Smith 方案）应用于多孔介质表述。通过利用质量连续性计算边界处的固有流体密度，并考虑了力项和孔隙率变化的影响，推导出了特定的程序。

核心贡献

一致的单域表述： 本文提出了首个能够恢复具有完全一致孔隙率缩放（基于表观速度）的 Brinkman 粘性应力张量的 LBM 表述。这使得模型可以处理空间变化的孔隙率，而无需引入人工界面条件。
压力与粘性项的修正： 作者推导并实现了特定的力项修正，以纠正在非均匀介质中出现的压力梯度和粘性应力项的不一致性，确保模型能够恢复正确的宏观极限。
独立参数调节： 该模型允许将声速和有效体粘度与剪切粘度独立指定。这增强了计算效率（通过允许更高的声速）和数值稳定性（通过调节体粘度）。
单体算法： 该方法在整个定义域内使用单一的动力学方程和求解器，简化了涉及自由流与多孔区域复杂几何结构的实现过程。

结果与验证
模型通过一系列基准测试和数值研究进行了验证：

一致性分析： 线性时间谱分析证实，该模型能够恢复正确的有效声速以及剪切模态和声学模态的耗散率。分析表明，带有所提修正的乘积形式平衡态保持了伽利略不变性和各向同性，这与以往的二阶平衡态表述不同。
基准流场：
- 库埃特流（Coute Flow）： 在重力作用下的均匀多孔通道中，模拟在不同雷诺数和达西数下，其速度和密度分布均与参考数值解高度吻合。
- 泊肃流（Poiseuille Flow）： 模型准确捕捉了均匀通道以及关键的非均匀孔隙率剖面通道中的流动。它正确预测了流体倾向于流经高孔隙率区域的流道效应（channeling effects）。
数值演示：
- 流经多孔圆柱： 模型成功模拟了受限流经多孔圆柱的过程，并根据雷诺数和核心孔隙率识别出五种不同的流动机制。它捕捉到了诸如尾迹稳定、脱离再循环以及远场尾迹不稳定性等现象，显示出良好的定性一致性。
- 带多孔障碍物的驱动顶盖流： 模型模拟了一个带有内部多孔障碍物的腔体，展示了渗透率如何影响涡旋形成、流动渗透和再循环模式。结果强调了随着孔隙率增加，流场从类固体行为向类渗透行为的转变。

意义与主张
作者声称，这项工作为模拟非均匀多孔介质中的传输提供了一个鲁棒且统一的框架。通过消除对显式界面边界条件的需求，并确保所有流体力学项的缩放一致性，该模型特别适用于：

工程多孔结构（具有空间变化的孔隙率）的早期、器件级设计研究。
模拟自由流与渗透体之间复杂的相互作用（例如燃料电池、过滤系统、多孔燃烧）。
捕捉自由流与多孔介质之间过渡连续或界面难以精确解析的系统中的传输现象。

论文结论指出，虽然目前的研究侧重于二维验证，但该表述本质上是多维的，并已准备好扩展到三维几何结构。未来的工作需包括引入能量平衡（热传递）、多组分混合物以及多相流，以全面解决现实工程系统中耦合传输的问题。

技术摘要：用于非均匀多孔介质系统的单体动力学算法

类似论文