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 1148 prodotti
-
103/2023 - 15/12/2023
Dimola N.; Kuchta M.; Mardal K.A.; Zunino P.
Robust Preconditioning of Mixed-Dimensional PDEs on 3d-1d domains coupled with Lagrange Multipliers | Abstract | | In the context of micro-circulation, the coexistence of two distinct length scales - the vascular radius and the tissue/organ scale - with a substantial difference in magnitude, poses significant challenges. To handle slender inclusions and simplify
the geometry involved, a technique called topological dimensionality reduction is used, which suppresses the manifold dimensions associated with the smaller characteristic length. However, the algebraic structure of the resulting discretized system presents a challenge in constructing efficient solution algorithms. This chapter addresses this challenge by developing a robust preconditioner for the 3d-1d problem using the operator preconditioning technique. The robustness of the preconditioner is demonstrated with respect to the problem parameters, except for the vascular radius. The vascular radius, as demonstrated, plays a fundamental role in the mathematical well-posedness of the problem and the effectiveness of the preconditioner. |
-
98/2023 - 14/12/2023
Lespagnol, F.; Grandmont, C.; Zunino, P.; Fernandez, M.A.
A mixed-dimensional formulation for the simulation of slender structures immersed in an incompressible flow | Abstract | | We consider the simulation of slender structures immersed in a three-dimensional (3D) flow. By exploiting the special geometric configuration of the slender structures, this particular problem can be modeled by mixed-dimensional coupled equations (3D for the fluid and 1D for the solid).
Several challenges must be faced when dealing with this type of problems. From a mathematical point of view, these include defining wellposed trace operators of codimension two. On the computational standpoint, the non-standard mathematical formulation makes it difficult to ensure the accuracy of the solutions obtained with the mixed-dimensional discrete formulation as compared to a fully resolved one.
We establish the continuous formulation using the Navier-Stokes equations for the fluid and a Timoshenko beam model for the structure. We complement these models with a mixed-dimensional version of the fluid-structure interface conditions, based on the projection of kinematic coupling conditions on a finite-dimensional Fourier space. Furthermore, we develop a discrete formulation within the framework of the finite element method, establish the energy stability of the scheme, provide extensive numerical evidence of the accuracy of the discrete formulation, notably with respect to a fully resolved (ALE based) model and a standard reduced modeling approach. |
-
100/2023 - 14/12/2023
Vitullo, P.; Cicci, L.; Possenti, L.; Coclite, A.; Costantino, M.L.; Zunino, P.
Sensitivity analysis of a multi-physics model for the vascular microenvironment | Abstract | | The vascular microenvironment is the scale at which microvascular transport, interstitial tissue properties and cell metabolism interact. The vascular microenvironment has been widely studied by means of quantitative approaches, including multi-physics mathematical models as it is a central system for the pathophysiology of many diseases, such as cancer.
The microvascular architecture is a key factor for the fluid balance and mass transfer in the vascular microenvironment, together with the physical parameters characterizing the vascular wall and the interstitial tissue. The scientific literature of this field has witnessed a long debate about which factor of this multifaceted system is the most relevant.
The purpose of this work is to combine the interpretative power of an advanced multi-physics model of the vascular microenvironment with state of the art, robust sensitivity analysis methods, in order to determine what factors affect the most some quantity of interest, related in particular to cancer treatment. We are particularly interested in comparing the factors related to the microvascular architecture with the ones affecting the physics of microvascular transport.
Ultimately, this work will provide further insight of how the vascular microenvironment affects cancer therapies, such as chemotherapy, radiotherapy or immunotherapy. |
-
96/2023 - 28/11/2023
Bonetti, S.; Botti, M.; Antonietti, P.F.
Robust discontinuous Galerkin-based scheme for the fully-coupled non-linear thermo-hydro-mechanical problem | Abstract | | We present and analyze a discontinuous Galerkin method for the numerical modeling of the non-linear fully-coupled thermo-hydro-mechanic problem. We propose a high-order symmetric weighted interior penalty scheme that supports general polytopal grids and is robust with respect to strong heteorgeneities in the model coefficients. We focus on the treatment of the non-linear convective transport term in the energy conservation equation and we propose suitable stabilization techniques that make the scheme robust for advection-dominated regimes. The stability analysis of the problem and the convergence of the fixed-point linearization strategy are addressed theoretically under mild requirements on the problem's data. A complete set of numerical simulations is presented in order to assess the convergence and robustness properties of the proposed method. |
-
95/2023 - 28/11/2023
Barnafi, N. A.; Regazzoni, F.; Riccobelli, D.
Reconstructing relaxed configurations in elastic bodies: mathematical formulation and numerical methods for cardiac modeling | Abstract | | Modeling the behavior of biological tissues and organs often necessitates the knowledge of their shape in the absence of external loads. However, when their geometry is acquired in-vivo through imaging techniques, bodies are typically subject to mechanical deformation due to the presence of external forces, and the load-free configuration needs to be reconstructed. This paper addresses this crucial and frequently overlooked topic, known as the inverse elasticity problem (IEP), by delving into both theoretical and numerical aspects, with a particular focus on cardiac mechanics. In this work, we extend Shield's seminal work to determine the structure of the IEP with arbitrary material inhomogeneities and in the presence of both body and active forces. These aspects are fundamental in computational cardiology, and we show that they may break the variational structure of the inverse problem. In addition, we show that the inverse problem might be ill-posed, even in the presence of constant Neumann boundary conditions and a polyconvex strain energy functional. We then present the results of extensive numerical tests to validate our theoretical framework, and to characterize the computational challenges associated with a direct numerical approximation of the IEP. Specifically, we show that this framework outperforms existing approaches both in terms of robustness and optimality, such as Sellier's iterative procedure, even when the latter is improved with acceleration techniques. A notable discovery is that multigrid preconditioners are, in contrast to standard elasticity, not efficient, and domain decomposition methods provide a much more reliable alternative. Finally, we successfully address the IEP for a full-heart geometry, demonstrating that the IEP formulation can compute the stress-free configuration in real-life scenarios where Sellier's algorithm proves inadequate. |
-
93/2023 - 25/11/2023
Andrini, D.; Magri, M.; Ciarletta, P.
Optimal surface clothing with elastic nets | Abstract | | The clothing problem aims at identifying the shape of a planar fabric for covering a target surface in the three-dimensional space. It poses significant challenges in various applications, ranging from fashion industry to digital manufacturing. Here, we propose a novel inverse design approach to the elastic clothing problem that is formulated as a constrained optimization problem. We assume that the textile behaves as an orthotropic, nonlinear elastic surface with fibers distributed along its warp and weft threads, and we enforce mechanical equilibrium as a variational problem. The target surface is frictionless, except at its boundary where the textile is pinned, imposing a unilateral obstacle condition for the reactive forces at the target surface. The constrained optimization problem also accounts for an elongation condition of the warp and weft fibers, possibly with bounded shearing angle. We numerically solve the resulting constrained optimization problem by means of a gradient descent algorithm. The numerical results are first validated against known clothing solutions for Chebyshev nets, taking the limit of inextensible fibers. We later unravel the interplay between thread and shear stiffness for driving the optimal cloth shape covering the hemisphere and the hemicatenoid. We show how the metric of these target surfaces strongly affects the resulting distribution of the reaction forces. When considering the limit of covering the full sphere, we show how clothing with elastic nets allows to avoid the onset of singularities in the corresponding Chebyshev net, by developing corners at the cloth boundary. |
-
92/2023 - 16/11/2023
Burzacchi, A.; Rossi, L.; Agasisti, T.; Paganoni, A. M.; Vantini, S.
Commuting time as a determinant of higher education students' performance: the case of Politecnico di Milano | Abstract | | Despite its crucial role in students' daily lives, commuting time remains an underexplored dimension in higher education research. To address this gap, this study focuses on challenges that students face in urban environments and investigates the impact of commuting time on the academic performance of first-year bachelor students of Politecnico di Milano, Italy.
This research employs an innovative two-step methodology. In the initial phase, machine learning algorithms trained on GPS data from anonymous users are used to construct accessibility maps to the university and to obtain an estimate of students' commuting times. In the subsequent phase, authors utilize polynomial linear mixed-effects models and investigate the factors influencing students' academic performance, with a particular emphasis on commuting time. Notably, this investigation incorporates a causal framework, which enables the establishment of causal relationships between commuting time and academic outcomes.
The findings underscore the significant impact of travel time on students' performance and may support policies and implications aiming at improving students' educational experience in metropolitan areas.
The study's innovation lies both in its exploration of a relatively uncharted factor and the novel methodologies applied in both phases. |
-
90/2023 - 10/11/2023
Gregorio, C.; Baj, G.; Barbati, G.; Ieva, F.
Dynamic treatment effect phenotyping through functional survival analysis | Abstract | | In recent years, research interest in personalised treatments has been growing. However, treatment effect heterogeneity and possibly time-varying treatment effects are still often overlooked in clinical studies. Statistical tools are needed for the identification of treatment response patterns, taking into account that treatment response is not constant over time. We aim to provide an innovative method to obtain dynamic treatment effect phenotypes on a time-to-event outcome, conditioned
on a set of relevant effect modifiers. The proposed method does not require the assumption of proportional hazards for the treatment effect, which is rarely realistic. We propose a spline-based survival neural network, inspired by the Royston-Parmar survival model, to estimate time-varying conditional treatment effects. We then exploit the functional nature of the resulting estimates to apply a functional clustering of the treatment effect curves in order to identify different patterns of treatment effects. The application that motivated this work is the discontinuation of treatment with Mineralocorticoid receptor Antagonists (MRAs) in patients with heart failure, where there is no clear evidence as to which patients it is the safest choice to discontinue treatment and, conversely, when it leads to a higher risk of adverse events. The data come from an electronic health record database. A simulation study was performed to assess the performance of the spline-based neural network and the stability of the treatment response phenotyping procedure. In light of the results, the suggested approach has the potential to support personalized medical choices by assessing treatment responses in various medical contexts over a period of time. |
|