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 1238 products
-
03/2017 - 01/13/2017
Fumagalli, I.; Parolini, N.; Verani, M.
On a free-surface problem with moving contact line: from variational principles to stable numerical approximations | Abstract | | We analyze a free-surface problem described by time-dependent Navier-Stokes equations. Surface tension, capillary effects and wall friction are taken into account in the evolution of the system, influencing the motion of the contact line - where the free surface hits the wall - and of the dynamics of the contact angle. The differential equations governing the phenomenon are first derived from the variational principle of minimum reduced dissipation, and then discretized by means of the ALE approach. The numerical properties of the resulting scheme are investigated, drawing a parallel with the physical properties holding at the continuous level. Some instability issues are addressed in detail, in the case of an explicit treatment of the geometry, and novel additional terms are introduced in the discrete formulation in order to damp the instabilities. Numerical tests assess the suitability of the approach, the influence of the parameters, and the effectiveness of the new stabilizing terms. |
-
02/2017 - 01/13/2017
Arena, M.; Calissano, A.; Vantini, S.
Monitoring Rare Categories in Sentiment and Opinion Analysis - Expo Milano 2015 on Twitter Platform. | Abstract | | This paper proposes a new aggregated classification scheme aimed
to support the implementation of text analysis methods in contexts
characterised by the presence of rare text categories. The proposed
approach starts from the aggregate supervised text classifier developed
by Hopkins and King and moves forward relying on rare event sampling
methods. In details, it enables the analyst to enlarge the number
of text categories whose proportions can be estimated preserving the
estimation accuracy of standard aggregate supervised algorithms and
reducing the working time w.r.t. to unconditionally increase the size
of the random training set. The approach is applied to study the daily
evolution of the web reputation of Expo Milano 2015, before, during
and after the event. The data set is constituted by about 900,000
tweets in Italian and 260,000 tweets in English, posted about the event
between March 2015 and December 2015. The analysis provides an
interesting portray of the evolution of Expo stakeholders’ opinions over
time and allow to identify the main drivers of Expo reputation. The
algorithm will be implemented as a running option of the next release
of R package ReadMe. |
-
01/2017 - 01/04/2017
Riccobelli, D.; Ciarletta, P.
Rayleigh-Taylor instability in soft elastic layers | Abstract | | This work investigates the morphological stability of a soft body composed of two heavy elastic layers, attached to a rigid surface and subjected only to the bulk gravity force. Using theoretical and computational tools, we characterize the selection of different patterns as well as their nonlinear evolution, unveiling the interplay between elastic and geometric effects for their formation.
Unlike similar gravity-induced shape transitions in fluids, as the Rayleigh-Taylor instability, we prove that the nonlinear elastic effects saturate the dynamic instability of the bifurcated solutions, displaying a rich morphological diagram where both digitations and stable wrinkling can emerge. The results of this work provide important guidelines for the design of novel soft systems with tunable shapes, with several applications in engineering sciences. |
-
58/2016 - 12/27/2016
Antonietti, P. F.; Bruggi, M. ; Scacchi, S.; Verani, M.
On the Virtual Element Method for Topology Optimization on polygonal meshes: a numerical study | Abstract | | It is well known that the solution of topology optimization problems may be affected both by the geometric properties of the computational mesh, which can steer the minimization process towards local (and non-physical) minima, and by the accuracy of the method employed to discretize the underlying differential problem, which may not be able to correctly capture the physics of the problem. In light of the above remarks, in this paper we consider polygonal meshes and employ the virtual element method (VEM) to solve two classes of paradigmatic topology optimization problems, one governed by nearly-incompressible and compressible linear elasticity and the other by Stokes equations. Several numerical results show the virtues of our polygonal VEM based approach with respect to more standard methods. |
-
57/2016 - 12/24/2016
Bassi, C.; Abbà, A.; Bonaventura, L.; Valdettaro, L.
Large Eddy Simulation of gravity currents with a high order DG method | Abstract | | This work deals with Direct Numerical Simulations (DNS) and Large Eddy Simulations (LES) of turbulent gravity currents, performed by means of a Discontinuous Galerkin (DG) Finite Element method.
In particular, a DG-LES approach in which the filter operation is built in the numerical discretization has been employed, similarly to VMS approaches.
Numerical simulations of non-Boussinesq lock-exchange benchmark problems show that, in the DNS case, the proposed method allows to correctly reproduce relevant features of variable density flows with gravity. Moreover LES results highlight the excessively high dissipation of the Smagorinsky model with respect to the Germano dynamic procedure, providing a first indication of the superiority of dynamic models in the context of gravity currents. |
-
56/2016 - 12/24/2016
Guerciotti, B.; Vergara, C.; Ippolito, S.; Quarteroni, A.; Antona, C.; Scrofani, R.
A computational fluid-structure interaction analysis of coronary Y-grafts | Abstract | | Coronary artery disease is one of the leading causes of death worldwide. The stenotic coronary vessels are generally treated with coronary artery bypass grafts (CABG), which can be either arterial (internal mammary
artery, radial artery) or venous (saphenous vein). However, the different
mechanical properties of the graft can influence the outcome of the procedure in terms of risk of restenosis and subsequent graft failure. In this paper, we perform a computational fluid-structure interaction (FSI) analysis of patient-specific multiple CABGs (Y-grafts) with the aim of better understanding the influence of the choice of bypass (arterial vs venous) on the risk of graft failure. Our results show that the use of a venous bypass results in a more disturbed flow field at the anastomosis and in higher stresses in the vessel wall with respect to the arterial one. This could explain the better long-term patency of the arterial bypasses experienced in the clinical practice. |
-
55/2016 - 12/15/2016
Antonietti, P. F.; Facciola' C.; Russo A.; Verani M.;
Discontinuous Galerkin approximation of flows in fractured porous media on polytopic grids | Abstract | | We present a numerical approximation of Darcy's flow through
a fractured porous medium which employs discontinuous Galerkin
methods on polytopic grids. For simplicity, we analyze the case of a
single fracture represented by a (d-1)-dimensional interface between
two d-dimensional subdomains, d = 2; 3. We propose a discontinuous
Galerkin finite element approximation for the
flow in the porous matrix which is coupled with a conforming finite element scheme for the flow in the fracture. Suitable (physically consistent) coupling conditions complete the model. We theoretically analyse the resulting formulation, prove its well-posedness, and derive optimal a priori error estimates in a suitable (mesh-dependent) energy norm. Two-dimensional numerical experiments are reported to assess the theoretical results. |
-
54/2016 - 12/13/2016
Vergara, C.; Le Van, D.; Quadrio, M.; Formaggia, L.; Domanin, M.
Large Eddy Simulations of blood dynamics in abdominal aortic aneurysms | Abstract | | In this work we address the study of transition to turbulence effects
in abdominal aortic aneurysms (AAA). The formation of transitional ef-
fects in such districts is caused by the heart pulsatility and the sudden
change of diameter, and has been recorded by means of clinical measures
and computational studies. Here we propose, for the first time, the use
of a large eddy simulation (LES) model to describe transition to turbu-
lence in realistic scenarios of AAA obtained from radiological images. To
this aim, we post-process the obtained numerical solutions to produce sig-
nificant quantities, such as the ensemble-averaged velocity and wall shear
stress, the standard deviation of the velocity field, the Q-criterion. The
results show the suitability of the considered LES model and the presence
of significant transitional effects during the mid-deceleration phase around the impingement region.
|
|
