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
-
13/2022 - 02/26/2022
Grasselli, M.; Parolini, N.; Poiatti, A.; Verani, M.
Non-isothermal non-Newtonian fluids: the stationary case | Abstract | | The stationary Navier-Stokes equations for a non-Newtonian incompressible fluid are coupled with the stationary heat equation and subject to Dirichlet type boundary conditions. The viscosity is supposed to depend on the temperature and the stress depends on the strain through a suit-able power law depending on p in (1,2) (shear thinning case). For this problem we establish the existence of a weak solution as well as we prove some regularity results both for the Navier-Stokes and the Stokes cases.Then, the latter case with the Carreau power law is approximated through a FEM scheme and some error estimates are obtained. Such estimates are then validated through some two-dimensional numerical experiments. |
-
12/2022 - 02/25/2022
Antonietti, P.F.; Dassi, F.; Manuzzi, E.
Machine Learning based refinement strategies for polyhedral grids with applications to Virtual Element and polyhedral Discontinuous Galerkin methods | Abstract | | We propose two new strategies based on Machine Learning techniques to handle polyhedral grid refinement, to be possibly employed within an adaptive framework. The first one employs the k-means clustering algorithm to partition the points of the polyhedron to be refined. This strategy is a variation of the well known Centroidal Voronoi Tessellation. The second one employs Convolutional Neural Networks to classify the “shape” of an element so that “ad-hoc” refinement criteria can be defined. This strategy can be used to enhance existing refinement strategies, including the k-means strategy, at a low online computational cost. We test the proposed algorithms considering two families of finite element methods that support arbitrarily shaped polyhedral elements, namely the Virtual Element Method (VEM) and the Polygonal Discontinuous Galerkin (PolyDG) method. We demonstrate that these strategies do preserve the structure and the quality of the underlaying grids, reducing the overall computational cost and mesh complexity. |
-
11/2022 - 02/22/2022
Sampaoli, S.; Agosti, A.; Pozzi, G.; Ciarletta,P.
A toy model of misfolded protein aggregation and neural damage propagation in neudegenerative diseases | Abstract | | Neurodegenerative diseases (NDs) result from the transformation and accumulation of misfolded proteins within the nervous system. They have common features, like the chronic nature and the progressive destruction of neurons in specific areas of the brain. Several mathematical models have been proposed to investigate the biological processes underlying NDs, focusing on the kinetics of polymerization and fragmentation at the microscale and on the spread of neural damage at a macroscopic level. The aim of this work is to bridge the gap between microscopic and macroscopic approaches proposing a toy partial differential model able to take into account both for the short-time dynamics of the misfolded proteins aggregating in plaques and the long-term evolution of tissue damage. Using the theoretical framework of mixtures theory, we considered the brain as a biphasic material made of misfolded protein aggregates and of healthy tissue. The resulting Cahn-Hilliard type equation for the misfolded proteins contains a growth term depending on the local availability of precursor proteins, that follow a reaction-diffusion equation. The misfolded proteins also posses a chemotactic mass flux driven by gradients of neural damage, that is caused by local accumulation of misfolded protein and that evolves slowly according to an Allen-Cahn equation. The partial differential model is solved numerically using the finite element method in a simple two-dimensional domain, evaluating the effects of the mobility of the misfolded protein and the diffusion of the neural damage. We considered both isotropic and anisotropic mobility coefficients, highlighting that the spreading front of the neural damage follows the direction of the largest eigenvalue of the mobility tensor. In both cases, we computed two biomarkers for quantifying the aggregation in plaques and the evolution of neural damage, that are in qualitative agreement with the characteristic Jack curves for many NDs. |
-
10/2022 - 02/22/2022
Fresca, S.; Manzoni, A.
Real-time simulation of parameter-dependent fluid flows through deep learning-based reduced order models | Abstract | | Simulating fluid flows in different virtual scenarios is of key importance in engineering applications. However, high-fidelity, full-order models relying, e.g., on the finite element method, are unaffordable whenever fluid flows must be simulated in almost real-time. Reduced order models (ROMs) relying, e.g., on proper orthogonal decomposition (POD) provide reliable approximations to parameter-dependent fluid dynamics problems in rapid times. However, they might require expensive hyper-reduction strategies for handling parameterized nonlinear terms, and enriched reduced spaces (or Petrov-Galerkin projections) if a mixed velocity-pressure formulation is considered, possibly hampering the evaluation of reliable solutions in real-time. Dealing with fluid-structure interactions entails even higher difficulties. The proposed deep learning (DL)-based ROMs overcome all these limitations by learning in a non-intrusive way both the nonlinear trial manifold and the reduced dynamics. To do so, they rely on deep neural networks, after performing a former dimensionality reduction through POD enhancing their training times substantially. The resulting POD-DL ROMs are shown to provide accurate results in almost real-time for the flow around a cylinder benchmark, the fluid-structure interaction between an elastic beam attached to a fixed, rigid block and a laminar incompressible flow, and the blood flow in a cerebral aneurysm. |
-
09/2022 - 02/22/2022
Corti, M.; Zingaro, A.; Dede', L.; Quarteroni, A.
Impact of Atrial Fibrillation on Left Atrium Haemodynamics: A Computational Fluid Dynamics Study | Abstract | | We analyze left atrium haemodynamics, highlighting differences among healthy individuals and patients affected by atrial fibrillation. The computational study is based on patient-specific geometries of the left atria to simulate blood flow dynamics. We devise a novel procedure aimed at recovering the boundary conditions for the 3D haemodynamics simulations, particularly useful in absence of specific ones provided by clinical measurements. With this aim, we introduce a parametric definition of the atria displacement, and we employ a closed-loop lumped parameter model of the whole cardiocirculatory system conveniently tuned on the basis of the patient characteristics. We evaluate a number of fluid dynamics indicators for the atrial haemodynamics, validating our numerical results in terms of several clinical measurements; we investigate the impact of geometrical and clinical features on the risk of thrombosis. To analyse the correlation of thrombus formation with atrial fibrillation, coherently with the medical evidence, we propose a novel indicator, which we call age stasis and that arises from the combination of Eulerian and Lagrangian quantities. This indicator identifies regions where the slow flow cannot rinse the chamber properly, accumulating stale blood particles and creating optimal conditions for clot formation.
|
-
08/2022 - 02/02/2022
Gobat, G.; Opreni, A.; Fresca, S.; Manzoni, A.; Frangi, A.
Reduced order modeling of nonlinear microstructures through Proper Orthogonal Decomposition | Abstract | | We apply the Proper Orthogonal Decomposition (POD) method for the efficient simulation of several scenarios undergone by Micro-Electro Mechanical-Systems, involving nonlinearites of geometric and electrostatic nature. The former type of nonlinearity, associated to the large displacements of the devices, leads to polynomial terms up to cubic order that are reduced through exact projection onto a low dimensional subspace spanned by the Proper Orthogonal Modes (POMs). On the contrary, electrostatic nonlinearities are modeled resorting to precomputed manifolds in terms of the amplitudes of the electrically active POMs. We extensively test the reliability of the assumed linear trial space in challenging applications focusing on resonators, micromirrors and arches also displaying internal resonances. We discuss several options to generate the matrix of snapshots using both classical time marching schemes and more advanced Harmonic Balance (HB) approaches. Furthermore, we propose a comparison between the periodic orbits computed with POD and the invariant manifold approximated with Direct Normal Form approaches, further stressing the reliability of the technique and its remarkable predictive capabilities, e.g., in terms of estimation of the frequency response function of selected output quantities of interest |
-
07/2022 - 02/02/2022
Sinigaglia, C.; Quadrelli, D.E.; Manzoni, A.; Braghin, F.
Fast active thermal cloaking through PDE-constrained optimization and reduced-order modeling | Abstract | | In this paper we show how to efficiently achieve thermal cloaking from a computational standpoint in several virtual scenarios by controlling a distribution of active heat sources. We frame this problem in the setting of PDE-constrained optimization, where the reference field is the solution of the time-dependent heat equation in the absence of the object to cloak. The optimal control problem then aims at actuating the space-time control field so that the thermal field outside the obstacle is indistinguishable from the reference field. In particular, we consider multiple scenarios where material’s thermal diffusivity, source intensity and obstacle’s temperature are allowed to vary within a user-defined range. To tackle the thermal cloaking problem in a rapid and reliable way, we rely on a parametrized reduced order model built through the reduced basis method, thus entailing huge computational speedups compared to high-fidelity, full-order model exploiting the finite element method while dealing both with complex target shapes and disconnected control domains. |
-
06/2022 - 02/02/2022
Pozzi, G.; Grammatica, B.; Chaabane, L.; Catucci, M.; Mondino, A.; Zunino, P.; Ciarletta, P.
T cell therapy against cancer: a predictive diffuse-interface mathematical model informed by pre-clinical studies | Abstract | | T cell therapy has become a new therapeutic opportunity against solid cancers. Predicting T cell behaviour and efficacy would help therapy optimization and clinical implementation. In this work, we model responsiveness of mouse prostate adenocarcinoma to T cell-based therapies. The mathematical model is based on a Cahn-Hilliard diffuse interface description of the tumour, coupled with Keller-Segel type equations describing immune components dynamics. The model is fed by pre-clinical magnetic resonance imaging data describing anatomical features of prostate adenocarcinoma developed in the context of the Transgenic Adenocarcinoma of the Mouse Prostate model. We perform computational simulations based on the finite element method to describe tumor growth dynamics in relation to local T cells concentrations. We report that when we include in the model the possibility to activate tumor-associated vessels and by that increase the number of T cells within the tumor mass, the model predicts higher therapeutic effects (tumor regression) shortly after therapy administration. The simulated results are found in agreement with reported experimental data. Thus, this diffuse-interface mathematical model well predicts T cell behavior in vivo and represents a proof-of-concept for the role such predictive strategies may play in optimization of immunotherapy against cancer. |
|
