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
-
18/2021 - 31/03/2021
Gigante, G.; Vergara, C.
On the choice of interface parameters in Robin-Robin loosely coupled schemes for fluid-structure interaction | Abstract | | We consider two loosely-coupled schemes for the solution of the fluid-structure interaction problem in presence of large added mass effect. In particular, we introduce the Robin-Robin and Robin-Neumann explicit schemes where suitable interface conditions of Robin type are used. For the estimate of interface Robin parameters which guarantee stability of the numerical solution, we propose to optimize the reduction factor of the corresponding strongly-coupled (implicit) scheme, by means of the Optimized Schwarz method. To check the suitabilty of our proposals, we show numerical results both in an ideal cylindrical domain and in a real human carotid. |
-
17/2021 - 31/03/2021
Chew, R.; Benacchio, T.; Hastermann, G.; Klein, R.
Balanced data assimilation with a blended numerical model | Abstract | | A challenge arising from the local Bayesian assimilation of data in an atmospheric flow simulation is the imbalances it may introduce. Fast-mode imbalances of the order of the slower dynamics can be negated by employing a blended numerical model with seamless access to the compressible and the soundproof pseudo-incompressible dynamics. Here, the blended modelling strategy by Benacchio et al. (2014) is upgraded in an advanced numerical framework and extended with a Bayesian local ensemble data assimilation method. Upon assimilation of data, the model configuration is switched to the pseudo-incompressible regime for one time-step. After that, the model configuration is switched back to the compressible model for the duration of the assimilation window. The switching between model regimes is repeated for each subsequent assimilation window. An improved blending strategy ensures that a single time-step in the pseudo-incompressible regime is sufficient to filter imbalances. This improvement is based on three innovations: (i) the association of pressure fields computed at different stages of the numerical integration with actual time levels; (ii) a conversion of pressure-related variables between the model regimes derived from low Mach number asymptotics; and (iii) a judicious selection of the pressure variables used in converting numerical model states when a switch of models occurs. Travelling vortex and bubble convection experiments show that the imbalance arising from assimilation of the momentum fields can be eliminated by using this blended model, thereby achieving balanced analysis fields. The leftover imbalance in the thermodynamics can be quantified by scale analysis. |
-
16/2021 - 31/03/2021
Salvador, M.; Dede', L.; Manzoni, A.
Non intrusive reduced order modeling of parametrized PDEs by kernel POD and neural networks | Abstract | | We propose a nonlinear reduced basis method for the efficient approximation of parametrized partial differential equations (PDEs), exploiting kernel proper orthogonal decomposition (KPOD) for the generation of a reduced-order space and neural networks for the evaluation of the reduced-order approximation. In particular, we use KPOD in place of the more classical POD, on a set of high-fidelity solutions of the problem at hand to extract a reduced basis. This method provides a more accurate approximation of the snapshots' set featuring a lower dimension, while maintaining the same efficiency as POD. A neural network (NN) is then used to find the coefficients of the reduced basis by following a supervised learning paradigm and shown to be effective in learning the map between the time/parameter values and the projection of the high-fidelity snapshots onto the reduced space. In this NN, both the number of hidden layers and the number of neurons vary according to the intrinsic dimension of the differential problem at hand and the size of the reduced space. This adaptively built NN attains good performances in both the learning and the testing phases. Our approach is then tested on two benchmark problems, a one-dimensional wave equation and a two-dimensional nonlinear lid-driven cavity problem. We finally compare the proposed KPOD-NN technique with a POD-NN strategy, showing that KPOD allows a reduction of the number of modes that must be retained to reach a given accuracy in the reduced basis approximation. For this reason, the NN built to find the coefficients of the KPOD expansion is smaller, easier and less computationally demanding to train than the one used in the POD-NN strategy. |
-
15/2021 - 10/03/2021
Fumagalli, A.; Patacchini, F.S.
Model adaptation in a discrete fracture network: existence of solutions and numerical strategies | Abstract | | Fractures are normally present in the underground and are, for some physical processes, of paramount importance. Their accurate description is fundamental to obtain reliable numerical outcomes useful, e.g., for energy management. Depending on the physical and geometrical properties of the fractures, fluid flow can behave differently, going from a slow Darcian regime to more complicated Brinkman or even Forchheimer regimes for high velocity. The main problem is to determine where in the fractures one regime is more adequate than others. In order to determine these low-speed and high-speed regions, this work proposes an adaptive strategy which is based on selecting the appropriate constitutive law linking velocity and pressure according to a threshold criterion on the magnitude of the fluid velocity itself. Both theoretical and numerical aspects are considered and investigated, showing the potentiality of the proposed approach. From the analytical viewpoint, we show existence of weak solutions to such model under reasonable hypotheses on the constitutive laws. To this end, we use a variational approach identifying solutions with minimizers of an underlying energy functional. From the numerical viewpoint, we propose a one-dimensional algorithm which tracks the interface between the low- and high-speed regions. By running numerical experiments using this algorithm, we illustrate some interesting behaviors of our adaptive model on a single fracture and small networks of intersecting fractures. |
-
14/2021 - 10/03/2021
Peli, R.; Menafoglio, A.; Cervino, M.; Dovera, L.; Secchi, P;
Physics-based Residual Kriging for dynamically evolving functional random fields | Abstract | | We present a novel approach named Physics-based Residual Kriging for the statistical prediction of spatially dependent functional data. It incorporates a physical model - expressed by a partial differential equation - within a Universal Kriging setting through a geostatistical modelization of the residuals with respect to the physical model.
The approach is extended to deal with sequential problems, where samples of functional data become available along consecutive time intervals, in a context where the physical and stochastic processes generating them evolve, as time intervals succeed one another.
An incremental modeling is used to account for both these dynamics and the misfit between previous predictions and actual observations.We apply Physics-based Residual Kriging to forecast production rates of wells operating in a mature reservoir according to a given drilling schedule. We evaluate the predictive errors of the method in two different case studies. The first deals with a single-phase reservoir where production is supported by fluid injection, while the second considers again a single-phase reservoir but the production is driven by rock compaction. |
-
13/2021 - 10/03/2021
Ferro, N.; Perotto, S.; Cangiani, A.
An anisotropic recovery-based error estimator for adaptive discontinuous Galerkin methods | Abstract | | We present a new recovery-based anisotropic error estimator for discontinuous Galerkin finite element approximations of advection-diffusion problems. We propose a metric-based algorithm for mesh adaptation which is driven by this error estimator. Numerical verification on several test cases, both in the steady and in the unsteady setting, shows the effectiveness of the algorithm in capturing the intrinsic directionalities of the solution. |
-
12/2021 - 02/03/2021
di Cristofaro, D.; Galimberti, C.; Bianchi, D.; Ferrante, R.; Ferro, N.; Mannisi, M.; Perotto, S
Adaptive topology optimization for innovative 3D printed metamaterials | Abstract | | An adaptive method for designing the infill pattern of 3D printed objects is proposed. In particular, new unit cells for metamaterials are designed in order to match prescribed mechanical specifications. To this aim, we resort to topology optimization at the microscale driven by an inverse homogenization to guarantee the desired properties at the macroscale. The whole procedure is additionally enriched with an anisotropic adaptive generation of the computational mesh. The proposed algorithm is first numerically verified both in a mono- and in a multi-objective context. Then, a mechanical validation and 3D manufacturing through fused-model-deposition are carried out to assess the feasibility of the proposed design workflow. |
-
11/2021 - 02/03/2021
Antonietti,P.F.; Manzini, G.; Mazzieri, I.; Scacchi, S.; Verani, M.
The conforming virtual element method for polyharmonic and elastodynamics problems: a review | Abstract | | In this paper we review recent results on the conforming virtual element approximation of polyharmonic and elastodynamics problems. The structure and the content of this review is motivated by three paradigmatic examples of applications: classical and anisotropic Cahn-Hilliard equation and phase field models for brittle fracture, that are briefly discussed in the first part of the paper. We present and discuss the mathematical details of the conforming virtual element approximation of linear polyharmonic problems, the classical Cahn-Hilliard equation and linear elastodynamics problems. |
|
