Lid-Driven Square Cavity Flow: A Benchmark Solution with an 8192x8192 Grid

2021 ◽  
Author(s):  
Carlos Marchi ◽  
Cosmo D. Santiago ◽  
Carlos Alberto Rezende de Carvalho Junior
Keyword(s):  
Reynolds Number ◽  
Stokes Equations ◽  
Accurate Determination ◽  
Stable Solution ◽  
Reynolds Numbers ◽  
Accurate Solution ◽  
Discretization Error ◽  
Square Cavity ◽  
Solution Vector ◽  
Wide Range

Abstract The incompressible steady-state fluid flow inside a lid-driven square cavity was simulated using the mass conservation and Navier-Stokes equations. This system of equations is solved for Reynolds numbers of up to 10,000 to the accuracy of the computational machine round-off error. The computational model used was the second-order accurate finite volume method. A stable solution is obtained using the iterative multigrid methodology with 8192 × 8192 volumes, while degree-10 interpolation and Richardson extrapolation were used to reduce the discretization error. The solution vector comprised five entries of velocities, pressure, and location. For comparison purposes, 65 different variables of interest were chosen, such as velocity profile, its extremum values and location, extremum values and location of the stream function. The discretization error for each variable of interest was estimated using two types of estimators and their apparent order of accuracy. The variations of the 11 selected variables are shown across 38 Reynolds number values between 0.0001 and 10,000. In this study, we provide a more accurate determination of the Reynolds number value at which the upper secondary vortex appears. The results of this study were compared with those of several other studies in the literature. The current solution methodology was observed to produce the most accurate solution till date for a wide range of Reynolds numbers.

scholarly journals Mixing Enhancement of Non-Newtonian Shear-Thinning Fluid for a Kenics Micromixer

Micromachines ◽  
2021 ◽  
Vol 12 (12) ◽  
pp. 1494
Author(s):  
Abdelkader Mahammedi ◽  
Naas Toufik Tayeb ◽  
Kwang-Yong Kim ◽  
Shakhawat Hossain
Keyword(s):  
Reynolds Number ◽  
Stokes Equations ◽  
Fluid Dynamic ◽  
Reynolds Numbers ◽  
Shear Thinning ◽  
Navier Stokes ◽  
Mixing Enhancement ◽  
Pressure Losses ◽  
Shear Thinning Fluids ◽  
Wide Range

In this work, a numerical investigation was analyzed to exhibit the mixing behaviors of non-Newtonian shear-thinning fluids in Kenics micromixers. The numerical analysis was performed using the computational fluid dynamic (CFD) tool to solve 3D Navier-Stokes equations with the species transport equations. The efficiency of mixing is estimated by the calculation of the mixing index for different cases of Reynolds number. The geometry of micro Kenics collected with a series of six helical elements twisted 180° and arranged alternately to achieve the higher level of chaotic mixing, inside a pipe with a Y-inlet. Under a wide range of Reynolds numbers between 0.1 to 500 and the carboxymethyl cellulose (CMC) solutions with power-law indices among 1 to 0.49, the micro-Kenics proves high mixing Performances at low and high Reynolds number. Moreover the pressure losses of the shear-thinning fluids for different Reynolds numbers was validated and represented.

The Effect of Pulsed Blowing on the Boundary Layer of a Highly Loaded Low Pressure Turbine Blade

2013 ◽  
Author(s):  
Marion Mack ◽  
Roland Brachmanski ◽  
Reinhard Niehuis
Keyword(s):  
Boundary Layer ◽  
Reynolds Number ◽  
Mach Number ◽  
Suction Side ◽  
Reynolds Numbers ◽  
Low Pressure ◽  
Trailing Edge ◽  
Low Pressure Turbine ◽  
Wide Range ◽  
Highly Loaded

The performance of the low pressure turbine (LPT) can vary appreciably, because this component operates under a wide range of Reynolds numbers. At higher Reynolds numbers, mid and aft loaded profiles have the advantage that transition of suction side boundary layer happens further downstream than at front loaded profiles, resulting in lower profile loss. At lower Reynolds numbers, aft loading of the blade can mean that if a suction side separation exists, it may remain open up to the trailing edge. This is especially the case when blade lift is increased via increased pitch to chord ratio. There is a trend in research towards exploring the effect of coupling boundary layer control with highly loaded turbine blades, in order to maximize performance over the full relevant Reynolds number range. In an earlier work, pulsed blowing with fluidic oscillators was shown to be effective in reducing the extent of the separated flow region and to significantly decrease the profile losses caused by separation over a wide range of Reynolds numbers. These experiments were carried out in the High-Speed Cascade Wind Tunnel of the German Federal Armed Forces University Munich, Germany, which allows to capture the effects of pulsed blowing at engine relevant conditions. The assumed control mechanism was the triggering of boundary layer transition by excitation of the Tollmien-Schlichting waves. The current work aims to gain further insight into the effects of pulsed blowing. It investigates the effect of a highly efficient configuration of pulsed blowing at a frequency of 9.5 kHz on the boundary layer at a Reynolds number of 70000 and exit Mach number of 0.6. The boundary layer profiles were measured at five positions between peak Mach number and the trailing edge with hot wire anemometry and pneumatic probes. Experiments were conducted with and without actuation under steady as well as periodically unsteady inflow conditions. The results show the development of the boundary layer and its interaction with incoming wakes. It is shown that pulsed blowing accelerates transition over the separation bubble and drastically reduces the boundary layer thickness.

scholarly journals Finite-amplitude steady waves in plane viscous shear flows

1985 ◽  
Vol 160 ◽  
pp. 281-295 ◽  
Author(s):  
F. A. Milinazzo ◽  
P. G. Saffman
Keyword(s):  
Boundary Layer ◽  
Reynolds Number ◽  
Stokes Equations ◽  
Reynolds Numbers ◽  
Finite Amplitude ◽  
Linear Instability ◽  
Navier Stokes ◽  
Unidirectional Flow ◽  
Viscous Shear ◽  
Finite Amplitude Waves

Computations of two-dimensional solutions of the Navier–Stokes equations are carried out for finite-amplitude waves on steady unidirectional flow. Several cases are considered. The numerical method employs pseudospectral techniques in the streamwise direction and finite differences on a stretched grid in the transverse direction, with matching to asymptotic solutions when unbounded. Earlier results for Poiseuille flow in a channel are re-obtained, except that attention is drawn to the dependence of the minimum Reynolds number on the physical constraint of constant flux or constant pressure gradient. Attempts to calculate waves in Couette flow by continuation in the velocity of a channel wall fail. The asymptotic suction boundary layer is shown to possess finite-amplitude waves at Reynolds numbers orders of magnitude less than the critical Reynolds number for linear instability. Waves in the Blasius boundary layer and unsteady Rayleigh profile are calculated by employing the artifice of adding a body force to cancel the spatial or temporal growth. The results are verified by comparison with perturbation analysis in the vicinity of the linear-instability critical Reynolds numbers.

scholarly journals Nonlinear amplification in hydrodynamic turbulence

2021 ◽  
Vol 930 ◽  
Author(s):  
Kartik P. Iyer ◽  
Katepalli R. Sreenivasan ◽  
P.K. Yeung
Keyword(s):  
Reynolds Number ◽  
Isotropic Turbulence ◽  
Stokes Equations ◽  
Three Dimensional ◽  
Reynolds Numbers ◽  
Navier Stokes ◽  
Navier Stokes Equations ◽  
Developed Turbulence ◽  
Advection Term ◽  
Nonlinear Amplification

Using direct numerical simulations performed on periodic cubes of various sizes, the largest being $8192^3$ , we examine the nonlinear advection term in the Navier–Stokes equations generating fully developed turbulence. We find significant dissipation even in flow regions where nonlinearity is locally absent. With increasing Reynolds number, the Navier–Stokes dynamics amplifies the nonlinearity in a global sense. This nonlinear amplification with increasing Reynolds number renders the vortex stretching mechanism more intermittent, with the global suppression of nonlinearity, reported previously, restricted to low Reynolds numbers. In regions where vortex stretching is absent, the angle and the ratio between the convective vorticity and solenoidal advection in three-dimensional isotropic turbulence are statistically similar to those in the two-dimensional case, despite the fundamental differences between them.

Jet-to-Plate Distance Effect on Heat Transfer From a Flat Plate to an Impinging Jet at Various Reynolds Numbers

2012 ◽  
Author(s):  
Patricia Streufert ◽  
Terry X. Yan ◽  
Mahdi G. Baygloo
Keyword(s):  
Heat Transfer ◽  
Experimental Data ◽  
Reynolds Number ◽  
Flat Plate ◽  
Numerical Solutions ◽  
Stokes Equations ◽  
Local Heat Transfer ◽  
Reynolds Numbers ◽  
Air Jet ◽  
Plate Distance

Local turbulent convective heat transfer from a flat plate to a circular impinging air jet is numerically investigated. The jet-to-plate distance (L/D) effect on local heat transfer is the main focus of this study. The eddy viscosity V2F turbulence model is used with a nonuniform structured mesh. Reynolds-Averaged Navier-Stokes equations (RANS) and the energy equation are solved for axisymmetric, three-dimensional flow. The numerical solutions obtained are compared with published experimental data. Four jet-to-plate distances, (L/D = 2, 4, 6 and 10) and seven Reynolds numbers (Re = 7,000, 15,000, 23,000, 50,000, 70,000, 100,000 and 120,000) were parametrically studied. Local and average heat transfer results are analyzed and correlated with Reynolds number and the jet-to-plate distance. Results show that the numerical solutions matched experimental data best at low jet-to-plate distances and lower Reynolds numbers, decreasing in ability to accurately predict the heat transfer as jet-to-plate distance and Reynolds number was increased.

Reconsideration of spanwise rotating turbulent channel flows via resolvent analysis

2018 ◽  
Vol 861 ◽  
pp. 200-222 ◽  
Author(s):  
Satoshi Nakashima ◽  
Mitul Luhar ◽  
Koji Fukagata
Keyword(s):  
Turbulent Flows ◽  
Stokes Equations ◽  
Turbulent Channel Flow ◽  
Reynolds Numbers ◽  
High Gain ◽  
Spectral Space ◽  
Rotation Rates ◽  
Wide Range ◽  
Effects Of Rotation ◽  
Spanwise Rotation

We study the effect of spanwise rotation in turbulent channel flow at both low and high Reynolds numbers by employing the resolvent formulation proposed by McKeon & Sharma (J. Fluid Mech., vol. 658, 2010, pp. 336–382). Under this formulation, the nonlinear terms in the Navier–Stokes equations are regarded as a forcing that acts upon the remaining linear dynamics to generate the turbulent velocity field in response. A gain-based decomposition of the forcing–response transfer function across spectral space yields models for highly amplified flow structures, or modes. Unlike linear stability analysis, this enables targeted analyses of the effects of rotation on high-gain modes that serve as useful low-order models for dynamically important coherent structures in wall-bounded turbulent flows. The present study examines a wide range of rotation rates. A posteriori comparisons at low Reynolds number ($\mathit{Re}_{\unicode[STIX]{x1D70F}}=180$) demonstrate that the resolvent formulation is able to quantitatively predict the effect of varying spanwise rotation rates on specific classes of flow structure (e.g. the near-wall cycle) as well as energy amplification across spectral space. For fixed inner-normalized rotation number, the effects of rotation at varying friction Reynolds numbers appear to be similar across spectral space, when scaled in outer units. We also consider the effects of rotation on modes with varying speed (i.e. modes that are localized in regions of varying mean shear), and provide suggestions for modelling the nonlinear forcing term.

Complex dynamics in a stratified lid-driven square cavity flow

2018 ◽  
Vol 855 ◽  
pp. 43-66 ◽  
Author(s):  
Ke Wu ◽  
Bruno D. Welfert ◽  
Juan M. Lopez
Keyword(s):  
Stokes Equations ◽  
Complex Dynamics ◽  
Boussinesq Approximation ◽  
Navier Stokes ◽  
Square Cavity ◽  
Navier Stokes Equations ◽  
Stable Temperature ◽  
Wide Range ◽  
Side Walls ◽  
Steady State Solutions

The dynamic response to shear of a fluid-filled square cavity with stable temperature stratification is investigated numerically. The shear is imposed by the constant translation of the top lid, and is quantified by the associated Reynolds number. The stratification, quantified by a Richardson number, is imposed by maintaining the temperature of the top lid at a higher constant temperature than that of the bottom, and the side walls are insulating. The Navier–Stokes equations under the Boussinesq approximation are solved, using a pseudospectral approximation, over a wide range of Reynolds and Richardson numbers. Particular attention is paid to the dynamical mechanisms associated with the onset of instability of steady state solutions, and to the complex and rich dynamics occurring beyond.

Prediction of Reynolds Number Effects on Low-Pressure Turbines Using a High-Order ILES Method

2019 ◽  
Author(s):  
Marc Bolinches-Gisbert ◽  
David Cadrecha Robles ◽  
Roque Corral ◽  
Fernando Gisbert
Keyword(s):  
Boundary Layer ◽  
Reynolds Number ◽  
Fourth Order ◽  
Stokes Equations ◽  
Reynolds Numbers ◽  
Low Pressure ◽  
High Order ◽  
Experimental Results ◽  
Numerical Code ◽  
Design System

Abstract A comprehensive comparison between Implicit Large Eddy Simulations (ILES) and experimental results of a modern highlift low-pressure turbine airfoil has been carried out for an array of Reynolds numbers (Re). Experimental data were obtained in a low-speed linear cascade at the Polithecnic University of Madrid using hot-wire anemometry and LDV. The numerical code is fourth order accurate, both in time and space. The spatial discretization of the compressible Navier-Stokes equations is based on a high-order Flux Reconstruction approach while a fourth order Runge-Kutta method is used to march in time the simulations. The losses, pressure coefficient distributions, and boundary layer and wake velocity profiles have been compared for an array of realistic Reynolds numbers. Moreover, boundary layer and wake velocity fluctuations are compared for the first time with experimental results. It is concluded that the accuracy of the numerical results is comparable to that of the experiments, especially for integral quantities such as the losses or exit angle. Turbulent fluctuations in the suction side boundary layer and the wakes are well predicted also. The elapsed time of the is about 140 hours on 40 Graphics Processor Units. The numerical tool is integrated within an industrial design system and reuses pre- and post-processing tools previously developed for another kind of applications. The trend of the losses with the Reynolds number has a sub-critical regime, where the losses scale with Re−1, and a supercrital regime, where the losses scale with Re−1/2. This trend can be seen both, in the simulations and the experiments.

Fluid Flow and Lateral Fluid Dispersion in Bounded Granular Beds

2002 ◽  
Author(s):  
Azita Soleymani ◽  
Eveliina Takasuo ◽  
Piroz Zamankhan ◽  
William Polashenski
Keyword(s):  
Reynolds Number ◽  
Porous Media ◽  
Numerical Solutions ◽  
Stokes Equations ◽  
Numerical Study ◽  
Fluid Velocity ◽  
Three Dimensional ◽  
Numerical Results ◽  
Transition To Turbulence ◽  
Wide Range

Results are presented from a numerical study examining the flow of a viscous, incompressible fluid through random packing of nonoverlapping spheres at moderate Reynolds numbers (based on pore permeability and interstitial fluid velocity), spanning a wide range of flow conditions for porous media. By using a laminar model including inertial terms and assuming rough walls, numerical solutions of the Navier-Stokes equations in three-dimensional porous packed beds resulted in dimensionless pressure drops in excellent agreement with those reported in a previous study (Fand et al., 1987). This observation suggests that no transition to turbulence could occur in the range of Reynolds number studied. For flows in the Forchheimer regime, numerical results are presented of the lateral dispersivity of solute continuously injected into a three-dimensional bounded granular bed at moderate Peclet numbers. Lateral fluid dispersion coefficients are calculated by comparing the concentration profiles obtained from numerical and analytical methods. Comparing the present numerical results with data available in the literature, no evidence has been found to support the speculations by others for a transition from laminar to turbulent regimes in porous media at a critical Reynolds number.

Supercritical Reynolds number simulation for two-dimensional flow over circular cylinders

1975 ◽  
Vol 70 (3) ◽  
pp. 529-542 ◽  
Author(s):  
Edmond Szechenyi
Keyword(s):  
Surface Roughness ◽  
Reynolds Number ◽  
Reynolds Numbers ◽  
Practical Significance ◽  
Circular Cylinders ◽  
Wind Tunnels ◽  
Bluff Bodies ◽  
Wide Range ◽  
Two Dimensional Flow ◽  
Different Diameters

In wind-tunnel tests on bluff bodies the Reynolds number is often limited to values that are very much smaller than those of the flows being simulated. In such cases the experiments may have no practical significance whatsoever since both the fluctuating and the steady aerodynamic phenomena can vary considerably with Reynolds number.This difficulty was encountered in an investigation of supercritical incompressible flow over cylinders, and an attempt at artificially increasing the Reynolds number by means of surface roughness was made. In order to evaluate this simulation technique, the influence of various grades of surface roughness on the aerodynamic forces acting on cylinders of different diameters was studied over a wide range of Reynolds numbers in two very different wind tunnels. The results allow very positive conclusions to be drawn.

Sign in / Sign up

Export Citation Format

Share Document