Inspired by the Equation-Free multiscale modeling approach, we demonstrate how the embed-learn-lift framework enables the construction of surrogate global normal-forms, namely minimal-dimensional reduced-order models (ROMs), from high-fidelity Navier-Stokes simulations. These surrogate models are then used for efficient and accurate bifurcation and stability analysis, thus dealing with the presence of continuous symmetries. The framework proceeds in four steps.

The current study aims to develop a non-intrusive Reduced Order Model (ROM) to reconstruct the full temperature field for a large-scale industrial application based on both numerical and experimental datasets. The proposed approach is validated against a domestic refrigerator.

In computational fluid dynamics, Large Eddy Simulation (LES) offers a compelling balance between accuracy and computational cost by resolving large-scale flow structures while modeling unresolved subgrid scales. However, its predictive capacity is critically dependent on the choice and calibration of subgrid-scale (SGS) models, which often involve problem-dependent parameters and exhibit intricate interactions with the numerical discretization.

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.

This work develops a rigorous numerical framework for solving time-dependent Optimal Control Problems (OCPs) governed by partial differential equations, with a particular focus on biomedical applications. The approach deals with adjoint-based Lagrangian methodology, which enables efficient gradient computation and systematic derivation of optimality conditions for both distributed and concentrated control formulations.

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

This work presents an overview of several nonlinear reduction strategies for data compression from various research fields, and a comparison of their performance when applied to problems characterized by diffusion and/or advection terms. We aim to create a common framework by unifying the notation referring to a common two-stage pipeline. At the same time, we underline their main differences and objectives by highlighting the diverse choices made for each stage.

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.

This work addresses the development of a physics-informed neural network (PINN) with a loss term derived from a discretized time-dependent full-order and reduced-order system.

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.

Numerical stabilization techniques are often employed in under-resolved simulations of convection-dominated flows to improve accuracy and mitigate spurious oscillations. Specifically, the evolve--filter--relax (EFR) algorithm is a framework which consists in evolving the solution, applying a filtering step to remove high-frequency noise, and relaxing through a convex combination of filtered and original solutions.

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.

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.

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.

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.

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.

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.

This chapter focuses on the combination of reduced order models and data-driven techniques applied to the study of turbulent flows in order to improve the pressure and velocity accuracy of standard reduced order methods. We focus on reduced order models constructed by means of Proper Orthogonal Decomposition with Galerkin approach, enhanced with two different stabilization techniques: (i) the supremizer enrichment approach, (ii) the pressure Poisson equation approach.

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.

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.

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.

This work discusses the correct modeling of the fully nonlinear free surface boundary conditions to be prescribed in water waves flow simulations based on potential flow theory. The main goal of such a discussion is that of identifying a mathematical formulation and a numerical treatment that can be used both to carry out transient simulations, and to compute steady solutions — for any flow admitting them. In the literature on numerical towing tank in fact, steady and unsteady fully nonlinear potential flow solvers are characterized by different mathematical formulations.

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 manuscript a POD-Galerkin based Reduced Order Model for unsteady Fluid-Structure Interaction problems is presented. The model is based on a partitioned algorithm, with semi-implicit treatment of the coupling conditions. A Chorin–Temam projection scheme is applied to the incompressible Navier–Stokes problem, and a Robin coupling condition is used for the coupling between the fluid and the solid. The coupled problem is based on an Arbitrary Lagrangian Eulerian formulation, and the Proper Orthogonal Decomposition procedure is used for the generation of the reduced basis.

As a major breakthrough in artificial intelligence and deep learning, Convolutional Neural Networks have achieved an impressive success in solving many problems in several fields including computer vision and image processing. Real-time performance, robustness of algorithms and fast training processes remain open problems in these contexts. In addition object recognition and detection are challenging tasks for resource-constrained embedded systems, commonly used in the industrial sector.

This contribution focuses on the development of Model Order Reduction (MOR) for one-way coupled steady state linear thermo-mechanical problems in a finite element setting. We apply Proper Orthogonal Decomposition (POD) for the computation of reduced basis space. On the other hand, for the evaluation of the modal coefficients, we use two different methodologies: the one based on the Galerkin projection (G) and the other one based on Artificial Neural Network (ANN). We aim to compare POD-G and POD-ANN in terms of relevant features including errors and computational efficiency.

Models with dominant advection always posed a difficult challenge for projection-based reduced order modelling. Many methodologies that have recently been proposed are based on the pre-processing of the full-order solutions to accelerate the Kolmogorov N−width decay thereby obtaining smaller linear subspaces with improved accuracy. These methods however must rely on the knowledge of the characteristic speeds in phase space of the solution, limiting their range of applicability to problems with explicit functional form for the advection field.

In this work we 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.

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.

Abstract This work explores the development and the analysis of an efficient reduced order model for the study of a bifurcating phenomenon, known as the Coand? effect, in a multi-physics setting involving fluid and solid media. Taking into consideration a fluid-structure interaction problem, we aim at generalizing previous works towards a more reliable description of the physics involved. In particular, we provide several insights on how the introduction of an elastic structure influences the bifurcating behavior.

We investigate various data-driven methods to enhance projection-based model reduction techniques with the aim of capturing bifurcating solutions. To show the effectiveness of the data-driven enhancements, we focus on the incompressible Navier-Stokes equations and different types of bifurcations. To recover solutions past a Hopf bifurcation, we propose an approach that combines proper orthogonal decomposition with Hankel dynamic mode decomposition. To approximate solutions close to a pitchfork bifurcation, we combine localized reduced models with artificial neural networks.

Reduced order modeling (ROM) techniques are numerical methods that approximate the solution of parametric partial differential equation (PDE) by properly combining the high-fidelity solutions of the problem obtained for several configurations, i.e. for several properly chosen values of the physical/geometrical parameters characterizing the problem. In this contribution, we propose an efficient non-intrusive data-driven framework involving ROM techniques in computational fluid dynamics (CFD) for hemodynamics applications.

Reduced basis approximations of Optimal Control Problems (OCPs) governed by steady partial differential equations (PDEs) with random parametric inputs are analyzed and constructed. Such approximations are based on a Reduced Order Model, which in this work is constructed using the method of weighted Proper Orthogonal Decomposition. This Reduced Order Model then is used to efficiently compute the reduced basis approximation for any outcome of the random parameter.

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 contribution describes the implementation of a data-driven shape optimization pipeline in a naval architecture application. We adopt reduced order models in order to improve the efficiency of the overall optimization, keeping a modular and equation-free nature to target the industrial demand. We applied the above mentioned pipeline to a realistic cruise ship in order to reduce the total drag. We begin by defining the design space, generated by deforming an initial shape in a parametric way using free form deformation.

In the field of parametric partial differential equations, shape optimization represents a challenging problem due to the required computational resources. In this contribution, a data-driven framework involving multiple reduction techniques is proposed to reduce such computational burden. Proper orthogonal decomposition (POD) and active subspace genetic algorithm (ASGA) are applied for a dimensional reduction of the original (high fidelity) model and for an efficient genetic optimization based on active subspace property.

The majority of the most common physical phenomena can be described using partial differential equations (PDEs). However, they are very often characterized by strong nonlinearities. Such features lead to the coexistence of multiple solutions studied by the bifurcation theory. Unfortunately, in practical scenarios, one has to exploit numerical methods to compute the solutions of systems of PDEs, even if the classical techniques are usually able to compute only a single solution for any value of a parameter when more branches exist.

It is well known in the Reduced Basis approximation of saddle point problems that the Galerkin projection on the reduced space does not guarantee the inf–sup approximation stability even if a stable high fidelity method was used to generate snapshots. For problems in computational fluid dynamics, the lack of inf–sup stability is reflected by the inability to accurately approximate the pressure field. In this context, inf–sup stability is usually recovered through the enrichment of the velocity space with suitable supremizer functions.

We investigate a projection-based reduced order model of the steady incompressible Navier–Stokes equations for moderate Reynolds numbers. In particular, we construct an “embedded” reduced basis space, by applying proper orthogonal decomposition to the Shifted Boundary Method, a high-fidelity embedded method recently developed. We focus on the geometrical parametrization through level-set geometries, using a fixed Cartesian background geometry and the associated mesh.

The present work focuses on the geometric parametrization and the reduced order modeling of the Stokes equation. We discuss the concept of a parametrized geometry and its application within a reduced order modeling technique. The full order model is based on the discontinuous Galerkin method with an interior penalty formulation. We introduce the broken Sobolev spaces as well as the weak formulation required for an affine parameter dependency. The operators are transformed from a fixed domain to a parameter dependent domain using the affine parameter dependency.

We provide first the functional analysis background required for reduced order modeling and present the underlying concepts of reduced basis model reduction. The projection-based model reduction framework under affinity assumptions, offline-online decomposition and error estimation is introduced. Several tools for geometry parametrizations, such as free form deformation, radial basis function interpolation and inverse distance weighting interpolation are explained.

We focus on steady and unsteady Navier–Stokes flow systems in a reduced-order modeling framework based on Proper Orthogonal Decomposition within a levelset geometry description and discretized by an unfitted mesh Finite Element Method. This work extends the approaches of [1], [2], [3] to nonlinear CutFEM discretization. We construct and investigate a unified and geometry independent reduced basis which overcomes many barriers and complications of the past, that may occur whenever geometrical morphings are taking place.

We discuss some recent results on biological and bio-inspired morphing, and use them to identify promising research directions for the future. In particular, we consider issues related to morphing at microscopic scales inspired by unicellular organisms. We focus on broad conceptual principles and, in particular, on morphing approaches based on the use of Gauss' theorema egregium (Gaussian morphing).

In this work, we study the mechanics of metamaterial sheets inspired by the pellicle of Euglenids. They are composed of interlocking elastic rods which can freely slide along their edges. We characterize the kinematics and the mechanics of these structures using the special Cosserat theory of rods and by assuming axisymmetric deformations of the tubular assembly. Through an asymptotic expansion, we investigate both structures that comprise a discrete number of rods and the limit case of a sheet composed by infinitely many rods.

A model order reduction technique is combined with an embedded boundary finite element method with a POD-Galerkin strategy. The proposed methodology is applied to parametrized heat transfer problems and we rely on a sufficiently refined shape-regular background mesh to account for parametrized geometries. In particular, the employed embedded boundary element method is the Shifted Boundary Method (SBM) recently proposed.

In this work, we present an approach for the efficient treatment of parametrized geometries in the context of POD-Galerkin reduced order methods based on Finite Volume full order approximations. On the contrary to what is normally done in the framework of finite element reduced order methods, different geometries are not mapped to a common reference domain: the method relies on basis functions defined on an average deformed configuration and makes use of the Discrete Empirical Interpolation Method (D-EIM) to handle together non-affinity of the parametrization and non-linearities.

The purpose of this work is to present different reduced order model strategies starting from full order simulations stabilized using a residual-based variational multiscale (VMS) approach. The focus is on flows with moderately high Reynolds numbers. The reduced order models (ROMs) presented in this manuscript are based on a POD-Galerkin approach. Two different reduced order models are presented, which differ on the stabilization used during the Galerkin projection. In the first case, the VMS stabilization method is used at both the full order and the reduced order levels.

Coronary artery bypass grafts (CABG) surgery is an invasive procedure performed to circumvent partial or complete blood flow blockage in coronary artery disease. In this work, we apply a numerical optimal flow control model to patient-specific geometries of CABG, reconstructed from clinical images of real-life surgical cases, in parameterized settings. The aim of these applications is to match known physiological data with numerical hemodynamics corresponding to different scenarios, arisen by tuning some parameters.

We present in this double contribution two different reduced order strategies for incompressible parameterized Navier-Stokes equations characterized by varying Reynolds numbers. The first strategy deals with low Reynolds number (laminar flow) and is based on a stabilized finite element method during the offline stage followed by a Galerkin projection on reduced basis spaces generated by a greedy algorithm. The second methodology is based on a full order finite volume discretization.

We consider the Navier-Stokes equations in a channel with a narrowing of varying height. The model is discretized with high-order spectral element ansatz functions, resulting in 6372 degrees of freedom. The steady-state snapshot solutions define a reduced order space through a standard POD procedure. The reduced order space allows to accurately and efficiently evaluate the steady-state solutions for different geometries. In particular, we detail different aspects of implementing the reduced order model in combination with a spectral element discretization.

In this work we deal with parametrized time dependent optimal control problems governed by partial differential equations. We aim at extending the standard saddle point framework of steady constraints to time dependent cases. We provide an analysis of the well-posedness of this formulation both for parametrized scalar parabolic constraint and Stokes governing equations and we propose reduced order methods as an effective strategy to solve them.

In this work we present a Reduced Order Model which is specifically designed to deal with turbulent flows in a finite volume setting. The method used to build the reduced order model is based on the idea of merging/combining projection-based techniques with data-driven reduction strategies. In particular, the work presents a mixed strategy that exploits a data-driven reduction method to approximate the eddy viscosity solution manifold and a classical POD-Galerkin projection approach for the velocity and the pressure fields, respectively.

The purpose of this work is to introduce a novel POD–Galerkin strategy for the semi-implicit hybrid high order finite volume/finite element solver introduced in Bermúdez et al. (2014) and Busto et al. (2018). The interest is into the incompressible Navier–Stokes equations coupled with an additional transport equation. The full order model employed in this article makes use of staggered meshes. This feature will be conveyed to the reduced order model leading to the definition of reduced basis spaces in both meshes.

This paper deals with a concept for a reconfigurable structure bio-inspired by the cell wall architecture of euglenids, a family of unicellular protists, and based on the relative sliding of adjacent strips. Uniform sliding turns a cylinder resulting from the assembly of straight and parallel strips into a cylinder of smaller height and larger radius, in which the strips are deformed into a family of parallel helices.

The purpose of this work is to introduce a novel POD-Galerkin strategy for the hybrid finite volume/finite element solver introduced in Bermúdez et al. 2014 and Busto et al. 2018. The interest is into the incompressible Navier-Stokes equations coupled with an additional transport equation. The full order model employed in this article makes use of staggered meshes. This feature will be conveyed to the reduced order model leading to the definition of reduced basis spaces in both meshes.

This work focuses on the computationally efficient detection of the buckling phenomena and bifurcation analysis of the parametric Von Kármán plate equations based on reduced order methods and spectral analysis. The computational complexity - due to the fourth order derivative terms, the non-linearity and the parameter dependence - provides an interesting benchmark to test the importance of the reduction strategies, during the construction of the bifurcation diagram by varying the parameter(s).

We consider the Navier-Stokes equations in a channel with varying Reynolds numbers. The model is discretized with high-order spectral element ansatz functions, resulting in 14 259 degrees of freedom. The steady-state snapshot solu- tions define a reduced order space, which allows to accurately evaluate the steady- state solutions for varying Reynolds number with a reduced order model within a fixed-point iteration. In particular, we compare different aspects of implementing the reduced order model with respect to the use of a spectral element discretization.

A parametric reduced order model based on proper orthogonal decom- position with Galerkin projection has been developed and applied for the modeling of heat transport in T-junction pipes which are widely found in nuclear power plants. Thermal mixing of different temperature coolants in T-junction pipes leads to tem- perature fluctuations and this could potentially cause thermal fatigue in the pipe walls.

In this work we propose and analyze a weighted proper orthogonal decomposition method to solve elliptic partial differential equations depending on random input data, for stochastic problems that can be transformed into parametric systems. The algorithm is introduced alongside the weighted greedy method. Our proposed method aims to minimize the error in a L2 norm and, in contrast to the weighted greedy approach, it does not require the availability of an error bound.

We propose a numerical pipeline for shape optimization in naval engineering involving two different non-intrusive reduced order method (ROM) techniques. Such methods are proper orthogonal decomposition with interpolation (PODI) and dynamic mode decomposition (DMD). The ROM proposed will be enhanced by active subspaces (AS) as a pre-processing tool that reduce the parameter space dimension and suggest better sampling of the input space. We will focus on geometrical parameters describing the perturbation of a reference bulbous bow through the free form deformation (FFD) technique.

In this manuscript we discuss weighted reduced order methods for stochastic partial differential equations. Random inputs (such as forcing terms, equation coefficients, boundary conditions) are considered as parameters of the equations. We take advantage of the resulting parametrized formulation to propose an efficient reduced order model; we also profit by the underlying stochastic assumption in the definition of suitable weights to drive to reduction process.

Efficient shape morphing techniques play a crucial role in the approximation of partial differential equations defined in parametrized domains, such as for fluid‐structure interaction or shape optimization problems. In this paper, we focus on inverse distance weighting (IDW) interpolation techniques, where a reference domain is morphed into a deformed one via the displacement of a set of control points. We aim at reducing the computational burden characterizing a standard IDW approach without significantly compromising the accuracy. To this aim, first we propose an improvement of IDW based on a geometric criterion that automatically selects a subset of the original set of control points. Then, we combine this new approach with a dimensionality reduction technique based on a proper orthogonal decomposition of the set of admissible displacements. This choice further reduces computational costs. We verify the performances of the new IDW techniques on several tests by investigating the trade‐off reached in terms of accuracy and efficiency.

We present the results of the first application in the naval architecture field of a methodology based on active subspaces properties for parameters space reduction. The physical problem considered is the one of the simulation of the hydrodynamic flow past the hull of a ship advancing in calm water.

In this work we consider (hierarchical, Lagrange) reduced basis approximation and a posteriori error estimation for elasticity problems in affinely parametrized geometries.

In 2010 IMO (International Maritime Organisation) introduced new rules in SOLAS with the aim of intrinsically increase the safety of passenger ships. This requirement is achieved by providing safe areas for passengers and essential services for allowing ship to Safely Return to Port (SRtP). The entry into force of these rules has changed the way to design passenger ships. In this respect big effort in the research has been done by industry to address design issues related to the impact on failure analysis of the complex interactions among systems.

Flutter instability in elastic structures subject to follower load, the most important cases being the famous Beck's and Pflüger's columns (two elastic rods in a cantilever configuration, with an additional concentrated mass at the end of the rod in the latter case), have attracted, and still attract, a thorough research interest. In this field, the most important issue is the validation of the model itself of follower force, a nonconservative action which was harshly criticized and never realized in practice for structures with diffused elasticity.

In this contribution, we present the results of the application of a parameter space reduction methodology based on active subspaces to the hull hydrodynamic design problem. Several parametric deformations of an initial hull shape are considered to assess the influence of the shape parameters considered on the hull total drag. The hull resistance is typically computed by means of numerical simulations of the hydrodynamic flow past the ship.

We present a model order reduction approach for parametrized laminar flow problems including viscous boundary layers. The viscous effects are captured by the incompressible Navier–Stokes equations in the vicinity of the boundary layer, whereas a potential flow model is used in the outer region. By this, we provide an accurate model that avoids imposing the Kutta condition for potential flows as well as an expensive numerical solution of a global viscous model. To account for the parametrized nature of the problem, we apply the reduced basis method.

The introduction of computational models in cardiovascular sciences has been progressively bringing new and unique tools for the investigation of the physiopathology. Together with the dramatic improvement of imaging and measuring devices on one side, and of computational architectures on the other one, mathematical and numerical models have provided a new, clearly noninvasive, approach for understanding not only basic mechanisms but also patient-specific conditions, and for supporting the design and the development of new therapeutic options.

This chapter presents an overview of model order reduction – a new paradigm in the field of simulation-based engineering sciences, and one that can tackle the challenges and leverage the opportunities of modern ICT technologies. Despite the impressive progress attained by simulation capabilities and techniques, a number of challenging problems remain intractable. These problems are of different nature, but are common to many branches of science and engineering.

POD–Galerkin reduced-order models (ROMs) for fluid-structure interaction problems (incompressible fluid and thin structure) are proposed in this paper. Both the high-fidelity and reduced-order methods are based on a Chorin-Temam operator-splitting approach. Two different reduced-order methods are proposed, which differ on velocity continuity condition, imposed weakly or strongly, respectively. The resulting ROMs are tested and compared on a representative haemodynamics test case characterized by wave propagation, in order to assess the capabilities of the proposed strategies.

In this work we apply reduced basis methods for parametric PDEs to an isogeometric formulation based on NURBS. The motivation for this work is an integrated and complete work pipeline from CAD to parametrization of domain geometry, then from full order to certified reduced basis solution. IsoGeometric Analysis (IGA) is a growing research theme in scientific computing and computational mechanics, as well as reduced basis methods for parametric PDEs.

Brittle materials fail catastrophically. In consequence of their limited flaw-tolerance, failure occurs by localized fracture and is typically a dynamic process. Recently, experiments on epithelial cell monolayers have revealed that this scenario can be significantly modified when the material susceptible to cracking is adhered to a hydrogel substrate. Thanks to the hydraulic coupling between the brittle layer and the poroelastic substrate, such a composite can develop a toughening mechanism that relies on the simultaneous growth of multiple cracks.

The flagellar swimming of euglenids, which are propelled by a single anterior flagellum, is characterized by a generalized helical motion. The 3D nature of this swimming motion, which lacks some of the symmetries enjoyed by more common model systems, and the complex flagellar beating shapes that power it make its quantitative description challenging. In this work, we provide a quantitative, 3D, highly resolved reconstruction of the swimming trajectories and flagellar shapes of specimens of Euglena gracilis.

In this work we present a reduced basis Smagorinsky turbulence model for steady flows. We approximate the non-linear eddy diffusion term using the Empirical Interpolation Method, and the velocity-pressure unknowns by an independent reduced-basis procedure. This model is based upon an a posteriori error estimation for Smagorinsky turbulence model. The theoretical development of the a posteriori error estimation is based on previous works, according to the Brezzi-Rappaz-Raviart stability theory, and adapted for the non-linear eddy diffusion term.

Numerical simulation of fluid flows requires important computational efforts but it is essential in engineering applications. Reduced Order Model (ROM) can be employed whenever fast simulations are required, or in general, whenever a trade-off between computational cost and solution accuracy is a preeminent issue as in process optimization and control.

We explore the Mode I fracture toughness of a polymer gel containing a semi-infinite, growing crack. First, an expression is derived for the energy release rate within the linearized, small-strain setting. This expression reveals a crack tip velocity-independent toughening that stems from the poroelastic nature of polymer gels. Then, we establish a poroelastic cohesive zone model that allows us to describe the micromechanics of fracture in gels by identifying the role of solvent pressure in promoting poroelastic toughening.

In this work a reduced-order computational framework for the study of haemodynamics in three-dimensional patient-specific configurations of coronary artery bypass grafts dealing with a wide range of scenarios is proposed. We combine several efficient algorithms to face at the same time both the geometrical complexity involved in the description of the vascular network and the huge computational cost entailed by time dependent patient-specific flow simulations. Medical imaging procedures allow to reconstruct patient-specific configurations from clinical data.

In this paper we propose a monolithic approach for reduced order modelling of parametrized fluid-structure interaction problems based on a proper orthogonal decomposition (POD)–Galerkin method. Parameters of the problem are related to constitutive properties of the fluid or structural problem, or to geometrical parameters related to the domain configuration at the initial time.

In this work, a Reduced Basis method, with basis functions sampled by a Proper Orthogonal Decomposition technique, has been employed to develop a reduced order model of a multi-physics parametrized Lead-cooled Fast Reactor single-channel. Being the first time that a reduced order model is developed in this context, the work focused on a methodological approach and the coupling between the neutronics and the heat transfer, where the thermal feedbacks on neutronics are explicitly taken into account, in time-invariant settings.

The coupling of a free flow with a flow through porous media has many potential applications in several fields related with computational science and engineering, such as blood flows, environmental problems or food technologies. We present a reduced basis method for such coupled problems. The reduced basis method is a model order reduction method applied in the context of parametrized systems. Our approach is based on a heterogeneous domain decomposition formulation, namely the Stokes-Darcy problem. Thanks to an offline/online-decomposition, computational times can be drastically reduced.

Brittle materials propagate opening cracks under tension. When stress increases beyond a critical magnitude, then quasistatic crack propagation becomes unstable. In the presence of several precracks, a brittle material always propagates only the weakest crack, leading to catastrophic failure. Here, we show that all these features of brittle fracture are fundamentally modified when the material susceptible to cracking is bonded to a hydrogel, a common situation in biological tissues.

In many engineering applications, the investigation of the vibro-acoustic response of structures is of great interest. Hence, great effort has been dedicated to improve methods in this field in the last twenty years. Classical techniques have the main drawback that they become unaffordable when high frequency impact waves are considered. In that sense, the Energy Finite Element Analysis (EFEA) is a good alternative to those methods. Based on an approximate model, EFEA gives time and locally space-averaged energy densities and has been proven to yield accurate results.

In this work we show some results about the reduced basis approximation of advection dominated parametrized problems, i.e. advection-diffusion problems with high Péclet number. These problems are of great importance in several engineering applications and it is well known that their numerical approximation can be affected by instability phenomena. In this work we compare two possible stabilization strategies in the framework of the reduced basis method, by showing numerical results obtained for a steady advection-diffusion problem.

In this work, a Reduced Order Model (ROM) for multigroup time-dependent parametrized reactor spatial kinetics is presented. The Reduced Basis method (built upon a high-fidelity ``truth'' finite element approximation) has been applied to model the neutronics behavior of a parametrized system composed by a control rod surrounded by fissile material. The neutron kinetics has been described by means of a parametrized multi-group diffusion equation where the height of the control rod (i.e., how much the rod is inserted) plays the role of the varying parameter.

In this paper, two modelling approaches based on a Modal Method (MM) and on the Proper Orthogonal Decomposition (POD) technique, for developing a control-oriented model of nuclear reactor spatial kinetics, are presented and compared. Both these methods allow developing neutronics description by means of a set of ordinary differential equations. The comparison of the outcomes provided by the two approaches focuses on the capability of evaluating the reactivity and the neutron flux shape in different reactor configurations, with reference to a TRIGA Mark II reactor.

We extend the classical empirical interpolation method [M. Barrault, Y. Maday, N.C. Nguyen and A.T. Patera, An empirical interpolation method: application to efficient reduced-basis discretization of partial differential equations. Compt. Rend. Math. Anal. Num. 339 (2004) 667-672] to a weighted empirical interpolation method in order to approximate nonlinear parametric functions with weighted parameters, e.g. random variables obeying various probability distributions.

The stochastic collocation method (Babuška et al. in SIAM J Numer Anal 45(3):1005-1034, 2007; Nobile et al. in SIAM J Numer Anal 46(5):2411-2442, 2008a; SIAM J Numer Anal 46(5):2309-2345, 2008b; Xiu and Hesthaven in SIAM J Sci Comput 27(3):1118-1139, 2005) has recently been applied to stochastic problems that can be transformed into parametric systems. Meanwhile, the reduced basis method (Maday et al.

Advection-dominated problems, which arise in many engineering situations, often require a fast and reliable approximation of the solution given some parameters as inputs. In this work we want to investigate the coupling of the reduced basis method - which guarantees rapidity and reliability - with some classical stabilization techiques to deal with the advection-dominated condition. We provide a numerical extension of the results presented in [1], focusing in particular on problems with curved boundary layers and inner fronts whose direction depends on the parameter.

This chapter reviews techniques of model reduction of fluid dynamics systems. Fluid systems are known to be difficult to reduce efficiently due to several reasons. First of all, they exhibit strong nonlinearities - which are mainly related either to nonlinear convection terms and/or some geometric variability - that often cannot be treated by simple linearization.

In this paper we review and we extend the reduced basis approximation and a posteriori error estimation for steady Stokes flows in affinely parametrized geometries, focusing on the role played by the Brezzi's and Babuška's stability constants. The crucial ingredients of the methodology are a Galerkin projection onto a low-dimensional space of basis functions properly selected, an affine parametric dependence enabling to perform competitive Offline-Online splitting in the computational procedure and a rigorous a posteriori error estimation on field variables.

The solution of inverse problems in cardiovascular mathematics is computationally expensive. In this paper, we apply a domain parametrization technique to reduce both the geometrical and computational complexities of the forward problem and replace the finite element solution of the incompressible Navier-Stokes equations by a computationally less-expensive reduced-basis approximation. This greatly reduces the cost of simulating the forward problem.

This work aims at identifying and quantifying uncertainties from various sources in human cardiovascular system based on stochastic simulation of a one-dimensional arterial network. A general analysis of different uncertainties and probability characterization with log-normal distribution of these uncertainties is introduced.

The set of solutions of a parameter-dependent linear partial differential equation with smooth coefficients typically forms a compact manifold in a Hilbert space. In this paper we review the generalized reduced basis method as a fast computational tool for the uniform approximation of the solution manifold We focus on operators showing an affine parametric dependence, expressed as a linear combination of parameter-independent operators through some smooth, parameter-dependent scalar functions.

In this work we deal with a stochastic optimal Robin boundary control problem constrained by an advection-diffusion-reaction elliptic equation with advection-dominated term. We assume that the uncertainty comes from the advection field and consider a stochastic Robin boundary condition as control function.