Modelling of density-dependent flow systems: Sensitivity to spatial and temporal discretizations and numerical schemes using seawat 2000

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

Research output: Chapter in Book/Report/Conference proceedingChapter

Abstract

Contaminant transport with significant density contrasts is of increasing interest in many subsurface-hydrology problems like seawater intrusion, saltwater up-coning in coastal aquifers, and dense contaminant migration. Several numerical codes have been developed to simulate these systems. Results from these codes are found to be sensitive to grid and time step sizes, and to the numerical schemes used to approximate the solution. Although the 3-D Finite Difference Method code, SEAWAT has been used over a wide range of problems, the sensitivity of the computed results to spatial and temporal discretization levels for different numerical schemes supported by the code has not been studied in detail. This study aims to: (1) provide an approximate guide for SEAWAT users to ensure proper selection of grid and time step sizes for a particular numerical scheme in order to minimize the numerical error (in 2-D) and (2) To investigate the ability of the SEAWAT code to simulate more complex 3-D systems. The Elder-Voss benchmark problem was selected for the study. Different levels of grid and time discretizations were found to produce significantly different results. A grid size with dx (0.38%) and dz each 0.6% of the total length and depth of the domain respectively is found to be fine enough to produce results with acceptable accuracy for most of the numerical schemes with Courant number (Cr) of 0.1. Some numerical schemes produce accurate results at coarser meshes compared with other schemes that produce similar accuracy at only an extremely fine mesh. For the frequently used Cr value of 1.0, the computed salinity-spread patterns varied significantly for different solution schemes compared with the case when Cr is 0.1 or 0.5. To ensure a high level of accuracy in the modeling results using SEAWAT, Cr should be ≤0.1 when the Peclet number is ≤ about 1. SEAWAT was able to capture the main physical features of the Elder-Voss convection pattern in 3-D at coarse mesh sizes compared with that produced by a Finite Element Method code “FEFLOW” at a very fine mesh. The 3-D modelling requires large numerical efforts compared with that of the 2-D case. Runtimes for the 2-D problem range from less than. hr for a coarser mesh to about 2.8 days in the case of a very fine mesh, whereas it takes around 6 days in the 3-D case with a relatively coarse resolution level. Similarity between 2-D and 3-D convection patterns was found to exist in terms of fingering, and up/down welling behaviours. Other factors like the level of accuracy required, computational expense and storage memory requirements should be considered along with the spatial and temporal resolution levels.

Original languageEnglish
Title of host publicationGroundwater for Sustainable Development
Subtitle of host publicationProblems, Perspectives and Challenges
PublisherCRC Press
Pages155-167
Number of pages13
ISBN (Electronic)9780203894569
ISBN (Print)9780415407762
Publication statusPublished - Jan 1 2008

Fingerprint

Salt water intrusion
Impurities
Peclet number
Hydrology
Aquifers
Finite difference method
modeling
Finite element method
Data storage equipment
convection
fingering
saline intrusion
coastal aquifer
mesh size
pollutant transport
finite difference method
finite element method
Convection
code
hydrology

ASJC Scopus subject areas

  • Engineering(all)
  • Environmental Science(all)

Cite this

Al Maktoumi, A., Lockington, D. A., & Volker, R. E. (2008). Modelling of density-dependent flow systems: Sensitivity to spatial and temporal discretizations and numerical schemes using seawat 2000. In Groundwater for Sustainable Development: Problems, Perspectives and Challenges (pp. 155-167). CRC Press.

Modelling of density-dependent flow systems : Sensitivity to spatial and temporal discretizations and numerical schemes using seawat 2000. / Al Maktoumi, Ali; Lockington, D. A.; Volker, R. E.

Groundwater for Sustainable Development: Problems, Perspectives and Challenges. CRC Press, 2008. p. 155-167.

Research output: Chapter in Book/Report/Conference proceedingChapter

Al Maktoumi, A, Lockington, DA & Volker, RE 2008, Modelling of density-dependent flow systems: Sensitivity to spatial and temporal discretizations and numerical schemes using seawat 2000. in Groundwater for Sustainable Development: Problems, Perspectives and Challenges. CRC Press, pp. 155-167.
Al Maktoumi A, Lockington DA, Volker RE. Modelling of density-dependent flow systems: Sensitivity to spatial and temporal discretizations and numerical schemes using seawat 2000. In Groundwater for Sustainable Development: Problems, Perspectives and Challenges. CRC Press. 2008. p. 155-167
Al Maktoumi, Ali ; Lockington, D. A. ; Volker, R. E. / Modelling of density-dependent flow systems : Sensitivity to spatial and temporal discretizations and numerical schemes using seawat 2000. Groundwater for Sustainable Development: Problems, Perspectives and Challenges. CRC Press, 2008. pp. 155-167
@inbook{fb03bc4203024cd493ef24d8410e82c9,
title = "Modelling of density-dependent flow systems: Sensitivity to spatial and temporal discretizations and numerical schemes using seawat 2000",
abstract = "Contaminant transport with significant density contrasts is of increasing interest in many subsurface-hydrology problems like seawater intrusion, saltwater up-coning in coastal aquifers, and dense contaminant migration. Several numerical codes have been developed to simulate these systems. Results from these codes are found to be sensitive to grid and time step sizes, and to the numerical schemes used to approximate the solution. Although the 3-D Finite Difference Method code, SEAWAT has been used over a wide range of problems, the sensitivity of the computed results to spatial and temporal discretization levels for different numerical schemes supported by the code has not been studied in detail. This study aims to: (1) provide an approximate guide for SEAWAT users to ensure proper selection of grid and time step sizes for a particular numerical scheme in order to minimize the numerical error (in 2-D) and (2) To investigate the ability of the SEAWAT code to simulate more complex 3-D systems. The Elder-Voss benchmark problem was selected for the study. Different levels of grid and time discretizations were found to produce significantly different results. A grid size with dx (0.38{\%}) and dz each 0.6{\%} of the total length and depth of the domain respectively is found to be fine enough to produce results with acceptable accuracy for most of the numerical schemes with Courant number (Cr) of 0.1. Some numerical schemes produce accurate results at coarser meshes compared with other schemes that produce similar accuracy at only an extremely fine mesh. For the frequently used Cr value of 1.0, the computed salinity-spread patterns varied significantly for different solution schemes compared with the case when Cr is 0.1 or 0.5. To ensure a high level of accuracy in the modeling results using SEAWAT, Cr should be ≤0.1 when the Peclet number is ≤ about 1. SEAWAT was able to capture the main physical features of the Elder-Voss convection pattern in 3-D at coarse mesh sizes compared with that produced by a Finite Element Method code “FEFLOW” at a very fine mesh. The 3-D modelling requires large numerical efforts compared with that of the 2-D case. Runtimes for the 2-D problem range from less than. hr for a coarser mesh to about 2.8 days in the case of a very fine mesh, whereas it takes around 6 days in the 3-D case with a relatively coarse resolution level. Similarity between 2-D and 3-D convection patterns was found to exist in terms of fingering, and up/down welling behaviours. Other factors like the level of accuracy required, computational expense and storage memory requirements should be considered along with the spatial and temporal resolution levels.",
author = "{Al Maktoumi}, Ali and Lockington, {D. A.} and Volker, {R. E.}",
year = "2008",
month = "1",
day = "1",
language = "English",
isbn = "9780415407762",
pages = "155--167",
booktitle = "Groundwater for Sustainable Development",
publisher = "CRC Press",

}

TY - CHAP

T1 - Modelling of density-dependent flow systems

T2 - Sensitivity to spatial and temporal discretizations and numerical schemes using seawat 2000

AU - Al Maktoumi, Ali

AU - Lockington, D. A.

AU - Volker, R. E.

PY - 2008/1/1

Y1 - 2008/1/1

N2 - Contaminant transport with significant density contrasts is of increasing interest in many subsurface-hydrology problems like seawater intrusion, saltwater up-coning in coastal aquifers, and dense contaminant migration. Several numerical codes have been developed to simulate these systems. Results from these codes are found to be sensitive to grid and time step sizes, and to the numerical schemes used to approximate the solution. Although the 3-D Finite Difference Method code, SEAWAT has been used over a wide range of problems, the sensitivity of the computed results to spatial and temporal discretization levels for different numerical schemes supported by the code has not been studied in detail. This study aims to: (1) provide an approximate guide for SEAWAT users to ensure proper selection of grid and time step sizes for a particular numerical scheme in order to minimize the numerical error (in 2-D) and (2) To investigate the ability of the SEAWAT code to simulate more complex 3-D systems. The Elder-Voss benchmark problem was selected for the study. Different levels of grid and time discretizations were found to produce significantly different results. A grid size with dx (0.38%) and dz each 0.6% of the total length and depth of the domain respectively is found to be fine enough to produce results with acceptable accuracy for most of the numerical schemes with Courant number (Cr) of 0.1. Some numerical schemes produce accurate results at coarser meshes compared with other schemes that produce similar accuracy at only an extremely fine mesh. For the frequently used Cr value of 1.0, the computed salinity-spread patterns varied significantly for different solution schemes compared with the case when Cr is 0.1 or 0.5. To ensure a high level of accuracy in the modeling results using SEAWAT, Cr should be ≤0.1 when the Peclet number is ≤ about 1. SEAWAT was able to capture the main physical features of the Elder-Voss convection pattern in 3-D at coarse mesh sizes compared with that produced by a Finite Element Method code “FEFLOW” at a very fine mesh. The 3-D modelling requires large numerical efforts compared with that of the 2-D case. Runtimes for the 2-D problem range from less than. hr for a coarser mesh to about 2.8 days in the case of a very fine mesh, whereas it takes around 6 days in the 3-D case with a relatively coarse resolution level. Similarity between 2-D and 3-D convection patterns was found to exist in terms of fingering, and up/down welling behaviours. Other factors like the level of accuracy required, computational expense and storage memory requirements should be considered along with the spatial and temporal resolution levels.

AB - Contaminant transport with significant density contrasts is of increasing interest in many subsurface-hydrology problems like seawater intrusion, saltwater up-coning in coastal aquifers, and dense contaminant migration. Several numerical codes have been developed to simulate these systems. Results from these codes are found to be sensitive to grid and time step sizes, and to the numerical schemes used to approximate the solution. Although the 3-D Finite Difference Method code, SEAWAT has been used over a wide range of problems, the sensitivity of the computed results to spatial and temporal discretization levels for different numerical schemes supported by the code has not been studied in detail. This study aims to: (1) provide an approximate guide for SEAWAT users to ensure proper selection of grid and time step sizes for a particular numerical scheme in order to minimize the numerical error (in 2-D) and (2) To investigate the ability of the SEAWAT code to simulate more complex 3-D systems. The Elder-Voss benchmark problem was selected for the study. Different levels of grid and time discretizations were found to produce significantly different results. A grid size with dx (0.38%) and dz each 0.6% of the total length and depth of the domain respectively is found to be fine enough to produce results with acceptable accuracy for most of the numerical schemes with Courant number (Cr) of 0.1. Some numerical schemes produce accurate results at coarser meshes compared with other schemes that produce similar accuracy at only an extremely fine mesh. For the frequently used Cr value of 1.0, the computed salinity-spread patterns varied significantly for different solution schemes compared with the case when Cr is 0.1 or 0.5. To ensure a high level of accuracy in the modeling results using SEAWAT, Cr should be ≤0.1 when the Peclet number is ≤ about 1. SEAWAT was able to capture the main physical features of the Elder-Voss convection pattern in 3-D at coarse mesh sizes compared with that produced by a Finite Element Method code “FEFLOW” at a very fine mesh. The 3-D modelling requires large numerical efforts compared with that of the 2-D case. Runtimes for the 2-D problem range from less than. hr for a coarser mesh to about 2.8 days in the case of a very fine mesh, whereas it takes around 6 days in the 3-D case with a relatively coarse resolution level. Similarity between 2-D and 3-D convection patterns was found to exist in terms of fingering, and up/down welling behaviours. Other factors like the level of accuracy required, computational expense and storage memory requirements should be considered along with the spatial and temporal resolution levels.

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

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

M3 - Chapter

AN - SCOPUS:85056041343

SN - 9780415407762

SP - 155

EP - 167

BT - Groundwater for Sustainable Development

PB - CRC Press

ER -