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 1239 prodotti
-
80/2021 - 02/12/2021
Sollini, M., Bartoli, F., Cavinato, L., Ieva, F., Ragni, A., Marciano, A., Zanca, R., Galli, L., Paiar, F., Pasqualetti , F. and Erba P. A.
[18F]FMCH PET/CT biomarkers and similarity analysis to refine the definition of oligometastatic prostate cancer | Abstract | | Background:The role of image-derived biomarkers in recurrent oligometastatic Prostate Cancer (PCa) is unexplored. This paper aimed to evaluate [18F]FMCH PET/CT radiomic analysis in patients with recurrent PCa after primary radical therapy. Specifically, we tested intra-patient lesions similarity in oligometastatic and plurimetastatic PCa, comparing the two most used definitions of oligometastatic disease.Methods:PCa patients eligible for [18F]FMCH PET/CT presenting biochemical failure after first-line curative treat-ments were invited to participate in this prospective observational trial. PET/CT images of 92 patients were visually and quantitatively analyzed. Each patient was classified as oligometastatic or plurimetastatic according to the total number of detected lesions (up to 3 and up to 5 or > 3 and > 5, respectively). Univariate and intra-patient lesions’ similarity analysis were performed.Results: [18F]FMCH PET/CT identified 370 lesions, anatomically classified as regional lymph nodes and distant metastases. Thirty-eight and 54 patients were designed oligometastatic and plurimetastatic, respectively, using a 3-lesion threshold. The number of oligometastic scaled up to 60 patients (thus 32 plurimetastatic patients) with a 5-lesion threshold. Similarity analysis showed high lesions’ heterogeneity. Grouping patients according to the number of metastases, patients with oligometastatic PCa defined with a 5-lesion threshold presented lesions heterogene-ity comparable to plurimetastic patients. Lesions within patients having a limited tumor burden as defined by three lesions were characterized by less heterogeneity.Conclusions:We found a comparable heterogeneity between patients with up to five lesions and plurimetastic patients, while patients with up to three lesions were less heterogeneous than plurimetastatic patients, featuring dif-ferent cells phenotypes in the two groups. Our results supported the use of a 3-lesion threshold to define oligometa-static PCa. |
-
79/2021 - 02/12/2021
Ferraccioli, F.; Sangalli, L.M.; Finos, L.
Some first inferential tools for spatial regression with differential regularization | Abstract | | Spatial regression with differential regularization is an innovative approach at the crossroad between functional data analysis and spatial data analysis. These models have been shown to be numerically efficient and capable to handle complex applied problems. On the other hand, their theoretical properties are still largely unexplored. Here we consider the discrete estimators in spatial regression models with differential regularization, obtained after numerical discretization, using an expansion on a finite element basis. We study the consistency and the asymptotic normality of these discrete estimators. We also propose a nonparametric test procedure for the linear part of the models, based on random sign-flipping of the score components. The test exploits an appropriate decomposition of the smoothing matrix, in order to reduce the effect of the spatial dependence, without any parametric assumption on the form of the correlation structure. The proposed test is shown to be superior to parametric alternatives. |
-
78/2021 - 02/12/2021
Bucelli, M.; Dede', L.; Quarteroni, A.; Vergara, C.
Partitioned and monolithic algorithms for the numerical solution of cardiac fluid-structure interaction | Abstract | | We review and compare different fluid-structure interaction (FSI) numerical methods in the context of heart modeling, aiming at assessing their computational efficiency for cardiac numerical simulations and selecting the most appropriate method for heart FSI. Blood dynamics within the human heart is characterized by active muscular action, during both contraction and relaxation phases of the heartbeat. The efficient solution of the FSI problem in this context is challenging, due to the added-mass effect (caused by the comparable densities of fluid and solid, typical of biomechanics) and to the complexity, nonlinearity and anisotropy of cardiac consitutive laws. In this work, we review existing numerical coupling schemes for FSI in the two classes of strongly-coupled partitioned and monolithic schemes. The schemes are compared on numerical tests that mimic the flow regime characterizing the heartbeat in a human ventricle, during both systole and diastole. Active mechanics is treated in both the active stress and active strain frameworks. Computational costs suggest the use of a monolithic method. We employ it to simulate a full heartbeat of a human ventricle, showing how it allows to efficiently obtain physiologically meaningful results. |
-
77/2021 - 28/11/2021
Guo, M.; Manzoni, A.; Amendt, M.; Conti, P.; Hesthaven, J.S.
Multi-fidelity regression using artificial neural networks: efficient approximation of parameter-dependent output quantities | Abstract | | Highly accurate numerical or physical experiments are often very time-consuming or expensive to obtain. When time or budget restrictions prohibit the generation of additional data, the amount of available samples may be too limited to provide satisfactory model results. Multi-fidelity methods deal with such problems by incorporating information from other sources, which are ideally well-correlated with the high-fidelity data, but can be obtained at a lower cost. By leveraging correlations between different data sets, multi-fidelity methods often yield superior generalization when compared to models based solely on a small amount of high fidelity data. In the current work, we present the use of artificial neural networks applied to multi- fidelity regression problems. By elaborating a few existing approaches, we propose new neural network architectures for multi-fidelity regression. The introduced models are compared against a traditional multi- fidelity regression scheme - co-kriging. A collection of artificial benchmarks are presented to measure the performance of the analyzed models. The results show that cross-validation in combination with Bayesian optimization leads to neural network models that outperform the co-kriging scheme. Additionally, we show an application of multi-fidelity regression to an engineering problem. The propagation of a pressure wave into an acoustic horn with parametrized shape and frequency is considered, and the index of reflection intensity is approximated using the proposed multi-fidelity models. A finite element, full-order model and a reduced-order model built through the reduced basis method are adopted as the high- and low-fidelity, respectively. It is shown that the multi-fidelity neural networks return outputs that achieve a comparable accuracy to those from the expensive, full-order model, using only very few full-order evaluations combined with a larger amount of inaccurate but cheap evaluations of the reduced order model. |
-
76/2021 - 28/11/2021
Ponti, L.; Perotto, S.; Sangalli, L.M.
A PDE-regularized smoothing method for space-time data over manifolds with application to medical data | Abstract | | We propose an innovative statistical-numerical method to model spatio-temporal data, observed over a generic two-dimensional Riemanian manifold.
The proposed approach consists of a regression model completed with a regularizing term based on the heat equation.
The model is discretized through a finite element scheme set on the manifold, and solved by resorting to a fixed point-based iterative algorithm.
This choice leads to a procedure which is highly efficient when compared with a monolithic approach, and which allows us to deal with
massive datasets. After a preliminary assessment on simulation study cases, we investigate the performance of the new estimation tool
in practical contexts by dealing with neuroimaging and hemodynamic data. |
-
75/2021 - 28/11/2021
Cicci, L.; Fresca, S.; Pagani, S.; Manzoni, A.; Quarteroni, A.
Projection-based reduced order models for parameterized nonlinear time-dependent problems arising in cardiac mechanics | Abstract | | The numerical simulation of several virtual scenarios arising in cardiac mechanics poses a computational challenge that can be alleviated if traditional full-order models (FOMs) are replaced by reduced order models (ROMs). For example, in the case of problems involving a vector of input parameters related, e.g., to material coefficients, projection-based ROMs provide mathematically rigorous physics-driven surrogate ROMs. In this work we demonstrate how, once trained, ROMs yield extremely accurate predictions (according to a prescribed tolerance) - yet cheaper than the ones provided by FOMs - of the structural deformation of the left ventricular tissue over an entire heartbeat, and of related output quantities of interest, such as the pressure-volume loop, for any desired input parameter values within a prescribed parameter range. However, the construction of ROM approximations for time-dependent cardiac mechanics is not straightforward, because of the highly nonlinear and multiscale nature of the problem, and almost never addressed. Our approach relies on the reduced basis method for parameterized partial differential equations. This technique performs a Galerkin projection onto a low-dimensional space for the displacement variable; the reduced space is built from a set of solution snapshots - obtained for different input parameter values and time instances - of the high-fidelity FOM, through the proper orthogonal decomposition technique. Then, suitable hyper-reduction techniques, such as the Discrete Empirical Interpolation Method, are exploited to efficiently handle nonlinear and parameter-dependent terms. In this work we show how a fast and reliable approximation of the time-dependent cardiac mechanical model can be achieved by a projection-based ROM, taking into account both passive and active mechanics for the left ventricle providing all the building blocks of the methodology, and highlighting those challenging aspects that are still open. |
-
74/2021 - 28/11/2021
Orlando,G.; Barbante, P. F.; Bonaventura, L.
An efficient IMEX-DG solver for the compressible Navier-Stokes equations with a general equation of state | Abstract | | We propose an efficient, accurate and robust IMEX solver for the compressible Navier-Stokes equation with general equation of state. The method, which is based on an h-adaptive Discontinuos Galerkin spatial discretization and on an Additive Runge Kutta IMEX method for time discretization, is tailored for low Mach number applications and allows to simulate low Mach regimes at a significantly reduced computational cost, while maintaining full second order accuracy also for higher Mach number regimes. The method has been implemented in the framework of the deal.II numerical library, whose adaptive mesh refinement capabilities are employed to enhance efficiency. Refinement indicators appropriate for real gas phenomena have been introduced. A number of numerical experiments on classical benchmarks for compressible flows and their extension to real gases demonstrate the properties of the proposed method. |
-
73/2021 - 28/11/2021
Marcinno, F.; Zingaro, A.; Fumagalli, I.; Dede', L.; Vergara, C.
A computational study of blood flow dynamics in the pulmonary arteries | Abstract | | In this work we study for the first time the blood dynamics in the pulmonary arteries
by means of a 3D-0D geometric multiscale approach, where a detailed 3D model for the pulmonary arteries is coupled with
a lumped parameters (0D) model of the cardiocirculatory system.
We propose to investigate two strategies for the numerical solution of the 3D-0D coupled problem: a Splitting Algorithm,
where information are exchanged between 3D and 0D models at each time step at the interfaces, and a One-Way Algorithm,
where the 0D is solved first off-line.
In our numerical experiments performed in a realistic patient-specific 3D domain with a physiologically calibrated 0D model,
we discuss first the issue on instabilities that may arise when not suitable connections are considered between 3D and 0D models;
second we compare the performance and accuracy of the two proposed numerical strategies. Finally, we report a
comparison between an healthy and an hypertensive case, providing a preliminary result highlighting how our method could
be used in future for clinical purposes. |
|
