Hierarchical Approaches for Efficient and Scalable Solution of Inverse Problems Governed by Partial Differential Equations
Inverse problems governed by partial differential equations (PDEs) is a means to learn, from data, unknown or uncertain aspects of PDE-based mathematical models. PDE-based models often are constructed from scientific principles and are pervasive in science and engineering. The improvement of such a...
Main Author: | |
---|---|
Other Authors: | |
Format: | Thesis |
Language: | English |
Published: |
eScholarship, University of California
2022
|
Subjects: | |
Online Access: | https://escholarship.org/uc/item/99x865w6 https://escholarship.org/content/qt99x865w6/qt99x865w6.pdf |
id |
ftcdlib:oai:escholarship.org:ark:/13030/qt99x865w6 |
---|---|
record_format |
openpolar |
spelling |
ftcdlib:oai:escholarship.org:ark:/13030/qt99x865w6 2024-09-15T18:12:34+00:00 Hierarchical Approaches for Efficient and Scalable Solution of Inverse Problems Governed by Partial Differential Equations Hartland, Tucker Andrew Petra, Noemi 2022-01-01 application/pdf https://escholarship.org/uc/item/99x865w6 https://escholarship.org/content/qt99x865w6/qt99x865w6.pdf en eng eScholarship, University of California qt99x865w6 https://escholarship.org/uc/item/99x865w6 https://escholarship.org/content/qt99x865w6/qt99x865w6.pdf public Mathematics Hierarchical Matrices Interior-point Method Inverse Problems Point Spread Function etd 2022 ftcdlib 2024-06-28T06:28:20Z Inverse problems governed by partial differential equations (PDEs) is a means to learn, from data, unknown or uncertain aspects of PDE-based mathematical models. PDE-based models often are constructed from scientific principles and are pervasive in science and engineering. The improvement of such a model can be accomplished by determining a ``best'' set of model parameters, or by reducing the uncertainty of the parameters that characterize the model. Ultimately, the goal is to improve the predictive capacity of such models and quantitatively understand the limitations of model-based predictions.In this dissertation, we describe algorithmic approaches for the Newton-based solution of large-scale computational inverse problems governed by PDEs. We present and motivate hierarchical matrix approximations of the Hessian, a key-component of Newton's method, as a means to exploit localized sensitivities of underlying elliptic PDE operators and which perform well for inverse problems with highly-informative data. To circumvent the computational challenges associated with generating hierarchical matrix approximations of matrix-free operators, such as the Hessian, we describe a local point spread function methodology, whose computational cost is independent of the problem discretization. It is numerically demonstrated that hierarchical matrix approximations of the Hessian can be effective for large-scale ice sheet inverse problems with highly informative data and thin ice sheets. Lastly, we present a scalable means to solve PDE- and bound-constrained optimization problems by an interior-point filter line-search strategy that leverages performant algebraic multigrid linear solvers. Bound constraints arise naturally in many inverse problems as a means to enforce sign-definiteness as dictated by e.g., physical principles. The inclusion of bound constraints does add particular computational challenges such as non-smooth complementarity conditions as part of the conditions for optimality. The Newton linear system solve ... Thesis Ice Sheet University of California: eScholarship |
institution |
Open Polar |
collection |
University of California: eScholarship |
op_collection_id |
ftcdlib |
language |
English |
topic |
Mathematics Hierarchical Matrices Interior-point Method Inverse Problems Point Spread Function |
spellingShingle |
Mathematics Hierarchical Matrices Interior-point Method Inverse Problems Point Spread Function Hartland, Tucker Andrew Hierarchical Approaches for Efficient and Scalable Solution of Inverse Problems Governed by Partial Differential Equations |
topic_facet |
Mathematics Hierarchical Matrices Interior-point Method Inverse Problems Point Spread Function |
description |
Inverse problems governed by partial differential equations (PDEs) is a means to learn, from data, unknown or uncertain aspects of PDE-based mathematical models. PDE-based models often are constructed from scientific principles and are pervasive in science and engineering. The improvement of such a model can be accomplished by determining a ``best'' set of model parameters, or by reducing the uncertainty of the parameters that characterize the model. Ultimately, the goal is to improve the predictive capacity of such models and quantitatively understand the limitations of model-based predictions.In this dissertation, we describe algorithmic approaches for the Newton-based solution of large-scale computational inverse problems governed by PDEs. We present and motivate hierarchical matrix approximations of the Hessian, a key-component of Newton's method, as a means to exploit localized sensitivities of underlying elliptic PDE operators and which perform well for inverse problems with highly-informative data. To circumvent the computational challenges associated with generating hierarchical matrix approximations of matrix-free operators, such as the Hessian, we describe a local point spread function methodology, whose computational cost is independent of the problem discretization. It is numerically demonstrated that hierarchical matrix approximations of the Hessian can be effective for large-scale ice sheet inverse problems with highly informative data and thin ice sheets. Lastly, we present a scalable means to solve PDE- and bound-constrained optimization problems by an interior-point filter line-search strategy that leverages performant algebraic multigrid linear solvers. Bound constraints arise naturally in many inverse problems as a means to enforce sign-definiteness as dictated by e.g., physical principles. The inclusion of bound constraints does add particular computational challenges such as non-smooth complementarity conditions as part of the conditions for optimality. The Newton linear system solve ... |
author2 |
Petra, Noemi |
format |
Thesis |
author |
Hartland, Tucker Andrew |
author_facet |
Hartland, Tucker Andrew |
author_sort |
Hartland, Tucker Andrew |
title |
Hierarchical Approaches for Efficient and Scalable Solution of Inverse Problems Governed by Partial Differential Equations |
title_short |
Hierarchical Approaches for Efficient and Scalable Solution of Inverse Problems Governed by Partial Differential Equations |
title_full |
Hierarchical Approaches for Efficient and Scalable Solution of Inverse Problems Governed by Partial Differential Equations |
title_fullStr |
Hierarchical Approaches for Efficient and Scalable Solution of Inverse Problems Governed by Partial Differential Equations |
title_full_unstemmed |
Hierarchical Approaches for Efficient and Scalable Solution of Inverse Problems Governed by Partial Differential Equations |
title_sort |
hierarchical approaches for efficient and scalable solution of inverse problems governed by partial differential equations |
publisher |
eScholarship, University of California |
publishDate |
2022 |
url |
https://escholarship.org/uc/item/99x865w6 https://escholarship.org/content/qt99x865w6/qt99x865w6.pdf |
genre |
Ice Sheet |
genre_facet |
Ice Sheet |
op_relation |
qt99x865w6 https://escholarship.org/uc/item/99x865w6 https://escholarship.org/content/qt99x865w6/qt99x865w6.pdf |
op_rights |
public |
_version_ |
1810450155616337920 |