SEAWAT 2000: Modelling unstable flow and sensitivity to discretization levels and numerical schemes

A. Al-Maktoumi, D. A. Lockington, R. E. Volker

Research output: Contribution to journalArticle

5 Citations (Scopus)

Abstract

A systematic analysis shows how results from the finite difference code SEAWAT are sensitive to choice of grid dimension, time step, and numerical scheme for unstable flow problems. Guidelines to assist in selecting appropriate combinations of these factors are suggested. While the SEAWAT code has been tested for a wide range of problems, the sensitivity of results to spatial and temporal discretization levels and numerical schemes has not been studied in detail for unstable flow problems. Here, the Elder-Voss-Souza benchmark problem has been used to systematically explore the sensitivity of SEAWAT output to spatio-temporal resolution and numerical solver choice. A grid size of 0.38 and 0.60% of the total domain length and depth respectively is found to be fine enough to deliver results with acceptable accuracy for most of the numerical schemes when Courant number (Cr) is 0.1. All numerical solvers produced similar results for extremely fine meshes; however, some schemes converged faster than others. For instance, the 3rd-order total variation-diminishing method (TVD3) scheme converged at a much coarser mesh than the standard finite difference methods (SFDM) upstream weighting (UW) scheme. The sensitivity of the results to Cr number depends on the numerical scheme as expected.

Original languageEnglish
Pages (from-to)1119-1129
Number of pages11
JournalHydrogeology Journal
Volume15
Issue number6
DOIs
Publication statusPublished - Sep 2007

Fingerprint

flow modeling
finite difference method
code
analysis
method

Keywords

  • Convection
  • Grid sensitivity
  • Numerical modeling
  • SEAWAT
  • Unstable flow

ASJC Scopus subject areas

  • Earth and Planetary Sciences (miscellaneous)
  • Water Science and Technology

Cite this

SEAWAT 2000 : Modelling unstable flow and sensitivity to discretization levels and numerical schemes. / Al-Maktoumi, A.; Lockington, D. A.; Volker, R. E.

In: Hydrogeology Journal, Vol. 15, No. 6, 09.2007, p. 1119-1129.

Research output: Contribution to journalArticle

@article{d884ff329d484ad28a13e8069d883144,
title = "SEAWAT 2000: Modelling unstable flow and sensitivity to discretization levels and numerical schemes",
abstract = "A systematic analysis shows how results from the finite difference code SEAWAT are sensitive to choice of grid dimension, time step, and numerical scheme for unstable flow problems. Guidelines to assist in selecting appropriate combinations of these factors are suggested. While the SEAWAT code has been tested for a wide range of problems, the sensitivity of results to spatial and temporal discretization levels and numerical schemes has not been studied in detail for unstable flow problems. Here, the Elder-Voss-Souza benchmark problem has been used to systematically explore the sensitivity of SEAWAT output to spatio-temporal resolution and numerical solver choice. A grid size of 0.38 and 0.60{\%} of the total domain length and depth respectively is found to be fine enough to deliver results with acceptable accuracy for most of the numerical schemes when Courant number (Cr) is 0.1. All numerical solvers produced similar results for extremely fine meshes; however, some schemes converged faster than others. For instance, the 3rd-order total variation-diminishing method (TVD3) scheme converged at a much coarser mesh than the standard finite difference methods (SFDM) upstream weighting (UW) scheme. The sensitivity of the results to Cr number depends on the numerical scheme as expected.",
keywords = "Convection, Grid sensitivity, Numerical modeling, SEAWAT, Unstable flow",
author = "A. Al-Maktoumi and Lockington, {D. A.} and Volker, {R. E.}",
year = "2007",
month = "9",
doi = "10.1007/s10040-007-0164-2",
language = "English",
volume = "15",
pages = "1119--1129",
journal = "Hydrogeology Journal",
issn = "1431-2174",
publisher = "Springer Heidelberg",
number = "6",

}

TY - JOUR

T1 - SEAWAT 2000

T2 - Modelling unstable flow and sensitivity to discretization levels and numerical schemes

AU - Al-Maktoumi, A.

AU - Lockington, D. A.

AU - Volker, R. E.

PY - 2007/9

Y1 - 2007/9

N2 - A systematic analysis shows how results from the finite difference code SEAWAT are sensitive to choice of grid dimension, time step, and numerical scheme for unstable flow problems. Guidelines to assist in selecting appropriate combinations of these factors are suggested. While the SEAWAT code has been tested for a wide range of problems, the sensitivity of results to spatial and temporal discretization levels and numerical schemes has not been studied in detail for unstable flow problems. Here, the Elder-Voss-Souza benchmark problem has been used to systematically explore the sensitivity of SEAWAT output to spatio-temporal resolution and numerical solver choice. A grid size of 0.38 and 0.60% of the total domain length and depth respectively is found to be fine enough to deliver results with acceptable accuracy for most of the numerical schemes when Courant number (Cr) is 0.1. All numerical solvers produced similar results for extremely fine meshes; however, some schemes converged faster than others. For instance, the 3rd-order total variation-diminishing method (TVD3) scheme converged at a much coarser mesh than the standard finite difference methods (SFDM) upstream weighting (UW) scheme. The sensitivity of the results to Cr number depends on the numerical scheme as expected.

AB - A systematic analysis shows how results from the finite difference code SEAWAT are sensitive to choice of grid dimension, time step, and numerical scheme for unstable flow problems. Guidelines to assist in selecting appropriate combinations of these factors are suggested. While the SEAWAT code has been tested for a wide range of problems, the sensitivity of results to spatial and temporal discretization levels and numerical schemes has not been studied in detail for unstable flow problems. Here, the Elder-Voss-Souza benchmark problem has been used to systematically explore the sensitivity of SEAWAT output to spatio-temporal resolution and numerical solver choice. A grid size of 0.38 and 0.60% of the total domain length and depth respectively is found to be fine enough to deliver results with acceptable accuracy for most of the numerical schemes when Courant number (Cr) is 0.1. All numerical solvers produced similar results for extremely fine meshes; however, some schemes converged faster than others. For instance, the 3rd-order total variation-diminishing method (TVD3) scheme converged at a much coarser mesh than the standard finite difference methods (SFDM) upstream weighting (UW) scheme. The sensitivity of the results to Cr number depends on the numerical scheme as expected.

KW - Convection

KW - Grid sensitivity

KW - Numerical modeling

KW - SEAWAT

KW - Unstable flow

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

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

U2 - 10.1007/s10040-007-0164-2

DO - 10.1007/s10040-007-0164-2

M3 - Article

AN - SCOPUS:34548490124

VL - 15

SP - 1119

EP - 1129

JO - Hydrogeology Journal

JF - Hydrogeology Journal

SN - 1431-2174

IS - 6

ER -