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

Kelbij Star,Giovanni Stabile,Gianluigi Rozza,Joris Degroote
DOI: https://doi.org/10.48550/arXiv.2003.01114
2020-08-04
Abstract: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.
Fluid Dynamics
What problem does this paper attempt to address?
The problem that this paper attempts to solve is: to develop a Reduced Order Model (ROM) for turbulent convection and buoyancy - driven flow of fluids with low Prandtl number (Pr) in order to achieve efficient and accurate numerical simulations. Specifically, the focus of the research is on the turbulent convection and buoyancy - driven flow of sodium over a backward - facing step. ### Main Problems and Challenges 1. **Complexity of Fluids with Low Prandtl Number** - The thermal conduction boundary layer of low - Pr fluids (such as liquid metals) is thicker, resulting in different spatial scale ranges for temperature and velocity. - The heat flux becomes more important in the near - wall region, and the traditional Reynolds analogy is no longer applicable. 2. **Difficulties in Numerical Modeling** - The turbulent heat flux in low - Pr fluids is difficult to model accurately, especially in buoyancy - assisted flows. - It is necessary to deal with complex turbulence models and heat transfer mechanisms. 3. **Computational Efficiency** - The computational cost of the Full Order Model (FOM), such as the Reynolds - Averaged Navier - Stokes (RANS) in the steady state, is high, especially in cases where real - time simulations or a large number of parameter configurations are required. ### Solutions The author adopts a reduced - order - model method based on Proper Orthogonal Decomposition (POD) and Galerkin projection, specifically including: - **POD - Galerkin Method**: Extract the main characteristic modes through POD, and use Galerkin projection to project the high - dimensional full - order model onto a low - dimensional subspace, thereby significantly reducing the computational cost. - **Turbulence Model**: Use the Linear Eddy Viscosity Model to calculate the Reynolds stress, and combine the Single Gradient Diffusion Hypothesis and local correlations to evaluate the turbulent Prandtl number. - **Data - Driven Interpolation**: For the influence of the eddy viscosity and turbulent heat diffusivity fields, a data - driven interpolation - based method is adopted. ### Experimental Verification This reduced - order model is applied to the buoyancy - assisted turbulent flow of sodium over a vertical backward - facing step, with the wall condition of a uniform heat flux applied. The results show that the ROM can well predict the velocity and temperature distributions in the same and new Richardson number (Ri) ranges, and is approximately 10^5 times faster than RANS simulations. ### Summary This research aims to solve the complex numerical - modeling problems of turbulent convection and buoyancy - driven flow of low - Pr fluids by developing an efficient reduced - order model, thereby providing a faster and more economical simulation tool for industrial applications.