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 reduced-order approaches are not taken into account. In particular, in this work we focus on a Proper Orthogonal Decomposition (POD)-based formulation and our goal is to build a closure or correction model, aimed to re-introduce the contribution of the discarded modes.

In the present work, we introduce a data-driven approach to enhance the accuracy of non-intrusive Reduced Order Models (ROMs). In particular, we focus on ROMs built using Proper Orthogonal Decomposition (POD) in an under-resolved and marginally-resolved regime, i.e. when the number of modes employed is not enough to capture the system dynamics. We propose a method to re-introduce the contribution of neglected modes through a quadratic correction term, given by the action of a quadratic operator on the POD coefficients.

The widespread adoption of embedded vision systems in industrial applications has highlighted the limitations of deep learning models, which are characterized by a high number of parameters. This is representing a significant concern within the scientific community due to the increased computational resources and memory required for training and inference of these models. Addressing this, we propose a flexible and effective methodology for neural network compression that integrates a pluggable dimensionality reduction layer with a Knowledge Distillation (KD) approach.

The conservation of cultural heritage increasingly relies on integrating technological innovation with domain expertise to ensure effective monitoring and predictive maintenance. This paper presents a novel framework to support the preservation of cultural assets, combining Internet of Things (IoT) and Artificial Intelligence (AI) technologies, enhanced with the physical knowledge of phenomena. The framework is structured into four functional layers that permit the analysis of 3D models of cultural assets and elaborate simulations based on the knowledge acquired from data and physics.

Reynolds-Averaged Navier-Stokes (RANS) models are widely used for turbulent flow simulations due to their computational efficiency, but their accuracy strongly depends on the selected turbulence closure and may vary across the flow domain. Space-dependent model aggregation has been shown to improve RANS predictions by combining multiple turbulence models, although at the cost of repeated high-fidelity simulations.

Numerical simulations of turbulent flows are well known to pose extreme computational challenges because of the huge number of dynamical degrees of freedom required to correctly describe the complex multiscale 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.

We present a novel approach to define the filter and relax steps in the evolve-filter-relax (EFR) framework for simulating turbulent flows. The EFR main advantages are its ease of implementation and computational efficiency. However, as it only contains two parameters (one for the filter step and one for the relax step) its flexibility is rather limited. In this work, we propose a data-driven approach in which the optimal filter is found based on DNS data in the frequency domain.

We propose a reinforcement learning (RL) framework for the dynamic selection of the filter parameter in Evolve-Filter (EF) regularization strategies for incompressible turbulent flows. Instead of prescribing the filter radius heuristically, the RL agent learns to adaptively control the filtering intensity in time, balancing numerical stability and physical accuracy. The methodology is assessed on two benchmark problems with fundamentally different dynamics: flow past a cylinder and decaying homogeneous turbulence.

In recent years, large-scale numerical simulations played an essential role in estimating the effects of explosion events in urban environments, for the purpose of ensuring the security and safety of cities. Such simulations are computationally expensive and, often, the time taken for one single computation is large and does not permit parametric studies. The aim of this work is therefore to facilitate real-time and multi-query calculations by employing a non-intrusive Reduced Order Method (ROM). 

Spatial filters have played a central role in large eddy simulation and, more recently, in reduced order model (ROM) stabilization for convection-dominated flows. Nevertheless, important open questions remain: in under-resolved regimes, which filter is most suitable for a given stabilization or closure model? Moreover, once a filter is selected, how should its parameters, such as the filter radius, be determined? Addressing these questions is essential for the reliable design and performance of filter-based stabilization strategies.

Graph Neural Networks (GNNs) are emerging as powerful tools for nonlinear Model Order Reduction (MOR) of time-dependent parameterized Partial Differential Equations (PDEs). However, existing methodologies struggle to combine geometric inductive biases with interpretable latent behavior, overlooking dynamics-driven features or disregarding spatial information.

We propose an efficient hyper-reduced order model (HROM) designed for segregated finite-volume solvers in geometrically parametrized problems. The method follows a discretize-then-project strategy: the full-order operators are first assembled using finite volume or finite element discretizations and then projected onto low-dimensional spaces using a small set of spatial sampling points, selected through hyper-reduction techniques such as DEIM. This approach removes the dependence of the online computational cost on the full mesh size.

Recently, general fractional calculus was introduced by Kochubei (2011) and Luchko (2021) as a further generalisation of fractional calculus, where the derivative and integral operator admits arbitrary kernel. Such a formalism will have many applications in physics and engineering, since the kernel is no longer restricted. We first extend the work of Al-Refai and Luchko (2023) on finite interval to arbitrary orders. Followed by, developing an efficient Petrov-Galerkin scheme by introducing Jacobi convolution series as basis functions.

The two-layer quasi-geostrophic equations (2QGE) serve as a simplified model for simulating wind-driven, stratified ocean flows. However, their numerical simulation remains computationally expensive due to the need for high-resolution meshes to capture a wide range of turbulent scales. This becomes especially problematic when several simulations need to be run because of, e.g., uncertainty in the parameter settings.

Real-world problems encountered in Computational Fluid Dynamics (CFD) are often governed by complex systems of parametrized partial differential equations. The resolution of such problems requires the employment of advanced numerical tools for simulation purposes. Classic numerical simulations, which aim to accurately replicate experimental data, may require thousands or even millions of degrees of freedom, resulting in time and memory-intensive processes.

This article presents a Galerkin projection-based reduced-order modeling (ROM) approach for segregated fluid–structure interaction (FSI) problems, formulated within an Arbitrary Lagrangian–Eulerian (ALE) framework at low Reynolds numbers using the Finite Volume Method (FVM). The ROM is constructed using Proper Orthogonal Decomposition (POD) and incorporates a data-driven technique that combines classical Galerkin projection with radial basis function (RBF) networks.

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.

Cardiovascular diseases are a leading cause of death in the world, driving the development of patient-specific and benchmark models for blood flow analysis. This chapter provides a theoretical overview of the main categories of Reduced Order Models (ROMs), focusing on both projection-based and data-driven approaches within a classical setup. We then present a hybrid ROM tailored for simulating blood flow in a patient-specific aortic geometry.

The current work aims to incorporate physics-based loss in Physics Informed Neural Network (PINN) directly using the numerical residual obtained from the governing equation in any dicretized forward solver. PINN's major difficulties in coupling with external forward solvers arise from the inability to access the discretized form (Finite difference, finite volume, finite element, etc.) of the governing equation directly through the network and to include them in its computational graph.

This paper presents a projection-based reduced order modelling (ROM) framework for unsteady parametrized optimal control problems (OCP s) arising from cardiovascular (CV) applications. In real-life scenarios, accurately defining outflow boundary conditions in patient-specific models poses significant challenges due to complex vascular morphologies, physiological conditions, and high computational demands. These challenges make it difficult to compute realistic and reliable CV hemodynamics by incorporating clinical data such as 4D magnetic resonance imaging.

The flow behavior in the continuous casting tundish plays a critical role in steel quality and is typically characterized via residence time distribution (RTD) curves. This study investigates the fluid flow behaviour in a single-strand tundish using numerical and experimental approaches. Full-order model (FOM) steady-state simulations were conducted under both isothermal and non-isothermal conditions to assess the influence of thermal buoyancy on the flow characteristics.

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.

We present a novel reduced-order Model (ROM) that leverages optimal transport (OT) theory and displacement interpolation to enhance the representation of nonlinear dynamics in complex systems. While traditional ROM techniques face challenges in this scenario, especially when data (i.e., observational snapshots) is limited, our method addresses these issues by introducing a data augmentation strategy based on OT principles. The proposed framework generates interpolated solutions tracing geodesic paths in the space of probability distributions, enriching the training dataset for the ROM.

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.

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.

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.

The two-layer quasi-geostrophic equations (2QGE) is a simplified model that describes the dynamics of a stratified, wind-driven ocean in terms of potential vorticity and stream function. Its numerical simulation is plagued by a high computational cost due to the size of the typical computational domain and the need for high resolution to capture the full spectrum of turbulent scales. In this paper, we present a data-driven reduced order model (ROM) for the 2QGE that drastically reduces the computational time to predict ocean dynamics, especially when there are variable physical parameters.

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.

Autoregressive and recurrent networks have achieved remarkable progress across various fields, from weather forecasting to molecular generation and Large Language Models. Despite their strong predictive capabilities, these models lack a rigorous framework for addressing uncertainty, which is key in scientific applications such as PDE solving, molecular generation and Machine Learning Force Fields. To address this shortcoming we present BARNN: a variational Bayesian Autoregressive and Recurrent Neural Network.

Traditional linear approximation methods, such as proper orthogonal decomposition and the reduced basis method, are ineffective for transport-dominated problems due to the slow decay of the Kolmogorov  -width. This results in reduced-order models that are both inefficient and inaccurate.

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.

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.

ATHENA is an open source Python package for reduction in parameter space. It implements several advanced numerical analysis techniques such as Active Subspaces (AS), Kernel-based Active Subspaces (KAS), and Nonlinear Level-set Learning (NLL) method. It is intended as a tool for regression, sensitivity analysis, and in general to enhance existing numerical simulations’ pipelines tackling the curse of dimensionality.

To shed light on the effect of the icing phenomenon on the vertical-axis wind turbine (VAWT) wake characteristics, we present a high-fidelity computational fluid dynamics simulation of the flow field of H-Darrieus turbine under the icing conditions. To address continuous geometry alteration due to the icing and predefined motion of the VAWT, a pseudo-steady approach proposed by Baizhuma et al. [“Numerical method to predict ice accretion shapes and performance penalties for rotating vertical axis wind turbines under icing conditions,” J. Wind Eng. Ind. Aerodyn.

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 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).

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.

Parameter space reduction has been proved to be a crucial tool to speed-up the execution of many numerical tasks such as optimization, inverse problems, sensitivity analysis, and surrogate models’ design, especially when in presence of high-dimensional parametrized systems. In this work we propose a new method called local active subspaces (LAS), which explores the synergies of active subspaces with supervised clustering techniques in order to carry out a more efficient dimension reduction in the parameter space.

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.

In this manuscript, we combine non-intrusive reduced-order models (ROMs) with space-dependent aggregation techniques to build a mixed-ROM, able to accurately capture the flow dynamics in different physical settings. The flow prediction obtained using the mixed formulation is derived from a convex combination of the predictions of several previously trained reduced-order models (ROMs), with each model assigned a space-dependent weight.

This research study explored 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 was conducted.

Non-affine parametric dependencies, nonlinearities and advection-dominated regimes of the model of interest can result in a slow Kolmogorov n-width decay, which precludes the realization of efficient reduced-order models based on linear subspace approximations. Among the possible solutions, there are purely data-driven methods that leverage autoencoders and their variants to learn a latent representation of the dynamical system, and then evolve it in time with another architecture.

Nowadays, the shipbuilding industry is facing a radical change toward solutions with a smaller environmental impact. This can be achieved with low emissions engines, optimized shape designs with lower wave resistance and noise generation, and by reducing the metal raw materials used during the manufacturing. This work focuses on the last aspect by presenting a complete structural optimization pipeline for modern passenger ship hulls which exploits advanced model order reduction techniques to reduce the dimensionality of both input parameters and outputs of interest.

In this work the development of a machine learning-based Reduced Order Model (ROM) for the investigation of hemodynamics in a patient-specific configuration of Coronary Artery Bypass Graft (CABG) is proposed. The computational domain is referred to left branches of coronary arteries when a stenosis of the Left Main Coronary Artery (LMCA) occurs. The method extracts a reduced basis space from a collection of high-fidelity solutions via a Proper Orthogonal Decomposition (POD) algorithm and employs Artificial Neural Networks (ANNs) for the computation of the modal coefficients.

The aim of this work is to present a model reduction technique in the framework of optimal control problems for partial differential equations. We combine two approaches used for reducing the computational cost of the mathematical numerical models: domain-decomposition (DD) methods and reduced-order modelling (ROM). In particular, we consider an optimisation-based domain-decomposition algorithm for the parameter-dependent stationary incompressible Navier-Stokes equations.

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.

In this work, we propose a model order reduction framework to deal with inverse problems in a non-intrusive setting. Inverse problems, especially in a partial differential equation context, require a huge computational load due to the iterative optimization process. To accelerate such a procedure, we apply a numerical pipeline that involves artificial neural networks to parametrize the boundary conditions of the problem in hand, compress the dimensionality of the (full-order) snapshots, and approximate the parametric solution manifold.

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.

We present a Large Eddy Simulation (LES) approach based on a nonlinear differential low-pass filter for the simulation of two-dimensional barotropic flows with under-refined meshes. For the implementation of such model, we choose a segregated three-step algorithm combined with a computationally efficient Finite Volume method. We assess the performance of our approach with the classical double-gyre wind forcing benchmark.

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.

Convolutional Neural Network (CNN) is one of the most important architectures in deep learning. The fundamental building block of a CNN is a trainable filter, represented as a discrete grid, used to perform convolution on discrete input data. In this work, we propose a continuous version of a trainable convolutional filter able to work also with unstructured data. This new framework allows exploring CNNs beyond discrete domains, enlarging the usage of this important learning technique for many more complex problems.

The focus of this work is on the application of classical Model Order Reduction techniques, such as Active Subspaces and Proper Orthogonal Decomposition, to Deep Neural Networks. We propose a generic methodology to reduce the number of layers in a pre-trained network by combining the aforementioned techniques for dimensionality reduction with input-output mappings, such as Polynomial Chaos Expansion and Feedforward Neural Networks.

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.

This work deals with the investigation of bifurcating fluid phenomena using a reduced order modelling setting aided by artificial neural networks. We discuss the POD-NN approach dealing with non-smooth solutions set of nonlinear parametrized PDEs. Thus, we study the Navier–Stokes equations describing: (i) the Coanda effect in a channel, and (ii) the lid driven triangular cavity flow, in a physical/geometrical multi-parametrized setting, considering the effects of the domain’s configuration on the position of the bifurcation points.

Numerical stabilization is often used to eliminate (alleviate) the spurious oscillations generally produced by full order models (FOMs) in under-resolved or marginally-resolved simulations of convection-dominated flows. In this paper, we investigate the role of numerical stabilization in reduced order models (ROMs) of marginally-resolved convection-dominated flows. Specifically, we investigate the FOM-ROM consistency, i.e., whether the numerical stabilization is beneficial both at the FOM and the ROM level.

In this paper, we present recent efforts to develop reduced order modeling (ROM) capabilities for spectral element methods (SEM). Namely, we detail the implementation of ROM for both continuous Galerkin and discontinuous Galerkin methods in the spectral/hp element library Nektar++. The ROM approaches adopted are intrusive methods based on the proper orthogonal decomposition (POD). They admit an offline-online decomposition, such that fast evaluations for parameter studies and many-queries are possible.

This work investigates the use of sparse polynomial interpolation as a model order reduction method for the parametrized incompressible Navier–Stokes equations. Numerical results are presented underscoring the validity of sparse polynomial approximations and comparing with established reduced basis techniques. Two numerical models serve to assess the accuracy of the reduced order models (ROMs), in particular parametric nonlinearities arising from curved geometries are investigated in detail.

In the study of micro-swimmers, both artificial and biological ones, many-query problems arise naturally. Even with the use of advanced high performance computing (HPC), it is not possible to solve this kind of problems in an acceptable amount of time. Various approximations of the Stokes equation have been considered in the past to ease such computational efforts but they introduce non- negligible errors that can easily make the solution of the problem inaccurate and unreliable.

This work deals with optimal control problems as a strategy to drive bifurcating solution of nonlinear parametrized partial differential equations towards a desired branch. Indeed, for these governing equations, multiple solution configurations can arise from the same parametric instance. We thus aim at describing how optimal control allows to change the solution profile and the stability of state solution branches.

In this work we present an integrated computational pipeline involving several model order reduction techniques for industrial and applied mathematics, as emerging technology for product and/or process design procedures. Its data-driven nature and its modularity allow an easy integration into existing pipelines. We describe a complete optimization framework with automated geometrical parameterization, reduction of the dimension of the parameter space, and non-intrusive model order reduction such as dynamic mode decomposition and proper orthogonal decomposition with interpolation.

We aim to reconstruct the latent space dynamics of high dimensional, quasi-stationary systems using model order reduction via the spectral proper orthogonal decomposition (SPOD).

We consider an optimal flow control problem in a patient-specific coronary artery bypass graft with the aim of matching the blood flow velocity with given measurements as the Reynolds number varies in a physiological range. Blood flow is modelled with the steady incompressible Navier-Stokes equations. The geometry consists in a stenosed left anterior descending artery where a single bypass is performed with the right internal thoracic artery.

We develop a Proper Orthogonal Decomposition (POD)-Galerkin based Reduced Order Model (ROM) for the efficient numerical simulation of the parametric Navier–Stokes equations in the stream function-vorticity formulation. Unlike previous works, we choose different reduced coefficients for the vorticity and stream function fields. In addition, for parametric studies we use a global POD basis space obtained from a database of time dependent full order snapshots related to sample points in the parameter space.

In this work we propose reduced order methods as a reliable strategy to efficiently solve parametrized optimal control problems governed by shallow waters equations in a solution tracking setting. The physical parametrized model we deal with is nonlinear and time dependent: this leads to very time consuming simulations which can be unbearable e.g. in a marine environmental monitoring plan application. Our aim is to show how reduced order modelling could help in studying different configurations and phenomena in a fast way.

This paper focuses on reduced-order models (ROMs) built for the efficient treatment of PDEs having solutions that bifurcate as the values of multiple input parameters change. First, we consider a method called local ROM that uses k-means algorithm to cluster snapshots and construct local POD bases, one for each cluster. We investigate one key ingredient of this approach: the local basis selection criterion.

In this work, we propose to efficiently solve time dependent parametrized optimal control problems governed by parabolic partial differential equations through the certified reduced basis method. In particular, we will exploit an error estimator procedure, based on easy-to-compute quantities which guarantee a rigorous and efficient bound for the error of the involved variables. First of all, we propose the analysis of the problem at hand, proving its well-posedness thanks to Nečas - Babuška theory for distributed and boundary controls in a space-time formulation.

Multi-fidelity models are of great importance due to their capability of fusing information coming from different numerical simulations, surrogates, and sensors. We focus on the approximation of high-dimensional scalar functions with low intrinsic dimensionality. By introducing a low dimensional bias we can fight the curse of dimensionality affecting these quantities of interest, especially for many-query applications. We seek a gradient-based reduction of the parameter space through linear active subspaces or a nonlinear transformation of the input space.

Geometrically parametrized Partial Differential Equations are nowadays widely used in many different fields as, for example, shape optimization processes or patient specific surgery studies. The focus of this work is on some advances for this topic, capable of increasing the accuracy with respect to previous approaches while relying on a high cost-benefit ratio performance.

We present a stabilized POD-Galerkin reduced order method (ROM) for a Leray model. For the implementation of the model, we combine a two-step algorithm called Evolve-Filter (EF) with a computationally efficient finite volume method. In both steps of the EF algorithm, velocity and pressure fields are approximated using different POD basis and coefficients. To achieve pressure stabilization, we consider and compare two strategies: the pressure Poisson equation and the supremizer enrichment of the velocity space.

We propose a Proper Orthogonal Decomposition (POD)-Galerkin based Reduced Order Model (ROM) for an implementation of the Leray model that combines a two-step algorithm called Evolve-Filter (EF) with a computationally efficient finite volume method. The main novelty of the proposed approach relies in applying spatial filtering both for the collection of the snapshots and in the reduced order model, as well as in considering the pressure field at reduced level. In both steps of the EF algorithm, velocity and pressure fields are approximated by using different POD basis and coefficients.

This work recasts time-dependent optimal control problems governed by partial differential equations in a Dynamic Mode Decomposition with control framework. Indeed, since the numerical solution of such problems requires a lot of computational effort, we rely on this specific data-driven technique, using both solution and desired state measurements to extract the underlying system dynamics.

Model reduction, Hydroacoustics, Large eddy simulation, Ffowcs Williams and Hawkings, Dynamic mode decomposition, Proper orthogonal decomposition},

In this work we focus on two different methods to deal with parametrized partial differential equations in an efficient and accurate way. Starting from high fidelity approximations built via the hierarchical model reduction discretization, we consider two approaches, both based on a projection model reduction technique. The two methods differ for the algorithm employed during the construction of the reduced basis. In particular, the former employs the proper orthogonal decomposition, while the latter relies on a greedy algorithm according to the certified reduced basis technique.

A Finite-Volume based POD-Galerkin reduced order modeling strategy for steady-state Reynolds averaged Navier--Stokes (RANS) simulation is extended for low-Prandtl number flow. The reduced order model is based on a full order model for which the effects of buoyancy on the flow and heat transfer are characterized by varying the Richardson number. The Reynolds stresses are computed with a linear eddy viscosity model. A single gradient diffusion hypothesis, together with a local correlation for the evaluation of the turbulent Prandtl number, is used to model the turbulent heat fluxes.

In this work we present an advanced computational pipeline for the approximation and prediction of the lift coefficient of a parametrized airfoil profile. The non-intrusive reduced order method is based on dynamic mode decomposition (DMD) and it is coupled with dynamic active subspaces (DyAS) to enhance the future state prediction of the target function and reduce the parameter space dimensionality.

In this work we present a Reduced Basis VMS-Smagorinsky Boussinesq model, applied to natural convection problems in a variable height cavity, in which the buoyancy forces are involved. We take into account in this problem both physical and geometrical parametrizations, considering the Rayleigh number as a parameter, so as the height of the cavity. We perform an Empirical Interpolation Method to approximate the sub-grid eddy viscosity term that lets us obtain an affine decomposition with respect to the parameters.

PyGeM is an open source Python package which allows to easily parametrize and deform 3D object described by CAD files or 3D meshes. It implements several morphing techniques such as free form deformation, radial basis function interpolation, and inverse distance weighting. Due to its versatility in dealing with different file formats it is particularly suited for researchers and practitioners both in academia and in industry interested in computational engineering simulations and optimization studies.

This work describes the implementation of a data-driven approach for the reduction of the complexity of parametrical partial differential equations (PDEs) employing Proper Orthogonal Decomposition (POD) and Gaussian Process Regression (GPR). This approach is applied initially to a literature case, the simulation of the stokes problems, and in the following to a real-world industrial problem, inside a shape optimization pipeline for a naval engineering problem.

In the present work, we investigate a cut finite element method for the parameterized system of second-order equations stemming from the splitting approach of a fourth order nonlinear geometrical PDE, namely the Cahn-Hilliard system. We manage to tackle the instability issues of such methods whenever strong nonlinearities appear and to utilize their flexibility of the fixed background geometry -- and mesh -- characteristic, through which, one can avoid e.g.

The J-integral is formulated in a direct manner for a gel consisting of a cross-linked polymer network and a mobile solvent. The form of the J-integral is given for a formulation that exploits the Helmholtz energy density of the gel and expressions are provided for it in both the unswollen reference configuration of the polymer network and in the current swollen configuration of the gel when small strains are superimposed on the swollen state.

Understanding the mechanics of brain embryogenesis can provide insights on pathologies related to brain development, such as lissencephaly, a genetic disease which causes a reduction of the number of cerebral sulci. Recent experiments on brain organoids have confirmed that gyrification, i.e. the formation of the folded structures of the brain, is triggered by the inhomogeneous growth of the peripheral region. However, the rheology of these cellular aggregates and the mechanics of lissencephaly are still matter of debate.

A parametric, hybrid reduced order model approach based on the Proper Orthogonal Decomposition with both Galerkin projection and interpolation based on Radial Basis Functions method is presented. This method is tested against a case of turbulent non-isothermal mixing in a T-junction pipe, a common ow arrangement found in nuclear reactor cooling systems. The reduced order model is derived from the 3D unsteady, incompressible Navier-Stokes equations weakly coupled with the energy equation.

We propose a computationally efficient framework to treat nonlinear partial differential equations having bifurcating solutions as one or more physical control parameters are varied. Our focus is on steady bifurcations. Plotting a bifurcation diagram entails computing multiple solutions of a parametrized, nonlinear problem, which can be extremely expensive in terms of computational time.

We introduce reduced order methods as an efficient strategy to solve parametrized non-linear and time dependent optimal flow control problems governed by partial differential equations. Indeed, optimal control problems require a huge computational effort in order to be solved, most of all in a physical and/or geometrical parametrized setting. Reduced order methods are a reliably suitable approach, increasingly gaining popularity, to achieve rapid and accurate optimal solutions in several fields, such as in biomedical and environmental sciences.

Reduced order modeling has gained considerable attention in recent decades owing to the advantages offered in reduced computational times and multiple solutions for parametric problems. The focus of this manuscript is the application of model order reduction techniques in various engineering and scientific applications including but not limited to mechanical, naval and aeronautical engineering. The focus here is kept limited to computational fluid mechanics and related applications.

In this contribution, we coupled the isogeometric analysis to a reduced order modelling technique in order to provide a computationally efficient solution in parametric domains. In details, we adopt the free-form deformation method to obtain the parametric formulation of the domain and proper orthogonal decomposition with interpolation for the computational reduction of the model.

In this work, Uncertainty Quantification (UQ) based on non-intrusive Polynomial Chaos Expansion (PCE) is applied to the CFD problem of the flow past an airfoil with parameterized angle of attack and inflow velocity. To limit the computational cost associated with each of the simulations required by the non-intrusive UQ algorithm used, we resort to a Reduced Order Model (ROM) based on Proper Orthogonal Decomposition (POD)-Galerkin approach.

We propose a computationally efficient framework to treat nonlinear partial differential equations having bifurcating solutions as one or more physical control parameters are varied. Our focus is on steady bifurcations. Plotting a bifurcation diagram entails computing multiple solutions of a parametrized, nonlinear problem, which can be extremely expensive in terms of computational time.

Reduced Order Models (ROMs), also known as Reduced Basis Methods (RBMs), have received considerable attention in recent years for their ability to drastically reduce CFD cost, particularly when dealing with parametrised problems in a multi-query setting.

This Special Issue gathers recent advances in ROM/RBM techniques for complex flow problems relevant to applications in mechanical and aerospace engineering, as well as medical and applied sciences.

This work presents a reduced order modeling technique built on a high fidelity embedded mesh finite element method. Such methods, and in particular the CutFEM method, are attractive in the generation of projection-based reduced order models thanks to their capabilities to seamlessly handle large deformations of parametrized domains and in general to handle topological changes.

In this work, Uncertainty Quantification (UQ) based on non-intrusive Polynomial Chaos Expansion (PCE) is applied to the CFD problem of the flow past an airfoil with parameterized angle of attack and inflow velocity. To limit the computational cost associated with each of the simulations required by the non-intrusive UQ algorithm used, we resort to a Reduced Order Model (ROM) based on Proper Orthogonal Decomposition (POD)-Galerkin approach.

Some euglenids, a family of aquatic unicellular organisms, can develop highly concerted, large-amplitude peristaltic body deformations. This remarkable behaviour has been known for centuries. Yet, its function remains controversial, and is even viewed as a functionless ancestral vestige. Here, by examining swimming Euglena gracilis in environments of controlled crowding and geometry, we show that this behaviour is triggered by confinement.

We propose a model order reduction technique integrating the Shifted Boundary Method (SBM) with a POD-Galerkin strategy. This approach allows to treat more complex parametrized domains in an efficient and straightforward way. The impact of the proposed approach is threefold.

First, problems involving parametrizations of complex geometrical shapes and/or large domain deformations can be efficiently solved at full-order by means of the SBM, an unfitted boundary method that avoids remeshing and the tedious handling of cut cells by introducing an approximate surrogate boundary.

Reduced-order modeling (ROM) commonly refers to the construction, based on a few solutions (referred to as snapshots) of an expensive discretized partial differential equation (PDE), and the subsequent application of low-dimensional discretizations of partial differential equations (PDEs) that can be used to more efficiently treat problems in control and optimization, uncertainty quantification, and other settings that require multiple approximate PDE solutions.

We consider a Leray model with a nonlinear differential low-pass filter for the simulation of incompressible fluid flow at moderately large Reynolds number (in the range of a few thousands) with under-refined meshes. For the implementation of the model, we adopt the three-step algorithm Evolve-Filter-Relax (EFR). The Leray model has been extensively applied within a Finite Element (FE) framework. Here, we propose to combine the EFR algorithm with a computationally efficient Finite Volume (FV) method.

A parametric Reduced Order Model (ROM) for buoyancy-driven flow is developed for which the Full Order Model (FOM) is based on the finite volume approximation and the Boussinesq approximation is used for modeling the buoyancy. Therefore, there exists a two-way coupling between the incompressible Boussinesq equations and the energy equation. The reduced basis is constructed with a Proper Orthogonal Decomposition (POD) approach and to obtain the Reduced Order Model, a Galerkin projection of the governing equations onto the reduced basis is performed.

Reduced order modeling (ROM) provides an efficient framework to compute solutions of parametric problems. Basically, it exploits a set of precomputed high-fidelity solutions—computed for properly chosen parameters, using a full-order model—in order to find the low dimensional space that contains the solution manifold. Using this space, an approximation of the numerical solution for new parameters can be computed in real-time response scenario, thanks to the reduced dimensionality of the problem.

In the reduced order modeling (ROM) framework, the solution of a parametric partial differential equation is approximated by combining the high-fidelity solutions of the problem at hand for several properly chosen configurations. Examples of the ROM application, in the naval field, can be found in [31, 24]. Mandatory ingredient for the ROM methods is the relation between the high-fidelity solutions and the parameters.

Marine propeller blade shape is constantly studied by engineers to obtain designs that allow for enhanced hydrodynamic performance while reducing vibrations and noise emissions. In such framework, shape parametrization and morphing algorithms are crucial elements of the numerical simulation and prototyping environment required for the evaluation of new blade geometries.

In this chapter we introduce a combined parameter and model reduction methodology and present its application to the efficient numerical estimation of a pressure drop in a set of deformed carotids. The aim is to simulate a wide range of possible occlusions after the bifurcation of the carotid. A parametric description of the admissible deformations, based on radial basis functions interpolation, is introduced. Since the parameter space may be very large, the first step in the combined reduction technique is to look for active subspaces in order to reduce the parameter space dimension.

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

In this work, we provide an integrated pipeline for the model-order reduction of turbulent flows around parametrised geometries in aerodynamics. In particular, free-form deformation is applied for geometry parametrisation, whereas two different reduced-order models based on proper orthogonal decomposition (POD) are employed in order to speed-up the full-order simulations: the first method exploits POD with interpolation, while the second one is based on domain decomposition.

In this work, we propose viable and eficient strategies for stabilized parametrized advection dominated problems, with random inputs. In particular, we investigate the combination of the wRB (weighted reduced basis) method for stochastic parametrized problems with the stabilized RB (reduced basis) method, which is the integration of classical stabilization methods (streamline/upwind Petrov-Galerkin (SUPG) in our case) in the ofine-online structure of the RB method.

An elastic bilayer, consisting of an equibiaxially pre-stretched sheet bonded to a stress-free one, spontaneously morphs into curved shapes in the absence of external loads or constraints. Using experiments and numerical simulations, we explore the role of geometry for square and rectangular samples in determining the equilibrium shape of the system, for a fixed pre-stretch. We classify the observed shapes over a wide range of aspect ratios according to their curvatures and compare measured and computed values, which show good agreement.

We present the results of the application of a parameter space reduction methodology based on active subspaces (AS) to the hull hydrodynamic design problem. Several parametric deformations of an initial hull shape are considered to assess the influence of the shape parameters on the hull wave resistance. Such problem is relevant at the preliminary stages of the ship design, when several flow simulations are carried out by the engineers to establish a certain sensibility with respect to the parameters, which might result in a high number of time consuming hydrodynamic simulations.

Shape optimization is a challenging task in many engineering fields, since the numerical solutions of parametric system may be computationally expensive. This work presents a novel optimization procedure based on reduced order modeling, applied to a naval hull design problem. The advantage introduced by this method is that the solution for a specific parameter can be expressed as the combination of few numerical solutions computed at properly chosen parametric points. The reduced model is built using the proper orthogonal decomposition with interpolation (PODI) method.

We propose reduced order methods as a suitable approach to face parametrized optimal control problems governed by partial differential equations, with applications in en- vironmental marine sciences and engineering. Environmental parametrized optimal control problems are usually studied for different configurations described by several physical and/or geometrical parameters representing different phenomena and structures. The solution of parametrized problems requires a demanding computational effort.

Vortex shedding around circular cylinders is a well known and studied phenomenon that appears in many engineering fields. In this work a Reduced Order Model (ROM) of the incompressible flow around a circular cylinder, built performing a Galerkin projection of the governing equations onto a lower dimensional space is presented. The reduced basis space is generated using a Proper Orthogonal Decomposition (POD) approach.

In this paper we apply a reduced basis framework for the computation of flow bifurcation (and stability) problems in fluid dynamics. The proposed method aims at reducing the complexity and the computational time required for the construction of bifurcation and stability diagrams. The method is quite general since it can in principle be specialized to a wide class of nonlinear problems, but in this work we focus on an application in incompressible fluid dynamics at low Reynolds numbers.

We focus on reducing the computational costs associated with the hydrodynamic stability of solutions of the incompressible Navier-Stokes equations for a Newtonian and viscous fluid in contraction-expansion channels. In particular, we are interested in studying steady bifurcations, occurring when non-unique stable solutions appear as physical and/or geometric control parameters are varied. The formulation of the stability problem requires solving an eigenvalue problem for a partial differential operator.

A fast computational framework is devised to the study of several configurations of patient-specific coronary artery bypass grafts. This is especially useful to perform a sensitivity analysis of the haemodynamics for different flow conditions occurring in native coronary arteries and bypass grafts, the investigation of the progression of the coronary artery disease and the choice of the most appropriate surgical procedure. A complete pipeline, from the acquisition of patientspecific medical images to fast parametrized computational simulations, is proposed.

Visualization and tracking of the facial whiskers is required in an increasing number of rodent studies. Though many approaches have been employed, only high-speed videography has proven adequate for measuring whisker motion and deformation during interaction with an object. However, whisker visualization and tracking is challenging for multiple reasons, primary among them the low contrast of the whisker against its background. Here we demonstrate a fluorescent dye method suitable for visualization of one or more rat whiskers.

In the control field, the study of the system dynamics is usually carried out relying on lumped-parameter or one-dimensional modelling. Even if these approaches are well suited for control purposes since they provide fast-running simulations and are easy to linearize, they may not be sufficient to deeply assess the complexity of the systems, in particular where spatial phenomena have a significant impact on dynamics. Reduced Order Methods (ROM) can offer the proper trade-off between computational cost and solution accuracy.

In this work we review a reduced basis method for the solution of uncertainty quantification problems. Based on the basic setting of an elliptic partial differential equation with random input, we introduce the key ingredients of the reduced basis method, including proper orthogonal decomposition and greedy algorithms for the construction of the reduced basis functions, a priori and a posteriori error estimates for the reduced basis approximations, as well as its computational advantages and weaknesses in comparison with a stochastic collocation method [I. Babuska, F. Nobile, and R.

In this work we provide a combination of isogeometric analysis with reduced order modelling techniques, based on proper orthogonal decomposition, to guarantee computational reduction for the numerical model, and with free-form deformation, for versatile geometrical parametrization. We apply it to computational fluid dynamics problems considering a Stokes flow model. The proposed reduced order model combines efficient shape deformation and accurate and stable velocity and pressure approximation for incompressible viscous flows, computed with a reduced order method.

In this work, two approaches, based on the certified Reduced Basis method, have been developed for simulating the movement of nuclear reactor control rods, in time-dependent non-coercive settings featuring a 3D geometrical framework. In particular, in a first approach, a piece-wise affine transformation based on subdomains division has been implemented for modelling the movement of one control rod. In the second approach, a staircase strategy has been adopted for simulating the movement of all the three rods featured by the nuclear reactor chosen as case study.

Several problems in applied sciences and engineering require reduction techniques in order to allow computational tools to be employed in the daily practice, especially in iterative procedures such as optimization or sensitivity analysis. Reduced order methods need to face increasingly complex problems in computational mechanics, especially into a multiphysics setting.

The aim of this work is to solve parametrized partial differential equations in computational domains represented by networks of repetitive geometries by combining reduced basis and domain decomposition techniques. The main idea behind this approach is to compute once, locally and for few reference shapes, some representative finite element solutions for different values of the parameters and with a set of different suitable boundary conditions on the boundaries: these functions will represent the basis of a reduced space where the global solution is sought for.

Even though capacity of modern computers keeps growing, it has become evident over the past decades that the complexity of practically relevant problems from science and engineering has grown even faster. Without highly efficient, reliable and robust numerical methods there will be no progress in many areas. The complexity of relevant problems in turn has lead to systems of equations with an extremely large number of unknowns that need to be handled numerically.

In this paper, we speculate on a possible application of Liquid Crystal Elastomers to the field of soft robotics. In particular, we study a concept for limbless locomotion that is amenable to miniaturisation. For this purpose, we formulate and solve the evolution equations for a strip of nematic elastomer, subject to directional frictional interactions with a flat solid substrate, and cyclically actuated by a spatially uniform, time-periodic stimulus (e.g., temperature change).

In this paper we develop and analyze a multilevel weighted reduced basis method for solving stochastic optimal control problems constrained by Stokes equations. We prove the analytic regularity of the optimal solution in the probability space under certain assumptions on the random input data. The finite element method and the stochastic collocation method are employed for the numerical approximation of the problem in the deterministic space and the probability space, respectively, resulting in many large-scale optimality systems to solve.

This paper extends the reduced basis method for the solution of parametrized optimal control problems presented in Negri et al. (2013) to the case of noncoercive (elliptic) equations, such as the Stokes equations. We discuss both the theoretical properties-with particular emphasis on the stability of the resulting double nested saddle-point problems and on aggregated error estimates-and the computational aspects of the method.

In this work, we present a stable proper orthogonal decomposition–Galerkin approximation for parametrized steady incompressible Navier–Stokes equations with low Reynolds number.

Shape optimization problems governed by PDEs result from many applications in computational fluid dynamics. These problems usually entail very large computational costs and require also a suitable approach for representing and deforming efficiently the shape of the underlying geometry, as well as for computing the shape gradient of the cost functional to be minimized. Several approaches based on the displacement of a set of control points have been developed in the last decades, such as the so-called free-form deformations.

We present some recent advances and improvements in shape parametrisation techniques of interfaces for reduced-order modelling with special attention to fluid–structure interaction problems and the management of structural deformations, namely, to represent them into a low-dimensional space (by control points). This allows to reduce the computational effort, and to significantly simplify the (geometrical) deformation procedure, leading to more efficient and fast reduced-order modelling applications in this kind of problems.

In this work, we propose viable and efficient strategies for the stabilization of the reduced basis approximation of an advection dominated problem. In particular, we investigate the combination of a classic stabilization method (SUPG) with the Offline-Online structure of the RB method. We explain why the stabilization is needed in both stages and we identify, analytically and numerically, which are the drawbacks of a stabilization performed only during the construction of the reduced basis (i.e. only in the Offline stage).

We present a method to generate a non-affine transfinite map from a given reference domain to a family of deformed domains. The map is a generalization of the Gordon-Hall transfinite interpolation approach. It is defined globally over the reference domain. Once we have computed some functions over the reference domain, the map can be generated by knowing the parametric expressions of the boundaries of the deformed domain. Being able to define a suitable map from a reference domain to a desired deformation is useful for the management of parameterized geometries.

In this paper we present a reduced order method for the solution of parametrized Stokes equations in domain composed by an arbitrary number of predefined shapes. The novelty of the proposed approach is the possibility to use a small set of precomputed bases to solve Stokes equations in very different computational domains, defined by combining one or more reference geometries. The selection of the basis functions is performed through an optimization greedy algorithm.

In this chapter we consider Reduced Basis (RB) approximations of parametrized Partial Differential Equations (PDEs). The the idea behind RB is to decouple the generation and projection stages (Offline/Online computational procedures) of the approximation process in order to solve parametrized PDEs in a fast, inexpensive and reliable way. The RB method, especially applied to 3D problems, allows great computational savings with respect to the classical Galerkin Finite Element (FE) Method.

We consider a method to efficiently evaluate in a real-time context an output based on the numerical solution of a partial differential equation depending on a large number of parameters. We state a result allowing to improve the computational performance of a three-step RB-ANOVA-RB method. This is a combination of the reduced basis (RB) method and the analysis of variations (ANOVA) expansion, aiming at compressing the parameter space without affecting the accuracy of the output.

We review the optimal design of an arterial bypass graft following either a (i) boundary optimal control approach, or a (ii) shape optimization formulation. The main focus is quantifying and treating the uncertainty in the residual flow when the hosting artery is not completely occluded, for which the worst-case in terms of recirculation effects is inferred to correspond to a strong orifice flow through near-complete occlusion.

This work deals with the development and application of reduction strategies for real-time and many query problems arising in fluid dynamics, such as shape optimization, shape registration (reconstruction), and shape parametrization. The proposed strategy is based on the coupling between reduced basis methods for the reduction of computational complexity and suitable shape parametrizations - such as free-form deformations or radial basis functions - for low-dimensional geometrical description.

In this work we propose and analyze a weighted reduced basis method to solve elliptic partial differential equations (PDEs) with random input data. The PDEs are first transformed into a weighted parametric elliptic problem depending on a finite number of parameters. Distinctive importance of the solution at different values of the parameters is taken into account by assigning different weights to the samples in the greedy sampling procedure. A priori convergence analysis is carried out by constructive approximation of the exact solution with respect to the weighted parameters.

We propose a suitable model reduction paradigm-the certified reduced basis method (RB)-for the rapid and reliable solution of parametrized optimal control problems governed by partial differential equations. In particular, we develop the methodology for parametrized quadratic optimization problems with elliptic equations as a constraint and infinite-dimensional control variable. First, we recast the optimal control problem in the framework of saddle-point problems in order to take advantage of the already developed RB theory for Stokes-type problems.