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 1249 prodotti
-
59/2022 - 29/08/2022
Boon, W. M.; Fumagalli, A.
A multipoint vorticity mixed finite element method for incompressible Stokes flow | Abstract | | We propose a mixed finite element method for Stokes flow with one degree of freedom per element and facet of simplicial grids. The method is derived by considering the vorticity-velocity-pressure formulation and eliminating the vorticity locally through the use of a quadrature rule. The discrete solution is pointwise divergence-free and the method is pressure robust. The theoretically derived convergence rates are confirmed by numerical experiments. |
-
58/2022 - 20/08/2022
Zingaro, A.; Bucelli, M.; Fumagalli, I.; Dede', L; Quarteroni, A.
Modeling isovolumetric phases in cardiac flows by an Augmented Resistive Immersed Implicit Surface Method | Abstract | | A major challenge in the computational fluid dynamics modeling of the heart function is the simulation of isovolumetric phases when the hemodynamics problem is driven by a prescribed boundary displacement.
During such phases, both atrioventricular and semilunar valves are closed: consequently, the ventricular pressure may not be uniquely defined, and spurious oscillations may arise in numerical simulations.}
In this paper, we propose a suitable modification of the Resistive Immersed Implicit Surface (RIIS) method (Fedele et al., 2017) by introducing a reaction term to correctly capture the pressure transients during isovolumetric phases. The method, that we call Augmented RIIS (ARIIS) method, extends the previously proposed ARIS method (This et al., 2020) to the case of a mesh which is not body-fitted to the valves. We test the proposed method on two different benchmark problems, including a new simplified problem that retains all the characteristics of a heart cycle. We apply the ARIIS method to a fluid dynamics simulation of a realistic left heart geometry, and we show that ARIIS allows to correctly simulate isovolumetric phases, differently from standard RIIS method. |
-
57/2022 - 15/08/2022
Ruffino, L.; Santoro, A.; Sparvieri, S.; Regazzoni, F.; Adebo, D.A.; Quarteroni, A.; Vergara, C.; Corno, A.F.
Computational analysis of cardiovascular effects of COVID- 19 infection in children | Abstract | | BACKGROUND. The COVID-19 disease can involve any body part; nevertheless, the most serious consequences affect respiratory and cardiocirculatory systems with variable symptoms. Although the effects of COVID-19 are not fully understood yet, clinical evidence has shown that the virus may cause acute myocardial injury and chronic damages to heart and
blood vessels. There is no or limited experience on pathophysiological effects of COVID-19 infection in children’s cardiovascular system.
OBJECTIVES. The aim of this work is to assess the effects of COVID-19 on the cardiovascular system in children, in terms, e.g., of increased pulmonary resistances, reduced cardiac contraction capacity.
METHODS. We used a computational model based on lumped parameters to describe the whole blood circulation. The model was calibrated to account for data coming from 5 child patients.
RESULTS. Our results highlighted that the effect of COVID-19 on the cardiovascular system in children was characterized by the reduction in cardiac blood pressure and volumes. In particular, we analyzed in detail two patients showing a correlation between myocardial compromise and severity of the infection.
CONCLUSIONS. This study demonstrates that COVID-19 infection causes a complex pathophysiological state to the cardiovascular system, both in asymptomatic and symptomatic pediatric patients. This information is very helpful to prevent long term cardiovascular complications of COVID-19 infection in children. A prospective study with regular cardiology follow-up is recommended. |
-
56/2022 - 11/08/2022
Africa, P.C.
lifex: a flexible, high performance library for the numerical solution of complex finite element problems | Abstract | | Numerical simulations are ubiquitous in mathematical and computational modeling, where many industrial and clinical applications are required to deal with multiphysics problems and with complex systems characterized by multiple spatial and temporal scales.
This document introduces the design and the capabilities of lifex, an open source C++ library for high performance finite element simulations of multiphysics, multiscale and multidomain problems. lifex offers a versatile solution to answer the emerging need for efficient computational tools that are also easily approachable by a wide community of users and developers. We showcase illustrative examples of use, benchmarks, advanced application scenarios and demonstrate its parallel performance up to thousands of cores. |
-
55/2022 - 11/08/2022
Cavinato, L.; Pegoraro, M.; Ragni, A.; Ieva, F.
Imaging-based representation and stratification of intra-tumor Heterogeneity via tree-edit distance | Abstract | | Personalized medicine is the future of medical practice. In oncology, tumor heterogeneity assessment represents a pivotal step for effective treatment planning and prognosis prediction. Despite new procedures for DNA sequencing and analysis, non-invasive methods for tumor characterization are needed to impact on daily routine. On purpose, imaging texture analysis is rapidly scaling, holding the promise to surrogate histopathological assessment of tumor lesions. In this work, we propose a tree-based representation strategy for describing intra-tumor heterogeneity of patients affected by metastatic cancer. We leverage radiomics information extracted from PET/CT imaging and we provide an exhaustive and easily readable summary of the disease spreading. We exploit this novel patient representation to perform cancer subtyping according to hierarchical clustering technique. To this purpose, a new heterogeneity-based distance between trees is defined and applied to a case study of Prostate Cancer (PCa). Clusters interpretation is explored in terms of concordance with severity status, tumor burden and biological characteristics. Results are promising, as the proposed method outperforms current literature approaches. Ultimately, the pro- posed methods draws a general analysis framework that would allow to extract knowledge from daily acquired imaging data of patients and provide insights for effective treatment planning. |
-
54/2022 - 11/08/2022
Bucelli, M.; Zingaro, A.; Africa, P. C.; Fumagalli, I.; Dede', L.; Quarteroni, A.
A mathematical model that integrates cardiac electrophysiology, mechanics and fluid dynamics: application to the human left heart | Abstract | | We propose a mathematical and numerical model for the simulation of the heart function that couples cardiac electrophysiology, active and passive mechanics and hemodynamics, and includes reduced models for cardiac valves and the circulatory system. Our model accounts for the major feedback effects among the different processes that characterize the heart function, including electro-mechanical and mechano-electrical feedback as well as force-strain and force-velocity relationships. Moreover, it provides a three-dimensional representation of both the cardiac muscle and the hemodynamics, coupled in a fluid-structure interaction (FSI) model. By leveraging the multiphysics nature of the problem, we discretize it in time with a segregated electrophysiology-force generation-FSI approach, allowing for efficiency and flexibility in the numerical solution. We employ a monolithic approach for the numerical discretization of the FSI problem. We use finite elements for the spatial discretization of those partial differential equations that contribute to the model. We carry out a numerical simulation on a realistic human left heart model, obtaining results that are qualitatively and quantitatively in agreement with physiological ranges and medical images. |
-
53/2022 - 28/07/2022
Antonietti, P.F; Cauzzi, C.; Mazzieri, I.; Melas L.; Stupazzini, M.
Numerical simulation of the Athens 1999 earthquake including simplified models of the Acropolis and the Parthenon: initial results and outlook | Abstract | | In this work we present a preliminary study of the seismic response of the Acropolis and of the Parthenon of Athens to the 1999 Mw 5.9 earthquake. The three-dimensional numerical model includes the surface topography of the Attica region, the seismogenic fault, and the most important geological units in the metropolitan area of Athens: the Acropolis hill and the Parthenon. The multiscale numerical model, designed in order to correctly propagate seismic waves up to $5Hz$, is solved through a discontinuous Galerkin spectral element method implemented in the open source library SPEED (speed.mox.polimi.it).
Numerical results show the effectiveness of this approach and highlight new challenges for dynamic soil-structure interaction problems at the regional scale. |
-
52/2022 - 28/07/2022
Fedele, M.; Piersanti, R.; Regazzoni, F.; Salvador, M.; Africa, P. C.; Bucelli, M.; Zingaro, A.; Dede', L.; Quarteroni, A.
A comprehensive and biophysically detailed computational model of the whole human heart electromechanics | Abstract | | While ventricular electromechanics is extensively studied in both physiological and pathological conditions, four-chamber heart models have only been addressed recently; most of these works however neglect atrial contraction. Indeed, as atria are characterized by a complex anatomy and a physiology that is strongly influenced by the ventricular function, developing computational models able to capture the physiological atrial function and atrioventricular interaction is very challenging. In this paper, we propose a biophysically detailed electromechanical model of the whole human heart that considers both atrial and ventricular contraction. Our model includes: i) an anatomically accurate whole-heart geometry; ii) a comprehensive myocardial fiber architecture; iii) a biophysically detailed microscale model for the active force generation; iv) a 0D closed-loop model of the circulatory system, fully-coupled with the mechanical model of the heart; v) the fundamental interactions among the different textit{core models}, such as the mechano-electric feedback or the fibers-stretch and fibers-stretch-rate feedbacks; vi) specific constitutive laws and model parameters for each cardiac region. Concerning the numerical discretization, we propose an efficient segregated-intergrid-staggered scheme and we employ recently developed stabilization techniques - regarding the circulation and the fibers-stretch-rate feedback - that are crucial to obtain a stable formulation in a four-chamber scenario. We are able to reproduce the healthy cardiac function for all the heart chambers, in terms of pressure-volume loops, time evolution of pressures, volumes and fluxes, and three-dimensional cardiac deformation, with unprecedented matching (to the best of our knowledge) with the expected physiology. We also show the importance of considering atrial contraction, fibers-stretch-rate feedback and suitable stabilization techniques, by comparing the results obtained with and without these features in the model. The proposed model represents the state-of-the-art electromechanical model of the iHEART ERC project - an Integrated Heart Model for the Simulation of the Cardiac Function - and is a fundamental step toward the building of physics-based digital twins of the human heart. |
|
