Further consideration of the shape-factor relationship for turbulent boundary layers

The lag-entrainment predictive scheme developed by Green et al. has been modified to include the pressure-gradient parameter Π₁. In the original model suggested by Green et al. the mass-flow shape factor H₁ is related to the common shape factor H, H₁ = f(H). In the present model H₁ is related to H, Reynolds number based on the local momentum thickness θ, and Π₁; thus H₁ = f(H, R_eθ, Π₁). The modified formula for H₁, is introduced into the original lag-entrainment integral model. Calculations are made to examine the present model for the predictions of the development of boundary layers approaching separation studied experimentally by the authors. Slightly improved predictions are obtained using the model developed by El Telbany et al. However, the present model proved to give an improved representation of the development of wall shear stress in cases the two-equation turbulence model proved to be unsuccessful.     

Flow in a cylinder driven by rotating split-disk endwalls

Flow of an incompressible viscous fluid contained in a cylindrical vessel (radius R, height H) is considered. Each of the cylinder endwalls is split into two parts which rotate steadily about the central axis with different rotation rates: the inner disk (r < r₁) rotating at Ω₁, and the outer annulus (r₁ < r < R) rotating at Ω₂. Numerical solutions to the axisymmetric Navier-Stokes equations are secured for small system Ekman numbers E ( v/(ΩH²)). In the linear regime, when the Rossby number Ro

Full-size image (1K)

, the numerical results are shown to be compatible with the theoretical prediction as well as the available experimental measurements. Emphasis is placed on the results in the nonlinear regime in which Ro is finite. Details of the structures of azimuthai and meridional flows are presented by the numerical results. For a fixed Ekman number, the gross features of the flow remain qualitatively unchanged as Ro increases. The meridional flows are characterized by two circulation cells. The shear layer is a region of intense axial flow toward the endwall and of vanishing radial velocity. The thicknesses of the shear layer near r = r₁ and the Ekman layer on the endwall scale with E and E , respectively. The numerical results are consistent with these scalings.     

A general correlation for the local loss coefficient in Newtonian axisymmetric sudden expansions

Results from numerical simulations and guidance from an approximated corrected-theory, developed by Oliveira and Pinho (1997), (Oliveira, P.J. and Pinho, F.T. 1997. Pressure drop coefficient of laminar Newtonian flow in axisymmetric sudden expansions. Int. J. Heat and Fluid flow 18, 518–529) have been used to arrive at a correlation expressing the irreversible loss coefficient for laminar Newtonian flow in axisymmetric sudden expansions. The correlation is valid for the ranges 1.5 < D₂/D₁ < 4 and 0.5 < Re < 200 with errors of less than 5%, except for 25 < Re < 100 where the error could be as much as 7%. The recirculation bubble length is also presented for the same range of conditions and the pressure recovery coefficient was calculated for Reynolds numbers above 15.     

Inertial instability of the Stewartson -layer

Inertial stability of a vertical shear layer (Stewartson E^1/4-layer) on the sidewall of a cylindrical tank with respect to stationary axisymmetric perturbations is inverstigated by means of a linear theory. The stability is determined by two non-dimensional parameters, the Rossby number Ro = U/2ΩL and Ekman number E = v/ΩH², where U and L = (E/4)1/4H are the characteristic velocity and width of the shear layer, respectively, Ω the angular velocity of the basic rotation, v the kinematic viscosity and H the depth of the tank.

For a given Ekman number, the flow is more unstable for larger values of the Rossby number. For E = 10⁻⁴, which is a typical value of the Ekman number realized in rotating tank experiments, the critical Rossby number Ro_c for instability and the critical axial wavenumber m_c non-dimensionalized by L⁻¹ are found to be 1.3670 and 8.97, respectively. The value of Ro_c increases and that of m_c decreases with increasing E.     

Flow regime transitions due to cavitation in the flow through an orifice

This paper presents both experimental and theoretical aspects of the flow regime transitions caused by cavitation when water is passing through an orifice. Cavitation inception marks the transition from single-phase to two-phase bubbly flow; choked cavitation marks the transition from two-phase bubbly flow to two-phase annular jet flow.

It has been found that the inception of cavitation does not necessarily require that the minimum static pressure at the vena contracta downstream of the orifice, be equal to the vapour pressure liquid. In fact, it is well above the vapour pressure at the point of inception. The cavitation number [σ = (P₃ − P_v)/(0.5 pV²); here P₃ is the downstream pressure, P_v is the vapour pressure of the liquid, ρ is the density of the liquid and V is the average liquid velocity at the orifice] at inception is independent of the liquid velocity but strongly dependent on the size of the geometry. Choked cavitation occurs when this minimum pressure approaches the vapour pressure. The cavitation number at the choked condition is a function of the ratio of the orifice diameter (d) to the pipe diameter (D) only. When super cavitation occurs, the dimensionless jet length [L/(D - d); where L is the dimensional length of the jet] can be correlated by using the cavitation number. The vaporization rate of the surface of the liquid jet in super cavitation has been evaluated based on the experiments.

Experiments have also been conducted in which air was deliberately introduced at the vena contracta to simulate the flow regime transition at choked cavitation. Correlations have been obtained to calculate the critical air flow rate required to cause the flow regime transition. By drawing an analogy with choked cavitation, where the air flow rate required to cause the transition is zero, the vapour and released gas flow rate can be predicted.     

Stress state of compound polygonal plate

The constructions made of bars and plates with holes, openings and bulges of various forms are widely used in modern industry. By loading these structural elements with different efforts, there appears concentration (accumulation) of stress whose values sometimes exceeds the admissible one. The durability of the given element is defined according to the quantity of these stresses. Since the failure of details and construction itself begins from the place where the stress concentration has the greatest value.

Therefore the exact determination of stress distribution in details (bars, plates, beams) is of great scientific and practical interest and is one of the important problems of the solid fracture.

Compound details (when the nucleus of different material is soldered to the hole) are often used to decrease the stress concentration.

In the present paper, we study a stress–strain state of polygonal plate weakened by a central elliptic hole with two linear cracks info which a rigid nucleus (elliptic cylinder with two linear bulges) of different material was put in (soldered) without preload.

The problem is solved by a complex variable functions theory stated in papers [Theory of Elasticity, Higher School, Moscow, 1976, p. 276; Plane Problem of Elasticity Theory of Plates with Holes, Cuts and Inclusions, Publishing House Highest School, Kiev, 1975, p. 228; Bidimensional Problem of Elasticity Theory, Stroyizdat, Moscow, 1991, p. 352; Science, Moscow (1996) 708; MSB AH USSR OTH 9 (1948) 1371].

Kolosov–Mushkelishvili complex potential (z) and ψ(z) satisfying the definite boundary conditions are sought in the form of sums of functional series.

After making several strict mathematical transformations, the problem is reduced to the solution of a system of linear algebraic equations with respect to the coefficients of expansions of functions (z) and ψ(z).

Determining the values of (z) and ψ(z), we can find the stress components σ_r, σ_θ and τ_rθ at any point of cross-section of the plate and nucleus on the basis of the known formulae. The obtained solution is illustrated by numerical example.

Changing the parameters A₁, m₁, e, A₂, and m₂ we can get the various contour plates.

For example, if we assume m₁=0, A₁=r, then the internal contour of L₁ becomes the circle of radius r with two rectilinear cracks (for the nucleus––a rectilinear bulges).

Further, if we assume a small semi-axis of the ellipse b₁ to be equal to zero (b₁=0), then a linear crack becomes the internal contour of L₁ (and the nucleus becomes the linear rigid inclusion made of other material). For m₂=0; A₂=R, the external contour L₂ turns into the circle of radius R.

The obtained method of solution may be applied and in other similar problems of elasticity theory; tension of compound polygonal plate, torsion and bending of compound prismatic beams, etc.     

Shear-free turbulent diffusion—classical and new scaling laws

We consider the problem of turbulence generation at a vibrating grid in the x₂–x₃ plane. Turbulence diffuses in the x₁-direction. Analyzing the multi-point correlation equation using Lie-group analysis, we find three different invariant solutions (scaling laws): classical diffusion-like solution (heat equation like), decelerating diffusion-wave solution and finite domain diffusion due to rotation. All solutions have been obtained using Lie-group (symmetry) methods. It is shown that if only one spatial dimension is considered, models based on Reynolds averaging are only capable to model either the diffusion-like solution or the decelerating diffusion-wave solution. The latter solution is only admitted under certain algebraic constraints on the model constants; e.g. in case of the K– model the model constants need to obey the relation c₂σ/σ_K=2. Turbulent diffusion on a finite domain induced by rotation is not admitted by any of the classical models. Finally, in the appendix it is shown that Lele's transformation (Phys. Fluids 28(1) (1985) 64) leads to a complete analytic solution of the steady diffusion problem modelled by the K– equation.     

A study of properties of vortex stretching and enstrophy generation in numerical and laboratory turbulence

Dynamically relevant alignments are used in order to show that regions with weak vorticity are not structureless, non-Gaussian and dynamically not passive. for example, the structure of vorticity in quasi-homogeneous/isotropic turbulent flows is associated with strong alignment between vorticity ω and the eigenvectors of the rate of strain tensor λ_i (especially — but not only — between ω and λ₂) rather than with intense vorticity only. Consequently, much larger regions of turbulent flow than just those with intense vorticity are spatially structured. The whole flow field — even with the weakest measurable enstrophy — is strongly non-Gaussian, which among other things is manifested in strong alignment between vorticity and the vortex stretching vector W_i ≡ ω_jS_ij. It is shown that the quasi-two-dimensional regions corresponding to large cos(ω, λ₂) are qualitatively different from purely two-dimensional ones, e.g. in that they possess essentially nonvanishing enstrophy generation, which is larger than its mean for the whole field.     

Electromagnetic propagation in periodic porous structures

A variational technique is employed to compute approximate propagation constants for electromagnetic waves in a dielectric structure which is periodic in the X−Y plane and translationally invariant in the Z-direction. The fundamental cell, in the periodic structure, is composed of a pore and the surrounding host media. The pore is a circle of radius R₀ filled with a dielectric ε₁ and the host dielectric characterized by ε₂. The size of the cell is characterized by the length A which is R₀.

Two limiting cases are considered. In the first, the pore size is assumed to be much smaller than the wavelength; this limit is motivated by microwave heating of porous material. The approximate propagation constants are explicitly computed for this case and are shown to depend upon the two dielectric constants, the relative areas of the two regions in the cell, and on a modal number. They are not given by a simple mixture formula.

In the second limit, the pore size is taken to be of the same order as the wavelength; this limit is motivated by the propagation of light in a holey fiber. In this case our argument directly yields the dispersion relationship recently derived by Ferrando et al. [Opt. Lett. 24 (1999) 276], using intuitive and physical reasoning. Thus, our method puts theirs into a mathematical framework from which other approximations might be deduced.     

An application of Nagumo''s lemma to some singularly perturbed systems

The existence and asymptotic behavior as ε → 0⁺ of periodic, almost periodic, and bounded solutions of the differential system x = f(t, x, y, ε), Ωy′ = g(t, x, y, ε), are considered where x, f; are n-vectors, y, g are m-vectors and Ω = diag{ε^h₁}…, ε^h_m for integral h_i, h₁ h₂ …, h_m. The principal tools are a lemma of Nagumo which allows the construction of appropriate upper and lower solutions and the asymptotic theory of singularly perturbed linear differential systems.     

Effusing core at the center of A potential vortex

Experiments were carried out to measure the base pressure distribution of a flow field induced by a potential vortex with its axis normal to a stationary disk. The center base pressure coefficient of the vortex, C₀(0), was found to be proportional to Reynolds number from Re = 2.0 × 10³ to Re > 2.5 × 10⁴, where Re is based on the disk radius and azimuthal velocity at the disk edge. This behavior of C₀(0) is at variance with the experimental results of Phillips (Phys. Fluids, 27, 2215, 1984) and Khoo (M. Eng. Thesis, Natl. Univ. Singapore, 1984), which showed vastly different trends depending on Re. Plausible reasons are suggested for the apparent discrepancies observed. Finally, the extent of the effusing core at the center, r₁ (taken to be the radial position where departure from the outer potential flow took place), was found to be proportional to Re^−1/2 for all Re values considered.     

Traveling disturbance appearing in boundary layer transition in a Yawed cylinder

Boundary layers that develop over a body in fluid flow are in most cases three-dimensional owing to the spin, yaw, or surface curvature of the body. Therefore, the study of three-dimensional (3D) boundary-layer transition is essential to work in practical aerodynamics. The present investigation is concerned with the problem of 3D boundary layers over a yawed body. A yawed cylinder model that represents the leading edge portion of a swept wing and the mechanism of crossflow instability are investigated in detail using hot-wire velocimetry and a flow visualization technique. As a result, traveling disturbances having frequencies f₁ and f₂, which differ by about one order of magnitude, are detected in the transition region. The phase velocities and directions of travel of those disturbances are measured. Results for the low-frequency disturbance f₁ show qualitative coincidence with results numerically predicted for a crossflow unsteady disturbance. Nameley, F₁ travels nearly spanwise to the yawed cylinder and very close to the cylinder wall. The results for the high-frequency disturbance f₂ good agreement with the existing experimental results. The ₂ disturbance is found to be the high-frequency inflectional secondary instability that appears in 3D boundary layer transition in general. A two-stage transition process, where stationary crossflow vortices appear as the primary instability and a traveling inflectional disturbance is generated as a secondary instability, was observed. Secondary instability seems to play a major role in turbulent transition.     

The interaction between a plane shear layer and a slender body

The interaction of a plane shear layer with a thin flat plate located in the nonlinear region of the shear layer has been investigated experimentally. The shear layer's velocity ratio is 0.375, and its Reynolds number is ΔUθ₀/μ = 625. It is found that for different angles of attack of the plate, the mixing layer is deflected toward the slower stream. In addition, the plate attenuates turbulent fluctuations of the shear layer structures.     

Flow measurements inside a heated multiple rotating cavity with axial throughflow

This paper discusses experimental results from a multiple cavity test rig representative of a high pressure compressor internal air system. Measurements of the axial, tangential and radial velocity components are presented. These were made using a two component, laser doppler anemometry (LDA) system for a range of non-dimensional parameters representative of engine conditions (Re up to 4 × 10⁶ and Re_z up to 1.8 × 10⁵). Tests were carried out for two different sizes of annular gap between the (non-rotating) drive shaft and the disc bores.

The axial and radial velocities inside the cavities are virtually zero. The size of the annular gap between disc bore and shaft has a significant effect on the radial distribution of tangential velocity. For the narrow annular gap (d_h/b = 0.092), there is an increase of non-dimensional tangential velocity V/Ωr with radial location from V/Ωr < 1 at the lower radii to solid body rotation V/Ωr = 1 further into the cavity. For the wider annular gap (d_h/b = 0.164), there is a decrease from V/Ωr > 1 at the lower radii to solid body rotation further into the cavity. An analysis of the frequency spectrum obtained from the tangential velocity measurements is consistent with a flow structure in the r– plane consisting of pairs of contra rotating vortices.     

Sound wave scattering by circular cylindrical shells

The symmetric (S₀) and antisymmetric (A₀) Lamb-type waves generated in a thin elastic cylindrical shell by normal incidence of an acoustic wave have been considered. The typical frequency dependencies (FD) of the backscattered acoustic pressure at on observetion point in the far field are presented. The spectra of the S₀ and A₀ waves are marked on them. It was found that if the A₀ wave is excited in the shell, then its phase velocity is greater than the sound velocity in the fluid surrounding the shell. The parameter which defines the center of the strong bending domain (SBD) is defined. It is shown that in this domain the A₀ wave is practically non-dispersive. Phase velocity data for the A₀ wave are given. Spectra and dispersion curves of the S₀ wave for shells which have different relative thickness s and which are made of different materials have been examined.     

圆盘近尾迹流场相关谱分析

对Re=22 000 时径厚比D/h=5 圆盘近尾迹开展大涡模拟数值研究. 通过对x/D=1, 2 和8 处脉动速度进行快速傅里叶变换（fast Fourier transformation, FFT）,发现3 个特征频率：斯特劳哈尔数St₂=0.123 为自然脱落频率,与文献结果相符;较小频率St₁=0.035,与回流区伸缩和剪切层附近涡旋脱落点的旋转相关;高频率St₃=1.3～1.7 则与剪切层湍流结构相关. 通过分析截面r/D=2.8 圆周上两点间流向速度相关系数、相干谱和相位谱,发现相关系数受涡旋脱落影响出现以30°或45°为周期的正负波动,表明轴面上涡旋脱落点具有随机性.     

Experimental and numerical determinations of the initial surface of phase transformation under biaxial loading in some polycrystalline shape-memory alloys

Biaxial proportional loading such as tension (compression)–internal pressure and bi-compression tests are performed on a Cu-Zn-Al and Cu-Al-Be shape memory polycrystals. These tests lead to the experimental determination of the initial surface of phase transformation (austenite→martensite) in the principal stress space (σ₁,σ₂). A first “micro–macro” modeling is performed as follows. Lattice measurements of the cubic austenite and the monoclinic martensite cells are used to determine the “nature” of the phase transformation, i.e. an exact interface between the parent phase and an untwinned martensite variant. The yield surface is obtained by a simple (Sachs constant stress) averaging procedure assuming random texture. A second modeling, performed in the context of the thermodynamics of irreversible processes, consists of a phenomenological approach at the scale of the polycrystal. These two models fit the experimental phase transformation surface well.     

On the dispersion of a solute in Poiseuille flow of a third-grade fluid in a channel

Gill and Sankarasubramanian's analysis of the dispersion of Newtonian fluids in laminar flow between two parallel walls are extended to the flow of non-Newtonian viscoelastic fluid (known as third-grade fluid) using a generalized dispersion model which is valid for all times after the solute injection. The exact expression is obtained for longitudinal convective coefficient K₁(Γ), which shows the effect of the added viscosity coefficient Γ on the convective coefficient. It is seen that the value of the K₁(Γ) for Γ≠0 is always smaller than the corresponding value for a Newtonian fluid. Also, the effect of the added viscosity coefficient on the K₂(t,Γ) (which is a measure of the longitudinal dispersion coefficient of the solute) is explored numerically. Finally, the axial distribution of the average concentration C_m of the solute over the channel cross-section is determined at a fixed instant after the solute injection for several values of the added viscosity coefficient.     

Effect of tube diameter on Preston tube calibration curves for the measurement of wall shear stress

The effect of tube diameter (d) on Preston tube calibration curves for the measurement of wall shear stress (τ_w) in a zero pressure gradient turbulent boundary layer has been investigated. Five different outside diameter tubes of 1.46, 1.82, 3.23, 4.76 and 5.54 mm, corresponding to (d/δ) of 0.022, 0.027, 0.048, 0.071 and 0.082 were used to measure τ_w at Reynolds numbers based on momentum thickness (R_θ) of 2800–4100. The calibration curves of Patel (V.C. Patel, J. Fluid Mech. 23 (part I) (1965) 185–208) and Bechert (D.W. Bechert, AIAA J. 34 (1) (1995) 205–206) are both dependent on the tube diameter. The maximum difference in the τ_w measurements from the different tubes using Patel's calibration is about 8%, while Bechert's calibration gives a maximum difference of approximately 18%.     

仿生波纹夹层结构耐撞性分析及优化

为提高薄壁夹层结构耐撞性,以虾螯为仿生原型,设计梯度分布的仿生波纹形夹层结构,包括单层、双层和三层波纹结构。以初始峰值载荷F_p、比吸能E_s为耐撞性指标,利用有限元法分析了单元高宽比γ（γ₁、γ₂和γ₃分别为单元第1层、第2层和第3层的高宽比）对波纹夹层结构耐撞性的影响,采用多目标粒子群优化方法得到了夹层结构最优参数。结果表明,单层波纹结构耐撞性随单元高宽比γ的增大逐渐变差,双层波纹结构下层结构单元高宽比γ对耐撞性的影响大于上层结构单元高宽比γ对耐撞性的影响,较小的γ值有利于提高三层波纹结构的比吸能。结构优化结果表明:单层结构最优尺寸γ₁为0.8;双层结构最优尺寸为γ₁ = 0.5和γ₂ = 1.2;三层结构最优组合为γ₁ = 0.6,γ₂ = 0.6和γ₃ = 0.9。上述结果可为薄壁夹层结构轻量化设计提供新思路。