Tuned Finite-Difference Diffusion Operators

TL;DR

This paper introduces tuned finite-difference diffusion operators that optimize stability and physical accuracy by acting at the diffusion scale rather than the grid scale, thereby improving simulation efficiency in fluid and magnetohydrodynamics.

Contribution

The authors develop finite-difference diffusion operators that minimize timestep restrictions while targeting the diffusion scale, enhancing stability and physical fidelity in simulations.

Findings

Operators peak at the diffusion scale, not the grid scale

They enable larger timesteps in high Schmidt and Prandtl number simulations

Operators behave as standard at larger scales, preserving accuracy

Abstract

Finite-difference simulations of fluid dynamics and magnetohydrodynamics generally require an explicit diffusion operator, either to maintain stability by attenuating grid-scale structure, or to implement physical diffusivities such as viscosity or resistivity. If the goal is stability only, the diffusion must act at the grid scale, but should affect structure at larger scales as little as possible. For physical diffusivities the diffusion scale depends on the problem, and diffusion may act at larger scales as well. Diffusivity undesirably limits the computational timestep in both cases. We construct tuned finite-difference diffusion operators that minimally limit the timestep while acting as desired near the diffusion scale. Such operators reach peak values at the diffusion scale rather than at the grid scale, but behave as standard operators at larger scales. We focus on the specific applications of hyperdiffusivity for numerical stabilization, and high Schmidt and high Prandtl number simulations where the diffusion scale greatly exceeds the grid scale.