Numerical methods for conjunctive two-dimensional surface and three-dimensional sub-surface flows

Masaru Morita, Ben Chie Yen

Research output: Contribution to journalArticle

23 Citations (Scopus)

Abstract

Sophisticated catchment runoff problems necessitate conjunctive modeling of overland flow and sub-surface flow. In this paper, finite difference numerical methods are studied for simulation of catchment runoff of two-dimensional surface flow interacting with three-dimensional unsaturated and saturated sub-surface flows. The equations representing the flows are mathematically classified as a type of heat diffusion equation. Therefore, two- and three-dimensional numerical methods for heat diffusion equations were investigated for applications to the surface and sub-surface flow sub-models in terms of accuracy, stability, and calculation time. The methods are the purely explicit method, Saul'yev's methods, the alternating direction explicit (ADE) methods, and the alternating direction implicit (ADI) methods. The methods are first examined on surface and sub-surface flows separately; subsequently, 12 selected combinations of methods were investigated for modeling the conjunctive flows. Saul'yev's downstream (S-d) method was found to be the preferred method for two-dimensional surface flow modeling, whereas the ADE method of Barakat and Clark is a less accurate, stable alternative. For the three-dimensional sub-surface flow model, the ADE method of Larkin (ADE-L) and Brian's ADI method are unconditionally stable and more accurate than the other methods. The calculations of the conjunctive models utilizing the S-d surface flow sub-model give excellent results and confirm the expectation that the errors of the surface and sub-surface sub-models interact. The surface sub-model dominates the accuracy and stability of the conjunctive model, whereas the sub-surface sub-model dominates the calculation time, suggesting the desirability of using a smaller time increment for the surface sub-model. Copyright (C) 2000 John Wiley and Sons, Ltd.

Original languageEnglish
Pages (from-to)921-957
Number of pages37
JournalInternational Journal for Numerical Methods in Fluids
Volume32
Issue number8
DOIs
Publication statusPublished - 2000 Apr

Fingerprint

Numerical methods
alternating direction implicit methods
drainage
Runoff
Catchments
heat

Keywords

  • Conjunctive flow model
  • Heat diffusion equation
  • Sub-surface flow
  • Surface flow
  • Transient flow

ASJC Scopus subject areas

  • Computational Theory and Mathematics
  • Computer Science Applications
  • Computational Mechanics
  • Mechanics of Materials
  • Safety, Risk, Reliability and Quality
  • Applied Mathematics
  • Condensed Matter Physics

Cite this

Numerical methods for conjunctive two-dimensional surface and three-dimensional sub-surface flows. / Morita, Masaru; Yen, Ben Chie.

In: International Journal for Numerical Methods in Fluids, Vol. 32, No. 8, 04.2000, p. 921-957.

Research output: Contribution to journalArticle

@article{a93bc7ee607c424891d26fac5efafd6c,
title = "Numerical methods for conjunctive two-dimensional surface and three-dimensional sub-surface flows",
abstract = "Sophisticated catchment runoff problems necessitate conjunctive modeling of overland flow and sub-surface flow. In this paper, finite difference numerical methods are studied for simulation of catchment runoff of two-dimensional surface flow interacting with three-dimensional unsaturated and saturated sub-surface flows. The equations representing the flows are mathematically classified as a type of heat diffusion equation. Therefore, two- and three-dimensional numerical methods for heat diffusion equations were investigated for applications to the surface and sub-surface flow sub-models in terms of accuracy, stability, and calculation time. The methods are the purely explicit method, Saul'yev's methods, the alternating direction explicit (ADE) methods, and the alternating direction implicit (ADI) methods. The methods are first examined on surface and sub-surface flows separately; subsequently, 12 selected combinations of methods were investigated for modeling the conjunctive flows. Saul'yev's downstream (S-d) method was found to be the preferred method for two-dimensional surface flow modeling, whereas the ADE method of Barakat and Clark is a less accurate, stable alternative. For the three-dimensional sub-surface flow model, the ADE method of Larkin (ADE-L) and Brian's ADI method are unconditionally stable and more accurate than the other methods. The calculations of the conjunctive models utilizing the S-d surface flow sub-model give excellent results and confirm the expectation that the errors of the surface and sub-surface sub-models interact. The surface sub-model dominates the accuracy and stability of the conjunctive model, whereas the sub-surface sub-model dominates the calculation time, suggesting the desirability of using a smaller time increment for the surface sub-model. Copyright (C) 2000 John Wiley and Sons, Ltd.",
keywords = "Conjunctive flow model, Heat diffusion equation, Sub-surface flow, Surface flow, Transient flow",
author = "Masaru Morita and Yen, {Ben Chie}",
year = "2000",
month = "4",
doi = "10.1002/(SICI)1097-0363(20000430)32:8<921::AID-FLD993>3.0.CO;2-3",
language = "English",
volume = "32",
pages = "921--957",
journal = "International Journal for Numerical Methods in Fluids",
issn = "0271-2091",
publisher = "John Wiley and Sons Ltd",
number = "8",

}

TY - JOUR

T1 - Numerical methods for conjunctive two-dimensional surface and three-dimensional sub-surface flows

AU - Morita, Masaru

AU - Yen, Ben Chie

PY - 2000/4

Y1 - 2000/4

N2 - Sophisticated catchment runoff problems necessitate conjunctive modeling of overland flow and sub-surface flow. In this paper, finite difference numerical methods are studied for simulation of catchment runoff of two-dimensional surface flow interacting with three-dimensional unsaturated and saturated sub-surface flows. The equations representing the flows are mathematically classified as a type of heat diffusion equation. Therefore, two- and three-dimensional numerical methods for heat diffusion equations were investigated for applications to the surface and sub-surface flow sub-models in terms of accuracy, stability, and calculation time. The methods are the purely explicit method, Saul'yev's methods, the alternating direction explicit (ADE) methods, and the alternating direction implicit (ADI) methods. The methods are first examined on surface and sub-surface flows separately; subsequently, 12 selected combinations of methods were investigated for modeling the conjunctive flows. Saul'yev's downstream (S-d) method was found to be the preferred method for two-dimensional surface flow modeling, whereas the ADE method of Barakat and Clark is a less accurate, stable alternative. For the three-dimensional sub-surface flow model, the ADE method of Larkin (ADE-L) and Brian's ADI method are unconditionally stable and more accurate than the other methods. The calculations of the conjunctive models utilizing the S-d surface flow sub-model give excellent results and confirm the expectation that the errors of the surface and sub-surface sub-models interact. The surface sub-model dominates the accuracy and stability of the conjunctive model, whereas the sub-surface sub-model dominates the calculation time, suggesting the desirability of using a smaller time increment for the surface sub-model. Copyright (C) 2000 John Wiley and Sons, Ltd.

AB - Sophisticated catchment runoff problems necessitate conjunctive modeling of overland flow and sub-surface flow. In this paper, finite difference numerical methods are studied for simulation of catchment runoff of two-dimensional surface flow interacting with three-dimensional unsaturated and saturated sub-surface flows. The equations representing the flows are mathematically classified as a type of heat diffusion equation. Therefore, two- and three-dimensional numerical methods for heat diffusion equations were investigated for applications to the surface and sub-surface flow sub-models in terms of accuracy, stability, and calculation time. The methods are the purely explicit method, Saul'yev's methods, the alternating direction explicit (ADE) methods, and the alternating direction implicit (ADI) methods. The methods are first examined on surface and sub-surface flows separately; subsequently, 12 selected combinations of methods were investigated for modeling the conjunctive flows. Saul'yev's downstream (S-d) method was found to be the preferred method for two-dimensional surface flow modeling, whereas the ADE method of Barakat and Clark is a less accurate, stable alternative. For the three-dimensional sub-surface flow model, the ADE method of Larkin (ADE-L) and Brian's ADI method are unconditionally stable and more accurate than the other methods. The calculations of the conjunctive models utilizing the S-d surface flow sub-model give excellent results and confirm the expectation that the errors of the surface and sub-surface sub-models interact. The surface sub-model dominates the accuracy and stability of the conjunctive model, whereas the sub-surface sub-model dominates the calculation time, suggesting the desirability of using a smaller time increment for the surface sub-model. Copyright (C) 2000 John Wiley and Sons, Ltd.

KW - Conjunctive flow model

KW - Heat diffusion equation

KW - Sub-surface flow

KW - Surface flow

KW - Transient flow

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

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

U2 - 10.1002/(SICI)1097-0363(20000430)32:8<921::AID-FLD993>3.0.CO;2-3

DO - 10.1002/(SICI)1097-0363(20000430)32:8<921::AID-FLD993>3.0.CO;2-3

M3 - Article

AN - SCOPUS:0034177403

VL - 32

SP - 921

EP - 957

JO - International Journal for Numerical Methods in Fluids

JF - International Journal for Numerical Methods in Fluids

SN - 0271-2091

IS - 8

ER -