Preprint / Version 1

Unconditionally stable symplectic integrators for the Navier-Stokes equations and other dissipative systems

##article.authors##

  • Sutthikiat Sungkeetanon The Citadel, The Military College of South Carolina https://orcid.org/0009-0007-0099-960X
  • Joseph Gaglione The Citadel, The Military College of South Carolina https://orcid.org/0009-0004-6210-0944
  • Robert Chapman The Citadel, The Military College of South Carolina
  • Tyler Kelly The Citadel, The Military College of South Carolina
  • Howard Cushman The Citadel, The Military College of South Carolina https://orcid.org/0009-0005-2681-0831
  • Blakeley Odom The Citadel, The Military College of South Carolina
  • Bryan MacGavin The Citadel, The Military College of South Carolina
  • Gafar Elamin The Citadel, The Military College of South Carolina https://orcid.org/0009-0009-1055-9688
  • Nathan Washuta The Citadel, The Military College of South Carolina
  • Jonathan Crosmer The Citadel, The Military College of South Carolina
  • Adam DeVoria The Citadel, The Military College of South Carolina https://orcid.org/0000-0001-5615-0807
  • John Sanders The Citadel, the Military College of South Carolina https://orcid.org/0000-0003-3059-3815

DOI:

https://doi.org/10.31224/4107

Keywords:

Navier-Stokes Equations, Symplectic Integrators, Variational Principles, Hamiltonian Mechanics

Abstract

Symplectic integrators offer vastly superior performance over traditional numerical techniques for conservative dynamical systems, but their application to \emph{dissipative} systems is inherently difficult due to dissipative systems' lack of symplectic structure. Leveraging the intrinsic variational structure of higher-order dynamics, this paper presents a general technique for applying existing symplectic integration schemes to dissipative systems, with particular emphasis on viscous fluids modeled by the Navier-Stokes equations. Two very simple such schemes are developed here. Not only are these schemes unconditionally stable for dissipative systems, they also outperform traditional methods with a similar degree of complexity in terms of accuracy for a given time step. For example, in the case of viscous flow between two infinite, flat plates, one of the schemes developed here is found to outperform both the implicit Euler method and the explicit fourth-order Runge-Kutta method in predicting the velocity profile. To the authors' knowledge, this is the very first time that a symplectic integration scheme has been applied successfully to the Navier-Stokes equations. We interpret the present success as direct empirical validation of the canonical Hamiltonian formulation of the Navier-Stokes problem recently published by Sanders~\emph{et al.} More sophisticated symplectic integration schemes are expected to exhibit even greater performance. It is hoped that these results will lead to improved numerical methods in computational fluid dynamics.

Downloads

Download data is not yet available.

Downloads

Posted

2024-11-12