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 1151 prodotti
-
23/2008 - 28/10/2008
Hysing, Shu-Ren; Turek, Stefan; Kuzmin, Dmitri; Parolini, Nicola; Burman, Erik; Ganesan, Sashikumaar; Tobiska, Lutz
Quantitative benchmark computations of two-dimensional bubble dynamics | Abstract | | Benchmark configurations for quantitative validation and comparison of incompressible interfacial flow
codes, which model two-dimensional bubbles rising in liquid columns, are proposed. The benchmark
quantities: circularity, center of mass, and mean rise velocity are defined and measured to monitor
convergence towards a reference solution. Comprehensive studies are undertaken by three independent
research groups, two representing Eulerian level set finite element codes, and one representing an ALE
moving grid approach.
The first benchmark test case considers a bubble with small density and viscosity ratios which
undergoes moderate shape deformation. The results from all codes agree very well allowing for target
reference values to be established. For the second test case, a bubble with a very low density compared
to that of the surrounding fluid, the results for all groups are in good agreement up to the point of
break up, after which all three codes predict different bubble shapes. This highlights the need for the
research community to invest more effort in obtaining reference solutions to problems involving break
up and coalescence.
Other research groups are encouraged to participate in these benchmarks by contacting the authors
and submitting their own data. The reference data for the computed benchmark quantities can also
be supplied for validation purposes |
-
22/2008 - 23/10/2008
Nobile, Fabio; Tempone, Raul
Analysis and implementation issues for the numerical approximation of parabolic equations with random coefficients | Abstract | | We consider the problem of numerically approximating statistical moments
of the solution of a time dependent linear parabolic partial differential
equation (PDE), whose coefficients and/or forcing terms are spatially
correlated random fields. The stochastic coefficients of the PDE are approximated
by truncated Karhunen-Lo eve expansions driven by a finite number
of uncorrelated random variables. After approximating the stochastic coefficients
the original stochastic PDE turns into a new deterministic parametric
PDE of the same type, the dimension of the parameter set being
equal to the number of random variables introduced.
After proving that the solution of the parametric PDE problem is analytic
with respect to the parameters, we consider global polynomial approximations
based on tensor product, total degree or sparse polynomial spaces
and constructed by either a Stochastic Galerkin or a Stochastic Collocation
approach. We derive convergence rates for the different cases and present
numerical results that show how these approaches are a valid alternative to
the more traditional Monte Carlo Method for this class of problems |
-
21/2008 - 29/09/2008
Antonietti, Paola; Suli, Endre
Domain Decomposition Preconditioning for Discontinuous Galerkin Approximation of Convection-Diffusion Problems | Abstract | | We study a class of nonoverlapping Schwarz preconditioners for DG ap- proximations of convection-di usion equations. The generalized minimal residual (GMRES) Krylov space-based iterative solver is accelerated with the proposed preconditioners. We discuss the issue of convergence of the re- sulting preconditioned iterative method, and demonstrate through numer- ical computations that the classical Schwarz convergence theory for non- symmetric and indefinite problems developed by Cai and Widlund [SIAM J. Sci. Statist. Comput. 13 (1992) 243--258], [SIAM J. Numer. Anal.
30 (1993) 936--952] cannot be applied to explain theoretically the converge observed numerically |
-
20/2008 - 24/09/2008
David, Filippo; Micheletti, Stefano; Perotto, Simona
Model adaptation enriched with an anisotropic mesh spacing for advection-diffusion-reaction systems | Abstract | | We propose a procedure aiming at reducing the computational cost involved
in the numerical approximation of (possibly nonlinear) advectiondiffusion-
reaction systems. The idea is to suitably combine a model with a
mesh adaptive procedure. In particular we first derive, separately, a model
error estimator and an anisotropic estimator for the discretization error,
suited for driving a model and a mesh adaptivity algorithm, respectively.
These two strategies are then properly combined, allowing for a merged
model-mesh control. The whole procedure is finally assessed on some numerical
test cases, essentially inspired by ecological and environmental application |
-
19/2008 - 03/09/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 - 24/07/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 - 17/07/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 - 11/07/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 |
|
