-
A time-adaptive algorithm for pressure dominated flows: a heuristic estimator
Authors:
Ivan Prusak,
Davide Torlo,
Monica Nonino,
Gianluigi Rozza
Abstract:
This work aims to introduce a heuristic timestep-adaptive algorithm for Computational Fluid Dynamics (CFD) and Fluid-Structure Interaction (FSI) problems where the flow is dominated by the pressure. In such scenarios, many time-adaptive algorithms based on the interplay of implicit and explicit time schemes fail to capture the fast transient dynamics of pressure fields. We present an algorithm tha…
▽ More
This work aims to introduce a heuristic timestep-adaptive algorithm for Computational Fluid Dynamics (CFD) and Fluid-Structure Interaction (FSI) problems where the flow is dominated by the pressure. In such scenarios, many time-adaptive algorithms based on the interplay of implicit and explicit time schemes fail to capture the fast transient dynamics of pressure fields. We present an algorithm that relies on a temporal error estimator using Backward Differentiation Formulae (BDF$k$) of order $k=2,3$. Specifically, we demonstrate that the implicit BDF$3$ solution can be well approximated by applying a single Newton-type nonlinear solver correction to the implicit BDF$2$ solution. The difference between these solutions determines our adaptive temporal error estimator. The effectiveness of our approach is confirmed by numerical experiments conducted on a backward-facing step flow CFD test case with Reynolds number $300$ and on a two-dimensional haemodynamics FSI benchmark.
△ Less
Submitted 29 June, 2024;
originally announced July 2024.
-
Calibration-Based ALE Model Order Reduction for Hyperbolic Problems with Self-Similar Travelling Discontinuities
Authors:
Monica Nonino,
Davide Torlo
Abstract:
We propose a novel Model Order Reduction framework that is able to handle solutions of hyperbolic problems characterized by multiple travelling discontinuities. By means of an optimization based approach, we introduce suitable calibration maps that allow us to transform the original solution manifold into a lower dimensional one. The novelty of the methodology is represented by the fact that the o…
▽ More
We propose a novel Model Order Reduction framework that is able to handle solutions of hyperbolic problems characterized by multiple travelling discontinuities. By means of an optimization based approach, we introduce suitable calibration maps that allow us to transform the original solution manifold into a lower dimensional one. The novelty of the methodology is represented by the fact that the optimization process does not require the knowledge of the discontinuities location. The optimization can be carried out simply by choosing some reference control points, thus avoiding the use of some implicit shock tracking techniques, which would translate into an increased computational effort during the offline phase. In the online phase, we rely on a non-intrusive approach, where the coefficients of the projection of the reduced order solution onto the reduced space are recovered by means of an Artificial Neural Network. To validate the methodology, we present numerical results for the 1D Sod shock tube problem, for the 2D double Mach reflection problem, also in the parametric case, and for the triple point problem.
△ Less
Submitted 20 August, 2024; v1 submitted 18 March, 2024;
originally announced March 2024.
-
Optimisation-Based Coupling of Finite Element Model and Reduced Order Model for Computational Fluid Dynamics
Authors:
Ivan Prusak,
Davide Torlo,
Monica Nonino,
Gianluigi Rozza
Abstract:
Using Domain Decomposition (DD) algorithm on non--overlapping domains, we compare couplings of different discretisation models, such as Finite Element (FEM) and Reduced Order (ROM) models for separate subcomponents. In particular, we consider an optimisation-based DD model where the coupling on the interface is performed using a control variable representing the normal flux. We use iterative gradi…
▽ More
Using Domain Decomposition (DD) algorithm on non--overlapping domains, we compare couplings of different discretisation models, such as Finite Element (FEM) and Reduced Order (ROM) models for separate subcomponents. In particular, we consider an optimisation-based DD model where the coupling on the interface is performed using a control variable representing the normal flux. We use iterative gradient-based optimisation algorithms to decouple the subdomain state solutions as well as to locally generate ROMs on each subdomain. Then, we consider FEM or ROM discretisation models for each of the DD problem components, namely, the triplet state1-state2-control. On the backward-facing step Navier-Stokes (NS) problem, we investigate the efficacy of the presented couplings in terms of optimisation iterations, optimal functional values and relative errors.
△ Less
Submitted 10 July, 2024; v1 submitted 16 February, 2024;
originally announced February 2024.
-
An optimisation-based domain-decomposition reduced order model for parameter-dependent non-stationary fluid dynamics problems
Authors:
Ivan Prusak,
Davide Torlo,
Monica Nonino,
Gianluigi Rozza
Abstract:
In this work, we address parametric non-stationary fluid dynamics problems within a model order reduction setting based on domain decomposition. Starting from the optimisation-based domain decomposition approach, we derive an optimal control problem, for which we present a convergence analysis in the case of non-stationary incompressible Navier-Stokes equations. We discretize the problem with the…
▽ More
In this work, we address parametric non-stationary fluid dynamics problems within a model order reduction setting based on domain decomposition. Starting from the optimisation-based domain decomposition approach, we derive an optimal control problem, for which we present a convergence analysis in the case of non-stationary incompressible Navier-Stokes equations. We discretize the problem with the finite element method and we compare different model order reduction techniques: POD-Galerkin and a non-intrusive neural network procedure. We show that the classical POD-Galerkin is more robust and accurate also in transient areas, while the neural network can obtain simulations very quickly though being less precise in the presence of discontinuities in time or parameter domain. We test the proposed methodologies on two fluid dynamics benchmarks with physical parameters and time dependency: the non-stationary backward-facing step and lid-driven cavity flow.
△ Less
Submitted 19 February, 2024; v1 submitted 3 August, 2023;
originally announced August 2023.
-
Geometry-based approximation of waves in complex domains
Authors:
Davide Pradovera,
Monica Nonino,
Ilaria Perugia
Abstract:
We consider wave propagation problems over 2-dimensional domains with piecewise-linear boundaries, possibly including scatterers. We assume that the wave speed is constant, and that the initial conditions and forcing terms are radially symmetric and compactly supported. We propose an approximation of the propagating wave as the sum of some special space-time functions. Each term in this sum identi…
▽ More
We consider wave propagation problems over 2-dimensional domains with piecewise-linear boundaries, possibly including scatterers. We assume that the wave speed is constant, and that the initial conditions and forcing terms are radially symmetric and compactly supported. We propose an approximation of the propagating wave as the sum of some special space-time functions. Each term in this sum identifies a particular field component, modeling the result of a single reflection or diffraction effect. We describe an algorithm for identifying such components automatically, based on the domain geometry. To showcase our proposed method, we present several numerical examples, such as waves scattering off wedges and waves propagating through a room in presence of obstacles. Software implementing our numerical algorithm is made available as open-source code.
△ Less
Submitted 26 June, 2024; v1 submitted 31 January, 2023;
originally announced January 2023.
-
An optimisation-based domain-decomposition reduced order model for the incompressible Navier-Stokes equations
Authors:
Ivan Prusak,
Monica Nonino,
Davide Torlo,
Francesco Ballarin,
Gianluigi Rozza
Abstract:
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 fo…
▽ More
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. Firstly, the problem is described on the subdomains coupled at the interface and solved through an optimal control problem, which leads to the complete separation of the subdomain problems in the DD method. On top of that, a reduced model for the obtained optimal-control problem is built; the procedure is based on the Proper Orthogonal Decomposition technique and a further Galerkin projection. The presented methodology is tested on two fluid dynamics benchmarks: the stationary backward-facing step and lid-driven cavity flow. The numerical tests show a significant reduction of the computational costs in terms of both the problem dimensions and the number of optimisation iterations in the domain-decomposition algorithm.
△ Less
Submitted 3 August, 2023; v1 submitted 26 November, 2022;
originally announced November 2022.
-
Projection based semi--implicit partitioned Reduced Basis Method for non parametrized and parametrized Fluid--Structure Interaction problems
Authors:
Monica Nonino,
Francesco Ballarin,
Gianluigi Rozza,
Yvon Maday
Abstract:
We present a partitioned Model Order Reduction method for multiphysics problems, that is based on a semi-implicit treatment of the coupling conditions, and on a projection scheme. The proposed Reduced Order Method is based on the Proper Orthogonal Decomposition and on a Galerkin projection onto the reduced basis spaces; we aim of addressing both time-dependent and time-dependent, parametrized Flui…
▽ More
We present a partitioned Model Order Reduction method for multiphysics problems, that is based on a semi-implicit treatment of the coupling conditions, and on a projection scheme. The proposed Reduced Order Method is based on the Proper Orthogonal Decomposition and on a Galerkin projection onto the reduced basis spaces; we aim of addressing both time-dependent and time-dependent, parametrized Fluid-Structure Interaction problems, where the fluid is incompressible and the structure is linear, elastic and two dimensional.
△ Less
Submitted 3 November, 2022; v1 submitted 10 January, 2022;
originally announced January 2022.
-
A monolithic and a partitioned Reduced Basis Method for Fluid-Structure Interaction problems
Authors:
Monica Nonino,
Francesco Ballarin,
Gianluigi Rozza
Abstract:
The aim of this work is to present an overview about the combination of the Reduced Basis Method (RBM) with two different approaches for Fluid-Structure Interaction (FSI) problems, namely a monolithic and a partitioned approach. We provide the details of implementation of two reduction procedures, and we then apply them to the same test case of interest. We first implement a reduction technique th…
▽ More
The aim of this work is to present an overview about the combination of the Reduced Basis Method (RBM) with two different approaches for Fluid-Structure Interaction (FSI) problems, namely a monolithic and a partitioned approach. We provide the details of implementation of two reduction procedures, and we then apply them to the same test case of interest. We first implement a reduction technique that is based on a monolithic procedure where we solve the fluid and the solid problems all at once. We then present another reduction technique that is based on a partitioned (or segregated) procedure: the fluid and the solid problems are solved separately and then coupled using a fixed point strategy. The toy problem that we consider is based on the Turek-Hron benchmark test case, with a fluid Reynolds number Re = 100.
△ Less
Submitted 21 June, 2021; v1 submitted 20 April, 2021;
originally announced April 2021.
-
A Reduced Order Cut Finite Element method for geometrically parameterized steady and unsteady Navier-Stokes problems
Authors:
Efthymios N. Karatzas,
Monica Nonino,
Francesco Ballarin,
Gianluigi Rozza
Abstract:
This work focuses on steady and unsteady Navier-Stokes equations in a reduced order modeling framework. The methodology proposed is based on a Proper Orthogonal Decomposition within a levelset geometry description and the problems of interest are discretized with an unfitted mesh Finite Element Method. We construct and investigate a unified and geometry independent reduced basis which overcomes ma…
▽ More
This work focuses on steady and unsteady Navier-Stokes equations in a reduced order modeling framework. The methodology proposed is based on a Proper Orthogonal Decomposition within a levelset geometry description and the problems of interest are discretized with an unfitted mesh Finite Element Method. 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. By employing a geometry independent reduced basis, we are able to avoid remeshing and transformation to reference configurations, and we are able to handle complex geometries. This combination of a fixed background mesh in a fixed extended background geometry with reduced order techniques appears beneficial and advantageous in many industrial and engineering applications, which could not be resolved efficiently in the past.
△ Less
Submitted 10 May, 2021; v1 submitted 10 October, 2020;
originally announced October 2020.
-
Overcoming slowly decaying Kolmogorov n-width by transport maps: application to model order reduction of fluid dynamics and fluid--structure interaction problems
Authors:
Monica Nonino,
Francesco Ballarin,
Gianluigi Rozza,
Yvon Maday
Abstract:
In this work we focus on reduced order modelling for problems for which the resulting reduced basis spaces show a slow decay of the Kolmogorov $n$-width, or, in practical calculations, its computational surrogate given by the magnitude of the eigenvalues returned by a proper orthogonal decomposition on the solution manifold. In particular, we employ an additional preprocessing during the offline p…
▽ More
In this work we focus on reduced order modelling for problems for which the resulting reduced basis spaces show a slow decay of the Kolmogorov $n$-width, or, in practical calculations, its computational surrogate given by the magnitude of the eigenvalues returned by a proper orthogonal decomposition on the solution manifold. In particular, we employ an additional preprocessing during the offline phase of the reduced basis method, in order to obtain smaller reduced basis spaces. Such preprocessing is based on the composition of the snapshots with a transport map, that is a family of smooth and invertible mappings that map the physical domain of the problem into itself. Two test cases are considered: a fluid moving in a domain with deforming walls, and a fluid past a rotating cylinder. Comparison between the results of the novel offline stage and the standard one is presented.
△ Less
Submitted 15 November, 2019;
originally announced November 2019.