-
Two Nitsche-based mixed finite element discretizations for the seepage problem in Richards' equation
Authors:
Federico Gatti,
Andrea Bressan,
Alessio Fumagalli,
Domenico Gallipoli,
Leonardo Maria Lalicata,
Simone Pittaluga,
Lorenzo Tamellini
Abstract:
This paper proposes two algorithms to impose seepage boundary conditions in the context of Richards' equation for groundwater flows in unsaturated media. Seepage conditions are non-linear boundary conditions, that can be formulated as a set of unilateral constraints on both the pressure head and the water flux at the ground surface, together with a complementarity condition: these conditions in pr…
▽ More
This paper proposes two algorithms to impose seepage boundary conditions in the context of Richards' equation for groundwater flows in unsaturated media. Seepage conditions are non-linear boundary conditions, that can be formulated as a set of unilateral constraints on both the pressure head and the water flux at the ground surface, together with a complementarity condition: these conditions in practice require switching between Neumann and Dirichlet boundary conditions on unknown portions on the boundary. Upon realizing the similarities of these conditions with unilateral contact problems in mechanics, we take inspiration from that literature to propose two approaches: the first method relies on a strongly consistent penalization term, whereas the second one is obtained by an hybridization approach, in which the value of the pressure on the surface is treated as a separate set of unknowns. The flow problem is discretized in mixed form with div-conforming elements so that the water mass is preserved. Numerical experiments show the validity of the proposed strategy in handling the seepage boundary conditions on geometries with increasing complexity.
△ Less
Submitted 10 July, 2024;
originally announced July 2024.
-
Lowering the Entry Bar to HPC-Scale Uncertainty Quantification
Authors:
Linus Seelinger,
Anne Reinarz,
Jean Benezech,
Mikkel Bue Lykkegaard,
Lorenzo Tamellini,
Robert Scheichl
Abstract:
Treating uncertainties in models is essential in many fields of science and engineering. Uncertainty quantification (UQ) on complex and computationally costly numerical models necessitates a combination of efficient model solvers, advanced UQ methods and HPC-scale resources. The resulting technical complexities as well as lack of separation of concerns between UQ and model experts is holding back…
▽ More
Treating uncertainties in models is essential in many fields of science and engineering. Uncertainty quantification (UQ) on complex and computationally costly numerical models necessitates a combination of efficient model solvers, advanced UQ methods and HPC-scale resources. The resulting technical complexities as well as lack of separation of concerns between UQ and model experts is holding back many interesting UQ applications.
The aim of this paper is to close the gap between advanced UQ methods and advanced models by removing the hurdle of complex software stack integration, which in turn will offer a straightforward way to scale even prototype-grade UQ applications to high-performance resources.
We achieve this goal by introducing a parallel software architecture based on UM-Bridge, a universal interface for linking UQ and models. We present three realistic applications from different areas of science and engineering, scaling from single machines to large clusters on the Google Cloud Platform.
△ Less
Submitted 27 April, 2023;
originally announced April 2023.
-
The Sparse Grids Matlab kit -- a Matlab implementation of sparse grids for high-dimensional function approximation and uncertainty quantification
Authors:
Chiara Piazzola,
Lorenzo Tamellini
Abstract:
The Sparse Grids Matlab Kit provides a Matlab implementation of sparse grids, and can be used for approximating high-dimensional functions and, in particular, for surrogate-model-based uncertainty quantification. It is lightweight, high-level and easy to use, good for quick prototyping and teaching; however, it is equipped with some features that allow its use also in realistic applications. The g…
▽ More
The Sparse Grids Matlab Kit provides a Matlab implementation of sparse grids, and can be used for approximating high-dimensional functions and, in particular, for surrogate-model-based uncertainty quantification. It is lightweight, high-level and easy to use, good for quick prototyping and teaching; however, it is equipped with some features that allow its use also in realistic applications. The goal of this paper is to provide an overview of the data structure and of the mathematical aspects forming the basis of the software, as well as comparing the current release of our package to similar available software.
△ Less
Submitted 9 October, 2023; v1 submitted 17 March, 2022;
originally announced March 2022.
-
Comparing Multi-Index Stochastic Collocation and Multi-Fidelity Stochastic Radial Basis Functions for Forward Uncertainty Quantification of Ship Resistance
Authors:
Chiara Piazzola,
Lorenzo Tamellini,
Riccardo Pellegrini,
Riccardo Broglia,
Andrea Serani,
Matteo Diez
Abstract:
This paper presents a comparison of two multi-fidelity methods for the forward uncertainty quantification of a naval engineering problem. Specifically, we consider the problem of quantifying the uncertainty of the hydrodynamic resistance of a roll-on/roll-off passengers ferry advancing in calm water and subject to two operational uncertainties (ship speed and payload). The first four statistical m…
▽ More
This paper presents a comparison of two multi-fidelity methods for the forward uncertainty quantification of a naval engineering problem. Specifically, we consider the problem of quantifying the uncertainty of the hydrodynamic resistance of a roll-on/roll-off passengers ferry advancing in calm water and subject to two operational uncertainties (ship speed and payload). The first four statistical moments (mean, variance, skewness, kurtosis), and the probability density function for such quantity of interest (QoI) are computed with two multi-fidelity methods, i.e., the Multi-Index Stochastic Collocation (MISC) method and an adaptive multi-fidelity Stochastic Radial Basis Functions (SRBF) algorithm. The QoI is evaluated via computational fluid dynamics simulations, which are performed with the in-house unsteady Reynolds-Averaged Navier-Stokes (RANS) multi-grid solver $χ$navis. The different fidelities employed by both methods are obtained by stopping the RANS solver at different grid levels of the multi-grid cycle. The performance of both methods are presented and discussed: in a nutshell, the findings suggest that, at least for the current implementations of both algorithms, MISC could be preferred whenever a limited computational budget is available, whereas for a larger computational budget SRBFs seem to be preferable, thanks to its robustness to the numerical noise in the evaluations of the QoI.
△ Less
Submitted 26 November, 2021; v1 submitted 1 June, 2021;
originally announced June 2021.
-
On the convergence of adaptive stochastic collocation for elliptic partial differential equations with affine diffusion
Authors:
Martin Eigel,
Oliver Ernst,
Björn Sprungk,
Lorenzo Tamellini
Abstract:
Convergence of an adaptive collocation method for the stationary parametric diffusion equation with finite-dimensional affine coefficient is shown. The adaptive algorithm relies on a recently introduced residual-based reliable a posteriori error estimator. For the convergence proof, a strategy recently used for a stochastic Galerkin method with an hierarchical error estimator is transferred to the…
▽ More
Convergence of an adaptive collocation method for the stationary parametric diffusion equation with finite-dimensional affine coefficient is shown. The adaptive algorithm relies on a recently introduced residual-based reliable a posteriori error estimator. For the convergence proof, a strategy recently used for a stochastic Galerkin method with an hierarchical error estimator is transferred to the collocation setting. Extensions to other variants of adaptive collocation methods (including the classical one proposed in the paper "Dimension-adaptive tensor-product quadratuture" Computing (2003) by T. Gerstner and M. Griebel) is explored.
△ Less
Submitted 16 June, 2021; v1 submitted 17 August, 2020;
originally announced August 2020.
-
Uncertainty Quantification of Ship Resistance via Multi-Index Stochastic Collocation and Radial Basis Function Surrogates: A Comparison
Authors:
Chiara Piazzola,
Lorenzo Tamellini,
Riccardo Pellegrini,
Riccardo Broglia,
Andrea Serani,
Matteo Diez
Abstract:
This paper presents a comparison of two methods for the forward uncertainty quantification (UQ) of complex industrial problems. Specifically, the performance of Multi-Index Stochastic Collocation (MISC) and adaptive multi-fidelity Stochastic Radial Basis Functions (SRBF) surrogates is assessed for the UQ of a roll-on/roll-off passengers ferry advancing in calm water and subject to two operational…
▽ More
This paper presents a comparison of two methods for the forward uncertainty quantification (UQ) of complex industrial problems. Specifically, the performance of Multi-Index Stochastic Collocation (MISC) and adaptive multi-fidelity Stochastic Radial Basis Functions (SRBF) surrogates is assessed for the UQ of a roll-on/roll-off passengers ferry advancing in calm water and subject to two operational uncertainties, namely the ship speed and draught. The estimation of expected value, standard deviation, and probability density function of the (model-scale) resistance is presented and discussed; the required simulations are obtained by the in-house unsteady multi-grid Reynolds Averaged Navier-Stokes (RANS) solver $χ$navis. Both MISC and SRBF use as multi-fidelity levels the evaluations on the different grid levels intrinsically employed by the RANS solver for multi-grid acceleration; four grid levels are used here, obtained as isotropic coarsening of the initial finest mesh. The results suggest that MISC could be preferred when only limited data sets are available. For larger data sets both MISC and SRBF represent a valid option, with a slight preference for SRBF, due to its robustness to noise.
△ Less
Submitted 4 November, 2020; v1 submitted 15 May, 2020;
originally announced May 2020.
-
Compressive Isogeometric Analysis
Authors:
Simone Brugiapaglia,
Lorenzo Tamellini,
Mattia Tani
Abstract:
This work is motivated by the difficulty in assembling the Galerkin matrix when solving Partial Differential Equations (PDEs) with Isogeometric Analysis (IGA) using B-splines of moderate-to-high polynomial degree. To mitigate this problem, we propose a novel methodology named CossIGA (COmpreSSive IsoGeometric Analysis), which combines the IGA principle with CORSING, a recently introduced sparse re…
▽ More
This work is motivated by the difficulty in assembling the Galerkin matrix when solving Partial Differential Equations (PDEs) with Isogeometric Analysis (IGA) using B-splines of moderate-to-high polynomial degree. To mitigate this problem, we propose a novel methodology named CossIGA (COmpreSSive IsoGeometric Analysis), which combines the IGA principle with CORSING, a recently introduced sparse recovery approach for PDEs based on compressive sensing. CossIGA assembles only a small portion of a suitable IGA Petrov-Galerkin discretization and is effective whenever the PDE solution is sufficiently sparse or compressible, i.e., when most of its coefficients are zero or negligible. The sparsity of the solution is promoted by employing a multilevel dictionary of B-splines as opposed to a basis. Thanks to sparsity and the fact that only a fraction of the full discretization matrix is assembled, the proposed technique has the potential to lead to significant computational savings. We show the effectiveness of CossIGA for the solution of the 2D and 3D Poisson equation over nontrivial geometries by means of an extensive numerical investigation.
△ Less
Submitted 28 October, 2020; v1 submitted 13 March, 2020;
originally announced March 2020.
-
Parametric shape optimization for combined additive-subtractive manufacturing
Authors:
Christian Altenhofen,
Marco Attene,
Oliver Barrowclough,
Michele Chiumenti,
Marco Livesu,
Federico Marini,
Massimiliano Martinelli,
Vibeke Skytt,
Lorenzo Tamellini
Abstract:
In the industrial practice, additive manufacturing processes are often followed by post-processing operations such as subtractive machining, milling, etc. to achieve the desired surface quality and dimensional accuracy. Hence, a given part must be 3D printed with extra material to enable such finishing phase. This combined additive/subtractive technique can be optimized to reduce manufacturing cos…
▽ More
In the industrial practice, additive manufacturing processes are often followed by post-processing operations such as subtractive machining, milling, etc. to achieve the desired surface quality and dimensional accuracy. Hence, a given part must be 3D printed with extra material to enable such finishing phase. This combined additive/subtractive technique can be optimized to reduce manufacturing costs by saving printing time and reducing material and energy usage. In this work, a numerical methodology based on parametric shape optimization is proposed for optimizing the thickness of the extra material, allowing for minimal machining operations while ensuring the finishing requirements. Moreover, the proposed approach is complemented by a novel algorithm for generating inner structures leading to reduced distortion and improved weight reduction. The computational effort induced by classical constrained optimization methods is alleviated by replacing both the objective and constraint functions by their sparse-grid surrogates. Numerical results showcase the effectiveness of the proposed approach.
△ Less
Submitted 3 July, 2019; v1 submitted 26 June, 2019;
originally announced July 2019.
-
On Expansions and Nodes for Sparse Grid Collocation of Lognormal Elliptic PDEs
Authors:
Oliver G. Ernst,
Björn Sprungk,
Lorenzo Tamellini
Abstract:
This work is a follow-up to our previous contribution ("Convergence of sparse collocation for functions of countably many Gaussian random variables (with application to elliptic PDEs)", SIAM J. Numer. Anal., 2018), and contains further insights on some aspects of the solution of elliptic PDEs with lognormal diffusion coefficients using sparse grids. Specifically, we first focus on the choice of un…
▽ More
This work is a follow-up to our previous contribution ("Convergence of sparse collocation for functions of countably many Gaussian random variables (with application to elliptic PDEs)", SIAM J. Numer. Anal., 2018), and contains further insights on some aspects of the solution of elliptic PDEs with lognormal diffusion coefficients using sparse grids. Specifically, we first focus on the choice of univariate interpolation rules, advocating the use of Gaussian Leja points as introduced by Narayan and Jakeman ("Adaptive Leja sparse grid constructions for stochastic collocation and high-dimensional approximation", SIAM J. Sci. Comput., 2014) and then discuss the possible computational advantages of replacing the standard Karhunen-Loève expansion of the diffusion coefficient with the Lévy-Ciesielski expansion, motivated by theoretical work of Bachmayr, Cohen, DeVore, and Migliorati ("Sparse polynomial approximation of parametric elliptic PDEs. part II: lognormal coefficients", ESAIM: M2AN, 2016). Our numerical results indicate that, for the problem under consideration, Gaussian Leja collocation points outperform Gauss-Hermite and Genz-Keister nodes for the sparse grid approximation and that the Karhunen-Loève expansion of the log diffusion coefficient is more appropriate than its Lévy-Ciesielski expansion for purpose of sparse grid collocation.
△ Less
Submitted 9 April, 2020; v1 submitted 4 June, 2019;
originally announced June 2019.
-
IGA-based Multi-Index Stochastic Collocation for random PDEs on arbitrary domains
Authors:
Joakim Beck,
Lorenzo Tamellini,
Raúl Tempone
Abstract:
This paper proposes an extension of the Multi-Index Stochastic Collocation (MISC) method for forward uncertainty quantification (UQ) problems in computational domains of shape other than a square or cube, by exploiting isogeometric analysis (IGA) techniques. Introducing IGA solvers to the MISC algorithm is very natural since they are tensor-based PDE solvers, which are precisely what is required b…
▽ More
This paper proposes an extension of the Multi-Index Stochastic Collocation (MISC) method for forward uncertainty quantification (UQ) problems in computational domains of shape other than a square or cube, by exploiting isogeometric analysis (IGA) techniques. Introducing IGA solvers to the MISC algorithm is very natural since they are tensor-based PDE solvers, which are precisely what is required by the MISC machinery. Moreover, the combination-technique formulation of MISC allows the straight-forward reuse of existing implementations of IGA solvers. We present numerical results to showcase the effectiveness of the proposed approach.
△ Less
Submitted 18 March, 2019; v1 submitted 3 October, 2018;
originally announced October 2018.
-
A sparse-grid isogeometric solver
Authors:
Joakim Beck,
Giancarlo Sangalli,
Lorenzo Tamellini
Abstract:
Isogeometric Analysis (IGA) typically adopts tensor-product splines and NURBS as a basis for the approximation of the solution of PDEs. In this work, we investigate to which extent IGA solvers can benefit from the so-called sparse-grids construction in its combination technique form, which was first introduced in the early 90s in the context of the approximation of high-dimensional PDEs. The tests…
▽ More
Isogeometric Analysis (IGA) typically adopts tensor-product splines and NURBS as a basis for the approximation of the solution of PDEs. In this work, we investigate to which extent IGA solvers can benefit from the so-called sparse-grids construction in its combination technique form, which was first introduced in the early 90s in the context of the approximation of high-dimensional PDEs. The tests that we report show that, in accordance to the literature, a sparse-grid construction can indeed be useful if the solution of the PDE at hand is sufficiently smooth. Sparse grids can also be useful in the case of non-smooth solutions when some a-priori knowledge on the location of the singularities of the solution can be exploited to devise suitable non-equispaced meshes. Finally, we remark that sparse grids can be seen as a simple way to parallelize pre-existing serial IGA solvers in a straightforward fashion, which can be beneficial in many practical situations.
△ Less
Submitted 17 February, 2018; v1 submitted 30 July, 2017;
originally announced July 2017.
-
Uncertainty Quantification of geochemical and mechanical compaction in layered sedimentary basins
Authors:
Ivo Colombo,
Fabio Nobile,
Giovanni Porta,
Anna Scotti,
Lorenzo Tamellini
Abstract:
In this work we propose an Uncertainty Quantification methodology for sedimentary basins evolution under mechanical and geochemical compaction processes, which we model as a coupled, time-dependent, non-linear, monodimensional (depth-only) system of PDEs with uncertain parameters. While in previous works (Formaggia et al. 2013, Porta et al., 2014) we assumed a simplified depositional history with…
▽ More
In this work we propose an Uncertainty Quantification methodology for sedimentary basins evolution under mechanical and geochemical compaction processes, which we model as a coupled, time-dependent, non-linear, monodimensional (depth-only) system of PDEs with uncertain parameters. While in previous works (Formaggia et al. 2013, Porta et al., 2014) we assumed a simplified depositional history with only one material, in this work we consider multi-layered basins, in which each layer is characterized by a different material, and hence by different properties. This setting requires several improvements with respect to our earlier works, both concerning the deterministic solver and the stochastic discretization. On the deterministic side, we replace the previous fixed-point iterative solver with a more efficient Newton solver at each step of the time-discretization. On the stochastic side, the multi-layered structure gives rise to discontinuities in the dependence of the state variables on the uncertain parameters, that need an appropriate treatment for surrogate modeling techniques, such as sparse grids, to be effective. We propose an innovative methodology to this end which relies on a change of coordinate system to align the discontinuities of the target function within the random parameter space. The reference coordinate system is built upon exploiting physical features of the problem at hand. We employ the locations of material interfaces, which display a smooth dependence on the random parameters and are therefore amenable to sparse grid polynomial approximations. We showcase the capabilities of our numerical methodologies through two synthetic test cases. In particular, we show that our methodology reproduces with high accuracy multi-modal probability density functions displayed by target state variables (e.g., porosity).
△ Less
Submitted 25 July, 2017; v1 submitted 10 March, 2017;
originally announced March 2017.
-
Convergence of Sparse Collocation for Functions of Countably Many Gaussian Random Variables (with Application to Elliptic PDEs)
Authors:
Oliver G. Ernst,
Björn Sprungk,
Lorenzo Tamellini
Abstract:
We give a convergence proof for the approximation by sparse collocation of Hilbert-space-valued functions depending on countably many Gaussian random variables. Such functions appear as solutions of elliptic PDEs with lognormal diffusion coefficients. We outline a general $L^2$-convergence theory based on previous work by Bachmayr et al. (2016) and Chen (2016) and establish an algebraic convergenc…
▽ More
We give a convergence proof for the approximation by sparse collocation of Hilbert-space-valued functions depending on countably many Gaussian random variables. Such functions appear as solutions of elliptic PDEs with lognormal diffusion coefficients. We outline a general $L^2$-convergence theory based on previous work by Bachmayr et al. (2016) and Chen (2016) and establish an algebraic convergence rate for sufficiently smooth functions assuming a mild growth bound for the univariate hierarchical surpluses of the interpolation scheme applied to Hermite polynomials. We verify specifically for Gauss-Hermite nodes that this assumption holds and also show algebraic convergence w.r.t. the resulting number of sparse grid points for this case. Numerical experiments illustrate the dimension-independent convergence rate.
△ Less
Submitted 28 March, 2017; v1 submitted 22 November, 2016;
originally announced November 2016.
-
Optimal-order isogeometric collocation at Galerkin superconvergent points
Authors:
Monica Montardini,
Giancarlo Sangalli,
Lorenzo Tamellini
Abstract:
In this paper we investigate numerically the order of convergence of an isogeometric collocation method that builds upon the least-squares collocation method presented in [1] and the variational collocation method presented in [2]. The focus is on smoothest B-splines/NURBS approximations, i.e, having global $C^{p-1}$ continuity for polynomial degree $p$. Within the framework of [2], we select as c…
▽ More
In this paper we investigate numerically the order of convergence of an isogeometric collocation method that builds upon the least-squares collocation method presented in [1] and the variational collocation method presented in [2]. The focus is on smoothest B-splines/NURBS approximations, i.e, having global $C^{p-1}$ continuity for polynomial degree $p$. Within the framework of [2], we select as collocation points a subset of those considered in [1], which are related to the Galerkin superconvergence theory. With our choice, that features local symmetry of the collocation stencil, we improve the convergence behaviour with respect to [2], achieving optimal $L^2$-convergence for odd degree B-splines/NURBS approximations. The same optimal order of convergence is seen in [1], where, however a least-squares formulation is adopted. Further careful study is needed, since the robustness of the method and its mathematical foundation are still unclear.
△ Less
Submitted 23 September, 2016; v1 submitted 7 September, 2016;
originally announced September 2016.
-
Multi-index Stochastic Collocation convergence rates for random PDEs with parametric regularity
Authors:
Abdul-Lateef Haji-Ali,
Fabio Nobile,
Lorenzo Tamellini,
Raul Tempone
Abstract:
We analyze the recent Multi-index Stochastic Collocation (MISC) method for computing statistics of the solution of a partial differential equation (PDEs) with random data, where the random coefficient is parametrized by means of a countable sequence of terms in a suitable expansion. MISC is a combination technique based on mixed differences of spatial approximations and quadratures over the space…
▽ More
We analyze the recent Multi-index Stochastic Collocation (MISC) method for computing statistics of the solution of a partial differential equation (PDEs) with random data, where the random coefficient is parametrized by means of a countable sequence of terms in a suitable expansion. MISC is a combination technique based on mixed differences of spatial approximations and quadratures over the space of random data and, naturally, the error analysis uses the joint regularity of the solution with respect to both the variables in the physical domain and parametric variables. In MISC, the number of problem solutions performed at each discretization level is not determined by balancing the spatial and stochastic components of the error, but rather by suitably extending the knapsack-problem approach employed in the construction of the quasi-optimal sparse-grids and Multi-index Monte Carlo methods. We use a greedy optimization procedure to select the most effective mixed differences to include in the MISC estimator. We apply our theoretical estimates to a linear elliptic PDEs in which the log-diffusion coefficient is modeled as a random field, with a covariance similar to a Matérn model, whose realizations have spatial regularity determined by a scalar parameter. We conduct a complexity analysis based on a summability argument showing algebraic rates of convergence with respect to the overall computational work. The rate of convergence depends on the smoothness parameter, the physical dimensionality and the efficiency of the linear solver. Numerical experiments show the effectiveness of MISC in this infinite-dimensional setting compared with the Multi-index Monte Carlo method and compare the convergence rate against the rates predicted in our theoretical analysis.
△ Less
Submitted 21 July, 2016; v1 submitted 17 November, 2015;
originally announced November 2015.
-
Multi-Index Stochastic Collocation for random PDEs
Authors:
Abdul-Lateef Haji-Ali,
Fabio Nobile,
Lorenzo Tamellini,
Raul Tempone
Abstract:
In this work we introduce the Multi-Index Stochastic Collocation method (MISC) for computing statistics of the solution of a PDE with random data. MISC is a combination technique based on mixed differences of spatial approximations and quadratures over the space of random data. We propose an optimization procedure to select the most effective mixed differences to include in the MISC estimator: suc…
▽ More
In this work we introduce the Multi-Index Stochastic Collocation method (MISC) for computing statistics of the solution of a PDE with random data. MISC is a combination technique based on mixed differences of spatial approximations and quadratures over the space of random data. We propose an optimization procedure to select the most effective mixed differences to include in the MISC estimator: such optimization is a crucial step and allows us to build a method that, provided with sufficient solution regularity, is potentially more effective than other multi-level collocation methods already available in literature. We then provide a complexity analysis that assumes decay rates of product type for such mixed differences, showing that in the optimal case the convergence rate of MISC is only dictated by the convergence of the deterministic solver applied to a one dimensional problem. We show the effectiveness of MISC with some computational tests, comparing it with other related methods available in the literature, such as the Multi-Index and Multilevel Monte Carlo, Multilevel Stochastic Collocation, Quasi Optimal Stochastic Collocation and Sparse Composite Collocation methods.
△ Less
Submitted 21 July, 2016; v1 submitted 29 August, 2015;
originally announced August 2015.
