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 1253 products
-
69/2023 - 09/20/2023
Ferro, N.; Micheletti, S.; Parolini, N.; Perotto, S.; Verani, M.; Antonietti, P. F.
Level set-fitted polytopal meshes with application to structural topology optimization | Abstract | | We propose a method to modify a polygonal mesh in order to fit the zero-isoline of a level set function by extending a standard body-fitted strategy to a tessellation with arbitrarily-shaped elements. The novel level set-fitted approach, in combination with a Discontinuous Galerkin finite element approximation, provides an ideal setting to model physical problems characterized by embedded or evolving complex geometries, since it allows skipping any mesh post-processing in terms of grid quality. The proposed methodology is firstly assessed on the linear elasticity equation, by verifying the approximation capability of the level set-fitted approach when dealing with configurations with heterogeneous material properties. Successively, we combine the level set-fitted methodology with a minimum compliance topology optimization technique, in order to deliver optimized layouts exhibiting crisp boundaries and reliable mechanical performances. An extensive numerical test campaign confirms the effectiveness of the proposed method. |
-
68/2023 - 09/14/2023
Vitullo, P.; Colombo, A.; Franco, N.R.; Manzoni, A.; Zunino, P.
Nonlinear model order reduction for problems with microstructure using mesh informed neural networks | Abstract | | Many applications in computational physics involve approximating problems with microstructure, characterized by multiple spatial scales in their data. However, these numerical solutions are often computationally expensive due to the need to capture fine details at small scales. As a result, simulating such phenomena becomes unaffordable for many-query applications, such as parametrized systems with multiple scale-dependent features. Traditional projection-based reduced order models (ROMs) fail to resolve these issues, even for second-order elliptic PDEs commonly found in engineering applications. To address this, we propose an alternative nonintrusive strategy to build a ROM, that combines classical proper orthogonal decomposition (POD) with a suitable neural network (NN) model to account for the small scales. Specifically, we employ sparse mesh-informed neural networks (MINNs), which handle both spatial dependencies in the solutions and model parameters simultaneously. We evaluate the performance of this strategy on benchmark problems and then apply it to approximate a real-life problem involving the impact of microcirculation in transport phenomena through the tissue microenvironment. |
-
66/2023 - 09/05/2023
Fresca, S.; Gobat, G.; Fedeli, P.; Frangi, A.; Manzoni, A.
Deep learning-based reduced order models for the real-time simulation of the nonlinear dynamics of microstructures | Abstract | | We propose a non-intrusive Deep Learning-based Reduced Order Model (DL-ROM) capable of capturing the complex dynamics of mechanical systems showing inertia and geometric nonlinearities. In the first phase, a limited number of high fidelity snapshots are used to generate a POD-Galerkin ROM which is subsequently exploited to generate the data, covering the whole parameter range, used in the training phase of the DL-ROM. A convolutional autoencoder is employed to map the system response onto a low-dimensional representation and, in parallel, to model the reduced nonlinear trial manifold. The system dynamics on the manifold is described by means of a deep feedforward neural network that is trained together with the autoencoder. The strategy is benchmarked against high fidelity solutions on a clamped-clamped beam and on a real micromirror with softening response and multiplicity of solutions. By comparing the different computational costs, we discuss the impressive gain in performance and show that the DL-ROM truly represents a real-time tool which can be profitably and efficiently employed in complex system-level simulation procedures for design and optimisation purposes. |
-
67/2023 - 09/05/2023
Conti, P.; Guo, M.; Manzoni, A.; Hesthaven, J.S.
Multi-fidelity surrogate modeling using long short-term memory networks | Abstract | | When evaluating quantities of interest that depend on the solutions to differential equations, we inevitably face the trade-off between accuracy and efficiency. Especially for parametrized, time-dependent problems in engineering computations, it is often the case that acceptable computational budgets limit the availability of high-fidelity, accurate simulation data. Multi-fidelity surrogate modeling has emerged as an effective strategy to overcome this difficulty. Its key idea is to leverage many low-fidelity simulation data, less accurate but much faster to compute, to improve the approximations with limited high-fidelity data. In this work, we introduce a novel data-driven framework of multi-fidelity surrogate modeling for parametrized, time-dependent problems using long short-term memory (LSTM) networks, to enhance output predictions both for unseen parameter values and forward in time simultaneously -- a task known to be particularly challenging for data-driven models. We demonstrate the wide applicability of the proposed approaches in a variety of engineering problems with high- and low-fidelity data generated through fine versus coarse meshes, small versus large time steps, or finite element full-order versus deep learning reduced-order models. Numerical results show that the proposed multi-fidelity LSTM networks not only improve single-fidelity regression significantly, but also outperform the multi-fidelity models based on feed-forward neural networks. |
-
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. |
|
