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.

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.

Reduced Order Models (ROMs) have been regarded as an efficient alternative to conventional high-fidelity Computational Fluid Dynamics (CFD) for accelerating the design and optimization processes in engineering applications. Many industrial geometries feature repeating subdomains or contain sub-regions governed by distinct physical phenomena, making them well-suited to Domain Decomposition (DD) techniques. The integration of ROM and DD is promising to further reduce computational costs by constructing local ROMs and assembling them into global solutions.

This special issue, dedicated to the inaugural edition of the ECCOMAS-IACM conference series Math 2 Product (M2P), contains eleven articles on frontier topics in Computational Science and Engineering, with a specific focus on next-generation methodologies for industrial problems and sustainability.

The complexity of the cardiovascular system needs to be accurately reproduced in order to promptly acknowledge health conditions; to this aim, advanced multifidelity and multiphysics numerical models are crucial. On one side, Full Order Models (FOMs) deliver accurate hemodynamic assessments, but their high computational demands hinder their real-time clinical application. In contrast, Reduced Order Models (ROMs) provide more efficient yet accurate solutions, essential for personalized healthcare and timely clinical decision-making.

Graph autoencoders have gained attention in nonlinear reduced-order modeling of parameterized partial differential equations defined on unstructured grids. Despite they provide a geometrically consistent way of treating complex domains, applying such architectures to parameterized dynamical systems for temporal prediction beyond the training data, i.e. the extrapolation regime, is still a challenging task due to the simultaneous need of temporal causality and generalizability in the parametric space.

In collaboration with industrial partners, this research leverages the concept of a “Digital Twin” to create a virtual replica of industrial machinery that monitors health, identifies anomalies, and predicts potential failures. Analysing time-series data from multiple sensors poses significant challenges, particularly as machines dynamically adjust their operating conditions to meet production demands, making traditional forecasting algorithms ineffective.

This work investigates model reduction techniques for nonlinear parameterized and time-dependent PDEs, specifically focusing on bifurcating phenomena in Computational Fluid Dynamics (CFD). We develop interpretable and non-intrusive Reduced Order Models (ROMs) capable of capturing dynamics associated with bifurcations by identifying a minimal set of coordinates. Our methodology combines the Sparse Identification of Nonlinear Dynamics (SINDy) method with a deep learning framework based on Autoencoder (AE) architectures.

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.

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.

Continuous casting is a widely adopted process in the steel industry, where maintaining high steel quality is paramount. Efficient prediction of grade intermixing during ladle changeover operations is critical for maintaining steel quality and minimizing material losses in the continuous casting process. Among various factors influencing grade intermixing, operating parameters play a significant role, in addition to tundish geometry and flow control devices.

Advection-dominated problems are predominantly noticed in nature, engineering systems, and various industrial processes. Traditional linear compression methods, such as proper orthogonal decomposition (POD) and reduced basis (RB) methods are ill-suited for these problems, due to slow Kolmogorov n-width decay. This results in inefficient and inaccurate reduced order models (ROMs). There are few non-linear approaches to accelerate the Kolmogorov n-width decay. In this work, we use a neural network shift augmented transformation technique that employs automatic shift detection.

This paper deals with the development of a Reduced-Order Model (ROM) to investigate haemodynamics in cardiovascular applications. It employs the use of Proper Orthogonal Decomposition (POD) for the computation of the basis functions and the Galerkin projection for the computation of the reduced coefficients.

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.

Information propagation characterizes how input correlations evolve across layers in deep neural networks. This framework has been well studied using mean-field theory, which assumes infinitely wide networks. However, these assumptions break down for practical, finite-size networks. In this work, we study information propagation in randomly initialized neural networks with finite width and reveal that the boundary between ordered and chaotic regimes exhibits a fractal structure.

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.

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.

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.

Forecasting atmospheric flows with traditional discretization methods, also called full order methods (e.g., finite element methods or finite volume methods), is computationally expensive. We propose to reduce the computational cost with a Reduced Order Model (ROM) that combines Extended Convolutional Autoencoders (E-CAE) and Reservoir Computing (RC).

It is well known that in the computational fluid dynamics simulations related to the cardiovascular system the enforcement of outflow boundary conditions is a crucial point. In fact, they highly affect the computed flow and a wrong setup could lead to unphysical results. In this chapter we discuss the main features of two different ways for the estimation of proper outlet boundary conditions in the context of hemodynamics simulations: on one side, a lumped parameter model of the downstream circulation and, on the other one, a technique based on optimal control.

Although the two-layer quasi-geostrophic equations (2QGE) are a simplified model for the dynamics of a stratified, wind-driven ocean, their numerical simulation is still plagued by the need for high resolution to capture the full spectrum of turbulent scales. Since such high resolution would lead to unreasonable computational times, it is typical to resort to coarse low-resolution meshes combined with the so-called eddy viscosity parameterization to account for the diffusion mechanisms that are not captured due to mesh under-resolution.

This study introduces a first step for constructing a hybrid reduced-order models (ROMs) for segregated fluid-structure interaction in an Arbitrary Lagrangian-Eulerian (ALE) approach at a high Reynolds number using the Finite Volume Method (FVM). The ROM is driven by proper orthogonal decomposition (POD) with hybrid techniques that combines the classical Galerkin projection and two data-driven methods (radial basis networks , and neural networks/ long short term memory). Results demonstrate the ROM ability to accurately capture the physics of fluid-structure interaction phenomena.

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

Patient-specific modeling of cardiovascular flows with high-fidelity is challenging due to its dependence on accurately estimated velocity boundary profiles, which are essential for precise simulations and directly influence wall shear stress calculations - key in predicting cardiovascular diseases like atherosclerosis. This data, often derived from in vivo modalities like 4D flow MRI, suffers from low resolution and noise.

Reduced Order Models (ROMs) have gained a great attention by the scientific community in the last years thanks to their capabilities of significantly reducing the computational cost of the numerical simulations, which is a crucial objective in applications like real time control and shape optimization. This contribution aims to provide a brief overview about such a topic.

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.

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.

The goal of this paper is to test solids4Foam, the fluid-structure interaction (FSI) toolbox developed for foam-extend (a branch of OpenFOAM), and assess its flexibility in handling more complex flows. For this purpose, we consider the interaction of an incompressible fluid described by a Leray model with a hyperelastic structure modeled as a Saint Venant-Kirchhoff material. We focus on a strongly coupled, partitioned fluid-structure interaction (FSI) solver in a finite volume environment, combined with an arbitrary Lagrangian-Eulerian approach to deal with the motion of the fluid domain.

Parallel to the need for new technologies and renewable energy resources to address sustainability, the emerging field of Artificial Intelligence (AI) has experienced continuous high-speed growth in the application of its capabilities of modelling, managing, processing, and making sense of data in the entire areas related to the production and management of energy.

In this work, we present GAROM, a new approach for reduced order modeling (ROM) based on generative adversarial networks (GANs). GANs attempt to learn to generate data with the same statistics of the underlying distribution of a dataset, using two neural networks, namely discriminator and generator. While widely applied in many areas of deep learning, little research is done on their application for ROM, i.e. approximating a high-fidelity model with a simpler one.

In this article, we propose a shape optimization pipeline for propeller blades, applied to naval applications. The geometrical features of a blade are exploited to parametrize it, allowing to obtain deformed blades by perturbating their parameters. The optimization is performed using a genetic algorithm that exploits the computational speed-up of reduced order models to maximize the efficiency of a given propeller. A standard offline–online procedure is exploited to construct the reduced-order model.

In this work, we address parametric non–stationary fluid dynamics problems within a model order reduction setting based on domain decomposition. Starting from the optimisation–based domain decomposition approach, we derive an optimal control problem, for which we present a convergence analysis in the case of non–stationary incompressible Navier–Stokes equations. We discretise the problem with the finite element method and we compare different model order reduction techniques: POD–Galerkin and a non–intrusive neural network procedures.

In this work, we address parametric non–stationary fluid dynamics problems within a model order reduction setting based on domain decomposition. Starting from the optimisation–based domain decomposition approach, we derive an optimal control problem, for which we present a convergence analysis in the case of non–stationary incompressible Navier–Stokes equations. We discretise the problem with the finite element method and we compare different model order reduction techniques: POD–Galerkin and a non–intrusive neural network procedures.

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.

Dynamic mode decomposition (DMD) has recently become a popular tool for the non-intrusive analysis of dynamical systems. Exploiting Proper Orthogonal Decomposition (POD) as a dimensionality reduction technique, DMD is able to approximate a dynamical system as a sum of spatial basis evolving linearly in time, thus enabling a better understanding of the physical phenomena and forecasting of future time instants. In this work we propose an extension of DMD to parameterized dynamical systems, focusing on the future forecasting of the output of interest in a parametric context.

In this paper, we propose hybrid data-driven ROM closures for fluid flows. These new ROM closures combine two fundamentally different strategies: (i) purely data-driven ROM closures, both for the velocity and the pressure; and (ii) physically based, eddy viscosity data-driven closures, which model the energy transfer in the system. The first strategy consists in the addition of closure/correction terms to the governing equations, which are built from the available data.

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.

Nonlinear extensions to the active subspaces method have brought remarkable results for dimension reduction in the parameter space and response surface design. We further develop a kernel-based nonlinear method. In particular, we introduce it in a broader mathematical framework that contemplates also the reduction in parameter space of multivariate objective functions. The implementation is thoroughly discussed and tested on more challenging benchmarks than the ones already present in the literature, for which dimension reduction with active subspaces produces already good results.

Left ventricular assist devices (LVADs) are used to provide haemodynamic support to patients with critical cardiac failure. Severe complications can occur because of the modifications of the blood flow in the aortic region. In this work, the effect of a continuous flow LVAD device on the aortic flow is investigated by means of a non-intrusive reduced order model (ROM) built using the proper orthogonal decomposition with interpolation (PODI) method based on radial basis functions (RBF).

In this contribution we propose reduced order methods to fast and reliably solve parametrized optimal control problems governed by time dependent nonlinear partial differential equations. Our goal is to provide a tool to deal with the time evolution of several nonlinear optimality systems in many-query context, where a system must be analysed for various physical and geometrical features. Optimal control can be used in order to fill the gap between collected data and mathematical model and it is usually related to very time consuming activities: inverse problems, statistics, etc.

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.

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.

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 develop a Reduced Order Model (ROM) for a Large Eddy Simulation (LES) approach that combines a three-step algorithm called Evolve-Filter-Relax (EFR) with a computationally efficient finite volume method. The main novelty of our ROM lies in the use within the EFR algorithm of a nonlinear, deconvolution-based indicator function that identifies the regions of the domain where the flow needs regularization.

A novel reduced order model (ROM) for incompressible flows is developed by performing a Galerkin projection based on a fully (space and time) discrete full order model (FOM) formulation. This 'discretize-then-project' approach requires no pressure stabilization technique (even though the pressure term is present in the ROM) nor a boundary control technique (to impose the boundary conditions at the ROM level). These are two main advantages compared to existing approaches.

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.

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.

In this work, we present an approach for the efficient treatment of parametrized geometries in the context of proper orthogonal decomposition (POD)-Galerkin reduced order methods based on finite-volume full order approximations.

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.

We consider the Navier–Stokes equations in a channel with a narrowing and walls of varying curvature. By applying the empirical interpolation method to generate an affine parameter dependency, the offline-online procedure can be used to compute reduced order solutions for parameter variations. The reduced-order space is computed from the steady-state snapshot solutions by a standard POD procedure. The model is discretised with high-order spectral element ansatz functions, resulting in 4752 degrees of freedom.

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.

Computational models are bringing novel perspectives into the investigation of cardiovascular issues. Together with imaging improvements, computational tools are increasingly adopted toward the understanding of patient-specific pathological situations, in the development of surgical planning, and in the support of interventional medical decisions. However, the complexity of patient-specific problems calls for an adequate and knowledgeable use of computational tools to obtain valuable results.

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.

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.

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.

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.

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

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.

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.

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.

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.

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.