Low-order finite elements for parameter identification in groundwater flow

Jiannan Xiang, Derek Elsworth

Research output: Contribution to journalArticle

4 Citations (Scopus)

Abstract

Improved evaluation and characterization of the geologic subsurface are necessary precursors to enhancing our ability to recover desirable resources and guard against undesirable contaminants. Mathematical models are increasingly used to assess the effects of these activities when the fidelity of material parameters describing the system control the reliability of the resulting prediction. Inverse models provide a formal means of space, or history, matching observed data to determine the unknown spatial distribution of the required parameters. The finite element method is one of the most popular numerical methods and is used in this study for transmissivity identification in steady groundwater flow. To optimize parameter identification, a comparison of different elements is carried out. It is well known that high-order elements usually result in improved accuracy in forward solution of engineering problems. However, this fact might not be true in the inverse solution. The study in this paper uses different element orders for Taylor's series analysis and in coding a finite element program. The analysis indicates that strong variation of transmissivity results in a larger residual error for high-order elements than for low-order elements. An example with an analytical solution is used for numerical comparison. The computed results show that low-order elements, rather than high-order elements, yield better results in parameter estimation. When the variation of unknown parameters is large, the error within the high-order elements is large. It is recommended that low-order finite elements or constant elements be used for an inverse solution. The comparison also illustrates that any minimization procedure may be used to minimize the residual error and to limit numerical difficulties in inverse solution. A two-dimensional example shows that the relative errors are very small when a constant element is used but that errors increase as the measured head distribution flattens.

Original languageEnglish (US)
Pages (from-to)256-266
Number of pages11
JournalApplied Mathematical Modelling
Volume15
Issue number5
DOIs
StatePublished - Jan 1 1991

Fingerprint

Higher-order Elements
Groundwater Flow
Groundwater flow
Parameter Identification
Element Order
Identification (control systems)
Finite Element
Flatten
History Matching
Inverse Model
Taylor series
Numerical Comparisons
Relative Error
Steady Flow
Spatial Distribution
Precursor
Unknown Parameters
Fidelity
Parameter Estimation
Analytical Solution

All Science Journal Classification (ASJC) codes

  • Modeling and Simulation
  • Applied Mathematics

Cite this

@article{cdc9e3ca041645c9b463a40db09ce516,
title = "Low-order finite elements for parameter identification in groundwater flow",
abstract = "Improved evaluation and characterization of the geologic subsurface are necessary precursors to enhancing our ability to recover desirable resources and guard against undesirable contaminants. Mathematical models are increasingly used to assess the effects of these activities when the fidelity of material parameters describing the system control the reliability of the resulting prediction. Inverse models provide a formal means of space, or history, matching observed data to determine the unknown spatial distribution of the required parameters. The finite element method is one of the most popular numerical methods and is used in this study for transmissivity identification in steady groundwater flow. To optimize parameter identification, a comparison of different elements is carried out. It is well known that high-order elements usually result in improved accuracy in forward solution of engineering problems. However, this fact might not be true in the inverse solution. The study in this paper uses different element orders for Taylor's series analysis and in coding a finite element program. The analysis indicates that strong variation of transmissivity results in a larger residual error for high-order elements than for low-order elements. An example with an analytical solution is used for numerical comparison. The computed results show that low-order elements, rather than high-order elements, yield better results in parameter estimation. When the variation of unknown parameters is large, the error within the high-order elements is large. It is recommended that low-order finite elements or constant elements be used for an inverse solution. The comparison also illustrates that any minimization procedure may be used to minimize the residual error and to limit numerical difficulties in inverse solution. A two-dimensional example shows that the relative errors are very small when a constant element is used but that errors increase as the measured head distribution flattens.",
author = "Jiannan Xiang and Derek Elsworth",
year = "1991",
month = "1",
day = "1",
doi = "10.1016/0307-904X(91)90003-8",
language = "English (US)",
volume = "15",
pages = "256--266",
journal = "Applied Mathematical Modelling",
issn = "0307-904X",
publisher = "Elsevier Inc.",
number = "5",

}

Low-order finite elements for parameter identification in groundwater flow. / Xiang, Jiannan; Elsworth, Derek.

In: Applied Mathematical Modelling, Vol. 15, No. 5, 01.01.1991, p. 256-266.

Research output: Contribution to journalArticle

TY - JOUR

T1 - Low-order finite elements for parameter identification in groundwater flow

AU - Xiang, Jiannan

AU - Elsworth, Derek

PY - 1991/1/1

Y1 - 1991/1/1

N2 - Improved evaluation and characterization of the geologic subsurface are necessary precursors to enhancing our ability to recover desirable resources and guard against undesirable contaminants. Mathematical models are increasingly used to assess the effects of these activities when the fidelity of material parameters describing the system control the reliability of the resulting prediction. Inverse models provide a formal means of space, or history, matching observed data to determine the unknown spatial distribution of the required parameters. The finite element method is one of the most popular numerical methods and is used in this study for transmissivity identification in steady groundwater flow. To optimize parameter identification, a comparison of different elements is carried out. It is well known that high-order elements usually result in improved accuracy in forward solution of engineering problems. However, this fact might not be true in the inverse solution. The study in this paper uses different element orders for Taylor's series analysis and in coding a finite element program. The analysis indicates that strong variation of transmissivity results in a larger residual error for high-order elements than for low-order elements. An example with an analytical solution is used for numerical comparison. The computed results show that low-order elements, rather than high-order elements, yield better results in parameter estimation. When the variation of unknown parameters is large, the error within the high-order elements is large. It is recommended that low-order finite elements or constant elements be used for an inverse solution. The comparison also illustrates that any minimization procedure may be used to minimize the residual error and to limit numerical difficulties in inverse solution. A two-dimensional example shows that the relative errors are very small when a constant element is used but that errors increase as the measured head distribution flattens.

AB - Improved evaluation and characterization of the geologic subsurface are necessary precursors to enhancing our ability to recover desirable resources and guard against undesirable contaminants. Mathematical models are increasingly used to assess the effects of these activities when the fidelity of material parameters describing the system control the reliability of the resulting prediction. Inverse models provide a formal means of space, or history, matching observed data to determine the unknown spatial distribution of the required parameters. The finite element method is one of the most popular numerical methods and is used in this study for transmissivity identification in steady groundwater flow. To optimize parameter identification, a comparison of different elements is carried out. It is well known that high-order elements usually result in improved accuracy in forward solution of engineering problems. However, this fact might not be true in the inverse solution. The study in this paper uses different element orders for Taylor's series analysis and in coding a finite element program. The analysis indicates that strong variation of transmissivity results in a larger residual error for high-order elements than for low-order elements. An example with an analytical solution is used for numerical comparison. The computed results show that low-order elements, rather than high-order elements, yield better results in parameter estimation. When the variation of unknown parameters is large, the error within the high-order elements is large. It is recommended that low-order finite elements or constant elements be used for an inverse solution. The comparison also illustrates that any minimization procedure may be used to minimize the residual error and to limit numerical difficulties in inverse solution. A two-dimensional example shows that the relative errors are very small when a constant element is used but that errors increase as the measured head distribution flattens.

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

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

U2 - 10.1016/0307-904X(91)90003-8

DO - 10.1016/0307-904X(91)90003-8

M3 - Article

VL - 15

SP - 256

EP - 266

JO - Applied Mathematical Modelling

JF - Applied Mathematical Modelling

SN - 0307-904X

IS - 5

ER -