COMPEST, a PEST-COMSOL interface for inverse multiphysics modelling: Development and application to isotopic fractionation of groundwater contaminants. (May 2019)

• Record Type: Journal Article
• Authors: Halloran, Landon J.S.; Brunner, Philip; Hunkeler, Daniel
• Is Part Of: Computers & geosciences. Volume 126(2019)
• Page Start: 107
• Page End: 119
• Publication Date: 2019-05
• Subjects: Inverse problems -- Parameter estimation -- Uncertainty analysis -- Numerical environmental models -- Hydrogeology -- Contaminant hydrology; Environmental policy -- Periodicals; 550.5
• DOI: 10.1016/j.cageo.2019.02.001
• Languages: English
• ISSNs: 0098-3004

Abstract

Abstract: In the geosciences, inverse problems, wherein observations corresponding to model outputs are known and model parameters are unknown, are commonplace. Many of these problems involve coupled physical, chemical, and other processes that can be modelled using forward finite-element models. Here, we present a novel interface, COMPEST, that connects the parameter estimation and uncertainty analysis package, PEST, with the finite-element modelling package, COMSOL Multiphysics. To demonstrate some of the capabilities of this approach, we also develop and present a novel model for the degradation and transport of chlorohydrocarbons in low-permeability units. This model integrates isotopic fractionation arising from degradation and diffusion. Three implementations of this model with increasing complexity are used to demonstrate the functionality of the developed interface. This linkage provides a means for parameter estimation, uncertainty analysis, and singular value decomposition to gain insight into the behaviour, identifiability, and interdependence of the various parameters in the model. COMPEST is written so as to be suited to a wide range of scientific and engineering applications and thus can be used to link any COMSOL model with PEST. This enables the use of advanced inverse modelling techniques previously unavailable to COMSOL users. Highlights: COMPEST, an interface connecting PEST and COMSOL Multiphysics, is presented. Forward finite-element models are used to solve inverse modelling problems. Novel isotope fractionation models for groundwater contaminants are developed. Parameter estimation and SVD are demonstrated with the developed models.