2,104 research outputs found
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 times faster than the RANS simulations that
are performed on eight cores.Comment: 26 pages, 15 figures, 3 table
Modeling Thermal Turbulence Using Implicit Large Eddy Simulation
A general description of a thermally coupled fluid flow is given by the incompressible Navier-Stokes equations coupled with the heat equation using Boussinesq approximation, whose mathematical structure is much well understood. A variational multiscale finite element approximation has been considered for the formulation of incompressible Navier-Stokes equation and heat equation. The complexity of these problems makes their numerical solution very difficult as the standard finite element method is unstable. In the incompressible Navier Stokes equations, two well known sources of numerical instabilities are the incompressibility constraint and the presence of the convective term.
Many stabilization techniques used nowadays are based on scale separation, splitting the unknown into a coarse part induced by the discretization of the domain and a fine subgrid part. The modeling of the subgrid scale and its influence leads to a modified coarse scale problem providing stability.
In convection-diffusion problem once global instabilities have been overcome by a stabilization method, there are still local oscillations near layers due to the lack of monotonicity of the method.
Shock capturing techniques are often employed to deal with them. Proper choice of stabilization and shock capturing techniques can eliminate the local instabilities near layers of convection-diffusion equation.
A very important issue of the formulation presented in this thermally coupled incompressible flow is the possibility to model turbulent flows. Some terms involving the velocity subgrid scale arise from the convective term in the Navier-Stokes equations which can be understood as the contribution from the Reynolds tensor of a LES approach and the contribution from the cross stress tensor. This opens the door of modeling thermal turbulence using LES automatically inherited by the formulation
used in this work.
Different classical benchmark problems are numerically solved in this thesis work for the convection-diffusion equation to show the capabilities of different combination of stabilization and shock capturing methods. In the case of thermally coupled incompressible flows some numerical and industrial examples are exhibited to check the performance of the different combination of stabilization
and shock capturing methods and to compare them. The objective is to conclude which method works better to approximate the exact solution and eliminate instabilities and local oscillations
- …