A robust numerical method for the potential vorticity based control variable transform in variational data assimilation

Sean Buckeridge, M J P Cullen, Robert Scheichl, M Wlasak

Research output: Contribution to journalArticle

  • 2 Citations

Abstract

The potential vorticity based control variable transformation for variational data assimilation, proposed in Cullen (2003), is a promising alternative to the currently more common vorticity based transformation. It leads to a better decorrelation of the control variables, but it involves solving a highly ill-conditioned elliptic partial differential equation (PDE), with a constraint. This PDE has so far been impossible to solve to any reasonable accuracy for realistic grid resolutions in finite difference formulations. Following on from the work in Buckeridge and Scheichl (2010) we propose a numerical method for it based on a Krylov subspace method with a multigrid preconditioner. The problem of interest includes a constraint in the form of two-dimensional elliptic solves embedded within the main three-dimensional problem. Thus the discretised problem cannot be formulated as a simple linear equation system with a sparse system matrix (as usual in elliptic PDEs). Therefore, in order to precondition the system we apply the multigrid method in Buckeridge and Scheichl (2010) to a simplified form of the three-dimensional operator (without the embedded two-dimensional problems) leading to an asymptotically optimal convergence of the preconditioned Krylov subspace method. The solvers used at the Met Office typically take over 100 iterations to converge to a residual tolerance of 0.1 and fail to converge to a tolerance of 10(-2). The method proposed in this paper, in contrast, can converge to a tolerance of 10(-2) within 15 iterations on all typical grid resolutions used at the Met Office, and is convergent to a tolerance of 10(-6). In addition, the method demonstrates almost optimal parallel scalability. Copyright
LanguageEnglish
Pages1083-1094
Number of pages12
JournalQuarterly Journal of the Royal Meteorological Society
Volume137
Issue number657
DOIs
StatusPublished - Apr 2011

Fingerprint

potential vorticity
data assimilation
numerical method
transform
tolerance
vorticity
method
matrix

Keywords

  • preconditioner
  • solver
  • multigrid
  • parallelisation
  • finite volume
  • potential vorticity

Cite this

A robust numerical method for the potential vorticity based control variable transform in variational data assimilation. / Buckeridge, Sean; Cullen, M J P; Scheichl, Robert; Wlasak, M.

In: Quarterly Journal of the Royal Meteorological Society, Vol. 137, No. 657, 04.2011, p. 1083-1094.

Research output: Contribution to journalArticle

@article{7e713df4c1524d349d4fa90e4c144ffa,
title = "A robust numerical method for the potential vorticity based control variable transform in variational data assimilation",
abstract = "The potential vorticity based control variable transformation for variational data assimilation, proposed in Cullen (2003), is a promising alternative to the currently more common vorticity based transformation. It leads to a better decorrelation of the control variables, but it involves solving a highly ill-conditioned elliptic partial differential equation (PDE), with a constraint. This PDE has so far been impossible to solve to any reasonable accuracy for realistic grid resolutions in finite difference formulations. Following on from the work in Buckeridge and Scheichl (2010) we propose a numerical method for it based on a Krylov subspace method with a multigrid preconditioner. The problem of interest includes a constraint in the form of two-dimensional elliptic solves embedded within the main three-dimensional problem. Thus the discretised problem cannot be formulated as a simple linear equation system with a sparse system matrix (as usual in elliptic PDEs). Therefore, in order to precondition the system we apply the multigrid method in Buckeridge and Scheichl (2010) to a simplified form of the three-dimensional operator (without the embedded two-dimensional problems) leading to an asymptotically optimal convergence of the preconditioned Krylov subspace method. The solvers used at the Met Office typically take over 100 iterations to converge to a residual tolerance of 0.1 and fail to converge to a tolerance of 10(-2). The method proposed in this paper, in contrast, can converge to a tolerance of 10(-2) within 15 iterations on all typical grid resolutions used at the Met Office, and is convergent to a tolerance of 10(-6). In addition, the method demonstrates almost optimal parallel scalability. Copyright",
keywords = "preconditioner, solver, multigrid, parallelisation, finite volume, potential vorticity",
author = "Sean Buckeridge and Cullen, {M J P} and Robert Scheichl and M Wlasak",
year = "2011",
month = "4",
doi = "10.1002/qj.826",
language = "English",
volume = "137",
pages = "1083--1094",
journal = "Quarterly Journal of the Royal Meteorological Society",
issn = "0035-9009",
publisher = "Wiley-Blackwell",
number = "657",

}

TY - JOUR

T1 - A robust numerical method for the potential vorticity based control variable transform in variational data assimilation

AU - Buckeridge,Sean

AU - Cullen,M J P

AU - Scheichl,Robert

AU - Wlasak,M

PY - 2011/4

Y1 - 2011/4

N2 - The potential vorticity based control variable transformation for variational data assimilation, proposed in Cullen (2003), is a promising alternative to the currently more common vorticity based transformation. It leads to a better decorrelation of the control variables, but it involves solving a highly ill-conditioned elliptic partial differential equation (PDE), with a constraint. This PDE has so far been impossible to solve to any reasonable accuracy for realistic grid resolutions in finite difference formulations. Following on from the work in Buckeridge and Scheichl (2010) we propose a numerical method for it based on a Krylov subspace method with a multigrid preconditioner. The problem of interest includes a constraint in the form of two-dimensional elliptic solves embedded within the main three-dimensional problem. Thus the discretised problem cannot be formulated as a simple linear equation system with a sparse system matrix (as usual in elliptic PDEs). Therefore, in order to precondition the system we apply the multigrid method in Buckeridge and Scheichl (2010) to a simplified form of the three-dimensional operator (without the embedded two-dimensional problems) leading to an asymptotically optimal convergence of the preconditioned Krylov subspace method. The solvers used at the Met Office typically take over 100 iterations to converge to a residual tolerance of 0.1 and fail to converge to a tolerance of 10(-2). The method proposed in this paper, in contrast, can converge to a tolerance of 10(-2) within 15 iterations on all typical grid resolutions used at the Met Office, and is convergent to a tolerance of 10(-6). In addition, the method demonstrates almost optimal parallel scalability. Copyright

AB - The potential vorticity based control variable transformation for variational data assimilation, proposed in Cullen (2003), is a promising alternative to the currently more common vorticity based transformation. It leads to a better decorrelation of the control variables, but it involves solving a highly ill-conditioned elliptic partial differential equation (PDE), with a constraint. This PDE has so far been impossible to solve to any reasonable accuracy for realistic grid resolutions in finite difference formulations. Following on from the work in Buckeridge and Scheichl (2010) we propose a numerical method for it based on a Krylov subspace method with a multigrid preconditioner. The problem of interest includes a constraint in the form of two-dimensional elliptic solves embedded within the main three-dimensional problem. Thus the discretised problem cannot be formulated as a simple linear equation system with a sparse system matrix (as usual in elliptic PDEs). Therefore, in order to precondition the system we apply the multigrid method in Buckeridge and Scheichl (2010) to a simplified form of the three-dimensional operator (without the embedded two-dimensional problems) leading to an asymptotically optimal convergence of the preconditioned Krylov subspace method. The solvers used at the Met Office typically take over 100 iterations to converge to a residual tolerance of 0.1 and fail to converge to a tolerance of 10(-2). The method proposed in this paper, in contrast, can converge to a tolerance of 10(-2) within 15 iterations on all typical grid resolutions used at the Met Office, and is convergent to a tolerance of 10(-6). In addition, the method demonstrates almost optimal parallel scalability. Copyright

KW - preconditioner

KW - solver

KW - multigrid

KW - parallelisation

KW - finite volume

KW - potential vorticity

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

UR - http://dx.doi.org/10.1002/qj.826

U2 - 10.1002/qj.826

DO - 10.1002/qj.826

M3 - Article

VL - 137

SP - 1083

EP - 1094

JO - Quarterly Journal of the Royal Meteorological Society

T2 - Quarterly Journal of the Royal Meteorological Society

JF - Quarterly Journal of the Royal Meteorological Society

SN - 0035-9009

IS - 657

ER -