关于“网格生成湍流再探究”的学术研究报告

托马斯·库兰（Thomas Kurian）与延斯·H·M·弗兰松（Jens H. M. Fransson）两位研究人员，来自瑞典斯德哥尔摩皇家理工学院（KTH）的Linné流动中心力学系，于2009年3月6日在学术期刊《流体动力学研究》（*Fluid Dynamics Research*， 卷41， 文章号021403）上发表了一篇题为“网格生成湍流再探究”（Grid-generated turbulence revisited）的受邀论文。该研究深入表征了新型湍流生成网格的特性，并对网格湍流的经典理论和方法进行了批判性审视与更新。

一、 学术背景 本研究主要隶属于流体力学，特别是实验湍流研究领域。网格生成湍流是湍流理论研究和模型验证的经典测试案例，已有超过七十年的研究历史。其实验室应用广泛，既可用于改善风洞流动品质，也可用于生成受控的湍流场。研究的核心动机源于对湍流特性如何影响后续物理过程（如边界层转捩）的深入理解需求。前人研究表明，自由来流湍流不仅其强度（Turbulence Intensity, Tu）对边界层转捩位置有重要影响，其特征长度尺度（如积分尺度、泰勒微尺度）同样至关重要。然而，在实验上，独立地改变湍流强度与特征长度尺度极为困难。通常，改变网格几何尺寸会同时改变这两个参数。因此，本研究旨在通过设计一组具有大致相同固体率（solidity， 即网格阻挡面积比例）但不同网眼宽度和杆直径的网格，尝试实现“在保持湍流强度大致不变的前提下，系统改变湍流特征长度尺度”这一实验目标，这为后续研究（如湍流诱导转捩）提供了关键工具。此外，论文还指出，网格生成的湍流常被错误地描述为各向同性，但实际上各向异性普遍存在。许多基于各向同性湍流假设推导的关系式（如耗散率、长度尺度计算公式）在实际应用时可能产生高达50%或更多的误差。因此，本研究另一个重要目标是利用详尽的实验数据，对比和验证各向同性与各向异性关系式的准确性，并系统地表征新设计网格的湍流衰变、长度尺度、能谱及耗散率等特性。最后，研究还初步探讨了顺压梯度（Favorable Pressure Gradient, FPG）对网格湍流特性的影响。

二、 详细研究流程 本研究是一项系统的实验研究，主要流程包括实验设施与网格准备、速度与压力测量、数据采集以及多层次的数据分析与验证。

1. 实验设施与研究对象：研究在斯德哥尔摩KTH的“最低湍流度”（Minimum Turbulence Level， MTL）风洞中进行。该风洞品质极高，背景扰动水平低，温控精确，并配备五自由度自动探针 traversing 系统。研究对象为七块湍流生成网格，其中五块是新设计的低湍流强度网格（LT1-LT5），由圆形细金属丝编织而成；另外两块（网格A和E）是较早研究的网格，用于对比和展示雷诺数效应。网格A为圆形杆编织网格，网格E为方形杆单平面网格。这些网格的关键设计在于其固体率（σ）大致相同（约0.38-0.44），但网眼宽度（m）和杆直径（d）不同（m/d 比值在4.0到6.0之间），这正是实现“变尺度、保强度”设想的核心。所有网格都安装在一个可精密移动的通用框架上，便于在风洞中更换和定位。

2. 测量方法与设备：速度测量采用热线风速仪（Hot-Wire Anemometry）。使用了两种探头：单丝探头用于获得更高的频率响应（36 kHz）；X型探头用于同时测量两个方向的脉动速度分量。通过一个特别设计的、可手动旋转90度的探头支架，研究者能够先后使用两个X型探头（一个测量u-v分量，旋转后另一个测量u-w分量）来近乎同步地获取所有三个方向（u， v， w）的脉动速度信息，这是获取完整湍流张量信息的关键。探头为自制，使用铂金丝，并进行了详细的 King‘s 律和角度校准。压力测量则使用微压差计测量网格前后的静压差，以确定网格的阻力系数。所有测量均在风洞测试段中心线进行，针对不同网格在不同来流速度下（LT网格：6， 9， 12 m/s； A/E网格：2， 4， 8， 12， 16， 30 m/s）采集了下游多个位置（以 x/m 表示）的数据。

3. 数据分析流程：数据分析是本研究工作量最大、也最具方法论贡献的部分。对于每个测量工况，研究者系统地计算了以下参数：

   • 各向异性度：通过比较三个方向脉动速度的均方根值（u_rms， v_rms， w_rms）来衡量。
   • 湍流强度衰变：分析湍流动能（k = 0.5*(u_rms^2 + v_rms^2 + w_rms^2)）沿下游的衰变。研究开发了一套客观的、与曲线拟合起始位置无关的流程来确定衰变幂律公式（t_u^2 = a * [(x - x0)/m]^b）中的系数a、b和虚拟原点x0。该流程通过系统性地改变拟合起始点，观察b、a、x0的收敛平台来确定最终值，避免了主观选择拟合区间带来的误差。
   • 湍流特征长度尺度：
     • 积分尺度（Λ）：通过两种独立方法计算：(i) 直接积分纵向空间自相关函数f®（利用泰勒冻结湍流假说将时间序列转换为空间序列）；(ii) 基于一维能谱E(n)的积分公式（Λ_x = u0/(2*u_rms^2) * ∫ E(n)/n dn），此方法隐含各向同性假设。横向积分尺度则通过移动一个热线探头测量两点空间互相关函数g®并积分得到。
     • 泰勒微尺度（λ）：采用 Hallbäck 等人（1989）的方法，通过分析计算出的λ值对时间差分步长Δt的依赖性，外推得到精确的λ_x。同时也与基于各向同性湍流衰变的理论公式（λ_x^2 = 20ν(x-x0)/u0）进行了对比。
     • 柯尔莫哥洛夫微尺度（η）：通过耗散率ε计算，η = (ν^3/ε)^{¹⁄₄}。耗散率ε又通过三种独立方法估算：(i) 基于大尺度动能衰变：ε = -u0 * dk/dx；(ii) 基于各向同性假设下的速度梯度方差：ε = 15ν ( (∂u/∂x)^2 )；(iii) 基于泰勒微尺度与各向同性假设：ε = 30ν (u_rms^2 / λ_x^2)。这种交叉验证至关重要。
   • 能谱分析：计算了u， v， w三个方向的一维能谱，并检查了惯性子区（-5/3律）的出现与雷诺数的关系。
   • 顺压梯度影响实验：在网格A和LT5的下游安装了一个非对称收缩段（在顶壁，与来流呈25度角），产生一个伴有微弱横向流的顺压梯度。测量了在此条件下的湍流特性，并与零压梯度（ZPG）情况进行对比。为了公平比较，引入了基于真实对流时间积分（t = ∫ dx/u(x)）的无量纲时间尺度 u0*t/m。

三、 主要研究结果 1. 压力阻力：测量的网格压力阻力系数与雷诺数（基于杆直径， Re_d）的关系符合经典理论。经孔隙率归一化后，数据与文献中的经验曲线趋势一致，但具体系数因网格而异，强调了针对特定网格进行实测标定的必要性。

1. 各向异性与雷诺数效应：所有网格生成的湍流在下游都表现出一定程度的各向异性，通常流向（u）分量最大。一个重要发现是，各向异性程度随着基于网眼宽度的雷诺数（Re_m）的增加而减小，即更高Re_m的流动更接近各向同性，这与 Schedvin 等人（1974）的结论一致。

2. 湍流衰变特性：应用新的客观拟合流程，明确了衰变指数b和系数a对雷诺数的依赖性。b值（绝对值）随着Re_m增加而增加（趋向于约-1.2），a值则减小。这表明不存在“普适”的网格湍流衰变状态，初始条件（由网格几何和雷诺数决定）至关重要。研究还试图确定进入幂律衰变区的临界位置（x/m）_crit，发现其同时依赖于Re_m和固体率σ，但未能找到一个普适的相关式，进一步印证了初始条件的影响。

3. 长度尺度与耗散率：这是本研究贡献最突出的部分。

   • 积分尺度：纵向积分尺度Λ_x的实验数据与基于各向同性假设的经典经验公式（Λ_x/m ∝ (x/m)^{¹⁄₂}）在趋势上相似，但系数差异很大。横向积分尺度Λ_z则能较好地仅用网眼宽度m来标度。研究发现，通过积分自相关函数（方法i）与通过积分能谱（方法ii， 基于各向同性假设）得到的Λ_x值基本重合。这是一个关键结果，意味着即使在非完全各向同性的流动中，基于能谱的积分公式（方法ii）仍然是有效的，这为实验分析提供了便利。
   • 泰勒微尺度：由精确的差分外推法得到的λ_x与经典各向同性衰变理论公式（λ_x^2 = 20ν(x-x0)/u0）在图形上吻合良好，但定量比较发现，对于某些网格，差异可达50%。因此，研究者建议该理论公式仅用于粗略估算。
   • 耗散率与柯尔莫哥洛夫尺度：最具警示性的发现之一是，基于各向同性假设的耗散率公式（ε = 15ν ( (∂u/∂x)^2 )）与基于动能直接衰变的公式（ε = -u0 * dk/dx）计算结果相差可达50%。这强烈提示，在明显各向异性的网格湍流中，盲目应用各向同性关系式计算耗散率会引入巨大误差。基于动能衰变（方法i）和基于泰勒尺度（方法iii）计算出的柯尔莫哥洛夫尺度η吻合良好。一个重要的新结果是，η能够非常漂亮地用m和Re_m进行标度，并得出了经验公式：η/m = 0.80 [ (x/m) / Re_m ]^{³⁄₄}。作者认为此式可推广用于其他实验设置中对η的估计。
   • 雷诺数转换关系：研究给出了宏观雷诺数（Re_Λ = u_rms Λ_x/ν）和微观雷诺数（Re_λ = u_rms λ_x/ν）与参数 [ (x/m) / Re_m ] 的经验关系。发现在 [ (x/m) / Re_m ] ≈ (1-2)×10^{-2} 处，关系曲线的斜率发生明显变化，这被认为是与雷诺数相关的现象。

4. 能谱：三个方向的能谱在用各自积分尺度标度后能很好地重合。对于高雷诺数情况（Re_Λ > 300），能谱中出现了明显的-5/3惯性子区，验证了 Kistler 和 Vrebalovich（1966）提出的最小Re_Λ准则。随着雷诺数增加，惯性子区的范围和-5/3斜率的符合度都得到改善。

5. 顺压梯度的影响：施加一个温和的顺压梯度后，对于网眼宽度较大的网格A，其积分尺度和泰勒微尺度分别减少了约20%和30%；然而，对于网眼宽度较小的网格LT5，长度尺度则没有明显变化。这表明，初始湍流尺度的大小决定了其对流动加速（收缩）的敏感程度。大尺度的湍流结构更能“感受”到平均流的拉伸效应，从而被阻尼。各向异性度的变化也可用流体元拉伸和环量守恒的简单理论定性解释。

四、 结论与意义 本研究成功设计并表征了一组新型低湍流强度生成网格。核心结论是：通过采用固定固体率、变化网眼宽度的网格设计策略，可以在实验上实现湍流强度大致不变而系统改变湍流特征长度尺度的目标。例如，在12 m/s下，网格LT5和网格A在x/m=100处都能产生约1.04%的流向湍流强度，但它们的积分尺度相差5倍，泰勒尺度相差3倍，柯尔莫哥洛夫尺度相差2倍。这为后续研究自由流湍流尺度对边界层转捩等物理过程的影响提供了不可或缺的实验手段。

此外，本研究对网格湍流研究的方法论做出了重要贡献： 1. 开发了确定湍流衰变参数的客观流程。 2. 通过交叉验证，系统评估了基于各向同性假设的经典公式（用于计算积分尺度、泰勒尺度、耗散率）在实际（各向异性）网格湍流中的适用性与误差范围，为实验工作者提供了重要的使用指南和警示。 3. 提出了关于柯尔莫哥洛夫尺度η和雷诺数转换关系的新标度律和经验公式。 4. 揭示了顺压梯度对湍流长度尺度的阻尼效应及其对初始尺度的依赖性。

五、 研究亮点 1. 研究目标新颖且具有高度应用价值：成功实现了“独立调控湍流强度与长度尺度”这一实验难题，直接服务于湍流-转捩相关性研究。 2. 方法论贡献突出：提出的客观衰变参数拟合方法、对多种经典关系式的系统性实验验证与误差量化，均具有很高的参考价值。 3. 数据详实、分析全面：对五个新网格和两个旧网格在宽雷诺数范围内进行了全方位（三维速度、多尺度、能谱、压力梯度影响）的测量与分析，数据集本身非常宝贵。 4. 批判性与建设性并存：不仅指出了广泛使用的各向同性关系式可能存在的巨大误差（批判性），还通过对比给出了更可靠的方法建议，并导出了新的实用经验公式（建设性）。

六、 其他有价值内容 论文最后部分关于顺压梯度影响的初步探索，虽然篇幅不长，但揭示了平均流变形与湍流尺度相互作用的一个有趣现象，为更复杂流动条件下的湍流演化研究提供了线索。同时，论文中关于风洞设施、定制化探头支架和网格安装框架的描述，也体现了高水平实验研究对细节的追求，对实验流体力学工作者有借鉴意义。