Analytical investigation of boundary layer growth and swirl intensity decay rate in a pipe

In this research, the developing turbulent swirling flow in the entrance region of a pipe is investigated analytically by using the boundary layer integral method. The governing equations are integrated through the boundary layer and obtained differential equations are solved with forth-order Adams predictor-corrector method. The general tangential velocity is applied at the inlet region to consider both free and forced vortex velocity profiles. The comparison between present model and available experimental data demonstrates the capability of the model in predicting boundary layer parameters (e.g. boundary layer growth, shear rate and swirl intensity decay rate). Analytical results showed that the free vortex velocity profile can better predict the boundary layer parameters in the entrance region than in the forced one. Also, effects of pressure gradient inside the boundary layer is investigated and showed that if pressure gradient is ignored inside the boundary layer, results deviate greatly from the experimental data.     

Asymptotic investigation of turbulent buoyant jets

Results of a theoretical investigation are presented for the averaged characteristics of turbulent slightly swirling jets discharged vertically upward into a stagnant homogeneous medium. A semi-empirical model of turbulence is used to close the system of initial turbulent boundary layer equations. Exact and approximate analytical solutions are suggested. The results of calculations are compared with experimental and numerical data of other authors.     

The flow structure of dilute gas-particle suspensions behind a shock wave moving along a flat surface

The asymptotic and numerical investigations of shock-induced boundary layers in gas-particle mixtures are presented.The Saffman lift force acting on a particle in a shear flow istaken into account.It is shown that particle migration across the boundary layer leads tointersections of particle trajectories.The corresponding modification of dusty gas model isproposed in this paper.The equations of two-phase sidewall boundary layer behind a shock wave moving at aconstant speed are obtained by using the method of matched asymptotic expansions.Themethod of the calculation of particle phase parameters in Lagrangian coordinates isdescribed in detail.Some numerical results for the case of small particle concentration aregiven.     

超声速平板边界层斜波失稳转捩过程研究

以5阶迎风和6阶对称紧致格式混合差分求解三维可压缩滤波Navier-Stokes方程,对Mach 数为4.5, Reynolds数为10000的空间发展平板边界层湍流进行了大涡模拟. 时间推进采用 紧致存储3阶Runge-Kutta方法,亚格子尺度模型为修正Smagorinsky涡黏性模型. 通过在 入口边界叠加一对线性最不稳定第一模态斜波扰动,数值模拟得到了平板层流边界层失稳转 捩直至湍流的演化过程. 对流场转捩过程中瞬时量及统计平均量的分析表明,数值模拟结果 与理论吻合,得到的Y型剪切层、交替\Lambda涡结构以及转捩后期的发卡涡结构的发展 变化与相关文献结果一致,湍流流谱定性合理.     

Numerical simulation of swirling flow in diffusers

A reduced form of Navier–Stokes equations is developed which does not have the usual minimum axial step size restriction. The equations are able to predict accurately turbulent swirling flow in diffusers. An efficient single sweep implicit scheme is developed in conjunction with a variable grid size domain-conforming co-ordinate system. The present scheme indicates good agreement with experimental results for (1) turbulent pipe flow, (2) turbulent diffuser flow, (3) turbulent swirling diffuser flow. The strong coupling between the swirl and the axial velocity profiles outside of the boundary layer region is demonstrated.     

Universal equations of three-dimensional boundary layer theory

We consider a parametric method for investigating three-dimensional laminar motion of an incompressible fluid in a boundary layer on a curved surface. It is found that the problem solution in the general case depends on four series of parameters, constructed from two components of the outer flow velocity and the two Lamé coefficients characterizing the shape of the immersed surface. From the general equations of the three-dimensional boundary layer we obtain a system of two universal equations which do not contain the characteristics of the outer flow. This system may be solved once and for all. As an example we consider the problem of the laminar boundary layer on the walls of an axisymmetric channel in the case of swirling outer flow. For this case we obtain numerical solutions of the system of universal equations in the local two-parameter approximation.     

Jeffrey-Hamel flow of Leslie-Ericksen nematic liquids

A similarity solution of the Leslie-Ericksen equations for nematic liquid crystals is obtained for flow between converging and diverging planar walls (Jeffrey-Hamel flow). There are three regions in the flow: extensional or compressional flow near the centerline, shear near the wall, and a wall boundary layer in which elastic stresses control the transition from the wall-induced orientation to the bulk behavior. The boundary layer thickness is obtained in closed form; the scaling with the Ericksen number depends on whether or not the boundary layer extends into the region of extensional flow. Imposition of a magnetic field with an azimuthal component in a converging flow can result in a Freedericksz-like transition from radial to transverse orientation at the center line at a critical field strength. This transition provides a new means to measure the irrotational viscosity λ₂.     

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.     

Three-dimensional coupled numerical model of creeping flow of viscous fluid

A three-dimensional coupled numerical model is developed to describe creeping flow in a computational domain that consists of a thick viscous layer overlaid with a thin multilayered viscous sheet. The density of the sheet is assumed to be lower than that of the layer. The model couples the Stokes equations describing the flow in the layer and the Reynolds equations describing the flow in the sheet. We investigate the long-time behavior of the flow in the sheet by using an asymptotic method and derive an ordinary differential equation for the sheet boundary displacements and the velocities at the interface between the sheet and the layer. The Stokes and Reynolds equations are coupled by applying the resulting equation as an internal boundary condition. Numerical implementation is based on a modified finite element method combined with the projection gradient method. The computational domain is discretized into rectangular hexahedra. Piecewise square basis functions are used. The model proposed enables different-type hydrodynamic equations to be coupled without any iterative improvements. As a result, the computational costs are reduced significantly in comparison with available coupled models. Numerical experiments confirm that the three-dimensional coupled model developed is of good accuracy.     

DNS and scaling laws from new symmetry groups of ZPG turbulent boundary layer flow

The Lie group, or symmetry approach, developed by Oberlack (see e.g. Oberlack [26] and references therein) is used to derive new scaling laws for various quantities of a zero pressure gradient turbulent boundary layer flow. The approach unifies and extends the work done by Oberlack for the mean velocity of stationary parallel turbulent shear flows. From the two-point correlation (TPC) equations the knowledge of the symmetries allows us to derive a variety of invariant solutions (scaling laws) for turbulent flows, one of which is the new exponential mean velocity profile that is found in the mid-wake region of flat-plate boundary layers. Further, a third scaling group was found in the TPC equations for the one-dimensional turbulent boundary layer. This is in contrast to the Navier–Stokes and Euler equations, which have one and two scaling groups, respectively. The present focus is on the exponential law in the outer region of turbulent boundary layer corresponding new scaling laws for one- and two-point correlation functions. A direct numerical simulation (DNS) of a flat plate turbulent boundary layer with zero pressure gradient was performed at two different Reynolds numbers Re=750,2240. The Navier–Stokes equations were numerically solved using a spectral method with up to 140 million grid points. The results of the numerical simulations are compared with the new scaling laws. TPC functions are presented. The numerical simulation shows good agreement with the theoretical results, however only for a limited range of applicability. PACS 02.20.-a, 47.11.+j, 47.27.Nz, 47.27.Eq     

LES study of global instability in annular jets

The paper presents results of numerical investigations of annular non–swirling jets. The research is performed with the use of large eddy simulations (LES) and a high–order code based on the compact differences and the Fourier approximation combined with a projection method for pressure–velocity coupling. We consider various Reynolds numbers and inner shear layer thicknesses and observe that both parameters affect the dynamics of the flow in such a way that the formation of the turbulent structures follows two different scenarios. For the particular group of the inlet conditions the flow pattern is found to be similar to low–swirling jets. The instantaneous flow–fields reveal the formation of the spiral structures located in the inner and outer mixing layers. The analysis shows that the structures are the result of the instability that leads to the precession of the recirculation region formed in the near–field. The winding sense of the spirals are clockwise or counterclockwise depending on the case and azimuthal distribution of perturbations.     

Nonlinear evolution analysis of T-S disturbance wave at finite amplitude in nonparallel boundary layers

ＩｎｔｒｏｄｕｃｔｉｏｎＮｏｎｌｉｎｅａｒｓｔａｂｉｌｉｔｙｐｒｏｂｌｅｍｉｎｎｏｎｐａｒａｌｌｅｌｂｏｕｎｄａｒｙｌａｙｅｒｓｎｅｅｄｓｍａｉｎｌｙｃｏｎｓｉｄｅｒｔｈｅｅｖｏｌｕｔｉｏｎｏｆｕｎｓｔａｂｉｌｉｔｙＴ＿Ｓｗａｖｅａｎｄｉｎｔｅｒａｃｔｉｏｎｂｅｔｗｅｅｎｗａｖｅｓｏｆｄｉｆｆｅｒｅｎｔｆｒｅｑｕｅｎｃｉｅｓ.Ｉｔｉｓｖｅｒｙｉｍｐｏｒｔａｎｔｔｏｓｉｍｕｌａｔｅａｃｃｕｒａｔｅｌｙｓｐａｔｉａｌｅｖｏｌｕｔｉｏｎｒｏｌｅｏｆｄｉｓｔｕｒｂａｎｃｅｉｎｂｏｕｎｄａｒｙｌａｙｅｒｓ,…     

Full coverage film cooling using compound angle

An experimental and numerical study is carried out on a cooling film issuing from a multiperforated wall of a simplified combustor. The objectives of this work are to achieve a better understanding of the dynamics of the film and to construct an experimental database on a simplified geometry in order to test numerical models. A parametric study of film cooling efficiency based on the direction of the cooling air injection is presented and shows that a swirling injection greatly enhances the cooling efficiency. As accounting for multiperforated walls in numerical simulations cannot be done at the jets scale because of computing resources, in this article are presented RANS computations performed using a uniform boundary condition to provide the injection of coolant. Two injection models are applied on this boundary and numerical results are compared to experimental data in the recovery region. The standard model is shown to be totally inappropriate while the multiperforation model delivers promising results although some weaknesses appear very close to the wall. To cite this article: B. Michel et al., C. R. Mecanique 337 (2009).     

采用不同黏性处理方法的宽无叶扩压器不稳定流动研究

径向无叶扩压器的全局稳定性可能受到核心主流失稳,出口回流与壁面边界层分离等因素影响,对于宽无叶扩压器,无黏核心主流与壁面边界层流动对不稳定扰动诱发的作用机理是当前研究的重点.本文首先通过数值计算获得了大宽度比孤立无叶扩压器平均流动,然后基于小扰动理论和周向均质假设,分别对欧拉方程与 Navier-Stokes 方程进行线性化,建立了基于无黏核心流动的稳定性分析方法,以及基于涡黏性与分子黏性的混合稳定性分析方法;通过与实验结果的对比,验证了混合稳定性分析方法预测所得流动失稳频率和全局直接模态的准确性;最后基于伴随方法获得了特征值的结构敏感性,揭示了不同黏性处理条件下宽无叶扩压器内全局不失稳扰动的源发区域.在只考虑核心主流的无黏条件下,宽无叶扩压器内流动不稳定扰动来源于流场中部,为二维的离心失稳;在同时考虑核心主流与边界层的作用时,宽无叶扩压器不稳定扰动不仅来源于扩压器流场中部的核心主流,壁面回流对于不稳定扰动的产生了重要影响.      

Improved methods for thin, surface boundary layer investigations

 New techniques are developed to improve the velocity flow-field measurement capability within a free-surface boundary layer region on which progressive capillary-gravity waves are present. Due to the extremely thin but rather vortical characteristics of the aforementioned boundary layer, conventional particle image velocimetry (PIV) methods fail to estimate velocity (and vorticity) vectors at an acceptable detection rate. This failure is a direct consequence of optimal PIV parameters that are difficult to achieve in practice for such flow situations. A new technique, Sub-pattern PIV, is developed. This method has features similar to both the super-resolution PIV (Keane et al. 1995) and the particle image distortion (PID) technique (Huang et al. 1993), but is predicated upon a very differential philosophy. Another difficulty that arises in experiments to investigate surface boundary layer flows is that the oscillating and deforming air–water interface has a mirror-like behavior that affects the images, and generates very noisy data. An alternative experimental setup that utilizes the Brewster angle phenomenon is adopted and the specular effects of the free-surface are removed successfully. This Brewster angle imaging, along with the Sub-pattern PIV technique, is used for the target application – a free-surface boundary layer investigation. It proved to be very effective. The methodology of both techniques is discussed, and the modified PIV procedure is validated by numerical probabilistic simulations. Application to the capillary-gravity wave boundary layer is presented in a subsequent paper. Received: 31 July 1997/Accepted: 4 February 1998     

Three-dimensional arbitrary Lagrangian–Eulerian numerical prediction method for non-linear free surface oscillation

A numerical prediction method has been proposed to predict non-linear free surface oscillation in an arbitrarily-shaped three-dimensional container. The liquid motions are described with Navier–Stokes equations rather than Laplace equations which are derived by assuming the velocity potential. The profile of a liquid surface is precisely represented with the three-dimensional curvilinear co-ordinates which are regenerated in each computational step on the basis of the arbitrary Lagrangian–Eulerian (ALE) formulation. In the transformed space, the governing equations are discretized on a Lagrangian scheme with sufficient numerical accuracy and the boundary conditions near the liquid surface are implemented in a complete manner. In order to confirm the applicability of the present computational technique, numerical simulations are conducted for the free oscillations of viscid and inviscid liquids and for highly non-linear oscillation. In addition, non-linear sloshing motions caused by horizontal and vertical excitations and a transition from non-linear sloshing to swirling are numerically predicted in three-dimensional cylindrical containers. Conclusively, it is shown that these sloshing motions associated with high non-linearity are reasonably predicted with the present numerical technique. © 1998 John Wiley & Sons, Ltd.     

Laminar boundary layer in a swirling flow

The spatial non-self-similar boundary layer in a compressible gas in a swirling flow is studied. Boundary-layer equations are written in variables ensuring constancy of the coefficients of first derivatives and are solved by the finite-difference method. Boundary-layer peculiarities in the presence of a return circulation region in the channel are clarified.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 43–49, January–February, 1976.     

DNS of a spatially developing turbulent boundary layer with passive scalar transport

A direct numerical simulation (DNS) of a spatially developing turbulent boundary layer over a flat plate under zero pressure gradient (ZPG) has been carried out. The evolution of several passive scalars with both isoscalar and isoflux wall boundary condition are computed during the simulation. The Navier–Stokes equations as well as the scalar transport equation are solved using a fully spectral method. The highest Reynolds number based on the free-stream velocity U_∞ and momentum thickness θ is Re_θ=830, and the molecular Prandtl numbers are 0.2, 0.71 and 2. To the authors’ knowledge, this Reynolds number is to date the highest with such a variety of scalars. A large number of turbulence statistics for both flow and scalar fields are obtained and compared when possible to existing experimental and numerical simulations at comparable Reynolds number. The main focus of the present paper is on the statistical behaviour of the scalars in the outer region of the boundary layer, distinctly different from the channel-flow simulations. Agreements as well as discrepancies are discussed while the influence of the molecular Prandtl number and wall boundary conditions is also highlighted. A Pr scaling for various quantities is proposed in outer scalings. In addition, spanwise two-point correlation and instantaneous fields are employed to investigate the near-wall streak spacing and the coherence between the velocity and the scalar fields. Probability density functions (PDF) and joint probability density functions (JPDF) are shown to identify the intermittency both near the wall and in the outer region of the boundary layer. The present simulation data will be available online for the research community.     

Large-scale spanwise periodicity in a turbulent boundary layer induced by highly ordered and directional surface roughness

The effect of converging–diverging riblet-type surface roughness (riblets arranged in a ‘herringbone’ pattern) are investigated experimentally in a zero pressure gradient turbulent boundary layer. For this initial parametric investigation three different parameters of the surface roughness are analysed in detail; the converging–diverging riblet yaw angle α, the streamwise fetch or development length over the rough surface F_x and the viscous-scaled riblet height h⁺. It is observed that this highly directional surface roughness pattern induces a large-scale spanwise periodicity onto the boundary layer, resulting in a pronounced spanwise modification of the boundary layer thickness. Hot-wire measurements reveal that above the diverging region, the local mean velocity increases while the turbulent intensity decreases, resulting in a thinner overall boundary layer thickness in these locations. The opposite situation occurs over the converging region, where the local mean velocity is decreased and the turbulent intensity increases, producing a locally thicker boundary layer. Increasing the converging–diverging angle or the viscous-scaled riblet height results in stronger spanwise perturbations. For the strongest convergent–divergent angle, the spanwise variation of the boundary layer thickness between the diverging and converging region is almost a factor of two. Such a large variation is remarkable considering that the riblet height is only 1% of the unperturbed boundary layer thickness. Increasing the fetch seems to cause the perturbations to grow further from the surface, while the overall strength of the induced high and low speed regions remain relatively unaltered. Further analysis of the pre-multiplied energy spectra suggests that the surface roughness has modified or redistributed the largest scale energetic structures.