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.     

Solutions of some classes of second order non-linear ordinary differential equations

A second order non-linear ordinary differential equation satisfied by a homogeneous function of u and v where u is a solution of the linear equation ÿ + p(t)ÿ + r(t)y = 0 and v = ωu, ω being an arbitrary function of t, is obtained. Defining ω suitably in two specific cases, solutions are obtained for a non-linear equation of the form ÿ + p(t)ÿ + q(t)y = μÿ²y^{−1 +} f(t)yⁿ where μ ≠ 1, n≠ 1. Applying our results, some classes of equations of the above type possessing solutions involving two or one or no arbitrary constants are derived. Some illustrative examples are also discussed.     

An unsymmetrical van der pol equation with stable subharmonics

Possible stable subharmonic solutions of the equation

ÿ − k(1 + 2cy − y²)ÿ + Y = bkμ cos μt, c > 0

, klarge, are discussed by the techniques used by J.E. Littlewood for van der Pol's equation in Acta Math. 97 (1957), that is the case of the above equation with c = 0 and

, k large. Their variation as c increases is also considered briefly.     

Experimental study of turbulence in isothermal and stably stratified homogeneous shear flows: Mean flow measurements

The evolution of freestream turbulence under the combined action of linear shear and stable linear temperature profile is investigated. The experiment is carried out in a small, open circuit, low-speed test cell that uses air as working fluid. The temperature gradient formed at the entrance to the test section by means of an array of 24 horizontal, differentially heated elements is varied to get a maximum Brunt-Vaisala frequency N_o[=({g/T_m}{∂T/∂y})^1/2] of 3.1⁻¹. Linear velocity profiles are produced using screens of variable mesh size. The Reynolds number Re_M based on centre-line velocity and mesh size is varied from 80 to 175. Isothermal studies are carried out in four different experiments with varying velocity gradients. The effect of inlet turbulence level on growth of turbulence is studied in these flows by keeping the shear parameter Sh (=(x/u)(∂u/∂y)) constant. The range of shear parameters considered is 2.5–7.0. Shear and stratification combined produce a maximum gradient Richardson number Ri_g (= N_o²/(∂u/∂y)²) of 0.0145. Results have been presented in terms of evolution of variance of velocity fluctuations, Reynolds shear stress and temperature fluctuations. Measurements show the following: In isothermal flows the growth rate of turbulence quantities depends on both shear parameter and inlet turbulence level. There are distinct stages in the evolution of the flow and that can be identified by the power-law exponent of growth of turbulence. Shear is seen to promote the growth of turbulence and accelerate it towards a fully developed equilibrium state. Stratification initially suppresses the growth of turbulence and, hence, enhances the degree of underdevelopment. Under these conditions shear becomes active and subsequently enhances the growth rate of turbulence quantities.     

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.     

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.     

Non-linear oscillations with multiple forcing terms

The problem of determining the transient response of a non-linear oscillator of the form ü + u = εƒ(u,u) + E(t) is studied by the method of multiple time scales, using the symbolic computation system MACSYMA. when the excitation E(t) consists of a finite number of harmonic forcing terms. Here ε is a small parameter and ƒ(u,u) is a non-linear function of its arguments. In particular, the Van der Pol and Duffing oscillators are studied in detail. It is found that when the forcing frequencies are not close to each other or close to the primary resonance of the system, then the response of the system is analogous to the behavior when only one forcing term is present. However, when the forcing frequencies are close to each other or close to the primary resonance, then the behavior is quite different, exhibiting certain oscillations not observed in the case of one forcing term.     

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.     

Instantaneous fluctuation velocity and skewness distributions upstream of transition onset

The development of streamwise orientated disturbances through the boundary layer thickness prior to transition onset for zero-pressure gradient boundary layer flow under the influence %Tu = 4.2 is presented. The analysis concentrates on the development of the maximum positive and negative of the fluctuation velocity in order to gain further insight into the transition process. The average location of the peak negative fluctuation velocity over a range of Reynolds numbers was measured in the upper portion of the boundary layer at y/δ ≈ 0.6, whereas the location of the peak positive value was measured at y/δ ≈ 0.3. The disturbance magnitude of the negative fluctuation velocity increased beyond that of the positive as transition onset approached. The distribution and disturbance magnitude of the maximum positive and negative fluctuation velocities indicate that the initiation of transition may occur on the low-speed components of the flow that are lifted up to the upper region of the boundary layer. This is in qualitative agreement with recent direct numerical simulations on the breakdown of the flow on the lifted low-speed streaks near the boundary layer edge. The results presented in this investigation also demonstrate the increased physical insight gained by examining the distributions of the maximum positive and negative of the streamwise fluctuation velocity component associated with the low- and high-speed streaks, compared to time-averaged values, in determining what structures cause the breakdown to turbulence.     

Fluid flow in a rotating cylindrical container with a rotating disk at the fluid surface

Fluid flow in a rotating cylindrical container of radius R_w and height H with a co-axially rotating disk of radius R_d at the fluid surface is numerically investigated. The container and the disk rotate with angular velocities Ω_w and Ω_d, respectively. We solve the axisymmetric Navier-Stokes equations using a finite-volume method. The effects of the relative directions and magnitudes of the disk and container rotations are studied. The calculations are carried out with various ratios of Ω_w and Ω_d for H/R_w = 2 and R_d/R_w = 0.7. Streamlines and velocity vectors in the meridional plane and azimuthal velocities are obtained. The flow fields in the meridional plane are discussed with relation to azimuthal velocities in the interior of the container. The numerical results are also compared with experimental data.     

Hierarchical structures in a turbulent pipe flow

A hierarchical structure (HS) analysis (β-test and γ-test) is applied to a fully developed turbulent pipe flow. Velocity signals are measured at two cross sections in the pipe and at a series of radial locations from the pipe wall. Particular attention is paid to the variation of turbulent statistics at wall units 10<y⁺<3000. It is shown that at all locations the velocity fluctuations satisfy the She–Leveque hierarchical symmetry (Phys. Rev. Lett. 72 (1994) 336). The measured HS parameters, β and γ, are interpreted in terms of the variation of fluid structures. Intense anisotropic fluid structures generated near the wall appear to be more singular than the most intermittent structures in isotropic turbulence and appear to be more outstanding compared to the background fluctuations; this yields a more intermittent velocity signal with smaller γ and β. As turbulence migrates into the logarithmic region, small-scale motions are generated by an energy cascade and large-scale organized structures emerge which are also less singular than the most intermittent structures of isotropic turbulence. At the center, turbulence is nearly isotropic, and β and γ are close to the 1994 She–Leveque predictions. A transition is observed from the logarithmic region to the center in which γ drops and the large-scale organized structures break down. We speculate that it is due to the growing eddy viscosity effects of widely spread turbulent fluctuations in a similar way as in the breakdown of the Taylor vortices in a turbulent Couette–Taylor flow at high Reynolds numbers.     

The effect of the Reynolds number on the Reynolds stresses and vorticity in a turbulent far-wake

Using two orthogonal arrays of 16 X-wires, eight in the (x,y)-plane and eight in the (x,z)-plane, the effect of the Reynolds number in a turbulent plane far-wake has been investigated for two values of Re_θ (based on the free stream velocity and the momentum thickness), i.e. 1350 and 4600. It is observed that as the Reynolds number increases the magnitudes of the measured Reynolds stresses increase, as does the size of two-point vorticity correlation iso-contours. Discernible differences are also observed in probability density function, spectra and three-dimensional topologies. The Reynolds number dependence seems to vanish when Re_θ5000.     

Correlations predicting frictional pressure drop and liquid holdup during horizontal gas-liquid pipe flow with a small liquid holdup

Experimental data and correlations available in the literature for the liquid holdup ε_L and the pressure gradient ΔP_TP/L for gas-liquid pipe flow, generally, do not cover the domain 0 < ε_L < 0.06. Reliable pressure-drop correlations for this holdup range are important for calculating flow rates of natural gas, containing traces of condensate. In the present paper attention is focused on reliable measurements of ε_L and ΔP_TPIL values and on the development of a phenomenological model for the liquid-holdup range 0 < ε_L < 0.06. This model is called the “apparent rough surface” model and is referred to as the ARS model. The experimental results presented in this paper refer to air-water and air-water + ethyleneglycol systems with varying transport properties in horizontal straight smooth glass tubes under steady-state conditions. The holdup and pressure gradient values predicted with the ARS model agree satisfactorily with both our experimental results and data obtained from the literature referring to small liquid-holdup values 0 < ε_L < 0.06. Further, it has been shown that in the domain 38 < < 72 mPa m the interfacial tension of the gas-liquid system has no significant effect on the liquid holdup. The pressure gradient, however, increases slightly with decreasing surface tension values.     

Automatic control system of a rear-wheel drive vehicle moving on a sloped weak sandy terrain

The general mechanism of tractive performance of a four-wheel vehicle with rear-wheel drive moving up and down a sloped sandy soil has been considered theoretically. For the given vehicle dimensions and terrain-wheel system constants, the relationships among the effective tractive or braking effort of the vehicle, the amount of sinkage of the front and rear wheels, and the slip ratio were analysed by simulation. The optimum eccentricity of the vehicle’s center of gravity and the optimum application height of the drawbar-pull for obtaining the largest value of maximum effective tractive or braking effort could be calculated by means of the analytical simulation program. For a 5.88 kN weight vehicle, it was found that the optimum eccentricity of the center of gravity e_opt was 1/6 for the range of slope angle—0βπ/24 rad during driving action of the rear wheel and e_opt was also 1/6 for the range of slope angle—π/24β0 rad during braking action of the rear wheel. The optimum application height H_opt was found to be 35 cm for the range of slope angle 0βπ/24 rad during driving action of the rear wheel and H_opt was 0 cm for the range of slope angle—π/24β0 rad during braking action of the rear wheel.     

Hole diameter effect on flow characteristics of wake behind porous fences having the same porosity

The hole diameter effect on the flow characteristics of wake behind porous fences has been investigated experimentally in a circulating water channel having a test section of 300^w×200^h×1200^l (mm). Three porous fences having different hole diameters of d=1.4,2.1,2.8 mm were tested in this study, but they have the same =38.5% geometric porosity. One thousand instantaneous velocity fields for each fence were measured consecutively by the hybrid PTV system employing a high-speed CCD camera. Free stream velocity was fixed at 10 cm/sec and the corresponding Reynolds number based on the fence height was Re=2,985. Consequently, the fence with the smallest hole diameter d=1.4 mm (d_1.4) decreases the streamwise velocity component and increases the vertical velocity component. Among the three hole diameters tested in this study, the d_1.4 fence has the largest turbulence intensity in the shear layer developed from the fence top. Regardless of the hole diameter, however, all three fences having the same porosity reduce the reduction of turbulent intensity in the lower region below the fence height (y/H<1).     

Experimental study of heat transfer characteristics due to confined impinging two-dimensional jets

Experimental results are presented for characteristics of impingement heat transfer caused by three slot jets. Experimental values were obtained for the dimensionless distance H = 0.5−3, dimensionless pitch P = 6−16, and Reynolds number Re = 500−8000. For laminar impinging flow, they were compared with numerical results. For turbulent impinging flow, two peaks of the local Nusselt number were obtained behind the second nozzle. The position of the second peak approached the nozzle as the space between nozzle and impinged surface decreased. The average Nusselt number between the central and second nozzles was determined from the ratio P/H and the Reynolds number based on the pitch of the nozzles.     

The instability of rhombic cell flows

Instability of two-dimensional periodic flows with rhombic cell structure represented by the stream function Ψ=cos kx+cosy is investigated. Stability characteristics are obtained for the Reynolds number R=1, 2, 3 and 4 and the ratio of the diagonals of the cell . Variation of the critical Reynolds number R_c with k is obtained, and the square cell flow (k=1) is found to be most stable (R_c=√2). It is found that R_c → 1 as k → 0, which leads to a finite gap between this limiting R_c and R_c=√2 for K=0 (Ψ=cos y).     

Equations differentielles non lineaires: Differentiabilite des solutions par rapport aux conditions initiales et problemes de la stabilite

In this paper, we examine the problem of the asymptotic stability of the solutions of a differential equation (E), Y = X(Y, t), from the point of view of the differentiability of the solutions with respect to the initial conditions. The method allows us to deal with cases in which the mapping X is not differentiable everywhere, and in which the variational equation of (E) is not defined in the usual meaning. This method can be carried on for the research of periodic solutions. We give two examples.

Zusammenfassung

In diesem Bericht wird das Problem der asymptotischen Stabilität der Lösungen einer Differential Gleichung (E), Y = X(Y, t), vom Standpunkt der Differenzierbarkeit der Lösungen in Bezug auf die Anfangsbedingungen untersucht. Das Verfahren erlaubt die Fälle zu behandeln, worin die Abbildung X nicht überall differenzierbar ist, und worin auch die Variations Gleichung nicht im üblichen Sinne definiert ist. Das Verfahren kann erweitert werden zur Untersuchung periodischer Lösungen. Zwei Beispiele werden vorgeführt.

Abstract

B aдннoй paбaтe иccлeдyeтcя вoпoc oб accимeтpиХecкoй ycтoйХибocти peшeний диффepeHциaлънoгo ypaвнeния (E), Y = X(Y, t), c тoХки зpeния диффepeнциpyeмocти peшeний пo oтнoшeнию к нaХaлъным ycлoвиям. Иcпoлъзyeмый мeтoд дaeт вoзмoжнocтъ иccлeдoвaтъ cлyХaи, кoгдa пpeoбpaзoвaниe X нe вcюдy диффepeнциaлнoe. Bapиaциoннoe ypaвнeниe для oпpeдeляeтcя нecкoлъкo инaХe Хeм oбыХнo. Дaнный мeтoд мoжeт иcпoлъзoвaтъcя для иccлeдoвa*ncy;ия пepиoдиХecкиx peшeний. Дaютcя двa пpимepa.     

On one-dimensional shock waves

In this paper uniform asymptotic expansions for the solutions of a system of differential equations are obtained in the domain containing a shock wave. It is shown, in particular, that the function θ(t,x)/ε contained in the expansions and describing the behavior of the solution in the neighborhood of the wave front has, generally speaking, a discontinuity of derivatives at the front. The results are applicable to one-dimensional problems in gas dynamics with low viscosity and heat-conductivity.     

Numerical prediction of locally forced turbulent separated and reattaching flow

An unsteady numerical simulation was performed for locally forced separated and reattaching flow over a backward-facing step. The local forcing was given to the separated and reattaching flow by means of a sinusoidally oscillating jet from a separation line. A version of the k––f_μ model was employed, in which the near-wall behavior without reference to distance and the nonequilibrium effect in the recirculation region were incorporated. The Reynolds number based on the step height (H) was fixed at Re_H=33 000, and the forcing frequency was varied in the range 0St_H2. The predicted results were compared and validated with the experimental data of Chun and Chun. It was shown that the unsteady locally forced separated and reattaching flows are predicted reasonably well with the k––f_μ model. To characterize the large-scale vortex evolution due to the local forcing, numerical flow visualizations were carried out.