Turbulent flow in converging nozzles,part one: boundary layer solution

The boundary layer integral method is used to investigate the development of the turbulent swirling flow at the entrance region of a conical nozzle. The governing equations in the spherical coordinate system are simplified with the boundary layer assumptions and integrated through the boundary layer. The resulting sets of differential equations are then solved by the fourth-order Adams predictor-corrector method. The free vortex and uniform velocity profiles are applied for the tangential and axial velocities at the inlet region, respectively. Due to the lack of experimental data for swirling flows in converging nozzles, the developed model is validated against the numerical simulations. The results of numerical simulations demonstrate the capability of the analytical model in predicting boundary layer parameters such as the boundary layer growth, the shear rate, the boundary layer thickness, and the swirl intensity decay rate for different cone angles. The proposed method introduces a simple and robust procedure to investigate the boundary layer parameters inside the converging geometries.     

高超声速层流尾迹的数值模拟

本文利用无波动、无自由参数、耗散的差分格式(NND格式),通过求解NS方程,数值模拟了高超声速层流尾迹的流动,清晰地给出了主激波、拐角膨胀波、迹激波及自由剪切层,所得流场物理量的分布与实验结果甚为一致。计算发现了底部迴流区由起始向定常的发展中,在瞬时流线图上经历了极限环形成、胀大、缩小、再胀大最后消失的演变过程。     

Flow field characterization of a jet and vortex actuator

 An actuator, which produces several different flow fields that may be used for active flow control, is characterized in still air using flow visualization and velocity measurements. The primary actuator-induced flow fields are: free jets, wall jets, and vortex flows. The non-dimensional parameters governing these actuator-induced flows are developed. For the vortex-flow regime, the operational range of the actuator increases as it’s size decreases without a significant decrease in either the actuator induced velocity or vortex core size. The velocity scaling is developed for the vortex flow and suggests that the optimum actuator efficiency occurs at a Stokes number of approximately 7.9 for the range of parameters surveyed. In a turbulent, zero pressure gradient boundary layer, measurements made just downstream of the actuator (when operated in the vortex mode) indicate a vortical disturbance is generated in the boundary layer. Received: 2 September 1998/Accepted: 9 January 1999     

A bypass wake induced laminar/turbulent transition

The process of laminar to turbulent transition induced by a von Karman vortex street wake, was studied for the case of a flat plate boundary layer. The boundary layer developed under zero pressure gradient conditions. The vortex street was generated by a cylinder positioned in the free stream. An X-type hot-wire probe located in the boundary layer, measured the streamwise and normal to the wall velocity components. The measurements covered two areas; the region of transition onset and development and the region where the wake and the boundary layer merged producing a turbulent flow. The evolution of Reynolds stresses and rms-values of velocity fluctuations along the transition region are presented and discussed. From the profiles of the Reynolds stress and the mean velocity profile, a ‘negative' energy production region along the transition region, was identified. A quadrant splitting analysis was applied to the instantaneous Reynolds stress signals. The contributions of the elementary coherent structures to the total Reynolds stress were evaluated, for several x-positions of the near wall region. Distinct regions in the streamwise and normal to the wall directions were identified during the transition.     

Specified discharge velocity models for numerical simulations of laminar vortex rings

We numerically and theoretically investigate the flow generated at the exit section of a piston/cylinder arrangement that is generally used in experiments to produce vortex rings. Accurate models for the velocity profile in this section (also called specified discharge velocity, SDV models) are necessary in (i) numerical simulations of laminar vortex rings that do not compute the flow inside the cylinder and (ii) in slug-models that provide a formula for the total circulation of the flow. Based on the theoretical and numerical analysis of the flow evolution in the entrance region of a pipe, we derive two new and easy to implement SDV models. A first model takes into account the unsteady evolution of the centerline velocity, while the second model also includes the time variation of the characteristics of the boundary layer at the exit plane of the vortex generator. The models are tested in axisymmetric direct numerical simulations of vortex rings. As distinguished from classical SDV model, the new models allow to accurately reproduce the characteristics of the flow. In particular, the time evolution of the total circulation is in good agreement with experimental results and previous numerical simulations including the vortex generator. The second model also provides a more realistic time evolution of the vortex ring circulation. Using the classical slug-model and the new correction for the centerline velocity, we finally derive a new and accurate analytical expression for the total circulation of the flow.     

Simulation of Entry-Region Flow in Open-Cell Metal Foam and Experimental Validation

Open-cell metal foam is distinguished from traditional porous media by its very high porosities (often greater than 90 %), and its web-like open structure and good permeability. As such, the foam is a very attractive core for many engineered systems, e.g., heat exchangers, filtration devices, catalysts, and reactors. The flow field inside the foam is rather complex due to flow reversal and vigorous mixing. This complexity is increased by the possible presence of an entry region. The entrance region in metal foam is usually underestimated and ignored, just like its counterpart in traditional porous media. In this paper, the actual entry length is determined by simulation and direct experiment on commercial open-cell aluminum foam. It is shown to be dependent on flow velocity and to reach a constant value for higher velocities. The complex and intrinsically random architecture of the foam is idealized using a unit geometrical model, in order to numerically investigate the flow field and pressure drop inside the foam. The Navier–Stokes equations are solved directly, and velocity and pressure fields are obtained for various approach velocities using a commercial numerical package. The entry length is ascertained from the behavior of the velocity field close to the entrance. Comparisons to experimental data were also carried out. The commercial foam that was used in the experiment had 10 ppi and porosity of 91.2 %. Air was forced to flow inside the foam using an open-loop wind tunnel. Good qualitative agreement between the modeling and experimental results are obtained. The agreement lends confidence to the modeling approach and the determined entry length.     

Numerical investigation of a spatially developing turbulent natural convection boundary layer along a vertical heated plate

Large-eddy simulation (LES) on a spatially developing natural convection boundary layer along a vertical heated plate was conducted. The heat transfer rate, friction velocity, mean velocity and temperature, and second-order turbulent properties both in the wall-normal and the stream-wise direction showed reasonable agreement with the findings of past experiments. The spectrum of velocity and temperature fluctuation showed a -2/3-power decay slope and -2-power decay slope respectively. Quadrant analysis revealed the inclination on Q1 and Q3 in the Reynolds stress and turbulent heat flux, changing their contribution along the distance from the plate surface. Following the convention, we defined the threshold region where the stream-wise mean velocity takes local maximum, the inner layer which is closer to the plate than the threshold region, the outer layer which is farther to the plate than the threshold region. The space correlation of stream-wise velocity tilted the head toward the wall in the propagating direction in the outer layer; on the other hand, the correlated motion had little inclination in the threshold region. The time history of the second invariant of gradient tensor Q revealed that the vortex strength oscillates both in the inner and the outer layers in between the laminar and the transition region. In the turbulent region, the vortex was often dominant in the outer layer. Instantaneous three-dimensional visualization of Q revealed the existence of high-speed fluid parcels associated with arch-shape vortices. These results were considered as an intrinsic structure in the outer layer, which is symmetrical to the structure of canonical smooth/rough wall bounded layer flow in forced convection.     

A transient natural convection of micropolar fluids over a vertical cylinder

A transient free convective boundary layer flow of micropolar fluids past a semi-infinite cylinder is analysed in the present study. The transformed dimensionless governing equations for the flow, microrotation and heat transfer are solved by using the implicit scheme. For the validation of the current numerical method heat transfer results for a Newtonian fluid case where the vortex viscosity is zero are compared with those available in the existing literature, and an excellent agreement is obtained. The obtained results concerning velocity, microrotation and temperature across the boundary layer are illustrated graphically for different values of various parameters and the dependence of the flow and temperature fields on these parameters is discussed. An increase in the vortex viscosity tends to increase the magnitude of microrotation and thus decreases the peak velocity of fluid flow. An increase in the vortex viscosity in micropolar fluids is shown to decrease the heat transfer rate.     

Velocity distributions in a hydrocyclone separator

 The internal three-dimensional flow field in a hydrocyclone was studied using laser velocimetry. Seven axial planes were investigated for three different inlet flow rates and three independent and different rejects rates. Results at each measurement plane showed that the measured tangential velocity profile behaves like a forced vortex at the region near the air core, and like a free vortex in the outer portion of the flow. The peak nondimensional tangential velocity decreases as the distance from the inlet region increases, however, the peak dimensional tangential velocity increases as the distance from the inlet region increases. The nondimensional peak tangential velocities are approximately equal for all of the flow rates. The magnitude of the tangential velocity increased in the inner forced vortex region as the rejects rate was increased. Backflows exist in the axial velocity profile near the inlet region, but these reversed flows disappear in the exit region. The dimensional vorticity is proportional to inlet flow rate and decreases with increasing rejects flow rates. Received: 27 February 2001/Accepted: 19 June 2001     

Asymptotic model of the development of separations inside a boundary layer under the action of a traveling pressure wave

A model of the nonlinear interaction between a pressure perturbation traveling at a constant velocity and an incompressible boundary layer is constructed when its near-wall part is described by the “inviscid boundary layer” equations. A steady-state solution is managed to obtain in the finite form under the assumption that it exists in a moving coordinate system. It is shown that the boundary layer can easily overcome pressure perturbations whose amplitude is not higher than the dynamic pressure calculated from the velocity of the pressure perturbation. At the higher pressure perturbation amplitudes a vortex sheet sheds from the body surface to the boundary layer. The vortex sheet represents an unstable surface of the tangential discontinuity which separates the regions of the direct and reverse separation flows. In the case of an arbitrary shape of the pressure perturbation the surface of the tangential discontinuity sheds from the body surface at a finite angle with the formation of a stagnation point. An example of the pressure perturbation in which the vortex sheet sheds from the body surface along the tangent is constructed.     

Laminar boundary layer flow of a nanofluid along a wedge in the presence of suction/injection

The behavior of an incompressible laminar boundary layer flow over a wedge in a nanofluid with suction or injection has been investigated. The model used for the nanofluid integrates the effects of the Brownian motion and thermophoresis parameters. The governing partial differential equations of this problem, subjected to their boundary conditions, are solved by the Runge-Kutta-Gill technique with the shooting method for finding the skin friction and the rate of heat and mass transfer. The result are presented in the form of velocity, temperature, and volume fraction profiles for different values of the suction/injection parameter, Brownian motion parameter, thermophoresis parameter, pressure gradient parameter, Prandtl number, and Lewis number. The conclusion is drawn that these parameters significantly affect the temperature and volume fraction profiles, but their influence on the velocity profile is comparatively smaller.     

Some Observations of Vortex Breakdown in a Confined Flow with Solid Body Rotation

The occurrence of breakdown in slender vortex flows as a ``bubble'' or ``spiral'' pattern depends on the degree of radial deflection of the vortex core from its original axis as shown in [1]. A smooth transition from a bubble to a spiral-type ``mode'' can be forced by inducing a small asymmetric disturbance which led to the conclusion, that the patterns do not represent different fundamental modes of breakdown. The subject presented herein addresses the following question: how does breakdown evolve in a swirling flow in which the vortex core is forced on a straight axis? In addition, what is the effect of turbulent inflow conditions? This type of vortex conditions is achieved in a spinning tube flow. The swirl is introduced at the entrance of the rotating tube with a honeycomb package and maintained by the viscous action in the boundary layer of the spinning tube. A diffuser at the end induces an adverse pressure gradient to force the breakdown. Flow visualization experiments are carried out to characterize the nature of breakdown over a range of different flow conditions. For some selected characteristic stages, detailed velocity fields were obtained using the method of Digital Particle-Image-Velocimetry (DPIV). The results show, that for the range of parameters investigated, breakdown is initiated at Rossby-numbers below a critical value of Ro ≈ 0.6 similar to those observed in other experiments. The bursted part of the vortex has a near axi-symmetric slender conical shape containing approximately stagnant flow. Its downstream end is characterized by a jump-like contraction where the flow evolves into a jet with enhanced swirl on the axis. It is only in this region downstream of the jump-like contraction that asymmetric instabilities and wavy flow patterns could be observed. Perturbations caused by them travel upstream but do not change the near-axisymmetric shape of the bursted part of the vortex.     

Experimental investigation of the unsteady structure of a transitional plane wall jet

 A laminar wall jet undergoing transition is investigated using the particle image velocimetry (PIV) technique. The plane wall jet is issued from a rectangular channel, with the jet-exit velocity profile being parabolic. The Reynolds number, based on the exit mean velocity and the channel width, is 1450. To aid the understanding of the global flow features, laser-sheet/smoke flow visualizations are performed along streamwise, spanwise, and cross-stream directions. Surface pressure measurements are made to correlate the instantaneous vorticity distribution with the surface pressure fluctuations. The instantaneous velocity and vorticity field measurements provide the basis for understanding the formation of the inner-region vortex and the subsequent interactions between the outer-region (free-shear-layer region) and inner-region (boundary-layer region) vortical structures. Results show that under the influence of the free-shear-layer vortex, the local boundary layer becomes detached from the surface and inviscidly unstable, and a vortex is formed in the inner region. Once this vortex has formed, the free-shear-layer vortex and the inner-region vortex form a vortex couple and convect downstream. The mutual interactions between these inner- and outer-region vortical structures dominate the transition process. Farther downstream, the emergence of the three-dimensional structure in the free shear layer initiates complete breakdown of the flow. Received: 8 November 1995/Accepted: 6 November 1996     

Direct Numerical Simulation of Self-Similar Turbulent Boundary Layers in Adverse Pressure Gradients

Direct numerical simulations of the Navier–Stokes equations have been carried out with the objective of studying turbulent boundary layers in adverse pressure gradients. The boundary layer flows concerned are of the equilibrium type which makes the analysis simpler and the results can be compared with earlier experiments and simulations. This type of turbulent boundary layers also permits an analysis of the equation of motion to predict separation. The linear analysis based on the assumption of asymptotically high Reynolds number gives results that are not applicable to finite Reynolds number flows. A different non-linear approach is presented to obtain a useful relation between the freestream variation and other mean flow parameters. Comparison of turbulent statistics from the zero pressure gradient case and two adverse pressure gradient cases shows the development of an outer peak in the turbulent energy in agreement with experiment. The turbulent flows have also been investigated using a differential Reynolds stress model. Profiles for velocity and turbulence quantities obtained from the direct numerical simulations were used as initial data. The initial transients in the model predictions vanished rapidly. The model predictions are compared with the direct simulations and low Reynolds number effects are investigated.     

Radiation effect on mixed convection along a vertical plate with uniform surface temperature

This paper investigates the effect of radiation on the forced and free convection flow of an optically dense viscous incompressible fluid along a heated vertical flat plate with uniform free stream and uniform surface temperature with Rosseland diffusion approximation. With appropriate transformations, the boundary layer equations governing the flow are reduced to local nonsimilarity equations valid in the forced convection regime as well as in the free convection regime. A group of transformation is, also, introduced to reduce the boundary layer equations to a set of local nonsimilarity equations valid in both the forced and free convection regimes. Solutions of the governing equations are obtained by employing the implicit finite difference methods together with Keller box scheme and are expressed in terms of local shear stress and local rate of heat transfer for a range of values of the pertinent parameters.     

Evolution of an initially columnar vortex terminating normal to a no-slip wall

The early evolution of an initially columnar vortex normal to a solid wall was examined. The vortex was generated by a pair of flaps in a water tank. Detrimental effects from the wall during the vortex generation were avoided by producing the vortex normal to a free surface and subsequently bringing a horizontal plate into contact with the surface. Digital particle image velocimetry (DPIV) measurements of the velocity and vorticity, together with laser induced fluorescence (LIF) visualizations, in a meridional plane revealed a toroidal structure with the appearance of an axisymmetric vortex breakdown bubble. Agreement was found between the measurements and numerical simulations of the axisymmetric Navier–Stokes equations. The results show that the flow in the effusive corner region is dominated by a Bödewadt-type spatially oscillatory boundary layer within the core region and a potential-like vortex boundary layer at large radii. The toroidal structure results from the interaction between these two boundary layers, leading to the roll up of a dominant shear layer within the Bödewadt structure, and does not develop from the columnar vortex itself.     

Experimental and numerical study of the stability of a pre-separated flow on an airfoil

The problem of the origin and evolution of two-dimensional waves of unstable disturbances in the boundary layer on an airfoil in the region of adverse pressure gradient in the preseparation flow region is solved numerically. The stability of the experimental velocity profiles, including the inflected profiles, is studied. As a result of the calculations, the boundaries of the instability region and the parameters of the maximally unstable disturbances (frequency, growth rate, wavelength, and propagation velocity) are determined for each velocity profile. The characteristics obtained in the present work are in good agreement with the real experimental parameters of instability waves. Institute of Theoretical and Applied Mechanics, Siberian Division, Russian Academy of Sciences, Novosibirsk 630090. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 40, No. 1, pp. 126–132, January–February, 1999.     

The axisymmetric boundary layer beneath a Rankine-like vortex

 An experimental investigation of the three-dimensional boundary layer induced by a Rankine-like vortex with its axis normal to a stationary disk is described. The velocity field through the boundary layer was measured for Reynolds number Re (based on the tangential velocity and radius at the disk edge) ranging from 10 000 to 25 000 at various radial distances by means of a 4-beam, 2-component Laser Doppler Anemometer. Our results show that the nature of the boundary layer is affected by two factors: an inflexional instability caused by the crossflow velocity profile and a stability factor caused by the favorable pressure gradient. At lower Reynolds number, the radial pressure gradient has a very strong stabilizing effect on the boundary layer and acts to revert it to its laminar state upstream of the effusing core. At higher Re the inflexional instability caused by the crossflow velocity dominates while the stabilizing influence of the favorable pressure gradient recedes. As such, laminar reversion likely occurs closer to the effusion core. Thus, the point of laminar reversion moves closer to the effusion core as the Reynolds number is increased. Received 23 May 1996 / Accepted 29 July 1996