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 1314 prodotti
-
34/2016 - 06/10/2016
Menafoglio, A.; Secchi, P.
Statistical analysis of complex and spatially dependent data: a review of Object Oriented Spatial Statistics | Abstract | | We review recent advances in Object Oriented Spatial Statistics, a system of ideas, algorithms and methods that allows the analysis of high dimensional and complex data when their spatial dependence is an important issue. At the intersection of different disciplines -- including mathematics, statistics, computer science and engineering -- Object Oriented Spatial Statistics provides the right perspective to address key problems in varied contexts, from Earth and life sciences to urban planning. We illustrate a few paradigmatic methods applied to problems of prediction, classification and smoothing, giving emphasis to the key ideas Object Oriented Spatial Statistics relies upon. |
-
32/2016 - 30/09/2016
Tarabelloni, N.; Schenone, E.; Collin, A.; Ieva, F.; Paganoni, A.M.; Gerbeau, J.-F.
Statistical Assessment and Calibration of Numerical ECG Models | Abstract | | Objective: Because of the inter-subject variability of ECGs in a healthy
population, it is not straightforward to assess the quality of synthetic ECGs
produced by deterministic mathematical models. We propose a statistical method
to address this question.
Methods: We use a dataset of 1588 healthy, real ECGs and we introduce a way to
calibrate the deterministic model so that its output fits the dataset. Our
approach is based on the concepts of spatial quantiles and spatial depths. These
notions are convenient to manipulate functional data since they provide a
non-parametric way to measure the discrepancy of the model output with a
distribution of data.
Results: The method is successfully applied to two very different models: a
phenomenological model based on ordinary differential equations, and a complex
biophysical model based on partial differential equations set on a
three-dimensional geometry of the heart and the torso. We show in particular
that the proposed calibration strategy allows us to improve the quality of the
ECG obtained with the biophysical model.
Significance: The proposed methodology is to our knowledge the first attempt to
assess the quality of synthetic ECGs with quantitative statistical arguments.
More generally it can be applied to other situations where a deterministic model
produces a functional output that has to be compared with a population of
measurements containing inter-subject variability. |
-
31/2016 - 09/09/2016
Antonietti, P.F.; Merlet, B.; Morgan, P.; Verani, M.
Convergence to equilibrium for a second-order time semi-discretization of the Cahn-Hilliard equation | Abstract | | We consider a second-order two-step time semi-discretization of the Cahn-Hilliard equation with an analytic nonlinearity. The time-step is chosen small enough so that the pseudo-energy associated with the discretization is non-increasing at every time iteration. We prove that the sequence generated by the scheme converges to a steady state as time tends to infinity. We also obtain convergence rates in the energy norm. The proof is based on the Lojasiewicz-Simon inequality. |
-
30/2016 - 09/09/2016
Abramowicz, K.; Häger, C.; Pini, A.; Schelin, L.; Sjöstedt de Luna, S.; Vantini, S.
Nonparametric inference for functional-on-scalar linear models applied to knee kinematic hop data after injury of the anterior cruciate ligament | Abstract | | Motivated by the analysis of the dependence of knee movement patterns during functional tasks on subject-speci c covariates, we introduce a distribution-free procedure for testing a functional-on-scalar linear model with fixed e ects. The procedure does not only test the global hypothesis on all the domain, but also selects the intervals where statistically signi ficant e ects are detected. We prove that the proposed tests are provided with an asymptotic control of the interval-wise error rate, i.e., the
probability of falsely rejecting any interval of true null hypotheses. The procedure is applied to one-leg hop data from a study on anterior cruciate ligament injury. We compare knee kinematics of three groups of individuals (two injured groups with di erent treatments, and one group of healthy controls), taking individual-speci c covariates into account. |
-
29/2016 - 30/07/2016
Miglio, E.; Parolini, N.; Penati, M.; Porcù, R.
GPU parallelization of brownout simulations with a non-interacting particles dynamic model | Abstract | | The term brownout refers to the uplift of sand particles in the air and is
generated when a helicopter is close to a dusty soil. When a brownout occurs the visibility area is remarkably restricted, thus the pilot may be disoriented and the helicopter may dangerously collide with the ground. Simulations of a brownout require tens of millions of particles in order to be significative, so that the execution of a serial program takes a very long time. In order to speedup the computation, the GPU-parallelization of a brownout simulation program is performed in order to obtain a notable speedup. The dynamics of the particles are considered in a Lagrangian way, under the effect of the gravity force and of a precomputed aerodynamic field. The particles are independent from each other since collisions between them are not taken into account. Thus trajectories are independent and the parallelization is very effective. In this paper we discuss in detail the impact of the techniques used for the GPU implementation of the parallel code on the performance. |
-
28/2016 - 30/07/2016
Antonietti, P.F.; Dal Santo, N.; Mazzieri, I.; Quarteroni, A.
A high-order discontinuous Galerkin approximation to ordinary differential equations with applications to elastodynamics | Abstract | | The aim of this work is to propose and analyze a new high order discontinuous Galerkin finite element method for the time integration of a Cauchy problem second order ordinary differential equations. These equations typically arise after space semi-discretization of second order hyperbolic-type differential problems, e.g., wave, elastodynamics and acoustics equation. After introducing the new method, we analyze its well-posedness and prove a-priori error estimates in a suitable (mesh-dependent) norm. Numerical results are also presented to verify the theoretical estimates. space-time finite elements, discontinuous Galerkin methods, second order hyperbolic equations. |
-
27/2016 - 25/07/2016
Repossi, E.; Rosso, R.; Verani, M.
A phase-field model for liquid-gas mixtures: mathematical modelling and Discontinuous Galerkin discretization | Abstract | | In this article we propose a phase-field approach to model a liquid-gas mixture that might also provide a description of the expansion stage of a metal foam inside a hollow mold. We conceive the mixture as a two-phase incompressible-compressible fluid governed by a Navier-Stokes-Cahn-Hilliard system of equations, and we adapt the Lowengrub-Truskinowsky model to take into account the expansion of the gaseous phase. The resulting system of equations is characterized by a velocity field that fails to be divergence-free, by a logarithmic term for the pressure that enters the Gibbs free-energy expression and by the viscosity that degenerates in the gas phase. In the second part of the article we propose an energy-based numerical scheme that, at the discrete level, preserves the mass conservation property and the energy dissipation law of the original system. We use a Discontinuous Galerkin approximation for the spatial approximation and a modified midpoint based scheme for the time approximation. |
-
26/2016 - 25/07/2016
Brunetto, D.; Calderoni, F.; Piccardi, C.
Communities in criminal networks: A case study | Abstract | | Criminal organizations tend to be clustered to reduce risks of detection and information leaks. Yet, the literature exploring the relevance of subgroups for their internal structure is so far very limited. The paper applies methods of community analysis to explore the structure of a criminal network representing the individuals' co-participation in meetings. It draws from a case study on a large law enforcement operation (``Operazione Infinito'') tackling the 'Ndrangheta, a mafia organization from Calabria, a southern Italian region. The results show that the network is indeed clustered and that communities are associated, in a non trivial way, with the internal organization of the 'Ndrangheta into different ``locali'' (similar to mafia families). Furthermore, the results of community analysis can improve the prediction of the ``locale'' membership of the criminals (up to two thirds of any random sample of nodes) and the leadership roles (above 90% precision in classifying nodes as either bosses or non-bosses). The implications of these findings on the interpretation of the structure and functioning of the criminal network are discussed. |
|
