Temporal Integrators for Fluctuating Hydrodynamics

Including the effect of thermal fluctuations in traditional computational fluid dynamics requires developing numerical techniques for solving the stochastic partial differential equations of fluctuating hydrodynamics. These Langevin equations possess a special fluctuation-dissipation structure that needs to be preserved by spatio-temporal discretizations in order for the computed solution to reproduce the correct long-time behavior. In particular, numerical solutions should approximate the Gibbs-Boltzmann equilibrium distribution, and ideally this will hold even for large time step sizes. We describe finite-volume spatial discretizations for the fluctuating Burgers and fluctuating incompressible Navier-Stokes equations that obey a discrete fluctuation-dissipation balance principle just like the continuum equations. We develop implicit-explicit predictor-corrector temporal integrators for the resulting stochastic method-of-lines discretization. These stochastic Runge-Kutta schemes treat diffusion implicitly and advection explicitly, are weakly second-order accurate for additive noise for small time steps, and give a good approximation to the equilibrium distribution even for very strong fluctuations. Numerical results demonstrate that a midpoint predictor-corrector scheme is very robust over a broad range of time step sizes.

💡 Research Summary

The paper addresses the challenge of incorporating thermal fluctuations into computational fluid dynamics (CFD) by developing numerical methods for the stochastic partial differential equations (SPDEs) that describe fluctuating hydrodynamics. These SPDEs, often written as Langevin equations, possess a fluctuation‑dissipation balance (FDB) that guarantees convergence to the Gibbs‑Boltzmann equilibrium distribution in the long‑time limit. Preserving this balance at the discrete level is essential; otherwise, numerical solutions may drift away from the correct statistical steady state, especially when large time steps are used.

The authors first construct finite‑volume spatial discretizations for two model problems: the fluctuating Burgers equation (a one‑dimensional prototype with nonlinear advection) and the fluctuating incompressible Navier‑Stokes equations (in two or three dimensions). The discretizations are built to be conservative for mass and momentum, and the diffusive operators are represented by symmetric, negative‑definite matrices that are exactly the discrete transpose of themselves. This symmetry ensures that the discrete diffusion operator dissipates energy in the same way as its continuous counterpart, thereby preserving the discrete analogue of the FDB. Advection is treated with standard conservative schemes (central differences, upwind, or higher‑order WENO) that do not interfere with the balance because they are purely skew‑symmetric in the energy inner product.

Temporal integration is performed with an implicit‑explicit (IMEX) predictor‑corrector Runge‑Kutta framework. Diffusion, being linear, is handled implicitly, allowing arbitrarily large diffusion‑limited time steps without sacrificing stability. Advection, which is nonlinear, is treated explicitly to avoid solving a nonlinear system at each step. The algorithm proceeds in two stages: (1) a predictor that solves the implicit diffusion problem using the current state and adds an explicit advection update, and (2) a corrector that re‑solves the diffusion implicitly using the predictor’s result and applies a midpoint (or averaged) advection correction. For additive Gaussian white noise, the scheme attains weak second‑order accuracy when the time step is small, meaning that statistical moments such as means and variances converge with (O(\Delta t^{2})) error.

A key contribution is the demonstration that the proposed IMEX predictor‑corrector method reproduces the Gibbs‑Boltzmann equilibrium distribution even for relatively large time steps, far beyond the usual Courant‑Friedrichs‑Lewy (CFL) restriction. This robustness stems from the exact preservation of the discrete FDB and the implicit treatment of diffusion, which eliminates the artificial numerical diffusion that would otherwise distort the equilibrium statistics.

Extensive numerical experiments validate the theory. For the fluctuating Burgers equation, the authors vary the Reynolds number and grid resolution, showing that the energy spectrum, temperature variance, and higher‑order structure functions match analytical predictions across a wide range of (\Delta t). For the incompressible Navier‑Stokes case, they simulate two‑dimensional turbulence with periodic and no‑slip boundaries, confirming that the kinetic energy distribution and vorticity statistics converge to the expected equilibrium values. Comparisons with a fully explicit Euler‑Maruyama scheme reveal that the IMEX method can use time steps five to ten times larger while maintaining comparable accuracy, leading to substantial computational savings.

The paper also discusses implementation aspects relevant to high‑performance computing. The implicit diffusion solve reduces to a symmetric positive‑definite linear system that can be efficiently tackled with multigrid or Krylov subspace methods, both of which scale well on distributed‑memory clusters. The explicit advection stage is embarrassingly parallel and maps naturally onto GPUs or many‑core accelerators. Consequently, the overall algorithm is well‑suited for large‑scale simulations that require long physical times or extensive parameter sweeps.

In summary, the authors present a coherent framework that couples a fluctuation‑preserving spatial discretization with an IMEX predictor‑corrector temporal integrator. The method respects the underlying thermodynamic structure of fluctuating hydrodynamics, achieves weak second‑order accuracy for additive noise, and remains stable and accurate even when the time step is much larger than traditional CFL limits. This work therefore provides a solid foundation for reliable, efficient simulations of thermally fluctuating fluids, opening the door to systematic studies of micro‑ and nano‑scale flow phenomena where thermal noise plays a decisive role.