MOX Reports
The preprint collection of the Laboratory for Modeling and Scientific Computation MOX. It mainly contains works on numerical
analysis and mathematical modeling applied to engineering problems. MOX web site is mox.polimi.it
Found 1152 products
-
25/2010 - 07/15/2010
Pena, G.; Prud'homme, C.; Quarteroni, A.
High Order Methods for the Approximation of the Incompressible Navier-Stokes Equations in a Moving Domain | Abstract | | In this paper we address the numerical approximation of the incompressible Navier-Stokes
equations in a moving domain by the spectral element method and high order time integrators. We present the Arbitrary Lagrangian Eulerian (ALE) formulation of the incompressible
Navier-Stokes equations and propose a numerical method based on the following kernels: a Lagrange basis associated with Fekete points in the spectral element method context, BDF time integrators, an ALE map of high degree, and a robust algebraic linear solver.
In particular, the high degree ALE map is appropriate to deal with a computational domain whose boundary is described with curved elements. Finally, we apply the proposed strategy to a
test case. |
-
24/2010 - 07/14/2010
Tamellini, Lorenzo; Formaggia, Luca; Miglio, Edie; Scotti, Anna
An Uzawa iterative scheme for the simulation of floating boats | Abstract | | The numerical simulation of a floating boat requires to solve a complex interaction between free surface flow and the dynamics of the boat. In this work we focus on solution schemes with an explicit representation of the fluid surface and on the specic problem of the imposition of the presence of a floating object. The floating object represents an unilateral constraint on
the fluid surface, which cannot rise above a the level dened by the floating object itself.
We propose an Uzawa type iterative scheme for this type of situations and we describe an efficient implementation in the context of a
particular finite element discretization and algebraic splitting procedure. |
-
23/2010 - 07/13/2010
Baeck, Joakim; Nobile, Fabio; Tamellini, Lorenzo; Tempone, Raul
Stochastic Spectral Galerkin and collocation methods for PDEs with random coefficients: a numerical comparison | Abstract | | Much attention has recently been devoted to the development of Stochastic Galerkin (SG) and Stochastic Collocation (SC) methods for uncertainty quantification. An open and relevant research topic is the comparison of these two methods. By introducing a suitable generalization of the classical sparse grid SC method, we are able to compare SG and SC on the same underlying multivariate polynomial space in terms of accuracy versus computational work.
The approximation spaces considered here include
isotropic and anisotropic versions of Tensor Product (TP), Total Degree(TD), Hyperbolic Cross (HC) and Smolyak (SM) polynomials.
Numerical results for linear elliptic SPDEs indicate a slight computational work advantage
of isotropic SC over SG, with SC-SM and SG-TD being the best choices of approximation spaces for each method. Finally, numerical results
corroborate the optimality of the theoretical estimate of anisotropy ratios introduced by the authors in a previous work for the construction of
anisotropic approximation spaces. |
-
22/2010 - 07/12/2010
D'Angelo, Carlo; Zunino, Paolo
Numerical approximation with Nitsche s coupling of transient Stokes /Darcy s flow problems applied to hemodynamics | Abstract | | In this work, we consider a time dependent coupled Stokes-Darcy flow problem and study an approximation method based on a unified finite element scheme complemented with implicit time stepping. Our finite element formulation relies on a weighing strategy in which the physical and discretization parameters are taken into account to robustly enforce interface and boundary conditions by means of the Nitsche s method. We prove unconditional absolute stability and optimal convergence of the scheme, and discuss the
algebraic properties of the associated discrete problem. Finally, we present numerical experiments confirming the predicted convergence behavior and algebraic properties, and report an application to the computational analysis of
blood flow and plasma filtration in arteries after the implantation of a vascular
graft. |
-
21/2010 - 06/30/2010
Grieco, Niccolò; Ieva, Francesca; Paganoni, Anna Maria
Provider Profiling Using Mixed Effects Models on a Case Study concerning STEMI Patients | Abstract | | Provider profiling is the evaluation process of the performance of hospitals, doctors, and other medical practitioners to enhance the quality of
medical care. In this work we describe statistical analyses conducted on MOMI2 (MOnth MOnitoring Myocardial Infarction in MIlan) survey, a collection of data concerning patients admitted with STEMI (ST-Elevation Myocardial Infarction) diagnosis in one of the hospitals belonging to the Milan urban area Network, with the aim of pointing out process indicators
to be used in health-care evaluation. An ffective statistical support to decisional process for clinical and organizational governance is then obtained analyzing and modelling data coming from clinical registers and administrative data banks. Typically this kind of data are grouped at first level by structure where patients are admitted to, so multilevel models have been considered and tted to catch and explain overdispersion phenomena.
Moreover, we compare three different techniques of hospitals classication based respectively on traditional survival rates comparison, on analysis of variance components in fitted generalized linear mixed eects models and
non parametric random eects estimation.
Key words:
Generalized Linear Mixed Models
In-hospital survival
ST-Elevation Myocardial Infarction
Provider Profiling |
-
20/2010 - 06/09/2010
Nobile, Fabio; Quarteroni, Alfio; Ruiz Baier, Ricardo
Numerical solution of an active strain formulation for the electromechanical activity in the heart | Abstract | | We propose a finite element method for solving a system of equations describing the coupling
between cardiac mechanics and electrical signalling. The model is based on a multiplicative
decomposition of the deformation tensor into a passive and active part, the latter carrying
the information of the electrical potential propagation and anisotropy of the cardiac tissue into the equations of incompressible or nearly incompressible nonlinear elasticity, governing the
mechanical response of the biological material. Moreover, by changing from an Eulerian to
a Lagrangian configuration, the underlying problem exhibits a nonlinear diffusion term in the
equations of the electrical propagation (namely, the bidomain and monodomain equations).
Piecewise quadratic finite elements are used to approximate the displacements field, while for
pressure, electrical potentials and ionic variables, we use piecewise linear elements. Various test cases show that the proposed method is able to capture some important features of the
studied phenomenon, and illustrate the behavior of the global model. |
-
19/2010 - 06/08/2010
Gaudio, Loredana; Quarteroni, Alfio
hN-adaptive spectral element discretization of optimal control problems for environmental applications | Abstract | | In this work, we propose an hN adaptive algorithm for optimal control problems whose space discretization is based on the spectral element method.
We focus on a special class of problems described by quadratic cost functional and linear advection-diffusion state equation. This model has been
employed to forecast and suitably reduce, by an optimization process, the pollutant concentration in air on a specic area of interest. Moreover, a
local adaptive algorithm is employed to achieve a desired accuracy on the cost functional by reducing the computational cost of the problem. |
-
18/2010 - 06/07/2010
Antonietti, P.F.; Fadel, N.A.; Verani, M.
Modelling and numerical simulation of the polymeric extrusion process in textile products | Abstract | | We present the mathematical modelling and the numerical simulation of the polymer extrusion process of hollow and non-hollow yarns. In particular, the polymeric material flow has been modelled as non-Newtonian, incompressible, isothermal and steady. We attempt to develop a possible strategy for effective die design in profile extrusion and we numerically asses the validity of the model on both academic and real industrial test cases. Keywords:Free boundary flow, Die-swell, Hollow yarn. 2000 Mathematics Subject Classification. 76D07, 35Q30. |
|