A numerical study of an operator splitting method for rotating flows with large ageostrophic initial data

Donald Jones, Alex Mahalov, Basil Nicolaenko

Research output: Contribution to journalArticle

5 Citations (Scopus)

Abstract

We propose an operator splitting method which is especially suitable for long-time integration of geophysical equations characterized by the presence of multiple-time scales and weak-operator splitting. The method is illustrated on the classical rotating shallow-water equations on a periodic domain with large ageostrophic (unprepared) initial data. The asymptotic splitting decomposes the solution into a first part which solves the quasigeostrophic equation; a second one which is the "slow" ageostrophic component of the flow; and a corrector. The particular decomposition we use ensures that the corrector is small for large rotation. By considering only the "slow" ageostrophic and quasigeostrophic components a numerical approximation to the shallow-water equations is derived that effectively removes the time-step restrictions caused by the presence of fast waves. The splitting is exact in the asymptotic limit of large rotation and includes the nonlinearity of the equations. Numerical examples are included. These examples demonstrate a significant reduction in the computational cost over direct numerical approximations of the shallow-water equations. We conclude with an outline of a general operator splitting method for more general primitive geophysical equations.

Original languageEnglish (US)
Pages (from-to)143-159
Number of pages17
JournalTheoretical and Computational Fluid Dynamics
Volume13
Issue number2
StatePublished - 1999

Fingerprint

operators
Water
shallow water
Decomposition
Costs
approximation
constrictions
nonlinearity
costs
decomposition

ASJC Scopus subject areas

  • Computational Mechanics
  • Mechanics of Materials
  • Physics and Astronomy(all)
  • Condensed Matter Physics

Cite this

@article{af779f9d8fbc44c2b8e5e41368fa9771,
title = "A numerical study of an operator splitting method for rotating flows with large ageostrophic initial data",
abstract = "We propose an operator splitting method which is especially suitable for long-time integration of geophysical equations characterized by the presence of multiple-time scales and weak-operator splitting. The method is illustrated on the classical rotating shallow-water equations on a periodic domain with large ageostrophic (unprepared) initial data. The asymptotic splitting decomposes the solution into a first part which solves the quasigeostrophic equation; a second one which is the {"}slow{"} ageostrophic component of the flow; and a corrector. The particular decomposition we use ensures that the corrector is small for large rotation. By considering only the {"}slow{"} ageostrophic and quasigeostrophic components a numerical approximation to the shallow-water equations is derived that effectively removes the time-step restrictions caused by the presence of fast waves. The splitting is exact in the asymptotic limit of large rotation and includes the nonlinearity of the equations. Numerical examples are included. These examples demonstrate a significant reduction in the computational cost over direct numerical approximations of the shallow-water equations. We conclude with an outline of a general operator splitting method for more general primitive geophysical equations.",
author = "Donald Jones and Alex Mahalov and Basil Nicolaenko",
year = "1999",
language = "English (US)",
volume = "13",
pages = "143--159",
journal = "Theoretical and Computational Fluid Dynamics",
issn = "0935-4964",
publisher = "Springer New York",
number = "2",

}

TY - JOUR

T1 - A numerical study of an operator splitting method for rotating flows with large ageostrophic initial data

AU - Jones, Donald

AU - Mahalov, Alex

AU - Nicolaenko, Basil

PY - 1999

Y1 - 1999

N2 - We propose an operator splitting method which is especially suitable for long-time integration of geophysical equations characterized by the presence of multiple-time scales and weak-operator splitting. The method is illustrated on the classical rotating shallow-water equations on a periodic domain with large ageostrophic (unprepared) initial data. The asymptotic splitting decomposes the solution into a first part which solves the quasigeostrophic equation; a second one which is the "slow" ageostrophic component of the flow; and a corrector. The particular decomposition we use ensures that the corrector is small for large rotation. By considering only the "slow" ageostrophic and quasigeostrophic components a numerical approximation to the shallow-water equations is derived that effectively removes the time-step restrictions caused by the presence of fast waves. The splitting is exact in the asymptotic limit of large rotation and includes the nonlinearity of the equations. Numerical examples are included. These examples demonstrate a significant reduction in the computational cost over direct numerical approximations of the shallow-water equations. We conclude with an outline of a general operator splitting method for more general primitive geophysical equations.

AB - We propose an operator splitting method which is especially suitable for long-time integration of geophysical equations characterized by the presence of multiple-time scales and weak-operator splitting. The method is illustrated on the classical rotating shallow-water equations on a periodic domain with large ageostrophic (unprepared) initial data. The asymptotic splitting decomposes the solution into a first part which solves the quasigeostrophic equation; a second one which is the "slow" ageostrophic component of the flow; and a corrector. The particular decomposition we use ensures that the corrector is small for large rotation. By considering only the "slow" ageostrophic and quasigeostrophic components a numerical approximation to the shallow-water equations is derived that effectively removes the time-step restrictions caused by the presence of fast waves. The splitting is exact in the asymptotic limit of large rotation and includes the nonlinearity of the equations. Numerical examples are included. These examples demonstrate a significant reduction in the computational cost over direct numerical approximations of the shallow-water equations. We conclude with an outline of a general operator splitting method for more general primitive geophysical equations.

UR - http://www.scopus.com/inward/record.url?scp=0033238477&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=0033238477&partnerID=8YFLogxK

M3 - Article

AN - SCOPUS:0033238477

VL - 13

SP - 143

EP - 159

JO - Theoretical and Computational Fluid Dynamics

JF - Theoretical and Computational Fluid Dynamics

SN - 0935-4964

IS - 2

ER -