The effect of inertia and vertical confinement on the flow past a circular cylinder in a Hele-Shaw configuration

2022 ◽  
Vol 934 ◽  
Author(s):  
C.A. Klettner ◽  
F.T. Smith

The Poiseuille flow (centreline velocity $U_c$ ) of a fluid (kinematic viscosity $\nu$ ) past a circular cylinder (radius $R$ ) in a Hele-Shaw cell (height $2h$ ) is traditionally characterised by a Stokes flow ( $\varLambda =(U_cR/\nu )(h/R)^2 \ll 1$ ) through a thin gap ( $\epsilon =h/R \ll 1$ ). In this work we use asymptotic methods and direct numerical simulations to explore the parameter space $\varLambda$ – $\epsilon$ when these conditions are not met. Starting with the Navier–Stokes equations and increasing $\varLambda$ (which corresponds to increasing inertial effects), four successive regimes are identified, namely the linear regime, nonlinear regimes I and II in the boundary layer (the ‘ inner’ region) and a nonlinear regime III in both the inner and outer region. Flow phenomena are studied with extensive comparisons made between reduced calculations, direct numerical simulations and previous analytical work. For $\epsilon =0.01$ , the limiting condition for a steady flow as $\varLambda$ is increased is the instability of the Poiseuille flow. However, for larger $\epsilon$ , this limit is at a much higher $\varLambda$ , resulting in a laminar separation bubble, of size ${O}(h)$ , forming for a certain range of $\epsilon$ at the back of the cylinder, where the azimuthal location was dependent on $\epsilon$ . As $\epsilon$ is increased to approximately 0.5, the secondary flow becomes increasingly confined adjacent to the sidewalls. The results of the analysis and numerical simulations are summarised in a plot of the parameter space $\varLambda$ – $\epsilon$ .

Author(s):  
Dhanush Vittal Shenoy ◽  
Mostafa Safdari Shadloo ◽  
Jorge Peixinho ◽  
Abdellah Hadjadj

Purpose Fluid flows in pipes whose cross-sectional area are increasing in the stream-wise direction are prone to separation of the recirculation region. This paper aims to investigate such fluid flow in expansion pipe systems using direct numerical simulations. The flow in circular diverging pipes with different diverging half angles, namely, 45, 26, 14, 7.2 and 4.7 degrees, are considered. The flow is fed by a fully developed laminar parabolic velocity profile at its inlet and is connected to a long straight circular pipe at its downstream to characterise recirculation zone and skin friction coefficient in the laminar regime. The flow is considered linearly stable for Reynolds numbers sufficiently below natural transition. A perturbation is added to the inlet fully developed laminar velocity profile to test the flow response to finite amplitude disturbances and to characterise sub-critical transition. Design/methodology/approach Direct numerical simulations of the Navier–Stokes equations have been solved using a spectral element method. Findings It is found that the onset of disordered motion and the dynamics of the localised turbulence patch are controlled by the Reynolds number, the perturbation amplitude and the half angle of the pipe. Originality/value The authors clarify different stages of flow behaviour under the finite amplitude perturbations and shed more light to flow physics such as existence of Kelvin–Helmholtz instabilities as well as mechanism of turbulent puff shedding in diverging pipe flows.


Author(s):  
Joshua R. Brinkerhoff ◽  
Metin I. Yaras

This paper describes numerical simulations of the instability mechanisms in a separation bubble subjected to a three-dimensional freestream pressure distribution. Two direct numerical simulations are performed of a separation bubble with laminar separation and turbulent reattachment under low freestream turbulence at flow Reynolds numbers and streamwise pressure distributions that approximate the conditions encountered on the suction side of typical low-pressure gas-turbine blades with blade sweep angles of 0° and 45°. The three-dimensional pressure field in the swept configuration produces a crossflow-velocity component in the laminar boundary layer upstream of the separation point that is unstable to a crossflow instability mode. The simulation results show that crossflow instability does not play a role in the development of the boundary layer upstream of separation. An increase in the amplification rate and most amplified disturbance frequency is observed in the separated-flow region of the swept configuration, and is attributed to boundary-layer conditions at the point of separation that are modified by the spanwise pressure gradient. This results in a slight upstream movement of the location where the shear layer breaks down to small-scale turbulence and modifies the turbulent mixing of the separated shear layer to yield a downstream shift in the time-averaged reattachment location. The results demonstrate that although crossflow instability does not appear to have a noticeable effect on the development of the transitional separation bubble, the 3D pressure field does indirectly alter the separation-bubble development by modifying the flow conditions at separation.


2010 ◽  
Vol 39 (6) ◽  
pp. 926-935 ◽  
Author(s):  
Per-Olov Åsén ◽  
Gunilla Kreiss ◽  
Dietmar Rempfer

1991 ◽  
Vol 227 ◽  
pp. 473-493 ◽  
Author(s):  
S. Sarkar ◽  
G. Erlebacher ◽  
M. Y. Hussaini ◽  
H. O. Kreiss

It is shown that the dilatational terms that need to be modelled in compressible turbulence include not only the pressure-dilatation term but also another term - the compressible dissipation. The nature of the compressible velocity field, which generates these dilatational terms, is explored by asymptotic analysis of the compressible Navier-Stokes equations in the case of homogeneous turbulence. The lowest-order equations for the compressible field are solved and explicit expressions for some of the associated one-point moments are obtained. For low Mach numbers, the compressible mode has a fast timescale relative to the incompressible mode. Therefore, it is proposed that, in moderate Mach number homogeneous turbulence, the compressible component of the turbulence is in quasi-equilibrium with respect to the incompressible turbulence. A non-dimensional parameter which characterizes this equilibrium structure of the compressible mode is identified. Direct numerical simulations (DNS) of isotropic, compressible turbulence are performed, and their results are found to be in agreement with the theoretical analysis. A model for the compressible dissipation is proposed; the model is based on the asymptotic analysis and the direct numerical simulations. This model is calibrated with reference to the DNS results regarding the influence of compressibility on the decay rate of isotropic turbulence. An application of the proposed model to the compressible mixing layer has shown that the model is able to predict the dramatically reduced growth rate of the compressible mixing layer.


1999 ◽  
Vol 393 ◽  
pp. 257-308 ◽  
Author(s):  
F. S. GODEFERD ◽  
L. LOLLINI

The goal of this work is to analyse how solid body rotation affects forced turbulence enclosed within solid boundaries, and to compare it to results of the experiment performed by Hopfinger et al. (1982). In order to identify various mechanisms associated with rotation, confinement, and forcing, a numerical pseudo-spectral code is used for performing direct numerical simulations. The geometry is simplified with respect to the experimental one. First, we are able to reproduce the linear regime, as propagating inertial waves that undergo reflections at the walls. Second, the Ekman pumping phenomenon, proportional to the rotation rate, is identified in freely decaying turbulence, for which the evolution of the flow bounded by walls is compared to the evolution of unbounded homogeneous turbulence. Finally we introduce a local forcing on a plane in physical space, for simulating the effect of an oscillating grid, so that diffusive turbulence is created, and we examine the structuring of the flow under the combination of the linear and nonlinear mechanisms. A transition to an almost two-dimensional state is shown to occur between the region close to the forcing and an outer region in which vortices appear, the number of which depends on the Reynolds and Rossby numbers. In this region, the anisotropy of turbulence is examined, and the numerical predictions are shown to reproduce many of the most important features present in the experimental flow.


2016 ◽  
Vol 798 ◽  
pp. 853-888 ◽  
Author(s):  
Andrew J. Hogg ◽  
Mohamad M. Nasr-Azadani ◽  
Marius Ungarish ◽  
Eckart Meiburg

Gravitationally driven motion arising from a sustained constant source of dense fluid in a horizontal channel is investigated theoretically using shallow-layer models and direct numerical simulations of the Navier–Stokes equations, coupled to an advection–diffusion model of the density field. The influxed dense fluid forms a flowing layer underneath the less dense fluid, which initially filled the channel, and in this study its speed of propagation is calculated; the outflux is at the end of the channel. The motion, under the assumption of hydrostatic balance, is modelled using a two-layer shallow-water model to account for the flow of both the dense and the overlying less dense fluids. When the relative density difference between the fluids is small (the Boussinesq regime), the governing shallow-layer equations are solved using analytical techniques. It is demonstrated that a variety of flow-field patterns are feasible, including those with constant height along the length of the current and those where the height varies continuously and discontinuously. The type of solution realised in any scenario is determined by the magnitude of the dimensionless flux issuing from the source and the source Froude number. Two important phenomena may occur: the flow may be choked, whereby the excess velocity due to the density difference is bounded and the height of the current may not exceed a determined maximum value, and it is also possible for the dense fluid to completely displace all of the less dense fluid originally in the channel in an expanding region close to the source. The onset and subsequent evolution of these types of motions are also calculated using analytical techniques. The same range of phenomena occurs for non-Boussinesq flows; in this scenario, the solutions of the model are calculated numerically. The results of direct numerical simulations of the Navier–Stokes equations are also reported for unsteady two-dimensional flows in which there is an inflow of dense fluid at one end of the channel and an outflow at the other end. These simulations reveal the detailed mechanics of the motion and the bulk properties are compared with the predictions of the shallow-layer model to demonstrate good agreement between the two modelling strategies.


2018 ◽  
Vol 848 ◽  
pp. 648-659 ◽  
Author(s):  
Yantao Yang ◽  
Roberto Verzicco ◽  
Detlef Lohse

We conduct direct numerical simulations for turbulent Rayleigh–Bénard (RB) convection, driven simultaneously by two scalar components (say, temperature and concentration) with different molecular diffusivities, and measure the respective fluxes and the Reynolds number. To account for the results, we generalize the Grossmann–Lohse theory for traditional RB convection (Grossmann & Lohse, J. Fluid Mech., vol. 407, 2000, pp. 27–56; Phys. Rev. Lett., vol. 86 (15), 2001, pp. 3316–3319; Stevens et al., J. Fluid Mech., vol. 730, 2013, pp. 295–308) to this two-scalar turbulent convection. Our numerical results suggest that the generalized theory can successfully capture the overall trends for the fluxes of two scalars and the Reynolds number without introducing any new free parameters. In fact, for most of the parameter space explored here, the theory can even predict the absolute values of the fluxes and the Reynolds number with good accuracy. The current study extends the generality of the Grossmann–Lohse theory in the area of buoyancy-driven convection flows.


2020 ◽  
Vol 142 (5) ◽  
Author(s):  
Chunning Ji ◽  
Zhimeng Zhang ◽  
Dong Xu ◽  
Narakorn Srinil

Abstract Understanding hydrodynamics of a free-spanning pipeline subjected to omni-directional flows is important to engineering design. In this study, horizontally oblique flows past a three-dimensional circular cylinder in the vicinity of a plane boundary are numerically investigated using direct numerical simulations. Parametric studies are carried out at the normal Reynolds number of 500, a fixed gap-to-diameter ratio of 0.8 and five flow inclination angles (α) ranging from 0 deg to 60 deg with an increment of 15 deg. Two distinct vortex-shedding modes are observed: parallel (α ≤ 15 deg) and oblique (α ≥ 30 deg) vortex shedding. The wake evolution is further divided into two or three stages depending on α. The occurrence of the oblique vortex shedding is accompanied by the base pressure gradient along the cylinder span and the resultant axial flows near the cylinder base. The total hydrodynamic drag and lift force coefficients decrease from being the parallel mode to the oblique mode, owing to the intensified three-dimensionality of wake flows and the phase differences in the spanwise vortex shedding. The independence principle (IP) is found to be valid in predicting hydrodynamic forces and wake patterns when α ≤ 15 deg. This IP might produce unacceptable errors when α > 15 deg. In comparison with the mean drag force, the fluctuating lift force is more sensitive to the inclination angle. The IP validity range is substantially smaller than that in the case of flow past a wall-free cylinder. Such finding would be practically useful for vortex-induced vibration prediction.


Sign in / Sign up

Export Citation Format

Share Document