Application of a Macroscopic Turbulence Model to Simulation of Flow in a Channel With Equally Spaced Porous Fins

Author(s):  
Marcelo J. S. de Lemos ◽  
Luzia A. Tofaneli

In this work, numerical solutions are presented for turbulent flow in a channel containing fins made with porous material. The condition of spatially periodic cell is applied longitudinally along the channel. A macroscopic two-equation turbulence model is employed in both the porous region and the clear fluid. The equations of mass continuity, momentum and turbulence transport equations are written for an elementary representative volume yielding a set of equations valid for the entire computational domain. Results are presented for the velocity field as a function of Reynolds, porosity and permeability of the fins. Pressure drop along the channel is compared with the case of solid material.

Author(s):  
Marcelo J. S. de Lemos ◽  
Luzia A. Tofaneli

In this work, numerical solutions are presented for turbulent flow in a channel containing fins made with porous material. The condition of spatially periodic cell is applied longitudinally along the channel. A macroscopic tow-equation turbulence model is employed in both the porous region and the clear fluid. The equations of momentum, mass continuity and turbulence transport equations are written for an elementary representative volume yielding a set of equations valid for the entire computational domain. These equations are discretized using the control volume method and the resulting systems of algebraic equations is relaxed with the SIMPLE method. Results are presented for the velocity field as a function of Reynolds number, porosity and permeability of the fins.


Author(s):  
Marcelo J. S. de Lemos ◽  
Renato A. Silva

A number of natural and engineering systems can be characterized by some sort of porous structure through which a working fluid permeates. Atmospheric boundary layers over tropical forests and vegetation can be modeled as flow over a porous layer of irregular surface. In addition, in engineering systems one can have components that make use of a working fluid flowing over irregular layers of porous material. This paper presents numerical solutions for such hybrid medium, considering here a channel partially filled with a sinusoidal porous layer saturated by a fluid flowing in laminar regime. One unique set of transport equations is applied to both regions. Effects of Reynolds number, porosity and permeability on mean and turbulence fields are investigated. For a fixed inlet mass flow rate, increase of either porosity or permeability reduced the strength of the recirculating motion over the porous layer.


Author(s):  
Jannette B. Frandsen ◽  
Alistair G. L. Borthwick

Nonlinear effects of standing wave motions in fixed and vertically excited tanks are numerically investigated. The present fully nonlinear model analyses two-dimensional waves in stable and unstable regions of the free-surface flow. Numerical solutions of the governing nonlinear potential flow equations are obtained using a finite-difference time-stepping scheme on adaptively mapped grids. A σ-transformation in the vertical direction that stretches directly between the free-surface and bed boundary is applied to map the moving free surface physical domain onto a fixed computational domain. A horizontal linear mapping is also applied, so that the resulting computational domain is rectangular, and consists of unit square cells. The small-amplitude free-surface predictions in the fixed and vertically excited tanks compare well with 2nd order small perturbation theory. For stable steep waves in the vertically excited tank, the free-surface exhibits nonlinear behaviour. Parametric resonance is evident in the instability zones, as the amplitudes grow exponentially, even for small forcing amplitudes. For steep initial amplitudes the predictions differ considerably from the small perturbation theory solution, demonstrating the importance of nonlinear effects. The present numerical model provides a simple way of simulating steep non-breaking waves. It is computationally quick and accurate, and there is no need for free surface smoothing because of the σ-transformation.


2021 ◽  
Author(s):  
Daniel Ferreira Corrêa Barbosa ◽  
Daniel da Silva Tonon ◽  
Luiz Henrique Lindquist Whitacker ◽  
Jesuino Takachi Tomita ◽  
Cleverson Bringhenti

Abstract The aim of this work is an evaluation of different turbulence models applied in Computational Fluid Dynamics (CFD) techniques in the turbomachinery area, in this case, in an axial turbine stage used in turbopump (TP) application. The tip clearance region was considered in this study because it has a high influence in turbomachinery performance. In this region, due to its geometry and the relative movement between the rotor row and casing, there are losses associated with vortices and secondary flow making the flowfield even more turbulent and complex. Moreover, the flow that leaks in the tip region does not participate in the energy transfer between the fluid and rotor blades, degradating the machine efficiency and performance. In this work, the usual flat tip rotor blade geometry was considered. The modeling of turbulent flow based on Reynolds Averaged Navier-Stokes (RANS) equations predicts the variation of turbine operational characteristics that is sufficient for the present turbomachine and flow analysis. Therefore, the appropriate choice of the turbulence model for the study of a given flow is essential to obtain adequate results using numerical approximations. This comparison become important due to the fact that there is no general turbulence model for all engineering applications that has fluid and flow. The turbomachine considered in the present work, is the first stage of the hydraulic axial turbine used in the Low Pressure Oxidizer Turbopump (LPOTP) of the Space Shuttle Main Engine (SSME), considering the 3.0% tip clearance configuration relative to rotor blade height. The turbulence models evaluated in this work were the SST (Shear Stress Transport), the k-ε Standard and the k-ε RNG. The computational domain was discretized in several control volumes based on unstructured mesh. All the simulations were performed using the commercial software developed by ANSYS, CFX v15.0 (ANSYS). All numerical settings and how the boundary conditions were imposed at different surfaces are explained in the work. The boundary conditions settings follow the same rule used in the test facility and needs some attention during the simulations to vary the Blade-Jet-Speed ratio parameter adequately. The results from numerical simulations, were synthesized and compared with the experimental data published by National Aeronautics and Space Administration (NASA), in which the turbine efficiency and its jet velocity parameter are analyzed for each turbulence model result. The work fluid considered in this work was water, the same fluid used in the NASA test facility.


2010 ◽  
Vol 132 (5) ◽  
Author(s):  
Seok-Ki Choi ◽  
Seong-O Kim ◽  
Hoon-Ki Choi

A numerical study for the evaluation of heat transfer correlations for sodium flows in a heat exchanger of a fast breeder nuclear reactor is performed. Three different types of flows such as parallel flow, cross flow, and two inclined flows are considered. Calculations are performed for these three typical flows in a heat exchanger changing turbulence models. The tested turbulence models are the shear stress transport (SST) model and the SSG-Reynolds stress turbulence model by Speziale, Sarkar, and Gaski (1991, “Modelling the Pressure-Strain Correlation of Turbulence: An Invariant Dynamical System Approach,” J. Fluid Mech., 227, pp. 245–272). The computational model for parallel flow is a flow past tubes inside a circular cylinder and those for the cross flow and inclined flows are flows past the perpendicular and inclined tube banks enclosed by a rectangular duct. The computational results show that the SST model produces the most reliable results that can distinguish the best heat transfer correlation from other correlations for the three different flows. It was also shown that the SSG-RSTM high-Reynolds number turbulence model does not deal with the low-Prandtl number effect properly when the Peclet number is small. According to the present calculations for a parallel flow, all the old correlations do not match with the present numerical solutions and a new correlation is proposed. The correlations by Dwyer (1966, “Recent Developments in Liquid-Metal Heat Transfer,” At. Energy Rev., 4, pp. 3–92) for a cross flow and its modified correlation that takes into account of flow inclination for inclined flows work best and are accurate enough to be used for the design of the heat exchanger.


2015 ◽  
Author(s):  
A. Idris ◽  
B. P. Huynh ◽  
Z. Abdullah

Ventilation is a process of changing air in an enclosed space. Air should continuously be withdrawn and replaced by fresh air from a clean external source to maintain internal good air quality, which may referred to air quality within and around the building structures. In natural ventilation the air flow is due through cracks in the building envelope or purposely installed openings. Its can save significant amount of fossil fuel based energy by reducing the needs for mechanical ventilation and air conditioning. Numerical predictions of air velocities and the flow patterns inside the building are determined. To achieve optimum efficiency of natural ventilation, the building design should start from the climatic conditions and orography of the construction to ensure the building permeability to the outside airflow to absorb heat from indoors to reduce temperatures. Effective ventilation in a building will affects the occupant health and productivity. In this work, computational simulation is performed on a real-sized box-room with dimensions 5 m × 5 m × 5 m. Single-sided ventilation is considered whereby openings are located only on the same wall. Two opening of the total area 4 m2 are differently arranged, resulting in 16 configurations to be investigated. A logarithmic wind profile upwind of the building is employed. A commercial Computational Fluid Dynamics (CFD) software package CFD-ACE of ESI group is used. A Reynolds Average Navier Stokes (RANS) turbulence model & LES turbulence model are used to predict the air’s flow rate and air flow pattern. The governing equations for large eddy motion were obtained by filtering the Navier-Stokes and continuity equations. The computational domain was constructed had a height of 4H, width of 9H and length of 13H (H=5m), sufficiently large to avoid disturbance of air flow around the building. From the overall results, the lowest and the highest ventilation rates were obtained with windward opening and leeward opening respectively. The location and arrangement of opening affects ventilation and air flow pattern.


1989 ◽  
Vol 111 (3) ◽  
pp. 333-340 ◽  
Author(s):  
J. F. Louis ◽  
A. Salhi

The turbulent flow between two rotating co-axial disks is driven by frictional forces. The prediction of the velocity field can be expected to be very sensitive to the turbulence model used to describe the viscosity close to the walls. Numerical solutions of the Navier–Stokes equations, using a k–ε turbulence model derived from Lam and Bremhorst, are presented and compared with experimental results obtained in two different configurations: a rotating cavity and the outflow between a rotating and stationary disk. The comparison shows good overall agreement with the experimental data and substantial improvements over the results of other analyses using the k–ε models. Based on this validation, the model is applied to the flow between counterrotating disks and it gives the dependence of the radial variation of the tangential wall shear stress on Rossby number.


Author(s):  
Dongmei Zhou ◽  
Kenneth S. Ball

This paper has two objectives, (1) to examine the effects of spatial resolution, (2) to examine the effects of computational box size, upon turbulence statistics and the amount of drag reduction with and without the control scheme of wall oscillation. Direct numerical simulation (DNS) of the fully developed turbulent channel flow was performed at Reynolds number of 200 based on the wall-shear velocity and the channel half-width by using spectral methods. For the first objective, four different grids were applied to the same computational domain and the biggest impact was observed on the logarithmic law of mean velocity profiles and on the amount of drag reduction with 28.3% for the coarsest mesh and 35.4% for the finest mesh. Other turbulence features such as RMS velocity fluctuations, RMS vorticity fluctuations, and bursting events were either overpredicted or underpredicted through coarse grids. For the second objective, two different minimal channels and one natural full channel were studied and 3% drag reduction difference was observed between the smallest minimal channel of 39.1% and the natural full channel of 36.2%. In the near-wall region, however, the minimal channel flow did not exhibit significant difference in the mean velocity profiles and other lower-order statistics. Finally, from this systematical study, it showed that the accuracy of DNS depends more on the spanwise resolution, and it also confirmed that a minimal channel model is able to catch key structures of turbulence in the near-wall region but is much less expensive.


Author(s):  
Alexandre Lamoureux ◽  
B. Rabi Baliga

A computational investigation of temporally- and spatially-periodic laminar two-dimensional fluid flow and heat transfer in staggered-plate arrays is presented in this paper. The objective and the novel aspect of this study is the investigation of the influence (on the numerical solutions) of including single and multiple representative geometric modules in the calculation domain, with spatially-periodic boundary conditions imposed on the instantaneous velocity and temperature fields in both the streamwise and the lateral directions. The following geometrical parameters, normalized with respect to a representative module height, were studied: a dimensionless plate length equal to 1, and a dimensionless plate thickness of 0.250. This relatively high value of dimensionless plate thickness, compared to those commonly encountered in rectangular offset-fin cores of compact heat exchangers, was deliberately chosen to induce and enhance the unsteady features of the fluid flow and heat transfer phenomena. Different specified values of the time-mean modular streamwise gradient of the reduced pressure were investigated, yielding values of Reynolds number (Kays and London definition) in the range of 100 to 625. The Prandtl number was fixed at 0.7. In the multiple-module simulations, for Reynolds number values exceeding 400, it was found that multiple solutions are possible: the particular solution which is obtained in any one simulation depends on the specified initial conditions. The results presented include time-mean modular friction factors, modular Colburn factors, and Strouhal numbers.


2015 ◽  
Vol 137 (9) ◽  
Author(s):  
Ya-Ling He ◽  
Wen-Quan Tao

In this paper, numerical simulation approaches for multiscale process of heat transfer and fluid flow are briefly reviewed, and the existing coupling algorithms are summarized. These molecular dynamics simulation (MDS)–finite volume method (FVM), MD–lattice Boltzmann method (LBM), and direct simulation of Monte Carlo method (DSMC)–FVM. The available reconstruction operators for LBM–FVM coupling are introduced. Four multiscale examples for fluid flow and heat transfer are presented by using these coupled methods. It is shown that by coupled method different resolution requirements in the computational domain can be satisfied successfully while computational time can be significantly saved. Further research needs for the study of multiscale heat transfer and fluid flow problems are proposed.


Sign in / Sign up

Export Citation Format

Share Document