MFEM: A modular finite element methods library
MFEM is an open-source, lightweight, flexible and scalable C++ library for modular finite
element methods that features arbitrary high-order finite element meshes and spaces …
element methods that features arbitrary high-order finite element meshes and spaces …
Analysis of boundary effects on PDE-based sampling of Whittle--Matérn random fields
We consider the generation of samples of a mean-zero Gaussian random field with Matérn
covariance function. Every sample requires the solution of a differential equation with …
covariance function. Every sample requires the solution of a differential equation with …
Multilevel approximation of Gaussian random fields: fast simulation
We propose and analyze several multilevel algorithms for the fast simulation of possibly
nonstationary Gaussian random fields (GRFs) indexed, for example, by the closure of a …
nonstationary Gaussian random fields (GRFs) indexed, for example, by the closure of a …
Efficient white noise sampling and coupling for multilevel Monte Carlo with nonnested meshes
When solving stochastic partial differential equations (SPDEs) driven by additive spatial
white noise, the efficient sampling of white noise realizations can be challenging. Here, we …
white noise, the efficient sampling of white noise realizations can be challenging. Here, we …
Fast sampling of parameterised Gaussian random fields
Gaussian random fields are popular models for spatially varying uncertainties, arising for
instance in geotechnical engineering, hydrology or image processing. A Gaussian random …
instance in geotechnical engineering, hydrology or image processing. A Gaussian random …
Multilevel hierarchical decomposition of finite element white noise with application to multilevel Markov chain Monte Carlo
HR Fairbanks, U Villa, PS Vassilevski - SIAM Journal on Scientific Computing, 2021 - SIAM
In this work we develop a new hierarchical multilevel approach to generate Gaussian
random field realizations in an algorithmically scalable manner that is well suited to …
random field realizations in an algorithmically scalable manner that is well suited to …
[HTML][HTML] Comparison and application of non-conforming mesh models for flow in fractured porous media using dual Lagrange multipliers
Geological settings, such as reservoirs, include fractures with different material properties
and geometric features. Hence, numerical simulations in applied geophysics demands for …
and geometric features. Hence, numerical simulations in applied geophysics demands for …
Accelerating Multilevel Markov Chain Monte Carlo Using Machine Learning Models
S Reddy, H Fairbanks - arXiv preprint arXiv:2405.11179, 2024 - arxiv.org
This work presents an efficient approach for accelerating multilevel Markov Chain Monte
Carlo (MCMC) sampling for large-scale problems using low-fidelity machine learning …
Carlo (MCMC) sampling for large-scale problems using low-fidelity machine learning …
Frontiers in Mortar methods for isogeometric analysis
Complex geometries as common in industrial applications consist of multiple patches, if
spline based parametrizations are used. The requirements for the generation of analysis …
spline based parametrizations are used. The requirements for the generation of analysis …
Estimating posterior quantity of interest expectations in a multilevel scalable framework
HR Fairbanks, S Osborn… - Numerical Linear Algebra …, 2021 - Wiley Online Library
Scalable approaches for uncertainty quantification are necessary for characterizing
prediction confidence in large‐scale subsurface flow simulations with uncertain …
prediction confidence in large‐scale subsurface flow simulations with uncertain …