关于自由射流中类孤立波相干结构的数值研究
一、 研究作者与发表信息
本研究的主要作者为胡子渊，其所属机构为北京大学力学与工程科学学院。该研究是其提交给北京大学的博士学位论文，导师为蔡庆东副教授。论文完成并提交于2025年12月。
二、 学术背景与研究目标
本研究隶属于流体力学领域，具体聚焦于湍流基础研究，特别是自由射流中的相干结构。自由射流作为一种广泛存在于自然界（如火山喷发）和工程应用（如喷气推进、燃烧）中的典型流动形态，其湍流结构的生成与演化机制一直是该领域的核心课题。以往的研究多从欧拉视角出发，重点关注射流中的涡结构，而相对忽视了其中的波状结构。
近年来，李存标等研究者在平板湍流边界层转捩过程中发现了一种新的流动结构——类孤立波相干结构（solitons-like coherent structures, SCS）。这种结构被认为在湍流的产生过程中可能具有普适性，但此前在自由射流这一重要流动形态中尚未得到系统研究。鉴于此，本研究旨在填补这一空白，其核心目标是：采用数值模拟与拉格朗日方法相结合的研究手段，系统探究不可压缩圆管自由射流中是否存在SCS，并深入分析其演化行为、运动学特性、生成机理以及对流动条件的响应规律，从而为理解自由射流乃至更广泛的湍流转捩与湍流维持机制提供新的视角和见解。
三、 详细研究流程与方法
本研究采用了严谨的计算流体力学研究范式，其工作流程主要包含以下几个关键步骤：
第一步：数值模型建立与验证
研究首先建立了圆管自由射流的数值模型。控制方程为不可压缩Navier-Stokes方程。求解器基于开源计算流体力学软件OpenFOAM，采用有限体积法进行数值求解。研究设置了合理的计算域和边界条件以模拟真实的自由射流发展过程。为确保模拟结果的可靠性，研究进行了严格的验证工作：将数值模拟得到的流场统计特性（如速度剖面）与理论解以及自由射流固有的自相似性特征进行对比。验证结果表明，模拟结果与理论和已知物理规律高度吻合，从而确认了所建立数值模型的准确性和可靠性，为后续的相干结构分析奠定了坚实基础。
第二步：流场特征分析与SCS初步识别
在获得可靠的流场数据后，研究首先从传统的欧拉视角进行分析。通过分析脉动速度场，研究者识别出射流中交替分布的高速区和低速区，这是湍流中相干结构的典型宏观表现。然而，仅凭欧拉速度场难以精确追踪结构的时空演化。因此，研究引入了拉格朗日示踪法。该方法通过追踪流场中特定“物质线”和“物质面”（即由相同流体质点构成的线和面）的变形与运动，能够更直观地揭示流动结构的拉格朗日演化过程。通过对物质线和物质面的演化进行追踪，研究首次在自由射流中确认了SCS的存在，并初步观察到了其三维波状形态。
第三步：SCS运动学特性的定量测量
为了深入理解SCS的本质，研究对其运动速度进行了定量测量。这里采用了相关性分析的方法。具体步骤是：在识别出的SCS区域，选取特定的流动变量（如速度）作为参考信号，通过计算该信号在不同空间位置、不同时间点之间的相关性（如皮尔逊相关系数或斯皮尔曼相关系数），可以确定波结构的传播速度和方向。这一方法的关键在于，如果结构是随流平动的涡，其传播速度应接近当地流速；而如果是具有独立传播特性的波，其波速则可能与当地流速存在显著差异。通过这一精细的定量分析，研究得以将SCS的运动学特征与传统的涡结构区分开来。
第四步：SCS演化机理与涡生成机制探究
在确认SCS存在及其波性质后，研究进一步探究其演化动力学和与旋涡生成的关系。首先，进行了通量分析。通过计算和分析SCS结构周围各个方向（特别是周向）的质量或动量通量，研究试图揭示驱动SCS径向（即垂直于射流轴线方向）拱起或发展的力学原因。其次，为了明确SCS与旋涡生成之间的内在联系，研究采用了拉格朗日平均涡量偏移（Lagrangian-averaged vorticity deviation, LAVD） 方法进行旋涡识别。LAVD是一种基于拉格朗日观点的涡识别准则，能够有效追踪流体微团旋转强度的历史积累，从而更准确地识别出旋涡核心区域。通过将LAVD识别的涡量聚集区与SCS引起的剪切层区域进行对比，可以直观地揭示两者在空间上的关联性。
第五步：雷诺数效应研究
为了考察SCS的普适性及其对流动条件的依赖性，研究系统性地改变了流动的雷诺数（Re），在2000至6600的范围内进行了多组数值模拟。对于每一组雷诺数下的模拟结果，重复上述第二至第四步的分析流程，特别是通过拉格朗日物质面追踪，观察和统计SCS的生成位置、数量、结构形态、失稳破碎过程以及波速等关键参数随雷诺数的变化规律。这一系列对比研究旨在揭示SCS演化行为中的普适性特征和受雷诺数影响的特异性规律。
四、 主要研究结果
本研究取得了一系列系统而深入的发现，各步骤结果环环相扣，共同支撑了最终结论。
1. SCS的确认与形态特征：
通过拉格朗日物质线与物质面的追踪，研究明确证实了在圆管自由射流中存在类孤立波相干结构（SCS）。研究识别出两种典型但形态相反的三维波结构：一种位于低速区，其在流向（沿射流轴线方向）减速的同时向射流中心轴方向（向内）拱起；另一种位于高速区，其在流向加速的同时向远离轴线的方向（向外）拱起。这两种结构都伴随着显著的剪切层和旋涡的生成。值得注意的是，自由射流中观测到的SCS演化规律，与之前在平板湍流边界层中观测到的结果在定性上一致。这一发现强烈暗示，SCS可能是流体从层流转捩到湍流过程中一种具有普适性的流动结构，而不仅限于边界层流动。
2. SCS的类孤立波性质：
通过相关性分析对SCS传播速度的定量测量得到了关键数据。结果显示，低速区SCS的波速约为当地中心线速度的5%至35%，而高速区SCS的波速约为当地流速的70%至90%（相当于中心线速度的25%至40%）。这一发现至关重要：SCS的传播速度与当地流体流速存在显著差异，且其传播方向具有优异的稳定性。这直接证实了SCS具有“类孤立波”的性质——即它是一种能够保持形状、以特定速度独立传播的波动结构，而非简单地被流体裹挟着运动的涡旋结构。这是本研究区分SCS与传统涡结构的核心证据。
3. SCS的演化机理与涡生成机制：
通量分析的结果表明，SCS在径向的发展（即拱起）主要是由其两侧的周向通量差所驱动。形象地说，是SCS两侧剪切层对结构的“挤压”效应，导致了其进一步的拱起变形。这揭示了SCS三维形态演化的一个内在动力来源。
更重要的发现来自LAVD分析。结果显示，涡量的高度聚集区与SCS所诱导产生的高剪切层区域在空间上高度重合。这清晰地揭示了SCS诱导产生旋涡的机理：SCS的演化过程会形成局部的高剪切层；这些高剪切层不仅会进一步驱动SCS结构本身的变形与发展，同时其不稳定性会直接导致旋涡的生成。这一发现将波状结构（SCS）与涡结构（湍流的核心要素）的生成在物理机制上联系了起来。
4. SCS的雷诺数效应：
对不同雷诺数下SCS演化的系统研究表明，其行为既存在普适特征，也受雷诺数显著影响。
- 普适特征：在不同雷诺数下，SCS都表现出一些共同规律。例如，在径向不同高度（r/d）的物质面上首次观测到SCS的流向位置（z/d）与该径向高度呈正相关（即越靠外，SCS出现得越晚）；而在同一物质面上观测到的SCS数量则与径向高度呈负相关（即越靠近外围，结构越少）。此外，SCS的径向扩展速率始终大于射流主体的径向扩展速率，且其波速在整个研究的雷诺数范围内都保持相对稳定的性质。
- 雷诺数影响：随着雷诺数增大，SCS的整体生成位置会向上游（更靠近喷口）移动，其数量增多，并且结构会更早地发生失稳和破碎。这表明更高的流动强度（更高Re）促进了SCS的生成，但也加速了其非线性演化和崩溃过程。一个特别有趣的发现是，在特定雷诺数（Re=3500）下，外部流场的卷吸效应会抑制SCS结构的生成，这揭示了复杂流动环境对SCS动力学的微妙影响。
五、 研究结论与价值
本研究通过高精度的直接数值模拟和创新的拉格朗日分析方法，首次系统性地在圆管自由射流中发现了类孤立波相干结构（SCS），并对其运动学、动力学特性及演化机理进行了深入研究。
科学价值：
- 拓展了相干结构的研究范畴：将SCS的存在从平板边界层流动证实并拓展至自由射流这一重要流动体系，增强了“SCS是转捩过程中一种普适性流动结构”这一论点的说服力。
- 揭示了新的湍流生成物理图像：研究明确了SCS具有独立于当地流速的传播特性（类孤立波性质），并阐明了其通过形成高剪切层来诱导旋涡生成的具体机理。这为理解湍流中“波”与“涡”的相互作用、以及湍流能量级串的初始环节提供了新的物理模型和见解。
- 提供了系统的参数化认识：通过研究雷诺数效应，量化了SCS的生成、演化与流动条件之间的关系，为预测和控制在类似流动中的相干结构行为提供了参考依据。
应用价值：对自由射流中SCS的深入理解，有助于发展更先进的流动控制策略。例如，通过主动干预SCS的生成或演化，可能实现对射流混合效率、噪声产生或燃烧稳定性的有效调控，在航空航天推进、化工混合、环境流体力学等领域具有潜在应用前景。
六、 研究亮点
- 研究视角新颖：突破了自由射流研究中长期侧重欧拉视角和涡结构的传统，率先从拉格朗日视角系统探究其中的波状相干结构，获得了全新的发现。
- 方法组合创新：创造性地将直接数值模拟（DNS）与拉格朗日示踪法相结合，并综合运用相关性分析、通量分析和LAVD涡识别等定量工具，形成了从识别、表征到机理分析的一套完整研究方法论。
- 发现具有突破性：首次在自由射流中定量证实了SCS的类孤立波性质，并清晰揭示了其诱导涡生成的物理链条（SCS → 高剪切层 → 旋涡），这是对湍流相干结构物理本质认识的重要深化。
- 研究系统全面：不仅进行了个例的深入分析，还通过系统的参数研究（雷诺数效应），揭示了SCS行为的普适规律和特异性，使结论更加坚实和全面。
七、 其他有价值内容
本研究的摘要、目录及绪论部分显示，论文对相干结构的研究历史、自由射流的基本特征与研究历程进行了详尽的文献综述，为该项研究的立项提供了坚实的学术背景支撑。此外，论文详细阐述了所使用的数值算法（如PIMPLE算法）、网格无关性验证过程以及数据后处理方法，体现了计算流体力学研究的高度严谨性和可重复性，其方法学本身也对后续相关研究具有参考价值。