%0 Journal Article %J Communications in Computational Physics %D 2021 %T A novel iterative penalty method to enforce boundary conditions in Finite Volume POD-Galerkin reduced order models for fluid dynamics problems %A Kelbij Star %A Giovanni Stabile %A Francesco Belloni %A Gianluigi Rozza %A Joris Degroote %X 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. %B Communications in Computational Physics %I Global Science Press %V 30 %P 34–66 %G eng %R https://doi.org/10.4208/cicp.OA-2020-0059 %0 Journal Article %J Applied Mathematical Modelling %D 2021 %T A POD-Galerkin reduced order model of a turbulent convective buoyant flow of sodium over a backward-facing step %A Kelbij Star %A Giovanni Stabile %A Gianluigi Rozza %A Joris Degroote %X

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.

%B Applied Mathematical Modelling %V 89 %P 486-503 %G eng %R 10.1016/j.apm.2020.07.029 %0 Journal Article %J International Journal for Numerical Methods in Fluids %D 2021 %T Reduced order models for the incompressible Navier-Stokes equations on collocated grids using a `discretize-then-project' approach %A Kelbij Star %A Benjamin Sanderse %A Giovanni Stabile %A Gianluigi Rozza %A Joris Degroote %B International Journal for Numerical Methods in Fluids %I Wiley %V 93 %P 2694–2722 %G eng %U https://doi.org/10.1002/fld.4994 %R 10.1002/fld.4994 %0 Journal Article %J Computers & Fluids %D 2020 %T A hybrid reduced order method for modelling turbulent heat transfer problems %A Sokratia Georgaka %A Giovanni Stabile %A Kelbij Star %A Gianluigi Rozza %A Michael J. Bluck %X

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.

%B Computers & Fluids %V 208 %P 104615 %G eng %U https://arxiv.org/abs/1906.08725 %R 10.1016/j.compfluid.2020.104615 %0 Conference Paper %B International Conference on Mathematics and Computational Methods Applied to Nuclear Science and Engineering, M and C 2019 %D 2019 %T POD-Galerkin Reduced Order Model of the Boussinesq Approximation for Buoyancy-Driven Enclosed Flows %A Kelbij Star %A Giovanni Stabile %A Sokratia Georgaka %A Francesco Belloni %A Gianluigi Rozza %A Joris Degroote %B International Conference on Mathematics and Computational Methods Applied to Nuclear Science and Engineering, M and C 2019 %@ 9780894487699 %G eng