2025-11-13T23:07:10.449786

Spatiotemporal statistics of the dissipation rate at the boundary of a turbulent flow using Diffusing-Wave Spectroscopy

Francisco, Lambret, AumaÃ®tre

We use Diffusing Wave Spectroscopy (DWS) to perform the first direct space- and time-resolved measurement of the dissipation rate~$Îµ$ at the boundary of a turbulent flow. We have shown in a previous publication that this technique provides maps of the dissipation rate of Newtonian fluids~\cite{Francisco}. Here, we apply the technique at the boundary of a turbulent flow generated in a square box by an impeller stirring the fluids. Although the measurement is made on a small region near the boundary, we show that the dissipation remains proportional to the injected power and follows the turbulent scaling $Îµ\propto \mathrm{Re}^3$, with Re being the Reynolds number ranging from $1.5 \times 10^4$ to $6 \times 10^5$. With this flow, there is no need for logarithmic corrections to reproduce the dissipation near the flat boundary. In addition, our setup allows us to measure the spatio-temporal fluctuations of the dissipation near the boundary. These fluctuations are quite large (the relative fluctuations are about 50\%) and are well described by a log-normal distribution, as expected for the dissipation rate in the bulk of homogeneous and isotropic turbulence (HIT) but Power Density Spectra (PDS) do not correspond to those expected for HIT \cite{Li07,Graham16,K62}

academic

Spatiotemporal statistics of the dissipation rate at the boundary of a turbulent flow using Diffusing-Wave Spectroscopy

基本信息

论文ID: 2511.03462
标题: Spatiotemporal statistics of the dissipation rate at the boundary of a turbulent flow using Diffusing-Wave Spectroscopy
作者: Enzo Francisco, Julien Lambret, Sébastien Aumâıtre
单位: Service de Physique de l'État Condensé, Université Paris-Saclay, CNRS, CEA, France
分类: physics.flu-dyn (流体动力学)
提交时间: 2025年11月5日
论文链接: https://arxiv.org/abs/2511.03462

摘要

本研究使用扩散波光谱技术(DWS)首次实现了湍流边界处耗散率ε的直接时空分辨测量。研究团队在方形容器中通过叶轮搅拌产生湍流，并在边界附近进行测量。尽管测量区域较小，但结果显示耗散率与注入功率成正比，并遵循湍流标度关系ε∝Re³（雷诺数范围为1.5×10⁴至6×10⁵）。与管道流动不同，该流动在平坦边界附近无需对数修正。此外，边界附近的耗散时空涨落很大（相对涨落约50%），符合对数正态分布，但功率密度谱与均匀各向同性湍流(HIT)的预期不符。

研究背景与动机

研究问题

湍流耗散率是流体力学中的关键物理量，但其直接测量极具挑战性，特别是在边界附近。本研究旨在开发一种能够直接测量湍流边界处耗散率时空分布的实验技术。

问题重要性

工程应用价值：耗散率决定了车辆的能量消耗和管道中的压力损失
地球物理意义：在气候模型中影响能量平衡和传输效率
基础科学意义：耗散结构可能驱动湍流间歇性
边界问题：许多实际应用关注表面附近的耗散和剪切率，如浸没物体的阻力、地球物理流动的侵蚀效率等

现有方法局限性

热线技术：需要假设泰勒冻结湍流假说，多线探针会产生扰动
粒子图像测速(PIV)：
- 精度受限于询问窗口网格尺寸
- 仅适用于中等雷诺数范围
- 双平面立体PIV技术复杂
直接数值模拟(DNS)：计算成本高，长时间模拟困难，可能低估罕见事件

研究动机

牛顿流体中的耗散率定义为： $\epsilon = \sum_{i,j} \frac{\nu}{2}(\partial_j u_i + \partial_i u_j)^2$

这需要测量所有方向上速度场的空间导数，即应变率张量的范数。DWS技术为解决这一难题提供了新途径。

核心贡献

首次实现：首次对湍流边界处耗散率进行直接的时空分辨测量
标度关系验证：证明边界耗散率与注入功率成正比，遵循ε∝Re³标度律，无需对数修正
涨落特性揭示：发现边界耗散涨落达50%，符合对数正态分布
HIT对比分析：发现边界耗散的统计特性部分符合HIT理论（PDF），但功率谱特性不同
技术拓展：将DWS技术从Taylor-Couette流动成功应用到完全湍流流动
高雷诺数测量：实现了Re高达6×10⁵的测量，远超传统方法

方法详解

任务定义

输入：湍流边界处的散射光强度时间序列
输出：边界附近耗散率ε的时空分布图
约束条件：

粒子半径r和光波长λ远小于光传输平均自由程l*
梯度特征长度大于光在流体中的光程
时间分辨率需小于流动演化时间尺度

DWS技术原理

基本方程

相干光照射含小颗粒的运动流体，产生散斑图案。背散射光强度的自相关函数为：

$g_2(\tau) = \frac{\langle I(t+\tau)I(t)\rangle}{\langle I(t+\tau)\rangle\langle I(t)\rangle} = \beta \exp\left(-\gamma\sqrt{6(\tau/\tau_o + \tau^2/\tau_V^2)}\right) + 1$

其中：

β：强度对比度
γ：仅依赖于光学几何和边界条件的参数
τₒ = 1/(Dk²)：布朗运动引起的时间尺度
τᵥ：流体运动引起的时间尺度

耗散率提取

τᵥ与耗散率的关系：

$\tau_V = \sqrt{30}/\left(l^*k\sqrt{\langle\epsilon\rangle/\nu}\right)$

由于τᵥ ∝ λ/l* ≪ 1/√(⟨ε⟩/ν)，可以捕获⟨ε⟩的时间演化。

实验装置

流动系统

容器：60 cm边长的方形水槽
叶轮：半径R=20 cm，包含1.2 cm厚圆盘和4个3.6 cm高叶片
转速范围：Ω = 0.02-0.8 转/秒
雷诺数定义：Re = 2πUL/ν，其中U = RΩ
雷诺数范围：1.5×10⁴ ≤ Re ≤ 6×10⁵
流体：170升去离子水 + 1.7 kg TiO₂颗粒（质量分数1%）

示踪粒子特性

平均水动力学半径：r = 200 nm
Stokes数：Sk ≤ 0.1（保证被动示踪行为）
沉降速度：< 10⁻⁷ μm/s（可忽略）
传输平均自由程：l* = 80 μm

特征长度尺度对比（见表1）

参数	数值范围
波长λ	0.532 μm
粒子半径r	0.2 μm
传输自由程l*	80 μm
容器尺寸L	60 cm
叶轮半径R	20 cm
Taylor微尺度Λₜ	5-0.8 mm
Kolmogorov尺度ηₖ	440-30 μm
粘性底层δᵥ	1mm-25 μm

光学系统

照明系统：2W Nd:YAG激光器（λ=532 nm），通过20×显微镜透镜扩束
两种探测方式：
- PMT+相关器：探测直径12 cm圆盘区域的平均耗散
- 高速相机：分辨5.1×5.1 cm²区域（128×128像素）的涨落
时间分辨率：
- 相关器采样时间：1.28 μs
- 相机帧率：430,000 fps
- 相关函数收敛需要：25,000帧（0.06 s）

测量协议

静止流体测量g₂(τ)，通过Stokes-Einstein公式已知τₒ，推导γ
运行叶轮至设定转速，从方程(3)确定τᵥ
静止测量验证γ值
大多数测量持续约2秒（相机RAM限制）
Re=2×10⁵时进行90秒长时间测量（45个2秒片段）

技术创新点

空间分辨优势：应变率张量估计精度不依赖于相机分辨率，可在相同配置下探索大雷诺数范围
非侵入性：光学测量，不扰动流场
直接测量：无需Taylor假说或复杂的速度场重构
时空分辨：同时获得空间分布和时间演化信息
高雷诺数适用性：突破传统方法的雷诺数限制

实验设置

流体系统参数

工作流体：去离子水（运动粘度ν ≈ 1×10⁻⁶ m²/s）
示踪粒子：TiO₂纳米颗粒
- 密度：ρₛ > ρw（但沉降可忽略）
- 浓度：1%质量分数
- 悬浮稳定性：通过去离子水保证
功率测量：扭矩计测量转矩C，计算注入功率P = C·Ω

光学参数验证

l*原位测量：空间分辨反射技术
- 使用另一束250 mW Nd:YAG激光
- CCD相机测量扩散光斑径向强度衰减
- 遵循∝(r/l*)⁻³规律
- 48小时实验期间l*保持稳定（连续搅拌条件下）

测量区域

PMT探测区域：直径12 cm圆盘，深度约5l* ≈ 400 μm
相机探测区域：5.1×5.1 cm²正方形，深度约5l* ≈ 400 μm
像素尺寸：约400 μm（大于l*，满足要求）
测量位置：叶轮高度处的容器壁面

实验条件控制

搅拌策略：高低转速交替，保持颗粒均匀混合
测量时长：
- 常规测量：约2秒
- 长时间测量：90秒（用于PDF统计）
宏观时间尺度：1/(4Ω) = 0.3-12.5秒

实验结果

主要结果

1. 耗散率标度关系（图2）

核心发现：

边界耗散率⟨ε⟩与注入功率P成正比
遵循湍流零次定律：⟨ε⟩ ∝ P ∝ (ν³/L⁴)·Re³
无需对数修正：与管道/通道流动不同

定量结果：

测量的边界耗散率约为平均注入功率的2.4倍
这种过量耗散符合叶轮附近边界过量耗散的预期

物理解释：

单叶轮在方形容器中产生的流动高度不均匀
耗散结构在空间分布不均匀
叶轮附近边界处存在过量耗散

2. 空间涨落特性（图3）

观测结果：

相对涨落：σε/⟨ε⟩ ≈ 50%（独立于Re）
标准差标度：σε ∝ (ν³/L⁴)·Re³
耗散结构尺度显著大于像素尺寸

结构特征：

耗散结构被大尺度流动平流输运
结构被拉伸变形
主流方向：环形流动（y方向）+ 上升极向分量

时间演化：

PMT测量（12 cm直径平均）：涨落约30%
未观察到与叶片通过频率相关的特征频率

3. 概率密度函数（图4）

统计分析：

90秒长时间测量（Re=2×10⁵）
变量：(log(ε) - ⟨log(ε)⟩)/σlog(ε)

关键发现：

PDF符合对数正态分布
与HIT理论预期一致
但测量位置在湍流边界层，非HIT区域

理论对比：

HIT中对数正态分布基于级联过程
边界层中存在复杂的层级结构
实验结果表明对数正态性在边界层仍然成立

4. 功率密度谱（图5）

实验谱特征：

小尺度涨落随Re增加而增强
像素分辨率足以捕获所有相关结构
在kηₖ = 1/ηₖ之前达到噪声水平（假设在惯性区）

与HIT理论对比：

精细Kolmogorov理论预测： $\langle\epsilon(x)\epsilon(x+r)\rangle \propto \langle\epsilon^2\rangle(L/r)^\mu$ 对应功率谱：k⁻⁽¹⁻μ⁾，μ ≈ 0.2-0.3
实验观测：谱衰减比k⁻⁰·⁸更陡

各向异性分析：

难以识别优先方向
使用径向波数k = √(k²ₓ + k²y)

DNS验证（Johns Hopkins数据库）

模拟参数

流动类型：湍流通道流
雷诺数：Re ≈ 4×10⁴
网格：2048×1536×512节点
计算域：8π × 3π × 2
边界条件：x、y方向周期，z方向无滑移
观测窗口：20×40节点（对应实验5.1×5.1 cm²）

DNS结果对比

耗散结构（图3c、3d）：

实验与DNS显示相似结构
DNS中各向异性更明显（平均流动方向明确）

PDF对比（图4）：

DNS涨落也符合对数正态分布
与实验结果一致

功率谱对比（图5插图）：

展向和流向功率谱
两个方向衰减均比HIT预测k⁻⁰·⁸更陡
与实验观测一致

重要结论： DNS结果证实实验发现不是测量误差，而是边界湍流的真实特性。

实验发现总结

标度律验证：边界耗散率遵循ε∝Re³，无需对数修正
局部过量耗散：叶轮高度处耗散率为平均值的2.4倍
大幅涨落：相对涨落50%，独立于Re
对数正态性：PDF符合HIT理论预期
谱特性差异：PDS衰减比HIT预测更快
DNS一致性：通道流DNS证实边界特性

方法	优点	局限性
热线技术12	高频响应	需Taylor假说；多线探针扰动流场
点测量PIV13	直接速度测量	单点测量；中等Re
双平面立体PIV14	平面分辨	技术复杂；网格限制精度
DNS11,15	完整信息	计算成本；低估罕见事件
DWS（本文）	直接测量；高Re；时空分辨	边界测量；需透明流体

结论与讨论

主要结论

技术突破：DWS首次实现湍流边界耗散率的直接时空分辨测量，在Re高达6×10⁵时仍有良好空间和时间分辨率
标度关系：边界耗散率与注入功率成正比，遵循纯幂律标度ε∝Re³，无需管道流中的对数修正
局部特性：叶轮高度处存在易于估计的过量耗散（2.4倍）
涨落特性：
- 耗散结构涨落强烈（相对涨落50%）
- 符合对数正态分布（与HIT一致）
- 功率谱衰减比HIT预测更快
普适性：通道流DNS证实测量特性代表湍流边界的普遍行为

物理意义

边界耗散机制

叶轮驱动流动的边界耗散不同于经典壁湍流
体耗散结构主导边界层耗散
粘性底层的统计权重随Re变化（93.6%→2.9%）

统计特性的双重性

对数正态性保持：级联过程的影响延伸到边界
谱特性改变：边界层复杂结构改变空间相关性

局限性

测量位置限制：
- 仅测量边界附近薄层（深度~400 μm）
- 无法直接测量体耗散
光学条件要求：
- 需要透明或半透明流体
- 需要合适的示踪粒子（r, λ ≪ l*）
时间分辨率约束：
- 相关函数收敛需要0.06秒
- 要求流动演化时间尺度足够长
- 限制了可研究的流动类型
空间分辨率：
- 像素尺寸~400 μm，大于最小尺度
- l*=80 μm在Re范围内穿越Kolmogorov尺度
- 高Re时可能失效（l* ≤ ηₖ）
流动类型：
- 本文流动高度不均匀
- 与HIT和经典壁湍流都有差异
- 普适性需进一步验证
统计收敛：
- 常规测量仅2秒
- PDF分析仅基于一次90秒测量
- 统计收敛性有待提高