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 1249 products
-
65/2023 - 09/03/2023
Gatti, F.; de Falco, C.; Perotto, S.; Formaggia, L.; Pastor, M.
A scalable well-balanced numerical scheme for the modelling of two-phase shallow granular landslide consolidation | Abstract | | We introduce a new method to efficiently solve a variant of the Pitman-Le two-phase depth-integrated system of equations, for the simulation of a fast landslide consolidation process. In particular, in order to cope with the loss of hyperbolicity typical of this system, we generalize Pelanti’s proposition for the Pitman-Le model to the case of a non-null excess pore water pressure configuration. The variant of the Pitman-Le model is numerically solved by relying on the approach the authors set to discretize the corresponding single-phase model, jointly with a fictitious inter-phase drag force which avoids arising the spurious numerical oscillations induced by the loss of hyperbolicity. To verify the reliability of the proposed simulation tool, we first assess the accuracy and efficiency of the new method in ideal scenarios. In particular, we investigate the well-balancing property and provide some relevant scaling results for the parallel implementation of the method. Successively, we challenge the procedure on real configurations from the available literature. |
-
64/2023 - 08/28/2023
Heltai, L.; Zunino, P.
Reduced Lagrange multiplier approach for non-matching coupling of mixed-dimensional domains | Abstract | | Many physical problems involving heterogeneous spatial scales, such as the flow through fractured porous media, the study of fiber-reinforced materials, or the modeling of the small circulation in living tissues -- just to mention a few examples -- can be described as coupled partial differential equations defined in domains of heterogeneous dimensions that are embedded into each other. This formulation is a consequence of geometric model reduction techniques that transform the original problems defined in complex three-dimensional domains into more tractable ones. The definition and the approximation of coupling operators suitable for this class of problems is still a challenge. We develop a general mathematical framework for the analysis and the approximation of partial differential equations coupled by non-matching constraints across different dimensions, focusing on their enforcement using Lagrange multipliers. In this context, we address in abstract and general terms the well-posedness, stability, and robustness of the problem with respect to the smallest characteristic length of the embedded domain. We also address the numerical approximation of the problem and we discuss the inf-sup stability of the proposed numerical scheme for some representative configuration of the embedded domain. The main message of this work is twofold: from the standpoint of the theory of mixed-dimensional problems, we provide general and abstract mathematical tools to formulate coupled problems across dimensions. From the practical standpoint of the numerical approximation, we show the interplay between the mesh characteristic size, the dimension of the Lagrange multiplier space, and the size of the inclusion in representative configurations interesting for applications. The latter analysis is complemented with illustrative numerical examples. |
-
63/2023 - 08/22/2023
Antonietti, P. F.; Matalon, P.; Verani M.
Iterative solution to the biharmonic equation in mixed form discretized by the Hybrid High-Order method | Abstract | | We consider the solution to the biharmonic equation in mixed form discretized by the Hybrid High-Order (HHO) methods. The two resulting second-order elliptic problems can be decoupled via the introduction of a new unknown, corresponding to the boundary value of the solution of the first Laplacian problem. This technique yields a global linear problem that can be solved iteratively via a Krylov-type method. More precisely, at each iteration of the scheme, two second-order elliptic problems have to be solved, and a normal derivative on the boundary has to be computed. In this work, we specialize this scheme for the HHO discretization. To this aim, an explicit technique to compute the discrete normal derivative of an HHO solution of a Laplacian problem is proposed. Moreover, we show that the resulting discrete scheme is well-posed. Finally, a new preconditioner is designed to speed up the convergence of the Krylov method. Numerical experiments assessing the performance of the proposed iterative algorithm on both two- and three-dimensional test cases are presented. |
-
62/2023 - 08/08/2023
Zappon, E.; Manzoni, A.; Quarteroni, A.
A staggered-in-time and non-conforming-in-space numerical framework for realistic cardiac electrophysiology outputs | Abstract | | Computer-based simulations of non-invasive cardiac electrical outputs, such as electrocardiograms and body surface potential maps, usually entail severe computational costs due to the need of capturing fine-scale processes and to the complexity of the heart-torso morphology. In this work, we model cardiac electrical outputs by employing a coupled model consisting of a reaction-diffusion model - either the bidomain model or the most efficient pseudo-bidomain model - on the heart, and an elliptic model in the torso. We then solve the coupled problem with a segregated and staggered in-time numerical scheme, that allows for independent and infrequent solution in the torso region. To further reduce the computational load, main novelty of this work is in introduction of an interpolation method at the interface between the heart and torso domains, enabling the use of non-conforming meshes, and the numerical framework application to realistic cardiac and torso geometries. The reliability and efficiency of the proposed scheme is tested against the corresponding state-of-the-art bidomain-torso model. Furthermore, we explore the impact of torso spatial discretization and geometrical non-conformity on the model solution and the corresponding clinical outputs. The investigation of the interface interpolation method provides insights into the influence of torso spatial discretization and of the geometrical non-conformity on the simulation results and their clinical relevance. |
-
61/2023 - 08/04/2023
Afriaca, P.C.A; Piersanti, R.; Regazzoni, F.; Bucelli, M.; Salvador, M.; Fedele, M.; Pagani, S.; Dede', L.; Quarteroni, A.
lifex-ep: a robust and efficient software for cardiac electrophysiology simulations | Abstract | | Simulating the cardiac function requires the numerical solution of multi-physics and multi-scale mathematical models. This underscores the need for streamlined, accurate, and high-performance computational tools. Despite the dedicated endeavors of various research teams, comprehensive and user-friendly software programs for cardiac simulations, capable of accurately replicating both physiological and pathological conditions, are still in the process of achieving full maturity within the scientific community. This work introduces lifex-ep, a publicly available software for numerical simulations of the electrophysiology activity of the cardiac muscle, under both physiological and pathological conditions. lifex-ep employs the monodomain equation to model the heart’s electrical activity. It incorporates both phenomenological and second-generation ionic models. These models are discretized using the Finite Element method on tetrahedral or hexahedral meshes. Additionally, lifex-ep integrates the generation of myocardial fibers based on Laplace-Dirichlet Rule-Based Methods, previously released in Africa et al., 2023, within lifex-fiber. As an alternative, users can also choose to import myofibers from a file. This paper provides a concise overview of the mathematical models and numerical methods underlying lifex-ep, along with comprehensive
implementation details and instructions for users. lifex-ep features
exceptional parallel speedup, scaling efficiently when using up to thousands of cores, and its implementation has been verified against an established benchmark problem for computational electrophysiology. We
showcase the key features of lifex-ep through various idealized and
realistic simulations conducted in both physiological and pathological
scenarios. Furthermore, the software offers a user-friendly and flexible
interface, simplifying the setup of simulations using self-documenting
parameter files. lifex-ep provides easy access to cardiac electrophysiology simulations for a wide user community. It offers a computational tool that integrates models and accurate methods for simulating cardiac electrophysiology within a high-performance framework, while maintaining a user-friendly interface. lifex-ep represents a valuable tool for conducting in silico patient-specific simulations. |
-
60/2023 - 08/02/2023
Negrini G.; Parolini N.; Verani M.
The Rhie-Chow stabilized Box Method for the Stokes problem | Abstract | | The Finite Volume method (FVM) is widely adopted in many different applications because of its built-in conservation properties, its ability to deal with arbitrary mesh and its computational efficiency. In this work, we consider the Rhie-Chow stabilized Box Method (RCBM) for the approximation of the Stokes problem. The Box Method (BM) is a piecewise linear Petrov-Galerkin formulation on the Voronoi dual mesh of a Delaunay triangulation, whereas the Rhie-Chow (RC) stabilization is a well known stabilization technique for FVM.
The first part of the paper provides a variational formulation of the RC stabilization and discusses the validity of crucial properties relevant for the well-posedeness and convergence of RCBM. Moreover, a numerical exploration of the convergence properties of the method on 2D and 3D test cases is presented. The last part of the paper considers the theoretically justification of the well-posedeness of RCBM and the experimentally observed convergence rates. This latter justification hinges upon suitable assumptions, whose validity is numerically explored. |
-
59/2023 - 08/01/2023
Corti, M.; Bonizzoni, F.; Antonietti, P.F.
Structure Preserving Polytopal Discontinuous Galerkin Methods for the Numerical Modeling of Neurodegenerative Diseases | Abstract | | Many neurodegenerative diseases are connected to the spreading of misfolded prionic proteins. In this paper, we analyse the process of misfolding and spreading of both alpha-synuclein and Amyloid-beta, related to Parkinson's and Alzheimer's diseases, respectively. We introduce and analyze a positivity-preserving numerical method for the discretization of the Fisher-Kolmogorov equation, modelling accumulation and spreading of prionic proteins. The proposed approximation method is based on the discontinuous Galerkin method on polygonal and polyhedral grids for space discretization and on theta-method time integration scheme. We prove the existence of the discrete solution and a convergence result where the Implicit Euler scheme is employed for time integration. We show that the proposed approach is structure-preserving, in the sense that it guaranteed that the discrete solution is non-negative, a feature that is of paramount importance in practical application. The numerical verification of our numerical model is performed both using a manufactured solution and considering wavefront propagation in two-dimensional polygonal grids. Next, we present a simulation of alpha-synuclein spreading in a two-dimensional brain slice in the sagittal plane. The polygonal mesh for this simulation is agglomerated maintaining the distinction of white and grey matter, taking advantage of the flexibility of PolyDG methods in the mesh construction. Finally, we simulate the spreading of Amyloid-beta in a patient-specific setting by using a three-dimensional geometry reconstructed from magnetic resonance images and an initial condition reconstructed from positron emission tomography. Our numerical simulations confirm that the proposed method is able to capture the evolution of Parkinson's and Alzheimer's diseases. |
-
57/2023 - 07/24/2023
Regazzoni, F.
An optimally convergent Fictitious Domain method for interface problems | Abstract | | We introduce a novel Fictitious Domain (FD) unfitted method for interface problems that achieves optimal convergence without the need for adaptive mesh refinements nor enrichments of the Finite Element spaces. The key aspect of the proposed method is that it extends the solution into the fictitious domain in a way that ensures high global regularity. Continuity of the solution across the interface is enforced through a boundary Lagrange multiplier. The subdomains coupling, however, is not achieved by means of the duality pairing with the Lagrange multiplier, but through an $L^2$ product with the $H^1$ Riesz representative of the latter, thus avoiding gradient jumps across the interface. Thanks to the enhanced regularity, the proposed method attains an increase, with respect to standard FD methods, of up to one order of convergence in energy norm. The Finite Element formulation of the method is presented, followed by its analysis. Numerical tests demonstrate its effectiveness. |
|
