Formation of a tripolar vortex in a stratified fluid

This paper reports on an experimental study of vortices in a stratified fluid. The vortices were generated by two different stirring devices, viz. a rotating sphere and a rotating bent rod. It was found that the vortices created with the rotating sphere are mostly axisymmetric and stable, whereas the vortices produced with the bent rod generally show instabilities, under certain conditions leading to the formation of a tripolar vortex. This report concentrates on this tripolar structure and presents quantitative information about the flow obtained through streak photography of tracer particles.     

Liquid and gas entrainment to a small break hole from a stratified two-phase region

A theoretical and experimental study has been conducted for liquid and gas entrainment to a small break hole from a stratified two-phase region. Theoretical correlations previously obtained for top, side and bottom entrainment were modified to express the relation between break flow rate, break quality and bulk water level so that they can be used easily for any break geometry. The modified correlations were assessed with experimental data obtained under room temperature and low pressure conditions using air and water. The experiment results were predicted well with the present model without using any adjustment coefficient when the upstream flow was symmetrical around the break. The effects of vortex, crossflow and wavy flow, observed in the experiment but not taken into account in the model, were empirically correlated based on the present correlation. By using the empirically modified correlations, the data in the literature, including high-pressure steam-water conditions, were reasonably predicted.     

Dust entrainment by means of a planar shock induced vortex over loose dust layers

Conclusions The recent study of Ben-Dor and Rayevsky (1994) regarding the interaction of planar shock waves with high density layers was reconsidered in order to demonstrate a possible dust entrainment mechanism which has not received appropriate attention so far. It was shown that, as a result of the interaction, a large vortex is generated. This vortex could in fact contribute to the entrainment of dust when planar shock waves interact with loose dusty layers.In addition, the effect of viscosity on the proposed dust entrainment mechanism was also considered.     

Numerical simulation of quasi-streamwise hairpin-like vortex generation in turbulent boundary layer

A mechanism for generation of near wall quasi-streamwise hairpin-like vortex （QHV） and secondary quasi-streamwise vortices （SQV） is presented. The conceptual model of resonant triad in the theory of hydrodynamic instability and direct numerical simulation of a turbulent boundary layer were applied to reveal the formation of QHV and SQV. The generation procedures and the characteristics of the vortex structures are obtained, which share some similarities with previous numerical simulations. The research using resonant triad conceptual model and numerical simulation provides a possibility for investigating and controling the vortex structures, which play a dominant role in the evolution of coherent structures in the near-wall region.     

Numerical calculations of three-dimensional turbulent vortex breakdown

Numerical solutions to the three-dimensional, unsteady, incompressible Reynolds-averaged Navier-Stokes equations have been obtained for bubble-type vortex breakdown. Two different turbulence models were employed: (1) standard K-ε and (2) an explicit, regularized algebraic Reynolds stress model. Results are computed at a Reynolds number of 10,000. The algebraic Reynolds stress model produced a breakdown bubble with a larger length-to-diameter ratio than did the K-ε model. Breakdown also occurred at lower levels of adverse pressure gradient for the algebraic stress model than for the K-ε model. In each case single-cell breakdown structures resulted. This is contrasted with numerical calculations for laminar breakdown which reveal the existence of complex multicell bubble breakdown structures.     

Simulating air entrainment and vortex dynamics in a hydraulic jump

The air entrainment characteristics of three separate Froude number hydraulic jumps are investigated numerically using an unsteady RANS, realizable k–ε turbulence model, with a Volume of Fluid treatment for the free surface. Mean velocity profiles, average void fraction, and Sauter mean diameter compare favorably with experimental data reported in literature. In all simulations, time-averaged void fraction profiles show good agreement with experimental values in the turbulent shear layer and an accurate representation of interfacial aeration at the free surface. Sauter mean diameter is well represented in the shear layer, and free surface entrainment results indicate bubble size remains relatively unchanged throughout the depth of the jump. Several different grid resolutions are tested in the simulations. Significant improvements in void fraction and bubble size comparison are seen when the diameter to grid size ratio of the largest bubbles in the shear layer surpasses eight. A three-dimensional simulation is carried out for one Froude number jump, showing an improvement in the prediction of entrained air and bubble size compared with two-dimensional results at a substantial increase in computation time. An analysis of three-dimensional vorticity shows a complex interaction between spanwise and streamwise vortical structures and entrained air bubbles. The jump is similar to a turbulent mixing layer, constrained by the free surface, with vortex pairing and subsequent fluctuations in free surface elevation. Downstream fluctuations of the toe are associated with a roll up of the primary spanwise vortex, fluctuations of the free surface, and counter-rotating streamwise vortex pairs. The action of these flow structures is likely responsible for the improvement in three-dimensional results.     

Hybrid turbulence simulation of spray impingement cooling: The effect of vortex motion on turbulent heat flux

Spray impingement has been a major interest of researchers in the areas of spray cooling, internal combustion, fire suppression and spray cooling, etc. for a long time. Numerous studies have been done in the area of spray cooling. Spray cooling with phase change takes advantage of relatively large amounts of latent heat and is capable of removing high heat fluxes from the surface, which has generated the interest of many researchers. In this paper, the turbulent characteristics of vapor formed during the spray impingement are studied. Water and gasoline are used in the numerical analysis of the two‐phase spray impingement on a heated wall. Hybrid turbulence modeling was used for the analysis where the subgrid scale model was employed away from the wall and k–ε model was used near the wall. Gasoline, at 298 K, was sprayed on the heated wall, kept constant at 650 K. The surrounding temperature was maintained at 400 K at the start of the simulation. In case of water and gasoline at Reynolds number 2750, the heated wall was kept constant at 400 K and the surrounding temperature was maintained at 298 K at the start of the simulations. The nozzle diameter of 100µm was used for this study, with the nozzle plate spacing ratio at 10. The spray was impinged on the flat plate at angles of 0, 15, and 30^°. Root mean‐squared velocities and turbulent heat flux were plotted in the water spray impingement for the different angles of impingement. The effect of turbulence on the heat transfer was observed. The effect of vortex motion on the turbulent heat flux values was analyzed using different Reynolds numbers of impingement and at different angles in case of gasoline. The turbulent heat flux attained the maximum values with high vortex formation. Upwash of fluid transported heat away from the wall, producing higher heat flux values in the region. Copyright © 2008 John Wiley & Sons, Ltd.     

Simulation of the Effect of the Freestream Turbulence Parameters on Heat Transfer in an Unsteady Boundary Layer

A variant of the two-parameter turbulence model which makes it possible continuously to calculate a flow region with laminar, transition and turbulent regimes is proposed for investigating the flow under conditions of high freestream turbulence intensity. It is shown that the properties of the thermal transition can be theoretically described using the quasi-steady turbulence model in the case of periodic freestream velocity distribution. The numerical results are compared with theoretical and experimental data. The approach proposed is developed for determining the combined effect of the parameters of harmonic fluctuations of the external velocity and freestream turbulence on the heat transfer characteristics on a flat plate with different boundary conditions for the enthalpy.     

Receptivity of the Blunt-Leading-Edge Plate Boundary Layer to Unsteady Vortex Perturbations

The response of the boundary layer on a plate with a blunt leading edge to frozen-in vortex perturbations whose vorticity is normal to the plate surface is found. It is shown that these vortices generate an inhomogeneity of the streamwise velocity component in the boundary layer. This inhomogeneity is analogous to the streaky structure developing as the degree of free-stream turbulence increases. The dependence of the amplitude and shape of the boundary layer inhomogeneity on the distance from the leading edge, the streamwise and spanwise scales, and other parameters is found for periodic and local initial perturbations. It is shown that the receptivity of the boundary layer decreases with increase in the frequency and with decrease in the streamwise perturbation scale.     

A short note on the entrainment and exit boundary conditions

An attempt is made to find out the suitable entrainment and exit boundary conditions in laminar flow situations. Streamfunction vorticity formulation of the Navier–Stokes equations are solved by ADI method. Two‐dimensional laminar plane wall jet flow is used to test different forms of the boundary conditions. Results are compared with the experimental and similarity solution and the proper boundary condition is suggested. The Kind 1 boundary condition is recommended. It consists of zero first derivative condition for velocity variable and for streamfunction equation, mixed derivative at the entrainment and exit boundaries. Copyright © 2005 John Wiley & Sons, Ltd.     

Transformation of hanging discontinuities into vortex systems in a stratified flow behind a cylinder

The mechanisms of formation of hanging discontinuities, vortex dipoles, and vortex arrays in the wave wake behind a cylinder moving at a constant velocity in a stratified fluid are investigated using various schlieren methods. The existence of discontinuities is attributable to the distortion of the internal-wave phase pattern in the shear flow and to the varying stratification and subsequent interaction of the waves with the appearing nonuniformities. Hanging discontinuities and vortex systems are low-velocity analogs of shock waves. An analysis of the internal-wave pattern indicates that the values of the normal velocity component differ on the upper and lower edges of the discontinuities. A regime diagram for flows of this kind is given.     

Computations of Interaction of an Incident Oblique Shock Wave with a Turbulent Boundary Layer on a Flat Plate

Results of numerical simulation of interaction between an oblique shock wave and a turbulent boundary layer formed in a supersonic (Mach number M =5) flow past a flat plate are presented. The computations are performed for three cases of interaction of different intensity, which result in an attached or detached flow. Numerical results are compared with experimental data. The effect of flow turbulence and shockwave unsteadiness on flow parameters is studied.     

A model for droplet entrainment in heated annular flow

The ability to accurately predict droplet entrainment in annular two-phase flow is required to effectively calculate the interfacial mass, momentum, and energy transfer, which characterizes nuclear reactor safety, system design, analysis, and performance. Most annular flow entrainment models in the open literature are formulated in terms of dimensionless groups, which do not directly account for interfacial instabilities. However, many researchers agree that there is a clear presence of interfacial instability phenomena having a direct impact on droplet entrainment. The present study proposes a model for droplet entrainment, based on the underlying physics of droplet entrainment from upward co-current annular film flow that is characteristic to light water reactor safety analysis. The model is developed based on a force balance and stability analysis that can be implemented into a transient three-field (continuous liquid, droplet, and vapor) two-phase heat transfer and fluid flow systems analysis computer code.     

Receptivity of the inclined-cylinder attachment-line boundary layer to vortex perturbations

The receptivity of the boundary layer in the neighborhood of the attachment line of a cylinder inclined to the flow with respect to periodic vortex perturbations frozen into the stream is investigated. The problem considered simulates the interaction between external turbulence and the leading-edge swept wing boundary layer. It is shown that if the direction of the external perturbation vector is almost parallel to the leading edge, then the external perturbations are considerably strengthened at the outer boundary layer edge. This effect can cause laminar-turbulent transition on the attachment line at subcritical Reynolds numbers.Translated from Izvestiya Rossiiskoi Academii Nauk, Mekhanika Zhidkosti i Gaza, No. 6, 2004, pp. 72–85. Original Russian Text Copyright © 2004 by Ustinov.     

Computation of shock wave/turbulent boundary layer interactions using a two-equation model with compressibility corrections

In this paper computational results for two different types of shock wave / turbulent boundary layer interaction flows are presented. It is shown that upstream effects of the shock induced separation cannot be reproduced by Wilcox's (1991) k--model, whereas downstream of the interaction, predictions of pressure distribution and skin friction are acceptable. The inclusion of the compressible part of the dissipation rate and the pressure dilatation in the model has noticeable, but not dramatic effects on wall pressure and skin friction in the selected flow cases.     

Large eddy simulation of turbulent statistical and transport properties in stably stratified flows

Three dimensional large eddy simulation （LES） is performed in the investigation of stably stratified turbulence with a sharp thermal interface. Main results are focused on the turbulent characteristic scale, statistical properties, transport properties, and temporal and spatial evolution of the scalar field. Results show that the buoyancy scale increases first, and then goes to a certain constant value. The stronger the mean shear, the larger the buoyancy scale. The overturning scale increases with the flow, and the mean shear improves the overturning scale. The flatness factor of temperature departs from the Gaussian distribution in a fairly large region, and its statistical properties are clearly different from those of the velocity fluctuations in strong stratified cases. Turbulent mixing starts from small scale motions, and then extends to large scale motions.     

湍流的相干结构:推演方法和结果

本文讨论壁面湍流发展的相干结构的观点.在简要的历史文献综述后,我们回顾一些基本观点, 并且介绍相干结构的思想.基于大量主要是由实验所得的结果,本文通过广泛运用的事件检 测技术,探讨湍流边界层内部和外部区域发生的现象.我们从边界层内部区域发生的现象、 边界层外部大尺度运动的发展和涡结构动力学的角度来描述流动的现象.在文章的 第2部分,介绍从背景流动中推演出湍流相干结构的各种方法以及在各种方法框架下所得到的结果, 讨论速度梯度张量不变量、压力的Hessian矩阵分析和本征正交分解等方法.每一个过程 都有``相干结构'的特定的定义,满足恰当的数学构架,并可以对湍流数据做相干结构动力 学分析.这一工作可能会对当前流体动力学家在湍流研究中用到的最新理论和技术的传播有 所贡献.     

用格子Boltzmann方法数值模拟混合层中旋涡的合并过程

用格子Boltzmann方法计算混合层中的流动问题。在流场的入口处加不同频率、振幅和相位的小扰动,观察混合层中旋涡的演进机理,模拟二维混合层中旋涡合并现象。在基本扰动波的基础上,又加入频率为基本波频率一半的亚谐波,得到了两个涡合并的计算结果,当加入的亚谐波频率为基本波频率的三分之一时,得到了三个涡合并的计算结果。这些计算结果与已有文献的结果基本一致,显示用格子Boltzmann方法模拟混合层问题是可行的。     

DNS OF 2D TURBULENT FLOW AROUND A SQUARE CYLINDER

Two-dimensional ‘turbulent’ flow around a rectangular cylinder has been simulated at Re=10,000 using a sixth- order-accurate finite volume method for the discretization of convection and diffusion. The spatial discretization consists of a combination of a seventh- order upwind-biased method for the convective terms and an eighth-order central method for the diffusive terms, discretized on a stretched and staggered grid. To cope with the stretching of the grid, Lagrange interpolations are used. The method applied to obtain a boundary condition for the velocity in the x-direction at the outflow boundary is shown not to affect the flow in the interior of the computational domain in a way that is visible in various snapshots of the vorticity field. The variation in the velocity in the x-direction with time is itself found to be relatively small near the outflow boundary. Several turbulence statistics have been gathered from a simulation of the flow developed during 77 dimensionless time units. Snapshots of the vorticity field of the developed flow show the presence of a vortex-street- like structure. Typical 2D turbulent behaviour, such as the appearance of monopolar, dipolar and tripolar vortices due to the amalgamation of vorticity in the wake and the x^−1/2 scaling of the velocity defect in the wake, has been obtained. © 1997 John Wiley & Sons, Ltd.