TY - JOUR T1 - A novel iterative penalty method to enforce boundary conditions in Finite Volume POD-Galerkin reduced order models for fluid dynamics problems JF - Communications in Computational Physics Y1 - 2021 A1 - Kelbij Star A1 - Giovanni Stabile A1 - Francesco Belloni A1 - Gianluigi Rozza A1 - Joris Degroote AB - A Finite-Volume based POD-Galerkin reduced order model is developed for fluid dynamic problems where the (time-dependent) boundary conditions are controlled using two different boundary control strategies: the control function method, whose aim is to obtain homogeneous basis functions for the reduced basis space and the penalty method where the boundary conditions are enforced in the reduced order model using a penalty factor. The penalty method is improved by using an iterative solver for the determination of the penalty factor rather than tuning the factor with a sensitivity analysis or numerical experimentation. The boundary control methods are compared and tested for two cases: the classical lid driven cavity benchmark problem and a Y-junction flow case with two inlet channels and one outlet channel. The results show that the boundaries of the reduced order model can be controlled with the boundary control methods and the same order of accuracy is achieved for the velocity and pressure fields. Finally, the speedup ratio between the reduced order models and the full order model is of the order 1000 for the lid driven cavity case and of the order 100 for the Y-junction test case. PB - Global Science Press VL - 30 ER - TY - JOUR T1 - A POD-Galerkin reduced order model of a turbulent convective buoyant flow of sodium over a backward-facing step JF - Applied Mathematical Modelling Y1 - 2021 A1 - Kelbij Star A1 - Giovanni Stabile A1 - Gianluigi Rozza A1 - Joris Degroote AB -

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 105 times faster than the RANS simulations that are performed on eight cores.

VL - 89 ER - TY - JOUR T1 - Reduced order models for the incompressible Navier-Stokes equations on collocated grids using a `discretize-then-project' approach JF - International Journal for Numerical Methods in Fluids Y1 - 2021 A1 - Kelbij Star A1 - Benjamin Sanderse A1 - Giovanni Stabile A1 - Gianluigi Rozza A1 - Joris Degroote PB - Wiley VL - 93 UR - https://doi.org/10.1002/fld.4994 ER - TY - JOUR T1 - A hybrid reduced order method for modelling turbulent heat transfer problems JF - Computers & Fluids Y1 - 2020 A1 - Sokratia Georgaka A1 - Giovanni Stabile A1 - Kelbij Star A1 - Gianluigi Rozza A1 - Michael J. Bluck AB -

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. For high Reynolds numbers, the eddy viscosity and eddy diffusivity are incorporated into the reduced order model with a Proper Orthogonal Decomposition (nested and standard) with Interpolation (PODI), where the interpolation is performed using Radial Basis Functions. The reduced order solver, obtained using a k-ω SST URANS full order model, is tested against the full order solver in a 3D T-junction pipe with parametric velocity inlet boundary conditions.

VL - 208 UR - https://arxiv.org/abs/1906.08725 ER - TY - CONF T1 - POD-Galerkin Reduced Order Model of the Boussinesq Approximation for Buoyancy-Driven Enclosed Flows T2 - International Conference on Mathematics and Computational Methods Applied to Nuclear Science and Engineering, M and C 2019 Y1 - 2019 A1 - Kelbij Star A1 - Giovanni Stabile A1 - Sokratia Georgaka A1 - Francesco Belloni A1 - Gianluigi Rozza A1 - Joris Degroote JF - International Conference on Mathematics and Computational Methods Applied to Nuclear Science and Engineering, M and C 2019 SN - 9780894487699 ER -