Showing papers on "Reynolds number published in 2005"

A scale-dependent Lagrangian dynamic model for large eddy simulation of complex turbulent flows

Abstract: A scale-dependent dynamic subgrid model based on Lagrangian time averaging is proposed and tested in large eddy simulations (LES) of high-Reynolds number boundary layer flows over homogeneous and heterogeneous rough surfaces. The model is based on the Lagrangian dynamic Smagorinsky model in which required averages are accumulated in time, following fluid trajectories of the resolved velocity field. The model allows for scale dependence of the coefficient by including a second test-filtering operation to determine how the coefficient changes as a function of scale. The model also uses the empirical observation that when scale dependence occurs (such as when the filter scale approaches the limits of the inertial range), the classic dynamic model yields the coefficient value appropriate for the test-filter scale. Validation tests in LES of high Reynolds number, rough wall, boundary layer flow are performed at various resolutions. Results are compared with other eddy-viscosity subgrid-scale models. Unlike the...

Analysis of the microfluid flow in an evaporating sessile droplet.

Abstract: The axisymmetric time-dependent flow field in an evaporating sessile droplet whose contact line is pinned is studied numerically and using an analytical lubrication theory with a zero-shear-stress boundary condition on the free surface of the droplet at low capillary and Reynolds numbers. A finite element algorithm is developed to solve simultaneously the vapor concentration and flow field in the droplet under conditions of slow evaporation. The finite element solution confirms the accuracy of the lubrication solution, especially when terms of higher order in the droplet flatness ratio (the ratio of droplet height to radius, h/R) are included in the lubrication theory to account more accurately for the singular flow near the contact line.

Numerical and experimental investigation of flow and scour around a circular pile

Abstract: The flow around a vertical circular pile exposed to a steady current is studied numerically and experimentally. The numerical model is a three-dimensional model. The model validation was achieved against new experimental data (which include two-component laser-Doppler anemometry (LDA) flow measurements and the hot-film bed shear stress measurements, and reported in the present paper) and the data of others, and a k-w turbulence model was used for closure. The model does not have a free-surface facility and therefore is applicable only to cases where the Froude number is small (Fr < 0(0.2)). The flow model was used to study the horseshoe vortex and lee-wake vortex flow processes around the pile. The influence on the horseshoe vortex of three parameters, namely the boundary-layer thickness, the Reynolds number and the bed roughness, was investigated. In the latter investigation, the steady solution of the model was chosen. A study of the influence of the unsteady solution on the previously mentioned flow processes was also carried out. The ranges of the parameters covered in the numerical simulations are: The boundary-layer-thickness-to-pile-diameter ratio is varied from 2 x 10 -2 to 10 2 , the pile Reynolds number from 10 2 to 2 × 10 6 , and the pile diameter-to-roughness ratio from 2 to about 10 3 . The amplification of the bed shear stress around the pile (including the areas under the horseshoe vortex and the lee-wake region) was obtained for various values of the previously mentioned parameters. The steady-state flow model was coupled with a morphologic model to calculate scour around a vertical circular pile exposed to a steady current in the case of non-cohesive sediment. The morphologic model includes (i) a two-dimensional bed load sediment-transport description, and (ii) a description of surface-layer sand slides for bed slopes exceeding the angle of repose. The results show that the present numerical simulation captures all the main features of the scour process. The equilibrium scour depth obtained from the simulation agrees well with the experiments for the upstream scour hole. Some discrepancy (up to 30 %) was observed, however, for the downstream scour hole. The calculations show that the amplification of the bed shear stress around the pile in the equilibrium state of the scour process is reduced considerably with respect to that experienced at the initial stage where the bed is plane.

Numerical solutions of 2‐D steady incompressible driven cavity flow at high Reynolds numbers

Abstract: SUMMARY Numerical calculations of the 2-D steady incompressible driven cavity flow are presented. The NavierStokes equations in streamfunction and vorticity formulation are solved numerically using a fine uniform grid mesh of 601 × 601. The steady driven cavity solutions are computed for Re ≤ 21,000 with a maximum absolute residuals of the governing equations that were less than 10 −10 . A new quaternary vortex at the bottom left corner and a new tertiary vortex at the top left corner of the cavity are observed in the flow field as the Reynolds number increases. Detailed results are presented and comparisons are made with benchmark solutions found in the literature.

Lattice-Boltzmann simulations of low-Reynolds-number flow past mono- and bidisperse arrays of spheres: results for the permeability and drag force

Abstract: We report on lattice-Boltzmann simulations of slow fluid flow past mono- and bidisperse random arrays of spheres. We have measured the drag force on the spheres for a range of diameter ratios, mass fractions and packing fractions; in total, we studied 58 different parameter sets. Our simulation data for the permeability agrees well with previous simulation results and the experimental findings. On the basis of our data for the individual drag force, we have formulated new drag force relations for both monodisperse and polydisperse systems, based on the Carman?Kozeny equations; the average deviation of our binary simulation data with the new relation is less than 5%. We expect that these new relations will significantly improve the numerical modelling of gas?solid systems with polydisperse particles, in particular with respect to mixing and segregation phenomena. For binary systems with large diameter ratios (1:4), the individual drag force on a particle, as calculated from our relations, can differ by up to a factor of five compared with predictions presently favoured in chemical engineering.

Componentwise energy amplification in channel flows

Abstract: We study the linearized Navier–Stokes (LNS) equations in channel flows from an input–output point of view by analysing their spatio-temporal frequency responses. Spatially distributed and temporally varying body force fields are considered as inputs, and components of the resulting velocity fields are considered as outputs into these equations. We show how the roles of Tollmien–Schlichting (TS) waves, oblique waves, and streamwise vortices and streaks in subcritical transition can be explained as input–output resonances of the spatio-temporal frequency responses. On the one hand, we demonstrate the effectiveness of input field components, and on the other, the energy content of velocity perturbation components. We establish that wall-normal and spanwise forces have much stronger influence on the velocity field than streamwise force, and that the impact of these forces is most powerful on the streamwise velocity component. We show this using the relative scaling of the different input–output system components with the Reynolds number. We further demonstrate that for the streamwise constant perturbations, the spanwise force localized near the lower wall has, by far, the strongest effect on the evolution of the velocity field. In this paper, we analyse the dynamical properties of the Navier–Stokes (NS) equations with spatially distributed and temporally varying body force fields. These fields are considered as inputs, and different combinations of the resulting velocity fields are considered as outputs. This input–output analysis can in principle be done in any geometry and for the full nonlinear NS equations. In such generality, however, it is difficult to obtain useful results. We therefore concentrate on the geometry of channel flows, and the input–output dynamics of the linearized Navier–Stokes (LNS)

Formation Criterion for Synthetic Jets

Abstract: Af ormation criterion for synthetic jets is proposed and validated. A synthetic jet actuator is a zero-net mass-flux device that imparts momentum to its surroundings. Je tf ormation is defined as the appearance of a time-averaged outward velocity along the jet axis and corresponds to the generation and subsequent convection or escape of a vortex ring. It is shown that over a wide range of operating conditions synthetic jet formation is governed by the jet Strouhal number Sr (or Reynolds number Re and Stokes number S). Both numerical simulations and experiments are performed to supplement available two-dimensional and axisymmetric synthetic jet formation data in the literature. The data support the jet formation criterion 1/Sr = Re/S 2 > K, where the constant K is approximately 1 and 0.16 for two-dimensional and axisymmetric synthetic jets, respectively. In addition, the dependence of the constant K on the normalized radius of curvature of a rounded orifice or slot is addressed. The criterion is expected to serve as a useful design guide for synthetic jet formation in flow control, heat transfer, and acoustic liner applications, in which a stronger jet is synonymous with increased momentum transfer, vorticity generation, and acoustic nonlinearities.

Topology optimization of channel flow problems

Abstract: This paper describes a topology design method for simple two-dimensional flow problems. We consider steady, incompressible laminar viscous flows at low-to-moderate Reynolds numbers. This makes the flow problem nonlinear and hence a nontrivial extension of the work of Borrvall and Petersson (2003).Further, the inclusion of inertia effects significantly alters the physics, enabling solutions of new classes of optimization problems, such as velocity-driven switches, that are not addressed by the earlier method. Specifically, we determine optimal layouts of channel flows that extremize a cost function which measures either some local aspect of the velocity field or a global quantity, such as the rate of energy dissipation. We use the finite element method to model the flow, and we solve the optimization problem with a gradient-based math-programming algorithm that is driven by analytical sensitivities. Our target application is optimal layout design of channels in fluid network systems. Using concepts borrowed from topology optimization of compliant mechanisms in solid mechanics, we introduce a method for the synthesis of fluidic components, such as switches, diodes, etc.

Highly resolved large-eddy simulation of separated flow in a channel with streamwise periodic constrictions

Abstract: High-resolution large-eddy simulation is used to investigate the mean and turbulence properties of a separated flow in a channel constricted by periodically distributed hill-shaped protrusions on one wall that obstruct the channel by 33% of its height and are arranged 9 hill heights apart. The geometry is a modification of an experimental configuration, the adaptation providing an extended region of post-reattachment recovery and allowing high-quality simulations to be performed at acceptable computing costs. The Reynolds number, based on the hill height and the bulk velocity above the crest is 10595. The simulated domain is streamwise as well as spanwise periodic, extending from one hill crest to the next in the streamwise direction and over 4.5 hill heights in the spanwise direction. This arrangement minimizes uncertainties associated with boundary conditions and makes the flow an especially attractive generic test case for validating turbulence closures for statistically two-dimensional separation. The emphasis of the study is on elucidating the turbulence mechanisms associated with separation, recirculation reattachment, acceleration and wall proximity. Hence, careful attention has been paid to resolution, and a body-fitted, low-aspect-ratio, nearly orthogonal numerical grid of close to 5 million nodes has been used. Unusually, the results of two entirely independent simulations with different codes for identical flow and numerical conditions are compared and shown to agree closely. Results are included for mean velocity, Reynolds stresses, anisotropy measures, spectra and budgets for the Reynolds stresses. Moreover, an analysis of structural characteristics is undertaken on the basis of instantaneous realizations, and links to features observed in the statistical results are identified and interpreted. Among a number of interesting features, a distinct ‘splatting’ of eddies on the windward hill side following reattachment is observed, which generates strong spanwise fluctuations that are reflected, statistically, by the spanwise normal stress near the wall exceeding that of the streamwise stress by a substantial margin, despite the absence of spanwise strain.

Large-Eddy Simulations of Turbulence

Abstract: Direct-numerical simulations of turbulence (DNS) consist in solving explicitly all the scales of motion, from the largest l i to the Kolmogorov dissipative scale l d. It is well known from the statistical theory of turbulence that l i/l d scales like R l 3/4, where R l is the large-scale Reynolds number u’l i/ν based upon the rms velocity fluctuation u′. Therefore, the total number of degrees of freedom necessary to represent the whole span of scales of a three-dimensional turbulent flow is of the order of R l 9/4 in three dimensions. In the presence of obstacles, around a wing or a fuselage for instance, and if one wants to simulate three-dimensionally all motions ranging from the viscous thickness δ v = ν/υ * ≈ 10-6 m up to 10 m, it would be necessary to put 1021 modes on the computer. Right now, the calculations done to the expense of not excessive computing times on the biggest machines take about 2 × 107 grid points, which is very far from the above estimation. Even with the unprecedented improvement of scientific computers, it may take several tenths of years (if it becomes ever possible) before DNS permit to simulate situations at Reynolds numbers comparable to those encountered in natural conditions.

Parameterization of turbulent fluxes and scales using homogeneous sheared stably stratified turbulence simulations

Abstract: Laboratory experiments on stably stratified grid turbulence have suggested that turbulent diffusivity , suggesting that such Reynolds–Richardson number or Reynolds–Froude number aggregates are more descriptive of stratified turbulent flow conditions than the conventional reliance on Richardson number alone.

The inertio-elastic planar entry flow of low-viscosity elastic fluids in micro-fabricated geometries

Abstract: The non-Newtonian flow of dilute aqueous polyethylene oxide (PEO) solutions through micro-fabricated planar abrupt contraction-expansions is investigated. The small lengthscales and high deformation rates in the contraction throat lead to significant extensional flow effects even with dilute polymer solutions having time constants on the order of milliseconds. By considering the definition of the elasticity number, El = Wi/Re, we show that the lengthscale of the geometry is key to the generation of strong viscoelastic effects, such that the same flow behaviour cannot be reproduced using the equivalent macro-scale geometry using the same fluid. We observe significant vortex growth upstream of the contraction plane, which is accompanied by an increase of more than 200% in the dimensionless extra pressure drop across the contraction. Streak photography and video-microscopy using epifluorescent particles shows that the flow ultimately becomes unstable and three-dimensional. The moderate Reynolds numbers (0.44 ≤ Re ≤ 64) associated with these high Weissenberg number (0 ≤ Wi ≤ 548) micro-fluidic flows results in the exploration of new regions of the Re-Wi parameter space in which the effects of both elasticity and inertia can be observed. Understanding such interactions will be increasingly important in micro-fluidic applications involving complex fluids and can best be interpreted in terms of the elasticity number, El = Wi/Re, which is independent of the flow kinematics and depends only on the fluid rheology and the characteristic size of the device.

An experimental investigation of swirling jets

Abstract: The 'plug' flow emerging from a long rotating tube into a large stationary reservoir has been used in an experimental investigation of centrifugally unstable swirling jets. A moderate Reynolds number, Re = 1000, was studied extensively, and swirl numbers, S, the ratio of nozzle exit rotating speed to the mean mass axial velocity, were in the range 0-1.1. Four regimes were covered: non-swirling jets with S=0, weakly swirling jets in the range 0 S c2 , where S c1 =0.6 and S c2 = 0.88. Particular attention was paid to the dominant role of the underlying vortical flow structures and their dynamic evolution. Kelvin-Helmholtz (K-H) instability in the axial shear layer, generated by the axial velocity, leading to vortex ring formation, dominated non-swirling and weakly swirling jets

Stokes first problem for an oldroyd-b fluid in a porous half space

Abstract: Based on a modified Darcy’s law for a viscoelastic fluid, Stokes’ first problem was extended to that for an Oldroyd-B fluid in a porous half space. By using Fourier sine transform, an exact solution was obtained. In contrast to the classical Stokes’ first problem for a clear fluid, there is a y-dependent steady state solution for an Oldroyd-B fluid in the porous half space, which is a damping exponential function with respect to the distance from the flat plate. The thickness of the boundary layer, which tends to be a limited value, is also different from that of a clear fluid. The effect of viscoelasticity on the unsteady flow in porous media is investigated. It was found if α>1∕4[(αt∕Re)+Re]2, oscillations in velocity occur obviously and the system exhibits viscoelastic behaviors, where α and αt are nondimensional relaxation and retardation times, respectively, Re is Reynold number in porous media. Some previous solutions of Stokes’ first problem corresponding to Maxwell fluid and Newtonian fluid in por...

Chaos and threshold for irreversibility in sheared suspensions

Abstract: Systems governed by time reversible equations of motion often give rise to irreversible behaviour. The transition from reversible to irreversible behaviour is fundamental to statistical physics, but has not been observed experimentally in many-body systems. The flow of a newtonian fluid at low Reynolds number can be reversible: for example, if the fluid between concentric cylinders is sheared by boundary motion that is subsequently reversed, then all fluid elements return to their starting positions. Similarly, slowly sheared suspensions of solid particles, which occur widely in nature and science, are governed by time reversible equations of motion. Here we report an experiment showing precisely how time reversibility6 fails for slowly sheared suspensions. We find that there is a concentration dependent threshold for the deformation or strain beyond which particles do not return to their starting configurations after one or more cycles. Instead, their displacements follow the statistics of an anisotropic random walk. By comparing the experimental results with numerical simulations, we demonstrate that the threshold strain is associated with a pronounced growth in the Lyapunov exponent (a measure of the strength of chaotic particle interactions). The comparison illuminates the connections between chaos, reversibility and predictability.

Experimental support for Townsend’s Reynolds number similarity hypothesis on rough walls

Abstract: The Reynolds number similarity hypothesis of Townsend [The Structure of Turbulent Shear Flow (Cambridge University Press, Cambridge, UK, 1976)] states that the turbulence beyond a few roughness heights from the wall is independent of the surface condition. The underlying assumption is that the boundary layer thickness δ is large compared to the roughness height k. This hypothesis was tested experimentally on two types of three-dimensional rough surfaces. Boundary layer measurements were made on flat plates covered with sand grain and woven mesh roughness in a closed return water tunnel at a momentum thickness Reynolds number Reθ of ∼14000. The boundary layers on the rough walls were in the fully rough flow regime (ks+⩾100) with the ratio of the boundary layer thickness to the equivalent sand roughness height δ∕ks greater than 40. The results show that the mean velocity profiles for rough and smooth walls collapse well in velocity defect form in the overlap and outer regions of the boundary layer. The Reynolds stresses for the two rough surfaces agree well throughout most of the boundary layer and collapse with smooth wall results outside of 3ks. Higher moment turbulence statistics and quadrant analysis also indicate the differences in the rough wall boundary layers are confined to y<5ks. The present results provide support for Townsend’s Reynolds number similarity hypothesis for uniform three-dimensional roughness in flows where δ∕ks⩾40.

Unsteady aerodynamics of fluttering and tumbling plates

Abstract: We investigate the aerodynamics of freely falling plates in a quasi-two-dimensional flow at Reynolds number of , which is typical for a leaf or business card falling in air. We quantify the trajectories experimentally using high-speed digital video at sufficient resolution to determine the instantaneous plate accelerations and thus to deduce the instantaneous fluid forces. We compare the measurements with direct numerical solutions of the two-dimensional Navier–Stokes equation. Using inviscid theory as a guide, we decompose the fluid forces into contributions due to acceleration, translation, and rotation of the plate. For both fluttering and tumbling we find that the fluid circulation is dominated by a rotational term proportional to the angular velocity of the plate, as opposed to the translational velocity for a glider with fixed angle of attack. We find that the torque on a freely falling plate is small, i.e. the torque is one to two orders of magnitude smaller than the torque on a glider with fixed angle of attack. Based on these results we revise the existing ODE models of freely falling plates. We get access to different kinds of dynamics by exploring the phase diagram spanned by the Reynolds number, the dimensionless moment of inertia, and the thickness-to-width ratio. In agreement with previous experiments, we find fluttering, tumbling, and apparently chaotic motion. We further investigate the dependence on initial conditions and find brief transients followed by periodic fluttering described by simple harmonics and tumbling with a pronounced period-two structure. Near the cusp-like turning points, the plates elevate, a feature which would be absent if the lift depended on the translational velocity alone.

Friction factors for pipe flow of xanthan-based concentrates of fire fighting foams

Abstract: In this paper we develop a friction factor correlation to predict the pressure drop during pumping and induction of concentrates of fire fighting foams containing around 1% of xanthan gum. Such concentrates are highly elastic, display small yield stress and exhibit significant thinning upon shearing. We demonstrate that in the turbulent regime, the Blasius equation normally used for Newtonian fluids seems to correlate well the friction factor with the Metzner-Reed Reynolds number. Our development provides an example of how the methodology used to develop the friction factor correlation can be applied to analyse a set of experimental data to verify its internal consistency. The friction factors developed in the paper can be applied to other foam concentrates that include a similar content of xanthan gum in their formulation, to predict pressure drop as a function of a flow rate and pipe diameter, provided that there exists an appropriate viscosity model. Subsequently, the paper presents experimental measurements of apparent viscosity for one foam concentrate and develops relevant viscosity models. We observe that the rheology of the concentrate is governed by the behaviour of xanthan gum. Although, the foam concentrate considered in the article is yield pseudo plastic (i.e., it follows the Herschel-Bulkley model), for the shear stresses normally encountered during pipe flow, the viscosity of the material can be described by a power law model. Over the temperature range of between 0 and 40 o C, the apparent viscosity displays only a weak dependence on temperature. Subsequent calculations of pressure drop with temperature demonstrate minor variation in pressure drop with temperature, but only in the turbulent flow regime. This suggests that induction systems intended to operate under widely varying temperature conditions should be designed to function in the turbulent flow regime.

The Development Lengths of Laminar Pipe and Channel Flows

Abstract: The authors’ research work into fully developed pulsating and oscillating laminar pipe and channel flows raised questions regarding the development length of the corresponding steady flow. For this development length, i.e., the distance from the entrance of the pipe to the axial position where the flow reaches the parabolic velocity profile of the Hagen-Poiseuille flow, a wide range of contradictory data exists. This is shown through a short review of the existing literature. Superimposed diffusion and convection, together with order of magnitude considerations, suggest that the normalized development length can be expressed as L∕D=C0+C1Re and for Re→0 one obtains C0=0.619, whereas for Re→∞ one obtains C1=0.0567. This relationship is given only once in the literature and it is presumed to be valid for all Reynolds numbers. Numerical studies show that it is only valid for Re→0 and Re→∞. The development length of laminar, plane channel flow was also investigated. The authors obtained similar results to those for the pipe flow: L∕D=C0′+C1′; Re, where C0′=0.631 and C1′=0.044. Finally, correlations are given to express L∕D analytically for the entire Re range for both laminar pipe and channel flows.

Wall-induced forces on a rigid sphere at finite Reynolds number

Abstract: We perform direct numerical simulations of a rigid sphere translating parallel to a flat wall in an otherwise quiescent ambient fluid. A spectral element method is employed to perform the simulations with high accuracy. For above about 100. Detailed analysis of the flow field around the sphere suggests that this increase is due to an imperfect bifurcation resulting in the formation of a double-threaded wake vortical structure. In addition to a non-rotating sphere, we also simulate a freely rotating sphere in order to assess the importance of free rotation on the translational motion of the sphere. We observe the effect of sphere rotation on lift and drag forces to be small. We also explore the effect of the wall on the onset of unsteadiness.

Experimental study of eddy structures in a turbulent boundary layer using particle image velocimetry

Abstract: Particle image velocimetry experiments have been performed in a turbulent boundary-layer wind tunnel in order to study the coherent structures taking part in the generation and preservation of wall turbulence. The particular wind tunnel used is suitable for high-resolution experiments ($delta gt 0.3$ m) at high Reynolds numbers (up to $R_{theta} = 19,000$ in the present results). Eddy structures were identified in instantaneous velocity maps in order to determine their mean characteristics and possible relationships between these structures. In the logarithmic region, the results show that the observed eddy structures appear to organize like elongated vortices, tilted downstream, mainly at an angle of about 45° and having a cane shape. The characteristics of these vortices appear here to be universal in wall units for $R_{theta},{leq},19,000$. They seem to find their origin at a wall distance of about 25 wall units as quasi-streamwise vortices and to migrate away from the wall while tilting to form a head and a leg. Away from the wall, their radius increases and their vorticity decreases very slowly so that their circulation is nearly constant. Near the wall, the picture obtained is in fair agreement with existing models. The analysis of the results indicates a universality of the buffer-layer mechanism, even at low Reynolds number, and a sensitivity of the logarithmic region to low-Reynolds-number effects.