Giovanni Stabile

Giovanni Stabile, Gianluigi Rozza

In this work a stabilised and reduced Galerkin projection of the incompressible unsteady Navier-Stokes equations for moderate Reynolds number is presented. The full-order model, on which the Galerkin projection is applied, is based on a finite volumes approximation. The reduced basis spaces are constructed with a POD approach. Two different pressure stabilisation strategies are proposed and compared: the former one is based on the supremizer enrichment of the velocity space, and the latter one is based on a pressure Poisson equation approach.

Gianluigi Rozza, Haris Malik, Nicola Demo et al.

Reduced order modeling has gained considerable attention in recent decades owing to the advantages offered in reduced computational times and multiple solutions for parametric problems. The focus of this manuscript is the application of model order reduction techniques in various engineering and scientific applications including but not limited to mechanical, naval and aeronautical engineering. The focus here is kept limited to computational fluid mechanics and related applications. The advances in the reduced order modeling with proper orthogonal decomposition and reduced basis method are presented as well as a brief discussion of dynamic mode decomposition and also some present advances in the parameter space reduction. Here, an overview of the challenges faced and possible solutions are presented with examples from various problems.

Francesco Romor, Giovanni Stabile, Gianluigi Rozza

Non-affine parametric dependencies, nonlinearities and advection-dominated regimes of the model of interest can result in a slow Kolmogorov n-width decay, which precludes the realization of efficient reduced-order models based on linear subspace approximations. Among the possible solutions, there are purely data-driven methods that leverage autoencoders and their variants to learn a latent representation of the dynamical system, and then evolve it in time with another architecture. Despite their success in many applications where standard linear techniques fail, more has to be done to increase the interpretability of the results, especially outside the training range and not in regimes characterized by an abundance of data. Not to mention that none of the knowledge on the physics of the model is exploited during the predictive phase. In order to overcome these weaknesses, we implement the non-linear manifold method introduced by Carlberg et al [37] with hyper-reduction achieved through reduced over-collocation and teacher-student training of a reduced decoder. We test the methodology on a 2d non-linear conservation law and a 2d shallow water models, and compare the results obtained with a purely data-driven method for which the dynamics is evolved in time with a long-short term memory network.

E. N. Karatzas, G. Stabile, N. Atallah et al.

A model order reduction technique is combined with an embedded boundary finite element method with a POD-Galerkin strategy. The proposed methodology is applied to parametrized heat transfer problems and we rely on a sufficiently refined shape-regular background mesh to account for parametrized geometries. In particular, the employed embedded boundary element method is the Shifted Boundary Method (SBM) recently proposed. This approach is based on the idea of shifting the location of true boundary conditions to a surrogate boundary, with the goal of avoiding cut cells near the boundary of the computational domain. This combination of methodologies has multiple advantages. In the first place, since the Shifted Boundary Method always relies on the same background mesh, there is no need to update the discretized parametric domain. Secondly, we avoid the treatment of cut cell elements, which usually need particular attention. Thirdly, since the whole background mesh is considered in the reduced basis construction, the SBM allows for a smooth transition of the reduced modes across the immersed domain boundary. The performances of the method are verified in two dimensional heat transfer numerical examples.

Dario Coscia, Laura Meneghetti, Nicola Demo et al.

Convolutional Neural Network (CNN) is one of the most important architectures in deep learning. The fundamental building block of a CNN is a trainable filter, represented as a discrete grid, used to perform convolution on discrete input data. In this work, we propose a continuous version of a trainable convolutional filter able to work also with unstructured data. This new framework allows exploring CNNs beyond discrete domains, enlarging the usage of this important learning technique for many more complex problems. Our experiments show that the continuous filter can achieve a level of accuracy comparable to the state-of-the-art discrete filter, and that it can be used in current deep learning architectures as a building block to solve problems with unstructured domains as well.

Sajad Salavatidezfouli, Giovanni Stabile, Gianluigi Rozza

This research study explores the applicability of Deep Reinforcement Learning (DRL) for thermal control based on Computational Fluid Dynamics. To accomplish that, the forced convection on a hot plate prone to a pulsating cooling jet with variable velocity has been investigated. We begin with evaluating the efficiency and viability of a vanilla Deep Q-Network (DQN) method for thermal control. Subsequently, a comprehensive comparison between different variants of DRL is conducted. Soft Double and Duel DQN achieved better thermal control performance among all the variants due to their efficient learning and action prioritization capabilities. Results demonstrate that the soft Double DQN outperforms the hard Double DQN. Moreover, soft Double and Duel can maintain the temperature in the desired threshold for more than 98% of the control cycle. These findings demonstrate the promising potential of DRL in effectively addressing thermal control systems.

Isabella Carla Gonnella, Martin W. Hess, Giovanni Stabile et al.

Parametric time-dependent systems are of a crucial importance in modeling real phenomena, often characterized by non-linear behaviors too. Those solutions are typically difficult to generalize in a sufficiently wide parameter space while counting on limited computational resources available. As such, we present a general two-stages deep learning framework able to perform that generalization with low computational effort in time. It consists in a separated training of two pipe-lined predictive models. At first, a certain number of independent neural networks are trained with data-sets taken from different subsets of the parameter space. Successively, a second predictive model is specialized to properly combine the first-stage guesses and compute the right predictions. Promising results are obtained applying the framework to incompressible Navier-Stokes equations in a cavity (Rayleigh-Bernard cavity), obtaining a 97% reduction in the computational time comparing with its numerical resolution for a new value of the Grashof number.

Sajad Salavatidezfouli, Saeid Rakhsha, Armin Sheidani et al.

This paper aims to comprehensively investigate the efficacy of various Model Order Reduction (MOR) and deep learning techniques in predicting heat transfer in a pulsed jet impinging on a concave surface. Expanding on the previous experimental and numerical research involving pulsed circular jets, this investigation extends to evaluate Predictive Surrogate Models (PSM) for heat transfer across various jet characteristics. To this end, this work introduces two predictive approaches, one employing a Fast Fourier Transformation augmented Artificial Neural Network (FFT-ANN) for predicting the average Nusselt number under constant-frequency scenarios. Moreover, the investigation introduces the Proper Orthogonal Decomposition and Long Short-Term Memory (POD-LSTM) approach for random-frequency impingement jets. The POD-LSTM method proves to be a robust solution for predicting the local heat transfer rate under random-frequency impingement scenarios, capturing both the trend and value of temporal modes. The comparison of these approaches highlights the versatility and efficacy of advanced machine learning techniques in modelling complex heat transfer phenomena.

Davide Papapicco, Nicola Demo, Michele Girfoglio et al.

Models with dominant advection always posed a difficult challenge for projection-based reduced order modelling. Many methodologies that have recently been proposed are based on the pre-processing of the full-order solutions to accelerate the Kolmogorov N-width decay thereby obtaining smaller linear subspaces with improved accuracy. These methods however must rely on the knowledge of the characteristic speeds in phase space of the solution, limiting their range of applicability to problems with explicit functional form for the advection field. In this work we approach the problem of automatically detecting the correct pre-processing transformation in a statistical learning framework by implementing a deep-learning architecture. The purely data-driven method allowed us to generalise the existing approaches of linear subspace manipulation to non-linear hyperbolic problems with unknown advection fields. The proposed algorithm has been validated against simple test cases to benchmark its performances and later successfully applied to a multiphase simulation.

Matteo Zancanaro, Markus Mrosek, Giovanni Stabile et al.

Geometrically parametrized Partial Differential Equations are nowadays widely used in many different fields as, for example, shape optimization processes or patient specific surgery studies. The focus of this work is on some advances for this topic, capable of increasing the accuracy with respect to previous approaches while relying on a high cost-benefit ratio performance. The main scope of this paper is the introduction of a new technique mixing up a classical Galerkin-projection approach together with a data-driven method to obtain a versatile and accurate algorithm for the resolution of geometrically parametrized incompressible turbulent Navier-Stokes problems. The effectiveness of this procedure is demonstrated on two different test cases: a classical academic back step problem and a shape deformation Ahmed body application. The results show into details the properties of the architecture we developed while exposing possible future perspectives for this work.

Saray Busto, Giovanni Stabile, Gianluigi Rozza et al.

The purpose of this work is to introduce a novel POD-Galerkin strategy for the hybrid finite volume/finite element solver introduced in Bermúdez et al. 2014 and Busto et al. 2018. The interest is into the incompressible Navier-Stokes equations coupled with an additional transport equation. The full order model employed in this article makes use of staggered meshes. This feature will be conveyed to the reduced order model leading to the definition of reduced basis spaces in both meshes. The reduced order model presented herein accounts for velocity, pressure, and a transport-related variable. The pressure term at both the full order and the reduced order level is reconstructed making use of a projection method. More precisely, a Poisson equation for pressure is considered within the reduced order model. Results are verified against three-dimensional manufactured test cases. Moreover a modified version of the classical cavity test benchmark including the transport of a species is analysed.