147 research outputs found
Regularity and sparse approximation of the recursive first moment equations for the lognormal Darcy problem
We study the Darcy boundary value problem with log-normal permeability field.
We adopt a perturbation approach, expanding the solution in Taylor series
around the nominal value of the coefficient, and approximating the expected
value of the stochastic solution of the PDE by the expected value of its Taylor
polynomial. The recursive deterministic equation satisfied by the expected
value of the Taylor polynomial (first moment equation) is formally derived.
Well-posedness and regularity results for the recursion are proved to hold in
Sobolev space-valued H\"older spaces with mixed regularity. The recursive first
moment equation is then discretized by means of a sparse approximation
technique, and the convergence rates are derived
Moment equations for the mixed formulation of the Hodge Laplacian with stochastic loading term
We study the mixed formulation of the stochastic Hodge-Laplace problem defined on an n-dimensional domain D (n≥1), with random forcing term. In particular, we focus on the magnetostatic problem and on the Darcy problem in the three-dimensional case. We derive and analyse the moment equations, that is, the deterministic equations solved by the mth moment (m≥1) of the unique stochastic solution of the stochastic problem. We find stable tensor product finite element discretizations, both full and sparse, and provide optimal order-of-convergence estimates. In particular, we prove the inf-sup condition for sparse tensor product finite element space
Fast Least-Squares Pad\'e approximation of problems with normal operators and meromorphic structure
In this work, we consider the approximation of Hilbert space-valued
meromorphic functions that arise as solution maps of parametric PDEs whose
operator is the shift of an operator with normal and compact resolvent, e.g.
the Helmholtz equation. In this restrictive setting, we propose a simplified
version of the Least-Squares Pad\'e approximation technique introduced in [6]
following [11]. In particular, the estimation of the poles of the target
function reduces to a low-dimensional eigenproblem for a Gramian matrix,
allowing for a robust and efficient numerical implementation (hence the "fast"
in the name). Moreover, we prove several theoretical results that improve and
extend those in [6], including the exponential decay of the error in the
approximation of the poles, and the convergence in measure of the approximant
to the target function. The latter result extends the classical one for scalar
Pad\'e approximation to our functional framework. We provide numerical results
that confirm the improved accuracy of the proposed method with respect to the
one introduced in [6] for differential operators with normal and compact
resolvent
MATHICSE Technical Report : Perturbation analysis for the Darcy problem with log-normal permeability
We study the single-phase flow in a saturated, bounded heterogeneous porous medium. We model the permeability as a log-normal random field. We perform a perturbation analysis, expanding the solution in Taylor series. The series is directly computable in the case of a random field parametrized by a finite number of random variables. On the other hand, in the case of an infinite dimensional random field, suitable equations satisfied by the derivatives of the stochastic solution can be derived. We give a theoretical upper bound for the norm of the residual of the Taylor expansion which predicts the divergence of the series as the polynomial degree goes to infinity. We provide a formula to compute the optimal degree for the Taylor polynomial and the maximum achievable accuracy of the perturbation approach. Our theoretical findings are confirmed by numerical experiments in the simple case where the permeability field is described using one random variable
Perturbation analysis for the Darcy problem with log-normal permeability
We study the single-phase flow in a saturated, bounded heterogeneous porous medium. We model the permeability as a log-normal random field. We perform a perturbation analysis, expanding the solution in Taylor series. The series is directly computable in the case of a random field parametrized by a finite number of random variables. On the other hand, in the case of an infinite dimensional random field, suitable equations satisfied by the derivatives of the stochastic solution can be derived. We give a theoretical upper bound for the norm of the residual of the Taylor expansion which predicts the divergence of the series as the polynomial degree goes to infinity. We provide a formula to compute the optimal degree for the Taylor polynomial and the maximum achievable accuracy of the perturbation approach. Our theoretical findings are confirmed by numerical experiments in the simple case where the permeability field is described using one random variabl
MATHICSE Technical Report : Moment equations for the mixed formulation of the Hodge Laplacian with stochastic data
We study the mixed formulation of the stochastic Hodge-Laplace problem defined on a n-dimensional domain D (), with random forcing term. In particular, we focus on the magnetostatic problem and on the Darcy problem in the three dimensional case. We derive and analyze the moment equations, that is the deterministic equations solved by the m-th moment () of the unique stochastic solution of the stochastic problem. We find stable tensor product finite element discretizations, both full and sparse, and provide optimal order of convergence estimates. In particular, we prove the inf-sup condition for sparse tensor product finite element spaces
MATHICSE Technical Report :Convergence analysis of Padé approximations for Helmholtz frequency response problems
The present work concerns the approximation of the solution map S associated to the parametric Helmholtz boundary value problem, i.e., the map which associates to each (real) wavenumber belonging to a given interval of interest the corresponding solution of the Helmholtz equation. We introduce a least squares rational Padé-type approximation technique applicable to any meromorphic Hilbert space-valued univariate map, and we prove the uniform convergence of the Padé approximation error on any compact subset of the interval of interest that excludes any pole. This general result is then applied to the Helmholtz solution map S, which is proven to be meromorphic in C, with a pole of order one in every (single or multiple) eigenvalue of the Laplace operator with the considered boundary conditions. Numerical tests are provided that confirm the theoretical upper bound on the Padé approximation error for the Helmholtz solution map
- …