Near-separating,self-similar turbulent boundary layer

The self-similar flow in a turbulent boundary layer, which is in a state close to separation as a result of the effect of adverse pressure gradient, is investigated. Such a boundary layer has a triple-deck asymptotic structure. Between outer and near-wall regions above the logarithmic sublayer, i. e. the constant-stress layer, an intermediate region — the gradient sublayer — is formed, where the shear stress varies linearly due to adverse longitudinal pressure gradient. In the external part of the gradient sublayer, the velocity profile obeys the square-root law. The velocity profile obtained from the solution for the outer region satisfies a slip condition on the wall. The slip value decreases as the similarity parameter increases and vanishes at the value of Ω = 0.0911, which corresponds to separation, here δ^* is the displacement thickness, and U and U′ are the free-stream velocity and its derivative with respect to the longitudinal coordinate. In this case, the exponent m in the law specifying the free-stream self-similar velocity distribution increases, with separation occurring not at the minimal value of m = −1/3, which corresponds to the strongest adverse pressure gradient, but at the value m = 0.228. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

The applicability of continuum models in the transitional regime of hypersonic flow over blunt bodies

Hypersonic rarefied gas flow over blunt bodies in the transitional flow regime (from continuum to free-molecule) is investigated. Asymptotically correct boundary conditions on the body surface are derived for the full and thin viscous shock layer models. The effect of taking into account the slip velocity and the temperature jump in the boundary condition along the surface on the extension of the limits of applicability of continuum models to high free-stream Knudsen numbers is investigated. Analytic relations are obtained, by an asymptotic method, for the heat transfer coefficient, the skin friction coefficient and the pressure as functions of the free-stream parameters and the geometry of the body in the flow field at low Reynolds number; the values of these coefficients approach their values in free-molecule flow (for unit accommodation coefficient) as the Reynolds number approaches zero. Numerical solutions of the thin viscous shock layer and full viscous shock layer equations, both with the no-slip boundary conditions and with boundary conditions taking into account the effects slip on the surface are obtained by the implicit finite-difference marching method of high accuracy of approximation. The asymptotic and numerical solutions are compared with the results of calculations by the Direct Simulation Monte Carlo method for flow over bodies of different shape and for the free-stream conditions corresponding to altitudes of 75–150 km of the trajectory of the Space Shuttle, and also with the known solutions for the free-molecule flow regine. The areas of applicability of the thin and full viscous shock layer models for calculating the pressure, skin friction and heat transfer on blunt bodies, in the hypersonic gas flow are estimated for various free-stream Knudsen numbers.     

Scaling laws for turbulent boundary layer flows in favorable and adverse pressure gradient

An analytical theory is proposed to describe incompressible plane and axisymmetric turbulent boundary layer flows in favorable and adverse pressure gradients for near-equilibrium conditions. Scaling laws for the mean velocity, the Reynolds stress components, and the skin friction have been established. A universal friction law makes it possible to represent the skin friction distributions corresponding to different Reynolds numbers and pressure gradients in terms of a function of one variable. The theory is based on general physical assumptions and does not involve any special hypotheses on the nature of the turbulent motion. (© 2009 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

微通道周期流动电位势及电粘性效应

求解了双电层的Poisson-Boltzmann方程和流体运动的Navier-Stokes方程,得到在周期压差作用下,二维微通道的周期流动电位势,流动诱导电场和液体流动速度的解析解．量纲分析表明,流体电粘性力与以下3个参数有关:1) 电粘性数,它表示定常流动时,通道最大电粘性力与压力梯度的比;2) 形状函数,它表示电粘性力在通道横截面的分布形态; 3) 耦合系数,它表示电粘性力的振幅衰减特征和相位差．分析结果表明,微通道周期流动诱导电场、流动速度与频率Reynolds数有关．在频率Reynolds数小于1时,流动诱导电场随频率Reynolds数变化很慢．在频率Reynolds数大于1时,流动诱导电场随频率Reynolds数的增加快速衰减．在通道宽度与双电层厚度比值较小情况下,电粘性效应对周期流动速度和流动诱导电场有重要影响．     

Hydrodynamic and thermal slip flow boundary layers over a flat plate with constant heat flux boundary condition

In this paper the boundary layer flow over a flat plat with slip flow and constant heat flux surface condition is studied. Because the plate surface temperature varies along the x direction, the momentum and energy equations are coupled due to the presence of the temperature gradient along the plate surface. This coupling, which is due to the presence of the thermal jump term in Maxwell slip condition, renders the momentum and energy equations non-similar. As a preliminary study, this paper ignores this coupling due to thermal jump condition so that the self-similar nature of the equations is preserved. Even this fundamental problem for the case of a constant heat flux boundary condition has remained unexplored in the literature. It was therefore chosen for study in this paper. For the hydrodynamic boundary layer, velocity and shear stress distributions are presented for a range of values of the parameter characterizing the slip flow. This slip parameter is a function of the local Reynolds number, the local Knudsen number, and the tangential momentum accommodation coefficient representing the fraction of the molecules reflected diffusively at the surface. As the slip parameter increases, the slip velocity increases and the wall shear stress decreases. These results confirm the conclusions reached in other recent studies. The energy equation is solved to determine the temperature distribution in the thermal boundary layer for a range of values for both the slip parameter as well as the fluid Prandtl number. The increase in Prandtl number and/or the slip parameter reduces the dimensionless surface temperature. The actual surface temperature at any location of x is a function of the local Knudsen number, the local Reynolds number, the momentum accommodation coefficient, Prandtl number, other flow properties, and the applied heat flux.     

Interaction between inner and outer layer in drag-reduced turbulent flows

Effects of wall-based skin-friction drag reduction strategies on the statistical properties of large-scale motions in moderate Reynolds number turbulent flows have been investigated by exploiting Direct Numerical Simulation of turbulent channels. To educe large scales, a new efficient parallel distributed memory algorithm has been implemented which delivers data-driven modes of increasing characteristic lengthscales: the Fast and Adaptive Bidimensional Empirical Mode Decomposition (FABEMD). The influence of wall-based skin friction reduction on large scales is studied by comparing single point statistics, such as r.m.s. fluctuations, and two-point statistics, as cross-correlation functions in controlled and uncontrolled channel flow fields at constant friction Reynolds number. The traditional way of observing large-scale footprinting at the wall, as cross-correlation of the streamwise velocity components at different wall distances, has been found to be unreliable when comparing drag-reduced flows, due to the arbitrary choice of a reference plane in the logarithmic layer. A more sound way of observing the footprinting via the correlation of the streamwise velocity with the friction velocity is addressed and shows an increase of the footprinting in drag-reduced flows. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Asymptotic description of incipient separation bubble bursting

The appearance of short laminar separation bubbles in high Reynolds number (Re) wall bounded flows due to appropriate adverse pressure gradient conditions is usually associated with minor effects on global flow properties (e.g. lift force). However, localized reverse flow regions are known to react very sensitively to perturbations and in further consequence may trigger the laminar-turbulent transition process or even cause global separation. The present investigation of marginally separated boundary layer flows is based on an asymptotic approach Re → ∞. Special emphasis is placed on solutions of the corresponding model equations which blow up within finite time indicating the ejection of a vortical structure and the emergence of shorter spatio-temporal scales reminiscent of the early transition scenario (‘ bubble bursting’ ). Within the framework of marginal separation theory, an alternative adjoint operator method is used to formulate evolution equations governing the viscous-inviscid interaction process in leading and higher order correction required for the study of later stages of the flow development. Their blow up structure specifies the initial condition of and the match to the subsequent triple deck stage. (© 2012 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Structure of weak laminar hydraulic jumps in single layer and two layer fluids

We consider the occurrence of hydraulic jumps in near critical single layer and two layer flows under the assumption that viscous effects are confined to a thin laminar boundary layer adjacent to the solid boundary. In the limit of large Reynolds number this leads to a structure problem formed by the classical triple deck equations supplemented with a novel nonlinear coupling condition which allows for the passage through the critical state. In the case of positive hydraulic jumps this passage is achieved by the local thickening of the boundary layer which acts as a viscous hump. Conversely, the pressure drop at the wall associated with negative hydraulic jumps causes the boundary layer to decrease locally thereby forming a local indentation required for the Froude number to pass through one in this case.     

壁面加热湍流耗散率标度指数测量的实验研究

对中等雷诺数下壁面常温和壁面加热的平板湍流边界层中速度和温度粗粒化的耗散率结构函数标度指数进行了实验测量.用热线风速仪测量了风洞中壁面常温和加热的平板湍流边界层中不同法向位置的流向速度分量和温度的时间序列信号,研究了由于湍流边界层近壁区域相干结构的存在而导致的非各向同性、非均匀性对湍流耗散率结构函数标度指数的影响,研究发现,中等雷诺数下壁面加热的边界条件和剪切湍流的平均速度梯度对速度和温度耗散率结构函数的标度指数没有影响,均匀各向同性湍流的耗散率结构函数标度指数的层次结构模型对壁面加热平板湍流边界层的速度和温度耗散率结构函数的标度指数也是适用的.     

Open-surface MHD flow over a curved wall in the 3-D thin-shear-layer approximation

3-D thin-shear-layer equations for flows of conducting fluids in a magnetic field have been derived in orthogonal body-oriented coordinates and then applied to the analysis of MHD open-surface flows over a curved wall. Unlike the classic boundary-layer-type equations, present ones permit information to be propagated upstream through the induced magnetic field. Another departure from the classic theory is that the normal momentum equation keeps the balance between the pressure gradient term, and those related to gravity, centrifugal forces, and Lorentz force. Thus, the normal pressure variations are allowed. The model describes basic 3-D effects due to the wall curvature and spatial variations of the applied magnetic field. As a particular case, equations for flows with rotational symmetry have been derived. Numerical calculations were performed for open-surface flows over a body of revolution under conditions relevant to a fusion reactor (Hartmann number is 8500). Some specific flow patterns, such as flow thickening and spiral-type flows, have been observed and discussed. A special attention has been paid to the analysis of the magnetic propulsion as a tool for the active flow control by applying an electric current. It has been shown that depending on the applied current, the axial pressure gradient can act as an adverse pressure gradient or propulsion force.     

Turbulent Boundary Layer Under the Influence of Adverse Pressure Gradient

The paper deals with the experimental analysis of turbulent boundary layer at the flat plate for large value of Reynolds number equal Re_θ ≈︁ 3000. The adverse pressure gradient generated by curvature of the upper wall corresponded to the case of pressure variation in axial compressor. The fully developed structure of turbulence was achieved by proper triggering of the boundary layer. The mean and turbulent flow-fields were investigated with the use of hot-wire technique while mean and instantaneous pressure fields were examined with piezoelectric transducers. The scaling and turbulence structure of fully developed turbulent boundary layer under the influence of adverse pressure gradient revealed the more pronounced contribution of outer region to the downstream development of turbulent boundary layer. (© 2009 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Dual solutions in mixed convection boundary layer flow of micropolar fluids

The steady mixed convection boundary layer flow over a vertical surface immersed in an incompressible micropolar fluid is considered in this paper. Employing suitable similarity transformations, the governing partial differential equations are transformed into ordinary differential equations, and the transformed equations are solved numerically by the Keller-box method. Numerical results are obtained for the skin friction coefficient and the local Nusselt number as well as the velocity, angular velocity and temperature profiles. Both cases of assisting and opposing buoyant flows are considered. It is found that dual solutions exist for the assisting flow, besides that usually reported in the literature for the opposing flow. Moreover, in contrast to the classical boundary layer theory, the separation point of the boundary layer is found to be distinct from the point of vanishing skin friction.     

Low-Reynolds-number modelling of flows over a backward-facing step

The complex turbulent flow behind a backward-facing step is modelled using a full Reynolds stress closure. In order to develop a closure model that can resolve the complex near-wall flow in the recirculation region and in the recovery region downstream of the reattachment point, the performance of a low and a high Reynolds number version of the full Reynolds stress closure is examined and compared. Furthermore, the effects of redistribution modelling on the calculated flow is studied by comparing the performance of three redistribution models: one return-to-isotropy model and two with mean-strain effects. The results are grid independent and show that the flow downstream of the step is best described by a low-Reynolds-number model that does not depend on the conventional wall function assumption. However, the skin friction behavior is correctly predicted by the stipulation of a wall function. Of the three redistribution models examined, the return-to-isotropy model gives results that are in excellent agreement with measurements. Finally, the calculated results are adversely affected by refining the redistribution models to include meanstrain effects.     

Asymptotic and numerical solutions of three‐dimensional boundary‐layer flow past a moving wedge

We consider a laminar boundary‐layer flow of a viscous and incompressible fluid past a moving wedge in which the wedge is moving either in the direction of the mainstream flow or opposite to it. The mainstream flows outside the boundary layer are approximated by a power of the distance from the leading boundary layer. The variable pressure gradient is imposed on the boundary layer so that the system admits similarity solutions. The model is described using 3‐dimensional boundary‐layer equations that contains 2 physical parameters: pressure gradient (β) and shear‐to‐strain‐rate ratio parameter (α). Two methods are used: a linear asymptotic analysis in the neighborhood of the edge of the boundary layer and the Keller‐box numerical method for the full nonlinear system. The results show that the flow field is divided into near‐field region (mainly dominated by viscous forces) and far‐field region (mainstream flows); the velocity profiles form through an interaction between 2 regions. Also, all simulations show that the subsequent dynamics involving overshoot and undershoot of the solutions for varying parameter characterizing 3‐dimensional flows. The pressure gradient (favorable) has a tendency of decreasing the boundary‐layer thickness in which the velocity profiles are benign. The wall shear stresses increase unboundedly for increasing α when the wedge is moving in the x‐direction, while the case is different when it is moving in the y‐direction. Further, both analysis show that 3‐dimensional boundary‐layer solutions exist in the range −1<α<∞. These are some interesting results linked to an important class of boundary‐layer flows.     

Laws of the wall for velocity and temperature in supersonic turbulent boundary layers

Scaling laws for velocity and temperature profiles in the near-wall region of sub- and supersonic turbulent boundary layers have been developed, which allow us to represent velocity and temperature profiles in compressible gas stream in terms of those in an incompressible boundary layer. They are obtained as asymptotic expansions of the solutions to the Reynolds equations in a small parameter — the Mach number based on the friction velocity and gas enthalpy on the wall. The leading term of the expansion for velocity corresponds to known Van Driest's formula. However, the obtained solution contains additional terms of order unity, which explains the contradiction between Van Driest's formula and experimental data. The law of the wall for temperature, which has been formulated for the first time, has an analogous structure. Besides the von Kármán constant and the turbulent Prandtl number in the logarithmic region, known for incompressible flow, the obtained relations contain three new universal constants, which do not depend on gas molecular properties and the specific heat ratio. (© 2017 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Modeling unsteady flow characteristics using smoothed particle hydrodynamics

Smoothed particle hydrodynamics (SPH) method has been extensively used to simulate unsteady free surface flows. The works dedicated to simulation of unsteady internal flows have been generally performed to study the transient start up of steady flows under constant driving forces and for low Reynolds number regimes. However, most of the fluid flow phenomena are unsteady by nature and at moderate to high Reynolds numbers. In this study, first a benchmark case (transient Poiseuille flow) is simulated to evaluate the ability of SPH to simulate internal transient flows at low and moderate Reynolds numbers (Re = 0.05, 500 and 1500). For this benchmark case, the performance of the two most commonly used formulations for viscous term modeling is investigated, as well as the effect of using the XSPH variant. Some points regarding using the symmetric form for pressure gradient modeling are also briefly discussed. Then, the application of SPH is extended to oscillating flows imposed by oscillating body force (Womersley type flow) and oscillating moving boundary (Stokes’ second problem) at different frequencies and amplitudes. There is a very good agreement between SPH results and exact solution even if there is a large phase lag between the oscillating pressure difference and moving boundary and the movement of the SPH particles generated. Finally, a modified formulation for wall shear stress calculations is suggested and verified against exact solutions. In all presented cases, the spatial convergence analysis is performed.     

壁湍流多尺度湍涡结构推广的自相似标度律

在水槽中测量了中等雷诺数下平板湍流边界层中的瞬时流向速度的时间序列,验证了Benzi提出的推广的自相似标度律,用子波变换将壁湍流脉动速度分解为多尺度湍涡结构的速度,研究了每一个尺度的湍涡速度结构函数的推广的自相似标度律。主要结论如下:湍流的统计性质是自相似的,这不仅适用于充分发展湍流,而且适用于中等雷诺数和低雷诺数湍流,而且具有相同的标度指数;推广的自相似标度律的适用的尺度范围远远大于惯性子区的范围,可以一直延伸至耗散区的尺度范围;推广的自相似标度律不仅适用于均匀各向同性湍流,也适用于剪切湍流如边界层湍流。     

Similarity solutions for mixed convection boundary layer flow over a permeable horizontal flat plate

The effects of suction and injection on steady laminar mixed convection boundary layer flow over a permeable horizontal flat plate in a viscous and incompressible fluid is investigated in this paper. The similarity solutions of the governing boundary layer equations are obtained for some values of the suction and injection parameter f₀, the constant exponent n of the wall temperature as well as the mixed convection parameter λ. The resulting system of nonlinear ordinary differential equations is solved numerically for both assisting and opposing flow regimes using a finite-difference scheme known as the Keller-box method. Numerical results for the reduced skin friction coefficient, the reduced local Nusselt number, and the velocity and temperature profiles are obtained for various values of the parameters considered. Dual solutions are found to exist for the opposing flow.     

The effect of variable viscosity on MHD viscoelastic fluid flow and heat transfer over a stretching sheet

An analysis has been carried out to study the momentum and heat transfer characteristics in an incompressible electrically conducting non-Newtonian boundary layer flow of a viscoelastic fluid over a stretching sheet. The partial differential equations governing the flow and heat transfer characteristics are converted into highly non-linear coupled ordinary differential equations by similarity transformations. The effect of variable fluid viscosity, Magnetic parameter, Prandtl number, variable thermal conductivity, heat source/sink parameter and thermal radiation parameter are analyzed for velocity, temperature fields, and wall temperature gradient. The resultant coupled highly non-linear ordinary differential equations are solved numerically by employing a shooting technique with fourth order Runge–Kutta integration scheme. The fluid viscosity and thermal conductivity, respectively, assumed to vary as an inverse and linear function of temperature. The analysis reveals that the wall temperature profile decreases significantly due to increase in magnetic field parameter. Further, it is noticed that the skin friction of the sheet decreases due to increase in the Magnetic parameter of the flow characteristics.     

On the initial phase of the triple deck stage in marginally separated flows

The method of matched asymptotic expansions is used to investigate marginally separated boundary layer flows (laminar or alternatively transitional separation bubbles) at high Reynolds numbers. Typical examples include, among others, the flow past slender airfoils at small to moderate angels of attack and channel flows with suction. As is well-known, classical (hierarchical) boundary layer computations usually break down under the action of an adverse pressure gradient on the flow, a scenario associated with the appearance of the Goldstein separation singularity. If, however, the parameter controlling the strength of the pressure gradient (the angle of attack or the relative suction rate in the examples mentioned above) is adjusted accordingly, the application of a local viscous-inviscid interaction strategy is capable of describing localized boundary layer separation. Moreover, taking into account unsteady effects and flow control devices allows the investigation of the conditions leading to forced or self-sustained vortex generation and the subsequent evolution process culminating in bubble bursting. Within the asymptotic formulation of this stage bubble bursting is associated with the formation of finite time singularities in the solution of the underlying equations and a corresponding break down. The distinct blow-up structure gives rise to a fully non-linear triple deck interaction stage featuring shorter spatio-temporal scales characteristic of the successive vortex evolution process. The paper will focus on the numerical treatment of the initial phase of the latter stage. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)