[論文レビュー] On the Accuracy of Finite-Volume Schemes for Fluctuating Hydrodynamics

A. Donev, Eric Vanden‐Eijnden+2 more|ArXiv.org|Jun 12, 2009

Fluid Dynamics and Turbulent Flows参考文献 53被引用数 4

ひとこと要約

本稿では、離散的平衡構造因子の波数および周波数依存解析を導入することで、ランダウ＝リフシッツ＝ナビエ＝ストークス（LLNS）方程式に対する高精度な有限体積スキームを開発する。特許的な3次ルンゲ＝クッタ時間積分スキームと、新たな確率的応力の離散化手法を提示し、3次元シミュレーションにおいて、安定性限界に近い状況でも、正確なフラクチュエーション・ディスシペーションバランスと最小限の長波長時間誤差を達成する。

ABSTRACT

This paper describes the development and analysis of finite-volume methods for the Landau-Lifshitz Navier-Stokes (LLNS) equations and related stochastic partial differential equations in fluid dynamics. The LLNS equations incorporate thermal fluctuations into macroscopic hydrodynamics by the addition of white-noise fluxes whose magnitudes are set by a fluctuation-dissipation relation. Originally derived for equilibrium fluctuations, the LLNS equations have also been shown to be accurate for non-equilibrium systems. Previous studies of numerical methods for the LLNS equations focused primarily on measuring variances and correlations computed at equilibrium and for selected non-equilibrium flows. In this paper, we introduce a more systematic approach based on studying discrete equilibrium structure factors for a broad class of explicit linear finite-volume schemes. This new approach provides a better characterization of the accuracy of a spatio-temporal discretization as a function of wavenumber and frequency, allowing us to distinguish between behavior at long wavelengths, where accuracy is a prime concern, and short wavelengths, where stability concerns are of greater importance. We use this analysis to develop a specialized third-order Runge Kutta scheme that minimizes the temporal integration error in the discrete structure factor at long wavelengths for the one-dimensional linearized LLNS equations. Together with a novel method for discretizing the stochastic stress tensor in dimension larger than one, our improved temporal integrator yields a scheme for the three-dimensional equations that satisfies a discrete fluctuation-dissipation balance for small time steps and is also sufficiently accurate even for time steps close to the stability limit.

研究の動機と目的

• 平衡分散および相関測定を超えた、LLNS方程式の有限体積スキームの精度を体系的に評価すること。
• 確率的流体力学シミュレーションにおいて、長波長での精度と短波長での安定性の両立という課題に対処すること。
• 1次元線形化LLNS方程式の離散的構造因子における長波長時間誤差を最小化する時間積分スキームを開発すること。
• 時間刻みが小さい場合に離散的フラクチュエーション・ディスシペーションバランスを保つ、3次元での新しいステージストレステンソルの離散化法を設計すること。
• 時間刻みが安定性限界に近づく状況でも精度を維持できる、完全に整合性のある3次元スキームを構築すること。

提案手法

• 波数および周波数にわたる空間的・時間的離散化の精度を評価するため、離散的平衡構造因子を診断指標として用いる。
• 1次元線形化LLNS方程式において、長波長での構造因子における時間積分誤差を最小化するように最適化された特許的な3次ルンゲ＝クッタスキームを導出する。
• 1次元を超える次元における確率的応力テンソルの新しい離散化法を導入し、一貫性のあるフラクチュエーション・ディスシペーションバランスを確保する。
• 波数ベクトルモードを分離するためにフーリエ空間解析を適用し、さまざまなスキーム下での離散的構造因子を評価する。
• 最適化された時間積分スキームと新しい応力離散化法を組み合わせ、3次元有限体積スキームを構築し、精度と安定性を向上させる。
• 時間刻みが小さい場合に正確な離散的フラクチュエーション・ディスシペーションバランスを示すことで、スキームの妥当性を検証する。また、安定性限界に近い状況でも高い精度を示すことを確認する。

実験結果

リサーチクエスチョン

• RQ1 LLNS方程式の有限体積スキームの精度を、さまざまな波数および周波数にわたって体系的にどのように定量的評価できるか？
• RQ2 1次元線形化LLNS方程式の離散的構造因子において、長波長時間誤差を最小化する時間積分スキームは何か？
• RQ3 3次元における確率的応力テンソルをどのように離散化すれば、離散的フラクチュエーション・ディスシペーション関係を保てるか？
• RQ4 フラクチュエーション・ディスシペーションバランスおよび構造因子の精度を維持しつつ、時間刻みをどの程度まで大きくできるか？
• RQ5 長波長での精度と時間刻みの安定性限界付近での安定性を両立できる3次元有限体積スキームを構築できるか？

主な発見

• 提案された3次ルンゲ＝クッタスキームは、1次元線形化LLNS方程式の離散的構造因子における長波長時間誤差を最小化する。
• 新規の3次元確率的応力離散化法により、時間刻みが小さい場合にスキームが離散的フラクチュエーション・ディスシペーション関係を満たすことが保証される。
• 組み合わせたスキームは、時間刻みが安定性限界に近づく状況でも、正確な構造因子とフラクチュエーション・ディスシペーションバランスを維持する。
• クランク＝ニコルソン半陽的スキームは、時間刻みに依存せずに正確な静的構造因子（S_k = 1）を達成するが、空間離散化誤差を有する。
• 対流拡散系において、予測子・修正子処理を施した半陽的スキームは、r = 0 のとき S_k = 1 を維持し、|r| > 0 のときには PC-2RNG の一次誤差と一致する。
• 構造因子を用いた解析により、長波長精度と短波長安定性の問題を明確に区別する優れた診断法が得られる。