Researcher profile

Gianluigi Rozza

Gianluigi Rozza contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
34works
0followers
10topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

34 published item(s)

preprint2026arXiv

An efficient hyper reduced-order model for segregated solvers for geometrical parametrization problems

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. The method is assessed on three benchmark problems: a linear transport equation, a nonlinear Burgers equation, and the incompressible Navier--Stokes equations. The results show that the hyper-reduced models closely match full-order solutions while achieving substantial reductions in computational time. Since only a sparse subset of mesh cells is evaluated during the online phase, the method is naturally parallelizable and scalable to very large meshes. These findings demonstrate that hyper-reduction can be effectively combined with segregated solvers and geometric parametrization to enable fast and accurate CFD simulations.

preprint2026arXiv

Data-driven reduced order model for residence time distribution analysis of an industrial-scale continuous casting tundish

The continuous casting tundish plays a critical role as a metallurgical reactor in the continuous casting process, with its flow characteristics serving as a key parameter in the production of high-quality steel. These characteristics are typically assessed through residence time distribution (RTD) curves. This study examines the flow behaviour in a single-strand continuous casting tundish through a combination of numerical simulations and experimental validation. Steady-state full order model (FOM) simulations are performed under both isothermal and non-isothermal conditions to evaluate the influence of thermal buoyancy on the velocity field, which is found to be negligible. The resulting flow fields are used to initialize transient tracer transport simulations for determining the RTD and flow volume partitioning. Subsequently, a data-driven reduced order model (ROM) is developed to predict the RTD response. Comparison of RTD curves obtained from experiments, FOM, and ROM shows excellent agreement, with the ROM accurately capturing the key flow characteristics at a fraction of the computational cost. These results highlight the potential of ROM techniques for efficient real-time analysis, design, and optimization of tundish operations in metallurgical processes.

preprint2026arXiv

Introduction to the Combination of Reduced Order Models and Domain Decomposition: State of the Art and Perspectives

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. Due to the complexity and necessity of coupling ROMs, many approaches have been proposed in recent years. This review provides a concise overview of existing methodologies combining ROM and DD. We categorize existing methods into intrusive (projection-based) and non-intrusive (data-driven) frameworks. Various strategies for generating local reduced bases and coupling them across subdomains are illustrated. Particular emphasis is placed on intrusive techniques, including equations, numerical algorithms, and practical implementations. The non-intrusive framework is also discussed, highlighting its general procedures, basic formulations, and underlying principles. Finally, we summarise the state of the literature, identify open challenges, and present perspectives on future implementation from an engineering viewpoint.

preprint2026arXiv

Latent Dynamics Graph Convolutional Networks for model order reduction of parameterized time-dependent PDEs

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. In this work, we address this gap by introducing Latent Dynamics Graph Convolutional Network (LD-GCN), a purely data-driven, encoder-free architecture that learns a global, low-dimensional representation of dynamical systems conditioned on external inputs and parameters. The temporal evolution is modeled in the latent space and advanced through time-stepping, allowing for time-extrapolation, and the trajectories are consistently decoded onto geometrically parameterized domains using a GNN. Our framework enhances interpretability by enabling the analysis of the reduced dynamics and supporting zero-shot prediction through latent interpolation. The methodology is mathematically validated via a universal approximation theorem for encoder-free architectures, and numerically tested on complex computational mechanics problems involving physical and geometric parameters, including the detection of bifurcating phenomena for Navier-Stokes equations. Code availability: https://github.com/lorenzotomada/ld-gcn-rom

preprint2026arXiv

Revisiting Deep Information Propagation: Fractal Frontier and Finite-size Effects

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 shows the fundamental complexity of neural network dynamics, in a setting that is independent of input data and optimization. To extend this analysis beyond multilayer perceptrons, we leverage recently introduced Fourier-based structured transforms, and show that information propagation in convolutional neural networks also follow the same behavior. In practice, our investigation highlights the importance of finite network depth with respect to the tradeoff between separation and robustness.

preprint2022arXiv

A data-driven partitioned approach for the resolution of time-dependent optimal control problems with dynamic mode decomposition

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. Thus, after the Dynamic Mode Decomposition operators construction, we reconstruct and perform future predictions for all the variables of interest at a lower computational cost with respect to the standard space-time discretized models. We test the methodology in terms of relative reconstruction and prediction errors on a boundary control for a Graetz flow and on a distributed control with Stokes constraints.

preprint2022arXiv

A data-driven Reduced Order Method for parametric optimal blood flow control: application to coronary bypass graft

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. The control variable is the unknown value of the normal stress at the outlet boundary, which is need for a correct set-up of the outlet boundary condition. For the numerical solution of the parametric optimal flow control problem, we develop a data-driven reduced order method that combines proper orthogonal decomposition (POD) with neural networks. We present numerical results showing that our data-driven approach leads to a substantial speed-up with respect to a more classical POD-Galerkin strategy proposed in [59], while having comparable accuracy.

preprint2022arXiv

A novel Large Eddy Simulation model for the Quasi-Geostrophic Equations in a Finite Volume setting

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 on the classical double-gyre wind forcing benchmark. The numerical experiments we present demonstrate that our nonlinear filter is an improvement over a linear filter since it is able to recover the four-gyre pattern of the time-averaged stream function even with extremely coarse meshes. In addition, our LES approach provides an average kinetic energy that compares well with the one computed with a Direct Numerical Simulation.

preprint2022arXiv

A POD-Galerkin reduced order model for the Navier-Stokes equations in stream function-vorticity formulation

We develop a Proper Orthogonal Decomposition (POD)-Galerkin based Reduced Order Model (ROM) for the efficient numerical simulation of the parametric Navier-Stokes equations in the stream function-vorticity formulation. Unlike previous works, we choose different reduced coefficients for the vorticity and stream function fields. In addition, for parametric studies we use a global POD basis space obtained from a database of time dependent full order snapshots related to sample points in the parameter space. We test the performance of our ROM strategy with the vortex merger benchmark. Accuracy and efficiency are assessed for both time reconstruction and physical parametrization.

preprint2022arXiv

An introduction to POD-Greedy-Galerkin reduced basis method

Partial differential equations can be used to model many problems in several fields of application including, e.g., fluid mechanics, heat and mass transfer, and electromagnetism. Accurate discretization methods (e.g., finite element or finite volume methods, the so-called full order models) are widely used to numerically solve these problems. However, when many physical and/or geometrical parameters are involved, the computational cost required by full order models becomes prohibitively expensive and this is not acceptable for real-time computations that are becoming more and more popular for rapid prototyping. Therefore, there is the need to introduce reduced order methods (also referred to as reduced basis methods) able to provide, as the input parameters change, fast and reliable solutions at a reduced computational cost.

preprint2022arXiv

Data-Driven Enhanced Model Reduction for Bifurcating Models in Computational Fluid Dynamics

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. Several numerical examples are shown to demonstrate the feasibility of the presented approaches.

preprint2022arXiv

Embedded domain Reduced Basis Models for the shallow water hyperbolic equations with the Shifted Boundary Method

We consider fully discrete embedded finite element approximations for a shallow water hyperbolic problem and its reduced-order model. Our approach is based on a fixed background mesh and an embedded reduced basis. The Shifted Boundary Method for spatial discretization is combined with an explicit predictor/multi-corrector time integration to integrate in time the numerical solutions to the shallow water equations, both for the full and reduced-order model. In order to improve the approximation of the solution manifold also for geometries that are untested during the offline stage, the snapshots have been pre-processed by means of an interpolation procedure that precedes the reduced basis computation. The methodology is tested on geometrically parametrized shapes with varying size and position.

preprint2022arXiv

Fast and accurate numerical simulations for the study of coronary artery bypass grafts by artificial neural network

In this work a machine learning-based Reduced Order Model (ROM) is developed to investigate in a rapid and reliable way the hemodynamic patterns in a patient-specific configuration of Coronary Artery Bypass Graft (CABG). The computational domain is composed by the left branches of coronary arteries when a stenosis of the Left Main Coronary Artery (LMCA) occurs. A reduced basis space is extracted from a collection of Finite Volume (FV) solutions of the incompressible Navier-Stokes equations by using the Proper Orthogonal Decomposition (POD) algorithm. Artificial Neural Networks (ANNs) are employed to compute the modal coefficients. Stenosis is introduced by morphing the volume meshes with a Free Form Deformation (FFD) by means of a Non-Uniform Rational Basis Spline (NURBS) volumetric parameterization.

preprint2022arXiv

Model Reduction Using Sparse Polynomial Interpolation for the Incompressible Navier-Stokes Equations

This work investigates the use of sparse polynomial interpolation as a model order reduction method for the 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 access the accuracy of the reduced order models (ROMs), in particular parametric nonlinearities arising from curved geometries are investigated in detail. Besides the accuracy of the ROMs, other important features of the method are covered, such as offline-online splitting, run time and ease of implementation. The findings establish sparse polynomial interpolation as another instrument in the toolbox of methods for breaking the curse of dimensionality.

preprint2022arXiv

Neural-network learning of SPOD latent dynamics

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). The proposed method is based on three fundamental steps: in the first, once that the mean flow field has been subtracted from the realizations (also referred to as snapshots), we compress the data from a high-dimensional representation to a lower dimensional one by constructing the SPOD latent space; in the second, we build the time-dependent coefficients by projecting the snapshots containing the fluctuations onto the SPOD basis and we learn their evolution in time with the aid of recurrent neural networks; in the third, we reconstruct the high-dimensional data from the learnt lower-dimensional representation. The proposed method is demonstrated on two different test cases, namely, a compressible jet flow, and a geophysical problem known as the Madden-Julian Oscillation. An extensive comparison between SPOD and the equivalent POD-based counterpart is provided and differences between the two approaches are highlighted. The numerical results suggest that the proposed model is able to provide low rank predictions of complex statistically stationary data and to provide insights into the evolution of phenomena characterized by specific range of frequencies. The comparison between POD and SPOD surrogate strategies highlights the need for further work on the characterization of the interplay of error between data reduction techniques and neural network forecasts.

preprint2022arXiv

Non-linear manifold ROM with Convolutional Autoencoders and Reduced Over-Collocation method

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

preprint2022arXiv

Reduced order modeling for spectral element methods: current developments in Nektar++ and further perspectives

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. An affine parameter dependency is exploited such that the reduced order model can be evaluated independent of the large-scale discretization size. The implementation in the context of SEM can be found in the open-source model reduction software ITHACA-SEM.

preprint2022arXiv

Scientific Machine Learning through Physics-Informed Neural Networks: Where we are and What's next

Physics-Informed Neural Networks (PINN) are neural networks (NNs) that encode model equations, like Partial Differential Equations (PDE), as a component of the neural network itself. PINNs are nowadays used to solve PDEs, fractional equations, integral-differential equations, and stochastic PDEs. This novel methodology has arisen as a multi-task learning framework in which a NN must fit observed data while reducing a PDE residual. This article provides a comprehensive review of the literature on PINNs: while the primary goal of the study was to characterize these networks and their related advantages and disadvantages. The review also attempts to incorporate publications on a broader range of collocation-based physics informed neural networks, which stars form the vanilla PINN, as well as many other variants, such as physics-constrained neural networks (PCNN), variational hp-VPINN, and conservative PINN (CPINN). The study indicates that most research has focused on customizing the PINN through different activation functions, gradient optimization techniques, neural network structures, and loss function structures. Despite the wide range of applications for which PINNs have been used, by demonstrating their ability to be more feasible in some contexts than classical numerical techniques like Finite Element Method (FEM), advancements are still possible, most notably theoretical issues that remain unresolved.

preprint2022arXiv

Thermomechanical modelling for industrial applications

In this work we briefly present a thermomechanical model that could serve as starting point for industrial applications. We address the non-linearity due to temperature dependence of material properties and heterogeneity due to presence of different materials. Finally a numerical example related to the simplified geometry of blast furnace hearth walls is shown with the aim of assessing the feasibility of the modelling framework. Keywords: nonlinear thermomechanical model, finite element method, heterogeneous material, blast furnace.

preprint2021arXiv

A Certified Reduced Basis Method for Linear Parametrized Parabolic Optimal Control Problems in Space-Time Formulation

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. Then, we derive error estimators to apply a Greedy method during the offline stage, in order to perform, during the online stage, a Galerkin projection onto a low-dimensional space spanned by properly chosen high-fidelity solutions. We tested the error estimators on two model problems governed by a Graetz flow: a physical parametrized distributed optimal control problem and a boundary optimal control problem with physical and geometrical parameters. The results have been compared to a previously proposed bound, based on the exact computation of the Babuška inf-sup constant, in terms of reliability and computational costs. We remark that our findings still hold in the steady setting and we propose a brief insight also for this simpler formulation.

preprint2021arXiv

A Comparison of Reduced-Order Modeling Approaches Using Artificial Neural Networks for PDEs with Bifurcating Solutions

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. Several criteria are compared and it is found that a criterion based on a regression artificial neural network (ANN) provides the most accurate results for a channel flow problem exhibiting a supercritical pitchfork bifurcation. The same benchmark test is then used to compare the local ROM approach with the regression ANN selection criterion to an established global projection-based ROM and a recently proposed ANN based method called POD-NN. We show that our local ROM approach gains more than an order of magnitude in accuracy over the global projection-based ROM. However, the POD-NN provides consistently more accurate approximations than the local projection-based ROM.

preprint2021arXiv

A numerical approach for heat flux estimation in thin slabs continuous casting molds using data assimilation

In the present work, we consider the industrial problem of estimating in real-time the mold-steel heat flux in continuous casting mold. We approach this problem by first considering the mold modeling problem (direct problem). Then, we plant the heat flux estimation problem as the inverse problem of estimating a Neumann boundary condition having as data pointwise temperature measurements in the interior of the mold domain. We also consider the case of having a total heat flux measurement together with the temperature measurements. We develop two methodologies for solving this inverse problem. The first one is the traditional Alifanov's regularization, the second one exploits the parameterization of the heat flux. We develop the latter method to have an offline-online decomposition with a computationally efficient online part to be performed in real-time. In the last part of this work, we test these methods on academic and industrial benchmarks. The results show that the parameterization method outclasses Alifanov's regularization both in performance and computational cost. Moreover, it proves to be robust with respect to the measurements noise. Finally, the tests confirm that the computational cost is suitable for real-time estimation of the heat flux.

preprint2021arXiv

ATHENA: Advanced Techniques for High Dimensional Parameter Spaces to Enhance Numerical Analysis

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. Source code, documentation, and several tutorials are available on GitHub at https://github.com/mathLab/ATHENA under the MIT license.

preprint2021arXiv

Hull shape design optimization with parameter space and model reductions, and self-learning mesh morphing

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 parameterization of the shape is applied directly to the computational mesh, propagating the generic deformation map applied to the surface (of the object to optimize) to the mesh nodes using a radial basis function (RBF) interpolation. Thus, topology and quality of the original mesh are preserved, enabling application of POD-based reduced order modeling techniques, and avoiding the necessity of additional meshing steps. Model order reduction is performed coupling POD and Gaussian process regression (GPR) in a data-driven fashion. The framework is validated on a benchmark ship.

preprint2021arXiv

Model order reduction for bifurcating phenomena in Fluid-Structure Interaction problems

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 behaviour. We have addressed the computational burden by developing a reduced order branch-wise algorithm based on a monolithic Proper Orthogonal Decomposition. We compared different constitutive relations for the solid, and we observed that a nonlinear hyper-elastic law delays the bifurcation w.r.t. the standard model, while the same effect is even magnified when considering linear elastic solid.

preprint2021arXiv

The Neural Network shifted-Proper Orthogonal Decomposition: a Machine Learning Approach for Non-linear Reduction of Hyperbolic Equations

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

preprint2020arXiv

A POD-Galerkin reduced order model of a turbulent convective buoyant flow of sodium over a backward-facing step

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. The contribution of the eddy viscosity and turbulent thermal diffusivity fields are considered in the reduced order model with an interpolation based data-driven method. The reduced order model is tested for buoyancy-aided turbulent liquid sodium flow over a vertical backward-facing step with a uniform heat flux applied on the wall downstream of the step. The wall heat flux is incorporated with a Neumann boundary condition in both the full order model and the reduced order model. The velocity and temperature profiles predicted with the reduced order model for the same and new Richardson numbers inside the range of parameter values are in good agreement with the RANS simulations. Also, the local Stanton number and skin friction distribution at the heated wall are qualitatively well captured. Finally, the reduced order simulations, performed on a single core, are about $10^5$ times faster than the RANS simulations that are performed on eight cores.

preprint2020arXiv

A Reduced Order Modeling technique to study bifurcating phenomena: application to the Gross-Pitaevskii 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. In order to reduce these demanding computational costs, our approach combines a continuation technique and Newton's method with a Reduced Order Modeling (ROM) technique, suitably supplemented with a hyper-reduction method. To demonstrate the effectiveness of our ROM approach, we trace the steady solution branches of a nonlinear Schrödinger equation, called Gross-Pitaevskii equation, as one or two physical parameters are varied. In the two parameter study, we show that our approach is 60 times faster in constructing a bifurcation diagram than a standard Full Order Method.

preprint2020arXiv

A Reduced-Order Shifted Boundary Method for Parametrized incompressible Navier-Stokes equations

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. This approach avoids both remeshing and the development of a reference domain formulation, as typically done in fitted mesh finite element formulations. Two-dimensional computational examples for one and three-parameter dimensions are presented to validate the convergence and the efficacy of the proposed approach.

preprint2020arXiv

Efficient computation of bifurcation diagrams with a deflated approach to reduced basis spectral element method

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. In this work we implemented an elaborated deflated continuation method, that relies on the spectral element method (SEM) and on the reduced basis (RB) one, to efficiently compute bifurcation diagrams with more parameters and more bifurcation points. The deflated continuation method can be obtained combining the classical continuation method and the deflation one: the former is used to entirely track each known branch of the diagram, while the latter is exploited to discover the new ones. Finally, when more than one parameter is considered, the efficiency of the computation is ensured by the fact that the diagrams can be computed during the online phase while, during the offline one, one only has to compute one-dimensional diagrams. In this work, after a more detailed description of the method, we will show the results that can be obtained using it to compute a bifurcation diagram associated with a problem governed by the Navier-Stokes equations.

preprint2020arXiv

Efficient Geometrical parametrization for Finite-Volume based Reduced Order Methods

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. In the first numerical example, different mesh motion strategies, based on a Laplacian smoothing technique and on a Radial Basis Function approach, are analyzed and compared on a heat transfer problem. Particular attention is devoted to the role of the non-orthogonal correction. In the second numerical example the methodology is tested on a geometrically parametrized incompressible Navier--Stokes problem. In this case, the reduced order model is constructed following the same segregated approach used at the full order level

preprint2020arXiv

MicroROM: An Efficient and Accurate Reduced Order Method to Solve Many-Query Problems in Micro-Motility

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. Reduced order modeling solves this issue by taking advantage of a proper subdivision between a computationally expensive offline phase and a fast and efficient online stage. This work presents the coupling of Boundary Element Method (BEM) and Reduced Basis (RB) Reduced Order Modeling (ROM) in two models of practical interest, obtaining accurate and reliable solutions to different many-query problems. Comparisons of standard reduced order modeling approaches in different simulation settings and a comparison to typical approximations to Stokes equations are also shown. Different couplings between a solver based on a HPC boundary element method for micro-motility problems and reduced order models are presented in detail. The methodology is tested on two different models: a robotic-bacterium-like and an Eukaryotic-like swimmer, and in each case two resolution strategies for the swimming problem, the split and monolithic one, are used as starting points for the ROM. An efficient and accurate reconstruction of the performance of interest is achieved in both cases proving the effectiveness of our strategy.

preprint2019arXiv

Reduced Basis Model Order Reduction for Navier-Stokes equations in domains with walls of varying curvature

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. The proposed reduced order model produces accurate approximations of steady-state solutions for a wide range of geometries and kinematic viscosity values. The application that motivated the present study is the onset of asymmetries (i.e., symmetry breaking bifurcation) in blood flow through a regurgitant mitral valve, depending on the Reynolds number and the valve shape. Through our computational study, we found that the critical Reynolds number for the symmetry breaking increases as the wall curvature increases.

preprint2018arXiv

An integrated data-driven computational pipeline with model order reduction for industrial and applied mathematics

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. Moreover several industrial examples are illustrated.