The computation of the eddy along the upper wall in the three‐dimensional flow over a backward‐facing step

A three‐dimensional laminar flow over a backward‐facing step is studied as a numerical experiment by solving the steady‐state, isothermal and incompressible Navier–Stokes equations using two different finite element codes. The Reynolds number ranges from 100 to 1050. The expansion ratio is 1:1.94, and the aspect ratio is 1:36.7. The numerical experiment reveals both eddies along the lower and upper walls downstream of the step. Results of computations regarding positions of detachment of the eddy along the upper wall and positions of reattachments of the eddies along both the lower and upper walls are tabulated along with positions and magnitudes of global extrema of shear rate within the eddies. The wall effects are shown by calculating streamlines along planes parallel/normal to the lateral walls of the domain and depicting how the streamlines are distorted close to the walls and how they assume a two‐dimensional configuration in the plane of symmetry. Comparisons are made with available numerical results and laboratory measurements. Copyright © 2011 John Wiley & Sons, Ltd.     

A numerical investigation of the influence of the aspect ratio in three‐dimensional separated flows

The influence of aspect ratio in three‐dimensional, numerical experiments of separated flows is studied in the case of the backward‐facing step at Reynolds numbers 600, 800, and 950. The computational domain is designed as an actual laboratory experiment. The governing equations are the steady state, isothermal, and incompressible Navier–Stokes equations. The expansion ratio of the computational domain is 1:2. The aspect ratio varies from 1:10 to 1:40. The results of the computations focus on the spanwise variations of the length and the strength of the two eddies along the lower and upper wall. It is concluded that both numerical and laboratory experiments should be designed with an aspect ratio of at least 1:20, if only the accuracy of the position of the detachment and the re‐attachment points matters. If the accuracy of the shear‐stress distributions is also taken into account, then an aspect ratio of at least 1:30 should be chosen. Finally, if the magnitudes of the vortex centers are also considered, then only the aspect ratio of 1:40 qualifies for a realization of two‐dimensional flow conditions in the plane of symmetry. This is contrary to the common practice in the field, at least from the side of laboratory experiments, where an aspect ratio of 1:10 is still considered adequate for this purpose. Copyright © 2011 John Wiley & Sons, Ltd.     

A Fourier–Chebyshev spectral collocation method for simulating flow past spheres and spheroids

An accurate Fourier–Chebyshev spectral collocation method has been developed for simulating flow past prolate spheroids. The incompressible Navier–Stokes equations are transformed to the prolate spheroidal co‐ordinate system and discretized on an orthogonal body fitted mesh. The infinite flow domain is truncated to a finite extent and a Chebyshev discretization is used in the wall‐normal direction. The azimuthal direction is periodic and a conventional Fourier expansion is used in this direction. The other wall‐tangential direction requires special treatment and a restricted Fourier expansion that satisfies the parity conditions across the poles is used. Issues including spatial and temporal discretization, efficient inversion of the pressure Poisson equation, outflow boundary condition and stability restriction at the pole are discussed. The solver has been validated primarily by simulating steady and unsteady flow past a sphere at various Reynolds numbers and comparing key quantities with corresponding data from experiments and other numerical simulations. Copyright © 1999 John Wiley & Sons, Ltd.     

Analytic Approximate Solution for a Flow of a Second-Grade Viscoelastic Fluid in a Converging Porous Channel

The problem of a two-dimensional steady flow of a second-grade fluid in a converging porous channel is considered. It is assumed that the fluid is injected into the channel through one wall and sucked from the channel through the other wall at the same velocity, which is inversely proportional to the distance along the wall from the channel origin. The equations governing the flow are reduced to ordinary differential equations. The boundary-value problem described by the latter equations is solved by the homotopy perturbation method. The effects of the Reynolds and crossflow Reynolds number on the flow characteristics are examined.     

Temporal and spatial transients in turbulent boundary layer flow over an oscillating wall

Direct numerical simulations have been performed to study the effect of an oscillating segment of the wall on a turbulent boundary layer flow. Two different oscillation amplitudes with equal oscillation period have been used, which allows a direct comparison between a relatively weak and strong forcing of the flow. The weaker forcing results in 18% drag reduction while the stronger forcing, with twice the amplitude, yields 29% drag reduction. The downstream development of the drag reduction is compared with earlier simulations and experiments. In addition, a simulation with identical oscillation parameters as in previous numerical and experimental investigations allows for an estimation of the effect of the Reynolds number on the drag reduction.Reductions in the Reynolds stresses and the important role that the edge of the Stokes layer has is explained.An estimation of the idealized power consumption shows that a positive energy budget is only possible for the weaker wall velocity case.Spatial and temporal transients are investigated and a transformation between spatial and temporal coordinates via a convection velocity is shown to facilitate a comparison between the two transients in a consistent manner. The streamwise shear exhibits a similar monotonic behavior in the spatial and temporal transients, while the non-monotinic temporal transient of the longitudinal Reynolds stress has no counterpart in the spatial development. Furthermore, the evolution in time of the spanwise Reynolds stress is very similar to previously reported channel flow data.The instantaneous spanwise velocity profile (only averaged in the homogeneous spanwise direction) will for the first time be presented from a boundary layer over an oscillating wall, and comparisons with the analytical solution to the laminar Navier–Stokes equations show very good agreement.     

Turbulence in a three‐dimensional wall‐bounded shear flow

A new turbulent flow with distinct three‐dimensional characteristics has been designed in order to study the impact of mean‐flow skewing on the turbulent coherent vortices and Reynolds‐averaged statistics. The skewing of a unidirectional plane Couette flow was achieved by means of a spanwise pressure gradient. Direct numerical simulations of the statistically steady Couette–Poiseuille flow enabled in‐depth explorations of the turbulence field in the skewed flow. The imposition of a modest spanwise gradient turned the mean flow about 8° away from the original Couette flow direction and this turning angle remained nearly the same over the entire cross section. Nevertheless, a substantial non‐alignment between the turbulent shear stress angle and the mean velocity gradient angle was observed. The structure parameter turned out to slightly exceed that in the pure Couette flow, contrary to the observations made in some other three‐dimensional shear flows. Coherent flow structures, which are known to be associated with the Reynolds shear stress in near‐wall regions, were identified by the λ₂‐criterion. Instantaneous and ensemble‐averaged vortices resembled those found in the unidirectional Couette flow. In the skewed flow, however, the vortex structures were turned to align with the local mean‐flow direction. The conventional symmetry between Case 1 and Case 2 vortices was broken due to the mean‐flow three‐dimensionality. The turning of the coherent vortices and the accompanying symmetry‐breaking gave rise to secondary and tertiary turbulent shear stress components. By averaging the already ensemble‐averaged shear stresses associated with Case 1 and Case 2 vortices in the homogeneous directions, a direct link between the educed near‐wall structures and the Reynolds‐averaged turbulent stresses was established. These observations provide evidence in support of the hypothesis that the structural model proposed for two‐dimensional turbulent boundary layers remains valid also in flows with moderate mean three‐dimensionality. Copyright © 2009 John Wiley & Sons, Ltd.     

Motion of a circular cylinder in a viscous liquid between parallel plates

The two-dimensional motion of a cylinder in a viscous fluid between two parallel walls of a vertical channel is studied. It is found that when the cylinder moves very closely along one of the channel walls, it always rotates in the direction opposite to that of contact rolling along the nearest wall. When the cylinder is away from the walls, its rotation depends on the Reynolds number of the flow. In this study two numerical methods were used. One is for the unsteady motion of a sedimenting cylinder initially released from a position close to one of the channel walls, where the Navier-Stokes equations are solved for the fluid and Newton's equations of motion are solved for the rigid cylinder. The other method is for the steady flow in which a cylinder is fixed in a uniform flow field where the channel walls are sliding past the cylinder at the speed of the approaching flow, or equivalently a cylinder is moving with a constant velocity in a quiescent fluid. The flow field, the drag, the side force (lift), and the torque experienced by the cylinder are studied in detail. The effects of the cylinder location in the channel, the size of the channel relative to the cylinder diameter, and the Reynolds number of the flow are examined. In the limit when the cylinder is translating very closely along one of the walls, the flow in the gap between the cylinder and the wall is solved analytically using lubrication theory, and the numerical solution in the other region is used to piece together the whole flow field.This research was supported by NSF DMR91-20668 through the Laboratory for Research on the Structure of Matter at the University of Pennsylvania and from the Research Foundation of the University of Pennsylvania.     

Numerical simulation of turbulent flows in trapezoidal meandering compound open channels

Three‐dimensional (3D) numerical study is presented to investigate the turbulent flow in meandering compound open channels with trapezoidal cross‐sections. The flow simulation is carried out by solving the 3D Reynolds‐averaged continuity and Navier–Stokes equations with Reynolds stress equation model (RSM) for steady‐state flow. Finite volume method (FVM) is applied to numerically solve the governing equations of fluid flow. The velocity magnitude, tangential velocity, transverse velocity and Reynolds stresses are calculated for various flow conditions. Good agreement between the simulated and available laboratory measurements was obtained, indicating that the RSM can accurately predict the complicated flow phenomenon. Comparison of the calculated secondary currents of four cases (one being inbank flow and other three being overbank flow) with different water depths reveals that (i) the inbank flow exhibits different flow behaviors from that of the overbank flow does and (ii) the water depth has significant effects on the magnitude and direction of secondary currents. Copyright © 2010 John Wiley & Sons, Ltd.     

Finite element Galerkin Approach for a Computational Study of Arterial Flow

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｍｅａｓｕｒｅｍｅｎｔｏｆｓｈｅａｒｓｔｒｅｓｓｉｎａｒｔｅｒｉａｌｆｌｏｗｈａｓｂｅｅｎｂｅｓｔｏｗｅｄｂｙｎａｔｕｒｅｂｅｃａｕｓｅｏｆｔｈｅｐｏｓｓｉｂｌｅｒｅｌａｔｉｏｎｂｅｔｗｅｅｎｔｈｅｓｈｅａｒｓｔｒｅｓｓａｔｔｈｅｗａｌｌａｎｄｅｘｉｓｔｅｎｃｅｏｆａｔｈｅｒｏｓｃｌｅｒｏｓｉｓ.Ｔｈｅｆｌｏｗｐｈｅｎｏｍｅｎａｉｎｌａｒｇｅａｎｄｍｅｄｉｕｍｓｉｚｅｄａｒｔｅｒｉｅｓｉｎｆｌｕｅｎｃｅｔｈｅｄｅｖｅｌｏｐｍｅｎｔｏｆａｔｈｅｒｏｓｃｌｅｒｏｓｉｃｄ…     

Reynolds stress modelling of rectangular open‐channel flow

A Reynolds stress model for the numerical simulation of uniform 3D turbulent open‐channel flows is described. The finite volume method is used for the numerical solution of the flow equations and transport equations of the Reynolds stress components. The overall solution strategy is the SIMPLER algorithm, and the power‐law scheme is used to discretize the convection and diffusion terms in the governing equations. The developed model is applied to a flow at a Reynolds number of 77000 in a rectangular channel with a width to depth ratio of 2. The simulated mean flow and turbulence structures are compared with measured and computed data from the literature. The computed flow vectors in the plane normal to the streamwise direction show a small vortex, called inner secondary currents, located at the juncture of the sidewall and the free surface as well as the free surface and bottom vortices. This small vortex causes a significant increase in the wall shear stress in the vicinity of the free surface. A budget analysis of the streamwise vorticity is carried out. It is found that both production terms by anisotropy of Reynolds normal stress and by Reynolds shear stress contribute to the generation of secondary currents. Copyright © 2006 John Wiley & Sons, Ltd.     

Finite element computations of two-dimensional arterial flow in the presence of a transverse magnetic field

A finite element solution of the Navier-Stokes equations for steady flow under the magnetic effect through a double-branched two-dimensional section of a three-dimensional model of the canine aorta is discussed. The numerical scheme involves transforming the physical co-ordinates to a curvilinear boundary-fitted co-ordinate system. The shear stress at the wall is calculated for a Reynolds number of 1000 with the branch-to-main aortic flow rate ratio as a parameter. The results are compared with earlier works involving experimental data and found to be in reasonable qualitative agreement. The steady flow, shear stress and branch flow under the effect of a magnetic field have been discussed in detail.     

Spatial characterization of the numerically simulated vorticity fields of a flow in a flume

The topology of large scale structures in a turbulent boundary layer is investigated numerically. Spatial characteristics of the large scale structure are presented through an original method, proper orthogonal decomposition (POD) of the three-dimensional vorticity fields. The DNS results, obtained by Tiselj et al. [23] for a fully developed turbulent flow in a flume, are used in the present work to analyze coherent structures with the proposed methodology. In contrast to the reconstruction methods that use instantaneous flow quantities, this approach utilizes the whole dataset of the numerical simulation. The analysis uses one thousand 3D vorticity fields from 50000 time steps of the simulation for the Reynolds number of 2600 (the turbulent Reynolds number Re^*=171). The computational domain is 2146×171×537 wall units and the grid resolution is 128×65×72 points (in streamwise, wall-normal and spanwise directions, respectively). Experimental results obtained by using particle image velocimetry (PIV) in a fully developed turbulent boundary layer in a flume, which were analyzed with the same statistical characterization method, are in agreement with the DNS analysis: the dominant vortical structure appears to have a longitudinal streamwise orientation, an inclination angle of about 8°, streamwise length of several hundred wall units, and a distance between the neighboring structures of about 100 wall units in the spanwise direction. PACS 47.27.Nz, 47.54+r     

Vortex shedding from a square cylinder in presence of a moving wall

A numerical study on the flow past a square cylinder placed parallel to a wall, which is moving at the speed of the far field has been made. Flow has been investigated in the laminar Reynolds number (based on the cylinder length) range. We have studied the flow field for different values of the cylinder to wall separation length. The governing unsteady Navier–Stokes equations are discretized through the finite volume method on a staggered grid system. A SIMPLE type of algorithm has been used to compute the discretized equations iteratively. A shear layer of negative vortex generates along the surface of the wall, which influences the vortex shedding behind the cylinder. The flow‐field is distinct from the flow in presence of a stationary wall. An alternate vortex shedding occurs for all values of gap height in the unsteady regime of the flow. The strong positive vortex pushes the negative vortex upwards in the wake. The gap flow in the undersurface of the cylinder is strong and the velocity profile overshoots. The cylinder experiences a downward force for certain values of the Reynolds number and gap height. The drag and lift are higher at lower values of the Reynolds number. Copyright © 2005 John Wiley & Sons, Ltd.     

Vortex reynolds number in turbulent boundary layers

The effects of vortex Reynolds number on the statistics of turbulence in a turbulent boundary layer have been investigated. Vortex Reynolds number is defined as the ratio of circulation around the vortex structure to the fluid viscosity. The vortex structure of the outer region was modeled and a full numerical simulation was then conducted using a high-order spectral method. A unit domain of the outer region of a turbulent boundary layer was assumed to be composed of essentially three elements: a wall, a Blasius mean shear, and an elliptic vortex inclined at 45° to the flow direction. The laminar base-flow Reynolds number is roughly in the same range as that of a turbulent boundary layer based on eddy viscosity, and the vortex-core diameter based on the boundary-layer thickness is nearly the same as the maximum mixing length in a turbulent boundary layer. The computational box size, namely, 500, 150, and 250 wall units in the streamwise, surface-normal, and spanwise directions, respectively, is approximately the same as the measured quasi-periodic spacings of the near-wall turbulence-producing events in a turbulent boundary layer. The effects of vortex Reynolds number and the signs of the circulation on the moments of turbulence were examined. The signs mimic the ejection and sweep types of organized motions of a turbulent boundary layer. A vortex Reynolds number of 200 describes the turbulence moments in the outer layer reasonably well.     

Computer-aided analysis of flow past a surface-mounted obstacle

The incompressible, laminar, isothermal flow of a Newtonian fluid at steady state past a surface-mounted obstacle (flow over a step) is studied in a two-dimensional numerical experiment using the Galerkin finite element method. The dimensionless Navier–Stokes equations are solved in the whole range of the laminar flow regime. The numerical predictions are compared with available experimental data. The emphasis in the discussion of the results is on the presentation of the streamlines for various Reynolds numbers, the pressure distribution over and downstream of the step, the shear stress distribution along the surface of the step and the length of the recirculation region as a function of the Reynolds number. This analysis may be used in numerous applications from agricultural to civil, mechanical and chemical engineering. © 1997 John Wiley & Sons, Ltd.     

Computational technique for flow in blood vessels with porous effects

IntroductionThehemodynamicsofflowsthroughbloodvesselsisofgreatinterest,becausethesevesselspresentasubstantialhealthriskandareamajorcauseofmortalityandmorbidityintheindustrializedworld .Researchpapersonthebloodflowhaveappearedbutmostofthemhaveneglectedtheporosityeffectsduetovesselwalls.Inthisstudyweareinterestedintheflowthroughabloodvesseltakingintoaccounttheporosityeffectsofthevessels.Fluidflowthroughaporousmediumisoffundamentalimportancetowiderangeofdisciplinesinthevariousbranchesofnaturalsci…     

Modeling of convective structures in a thin layer of silicone oil in air flow under heating

Numerical simulation was performed to study convective structures in a thin silicone oil layer heated from below, whose free surface is exposed to air flow generating drift flow. The basic equations are transformation to a form suitable for spectral methods. The steady flow velocity profile obtained in a laboratory experiment is calculated. It is shown that increasing the Reynolds number leads to the transition from polygonal convective cells to longitudinal rolls (elongated along the flow). The dependence of the transition Reynolds number on the temperature on the lower boundary of the layer is obtained. The calculation results are compared with experimental data.     

Analysis of Linear Viscous Flow with a Free Surface in a Circular Cylinder Subjected to Small-Amplitude Circular Oscillations

A fluid flow with a free surface inside a circular cylinder subjected to horizontal, circular oscillation was analyzed theoretically and numerically under the assumption of small-amplitude oscillation and high Reynolds number. It was shown that the nature of the oscillatory flow is of a standing-wave type when projected onto the axial plane and of a progressive-wave type when projected onto the azimuthal plane. The Stokes drift motion in the azimuthal direction and the steady streaming velocity at the edge of the bottom- and side-wall boundary layers are then used in the numerical computation for the steady axisymmetric recirculatory flow outside the boundary layers. We have found that the solutions can be well predicted by asymptotic analysis for the full Navier–Stokes equations in the low streaming-Reynolds-number limit. A simple experiment on flow visualization revealed a good agreement in the surface flow pattern on the bottom wall. It also provided steady recirculatory flows that were not much different from the numerical results on the whole. Received 16 August 1999 and accepted 14 January 2000     

Unsteady laminar flow developing in a curved duct

The flow developing in a tightly curved U-bend of square cross section has been investigated experimentally and via numerical simulation. Both long-time averages and time histories of the longitudinal (streamwise) component of velocity were measured using a laser-Doppler velocimeter. The Reynolds number investigated was Re = 1400. The data were obtained at different bend angles, θ, and were confined to the symmetry plane of the bend. At Re = 1400, the flow entering the bend is steady, but by θ = 90° it develops an oscillatory component of motion along the outer-radius wall. Autocorrelations and energy spectra derived from the time histories yield a base frequency of approximately 0.1 Hz for these oscillations. Flow-visualization studies showed that the proximity of the outer-radius wall served to damp the amplitude of the spanwise oscillations.

Numerical simulations of the flow were performed using both steady and unsteady version of the finite-difference elliptic calculation procedure of Humphrey et al. (1977). Although the unsteadiness observed experimentally does not arise spontaneously in the calculations, numerical experiments involving the imposition of a periodic time-dependent perturbation at the inlet plane suggest that the U-bend acts upon the incoming flow so as to damp the amplitude of the imposed oscillation while altering its frequency.

The oscillations observed experimentally, and numerically as a result of the periodic perturbation, have been linked to the formation of Goertler-type vortices of the outer-radius wall in the developing flow. The vortices, which develop as a result of the centrifugal instability of the flow on the outer-radius wall, undergo a further transition to an unsteady regime at higher flow rates.     

The Tilting Mechanism of a Longitudinal Vortical Structure in a Homogeneous Shear Flow with and Without Spanwise Rotation

A longitudinal vortical structure is typically observed in near-wall turbulence. This vortical structure is elongated in the streamwise direction, though it is also tilted in the spanwise direction. The sense of this spanwise tilting is determined by the sign of the streamwise vorticity associated with the vortex, and longitudinal vortical structures with a different streamwise vorticity become asymmetric (mirror symmetric). The tilting must be due to the combined effects of the non-linear terms and mean spanwise vorticity associated with the mean shear. However, the detailed mechanism of the tilting is not well known. To study the tilting in detail, we performed direct numerical simulations of a homogeneous shear flow where the longitudinal vortical structures similar to those in the near-wall region are observed. In particular, the effects of spanwise system rotation as well as the Reynolds number on the vortical structure are studied. As a result, we found that spanwise system rotation has more marked effects on the vortical structure than the Reynolds number. When the system rotation is imposed in the same direction as the mean spanwise vorticity, the tilting is enhanced, while the system rotation of the opposite direction attenuates it. We also found that when the longitudinal vortical structure is tilted in the spanwise direction, it is sandwiched between the streamwise vorticity of the opposite sign. The cyclonic rotation enhances the streamwise vorticity of the opposite sign, though the longitudinal vortical structure at the center is attenuated. In the anticyclonic case, the streamwise vorticity of the opposite sign almost disappears and the longitudinal vortical structure is isolated from the surrounding flow.