《Journal of Fluid Mechanics》：Taylor dispersion in long straight channels of arbitrary and gradually changing size and shape
编辑推荐：
中性溶质在流动中的平流分散与许多物理系统的设计相关，包括微流控装置和化学分离系统。这些装置中通道的横截面形状通常很复杂，并且可能沿流动方向变化（形状或面积）。然而，先前对变几何通道中泰勒-阿瑞斯分散（Taylor–Aris dispersion）的分析通常集中
中性溶质在流动中的平流分散与许多物理系统的设计相关，包括微流控装置和化学分离系统。这些装置中通道的横截面形状通常很复杂，并且可能沿流动方向变化（形状或面积）。然而，先前对变几何通道中泰勒-阿瑞斯分散（Taylor–Aris dispersion）的分析通常集中在周期性通道或具有简单横截面的通道上。据研究人员所知，目前尚不存在用于预测任意横截面（其形状或尺寸沿轴向缓慢变化）通道中分散的降阶理论。在当前研究中，研究人员开发了一个适用于分析具有此类形状通道中分散的降阶模型理论框架，适用于稳态或时变流动条件。研究人员首先推导了一个关于面积平均溶质浓度的一维偏微分方程（PDE），以及完整的三维（3-D）示踪剂分布的相关表达式。接下来，研究人员建立了两个耦合的常微分方程（ODE）来控制溶质轴向平均位置和方差。此外，研究人员推导了一个常微分方程，用于设计通道以产生先验指定的溶质方差分布。研究人员将模型应用于变半径轴对称通道、变“哑铃”横截面通道和变宽度矩形通道的示例几何形状。研究人员将模型性能与所有测试几何形状的布朗动力学模拟进行了基准比较。总体而言，研究人员的框架易于应用于不同的几何形状和流动条件，并可能为微流控系统的设计和优化提供关键见解。
**论文解读：在任意缓慢变化截面长直通道中的泰勒分散降阶模型**
**研究背景与问题**
泰勒-阿瑞斯分散（Taylor–Aris dispersion）是描述中性溶质在流动通道中因平流与横向扩散耦合而导致的轴向展宽现象，对于微流控装置、化学分离系统及环境流动等物理系统的设计具有基础重要性。传统泰勒分散理论（G.I. Taylor, 1953；Aris, 1956）适用于均匀截面通道，而实际装置中的通道截面往往形状复杂，并可能沿流向变化。尽管已有研究将泰勒分散分析扩展至变几何通道，但大多局限于周期性通道或简单截面（如轴对称变半径管道）。同时，广义泰勒分散理论（macrotransport paradigm）虽能处理周期性系统，但局限于长时间渐近行为。目前尚缺乏一个统一的降阶模型，能够预测任意缓慢变化截面（形状、尺寸或两者兼有）通道中的分散行为，并允许在有限时间尺度内进行工程设计。因此，研究人员开展了本研究，旨在构建一个普适且易于应用的简化理论框架，以填补这一空白。该论文发表在《Journal of Fluid Mechanics》。
**主要技术方法**
研究人员采用面积平均与偏差量方法，引入一个新的“溶质偏差源场”（solute deviation source field, Ψ），该场通过求解通道横截面上的泊松方程（Poisson equation）确定，其源项由轴向速度偏差和截面面积的局部相对变化率构成。基于此，研究人员推导出描述面积平均浓度的一维偏微分方程（PDE），以及用于重构完整三维浓度场的表达式。在溶质云团宽度远小于通道几何变化轴向尺度（γ?1）的极限下，进一步推导出两个耦合的常微分方程（ODE），控制溶质轴向平均位置和方差，并得到一个用于工程化设计通道以产生指定方差分布的非线性一阶ODE。模型验证采用布朗动力学模拟（Brownian dynamics simulations），对三种不同几何形状的通道进行基准比较。
**研究结果**
**3.1 变半径轴对称通道**
研究人员将模型应用于压力驱动流（PDF）中半径缓慢变化的轴对称通道。通过求解Ψ场，得到一维PDE（方程3.7）和ODE（方程3.9-3.10），并证明该框架在γ?1极限下可简化为Chang & Santiago (2023)的模型。在Péclet数（Pe）为10、100和1000时，将PDE模型与布朗动力学模拟进行对比，发现两者在溶质方差随平均位置变化上表现出极好的一致性。即使当γ>1时（溶质云团跨越多个周期性几何周期），PDE模型仍保持准确，因为溶质在多个周期内采样了局部分散率，使得方差增长近似线性。
**3.2 哑铃横截面几何**
研究人员分析了一个沿流向逐渐变细的哑铃形截面通道（手柄高度沿轴向线性减小）。通过计算流体动力学得到三维速度场，数值求解Ψ场，并利用PDE和ODE模型预测分散。布朗动力学模拟显示，在γ最大约0.15的条件下，ODE和PDE模型均与模拟结果高度吻合。值得注意的是，尽管局部ηPe（η为扩散长度与溶质云宽度之比）在通道末端超过1，模型仍准确预测分散。进一步分析三维浓度场发现，溶质在哑铃手柄区域速度显著滞后，导致Ψu′乘积为较大负值，从而增强有效扩散系数（D*eff = 1 - ?Ψu′?/D），表明狭窄连接结构对剪切分散有显著贡献。
**3.3 变宽度矩形通道**
研究人员将模型应用于恒定高度、宽度（即长宽比）沿轴向变化的矩形通道，并仅考虑γ?1的极限。通过解析求解Ψ场，得到有效扩散系数表达式。首先，绘制了零方差增长曲面（方程3.24），表明在固定Péclet数下，使方差增长为负所需的局部壁面斜率随通道长宽比增大而减小。其次，利用逆向工程ODE（方程2.31）设计矩形通道，以实现恒定的、正弦变化的以及正弦加线性增长的溶质方差分布。布朗动力学模拟验证了这些工程化通道的方差预测，结果与目标值吻合良好。
**总结与结论**
研究人员开发并展示了一个通用的泰勒分散分析框架，适用于任意截面、沿轴向缓慢变化的长直通道。通过引入溶质偏差源场Ψ，将三维对流-扩散方程简化为一个一维PDE，并进一步在溶质云团局部化假设下得到两个耦合ODE，用于预测溶质平均位置和方差。该框架还允许逆向工程，即通过求解ODE设计通道以产生指定的溶质方差分布。模型在三个截然不同的几何形状（轴对称、哑铃形、矩形）中均与布朗动力学模拟验证一致，并成功预测了三维浓度场。主要假设包括η?1、c0′/c0?1、γ(c0′/c0)?1、ηPe?1及St(c0′/c0)?1，其中ηPe是关键参数。该模型适用于多种流动条件（如压力驱动流、电渗流等），可望为微流控系统的设计与优化提供关键见解。
