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 1362 prodotti
-
71/2026 - 24/09/2026
Radisic, I.; Regazzoni, F.; Dede', L.; Quarteroni, A.
A mathematical model for irreversible damage of the collagen scaffold in the myocardium | Abstract | | We propose a dissipative, irreversible damage model for anisotropic media in large deformations to address the damage process of the myocardium following a cardiac infarction. We model damage to the collagen in the cleavage planes resulting from an increased load to the passive tissue. Starting from variational principles, we derive a quasi-static differential model governing the irreversible evolution of the damage. We apply our model to two test cases. First, a rectangular passive slab geometry to which an indenter-like load is applied. Then, we couple our model with a model for the numerical simulation of left ventricular electromechanics. The computational results obtained with our novel mathematical model show that a redistribution of the load within a left ventricle occurring following a cardiac infarction can produce the damage to the collagen scaffold in the infarcted area consistent with experimental results. Our model and results support the evidence that passive load-dependent dissipative phenomena are implicated in post-infarction remodelling. |
-
70/2026 - 24/09/2026
Micheletti, S.
Longitudinal versus quasi-P waves in weakly anisotropic elasticity: operator comparison and long-time dephasing | Abstract | | In isotropic elasticity, P-wave displacement is longitudinal and curl-free. In anisotropic elasticity, the physical quasi-P (qP) mode is instead selected by the principal eigenvector of the Christoffel matrix and is generally not parallel to the wave normal. We compare, in homogeneous two-dimensional vertical transverse isotropy (VTI), the exact elastic qP propagator with a longitudinally constrained propagator that suppresses transverse displacement by construction. Both dynamics are placed in a periodic operator and variational framework. The qP space is an invariant spectral subspace of the elastic operator, whereas the longitudinal space yields a Galerkin restriction that is generally not invariant. The corresponding forms are coercive on mean-zero $H^1$, their generators are positive self-adjoint operators of order two, and both Cauchy problems are well posed and energy conservative. A two-dimensional Rayleigh-quotient identity expresses the longitudinal squared-speed deficit as the qP--qSV spectral gap multiplied by the squared polarization mismatch. Building on classical weak-anisotropy perturbation theory, we derive an explicit VTI specialization of the qP projector and phase symbol. The projector error and qSV contamination are $O(eta)$ in amplitude, whereas the qP phase-speed error is $O(eta^2)$. For longitudinal initial data, the field error admits an orthogonal decomposition into modal contamination and qP dephasing. A uniform estimate bounds it by $C|eta||u_0|+Cteta^2|u_0|_{H^1}$ on bounded times. After removal of the common carrier, an abstract two-branch multiplier theorem gives field-level limits on the intermediate scale $t=tau/eta$ and the diffractive scale $t=tau/eta^2$. The first yields a finite polarization--phase crossover, whereas the second yields a generally order-one dephasing profile. Exact-in-time Fourier experiments verify the identities, asymptotic coefficients, convergence rates, and two long-time limits. The main contribution is therefore not a new local qP approximation, but a quantitative error theory describing when a longitudinal elastic model approximates the exact qP branch and when its second-order phase accuracy is exhausted by long-time dephasing. |
-
69/2026 - 04/09/2026
Zecchi, A. A.; Ferro, N.; Perotto, S.
Complexity-controlled QUBO formulation for topology optimization via anisotropic mesh adaptation | Abstract | | QUBO formulation provides a natural route to address topology optimization of continuum structures with classical heuristic, quantum-inspired, and quantum annealing solvers. Their practical applicability, however, is still limited by the number of binary variables that can be handled, especially on near-term quantum hardware. In this work, we tackle this bottleneck by proposing QUBOSIMPATY algorithm which combines a QUBO formulation of the SIMP minimum compliance problem with anisotropic recovery-based mesh adaptation, enriched with metric-based complexity control.
Specifically, we exploit anisotropic mesh adaptation to reduce the number of density degrees of freedom, while retaining directional resolution along the sharp material-void interfaces that characterize optimized layouts. A metric rescaling procedure is applied to prescribe the number of mesh vertices, and hence the number of binary variables, so as to comply with the severe size constraints imposed by quantum annealing hardware.
Numerical tests in two and three dimensions are carried out with classical QUBO heuristics, including simulated annealing and tabu search. Moreover, a proof-of-concept experiment is performed on quantum annealing hardware. The results show that the proposed strategy improves the suppression of intermediate densities while keeping the optimization problem within the admissible size range. |
-
68/2026 - 04/09/2026
Patanè, G.; Greven, S.; Menafoglio, A.
Random mixtures in Bayes Hilbert spaces | Abstract | | We present a framework for the analysis and unmixing of random density mixtures in the Bayes Hilbert space. General identifiability results for mixtures in Hilbert spaces are established and applied to the Bayes Hilbert space setting. Building on these results, we propose a penalised maximum likelihood approach for the unmixing of Bayes Hilbert mixtures aimed at recovering the statistically space-efficient representation, together with a computationally efficient coordinate-wise maximisation algorithm for its implementation. The methodology is illustrated through a hyperspectral data application, where observations can be naturally embedded in the Bayes Hilbert space and analyzed in terms of distributional shape rather than amplitude. A complementary simulation study demonstrates the interpretability and practical performance of the proposed approach. |
-
67/2026 - 28/08/2026
Ciardulli, S.; Fontana, N.; Vantini, S.; Ieva F.
Generalized propensity score weighting for functional causal inference framework | Abstract | | Estimating causal effects in observational studies requires adjustment for confounding, a task that becomes challenging when the exposure is a function observed over a continuous domain rather than a scalar variable. We develop a functional propensity score weighting framework that achieves covariate balance by removing dependence between time-varying treatments and observed confounders, thereby enabling estimation of marginal causal effects in settings with functional treatments, covariates, and outcomes.
We propose a dual formulation of the weight estimation problem that yields a smooth unconstrained optimization and improves computational scalability. The proposed framework extends naturally to settings with time-varying covariates and to longitudinal outcomes via a function-on-function marginal structural model, allowing estimation of causal effect surfaces. The proposed method improves covariate balance, estimation accuracy, and computational efficiency compared to the existing approach and retains these properties when extended to functional covariates and outcomes. We apply the method to data from the UK Biobank to estimate the causal effect of body mass index trajectories on the risk of Type 2 Diabetes and on subsequent glycated hemoglobin trajectories, a functional measure of metabolic status. |
-
66/2026 - 28/07/2026
Antonietti, P.F.;Bonizzoni, F.;Corti, M.; De March, N.;Di Noto, S.;Regazzoni, F.
A stability-preserving polytopal discontinuous Galerkin method for the Fisher-Kolmogorov model with applications to neurodegenerative diseases | Abstract | | The Fisher-Kolmogorov model is one of the most widely used models in the study of neurodegenerative diseases, owing to its simple structure as a nonlinear reaction-diffusion equation. In particular, it is commonly employed to describe proteinopathies such as Alzheimer's and Parkinson's diseases. Under suitable assumptions, non-negativity of the solution is guaranteed at the continuous level, which is physically relevant since the solution represents a relative concentration. However, this property is not generally preserved at the discrete level, potentially leading to unphysical and unstable numerical approximations.
In this work, we analyze a modified version of the Fisher-Kolmogorov model that stabilizes the dynamics around the unstable equilibrium c=0. For the spatial discretization, we adopt a discontinuous Galerkin method on polygonal and polyhedral meshes, coupled with the Crank--Nicolson scheme for time integration. We derive stability and a-priori error estimates for the semi-discrete problem.
The theoretical findings are supported by numerical experiments, including convergence studies in both two and three dimensions. Finally, we validate the model through simulations of alpha-synuclein diffusion in a two-dimensional agglomerated brain section, demonstrating the high-order accuracy and robustness of the proposed method. |
-
65/2026 - 17/07/2026
De Sanctis, M.F.; Arnone, E.; Ieva, F.; Sangalli, L.M.
Modeling group heterogeneity in spatio-temporal data via physics-informed regression | Abstract | | We propose a physics-informed semiparametric framework for modeling spatiotemporal data with group structure. The approach extends classical mixed effects regression by incorporating a nonparametric space–time component, regularized through a partial differential equation to encode the underlying physical dynamics, while random effects capture group-specific variability. Estimation is carried out via a two-step procedure based on a functional extension of the Iteratively Reweighted Least Squares algorithm. We establish asymptotic properties of both fixed and random effect estimators and we assess the performance of the method through simulation studies against state-of-the-art alternatives. The proposed framework is applied to hourly nitrogen dioxide data over Lombardy (Italy), where random effects account for measurement heterogeneity across monitoring stations with different sensor technologies, demonstrating its effectiveness in capturing both physical dynamics and group heterogeneity. |
-
64/2026 - 15/07/2026
Radisic, I.; Tirotta, R.; Zingaro, A.; Pagani, S.; Dede', L.
Physics-Informed Neural Networks for the High-Resolution Reconstruction of Flow Measurement Indicators in Fluid Dynamics | Abstract | | Accurate, spatially resolved flow field measurements are essential for the reliable assessment of hemodynamic quantities in cardiovascular research and clinical practice. Experimental techniques, such as 4D flow MRI, PIV, or Doppler ultrasound, often yield data that are sparse, noisy, or under-resolved, particularly near vessel walls and in regions of complex flow. This limits the fidelity of distributed or derived hemodynamic indicators such as the wall shear stress and the clinical utility of such measurements. To address these challenges, we propose a physics-informed neural network (PINN) framework that integrates the incompressible Navier-Stokes equations with velocity measurements coming from experimental flow field data. By embedding physical laws into data, PINN enhances the reconstruction of velocity fields, enables the estimation of unmeasured quantities such as pressure and wall shear stress, and improves the spatial resolution of hemodynamic indicators. We show the effectiveness of our approach using both in silico and experimental data. First, we apply our method to the FDA nozzle benchmark, leveraging both control particle image velocimetry (PIV) measurements and computational fluid dynamics (CFD) simulations. Next, we apply our method to the more complex case of blood flow in an aneurysm model, exploiting in vitro 4D flow MRI data. In both cases, the synergy between data-driven learning and physics-based regularization yields results that align more closely with ground truth observations than standard CFD or pure data-driven approaches. Our findings highlight the potential of PINNs to improve the fidelity of under-resolved flow field measurements and yield spatially resolved hemodynamic indicators. |
|
