2025-11-19T00:46:13.540157

Techniques for improved statistical convergence in quantification of eddy diffusivity moments

Lavacot, Liu, Morgan et al.
While recent approaches, such as the macroscopic forcing method (MFM) or Green's function-based approaches, can be used to compute Reynolds-averaged Navier--Stokes closure operators using forced direct numerical simulations, MFM can also be used to directly compute moments of the effective nonlocal and anisotropic eddy diffusivities. The low-order spatial and temporal moments contain limited information about the eddy diffusivity but are often sufficient for quantification and modeling of nonlocal and anisotropic effects. However, when using MFM to compute eddy diffusivity moments, the statistical convergence can be slow for higher-order moments. In this work, we demonstrate that using the same direct numerical simulation (DNS) for all forced MFM simulations improves statistical convergence of the eddy diffusivity moments. We present its implementation in conjunction with a decomposition method that handles the MFM forcing semi-analytically and allows for consistent boundary condition treatment, which we develop for both scalar and momentum transport. We demonstrate that for a two-dimensional Rayleigh--Taylor instability case study, using the same DNS for all forced MFM simulations results in convergence with O(100) simulations rather than O(1000) simulations. We then demonstrate the impacts of improved convergence on the quantification of the eddy diffusivity.
academic

Techniques for improved statistical convergence in quantification of eddy diffusivity moments

基本信息

  • 论文ID: 2503.06418
  • 标题: Techniques for improved statistical convergence in quantification of eddy diffusivity moments
  • 作者: Dana L. O.-L. Lavacot (Stanford University), Jessie Liu (Stanford University), Brandon E. Morgan (Lawrence Livermore National Laboratory), Ali Mani (Stanford University)
  • 分类: physics.flu-dyn
  • 发表时间: October 15, 2025 (arXiv预印本)
  • 论文链接: https://arxiv.org/abs/2503.06418v2

摘要

本文针对宏观强迫方法(MFM)在计算涡扩散系数矩的统计收敛性问题，提出了改进技术。研究表明，对所有强迫MFM模拟使用相同的直接数值模拟(DNS)可以显著改善涡扩散系数矩的统计收敛性。作者开发了结合分解方法的实现方案，能够半解析处理MFM强迫并确保边界条件的一致性处理。在二维Rayleigh-Taylor不稳定性案例研究中，该方法将收敛所需的模拟次数从O(1000)降低到O(100)。

研究背景与动机

问题背景

在湍流建模中，Reynolds平均Navier-Stokes (RANS)方程需要闭合模型来处理未闭合项，如标量通量或Reynolds应力。传统的Boussinesq近似假设局部和各向同性的涡扩散系数，但对于复杂流动往往无效。

核心问题

  1. 非局部和各向异性效应：真实的涡扩散系数通常是非局部和各向异性的，需要更精确的描述
  2. 统计收敛缓慢：使用MFM计算高阶涡扩散系数矩时，统计收敛性较差
  3. 计算成本高昂：传统方法需要大量的DNS模拟才能达到统计收敛

研究动机

现有的宏观强迫方法(MFM)和Green函数方法可以计算非局部和各向异性的涡扩散系数，但在计算高阶矩时存在统计收敛缓慢的问题。本文旨在开发改进技术来加速统计收敛，降低计算成本。

核心贡献

  1. 识别并解决统计误差放大问题：发现使用独立donor模拟会导致高阶矩的统计误差放大，提出使用单一donor的解决方案
  2. 开发分解MFM方法：扩展了分解方法用于一般非定常和混沌流动，以及动量输运
  3. 显著提高计算效率：在2D Rayleigh-Taylor不稳定性中实现了约10倍的计算加速
  4. 提供完整的理论和实现框架：为标量和动量输运提供了完整的数学框架和实现细节

方法详解

任务定义

本文研究如何高效计算非局部和各向异性涡扩散系数的时空矩，这些矩包含了涡扩散系数核的有限但重要的信息，足以量化和建模非局部和各向异性效应。

涡扩散系数矩的数学描述

涡扩散系数可以通过Kramers-Moyal展开表示为矩的形式：

−⟨ui′c′⟩(x,t)=[Dij00(x,t)+Dijk10(x,t)∂∂xk+⋯+Dij01(x,t)∂∂t+⋯ ]∂⟨c⟩∂xj-\langle u'_i c' \rangle(x,t) = \left[D^{00}_{ij}(x,t) + D^{10}_{ijk}(x,t)\frac{\partial}{\partial x_k} + \cdots + D^{01}_{ij}(x,t)\frac{\partial}{\partial t} + \cdots\right]\frac{\partial\langle c \rangle}{\partial x_j}

其中各阶矩定义为：

  • Dij00(x,t)=∫∫Dij(x,x′,t,t′)dx′dt′D^{00}_{ij}(x,t) = \int\int D_{ij}(x,x',t,t')dx'dt' (零阶时空矩)
  • Dijk10(x,t)=∫∫(xk′−xk)Dij(x,x′,t,t′)dx′dt′D^{10}_{ijk}(x,t) = \int\int (x'_k-x_k)D_{ij}(x,x',t,t')dx'dt' (一阶空间矩)
  • Dij01(x,t)=∫∫(t′−t)Dij(x,x′,t,t′)dx′dt′D^{01}_{ij}(x,t) = \int\int (t'-t)D_{ij}(x,x',t,t')dx'dt' (一阶时间矩)

宏观强迫方法(MFM)

MFM通过在标量输运方程中添加宏观强迫来计算涡扩散系数矩：

∂c∂t+∂∂xi(uic)=DM∂2c∂xi∂xi+s\frac{\partial c}{\partial t} + \frac{\partial}{\partial x_i}(u_i c) = D_M \frac{\partial^2 c}{\partial x_i \partial x_i} + s

通过指定不同的平均标量梯度（如多项式形式），可以直接获得相应的涡扩散系数矩。

分解MFM方法

分解方法将标量脉动c′c'展开为：

c′(x,t)=[cj00(x,t)+cjk10(x,t)∂∂xk+⋯+cj01(x,t)∂∂t+⋯ ]∂⟨c⟩∂xjc'(x,t) = \left[c^{00}_j(x,t) + c^{10}_{jk}(x,t)\frac{\partial}{\partial x_k} + \cdots + c^{01}_j(x,t)\frac{\partial}{\partial t} + \cdots\right]\frac{\partial\langle c \rangle}{\partial x_j}

这允许半解析处理MFM强迫，解决了周期边界条件的兼容性问题。

技术创新点

  1. 单一donor策略：通过使用相同的DNS模拟为所有receiver方程提供流场，避免了独立donor模拟间的差异导致的误差放大
  2. 误差传播分析：理论分析表明，使用独立donor时，高阶矩的统计误差会按空间坐标的幂次放大： D1=F1−x1F0+x1(ε01−ε00)+ε11D^1 = F^1 - x_1 F^0 + x_1(\varepsilon^{01} - \varepsilon^{00}) + \varepsilon^{11}
  3. 分解方法的扩展：将原本用于定常层流的分解方法扩展到一般非定常混沌流动和动量输运

实验设置

案例研究：二维Rayleigh-Taylor不稳定性

  • 计算域：2049×2049网格的正方形域
  • 边界条件：x₁方向周期性，x₂方向无滑移和无穿透
  • 物理参数：
    • Atwood数：A = 0.05
    • 最大Mach数：Ma_max = 0.05
    • Grashof数：Gr = 1
    • Schmidt数：Sc = 1
    • Reynolds数：Re_T = 54, Re_L = 8,000

评价指标

  • 统计收敛性：通过不同实现次数下的均方根误差(RMSE)评估
  • 计算效率：比较达到相同收敛精度所需的模拟次数
  • 建模性能：通过匹配矩逆(MMI)方法构建的闭合算子性能评估

对比方法

  • 标准MFM（独立donor）：每个矩使用独立的donor-receiver对
  • 分解MFM（单一donor）：所有receiver方程使用同一个donor模拟

实验结果

主要结果

统计收敛性对比

使用单一donor的分解MFM相比独立donor的标准MFM在统计收敛性上有显著改善：

  1. 零阶矩D⁰⁰：两种方法表现相似，因为计算方式数学上等价
  2. 一阶时间矩D⁰¹：分解MFM在100个实现下即可达到良好收敛，而标准MFM需要1000个实现
  3. 一阶空间矩D¹⁰：分解MFM在100个实现下的误差明显低于标准MFM

计算效率提升

  • 加速比：约10倍的计算加速（从O(1000)降至O(100)个模拟）
  • 误差缩放：理论分析和数值结果都证实了误差按时间和空间坐标放大的现象

消融实验：湍流槽道流验证

在Re_τ = 180的湍流槽道流中验证了方法的有效性：

  • 标准MFM（单一donor）与分解MFM结果几乎完全一致
  • 证明了改进主要来自donor的统一而非分解方法本身

建模性能评估

使用改进收敛的涡扩散系数矩构建的MMI闭合算子：

  • 分解MFM：在200个实现下构建的模型与DNS结果吻合良好
  • 标准MFM：相同实现数下构建的模型出现显著偏差，可能导致错误的建模结论

相关工作

涡扩散系数计算方法

  1. Green函数方法：Hamba等发展的精确但计算昂贵的方法
  2. 快速MFM：Bryngelson等开发的利用隐藏稀疏性的加速方法
  3. 伴随MFM：Liu等提出的基于伴随方法的目标计算技术

非局部建模理论

  • Kraichnan的精确非局部表达式
  • 矩-梯度方法在扩散建模中的应用
  • Kramers-Moyal展开在湍流建模中的应用

结论与讨论

主要结论

  1. 统计误差源识别：证明了独立donor模拟间的微小差异会在高阶矩计算中放大为O(1)的误差
  2. 有效解决方案：使用单一donor可以消除误差放大，显著改善统计收敛性
  3. 方法通用性：分解MFM适用于标量和动量输运，能处理周期边界条件问题

局限性

  1. 流动类型限制：主要在2D Rayleigh-Taylor不稳定性中验证，需要更多流动类型的测试
  2. 数值实现依赖：改进程度可能依赖于具体的数值求解器和并行实现
  3. 理论分析简化：误差分析基于简化假设，实际情况可能更复杂

未来方向

  1. 扩展应用：在更多复杂流动中验证方法有效性
  2. 理论完善：发展更完整的误差传播理论
  3. 算法优化：结合其他加速技术进一步提高计算效率

深度评价

优点

  1. 问题识别准确：准确识别了MFM中的统计误差放大问题
  2. 解决方案有效：提出的单一donor策略简单而有效
  3. 理论分析深入：提供了误差传播的数学分析
  4. 实验验证充分：通过多个案例验证了方法的有效性
  5. 实用价值高：显著降低了计算成本，提高了方法的实用性

不足

  1. 案例有限：主要在2D RT不稳定性中验证，缺乏3D和其他类型流动的验证
  2. 理论局限：误差分析基于一维简化情况，实际多维情况更复杂
  3. 依赖性问题：改进效果可能依赖于特定的数值实现和并行策略

影响力

  1. 学术贡献：为MFM方法的改进提供了重要技术路径
  2. 实用价值：显著提高了涡扩散系数矩计算的效率
  3. 方法推广：分解MFM方法具有良好的推广前景

适用场景

  1. 湍流建模：需要精确涡扩散系数的RANS建模
  2. 混合问题：Rayleigh-Taylor、Richtmyer-Meshkov等不稳定性研究
  3. 非局部效应研究：需要量化非局部和各向异性效应的流动问题

参考文献

本文引用了37篇重要文献，涵盖了湍流理论、数值方法、和MFM相关研究，为研究提供了坚实的理论基础。

总体评价：这是一篇高质量的计算流体力学论文，针对实际计算问题提出了有效的解决方案。方法简单实用，理论分析深入，实验验证充分，具有重要的学术价值和实用意义。