如何借助“粒子追踪模块”，开发分析静态混合器性能的App-电子发烧友网

静态混合器因其高效、成本低廉、易于安装且维护要求低等优点，成为各类工程领域的常用工具。在评估混合器能否满足某种使用目的时，一个重要的判断指标是得到的混合物是否足够均匀。在本篇文章中，我们将介绍如何借助“粒子追踪模块”，开发一款能定量和定性分析静态混合器性能的 App。

层流静态粒子混合器设计器 App 的设计基础

在创建“层流静态粒子混合器设计器”App 前，我们首先前往官网的“案例下载”页面下载层流混合器教程。此模型通过计算贯穿设备的粒子轨迹对静态混合器的混合性能进行了评估。如果您希望深入了解此教程及学习基本的混合器建模知识，我们推荐您阅读下列文章：

如何模拟层流静态混合器内的粒子追踪

层流静态混合器建模

静态混合器建模

上述教程使用了与我们的 App 相同的几何模型。如下图所示，模型由一段装有三个螺旋叶片的管道构成，叶片可交替旋转。搅拌叶片的表面呈灰色，周围勾勒出管道的轮廓线。流体运送粒子通过管道时，静态搅拌叶片会使粒子混合在一起。

层流静态混合器模型的几何结构。

用于研究静态混合器性能的 App

借助下图中的层流静态粒子混合器设计器 App，我们可以先计算出通过混合器的粒子轨迹，然后利用一些内置后处理工具定量及定性地对混合器的性能进行评估。

“层流静态粒子混合器设计器”的用户界面（UI）截图。

此 App 中包含大量的几何参数和材料属性，在创建混合器模型时，可选择使用一个、两个或三个螺旋式搅拌叶片。除此之外，通过 App 的高级设置，您还能更改粒子数量和后处理参数。

为了更准确地绘制静态混合器内不同物质的分布及变化，我们可以释放粒子，并使用流体流动粒子追踪接口对其轨迹进行计算。粒子位置是通过牛顿运动方程公式计算出来的，而计算位置矢量的分量则需求解一组二阶方程：

其中

在模型中，粒子只受到曳力的作用，并使用 Stokes 曳力定律计算得到：

其中应用了下列物理量：

Stokes 曳力定律适用于相对雷诺数远小于 1 的粒子；也就是，

其中

静态混合器内粒子轨迹的绘图。

借助 App 开发器量化静态混合器的性能

肉眼观察在一定程度上可以判断混合物的均匀程度，针对本文的案例，混合性能的可视化能够通过创建粒子位置的相图来实现。在相图中，我们可以在任意的二维相空间中绘制粒子轨迹——也就是说，可以通过用户自定义的表达式在二维绘图的坐标轴上表示粒子轨迹。相图常用于绘制特定方向上粒子位置随动量的变化，即相空间分布。

在下方动画中，相图用于观察每个粒子在通过混合器的过程中横向位置的变化。管道设为 y 方向，因此横向为 x 和 z 方向。我们使用不同的颜色来表示初始状态下每种粒子所占据的象限。从图中可以看出，初始时释放的深蓝色粒子位于x和z轴的正方向。

图显示了粒子在通过混合器过程中的横向位置。

相位图定性地表面出口处的粒子没有完全均匀地混合。图中仍存在粒子数量密度偏高或偏低的区域，同时还能看到相同颜色的粒子（来自同一象限）聚集在了一起。

相图存在一个潜在的缺陷，这是因为它绘制的是同一时刻相空间内的粒子，而非基于相同y坐标的粒子。生成的混合器可视化绘图也因此带有一定的误导性，这因为一些更靠近叶片的粒子可能会晚于其他粒子到达出口。因此我们提出了另一种方法：创建庞加莱映射，该映射可绘制出粒子轨迹与特定位置上平面的交点。

下图中，每个截面上粒子的颜色取决于释放时该粒子的初始x坐标，蓝色表示初始坐标为正方向，红色表示为负方向。接下来我们再一次观察出口处红蓝粒子的聚集情况。

庞加莱映射显示了二维绘图中粒子的位置。

我们可以从相图与庞加莱映射中获取大量有关混和器性能的信息，然而对于精密的工业应用而言，多数信息过于主观。肉眼观察可以大概判断不同物质处于完全分离、部分混合或充分混合的状态，但这些说法太模糊，难以量化。举例来说，任何一个观察者都能看出上图中聚集着一小群同色粒子，但要指定数值来描述粒子混合程度就十分困难了。

幸运的是，“App 开发器”和“方法编辑器”提供的工具可以帮助我们创建专业的高端后处理程序，这些程序可以对特定混合器几何模型的性能进行赋值。分布指数是一个常用于评估粒子空间分布均匀性的度量标准，它被定义为方差和平均数之比：