Quaderni MOX
Pubblicazioni
del Laboratorio di Modellistica e Calcolo Scientifico MOX. I lavori riguardano prevalentemente il campo dell'analisi numerica, della statistica e della modellistica matematica applicata a problemi di interesse ingegneristico. Il sito del Laboratorio MOX è raggiungibile
all'indirizzo mox.polimi.it
Trovati 1288 prodotti
-
39/2012 - 16/09/2012
Ieva, F.; Paganoni, A.M.; Ziller, S.
Operational risk management: a statistical perspective | Abstract | | This work presents a statistical model for operational risk management. We distinguish different types of operational Event, we model the probability of event occurrence (the frequency distribution) and the economic impact of the single event (the severity distribution), and then the aggregated distribution is obtained through convolution of frequency and severity, for each event type. The main problem is the parameters estimation of the severity distribution above a suitable threshold, that we consider as an unknown parameter to be estimated as well. An application to a case study is also presented.
|
-
38/2012 - 15/09/2012
Antonietti, P.F.; Bigoni, N.; Verani, M.
Mimetic finite difference approximation of quasilinear elliptic problems | Abstract | | In this work we approximate the solution of a quasilinear elliptic problem of monotone type by using the Mimetic Finite Difference (MFD) method. Under a suitable approximation assumption, we prove that the MFD approximate solution converges, with optimal rate, to the exact solution in a mesh-dependent energy norm. The resulting nonlinear discrete problem is then solved iteratively via linearization by applying the Kacanov method. The convergence of the Kacanov algorithm in the discrete mimetic framework is also proved. Several numerical experiments confirm the theoretical analysis. |
-
37/2012 - 14/09/2012
Nobile, F.; Pozzoli, M.; Vergara, C.
Exact and inexact partitioned algorithms for fluid-structure interaction problems with finite elasticity in haemodynamics | Abstract | | In this paper we consider the numerical solution of the three-dimensional
(3D) fluid-structure interaction problem in haemodynamics, in the case of
physiological geometries and data, and finite elasticity vessel deformations.
We introduce new partitioned algorithms and compare their efficiency with
that of existing ones. We also study some new inexact variants, obtained
from semi-implicit approximations, and show that they allow to improve the
efficiency while preserving the accuracy of the related exact (implicit) scheme. |
-
36/2012 - 13/09/2012
Canuto, C.; Verani, M.
On the Numerical Analysis of Adaptive Spectral/hp Methods for Elliptic Problems | Abstract | | We provide an overview of the state of the art of adaptive strategies for high-order $hp$ discretizations of partial differential
equations; at the same time, we draw attention on some recent results of ours concerning the convergence and complexity analysis of adaptive algorithm of spectral and spectral-element type. Complexity is studied under the assumption that the solution belongs to a sparsity class of exponential type, which means that its best $N$-term approximation error in the chosen piecewise polynomial basis decays at an exponential rate with respect to $N$. |
-
35/2012 - 10/09/2012
Pigoli, D.; Aston, J.A.D.; Dryden, I.L.; Secchi, P.
Distances and Inference for Covariance Functions | Abstract | | A framework is developed for inference concerning the covariance operator of a functional random process, where the covariance operator itself is an object of interest for the statistical analysis. Distances for comparing positive definite covariance matrices are either extended or shown to be inapplicable for functional data. In particular, an infinite dimensional analogue of
the Procrustes size and shape distance is developed. The convergence of the finite dimensional approximations to the infinite dimensional distance metrics is also shown. To perform inference, a Fréchet estimator for the average covariance function is introduced, and a
permutation procedure to test the equality of the covariance operator between two groups is then considered. The proposed techniques are applied to two problems where inference concerning the covariance is of interest. Firstly, in data arising from a study into cerebral aneurysms, it is of interest to determine whether two groups of data can be combined when comparing with a third group. For this to be done, it is necessary to assess whether the covariance structures of the two groups are the same or different. Secondly, in a philological study of cross-linguistic dependence, the use of covariance operators has been suggested as a way to incorporate quantitative phonetic information. It is shown
that distances between languages derived from phonetic covariance functions can provide insight into relationships between the Romance languages. |
-
34/2012 - 04/09/2012
Menafoglio, A.; Dalla Rosa, M.; Secchi, P.
A Universal Kriging predictor for spatially dependent functional data of a Hilbert Space | Abstract | | We address the problem of predicting spatially dependent functional data belonging to a Hilbert space, with a Functional Data Analysis approach. Having defined new global measures of spatial variability for functional random processes, we derive a Universal Kriging predictor for functional data.
Consistently with the new established theoretical results, we develop a two-step procedure for predicting georeferenced functional data: first model selection and estimation of the spatial mean (drift), then Universal Kriging prediction on the basis of the identified dichotomy model, sum of deterministic drift and stochastic residuals.
The proposed methodology is tested by means of a simulation study and finally applied to daily mean temperatures curves aiming at reconstructing the space-time field of temperatures of Canada s Maritimes Provinces. |
-
33/2012 - 24/08/2012
Motamed, M.; Nobile, F.; Tempone, R.
Analysis and computation of the elastic wave equation with random coefficients | Abstract | | We analyze the stochastic initial-boundary value problem for the elastic wave equation with random coefficients and deterministic data. We propose a stochastic collocation method for computing statistical moments of the solution or statistics of some given quantities of interest. We study the convergence rate of the error in the stochastic collocation method. In particular, we show that, the rate of convergence depends on the regularity of the solution or the quantity of interest in the stochastic space, which is in turn related to the regularity of the deterministic data in the physical space and the type of the quantity of interest. We demonstrate that a fast rate of convergence is possible in two cases: for the elastic wave solutions with high regular data; and for some high regular quantities of interest even in the presence of low regular data. We perform numerical examples, including a simplified earthquake, which confirm the analysis and show that the collocation method is a valid alternative to the more traditional Monte Carlo sampling method for problems with high stochastic regularity. |
-
32/2012 - 22/08/2012
Formaggia, L.; Fumagalli, A.; Scotti A.; Ruffo, P
A reduced model for Darcy s problem in networks of fractures | Abstract | | Subsurface flows are influenced by the presence of faults and large fractures which act as preferential paths or barriers for the flow. In
literature models were proposed to handle fractures in a porous medium as objects of codimension 1. In this work we consider the case of a network of intersecting fractures, with the aim of deriving physically consistent and
effective interface conditions to impose at the intersection between fractures. This new model accounts for the angle between fractures at the intersections and allows for jumps of pressure across the intersection. This latter property permits to describe more accurately the flow when fractures are characterised by different properties, than other models that impose
pressure continuity. The main mathematical properties of the model, derived in the two-dimensional setting, are analysed. As concerns the numerical discretization we allow the grids of the fractures to be independent, thus in general non-matching at the intersection, by means of the
extended finite element method (XFEM), to increase the flexibility of the method in the case of complex geometries characterized by a high number of fractures. |
|
