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

渦拡散係数モーメント定量化における統計収束改善技術

基本情報

論文ID: 2503.06418
タイトル: Techniques for improved statistical convergence in quantification of eddy diffusivity moments
著者: Dana L. O.-L. Lavacot (スタンフォード大学)、Jessie Liu (スタンフォード大学)、Brandon E. Morgan (ローレンス・リバモア国立研究所)、Ali Mani (スタンフォード大学)
分類: physics.flu-dyn
発表日: 2025年10月15日 (arXiv プレプリント)
論文リンク: https://arxiv.org/abs/2503.06418v2

要旨

本論文は、巨視的強制法(MFM)における渦拡散係数モーメント計算の統計収束性問題に対して、改善技術を提案している。研究により、すべての強制MFM シミュレーションに対して同一の直接数値シミュレーション(DNS)を使用することで、渦拡散係数モーメントの統計収束性が著しく改善されることが示された。著者らは分解法を組み合わせた実装方案を開発し、MFM強制を半解析的に処理し、境界条件の一貫性を確保することができる。二次元Rayleigh-Taylor不安定性のケーススタディでは、本手法により収束に必要なシミュレーション回数をO(1000)からO(100)に削減した。

研究背景と動機

問題背景

乱流モデリングにおいて、Reynolds平均Navier-Stokes (RANS)方程式は、スカラー通量やReynolds応力などの未閉包項を処理するための閉包モデルを必要とする。従来のBoussinesq近似は局所的かつ等方的な渦拡散係数を仮定しているが、複雑な流れに対しては往々にして無効である。

中心的課題

非局所性と異方性効果: 実際の渦拡散係数は通常、非局所的かつ異方的であり、より正確な記述が必要である
統計収束の遅さ: MFMを使用して高次の渦拡散係数モーメントを計算する際、統計収束性が不良である
計算コストの高さ: 従来の手法では統計収束を達成するために多数のDNSシミュレーションが必要である

研究動機

既存の巨視的強制法(MFM)とGreen関数法は非局所的かつ異方的な渦拡散係数を計算できるが、高次モーメント計算時に統計収束が遅い問題が存在する。本論文は統計収束を加速し、計算コストを削減する改善技術の開発を目指している。

中核的貢献

統計誤差増幅問題の識別と解決: 独立したドナーシミュレーションの使用が高次モーメントの統計誤差増幅を招くことを発見し、単一ドナー使用による解決策を提案
分解MFM法の開発: 一般的な非定常かつカオス的流れ、および運動量輸送に対する分解法の拡張
計算効率の大幅な向上: 2D Rayleigh-Taylor不安定性において約10倍の計算加速を実現
完全な理論および実装フレームワークの提供: スカラーおよび運動量輸送に対する完全な数学的フレームワークと実装詳細を提供

方法論の詳細

タスク定義

本論文は、非局所的かつ異方的な渦拡散係数の時空間モーメントを効率的に計算する方法を研究している。これらのモーメントは渦拡散係数カーネルの有限ながら重要な情報を含み、非局所性および異方性効果の定量化とモデリングに十分である。

渦拡散係数モーメントの数学的記述

渦拡散係数はKramers-Moyal展開によってモーメント形式で表現できる:

$-\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}$

ここで各次数のモーメントは以下のように定義される:

$D^{00}_{ij}(x,t) = \int\int D_{ij}(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'$ (一次空間モーメント)
$D^{01}_{ij}(x,t) = \int\int (t'-t)D_{ij}(x,x',t,t')dx'dt'$ (一次時間モーメント)

巨視的強制法(MFM)

MFMはスカラー輸送方程式に巨視的強制を追加することで渦拡散係数モーメントを計算する:

$\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'(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強制の半解析的処理が可能となり、周期境界条件の適合性問題が解決される。

技術的革新点

単一ドナー戦略: すべてのレシーバー方程式に対して同一のDNSシミュレーションを流場供給源として使用することで、独立したドナーシミュレーション間の差異に起因する誤差増幅を回避
誤差伝播分析: 理論分析により、独立したドナー使用時に高次モーメントの統計誤差が空間座標の冪乗に従って増幅されることが示される: $D^1 = F^1 - x_1 F^0 + x_1(\varepsilon^{01} - \varepsilon^{00}) + \varepsilon^{11}$
分解法の拡張: 元々定常層流に適用されていた分解法を、一般的な非定常カオス流れおよび運動量輸送に拡張

実験設定

ケーススタディ: 二次元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