On the influence of the exit conditions on the entrainment rate in the development region of a free, round, turbulent jet

The development region of a free round turbulent jet was experimentally investigated in order to determine the dependence of the jet development on the exit conditions. In particular the influence of the exit conditions on the entrainment rate in both the axisymmetrix mixing layer region and at a distance of 20D from the nozzle was investigated.The entrainment rate at 20D depends noticeably on the exit boundary layer state, i.e. laminar or turbulent, and on the turbulence intensity of the core. A laminar exit boundary layer showed the highest value of entrainment rate. A turbulent one reduced the entrainment rate at 20D by 15%. A grid placed at the nozzle exit, in order to increase the turbulence intensity of the core, reduced the entrainment rate at 20D by 40%. It must be mentioned that this grid also disturbed the exit boundary layer.An empirical relation for the entrainment rate at the end of the mixing layer region as a function of the mixing layer parameters, is successful in predicting the entrainment rate at 20D within an accuracy of 5%. This result suggests a strong dependence of the jet development region on mixing layer processes, i.e. the development of organized structures and their coalescence. Possibly there is no asymptotic region which is fully independent of the mixing layer processes, and which is not more or less influenced by the exit conditions via the influence on these processes.     

Large eddy simulation in a turbulent channel flow using exit boundary conditions

The influence of the exit boundary conditions on the vanishing first derivative of the velocity components and constant pressure on the large eddy simulation of the fully developed turbulent channel flow has been investigated for equidistant and stretched grids at the channel exit. Results show that the chosen exit boundary conditions introduce some small disturbances that are mostly damped by the grid stretching. The difference of rms values between the fully developed turbulent channel flow with periodicity conditions and the fully developed channel flow using inlet and the exit boundary conditions is less than 10% for the equidistant grids and less than 5% for the stretched grids. The chosen boundary conditions are of interest because they may be used in complex problems with back flow at the exit. Copyright © 1999 John Wiley & Sons, Ltd.     

Shock stability in magnetogasdynamic channel flows

The question of shock stability in a perfect-gas channel flow was examined in [1] in the onedimensional approximation under various assumptions: the disturbances are not reflected from the channel exit section, weak shock, etc. The results were found to coincide for two specific forms of the boundary conditions at the channel exit, from which it was concluded that the shock was not sensitive to the exit boundary condition. In [2] the question of shock stability was studied numerically in relation to a conducting-gas flow in a flat channel of constant cross section in the presence of a magnetic field (zero electric field intensity). It was established that the shock stability is significantly affected by the form of the conductivity law. A condition for the limiting regime between the stable and unstable regions was also given for flow with a shock wave. It was assumed that the pressure in the channel exit section is given. In this paper the effect of the exit boundary condition on shock stability in gasdynamic and magnetogasdynamic flows is demonstrated for small magnetic Reynolds numbers. Stability criteria are obtained for shocks near the channel exit for a specific exit condition. The influence of electromagnetic effects (conductivity law, electric load factor) on shock stability is investigated.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, Vol. 11, No. 1, pp. 16–23, January–February, 1970.The author is grateful to A. G. Kulikovskii for discussing his work.     

Local similarity solution for the flow of a “second-grade” viscoelastic fluid above a moving plate

The boundary layer flow of a viscoelastic fluid of the second-grade type over a rigid continuous plate moving through an otherwise quiescent fluid with constant velocity U is studied. Assuming the flow to be laminar and two-dimensional, local similarity solution is found with fluid's elasticity and plate's withdrawal speed as the main variables. Results are presented for velocity profiles, boundary layer thickness, wall skin friction coefficient and fluid entrainment in terms of the local Deborah number. A marked formation of boundary layer is predicted, even at low Reynolds numbers, provided the Deborah number is sufficiently large. The boundary layer thickness and the wall skin friction coefficient are found to scale with fluid's elasticity—both decreasing the higher the fluid's elasticity. The amount of fluid entrained is also predicted to decrease whenever a fluid exhibits elastic behavior.     

Consistent inlet and outlet boundary conditions for particle methods

In this paper, simple and consistent open boundary conditions are presented for the numerical simulation of viscous incompressible laminar flows. The present approach is based on an arbitrary Lagrangian-Eulerian particle method using upwind interpolation. Three kinds of inlet/outlet boundary conditions are proposed for particle methods, a pressure specified inlet/outlet condition, a velocity profile specified inlet/outlet condition, and a fully developed flow outlet condition. These inlet/outlet conditions are realized by using boundary particles and modification to the physical value such as velocity. Poiseuille flows, flows over a backward-facing step, and flows in a T-shape branch are calculated. The results are compared with those of mesh-based methods such as the finite volume method. The method presented herein exhibits accuracy and numerical stability.     

Existence and Stability of Multidimensional Transonic Flows through an Infinite Nozzle of Arbitrary Cross-Sections

We establish the existence and stability of multidimensional steady transonic flows with transonic shocks through an infinite nozzle of arbitrary cross-sections, including a slowly varying de Laval nozzle. The transonic flow is governed by the inviscid potential flow equation with supersonic upstream flow at the entrance, uniform subsonic downstream flow at the exit at infinity, and the slip boundary condition on the nozzle boundary. Our results indicate that, if the supersonic upstream flow at the entrance is sufficiently close to a uniform flow, there exists a solution that consists of a C ^1,α subsonic flow in the unbounded downstream region, converging to a uniform velocity state at infinity, and a C ^1,α multidimensional transonic shock separating the subsonic flow from the supersonic upstream flow; the uniform velocity state at the exit at infinity in the downstream direction is uniquely determined by the supersonic upstream flow; and the shock is orthogonal to the nozzle boundary at every point of their intersection. In order to construct such a transonic flow, we reformulate the multidimensional transonic nozzle problem into a free boundary problem for the subsonic phase, in which the equation is elliptic and the free boundary is a transonic shock. The free boundary conditions are determined by the Rankine–Hugoniot conditions along the shock. We further develop a nonlinear iteration approach and employ its advantages to deal with such a free boundary problem in the unbounded domain. We also prove that the transonic flow with a transonic shock is unique and stable with respect to the nozzle boundary and the smooth supersonic upstream flow at the entrance.     

Laser-Doppler measurements of laminar and turbulent flow in a pipe bend

Laser-Doppler measurements are reported for laminar and turbulent flow through a 90° bend of circular cross-section with mean radius of curvature equal to 2.8 times the diameter. The measurements were made in cross-stream planes 0.58 diameters upstream of the bend inlet plane, in 30, 60 and 75° planes in the bend and in planes one and six diameters downstream of the exit plane. Three sets of data were obtained: for laminar flow at Reynolds numbers of 500 and 1093 and for turbulent flow at the maximum obtainable Reynolds number of 43 000. The results show the development of strong pressure-driven secondary flows in the form of a pair of counter-rotating vortices in the streamwise direction. The strength and character of the secondary flows were found to depend on the thickness and nature of the inlet boundary layers, inlet conditions which could not be varied independently of Reynolds number. The quantitative anemometer measurements are supported by flow visualization studies. Refractive index matching at the fluid-wall interface was not used; the measurements consist, therefore, of streamwise components of mean and fluctuating velocities only, supplemented by wall pressure measurements for the turbulent flow. The displacement of the laser measurement volume due to refraction is allowed for in simple geometrical calculations. The results are intenden for use as benchmark data for calibrating flow calculation methods.     

CABARET solutions on graphics processing units for NASA jets: Grid sensitivity and unsteady inflow condition effect

The GPU CABARET method for solving the Navier–Stokes equations coupled with the Ffowcs Williams–Hawkings scheme for far-field noise predictions is applied for conditions of the NASA SHJAR experiment corresponding to Set Point 3 and 7 in accordance with Tanna's classification. The questions addressed include the sensitivity of the flow and noise spectra solutions to the grid resolution and the inflow condition at the nozzle exit. To study the grid sensitivity, several “hand-made” multi-block curvilinear grids are considered along with a simple hanging-nodes-type grid that was automatically generated with OpenFOAM, whose solutions are cross-verified. To study the effect of the inflow jet condition, the flow and noise solutions based on the laminar inflow condition for Set Point 7 case are compared with the same based on modifying the interior nozzle geometry with a turbulence grid to generate the initial unsteadiness inside the nozzle so that both the centerline velocity fluctuations and the jet Mach number at the nozzle exit are preserved in accordance with the experiment. The numerical solutions obtained are compared with the experimental data and reference LES solutions available in the literature.     

Uniqueness of the self-similar solutions of three-dimensional laminar boundary-layer equations

The equations of the three-dimensional laminar boundary layer on lines of flow outflow and inflow are studied for conical outer flow under the assumption that the Prandtl number and the productρμ are constant. It is shown that in the case of a positive velocity gradient of the secondary flow (α₁>0) the additional conditions which result from the physical flow pattern determine a unique solution of the system of boundary-layer equations. For a negative velocity gradient of the secondary flow (α₁≤0) these conditions are satisfied by two solutions. An approximate solution is obtained for the boundary layer equations which is in rather good agreement with the numerical integration results. Compressible gas flow in a three-dimensional laminar boundary layer is described by a system of nonlinear differential equations whose solution is not unique for given boundary conditions. Therefore additional conditions resulting from the physical pattern of the gas flow are imposed on the resulting solution. In the solution of problems with a negative pressure gradient these additional conditions are sufficient for a unique selection of the solution of the boundary-layer equations. However, in the case of a positive pressure gradient the solution of the boundary-layer equations satisfying the boundary and additional conditions may not be unique. In particular, in [1] in a study of a three-dimensional laminar boundary layer in the vicinity of the stagnation point it was shown that for $$c = {{\frac{{\partial v_e }}{{\partial y}}} \mathord{\left/ {\vphantom {{\frac{{\partial v_e }}{{\partial y}}} {\frac{{\partial u_e }}{{\partial x}}}}} \right. \kern-\nulldelimiterspace} {\frac{{\partial u_e }}{{\partial x}}}} > 0$$ the solution is unique, while for c<0 there are two solutions. In the present paper we study the question of the uniqueness of the self-similar solution of the three-dimensional laminar boundary-layer equations on lines of flow outflow and inflow for a conical outer flow.     

Numerical investigations of radial jet reattachment flows

A finite volume computational scheme to solve the Navier-Stokes equations for the laminar flow fields of partially enclosed axial and radial jets impinging on a flat plate has been devised and tested. This scheme is based on the SIMPLEC technique. However, because of the backflow at the ‘outflow’ boundary, the SIMPLEC pressure correction technique has to be modified. The need for this modification, necessitated by the convergence failure, showed the ‘hidden’ pressure boundary condition of SIMPLE-type techniques. Test computations with the present scheme for flows in a channel with a built-in cylinder show that the location of the exit boundary affects very slightly the separated flow behind the cylinder. Computed Squire jet flows compare quite well with the available analytical solution. Finally, impinging radial jets have been computed for different Reynolds numbers. The results show the critical Reynolds number below which a steady solution is obtained and above which periodic and eventually chaotic flows result.     

Prediction of amount of entrained droplets in vertical annular two-phase flow

Prediction of amount of entrained droplets or entrainment fraction in annular two-phase flow is essential for the estimation of dryout condition and analysis of post dryout heat transfer in light water nuclear reactors and steam boilers. In this study, air–water and organic fluid (Freon-113) annular flow entrainment experiments have been carried out in 9.4 and 10.2 mm diameter test sections, respectively. Both the experiments covered three distinct pressure conditions and wide range of liquid and gas flow conditions. The organic fluid experiments simulated high pressure steam–water annular flow conditions. In each experiment, measurements of entrainment fraction, droplet entrainment rate and droplet deposition rate have been performed by using the liquid film extraction method. A simple, explicit and non-dimensional correlation developed by Sawant [Sawant, P.H., Ishii, M., Mori, M., 2008. Droplet entrainment correlation in vertical upward co-current annular two-phase flow. Nucl. Eng. Des. 238 (6), 1342–1352] for the prediction of entrainment fraction is further improved in this study in order to account for the existence of critical gas and liquid flow rates below which no entrainment is possible.Additionally, a new correlation is proposed for the estimation of minimum liquid film flow rate at the maximum entrainment fraction condition. The improved correlation successfully predicted the newly collected air–water and Freon-113 entrainment fraction data. Furthermore, the correlations satisfactorily compared with the air–water, helium–water and air–genklene experimental data measured by Willetts [Willetts, I.P., 1987. Non-aqueous annular two-phase flow. D.Phil. Thesis, University of Oxford]. However, comparison of the correlations with the steam–water data available in literature showed significant discrepancies. It is proposed that these discrepancies might have been caused due to the inadequacy of the liquid film extraction method used to measure the entrainment fraction or due to the change in mechanism of entrainment under high liquid flow conditions.     

Study of parameters and entrainment of a jet in cross-flow arrangement with transition at two low Reynolds numbers

Investigation of the mixing process is one of the main issues in chemical engineering and combustion and the configuration of a jet into a cross-flow (JCF) is often employed for this purpose. Experimental data are gained for the symmetry plane in a JCF-arrangement of an air flow using a combination of particle image velocimetry (PIV) with laser-induced fluorescence (LIF). The experimental data with thoroughly measured boundary conditions are complemented with direct numerical simulations, which are based on idealized boundary conditions. Two similar cases are studied with a fixed jet-to-cross-flow velocity ratio of 3.5 and variable cross-flow Reynolds numbers equal to 4,120 and 8,240; in both cases the jet issues from the pipe at laminar conditions. This leads to a laminar-to-turbulent transition, which depends on the Reynolds number and occurs quicker for the case with higher Reynolds number in both experiments and simulations as well. It was found that the Reynolds number only slightly affects the jet trajectory, which in the case with the higher Reynolds number is slightly deeper. It is attributed to the changed boundary layer shape of the cross-flow. Leeward streamlines bend toward the jet and are responsible for the strong entrainment of cross-flow fluid into the jet. Velocity components are compared for the two Reynolds numbers at the leeward side at positions where strongest entrainment is present and a pressure minimum near the jet trajectory is found. The numerical simulations showed that entrainment is higher for the case with the higher Reynolds number. The latter is attributed to the earlier transition in this case. Fluid entrainment of the jet in cross-flow is more than twice stronger than for a similar flow of a jet issuing into a co-flowing stream. This comparison is made along the trajectory of the two jets at a distance of 5.5 jet diameters downstream and is based on the results from the direct numerical simulations and recently published experiments of a straight jet into a co-flow. Mixing is further studied by means of second-order statistics of the passive scalar variance and the Reynolds fluxes. Windward and leeward sides of the jet exhibit different signs for the time-averaged streamwise Reynolds flux 〈v _x ′c′〉. The large coherent structures which contribute to this effect are investigated by means of timely correlated instantaneous PIV-LIF camera snapshots and their contribution to the average statistics of 〈v _x ′c′〉 are discussed. The discussion on mixing capabilities of the jet in cross-flow is supported by simulation results showing the instantaneous three-dimensional coherent structures defined in terms of the pressure fluctuations.     

强吸气旋转圆筒壁面湍流边界层建模及计算

提出了湍流边界层的一种简单、快速计算方法, 用以求解强吸气作用下旋转圆筒表面边界层流动. 首先, 理论分析了同心圆筒间的旋转流体运动, 外筒静止、内筒旋转且为多孔吸气条件. 强吸气情况下旋转流动主要表现为内筒壁面附近的边界层流动, 基于这一事实得到了周向速度分布的解析表达式. 其次, 通过引入新参数扩展Cebeci-Smith代数湍流模型, 使其能考虑流线曲率、壁面吸气、低Reynolds数效应等因素. 针对这些因素的综合影响, 采用解析修正和经验参数对模型进行调整. 同时, 基于Reynolds应力湍流模型的仿真结果, 校准代数湍流模型中的经验参数. 最后, 给出基于广义Cebeci-Smith湍流模型的旋转壁面边界层流动的迭代算法, 该算法适用于需要特殊迭代过程的轴向及周向流动均匀情况. 计算了不同旋转速度和吸气强度组合工况下的边界层流动, 其周向速度和湍流强度分布与基于Reynolds应力湍流模型的计算结果非常接近. 并且表明, 当Reynolds应力湍流模型数值模拟预测内筒边界层为稳定层流时, 该方法也再现了相同初始条件下的层流边界层.      

Einfluß der thermischen Randbedingung auf den laminaren Wärmeübergang im Kreisrohr bei hydrodynamischer Einlaufströmung

Heat transfer characteristics of hydrodynamically developing laminar flow in a circular duct with different thermal boundary conditions were calculated by solving the equations of continuity, motion and energy in finite difference form. Results are presented for linear, sinusoidal and exponential variations of the prescribed wall heat flux along the duct length. A comparison shows that the influence of the thermal boundary condition on heat transfer increases with increasing development of velocity and temperature profiles. As a side result an improved correlation for heat transfer with constant wall heat flux in hydrodynamically developing flow is presented.     

Mathematical model of heat transfer in a fluid with account for its relaxation properties

Using the terms that take account for the temporal and spatial nonlocality (time variation of the heat flux and the temperature gradient) in the formula of Fourier’s law for the heat flux a differential equation for a fluid in motion is derived that contains the second time derivative and themixed derivative with respect to the spatial and temporal variables. Numerical solution of the problem of heat transfer in the laminar fluid flow in a plane channel demonstrates that, in view of the lag in the time variation of the heat flux from zero to a certain maximum value, the boundary condition of the first kind (thermal shock) cannot be instantaneously realized. The process of its stabilization on the wall is characterized by a certain time interval, whose duration is determined by the relaxation properties of the fluid. At large values of the dimensionless coefficients of the heat flux relaxation and the temperature gradient the boundary condition of the first kind can be realized only as the steady state is attainted, as Fo→∞. In this case, the flow does not contain temperature jumps and negative temperature values.     

Radiation effects on the thermal boundary layer flow over a moving plate with convective boundary condition

The steady laminar boundary layer flow over a moving plate in a moving fluid with convective surface boundary condition and in the presence of thermal radiation is investigated in this paper. Under certain conditions, the present problem reduces to the classical Blasius and Sakiadis problems. The effects of radiation and convective parameters on the thermal field are thoroughly examined and discussed. Dual solutions are found to exist when the plate and the fluid move in the opposite directions.     

Is Helmholtz Resonance a Problem for Micro-jet Actuators?

A theoretical analysis is described that determines the conditions for Helmholtz resonance for a popular class of self-contained microjet actuator used in both synthetic- and pressure-jump (pulse-jet) mode. It was previously shown that the conditions for Helmholtz resonance are identical to those for optimizing actuator performance for maximum mass flux. The methodology is described for numerical-simulation studies on how Helmholtz resonance affects the interaction of active and nominally inactive micro-jet actuators with a laminar boundary layer. Two sets of numerical simulations were carried out. The first set models the interaction of an active actuator with the boundary layer. These simulations confirm that our criterion for Helmholtz resonance is broadly correct. When it is satisfied we find that the actuator cannot be treated as a predetermined wall boundary condition because the interaction with the boundary layer changes the pressure difference across the exit orifice thereby affecting the outflow from the actuator. We further show that strong inflow cannot be avoided even when the actuator is used in pressure-jump mode. In the second set of simulations two-dimensional Tollmien–Schlichting waves, with frequency comparable with, but not particularly close to, the Helmholtz resonant frequency, are incident on a nominally inactive micro-jet actuator. The simulations show that under these circumstances the actuators act as strong sources of 3D Tollmien–Schlichting waves. It is surmised that in the real-life aeronautical applications with turbulent boundary layers broadband disturbances of the pressure field, including acoustic waves, would cause nominally inactive actuators, possibly including pulsed jets, to act as strong disturbance sources. Should this be true it would probably be disastrous for engineering applications of such massless microjet actuators for flow control.     

MODELLING AND CALCULATION OF THE TURBULENT BOUNDARY LAYER ON A ROTATING CYLINDER SURFACE WITH STRONG SUCTION 1)

提出了湍流边界层的一种简单、快速计算方法, 用以求解强吸气作用下旋转圆筒表面边界层流动. 首先, 理论分析了同心圆筒间的旋转流体运动, 外筒静止、内筒旋转且为多孔吸气条件. 强吸气情况下旋转流动主要表现为内筒壁面附近的边界层流动, 基于这一事实得到了周向速度分布的解析表达式. 其次, 通过引入新参数扩展Cebeci-Smith代数湍流模型, 使其能考虑流线曲率、壁面吸气、低Reynolds数效应等因素. 针对这些因素的综合影响, 采用解析修正和经验参数对模型进行调整. 同时, 基于Reynolds应力湍流模型的仿真结果, 校准代数湍流模型中的经验参数. 最后, 给出基于广义Cebeci-Smith湍流模型的旋转壁面边界层流动的迭代算法, 该算法适用于需要特殊迭代过程的轴向及周向流动均匀情况. 计算了不同旋转速度和吸气强度组合工况下的边界层流动, 其周向速度和湍流强度分布与基于Reynolds应力湍流模型的计算结果非常接近. 并且表明, 当Reynolds应力湍流模型数值模拟预测内筒边界层为稳定层流时, 该方法也再现了相同初始条件下的层流边界层.     

Simultaneous velocity and passive scalar concentration measurements in low Reynolds number neutrally buoyant turbulent round jets

Velocity and scalar concentration characteristics of low Reynolds number (Re) neutrally buoyant turbulent round jets were studied using coupled particle image velocimetry and laser induced fluorescence. Experiments were conducted on a jet with a fully developed pipe exit profile at Re = 1,500 and Re = 4,000. Measurements were made in the far field (60 < x/D < 80). Results show that the centerline velocity decay constant increases with Re and the virtual origin is located far from the jet exit. The Re does not have such an influence on the scalar concentration decay constant and scalar field virtual origin. Centerline turbulent intensities show evidence of becoming self-similar past x/D ≈ 75 but the same is not observed in the scalar fluctuations. The latter are strongly influenced by the Re, increasing as entrainment of ambient flow on the jet axis increases as the flow approaches laminar conditions. Scalar fluxes showed a marked Re dependence as well. Results were used to estimate the turbulent Schmidt number which decreases as the Re decreases.     

Flow and Mixing Fields for Transported Scalar PDF Simulations of a Piloted Jet Diffusion Flame (‘Delft Flame III’)

Numerical simulation results are presented for ‘Delft Flame III’, a piloted jet diffusion flame with strong turbulence–chemistry interaction. While pilot flames emerge from 12 separate holes in the experiments, the simulations are performed on a rectangular grid, under the assumption of axisymmetry. In the first part of the paper, flow and mixing field results are presented with a non-linear first order k–ε model, with the transport equation for ε based on a modeled enstrophy transport equation, for cold and reactive flows. For the latter, the turbulence model is applied in combination with pre-assumed β-PDF modeling for the turbulence–chemistry interaction. The mixture fraction serves as conserved scalar. Two chemistry models are considered: chemical equilibrium and a steady laminar flamelet model. The importance of the turbulence model is highlighted. The influence of the chemistry model is noticeable too. A procedure is described to construct appropriate inlet boundary conditions. Still, the generation of accurate inlet boundary conditions is shown to be far less important, their effect being local, close to the nozzle exit. In the second part of the paper, results are presented with the transported scalar PDF approach as turbulence–chemistry interaction model. A C₁ skeletal scheme serves as chemistry model, while the EMST method is applied as micro-mixing model. For the transported PDF simulations, the model for the pilot flames, as an energy source term in the mean enthalpy transport equation, is important with respect to the accuracy of the flow field predictions. It is explained that the strong influence on the flow and mixing field is through the turbulent shear stress force in the region, close to the nozzle exit.