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 1251 products
-
19/2008 - 09/03/2008
Badia, Santiago; Nobile, Fabio; Vergara, Christian
Robin-Robin preconditioned Krylov methods for fluid-structure interaction problems | Abstract | | In this work we propose a Robin-Robin preconditioner combined
with Krylov iterations for the solution of the interface system
arising in fluid-structure interaction (FSI) problems. It can be
seen as a partitioned FSI procedure and in this respect it
generalizes the ideas introduced in [{ sl Badia, Nobile and Vergara,
J. Comput. Phys. { bf{227}} (2008) 7027 --7051]}.
We analyze the convergence of GMRES iterations with the Robin-Robin
preconditioner on a model problem and compare its efficiency with some
existing algorithms. The method is shown to be very efficient for many
challenging fluid-structure interaction problems, such as those
characterized by a large added-mass effect or by enclosed fluids. In
particular, the possibility to solve balloon-type problems without any
special treatment makes this algorithm very appealing compared to the
computationally intensive existing approaches |
-
18/2008 - 07/24/2008
Bonaventura, Luca; Castruccio, Stefano; Crippa, Paola; Lonati, Giovanni
Geostatistical estimate of PM10 concentrations in Northern Italy: validation of kriging reconstructions with classical and flexible variogram models | Abstract | | The applicability of classical geostatistical tools to the reconstruction
of PM10 concentration fields over the entire Po Valley has been assessed,
based on a large dataset of daily PM10 data spanning the period 2003-
2006. The impact of data detrending by the median polish procedure and
of the variogram model chosen for the geostatistical estimates have been
investigated, by comparison of the results obtained with several isotropic
variogram models as well as with anisotropic
exible variogram models.
The relative merits of the different approaches were evaluated by crossvalidating
the resulting reconstructions and performing normality tests on
the corresponding residuals. Although exponential and linear variograms
yield reliable reconstructions in most of the cases, the analysis has highlighted
significant seasonal and interannual variations in the basic features
of the estimated concentration fields and residual correlation structure. As
a consequence, none of the classical models is able to cope with all the
different situations encountered, while the anisotropic
exible variogram
models appear to provide a more robust tool for automatic reconstruction
of the PM10 concentration fields without expert user intervention. |
-
17/2008 - 07/17/2008
Ern, Alexandre; Perotto, Simona; Veneziani, Alessandro
Hierarchical model reduction for advection-diffusion-reaction problems | Abstract | | Some engineering problems ranging from blood flow to river flow, from
internal combustion engines to electronic devices have been recently modelled by coupling problems with different space dimensions (geometrical
multiscale method). In this paper we focus on a new approch, where different levels of detail of the problem at hand stem from a different selection of the dimension of a suitable function space. The coarse and fine models
are thus identified in a straightforward way. Moreover this approach lends
itself to an automatic model adaptive strategy. The approach is addressed
on a 2D linear advection-diffusion reaction problem |
-
16/2008 - 07/11/2008
Formaggia, Luca; Miglio, Edie; Mola, Andrea; Scotti, Anna
Numerical simulation of the dynamics of boats by a variational inequality approach | Abstract | | In this paper we present some recent numerical studies on fluidstructure
interaction problems in the presence of free surface flow. We
consider the dynamics of a rowing boat, simulated as a rigid body. We
focus on an approach based on formulating the floating body problem
as an inequality constraint on the water elevation. A splitting
procedure is used to develop an efficient numerical scheme where the
inequality constraint is imposed only on a wave like equation representing
an hydrostatic approximation of the hydrodynamic equations.
Numerical tests demostrate the effectiveness of the proposed procedure |
-
15/2008 - 07/02/2008
Micheletti, Stefano; Perotto, Simona
An anisotropic mesh adaptation procedure for an optimal control problem of the advection-diffusion-reaction equation | Abstract | | We derive an anisotropic a posteriori error estimator for a PDE-con-
strained optimal control problem, governed by the scalar advection-diffusion-
reaction equation. With a view to the advection dominated case, a strongly
consistent symmetric stabilization is employed so that the “optimize-then-
discretize” and “discretize-then-optimize” philosophies coincide and lead
to the same discrete problem. The estimator is turned into an anisotropic
mesh adaptation procedure which allows us to approximate the cost func-
tional within to a given tolerance. Both an academic and a realistic test,
inspired by an environmental application, assess the performance of the
proposed approach. |
-
14/2008 - 07/01/2008
D'Angelo, Carlo; Zunino, Paolo
A finite element method based on weighted interior penalties for heterogeneous incompressible flows | Abstract | | We propose a finite element scheme for the approximation of multidomain
heterogeneous problems arising in the general framework of linear
incompressible flows (e.g. Stokes’ and Darcy’s equations). We exploit stabilized
mixed finite elements together with Nitsche type matching conditions
that automatically adapt to the coupling of different subproblem combinations.
Optimal error estimates are derived for the coupled problem. Finally,
we propose and analyze an iterative splitting strategy for the approximation
of the multidomain solution by means of a sequence of independent
and local subproblems. Thanks to the introduction of a suitable relaxation
strategy, the iterative method turns out to be convergent for any possible
coupling between subproblems |
-
13/2008 - 06/30/2008
Sangalli, Laura Maria; Secchi, Piercesare; Vantini, Simone; Vitelli, Valeria
K-means alignment for curve clustering | Abstract | | We deal with the problem of curve clustering when curves are mis- aligned. We propose a k-means alignment algorithm which jointly cluster and align the curves. We illustrate the procedure via simulation studies and applications to real data.
|
-
12/2008 - 06/26/2008
Passerini, Tiziano; de Luca, Maria Rita; Formaggia, Luca; Quarteroni, Alfio; Veneziani, Alessandro
A 3D/1D geometrical multiscale model of cerebral vasculature | Abstract | | Geometrical multiscale modeling is a strategy advocated in computational
hemodynamics for representing in a single numerical model dynamics that
involve different space scales. This approach is particularly useful to describe complex
networks such as the circle of Willis in the cerebral vasculature. In this paper
we present a multiscale model of the cerebral circulation where a one dimensional
description of the circle of Willis, relying on the one-dimensional Euler equations, is
coupled to a fully three dimensional model of a carotid artery, based on the solution
of the incompressible Navier-Stokes equations. Even if vascular compliance is often
not relevant to the meaningfulness of 3D results, it is crucial in the multiscale model,
since it is the driving mechanism of pressure wave propagation. Unfortunately, 3D
simulations in compliant domains still demand computational costs significantly
higher than the rigid case. Appropriate matching conditions between the two models
have been devised to concentrate the effects of the compliance at the interfaces and
to obtain reliable results still solving a 3D problems on rigid vessels |
|