Euler-Lagrange equations for the spectral element shallow water system

Julia C. Levin, Dale B. Haidvogel, B. Chua, A. F. Bennett, M. Iskandarani

Research output: Contribution to journalArticle

5 Citations (Scopus)

Abstract

We present the derivation of the discrete Euler-Lagrange equations for an inverse spectral element ocean model based on the shallow water equations. We show that the discrete Euler-Lagrange equations can be obtained from the continuous Euler-Lagrange equations by using a correct combination of the weak and the strong forms of derivatives in the Galerkin integrals, and by changing the order with which elemental assembly and mass averaging are applied in the forward and in the adjoint systems. Our derivation can be extended to obtain an adjoint for any Galerkin finite element and spectral element system. We begin the derivations using a linear wave equation in one dimension. We then apply our technique to a two-dimensional shallow water ocean model and test it on a classic double-gyre problem. The spectral element forward and adjoint ocean models can be used in a variety of inverse applications, ranging from traditional data assimilation and parameter estimation, to the less traditional model sensitivity and stability analyses, and ensemble prediction. Here the Euler-Lagrange equations are solved by an indirect representer algorithm.

Original languageEnglish (US)
Pages (from-to)348-377
Number of pages30
JournalOcean Modelling
Volume12
Issue number3-4
DOIs
StatePublished - Jan 12 2006

Fingerprint

shallow water
Water
ocean
shallow-water equation
wave equation
Wave equations
gyre
data assimilation
Parameter estimation
Derivatives
prediction

All Science Journal Classification (ASJC) codes

  • Computer Science (miscellaneous)
  • Oceanography
  • Geotechnical Engineering and Engineering Geology
  • Atmospheric Science

Keywords

  • 4DVAR variational data assimilation
  • Euler-Lagrange equations
  • Inverse ocean model
  • Spectral element method
  • Twin experiment

Cite this

Levin, Julia C. ; Haidvogel, Dale B. ; Chua, B. ; Bennett, A. F. ; Iskandarani, M. / Euler-Lagrange equations for the spectral element shallow water system. In: Ocean Modelling. 2006 ; Vol. 12, No. 3-4. pp. 348-377.
@article{658d9674c42146a98283aabcd3324382,
title = "Euler-Lagrange equations for the spectral element shallow water system",
abstract = "We present the derivation of the discrete Euler-Lagrange equations for an inverse spectral element ocean model based on the shallow water equations. We show that the discrete Euler-Lagrange equations can be obtained from the continuous Euler-Lagrange equations by using a correct combination of the weak and the strong forms of derivatives in the Galerkin integrals, and by changing the order with which elemental assembly and mass averaging are applied in the forward and in the adjoint systems. Our derivation can be extended to obtain an adjoint for any Galerkin finite element and spectral element system. We begin the derivations using a linear wave equation in one dimension. We then apply our technique to a two-dimensional shallow water ocean model and test it on a classic double-gyre problem. The spectral element forward and adjoint ocean models can be used in a variety of inverse applications, ranging from traditional data assimilation and parameter estimation, to the less traditional model sensitivity and stability analyses, and ensemble prediction. Here the Euler-Lagrange equations are solved by an indirect representer algorithm.",
keywords = "4DVAR variational data assimilation, Euler-Lagrange equations, Inverse ocean model, Spectral element method, Twin experiment",
author = "Levin, {Julia C.} and Haidvogel, {Dale B.} and B. Chua and Bennett, {A. F.} and M. Iskandarani",
year = "2006",
month = "1",
day = "12",
doi = "10.1016/j.ocemod.2005.06.002",
language = "English (US)",
volume = "12",
pages = "348--377",
journal = "Ocean Modelling",
issn = "1463-5003",
publisher = "Elsevier BV",
number = "3-4",

}

Euler-Lagrange equations for the spectral element shallow water system. / Levin, Julia C.; Haidvogel, Dale B.; Chua, B.; Bennett, A. F.; Iskandarani, M.

In: Ocean Modelling, Vol. 12, No. 3-4, 12.01.2006, p. 348-377.

Research output: Contribution to journalArticle

TY - JOUR

T1 - Euler-Lagrange equations for the spectral element shallow water system

AU - Levin, Julia C.

AU - Haidvogel, Dale B.

AU - Chua, B.

AU - Bennett, A. F.

AU - Iskandarani, M.

PY - 2006/1/12

Y1 - 2006/1/12

N2 - We present the derivation of the discrete Euler-Lagrange equations for an inverse spectral element ocean model based on the shallow water equations. We show that the discrete Euler-Lagrange equations can be obtained from the continuous Euler-Lagrange equations by using a correct combination of the weak and the strong forms of derivatives in the Galerkin integrals, and by changing the order with which elemental assembly and mass averaging are applied in the forward and in the adjoint systems. Our derivation can be extended to obtain an adjoint for any Galerkin finite element and spectral element system. We begin the derivations using a linear wave equation in one dimension. We then apply our technique to a two-dimensional shallow water ocean model and test it on a classic double-gyre problem. The spectral element forward and adjoint ocean models can be used in a variety of inverse applications, ranging from traditional data assimilation and parameter estimation, to the less traditional model sensitivity and stability analyses, and ensemble prediction. Here the Euler-Lagrange equations are solved by an indirect representer algorithm.

AB - We present the derivation of the discrete Euler-Lagrange equations for an inverse spectral element ocean model based on the shallow water equations. We show that the discrete Euler-Lagrange equations can be obtained from the continuous Euler-Lagrange equations by using a correct combination of the weak and the strong forms of derivatives in the Galerkin integrals, and by changing the order with which elemental assembly and mass averaging are applied in the forward and in the adjoint systems. Our derivation can be extended to obtain an adjoint for any Galerkin finite element and spectral element system. We begin the derivations using a linear wave equation in one dimension. We then apply our technique to a two-dimensional shallow water ocean model and test it on a classic double-gyre problem. The spectral element forward and adjoint ocean models can be used in a variety of inverse applications, ranging from traditional data assimilation and parameter estimation, to the less traditional model sensitivity and stability analyses, and ensemble prediction. Here the Euler-Lagrange equations are solved by an indirect representer algorithm.

KW - 4DVAR variational data assimilation

KW - Euler-Lagrange equations

KW - Inverse ocean model

KW - Spectral element method

KW - Twin experiment

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

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

U2 - 10.1016/j.ocemod.2005.06.002

DO - 10.1016/j.ocemod.2005.06.002

M3 - Article

AN - SCOPUS:30044434356

VL - 12

SP - 348

EP - 377

JO - Ocean Modelling

JF - Ocean Modelling

SN - 1463-5003

IS - 3-4

ER -