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

Fuente: arXiv
Saved in:
Bibliographic Details
Main Authors: Sungkeetanon, Sutthikiat, Gaglione, Joseph S., Chapman, Robert L., Kelly, Tyler M., Cushman, Howard A., Odom, Blakeley H., MacGavin, Bryan, Elamin, Gafar A., Washuta, Nathan J., Crosmer, Jonathan E., DeVoria, Adam C., Sanders, John W.
Format: Preprint
Published: 2024
Subjects:
Online Access:
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866915028729331712
author Sungkeetanon, Sutthikiat
Gaglione, Joseph S.
Chapman, Robert L.
Kelly, Tyler M.
Cushman, Howard A.
Odom, Blakeley H.
MacGavin, Bryan
Elamin, Gafar A.
Washuta, Nathan J.
Crosmer, Jonathan E.
DeVoria, Adam C.
Sanders, John W.
author_facet Sungkeetanon, Sutthikiat
Gaglione, Joseph S.
Chapman, Robert L.
Kelly, Tyler M.
Cushman, Howard A.
Odom, Blakeley H.
MacGavin, Bryan
Elamin, Gafar A.
Washuta, Nathan J.
Crosmer, Jonathan E.
DeVoria, Adam C.
Sanders, John W.
contents 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.
format Preprint
id arxiv_https___arxiv_org_abs_2411_13569
institution arXiv
publishDate 2024
record_format arxiv
spellingShingle Unconditionally stable symplectic integrators for the Navier-Stokes equations and other dissipative systems
Sungkeetanon, Sutthikiat
Gaglione, Joseph S.
Chapman, Robert L.
Kelly, Tyler M.
Cushman, Howard A.
Odom, Blakeley H.
MacGavin, Bryan
Elamin, Gafar A.
Washuta, Nathan J.
Crosmer, Jonathan E.
DeVoria, Adam C.
Sanders, John W.
Numerical Analysis
Dynamical Systems
Fluid Dynamics
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.
title Unconditionally stable symplectic integrators for the Navier-Stokes equations and other dissipative systems
topic Numerical Analysis
Dynamical Systems
Fluid Dynamics
url https://arxiv.org/abs/2411.13569