Finite elements for Matérn-type random fields: Uncertainty in computational mechanics and design optimization

Abstract

This work highlights an approach for incorporating realistic uncertainties into scientific computing workflows based on finite elements, focusing on applications in computational mechanics and design optimization. We leverage Matérn-type Gaussian random fields (GRFs) generated using the SPDE method to model aleatoric uncertainties, including environmental influences, variating material properties, and geometric ambiguities. Our focus lies on delivering practical GRF realizations that accurately capture imperfections and variations and understanding how they impact the predictions of computational models and the topology of optimized designs. We describe a numerical algorithm based on solving a generalized SPDE to sample GRFs on arbitrary meshed domains. The algorithm leverages established techniques and integrates seamlessly with the open-source finite element library MFEM and associated scientific computing workflows, like those found in industrial and national laboratory settings. Our solver scales efficiently for large-scale problems and supports various domain types, including surfaces and embedded manifolds. We showcase its versatility through biomechanics and topology optimization applications. The flexibility and efficiency of SPDE-based GRF generation empower us to run large-scale optimization problems on 2D and 3D domains, including finding optimized designs on embedded surfaces, and to generate topologies beyond the reach of conventional techniques. Moreover, these capabilities allow us to model geometric uncertainties of reconstructed submanifolds, such as the surfaces of cerebral aneurysms. In addition to offering benefits in these specific domains, the proposed techniques transcend specific applications and generalize to arbitrary forward and backward problems in uncertainty quantification involving finite elements.

arXiv comment38 pages, 21 figures

References in corpus (20)

• An extension problem related to the fractional Laplacian
• MFEM: a modular finite element methods library
• The AAA algorithm for rational approximation
• Geometrical structure of Laplacian eigenfunctions
• Numerical Methods for Fractional Diffusion
• Isotropic Gaussian random fields on the sphere: Regularity, fast simulation and stochastic partial differential equations
• Numerical solving the boundary value problem for fractional powers of elliptic operators
• The rational SPDE approach for Gaussian random fields with general smoothness
• Fast simulation of Gaussian random fields
• Tensor FEM for spectral fractional diffusion
• $hp$-Finite Elements for Fractional Diffusion
• High-performance finite elements with MFEM
• A fractional PDE model for turbulent velocity fields near solid walls
• Adaptive sampling strategies for risk-averse stochastic optimization with constraints
• Statistically equivalent surrogate material models and the impact of random imperfections on elasto-plastic response
• A Block Circulant Embedding Method for Simulation of Stationary Gaussian Random Fields on Block-regular Grids
• Learning the structure of wind: A data-driven nonlocal turbulence model for the atmospheric boundary layer
• Surface finite element approximation of spherical Whittle--Matérn Gaussian random fields
• Numerical simulation of endovascular treatment options for cerebral aneurysms
• AdAdaGrad: Adaptive Batch Size Schemes for Adaptive Gradient Methods

Cited by in corpus (1)

• Numerical simulation of dilute polymeric fluids with memory effects in the turbulent flow regime