[论文解读] 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）方程。该研究提出了一种专用的三阶龙格-库塔时间积分格式，并设计了一种新颖的随机应力离散化方法，在3D模拟中实现了精确的涨落-耗散平衡与最小的长波长时间误差，即使在接近稳定性极限时亦保持稳定。

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方程中的精度，超越平衡方差与关联函数的测量范围。
• 解决在随机流体动力学模拟中，长波长精度与短波长稳定性之间的平衡挑战。
• 开发一种时间积分格式，以最小化1D线性化LLNS方程离散结构因子中的长波长时间误差。
• 设计一种新颖的3D随机应力张量离散化方法，确保在小时间步长下保持离散涨落-耗散平衡。
• 构建一个完全一致的3D格式，即使在时间步长接近稳定性极限时仍能保持精度。

提出的方法

• 使用离散平衡结构因子作为诊断工具，评估在不同波数和频率下的空间-时间离散化精度。
• 推导一种专用的三阶龙格-库塔格式，专为最小化1D线性化LLNS方程在长波长下的时间积分误差而优化。
• 提出一种在大于一维情况下离散随机应力张量的新方法，确保涨落-耗散平衡的一致性。
• 通过傅里叶空间分析解耦波矢量模态，并在不同格式下评估离散结构因子。
• 将优化的时间积分格式与新的应力离散化方法结合，构建一种具有更高精度与稳定性的3D有限体积格式。
• 通过验证表明，该格式在小时间步长下实现了精确的离散涨落-耗散平衡，并在接近稳定性极限时仍保持鲁棒的精度。

实验结果

研究问题

• RQ1 如何系统量化有限体积格式在LLNS方程中不同波数和频率下的精度？
• RQ2 何种时间积分格式能最小化1D线性化LLNS方程离散结构因子中的长波长时间误差？
• RQ3 如何在3D中对随机应力张量进行离散化，以保持离散涨落-耗散关系？
• RQ4 在保持精确涨落-耗散平衡与结构因子精度的前提下，时间步长可被增大到何种程度？
• RQ5 能否构建一种3D有限体积格式，使其在长波长下保持精度，同时在时间步长极限附近保持稳定？

主要发现

• 所提出的三阶龙格-库塔格式能最小化1D线性化LLNS方程在长波长下的离散结构因子中的时间积分误差。
• 新的3D随机应力离散化方法确保该格式在小时间步长下满足离散涨落-耗散关系。
• 组合格式即使在时间步长接近稳定性极限时，仍能保持准确的结构因子与涨落-耗散平衡。
• Crank-Nicolson半隐式方法能实现与时间步长无关的完美静态结构因子（S_k = 1），但存在空间离散化误差。
• 对于对流-扩散系统，采用预测-校正处理的半隐式格式在r = 0时保持S_k = 1，并在|r| > 0时与PC-2RNG的主导误差阶数一致。
• 通过结构因子进行的分析为区分长波长精度与短波长稳定性问题提供了更优的诊断手段。