🔢 mathematics

Development and application of a multiphase Lagrangian structure function model in anisotropic turbulence

本研究利用直接数值模拟并推导出一个精确的守恒律，旨在研究壁面湍流中的流动各向异性如何影响沉降惯性颗粒的尺度级能量学，揭示了颗粒惯性和沉降如何将大尺度各向异性保持至各向同性范围内，并为未来的连续介质模型提供参考。

原作者： Andrew P. Grace, David Richter

发布于 2026-07-07

📖 1 分钟阅读🧠 深度阅读

原作者： Andrew P. Grace, David Richter

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

大局观：风中的尘埃

想象一个充满旋转气流（湍流）的房间。现在，想象你扔进去一堆沉重的弹珠（颗粒）和一把轻飘飘的小羽毛。

空气： 空气以混乱、旋转的方式运动。靠近地板的地方，空气的表现与靠近天花板的地方不同，因为地板阻碍了流动。这被称为各向异性（意味着规则会根据你观察的方向而改变）。
颗粒： 重弹珠具有惯性。它们不会瞬间跟随空气运动；它们会保持原有的运动方向，就像汽车试图转弯但由于惯性而发生侧滑一样。而轻盈的羽毛则会完美地跟随空气。

这篇论文提出了一个具体的问题：如果靠近地板的空气是“被挤压”且具有方向性的（各向异性的），那么当这些弹珠进入那些本该是完全圆润且随机的微小、混乱旋流时，它们还会记得那种形状吗？

主要发现：“记忆”效应

研究人员发现：是的，重弹珠确实记得。

即使当空气旋流变得如此之小，以至于看起来是完全圆润且随机的（各向同性的）时，重弹珠的运动方式仍然反映了它们所来自的较大规模气流的“挤压”形状。这就像弹珠带着来自地板处那种大尺度、有方向性风力的“记忆”，一直向上进入了微小的随机旋流中。

轻颗粒（羽毛）： 它们会迅速忘记大局。如果空气是随机的，它们的运动也是随机的。
重颗粒（弹珠）： 它们保留了“大风”的形状。如果大风是水平拉伸的，那么即使在微小的随机旋流中，弹珠也会保持水平运动。

他们是如何研究的

作者们并非仅仅靠猜测；他们构建了一个超级计算机模拟（一个“虚拟风洞”）来观察这一现象。

设置： 他们模拟了一面墙（地板）以及流经它的气流。他们将颗粒丢入这种流体中，并让重力将它们向下牵引。
测量： 他们使用了一种叫做结构函数的工具。你可以把它想象成一个“能量距离尺”。
- 他们观察间隔一定距离的一对颗粒。
- 他们问道：“如果我观察两颗靠得很近的颗粒，它们的运动差异有多大？那么两颗离得很远的颗粒呢？”
- 通过在不同距离下进行测量，他们可以观察能量是如何从大旋流传递到小旋流的。

两个关键因素

论文强调了改变颗粒行为的两个主要因素：

惯性（颗粒有多重）：
- 如果颗粒很轻，它的行为就和空气一致。
- 如果颗粒很重，它就会忽略空气中微小、快速的变化，而是继续做它之前正在做的事情。这就是为什么它们能保留“大风”的记忆。
沉降（重力）：
- 颗粒受重力作用而下落。研究人员观察了在重力较弱与重力较强的情况下会发生什么。
- 令人惊讶的是： 当重力非常强时，它实际上有助于颗粒“重置”并重新变得更加随机（各向同性），尤其是在流体中间区域。但在靠近地板（墙壁）的地方，重颗粒仍然难以摆脱墙壁带来的方向性影响。

数学“配方”

作者还写出了一套新的规则（方程）来描述这一过程。

想象你在预测天气。你有一套完美的空气方程，但当你加入重弹珠时，方程变得复杂了，因为弹珠不仅仅是跟随空气，它们还带着自己的历史信息。
作者推导出了一个“守恒定律”。你可以把它看作是一个“能量收支表”。它追踪不同尺寸下的能量情况。
他们发现，颗粒的“能量预算”与空气不同，因为颗粒具有这种记忆效应。它们不仅对身边的空气做出反应，还会对刚才经过的空气做出反应。

核心结论

在靠近墙壁的湍流中：

空气在极小的尺度上最终会变得完全随机且圆润。
重颗粒并不会变得完全随机。由于惯性的作用，它们会保持被“拉伸”的状态，并朝着大风的方向运动，从而携带了大尺度流动的记忆。

这很重要，因为如果你想预测现实世界中（存在墙壁和重力的情况下）灰尘、雨滴或工业颗粒的运动，你不能仅仅假设它们的行为与空气一致。你必须考虑到它们对大尺度、有方向性风力的“记忆”。

技术摘要：各向异性湍流中多相拉格朗日结构函数模型的开发与应用

问题陈述
惯性颗粒对湍流的能量响应对于理解多相动力学以及开发颗粒扩散的降阶模型至关重要。虽然已有大量研究探讨了均匀各向同性湍流（HIT）中的颗粒能量学，但实际流动通常表现出空间非均匀性和各向异性，特别是在靠近固体边界处。在壁面湍流（如湍流边界层，TBL）中，大尺度各向异性是由平均剪切和边界条件产生的。一个关键的未解问题是，这种大尺度各向异性是如何传递到颗粒相的，以及它是否会在形式上各向同性的小尺度范围内持续存在。此外，在这些环境中，颗粒惯性、重力沉降与流体各向异性之间复杂的耦合关系需要进一步研究，以准确模拟惯性沉降颗粒。

研究方法
作者采用结合理论与数值的方法，研究了在净沉积条件下湍流边界层中的尺度相关颗粒能量学。

理论框架：
- 研究利用动力学理论推导了任意尺度下颗粒相二阶结构函数（SoSF）的精确（尽管是不封闭的）守恒定律。
- 颗粒动力学使用 Maxey-Riley 方程进行建模，针对大颗粒-流体密度比极限下的微小、球形、无相互作用颗粒。
- 作者定义了一个基于纵向与垂直二阶结构函数差异的颗粒相各向异性度量 $\eta^{(p)}$ 。
- 通过对推导出的连续性方程进行渐近分析，在假设相对速度矩和径向分布函数遵循幂律行为的前提下，阐明了控制能量通量和各向异性的各项之间的平衡。
数值模拟：
- 对单向耦合颗粒载流壁面湍流进行直接数值模拟（DNS），摩擦雷诺数取中等水平（ $Re_\tau = 315$ ）。
- 模拟范围涵盖了湍流边界层的对数层和缓冲层。
- 研究通过改变 Stokes 数（ $St^+$ ）来评估颗粒惯性，并通过沉降参数（ $Sv^+$ ）来评估重力漂移的影响。
- 同时利用欧拉（流体）和拉格朗日（颗粒）视角，并特别关注“片状”（slab-wise）统计，即垂直分离度为零的情况，以隔离局部垂直动力学。

主要贡献

守恒定律的推导： 本文推导出了任意尺度下颗粒相能量学的精确守恒方程，明确识别了传输项、与一阶矩的耦合项以及能量通量项。这为分析各向异性多相流中的尺度相关能量传递提供了一个连续介质框架。
各向异性传递的量化： 本工作系统地量化了载流流体的大尺度各向异性是如何“捐赠”给颗粒相的，以及颗粒惯性和沉降如何修正这一过程。
渐近分析： 作者通过渐近分析来解释 DNS 数据，将观察到的各向异性度量的非单调行为与大尺度流体各向异性与小尺度内尺度间能量传递之间的竞争联系起来。

结果

惯性与各向异性： 在小 Stokes 数时，颗粒相在小尺度处表现出接近各向同性的特征，与流体一致。然而，随着惯性增加（ $St^+ = 10$ ），颗粒相即使在小尺度下仍保持显著的各向异性，表明颗粒惯性保留了大尺度流体各向异性的记忆。
沉降效应： $Sv^+$ 的影响高度依赖于流场区域。在缓冲层（ $z^+ = 30$ ，其底层流体本质上是各向异性的）中，增加 $Sv^+$ 会导致结构函数的复杂变化，表现为垂直分量增加而纵向分量减少。在对数层（ $z^+ = 90$ ，其小尺度流体是各向同性的）中，高沉降率（ $Sv^+$ ）倾向于驱动颗粒相在小尺度下趋向于各向同性状态。
尺度非局域性： 与预期的流体 $r^2$ 标度相比，颗粒结构函数在耗散区表现出较浅的标度（趋向于 $r^{2/3}$ ）。这反映了惯性颗粒的尺度非局域性，该特性主要由惯性驱动，并受沉降调节。
各向异性度量行为： 颗粒各向异性度量 $\eta^{(p)}$ 随分离尺度 $r$ 的减小呈现非单调行为。在 $r/z \to 0$ 时，它在达到一个局部极小值后再次上升。这表明大尺度流体施加的各向异性与颗粒相内部尺度间能量传递产生的各向异性之间存在竞争。

意义与主张
本文声称提供了关于各向异性湍流中惯性颗粒尺度相关能量通量的关键见解。通过证明颗粒惯性可以将大尺度各向异性保留至形式上各向同性的小尺度，作者指出了假设颗粒相具有局部各向同性的模型的局限性。推导出的守恒定律及识别出的各向异性传递机制，对于开发惯性沉降颗粒的连续介质模型具有直接意义，特别是关于双向耦合动力学如何受到这些非局域各向异性效应影响的问题。本研究强调，为了实现准确建模，必须同时考虑颗粒惯性和重力漂移，因为它们的影响并非简单的线性叠加，而是在颗粒动力学中存在隐式耦合。