Showing 1–50 of 182 results for author: Rozza, G

Data-driven Discovery of Delay Differential Equations with Discrete Delays

Authors: Alessandro Pecile, Nicola Demo, Marco Tezzele, Gianluigi Rozza, Dimitri Breda

Abstract: The Sparse Identification of Nonlinear Dynamics (SINDy) framework is a robust method for identifying governing equations, successfully applied to ordinary, partial, and stochastic differential equations. In this work we extend SINDy to identify delay differential equations by using an augmented library that includes delayed samples and Bayesian optimization. To identify a possibly unknown delay we… ▽ More The Sparse Identification of Nonlinear Dynamics (SINDy) framework is a robust method for identifying governing equations, successfully applied to ordinary, partial, and stochastic differential equations. In this work we extend SINDy to identify delay differential equations by using an augmented library that includes delayed samples and Bayesian optimization. To identify a possibly unknown delay we minimize the reconstruction error over a set of candidates. The resulting methodology improves the overall performance by remarkably reducing the number of calls to SINDy with respect to a brute force approach. We also address a multivariate setting to identify multiple unknown delays and (non-multiplicative) parameters. Several numerical tests on delay differential equations with different long-term behavior, number of variables, delays, and parameters support the use of Bayesian optimization highlighting both the efficacy of the proposed methodology and its computational advantages. As a consequence, the class of discoverable models is significantly expanded. △ Less

Submitted 28 July, 2024; originally announced July 2024.

Non-intrusive model reduction of advection-dominated hyperbolic problems using neural network shift augmented manifold transformations

Authors: Harshith Gowrachari, Nicola Demo, Giovanni Stabile, Gianluigi Rozza

Abstract: Advection-dominated problems are commonly noticed in nature, engineering systems, and a wide range of industrial processes. For these problems, linear approximation methods (proper orthogonal decomposition and reduced basis method) are not suitable, as the Kolmogorov $n$-width decay is slow, leading to inefficient and inaccurate reduced order models. There are few non-linear approaches to accelera… ▽ More Advection-dominated problems are commonly noticed in nature, engineering systems, and a wide range of industrial processes. For these problems, linear approximation methods (proper orthogonal decomposition and reduced basis method) are not suitable, as the Kolmogorov $n$-width decay is slow, leading to inefficient and inaccurate reduced order models. There are few non-linear approaches to accelerate the Kolmogorov $n$-width decay. In this work, we use a neural-network shift augmented transformation technique, that employs automatic-shit detection and detects the optimal non-linear transformation of the full-order model solution manifold $\mathcal{M}$. We exploit a deep-learning framework to derive parameter-dependent bijective mapping between the manifold $\mathcal{M}$ and the transformed manifold $\tilde{\mathcal{M}}$. It consists of two neural networks, 1) ShiftNet, to employ automatic-shift detection by learning the shift-operator, which finds the optimal shifts for numerous snapshots of the full-order solution manifold, to accelerate the Kolmogorov $n$-width decay, and 2) InterpNet, which learns the reference configuration and can reconstruct the field values of the same, for each shifted grid distribution. We construct non-intrusive reduced order models on the resulting transformed linear subspaces and employ automatic-shift detection for predictions. We test our methodology on advection-dominated problems, such as 1D travelling waves, 2D isentropic convective vortex and 2D two-phase flow test cases. This work leads to the development of the complete NNsPOD-ROM algorithm for model reduction of advection-dominated problems, comprising both offline-online stages. △ Less

Submitted 25 July, 2024; originally announced July 2024.

Comments: 22 pages, 19 Figures

On the accuracy and efficiency of reduced order models: towards real-world applications

Authors: Pierfrancesco Siena, Paquale Claudio Africa, Michele Girfoglio, Gianluigi Rozza

Abstract: This chapter provides an extended overview about Reduced Order Models (ROMs), with a focus on their features in terms of efficiency and accuracy. In particular, the aim is to browse the more common ROM frameworks, considering both intrusive and data-driven approaches. We present the validation of such techniques against several test cases. The first one is an academic benchmark, the thermal block… ▽ More This chapter provides an extended overview about Reduced Order Models (ROMs), with a focus on their features in terms of efficiency and accuracy. In particular, the aim is to browse the more common ROM frameworks, considering both intrusive and data-driven approaches. We present the validation of such techniques against several test cases. The first one is an academic benchmark, the thermal block problem, where a Poisson equation is considered. Here a classic intrusive ROM framework based on a Galerkin projection scheme is employed. The second and third test cases come from real-world applications, the one related to the investigation of the blood flow patterns in a patient specific coronary arteries configuration where the Navier Stokes equations are addressed and the other one concerning the granulation process within pharmaceutical industry where a fluid-particle system is considered. Here we employ two data-driven ROM approaches showing a very relevant trade-off between accuracy and efficiency. In the last part of the contribution, two novel technological platforms, ARGOS and ATLAS, are presented. They are designed to provide a user-friendly access to data-driven models for real-time predictions for complex biomedical and industrial problems. △ Less

Submitted 30 April, 2024; originally announced July 2024.

A time-adaptive algorithm for pressure dominated flows: a heuristic estimator

Authors: Ivan Prusak, Davide Torlo, Monica Nonino, Gianluigi Rozza

Abstract: This work aims to introduce a heuristic timestep-adaptive algorithm for Computational Fluid Dynamics (CFD) and Fluid-Structure Interaction (FSI) problems where the flow is dominated by the pressure. In such scenarios, many time-adaptive algorithms based on the interplay of implicit and explicit time schemes fail to capture the fast transient dynamics of pressure fields. We present an algorithm tha… ▽ More This work aims to introduce a heuristic timestep-adaptive algorithm for Computational Fluid Dynamics (CFD) and Fluid-Structure Interaction (FSI) problems where the flow is dominated by the pressure. In such scenarios, many time-adaptive algorithms based on the interplay of implicit and explicit time schemes fail to capture the fast transient dynamics of pressure fields. We present an algorithm that relies on a temporal error estimator using Backward Differentiation Formulae (BDF$k$) of order $k=2,3$. Specifically, we demonstrate that the implicit BDF$3$ solution can be well approximated by applying a single Newton-type nonlinear solver correction to the implicit BDF$2$ solution. The difference between these solutions determines our adaptive temporal error estimator. The effectiveness of our approach is confirmed by numerical experiments conducted on a backward-facing step flow CFD test case with Reynolds number $300$ and on a two-dimensional haemodynamics FSI benchmark. △ Less

Submitted 29 June, 2024; originally announced July 2024.

A hybrid reduced-order model for segregated fluid-structure interaction solvers in an ALE approach at high Reynolds number

Authors: Valentin Nkana Ngan, Giovanni Stabile, Andrea Mola, Gianluigi Rozza

Abstract: This study introduces a first step for constructing a hybrid reduced-order models (ROMs) for segregated fluid-structure interaction in an Arbitrary Lagrangian-Eulerian (ALE) approach at a high Reynolds number using the Finite Volume Method (FVM). The ROM is driven by proper orthogonal decomposition (POD) with hybrid techniques that combines the classical Galerkin projection and two data-driven met… ▽ More This study introduces a first step for constructing a hybrid reduced-order models (ROMs) for segregated fluid-structure interaction in an Arbitrary Lagrangian-Eulerian (ALE) approach at a high Reynolds number using the Finite Volume Method (FVM). The ROM is driven by proper orthogonal decomposition (POD) with hybrid techniques that combines the classical Galerkin projection and two data-driven methods (radial basis networks , and neural networks/ long short term memory). Results demonstrate the ROM ability to accurately capture the physics of fluid-structure interaction phenomena. This approach is validated through a case study focusing on flow-induced vibration (FIV) of a pitch-plunge airfoil at a high Reynolds number 10000000. △ Less

Submitted 18 June, 2024; originally announced June 2024.

Parametric Intrusive Reduced Order Models enhanced with Machine Learning Correction Terms

Authors: Anna Ivagnes, Giovanni Stabile, Gianluigi Rozza

Abstract: In this paper, we propose an equation-based parametric Reduced Order Model (ROM), whose accuracy is improved with data-driven terms added into the reduced equations. These additions have the aim of reintroducing contributions that in standard ROMs are not taken into account. In particular, in this work we consider two types of contributions: the turbulence modeling, added through a reduced-order a… ▽ More In this paper, we propose an equation-based parametric Reduced Order Model (ROM), whose accuracy is improved with data-driven terms added into the reduced equations. These additions have the aim of reintroducing contributions that in standard ROMs are not taken into account. In particular, in this work we consider two types of contributions: the turbulence modeling, added through a reduced-order approximation of the eddy viscosity field, and the correction model, aimed to re-introduce the contribution of the discarded modes. Both approaches have been investigated in previous works and the goal of this paper is to extend the model to a parametric setting making use of ad-hoc machine learning procedures. More in detail, we investigate different neural networks' architectures, from simple dense feed-forward to Long-Short Term Memory neural networks, in order to find the most suitable model for the re-introduced contributions. We tested the methods on two test cases with different behaviors: the periodic turbulent flow past a circular cylinder and the unsteady turbulent flow in a channel-driven cavity. In both cases, the parameter considered is the Reynolds number and the machine learning-enhanced ROM considerably improved the pressure and velocity accuracy with respect to the standard ROM. △ Less

Submitted 6 June, 2024; originally announced June 2024.

A brief review of Reduced Order Models using intrusive and non-intrusive techniques

Authors: Guglielmo Padula, Michele Girfoglio, Gianluigi Rozza

Abstract: Reduced Order Models (ROMs) have gained a great attention by the scientific community in the last years thanks to their capabilities of significantly reducing the computational cost of the numerical simulations, which is a crucial objective in applications like real time control and shape optimization. This contribution aims to provide a brief overview about such a topic. We discuss both an intrus… ▽ More Reduced Order Models (ROMs) have gained a great attention by the scientific community in the last years thanks to their capabilities of significantly reducing the computational cost of the numerical simulations, which is a crucial objective in applications like real time control and shape optimization. This contribution aims to provide a brief overview about such a topic. We discuss both an intrusive framework based on a Galerkin projection technique and non-intrusive approaches, including Physics Informed Neural Networks (PINN), purely Data-Driven Neural Networks (DDNN), Radial Basis Functions (RBF), Dynamic Mode Decomposition (DMD) and Gaussian Process Regression (GPR). We also briefly mention geometrical parametrization and dimensionality reduction methods like Active Subspaces (AS). Then we present some results related to academic test cases as well as a preliminary investigation related to an industrial application. △ Less

Submitted 4 June, 2024; v1 submitted 1 June, 2024; originally announced June 2024.

Stabilized POD Reduced Order Models for convection-dominated incompressible flows

Authors: Pierfrancesco Siena, Michele Girfoglio, Annalisa Quaini, Gianluigi Rozza

Abstract: We present a comparative computational study of two stabilized Reduced Order Models (ROMs) for the simulation of convection-dominated incompressible flow (Reynolds number of the order of a few thousands). Representative solutions in the parameter space, which includes either time only or time and Reynolds number, are computed with a Finite Volume method and used to generate a reduced basis via Pro… ▽ More We present a comparative computational study of two stabilized Reduced Order Models (ROMs) for the simulation of convection-dominated incompressible flow (Reynolds number of the order of a few thousands). Representative solutions in the parameter space, which includes either time only or time and Reynolds number, are computed with a Finite Volume method and used to generate a reduced basis via Proper Orthogonal Decomposition (POD). Galerkin projection of the Navier-Stokes equations onto the reduced space is used to compute the ROM solution. To ensure computational efficiency, the number of POD modes is truncated and ROM solution accuracy is recovered through two stabilization methods: i) adding a global constant artificial viscosity to the reduced dimensional model, and ii) adding a different value of artificial viscosity for the different POD modes. We test the stabilized ROMs for fluid flow in an idealized medical device consisting of a conical convergent, a narrow throat, and a sudden expansion. Both stabilization methods significantly improve the ROM solution accuracy over a standard (non-stabilized) POD-Galerkin model. △ Less

Submitted 30 April, 2024; originally announced April 2024.

Computational study of numerical flux schemes for mesoscale atmospheric flows in a Finite Volume framework

Authors: Nicola Clinco, Michele Girfoglio, Annalisa Quaini, Gianluigi Rozza

Abstract: We develop, and implement in a Finite Volume environment, a density-based approach for the Euler equations written in conservative form using density, momentum, and total energy as variables. Under simplifying assumptions, these equations are used to describe non-hydrostatic atmospheric flow. The well-balancing of the approach is ensured by a local hydrostatic reconstruction updated in runtime dur… ▽ More We develop, and implement in a Finite Volume environment, a density-based approach for the Euler equations written in conservative form using density, momentum, and total energy as variables. Under simplifying assumptions, these equations are used to describe non-hydrostatic atmospheric flow. The well-balancing of the approach is ensured by a local hydrostatic reconstruction updated in runtime during the simulation to keep the numerical error under control. To approximate the solution of the Riemann problem, we consider four methods: Roe-Pike, HLLC, AUSM+-up and HLLC-AUSM. We assess our density-based approach and compare the accuracy of these four approximated Riemann solvers using two two classical benchmarks, namely the smooth rising thermal bubble and the density current. △ Less

Submitted 30 April, 2024; originally announced April 2024.

Comments: 22 pages

Linear and nonlinear filtering for a two-layer quasi-geostrophic ocean model

Authors: Lander Besabe, Michele Girfoglio, Annalisa Quaini, Gianluigi Rozza

Abstract: Although the two-layer quasi-geostrophic equations (2QGE) are a simplified model for the dynamics of a stratified, wind-driven ocean, their numerical simulation is still plagued by the need for high resolution to capture the full spectrum of turbulent scales. Since such high resolution would lead to unreasonable computational times, it is typical to resort to coarse low-resolution meshes combined… ▽ More Although the two-layer quasi-geostrophic equations (2QGE) are a simplified model for the dynamics of a stratified, wind-driven ocean, their numerical simulation is still plagued by the need for high resolution to capture the full spectrum of turbulent scales. Since such high resolution would lead to unreasonable computational times, it is typical to resort to coarse low-resolution meshes combined with the so-called eddy viscosity parameterization to account for the diffusion mechanisms that are not captured due to mesh under-resolution. We propose to enable the use of further coarsened meshes by adding a (linear or nonlinear) differential low-pass to the 2QGE, without changing the eddy viscosity coefficient. While the linear filter introduces constant (additional) artificial viscosity everywhere in the domain, the nonlinear filter relies on an indicator function to determine where and how much artificial viscosity is needed. Through several numerical results for a double-gyre wind forcing benchmark, we show that with the nonlinear filter we obtain accurate results with very coarse meshes, thereby drastically reducing the computational time (speed up ranging from 30 to 300). △ Less

Submitted 17 April, 2024; originally announced April 2024.

Comments: 34 pages, 13 figures

Kinetic data-driven approach to turbulence subgrid modeling

Authors: Giulio Ortali, Alessandro Gabbana, Nicola Demo, Gianluigi Rozza, Federico Toschi

Abstract: Numerical simulations of turbulent flows are well known to pose extreme computational challenges due to the huge number of dynamical degrees of freedom required to correctly describe the complex multi-scale statistical correlations of the velocity. On the other hand, kinetic mesoscale approaches based on the Boltzmann equation, have the potential to describe a broad range of flows, stretching well… ▽ More Numerical simulations of turbulent flows are well known to pose extreme computational challenges due to the huge number of dynamical degrees of freedom required to correctly describe the complex multi-scale statistical correlations of the velocity. On the other hand, kinetic mesoscale approaches based on the Boltzmann equation, have the potential to describe a broad range of flows, stretching well beyond the special case of gases close to equilibrium, which results in the ordinary Navier-Stokes dynamics. Here we demonstrate that, by properly tuning, a kinetic approach can statistically reproduce the quantitative dynamics of the larger scales in turbulence, thereby providing an alternative, computationally efficient and physically rooted approach towards subgrid scale (SGS) modeling in turbulence. More specifically we show that by leveraging on data from fully resolved Direct Numerical Simulation (DNS) data we can learn a collision operator for the discretized Boltzmann equation solver (the lattice Boltzmann method), which effectively implies a turbulence subgrid closure model. The mesoscopic nature of our formulation makes the learning problem fully local in both space and time, leading to reduced computational costs and enhanced generalization capabilities. We show that the model offers superior performance compared to traditional methods, such as the Smagorinsky model, being less dissipative and, therefore, being able to more closely capture the intermittency of higher-order velocity correlations. △ Less

Submitted 27 March, 2024; originally announced March 2024.

A LSTM-enhanced surrogate model to simulate the dynamics of particle-laden fluid systems

Authors: Arash Hajisharifi, Rahul Halder, Michele Girfoglio, Andrea Beccari, Domenico Bonanni, Gianluigi Rozza

Abstract: The numerical treatment of fluid-particle systems is a very challenging problem because of the complex coupling phenomena occurring between the two phases. Although accurate mathematical modelling is available to address this kind of application, the computational cost of the numerical simulations is very expensive. The use of the most modern high-performance computing infrastructures could help t… ▽ More The numerical treatment of fluid-particle systems is a very challenging problem because of the complex coupling phenomena occurring between the two phases. Although accurate mathematical modelling is available to address this kind of application, the computational cost of the numerical simulations is very expensive. The use of the most modern high-performance computing infrastructures could help to mitigate such an issue but not completely fix it. In this work, we develop a non-intrusive data-driven reduced order model (ROM) for Computational Fluid Dynamics (CFD) - Discrete Element Method (DEM) simulations. The ROM is built using the proper orthogonal decomposition (POD) for the computation of the reduced basis space and the Long Short-Term Memory (LSTM) network for the computation of the reduced coefficients. We are interested in dealing both with system identification and prediction. The most relevant novelties rely on (i) a filtering procedure of the full-order snapshots to reduce the dimensionality of the reduced problem and (ii) a preliminary treatment of the particle phase. The accuracy of our ROM approach is assessed against the classic Goldschmidt fluidized bed benchmark problem. Finally, we also provide some insights about the efficiency of our ROM approach. △ Less

Submitted 21 March, 2024; originally announced March 2024.

Comments: 13 Figures

Enhancing non-intrusive Reduced Order Models with space-dependent aggregation methods

Authors: Anna Ivagnes, Niccolò Tonicello, Paola Cinnella, Gianluigi Rozza

Abstract: In this manuscript, we combine non-intrusive reduced order models (ROMs) with space-dependent aggregation techniques to build a mixed-ROM. The prediction of the mixed formulation is given by a convex linear combination of the predictions of some previously-trained ROMs, where we assign to each model a space-dependent weight. The ROMs taken into account to build the mixed model exploit different re… ▽ More In this manuscript, we combine non-intrusive reduced order models (ROMs) with space-dependent aggregation techniques to build a mixed-ROM. The prediction of the mixed formulation is given by a convex linear combination of the predictions of some previously-trained ROMs, where we assign to each model a space-dependent weight. The ROMs taken into account to build the mixed model exploit different reduction techniques, such as Proper Orthogonal Decomposition (POD) and AutoEncoders (AE), and/or different approximation techniques, namely a Radial Basis Function Interpolation (RBF), a Gaussian Process Regression (GPR) or a feed-forward Artificial Neural Network (ANN). The contribution of each model is retained with higher weights in the regions where the model performs best, and, vice versa, with smaller weights where the model has a lower accuracy with respect to the other models. Finally, a regression technique, namely a Random Forest, is exploited to evaluate the weights for unseen conditions. The performance of the aggregated model is evaluated on two different test cases: the 2D flow past a NACA 4412 airfoil, with an angle of attack of 5 degrees, having as parameter the Reynolds number varying between 1e5 and 1e6 and a transonic flow over a NACA 0012 airfoil, considering as parameter the angle of attack. In both cases, the mixed-ROM has provided improved accuracy with respect to each individual ROM technique. △ Less

Submitted 8 March, 2024; originally announced March 2024.

Optimized Bayesian Framework for Inverse Heat Transfer Problems Using Reduced Order Methods

Authors: Kabir Bakhshaei, Umberto Emil Morelli, Giovanni Stabile, Gianluigi Rozza

Abstract: A stochastic inverse heat transfer problem is formulated to infer the transient heat flux, treated as an unknown Neumann boundary condition. Therefore, an Ensemble-based Simultaneous Input and State Filtering as a Data Assimilation technique is utilized for simultaneous temperature distribution prediction and heat flux estimation. This approach is incorporated with Radial Basis Functions not only… ▽ More A stochastic inverse heat transfer problem is formulated to infer the transient heat flux, treated as an unknown Neumann boundary condition. Therefore, an Ensemble-based Simultaneous Input and State Filtering as a Data Assimilation technique is utilized for simultaneous temperature distribution prediction and heat flux estimation. This approach is incorporated with Radial Basis Functions not only to lessen the size of unknown inputs but also to mitigate the computational burden of this technique. The procedure applies to the specific case of a mold used in Continuous Casting machinery, and it is based on the sequential availability of temperature provided by thermocouples inside the mold. Our research represents a significant contribution to achieving probabilistic boundary condition estimation in real-time handling with noisy measurements and errors in the model. We additionally demonstrate the procedure's dependence on some hyperparameters that are not documented in the existing literature. Accurate real-time prediction of the heat flux is imperative for the smooth operation of Continuous Casting machinery at the boundary region where the Continuous Casting mold and the molten steel meet which is not also physically measurable. Thus, this paves the way for efficient real-time monitoring and control, which is critical for preventing caster shutdowns. △ Less

Submitted 29 February, 2024; originally announced February 2024.

A comparative computational study of different formulations of the compressible Euler equations for mesoscale atmospheric flows in a finite volume framework

Authors: Michele Girfoglio, Annalisa Quaini, Gianluigi Rozza

Abstract: We consider three conservative forms of the mildly compressible Euler equations, called CE1, CE2 and CE3, with the goal of understanding which leads to the most accurate and robust pressure-based solver in a finite volume environment. Forms CE1 and CE2 are both written in density, momentum, and specific enthalpy, but employ two different treatments of the buoyancy and pressure gradient terms: for… ▽ More We consider three conservative forms of the mildly compressible Euler equations, called CE1, CE2 and CE3, with the goal of understanding which leads to the most accurate and robust pressure-based solver in a finite volume environment. Forms CE1 and CE2 are both written in density, momentum, and specific enthalpy, but employ two different treatments of the buoyancy and pressure gradient terms: for CE1 it is the standard pressure splitting implemented in open-source finite volume solvers (e.g., OpenFOAM), while for CE2 it is the typical pressure splitting found in computational atmospheric studies. Form CE3 is written in density, momentum, and potential temperature, with the buoyancy and pressure terms addressed as in CE2. For each formulation, we adopt a computationally efficient splitting approach. The three formulations are thoroughly assessed and compared through six benchmark tests involving dry air flow over a flat terrain or orography. We found that all three models are able to provide accurate results for the tests with a flat terrain, although the solvers based on the CE2 and CE3 forms are more robust. As for the mountain tests, CE1 solutions become unstable, while the CE2 and CE3 models provide results in very good agreement with data in the literature, the CE3 model being the most accurate. Hence, the CE3 model is the most accurate, reliable, and robust for the simulation of mesoscale atmospheric flows when using a pressure-based approach and space discretization by a finite volume method. △ Less

Submitted 28 February, 2024; originally announced February 2024.

Comments: 23 pages, 15 figures, 3 tables

A stochastic perturbation approach to nonlinear bifurcating problems

Authors: Isabella Carla Gonnella, Moaad Khamlich, Federico Pichi, Gianluigi Rozza

Abstract: Incorporating probabilistic terms in mathematical models is crucial for capturing and quantifying uncertainties of real-world systems. However, stochastic models typically require large computational resources to produce meaningful statistics. For such reason, the development of reduction techniques becomes essential for enabling efficient and scalable simulations of complex scenarios while quanti… ▽ More Incorporating probabilistic terms in mathematical models is crucial for capturing and quantifying uncertainties of real-world systems. However, stochastic models typically require large computational resources to produce meaningful statistics. For such reason, the development of reduction techniques becomes essential for enabling efficient and scalable simulations of complex scenarios while quantifying the underlying uncertainties. In this work, we study the accuracy of Polynomial Chaos (PC) surrogate expansion of the probability space on a bifurcating phenomena in fluid dynamics, namely the Coandă effect. In particular, we propose a novel non-deterministic approach to generic bifurcation problems, where the stochastic setting gives a different perspective on the non-uniqueness of the solution, also avoiding expensive simulations for many instances of the parameter. Thus, starting from the formulation of the Spectral Stochastic Finite Element Method (SSFEM), we extend the methodology to deal with solutions of a bifurcating problem, by working with a perturbed version of the deterministic model. We discuss the link between the deterministic and the stochastic bifurcation diagram, highlighting the surprising capability of PC polynomials coefficients of giving insights on the deterministic solution manifold. △ Less

Submitted 30 July, 2024; v1 submitted 26 February, 2024; originally announced February 2024.

A Predictive Surrogate Model for Heat Transfer of an Impinging Jet on a Concave Surface

Authors: Sajad Salavatidezfouli, Saeid Rakhsha, Armin Sheidani, Giovanni Stabile, Gianluigi Rozza

Abstract: 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 vari… ▽ More 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. △ Less

Submitted 16 February, 2024; originally announced February 2024.

Optimisation-Based Coupling of Finite Element Model and Reduced Order Model for Computational Fluid Dynamics

Authors: Ivan Prusak, Davide Torlo, Monica Nonino, Gianluigi Rozza

Abstract: Using Domain Decomposition (DD) algorithm on non--overlapping domains, we compare couplings of different discretisation models, such as Finite Element (FEM) and Reduced Order (ROM) models for separate subcomponents. In particular, we consider an optimisation-based DD model where the coupling on the interface is performed using a control variable representing the normal flux. We use iterative gradi… ▽ More Using Domain Decomposition (DD) algorithm on non--overlapping domains, we compare couplings of different discretisation models, such as Finite Element (FEM) and Reduced Order (ROM) models for separate subcomponents. In particular, we consider an optimisation-based DD model where the coupling on the interface is performed using a control variable representing the normal flux. We use iterative gradient-based optimisation algorithms to decouple the subdomain state solutions as well as to locally generate ROMs on each subdomain. Then, we consider FEM or ROM discretisation models for each of the DD problem components, namely, the triplet state1-state2-control. On the backward-facing step Navier-Stokes (NS) problem, we investigate the efficacy of the presented couplings in terms of optimisation iterations, optimal functional values and relative errors. △ Less

Submitted 10 July, 2024; v1 submitted 16 February, 2024; originally announced February 2024.

Comments: Revised version

PyDMD: A Python package for robust dynamic mode decomposition

Authors: Sara M. Ichinaga, Francesco Andreuzzi, Nicola Demo, Marco Tezzele, Karl Lapo, Gianluigi Rozza, Steven L. Brunton, J. Nathan Kutz

Abstract: The dynamic mode decomposition (DMD) is a simple and powerful data-driven modeling technique that is capable of revealing coherent spatiotemporal patterns from data. The method's linear algebra-based formulation additionally allows for a variety of optimizations and extensions that make the algorithm practical and viable for real-world data analysis. As a result, DMD has grown to become a leading… ▽ More The dynamic mode decomposition (DMD) is a simple and powerful data-driven modeling technique that is capable of revealing coherent spatiotemporal patterns from data. The method's linear algebra-based formulation additionally allows for a variety of optimizations and extensions that make the algorithm practical and viable for real-world data analysis. As a result, DMD has grown to become a leading method for dynamical system analysis across multiple scientific disciplines. PyDMD is a Python package that implements DMD and several of its major variants. In this work, we expand the PyDMD package to include a number of cutting-edge DMD methods and tools specifically designed to handle dynamics that are noisy, multiscale, parameterized, prohibitively high-dimensional, or even strongly nonlinear. We provide a complete overview of the features available in PyDMD as of version 1.0, along with a brief overview of the theory behind the DMD algorithm, information for developers, tips regarding practical DMD usage, and introductory coding examples. All code is available at https://github.com/PyDMD/PyDMD . △ Less

Submitted 12 February, 2024; originally announced February 2024.

Physics Informed Neural Network Framework for Unsteady Discretized Reduced Order System

Authors: Rahul Halder, Giovanni Stabile, Gianluigi Rozza

Abstract: This work addresses the development of a physics-informed neural network (PINN) with a loss term derived from a discretized time-dependent reduced-order system. In this work, first, the governing equations are discretized using a finite difference scheme (whereas, any other discretization technique can be adopted), then projected on a reduced or latent space using the Proper Orthogonal Decompositi… ▽ More This work addresses the development of a physics-informed neural network (PINN) with a loss term derived from a discretized time-dependent reduced-order system. In this work, first, the governing equations are discretized using a finite difference scheme (whereas, any other discretization technique can be adopted), then projected on a reduced or latent space using the Proper Orthogonal Decomposition (POD)-Galerkin approach and next, the residual arising from discretized reduced order equation is considered as an additional loss penalty term alongside the data-driven loss term using different variants of deep learning method such as Artificial neural network (ANN), Long Short-Term Memory based neural network (LSTM). The LSTM neural network has been proven to be very effective for time-dependent problems in a purely data-driven environment. The current work demonstrates the LSTM network's potential over ANN networks in physics-informed neural networks (PINN) as well. The potential of using discretized governing equations instead of continuous form lies in the flexibility of input to the PINN. Different sizes of data ranging from small, medium to big datasets are used to assess the potential of discretized-physics-informed neural networks when there is very sparse or no data available. The proposed methods are applied to a pitch-plunge airfoil motion governed by rigid-body dynamics and a one-dimensional viscous Burgers' equation. The current work also demonstrates the prediction capability of various discretized-physics-informed neural networks outside the domain where the data is available or governing equation-based residuals are minimized. △ Less

Submitted 29 January, 2024; v1 submitted 23 November, 2023; originally announced November 2023.

Comments: 22 pages, 14 figures

Modal Analysis of the Wake Shed Behind a Horizontal Axis Wind Turbine with Flexible Blades

Authors: Sajad Salavatidezfouli, Armin Sheidani, Kabir Bakhshaei, Ali Safari, Arash Hajisharifi, Giovanni Stabile, Gianluigi Rozza

Abstract: The proper orthogonal decomposition has been applied on a full-scale horizontal-axis wind turbine to shed light on the wake characteristics behind the wind turbine. In reality, the blade tip experiences high deflections even at the rated conditions which definitely alter the wake flow field, and in the case of a wind farm, may complicate the inlet conditions of the downstream wind turbine. The tur… ▽ More The proper orthogonal decomposition has been applied on a full-scale horizontal-axis wind turbine to shed light on the wake characteristics behind the wind turbine. In reality, the blade tip experiences high deflections even at the rated conditions which definitely alter the wake flow field, and in the case of a wind farm, may complicate the inlet conditions of the downstream wind turbine. The turbine under consideration is the full-scale model of the National Renewable Energy Laboratory 5MW onshore wind turbine which is accompanied by several simulation complexities including turbulence, mesh motion and fluid-structure interaction. Results indicated an almost similar modal behaviour for the rigid and flexible turbines at the wake region. In addition, more flow structures in terms of local vortices and fluctuating velocity fields take place at the far wake region. The flow structures due to the wake shed from the tower tend to move towards the center and merge with that of the nacelle leading to an integral vortical structure 2.5 diameter away from the rotor. Also, it is concluded that the exclusion of the tower leads to missing a major part of the wake structures, especially at far-wake positions. △ Less

Submitted 8 May, 2024; v1 submitted 14 November, 2023; originally announced November 2023.

Deep Reinforcement Learning for the Heat Transfer Control of Pulsating Impinging Jets

Authors: Sajad Salavatidezfouli, Giovanni Stabile, Gianluigi Rozza

Abstract: 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. Subse… ▽ More 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. △ Less

Submitted 25 September, 2023; originally announced September 2023.

Journal ref: Advances in Computational Science and Engineering 2024

A Reduced Order Model formulation for left atrium flow: an Atrial Fibrillation case

Authors: Caterina Balzotti, Pierfrancesco Siena, Michele Girfoglio, Giovanni Stabile, Jorge Dueñas-Pamplona, José Sierra-Pallares, Ignacio Amat-Santos, Gianluigi Rozza

Abstract: A data-driven Reduced Order Model (ROM) based on a Proper Orthogonal Decomposition - Radial Basis Function (POD-RBF) approach is adopted in this paper for the analysis of blood flow dynamics in a patient-specific case of Atrial Fibrillation (AF). The Full Order Model (FOM) is represented by incompressible Navier-Stokes equations, discretized with a Finite Volume (FV) approach. Both the Newtonian a… ▽ More A data-driven Reduced Order Model (ROM) based on a Proper Orthogonal Decomposition - Radial Basis Function (POD-RBF) approach is adopted in this paper for the analysis of blood flow dynamics in a patient-specific case of Atrial Fibrillation (AF). The Full Order Model (FOM) is represented by incompressible Navier-Stokes equations, discretized with a Finite Volume (FV) approach. Both the Newtonian and the Casson's constitutive laws are employed. The aim is to build a computational tool able to efficiently and accurately reconstruct the patterns of relevant hemodynamics indices related to the stasis of the blood in a physical parametrization framework including the cardiac output in the Newtonian case and also the plasma viscosity and the hematocrit in the non-Newtonian one. Many FOM-ROM comparisons are shown to analyze the performance of our approach as regards errors and computational speed-up. △ Less

Submitted 19 September, 2023; originally announced September 2023.

Comments: 21 pages, 14 figures

Optimal Transport-inspired Deep Learning Framework for Slow-Decaying Problems: Exploiting Sinkhorn Loss and Wasserstein Kernel

Authors: Moaad Khamlich, Federico Pichi, Gianluigi Rozza

Abstract: Reduced order models (ROMs) are widely used in scientific computing to tackle high-dimensional systems. However, traditional ROM methods may only partially capture the intrinsic geometric characteristics of the data. These characteristics encompass the underlying structure, relationships, and essential features crucial for accurate modeling. To overcome this limitation, we propose a novel ROM fr… ▽ More Reduced order models (ROMs) are widely used in scientific computing to tackle high-dimensional systems. However, traditional ROM methods may only partially capture the intrinsic geometric characteristics of the data. These characteristics encompass the underlying structure, relationships, and essential features crucial for accurate modeling. To overcome this limitation, we propose a novel ROM framework that integrates optimal transport (OT) theory and neural network-based methods. Specifically, we investigate the Kernel Proper Orthogonal Decomposition (kPOD) method exploiting the Wasserstein distance as the custom kernel, and we efficiently train the resulting neural network (NN) employing the Sinkhorn algorithm. By leveraging an OT-based nonlinear reduction, the presented framework can capture the geometric structure of the data, which is crucial for accurate learning of the reduced solution manifold. When compared with traditional metrics such as mean squared error or cross-entropy, exploiting the Sinkhorn divergence as the loss function enhances stability during training, robustness against overfitting and noise, and accelerates convergence. To showcase the approach's effectiveness, we conduct experiments on a set of challenging test cases exhibiting a slow decay of the Kolmogorov n-width. The results show that our framework outperforms traditional ROM methods in terms of accuracy and computational efficiency. △ Less

Submitted 26 August, 2023; originally announced August 2023.

MSC Class: 68T05; 65D99; 41A05; 65N30; 41A46; 90C25

Generative Models for the Deformation of Industrial Shapes with Linear Geometric Constraints: model order and parameter space reductions

Authors: Guglielmo Padula, Francesco Romor, Giovanni Stabile, Gianluigi Rozza

Abstract: Real-world applications of computational fluid dynamics often involve the evaluation of quantities of interest for several distinct geometries that define the computational domain or are embedded inside it. For example, design optimization studies require the realization of response surfaces from the parameters that determine the geometrical deformations to relevant outputs to be optimized. In thi… ▽ More Real-world applications of computational fluid dynamics often involve the evaluation of quantities of interest for several distinct geometries that define the computational domain or are embedded inside it. For example, design optimization studies require the realization of response surfaces from the parameters that determine the geometrical deformations to relevant outputs to be optimized. In this context, a crucial aspect to be addressed are the limited resources at disposal to computationally generate different geometries or to physically obtain them from direct measurements. This is the case for patient-specific biomedical applications for example. When additional linear geometrical constraints need to be imposed, the computational costs increase substantially. Such constraints include total volume conservation, barycenter location and fixed moments of inertia. We develop a new paradigm that employs generative models from machine learning to efficiently sample new geometries with linear constraints. A consequence of our approach is the reduction of the parameter space from the original geometrical parametrization to a low-dimensional latent space of the generative models. Crucial is the assessment of the quality of the distribution of the constrained geometries obtained with respect to physical and geometrical quantities of interest. Non-intrusive model order reduction is enhanced since smaller parametric spaces are considered. We test our methodology on two academic test cases: a mixed Poisson problem on the 3d Stanford bunny with fixed barycenter deformations and the multiphase turbulent incompressible Navier-Stokes equations for the Duisburg test case with fixed volume deformations of the naval hull. △ Less

Submitted 7 August, 2023; originally announced August 2023.

Explicable hyper-reduced order models on nonlinearly approximated solution manifolds of compressible and incompressible Navier-Stokes equations

Authors: Francesco Romor, Giovanni Stabile, Gianluigi Rozza

Abstract: A slow decaying Kolmogorov n-width of the solution manifold of a parametric partial differential equation precludes the realization of efficient linear projection-based reduced-order models. This is due to the high dimensionality of the reduced space needed to approximate with sufficient accuracy the solution manifold. To solve this problem, neural networks, in the form of different architectures,… ▽ More A slow decaying Kolmogorov n-width of the solution manifold of a parametric partial differential equation precludes the realization of efficient linear projection-based reduced-order models. This is due to the high dimensionality of the reduced space needed to approximate with sufficient accuracy the solution manifold. To solve this problem, neural networks, in the form of different architectures, have been employed to build accurate nonlinear regressions of the solution manifolds. However, the majority of the implementations are non-intrusive black-box surrogate models, and only a part of them perform dimension reduction from the number of degrees of freedom of the discretized parametric models to a latent dimension. We present a new intrusive and explicable methodology for reduced-order modelling that employs neural networks for solution manifold approximation but that does not discard the physical and numerical models underneath in the predictive/online stage. We will focus on autoencoders used to compress further the dimensionality of linear approximants of solution manifolds, achieving in the end a nonlinear dimension reduction. After having obtained an accurate nonlinear approximant, we seek for the solutions on the latent manifold with the residual-based nonlinear least-squares Petrov-Galerkin method, opportunely hyper-reduced in order to be independent from the number of degrees of freedom. New adaptive hyper-reduction strategies are developed along with the employment of local nonlinear approximants. We test our methodology on two nonlinear time-dependent parametric benchmarks involving a supersonic flow past a NACA airfoil with changing Mach number and an incompressible turbulent flow around the Ahmed body with changing slant angle. △ Less

Submitted 7 August, 2023; originally announced August 2023.

Friedrichs' systems discretized with the Discontinuous Galerkin method: domain decomposable model order reduction and Graph Neural Networks approximating vanishing viscosity solutions

Authors: Francesco Romor, Davide Torlo, Gianluigi Rozza

Abstract: Friedrichs' systems (FS) are symmetric positive linear systems of first-order partial differential equations (PDEs), which provide a unified framework for describing various elliptic, parabolic and hyperbolic semi-linear PDEs such as the linearized Euler equations of gas dynamics, the equations of compressible linear elasticity and the Dirac-Klein-Gordon system. FS were studied to approximate PDEs… ▽ More Friedrichs' systems (FS) are symmetric positive linear systems of first-order partial differential equations (PDEs), which provide a unified framework for describing various elliptic, parabolic and hyperbolic semi-linear PDEs such as the linearized Euler equations of gas dynamics, the equations of compressible linear elasticity and the Dirac-Klein-Gordon system. FS were studied to approximate PDEs of mixed elliptic and hyperbolic type in the same domain. For this and other reasons, the versatility of the discontinuous Galerkin method (DGM) represents the best approximation space for FS. We implement a distributed memory solver for stationary FS in deal.II. Our focus is model order reduction. Since FS model hyperbolic PDEs, they often suffer from a slow Kolmogorov n-width decay. We develop two approaches to tackle this problem. The first is domain decomposable reduced-order models (DD-ROMs). We will show that the DGM offers a natural formulation of DD-ROMs, in particular regarding interface penalties, compared to the continuous finite element method. We also develop new repartitioning strategies to obtain more efficient local approximations of the solution manifold. The second approach involves graph neural networks used to infer the limit of a succession of projection-based linear ROMs corresponding to lower viscosity constants: the heuristic behind is to develop a multi-fidelity super-resolution paradigm to mimic the mathematical convergence to vanishing viscosity solutions while exploiting to the most interpretable and certified projection-based ROMs. △ Less

Submitted 7 August, 2023; originally announced August 2023.

An optimisation-based domain-decomposition reduced order model for parameter-dependent non-stationary fluid dynamics problems

Authors: Ivan Prusak, Davide Torlo, Monica Nonino, Gianluigi Rozza

Abstract: In this work, we address parametric non-stationary fluid dynamics problems within a model order reduction setting based on domain decomposition. Starting from the optimisation-based domain decomposition approach, we derive an optimal control problem, for which we present a convergence analysis in the case of non-stationary incompressible Navier-Stokes equations. We discretize the problem with the… ▽ More In this work, we address parametric non-stationary fluid dynamics problems within a model order reduction setting based on domain decomposition. Starting from the optimisation-based domain decomposition approach, we derive an optimal control problem, for which we present a convergence analysis in the case of non-stationary incompressible Navier-Stokes equations. We discretize the problem with the finite element method and we compare different model order reduction techniques: POD-Galerkin and a non-intrusive neural network procedure. We show that the classical POD-Galerkin is more robust and accurate also in transient areas, while the neural network can obtain simulations very quickly though being less precise in the presence of discontinuities in time or parameter domain. We test the proposed methodologies on two fluid dynamics benchmarks with physical parameters and time dependency: the non-stationary backward-facing step and lid-driven cavity flow. △ Less

Submitted 19 February, 2024; v1 submitted 3 August, 2023; originally announced August 2023.

Comments: arXiv admin note: substantial text overlap with arXiv:2211.14528

Mathematical modelling and computational reduction of molten glass fluid flow in a furnace melting basin

Authors: Francesco Ballarin, Enrique Delgado Ávila, Andrea Mola, Gianluigi Rozza

Abstract: In this work, we present the modelling and numerical simulation of a molten glass fluid flow in a furnace melting basin. We first derive a model for a molten glass fluid flow and present numerical simulations based on the Finite Element Method (FEM). We further discuss and validate the results obtained from the simulations by comparing them with experimental results. Finally, we also present a non… ▽ More In this work, we present the modelling and numerical simulation of a molten glass fluid flow in a furnace melting basin. We first derive a model for a molten glass fluid flow and present numerical simulations based on the Finite Element Method (FEM). We further discuss and validate the results obtained from the simulations by comparing them with experimental results. Finally, we also present a non-intrusive Proper Orthogonal Decomposition (POD) based on Artificial Neural Networks (ANN) to efficiently handle scenarios which require multiple simulations of the fluid flow upon changing parameters of relevant industrial interest. This approach lets us obtain solutions of a complex 3D model, with good accuracy with respect to the FEM solution, yet with negligible associated computational times. △ Less

Submitted 23 July, 2024; v1 submitted 27 July, 2023; originally announced July 2023.

A Comparison of Data-Driven Reduced Order Models for the Simulation of Mesoscale Atmospheric Flow

Authors: Arash Hajisharifi, Michele Girfoglio, Annalisa Quaini, Gianluigi Rozza

Abstract: The simulation of atmospheric flows by means of traditional discretization methods remains computationally intensive, hindering the achievement of high forecasting accuracy in short time frames. In this paper, we apply three reduced order models that have successfully reduced the computational time for different applications in computational fluid dynamics while preserving accuracy: Dynamic Mode D… ▽ More The simulation of atmospheric flows by means of traditional discretization methods remains computationally intensive, hindering the achievement of high forecasting accuracy in short time frames. In this paper, we apply three reduced order models that have successfully reduced the computational time for different applications in computational fluid dynamics while preserving accuracy: Dynamic Mode Decomposition (DMD), Hankel Dynamic Mode Decomposition (HDMD), and Proper Orthogonal Decomposition with Interpolation (PODI). The three methods are compared in terms of computational time and accuracy in the simulation of two well-known benchmarks for mesoscale flow. The accuracy of the DMD and HDMD solutions deteriorates rather quickly as the forecast time window expands, although these methods are designed to predict the dynamics of a system. The reason is likely the strong nonlinearity in the benchmark flows. The PODI solution is accurate for the entire duration of the time interval of interest thanks to the use of interpolation with radial basis functions. This holds true also when the model features a physical parameter expected to vary in a given range, as is typically the case in weather prediction. △ Less

Submitted 17 July, 2023; originally announced July 2023.

Geometrically Parametrised Reduced Order Models for the Study of Hysteresis of the Coanda Effect in Finite Element-based Incompressible Fluid Dynamics

Authors: J. R. Bravo, G. Stabile, M. Hess, J. A. Hernandez, R. Rossi, G. Rozza

Abstract: This article presents a general reduced order model (ROM) framework for addressing fluid dynamics problems involving time-dependent geometric parametrisations. The framework integrates Proper Orthogonal Decomposition (POD) and Empirical Cubature Method (ECM) hyper-reduction techniques to effectively approximate incompressible computational fluid dynamics simulations. To demonstrate the applicabili… ▽ More This article presents a general reduced order model (ROM) framework for addressing fluid dynamics problems involving time-dependent geometric parametrisations. The framework integrates Proper Orthogonal Decomposition (POD) and Empirical Cubature Method (ECM) hyper-reduction techniques to effectively approximate incompressible computational fluid dynamics simulations. To demonstrate the applicability of this framework, we investigate the behavior of a planar contraction-expansion channel geometry exhibiting bifurcating solutions known as the Coanda effect. By introducing time-dependent deformations to the channel geometry, we observe hysteresis phenomena in the solution. The paper provides a detailed formulation of the framework, including the stabilised finite elements full order model (FOM) and ROM, with a particular focus on the considerations related to geometric parametrisation. Subsequently, we present the results obtained from the simulations, analysing the solution behavior in a phase-space for the fluid velocity at a probe point, considered as the Quantity of Interest (QoI). Through qualitative and quantitative evaluations of the ROMs and hyper-reduced order models (HROMs), we demonstrate their ability to accurately reproduce the complete solution field and the QoI. While HROMs offer significant computational speedup, enabling efficient simulations, they do exhibit some errors, particularly for testing trajectories. However, their value lies in applications where the detection of the Coanda effect holds paramount importance, even if the selected bifurcation branch is incorrect. Alternatively, for more precise results, HROMs with lower speedups can be employed. △ Less

Submitted 2 May, 2024; v1 submitted 11 July, 2023; originally announced July 2023.

Reduced order models for the buckling of hyperelastic beams

Authors: Federico Pichi, Gianluigi Rozza

Abstract: In this paper, we discuss reduced order modelling approaches to bifurcating systems arising from continuum mechanics benchmarks. The investigation of the beam's deflection is a relevant topic of investigation with fundamental implications on their design for structural analysis and health. When the beams are exposed to external forces, their equilibrium state can undergo to a sudden variation. Thi… ▽ More In this paper, we discuss reduced order modelling approaches to bifurcating systems arising from continuum mechanics benchmarks. The investigation of the beam's deflection is a relevant topic of investigation with fundamental implications on their design for structural analysis and health. When the beams are exposed to external forces, their equilibrium state can undergo to a sudden variation. This happens when a compression, acting along the axial boundaries, exceeds a certain critical value. Linear elasticity models are not complex enough to capture the so-called beam's buckling, and nonlinear constitutive relations, as the hyperelastic laws, are required to investigate this behavior, whose mathematical counterpart is represented by bifurcating phenomena. The numerical analysis of the bifurcating modes and the post-buckling behavior, is usually unaffordable by means of standard high-fidelity techniques such (as the Finite Element method) and the efficiency of Reduced Order Models (ROMs), e.g.\ based on Proper Orthogonal Decomposition (POD), are necessary to obtain consistent speed-up in the reconstruction of the bifurcation diagram. The aim of this work is to provide insights regarding the application of POD-based ROMs for buckling phenomena occurring for 2-D and 3-D beams governed by different constitutive relations. The benchmarks will involve multi-parametric settings with geometrically parametrized domains, where the buckling's location depends on the material and geometrical properties induced by the parameter. Finally, we exploit the acquired notions from these toy problems, to simulate a real case scenario coming from the Norwegian petroleum industry. △ Less

Submitted 31 May, 2023; originally announced May 2023.

Generative Adversarial Reduced Order Modelling

Authors: Dario Coscia, Nicola Demo, Gianluigi Rozza

Abstract: In this work, we present GAROM, a new approach for reduced order modelling (ROM) based on generative adversarial networks (GANs). GANs have the potential to learn data distribution and generate more realistic data. While widely applied in many areas of deep learning, little research is done on their application for ROM, i.e. approximating a high-fidelity model with a simpler one. In this work, we… ▽ More In this work, we present GAROM, a new approach for reduced order modelling (ROM) based on generative adversarial networks (GANs). GANs have the potential to learn data distribution and generate more realistic data. While widely applied in many areas of deep learning, little research is done on their application for ROM, i.e. approximating a high-fidelity model with a simpler one. In this work, we combine the GAN and ROM framework, by introducing a data-driven generative adversarial model able to learn solutions to parametric differential equations. The latter is achieved by modelling the discriminator network as an autoencoder, extracting relevant features of the input, and applying a conditioning mechanism to the generator and discriminator networks specifying the differential equation parameters. We show how to apply our methodology for inference, provide experimental evidence of the model generalisation, and perform a convergence study of the method. △ Less

Submitted 25 May, 2023; originally announced May 2023.

A reduced-order model for segregated fluid-structure interaction solvers based on an ALE approach

Authors: Valentin Nkana Ngan, Giovanni Stabile, Andrea Mola, Gianluigi Rozza

Abstract: This article presents a Galerkin projection model order reduction approach for fluid-structure interaction (FSI) problems in the Finite Volume context. The reduced-order model (ROM) is based on proper orthogonal decomposition (POD), where a reduced basis is formed using energy-dominant POD modes. The reduced basis also consists of characteristics of the POD time modes derived from the POD time mod… ▽ More This article presents a Galerkin projection model order reduction approach for fluid-structure interaction (FSI) problems in the Finite Volume context. The reduced-order model (ROM) is based on proper orthogonal decomposition (POD), where a reduced basis is formed using energy-dominant POD modes. The reduced basis also consists of characteristics of the POD time modes derived from the POD time modes coefficients. In addition, the solution state vector comprises the mesh deformation, considering the structural motion in FSI. The results are obtained by applying the proposed method to time-dependent problems governed by the 2D incompressible Navier-Stokes equations. The main objective of this work is to introduce a hybrid technique mixing up the classical Galerkin-projection approach with a data-driven method to obtain a versatile and accurate algorithm for resolving FSI problems with moving meshes. The effectiveness of this approach is demonstrated in the case study of vortex-induced vibrations (VIV) of a cylinder at Reynolds number Re = 200. The results show the stability and accuracy of the proposed method with respect to the high-dimensional model by capturing transient flow fields and, more importantly, the forces acting on the moving objects. △ Less

Submitted 22 May, 2023; originally announced May 2023.

Filter stabilization for the mildly compressible Euler equations with application to atmosphere dynamics simulations

Authors: Nicola Clinco, Michele Girfoglio, Annalisa Quaini, Gianluigi Rozza

Abstract: We present a filter stabilization technique for the mildly compressible Euler equations that relies on a linear or nonlinear indicator function to identify the regions of the domain where artificial viscosity is needed and determine its amount. For the realization of this technique, we adopt a three step algorithm called Evolve-Filter-Relax (EFR), which at every time step evolves the solution (i.e… ▽ More We present a filter stabilization technique for the mildly compressible Euler equations that relies on a linear or nonlinear indicator function to identify the regions of the domain where artificial viscosity is needed and determine its amount. For the realization of this technique, we adopt a three step algorithm called Evolve-Filter-Relax (EFR), which at every time step evolves the solution (i.e., solves the Euler equations on a coarse mesh), then filters the computed solution, and finally performs a relaxation step to combine the filtered and non-filtered solutions. We show that the EFR algorithm is equivalent to an eddy-viscosity model in Large Eddy Simulation. Three indicator functions are considered: a constant function (leading to a linear filter), a function proportional to the norm of the velocity gradient (recovering a Smagorinsky-like model), and a function based on approximate deconvolution operators. Through well-known benchmarks for atmospheric flow, we show that the deconvolution-based filter yields stable solutions that are much less dissipative than the linear filter and the Samgorinsky-like model and we highlight the efficiency of the EFR algorithm. △ Less

Submitted 22 May, 2023; originally announced May 2023.

A Shape Optimization Pipeline for Marine Propellers by means of Reduced Order Modeling Techniques

Authors: Anna Ivagnes, Nicola Demo, Gianluigi Rozza

Abstract: In this paper, we propose a shape optimization pipeline for propeller blades, applied to naval applications. The geometrical features of a blade are exploited to parametrize it, allowing to obtain deformed blades by perturbating their parameters. The optimization is performed using a genetic algorithm that exploits the computational speed-up of reduced order models to maximize the efficiency of a… ▽ More In this paper, we propose a shape optimization pipeline for propeller blades, applied to naval applications. The geometrical features of a blade are exploited to parametrize it, allowing to obtain deformed blades by perturbating their parameters. The optimization is performed using a genetic algorithm that exploits the computational speed-up of reduced order models to maximize the efficiency of a given propeller. A standard offline-online procedure is exploited to construct the reduced-order model. In an expensive offline phase, the full order model, which reproduces an open water test, is set up in the open-source software OpenFOAM and the same full order setting is used to run the CFD simulations for all the deformed propellers. The collected high-fidelity snapshots and the deformed parameters are used in the online stage to build the non-intrusive reduced-order model. This paper provides a proof of concept of the pipeline proposed, where the optimized propeller improves the efficiency of the original propeller. △ Less

Submitted 15 January, 2024; v1 submitted 12 May, 2023; originally announced May 2023.

Journal ref: International Journal for Numerical Methods in Engineering (2024)

A physics-based reduced order model for urban air pollution prediction

Authors: Moaad Khamlich, Giovanni Stabile, Gianluigi Rozza, László Környei, Zoltán Horváth

Abstract: This article presents an innovative approach for developing an efficient reduced-order model to study the dispersion of urban air pollutants. The need for real-time air quality monitoring has become increasingly important, given the rise in pollutant emissions due to urbanization and its adverse effects on human health. The proposed methodology involves solving the linear advection-diffusion probl… ▽ More This article presents an innovative approach for developing an efficient reduced-order model to study the dispersion of urban air pollutants. The need for real-time air quality monitoring has become increasingly important, given the rise in pollutant emissions due to urbanization and its adverse effects on human health. The proposed methodology involves solving the linear advection-diffusion problem, where the solution of the Reynolds-averaged Navier-Stokes equations gives the convective field. At the same time, the source term consists of an empirical time series. However, the computational requirements of this approach, including microscale spatial resolution, repeated evaluation, and low time scale, necessitate the use of high-performance computing facilities, which can be a bottleneck for real-time monitoring. To address this challenge, a problem-specific methodology was developed that leverages a data-driven approach based on Proper Orthogonal Decomposition with regression (POD-R) coupled with Galerkin projection (POD-G) endorsed with the discrete empirical interpolation method (DEIM). The proposed method employs a feedforward neural network to non-intrusively retrieve the reduced-order convective operator required for online evaluation. The numerical framework was validated on synthetic emissions and real wind measurements. The results demonstrate that the proposed approach significantly reduces the computational burden of the traditional approach and is suitable for real-time air quality monitoring. Overall, the study advances the field of reduced order modeling and highlights the potential of data-driven approaches in environmental modeling and large-scale simulations. △ Less

Submitted 26 May, 2023; v1 submitted 8 May, 2023; originally announced May 2023.

Applicable Methodologies for the Mass Transfer Phenomenon in Tumble Dryers: A Review

Authors: Sajad Salavatidezfouli, Arash Hajisharifi, Michele Girfoglio, Giovanni Stabile, Gianluigi Rozza

Abstract: Tumble dryers offer a fast and convenient way of drying textiles independent of weather conditions and therefore are frequently used in ordinary households. However, artificial drying of textiles consumes considerable amounts of energy, approximately 8.2 percent of the residential electricity consumption is for drying of textiles in northern European countries (Cranston et al., 2019). Several auth… ▽ More Tumble dryers offer a fast and convenient way of drying textiles independent of weather conditions and therefore are frequently used in ordinary households. However, artificial drying of textiles consumes considerable amounts of energy, approximately 8.2 percent of the residential electricity consumption is for drying of textiles in northern European countries (Cranston et al., 2019). Several authors have investigated the aspects of the clothes drying cycle with experimental and numerical methods to understand and improve the process. The first turning point study on understanding the physics of evaporation for tumble dryers was presented by Lambert et al. (1991) in the early 90s. With the aid of Chilton_Colburn analogy, they introduced the concept of area-mass transfer coefficient to address evaporation rate. Afterwards, several experimental or numerical studies were published based on this concept, and furthermore, the model was then developed into 0-dimensional (Deans, 2001) and 1-dimensional (Wei et al., 2017) to gain more accuracy. The evaporation rate is considered to be the main system parameter for dryers with which other performance parameters including drying time, effectiveness, moisture content and efficiency can be estimated. More recent literature focused on utilizing dimensional analysis or image processing techniques to correlate drying indices with system parameters. However, the validity of these regressed models is machine-specific, and hence, cannot be generalized yet. All the previous models for estimating the evaporation rate in tumble dryers are discussed. The review of the related literature showed that all of the previous models for the prediction of the evaporation rate in the clothes dryers have some limitations in terms of accuracy and applicability. △ Less

Submitted 7 April, 2023; originally announced April 2023.

Comments: 25 Pages

Model Order Reduction for Deforming Domain Problems in a Time-Continuous Space-Time Setting

Authors: Fabian Key, Max von Danwitz, Francesco Ballarin, Gianluigi Rozza

Abstract: In the context of simulation-based methods, multiple challenges arise, two of which are considered in this work. As a first challenge, problems including time-dependent phenomena with complex domain deformations, potentially even with changes in the domain topology, need to be tackled appropriately. The second challenge arises when computational resources and the time for evaluating the model beco… ▽ More In the context of simulation-based methods, multiple challenges arise, two of which are considered in this work. As a first challenge, problems including time-dependent phenomena with complex domain deformations, potentially even with changes in the domain topology, need to be tackled appropriately. The second challenge arises when computational resources and the time for evaluating the model become critical in so-called many query scenarios for parametric problems. For example, these problems occur in optimization, uncertainty quantification (UQ), or automatic control and using highly resolved full-order models (FOMs) may become impractical. To address both types of complexity, we present a novel projection-based model order reduction (MOR) approach for deforming domain problems that takes advantage of the time-continuous space-time formulation. We apply it to two examples that are relevant for engineering or biomedical applications and conduct an error and performance analysis. In both cases, we are able to drastically reduce the computational expense for a model evaluation and, at the same time, to maintain an adequate accuracy level. All in all, this work indicates the effectiveness of the presented MOR approach for deforming domain problems taking advantage of a time-continuous space-time setting. △ Less

Submitted 29 March, 2023; originally announced March 2023.

Weighted reduced order methods for uncertainty quantification in computational fluid dynamics

Authors: Julien Genovese, Francesco Ballarin, Gianluigi Rozza, Claudio Canuto

Abstract: In this manuscript we propose and analyze weighted reduced order methods for stochastic Stokes and Navier-Stokes problems depending on random input data (such as forcing terms, physical or geometrical coefficients, boundary conditions). We will compare weighted methods such as weighted greedy and weighted POD with non-weighted ones in case of stochastic parameters. In addition we will analyze diff… ▽ More In this manuscript we propose and analyze weighted reduced order methods for stochastic Stokes and Navier-Stokes problems depending on random input data (such as forcing terms, physical or geometrical coefficients, boundary conditions). We will compare weighted methods such as weighted greedy and weighted POD with non-weighted ones in case of stochastic parameters. In addition we will analyze different sampling and weighting choices to overcome the curse of dimensionality with high dimensional parameter spaces. △ Less

Submitted 25 March, 2023; originally announced March 2023.

GEA: a new finite volume-based open source code for the numerical simulation of atmospheric and ocean flows

Authors: Michele Girfoglio, Annalisa Quaini, Gianluigi Rozza

Abstract: We introduce GEA (Geophysical and Environmental Applications), a new open-source atmosphere and ocean modeling framework within the finite volume C++ library OpenFOAM. Here, we present the development of a non-hydrostatic atmospheric model consisting of a pressure-based solver for the Euler equations written in conservative form using density, momentum, and total energy as variables. We validate t… ▽ More We introduce GEA (Geophysical and Environmental Applications), a new open-source atmosphere and ocean modeling framework within the finite volume C++ library OpenFOAM. Here, we present the development of a non-hydrostatic atmospheric model consisting of a pressure-based solver for the Euler equations written in conservative form using density, momentum, and total energy as variables. We validate the solver for two idealized test cases involving buoyancy driven flows: smooth and non-smooth rising thermal bubble. Through qualitative and quantitative comparisons against numerical data available in the literature, we show that our approach is accurate. △ Less

Submitted 18 March, 2023; originally announced March 2023.

Comments: 8 pages, 2 figures, 2 tables. arXiv admin note: substantial text overlap with arXiv:2302.04836

A DeepONet multi-fidelity approach for residual learning in reduced order modeling

Authors: Nicola Demo, Marco Tezzele, Gianluigi Rozza

Abstract: In the present work, we introduce a novel approach to enhance the precision of reduced order models by exploiting a multi-fidelity perspective and DeepONets. Reduced models provide a real-time numerical approximation by simplifying the original model. The error introduced by the such operation is usually neglected and sacrificed in order to reach a fast computation. We propose to couple the model… ▽ More In the present work, we introduce a novel approach to enhance the precision of reduced order models by exploiting a multi-fidelity perspective and DeepONets. Reduced models provide a real-time numerical approximation by simplifying the original model. The error introduced by the such operation is usually neglected and sacrificed in order to reach a fast computation. We propose to couple the model reduction to a machine learning residual learning, such that the above-mentioned error can be learned by a neural network and inferred for new predictions. We emphasize that the framework maximizes the exploitation of high-fidelity information, using it for building the reduced order model and for learning the residual. In this work, we explore the integration of proper orthogonal decomposition (POD), and gappy POD for sensors data, with the recent DeepONet architecture. Numerical investigations for a parametric benchmark function and a nonlinear parametric Navier-Stokes problem are presented. △ Less

Submitted 17 November, 2023; v1 submitted 24 February, 2023; originally announced February 2023.

A Non-Intrusive Data-Driven Reduced Order Model for Parametrized CFD-DEM Numerical Simulations

Authors: Arash Hajisharifi, Francesco Romano`, Michele Girfoglio, Andrea Beccari, Domenico Bonanni, Gianluigi Rozza

Abstract: The investigation of fluid-solid systems is very important in a lot of industrial processes. From a computational point of view, the simulation of such systems is very expensive, especially when a huge number of parametric configurations needs to be studied. In this context, we develop a non-intrusive data-driven reduced order model (ROM) built using the proper orthogonal decomposition with interp… ▽ More The investigation of fluid-solid systems is very important in a lot of industrial processes. From a computational point of view, the simulation of such systems is very expensive, especially when a huge number of parametric configurations needs to be studied. In this context, we develop a non-intrusive data-driven reduced order model (ROM) built using the proper orthogonal decomposition with interpolation (PODI) method for Computational Fluid Dynamics (CFD) -- Discrete Element Method (DEM) simulations. The main novelties of the proposed approach rely in (i) the combination of ROM and FV methods, (ii) a numerical sensitivity analysis of the ROM accuracy with respect to the number of POD modes and to the cardinality of the training set and (iii) a parametric study with respect to the Stokes number. We test our ROM on the fluidized bed benchmark problem. The accuracy of the ROM is assessed against results obtained with the FOM both for Eulerian (the fluid volume fraction) and Lagrangian (position and velocity of the particles) quantities. We also discuss the efficiency of our ROM approach. △ Less

Submitted 24 February, 2023; originally announced February 2023.

Journal ref: July 2023-Journal of Computational Physics

A Graph-based Framework for Complex System Simulating and Diagnosis with Automatic Reconfiguration

Authors: Martina Teruzzi, Nicola Demo, Gianluigi Rozza

Abstract: Fault detection has a long tradition: the necessity to provide the most accurate diagnosis possible for a process plant criticality is somehow intrinsic in its functioning. Continuous monitoring is a possible way for early detection. However, it is somehow fundamental to be able to actually simulate failures. Reproducing the issues remotely allows to quantify in advance their consequences, causing… ▽ More Fault detection has a long tradition: the necessity to provide the most accurate diagnosis possible for a process plant criticality is somehow intrinsic in its functioning. Continuous monitoring is a possible way for early detection. However, it is somehow fundamental to be able to actually simulate failures. Reproducing the issues remotely allows to quantify in advance their consequences, causing literally no real damage. Within this context, signed directed graphs have played an essential role within the years, managing to model with a relatively simple theory diverse elements of an industrial network, as well as the logic relations between them.\\ In this work we present a quantitative approach, employing directed graphs to the simulation and automatic reconfiguration of a fault in a network. To model the typical operation of industrial plants, we propose several additions with respect to the standard graphs: 1. a quantitative measure to control the overall residual capacity, 2. nodes of different categories - and then different behaviors - and 3. a fault propagation procedure based on the predecessors and the redundancy of the system. The obtained graph is able to mimic the behaviour of the real target plant when one or more faults occur. Additionally, we also implement a generative approach capable to activate a particular category of nodes in order to contain the issue propagation, equipping the network with the capability of reconfigure itself and resulting then in a mathematical tool useful not only for simulating and monitoring, but also to design and optimize complex plants. The final asset of the system is provided in output with its complete diagnostics, and a detailed description of the steps that have been carried out to obtain the final realization. △ Less

Submitted 10 February, 2023; originally announced February 2023.

Journal ref: Mathematics in Engineering (2024)

Validation of an OpenFOAM-based solver for the Euler equations with benchmarks for mesoscale atmospheric modeling

Authors: Michele Girfoglio, Annalisa Quaini, Gianluigi Rozza

Abstract: Within OpenFOAM, we develop a pressure-based solver for the Euler equations written in conservative form using density, momentum, and total energy as variables. Under simplifying assumptions, these equations are used to describe non-hydrostatic atmospheric flow. For the stabilization of the Euler equations and to capture sub-grid processes, we consider two Large Eddy Simulation models: the classic… ▽ More Within OpenFOAM, we develop a pressure-based solver for the Euler equations written in conservative form using density, momentum, and total energy as variables. Under simplifying assumptions, these equations are used to describe non-hydrostatic atmospheric flow. For the stabilization of the Euler equations and to capture sub-grid processes, we consider two Large Eddy Simulation models: the classical Smagorinsky model and the one equation eddy-viscosity model. To achieve high computational efficiency, our solver uses a splitting scheme that decouples the computation of each variable. The numerical results obtained with our solver are validated against numerical data available in the literature for two classical benchmarks: the rising thermal bubble and the density current. Through qualitative and quantitative comparisons, we show that our approach is accurate. This work is meant to lay the foundation for a new open source package specifically created for quick assessment of new computational approaches for the simulation of atmospheric flows at the mesoscale level △ Less

Submitted 9 February, 2023; originally announced February 2023.

Reduced Basis, Embedded Methods and Parametrized Levelset Geometry

Authors: Efthymios N. Karatzas, Giovanni Stabile, Francesco Ballarin, Gianluigi Rozza

Abstract: In this chapter we examine reduced order techniques for geometrical parametrized heat exchange systems, Poisson, and flows based on Stokes, steady and unsteady incompressible Navier-Stokes and Cahn-Hilliard problems. The full order finite element methods, employed in an embedded and/or immersed geometry framework, are the Shifted Boundary (SBM) and the Cut elements (CutFEM) methodologies, with app… ▽ More In this chapter we examine reduced order techniques for geometrical parametrized heat exchange systems, Poisson, and flows based on Stokes, steady and unsteady incompressible Navier-Stokes and Cahn-Hilliard problems. The full order finite element methods, employed in an embedded and/or immersed geometry framework, are the Shifted Boundary (SBM) and the Cut elements (CutFEM) methodologies, with applications mainly focused in fluids. We start by introducing the Nitsche's method, for both SBM/CutFEM and parametrized physical problems as well as the high fidelity approximation. We continue with the full order parameterized Nitsche shifted boundary variational weak formulation, and the reduced order modeling ideas based on a Proper Orthogonal Decomposition Galerkin method and geometrical parametrization, quoting the main differences and advantages with respect to a reference domain approach used for classical finite element methods, while stability issues may overcome employing supremizer enrichment methodologies. Numerical experiments verify the efficiency of the introduced ``hello world'' problems considering reduced order results in several cases for one, two, three and four dimensional geometrical kind of parametrization. We investigate execution times, and we illustrate transport methods and improvements. A list of important references related to unfitted methods and reduced order modeling are [11, 8, 9, 10, 7, 6, 12]. △ Less

Submitted 29 January, 2023; originally announced January 2023.

MSC Class: 78M34; 97N40; 35Q35

A two stages Deep Learning Architecture for Model Reduction of Parametric Time-Dependent Problems

Authors: Isabella Carla Gonnella, Martin W. Hess, Giovanni Stabile, Gianluigi Rozza

Abstract: 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 l… ▽ More 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. △ Less

Submitted 25 January, 2023; v1 submitted 24 January, 2023; originally announced January 2023.

Stabilized Weighted Reduced Order Methods for Parametrized Advection-Dominated Optimal Control Problems governed by Partial Differential Equations with Random Inputs

Authors: Fabio Zoccolan, Maria Strazzullo, Gianluigi Rozza

Abstract: In this work, we analyze Parametrized Advection-Dominated distributed Optimal Control Problems with random inputs in a Reduced Order Model (ROM) context. All the simulations are initially based on a finite element method (FEM) discretization; moreover, a space-time approach is considered when dealing with unsteady cases. To overcome numerical instabilities that can occur in the optimality system f… ▽ More In this work, we analyze Parametrized Advection-Dominated distributed Optimal Control Problems with random inputs in a Reduced Order Model (ROM) context. All the simulations are initially based on a finite element method (FEM) discretization; moreover, a space-time approach is considered when dealing with unsteady cases. To overcome numerical instabilities that can occur in the optimality system for high values of the Péclet number, we consider a Streamline Upwind Petrov-Galerkin technique applied in an optimize-then-discretize approach. We combine this method with the ROM framework in order to consider two possibilities of stabilization: Offline-Only stabilization and Offline-Online stabilization. Moreover, we consider random parameters and we use a weighted Proper Orthogonal Decomposition algorithm in a partitioned approach to deal with the issue of uncertainty quantification. Several quadrature techniques are used to derive weighted ROMs: tensor rules, isotropic sparse grids, Monte-Carlo and quasi Monte-Carlo methods. We compare all the approaches analyzing relative errors between the FEM and ROM solutions and the computational efficiency based on the speedup-index. △ Less

Submitted 5 January, 2023; originally announced January 2023.

Comments: 27 pages, 42 figures, 4 tables

MSC Class: 49J20; 49M41; 65M60; 60H25; 60H35

A Streamline upwind Petrov-Galerkin Reduced Order Method for Advection-Dominated Partial Differential Equations under Optimal Control

Authors: Fabio Zoccolan, Maria Strazzullo, Gianluigi Rozza

Abstract: In this paper we will consider distributed Linear-Quadratic Optimal Control Problems dealing with Advection-Diffusion PDEs for high values of the Péclet number. In this situation, computational instabilities occur, both for steady and unsteady cases. A Streamline Upwind Petrov-Galerkin technique is used in the optimality system to overcome these unpleasant effects. We will apply a finite element m… ▽ More In this paper we will consider distributed Linear-Quadratic Optimal Control Problems dealing with Advection-Diffusion PDEs for high values of the Péclet number. In this situation, computational instabilities occur, both for steady and unsteady cases. A Streamline Upwind Petrov-Galerkin technique is used in the optimality system to overcome these unpleasant effects. We will apply a finite element method discretization in an optimize-then-discretize approach. Concerning the parabolic case, a stabilized space-time framework will be considered and stabilization will also occur in both bilinear forms involving time derivatives. Then we will build Reduced Order Models on this discretization procedure and two possible settings can be analyzed: whether or not stabilization is needed in the online phase, too. In order to build the reduced bases for state, control, and adjoint variables we will consider a Proper Orthogonal Decomposition algorithm in a partitioned approach. It is the first time that Reduced Order Models are applied to stabilized parabolic problems in this setting. The discussion is supported by computational experiments, where relative errors between the FEM and ROM solutions are studied together with the respective computational times. △ Less

Submitted 1 April, 2024; v1 submitted 5 January, 2023; originally announced January 2023.

Comments: 27 pages, 36 figures, 4 tables

MSC Class: 49J20; 49M41; 65M60

Assessment of URANS and LES Methods in Predicting Wake Shed Behind a Vertical Axis Wind Turbine

Authors: Armin Sheidani, Sajad Salavatidezfouli, Giovanni Stabile, Gianluigi Rozza

Abstract: In order to shed light on the Vertical-Axis Wind Turbines (VAWT) wake characteristics, in this paper we present high-fidelity CFD simulations of the flow around an exemplary H-shaped VAWT turbine, and we propose to apply Proper Orthogonal Decomposition (POD) to the computed flow field in the near wake of the rotor. The turbine under consideration was widely studied in previous experimental and com… ▽ More In order to shed light on the Vertical-Axis Wind Turbines (VAWT) wake characteristics, in this paper we present high-fidelity CFD simulations of the flow around an exemplary H-shaped VAWT turbine, and we propose to apply Proper Orthogonal Decomposition (POD) to the computed flow field in the near wake of the rotor. The turbine under consideration was widely studied in previous experimental and computational investigations. In the first part of the study, multiple Reynolds-Averaged Navier-Stokes (RANS) simulations were performed at the Tip Speed Ratio (TSR) of peak power coefficient, to select the most accurate turbulence model with respect to available data. In the following step, further RANS numerical simulations were performed at different TSRs to compare the power coefficient against experimental data. Then, Large Eddy Simulation (LES) was applied for multiple TSR conditions. The spatial and temporal POD modes along with modal energy for the RANS and LES results were extracted, and the performance of the turbulence models was assessed. Also, an interpretation of the POD modes with respect to the flow structures was given to highlight the most significant time and length scales of the predictions considering the different dynamical levels of approximations of the computational models. △ Less

Submitted 26 December, 2022; originally announced December 2022.

Comments: 15 Pages