Influence of unsteady wake on a turbulent separation bubble

 An experimental study was made of turbulent separated and reattaching flow over a blunt body, where unsteady wake was generated by a spoked-wheel type of wake generator with cylindrical rods. The influence of unsteady wake was scrutinized by altering the rotation direction (clockwise and counter-clockwise) and the normalized passing frequency (0 ≤ St _H  ≤ 0.20). The Reynolds number based on the cylindrical rod was Re _d =375. A phase-averaging technique was employed to characterize the unsteady wake. The effect of different rotation directions, which gave a significant reduction of x _R , was examined in detail. The wall pressure fluctuations on the blunt body were analyzed in terms of the spectrum and the coherence. Received: 15 January 2001 / Accepted: 17 July 2001     

Heat transfer by laminar flow in a vertical and horizontal pipe with twisted-tape inserts

 This article presents the results of laboratory research on heat exchange while heating water in horizontal and vertical tubes with twisted-tape inserts. The scope of the research: 70 ≤ Re ≤ 4000 3.6 ≤ Pr ≤ 5.9 8.6 ≤ Gz ≤ 540 The research was held for three cases: – horizontal experimental tube – vertical experimental tube, the direction of flow according to the free convection vector – vertical experimental tube, the direction of flow not in accordance with the free convection vector For such cases the correlation equation was defined Nu_T=f(Gz; y), Nu = f(Gz) and the proportion Nu_T/Nu was analysed. Received on 30 March 2000     

Second-order nonlinear differential equations with infinite set of periodic solutions

For the differential equation u″ = f(t, u, u′), where the function f: R × R ² → R is periodic in the first variable and f (t, x, 0) ≡ 0, sufficient conditions for the existence of a continuum of nonconstant periodic solutions are found. Published in Neliniini Kolyvannya, Vol. 11, No. 4, pp. 495–500, October–December, 2008.     

Experimental investigation of mixed convection heat transfer in a horizontal and inclined rectangular channel

 Mixed convection heat transfer in rectangular channels has been investigated experimentally under various operating conditions. The lower surface of the channel is subjected to a uniform heat flux, sidewalls are insulated and adiabatic, and the upper surface is exposed to the surrounding fluid. Experiments were conducted for Pr=0.7, aspect ratios AR=5 and 10, inclination angles 0° ≤ θ ≤ 30°, Reynolds numbers 50 ≤ Re ≤ 1000, and modified Grashof numbers Gr*=7.0 × 10⁵ to 4.0 × 10⁷. From the parametric study, local Nusselt number distributions were obtained and effects of channel inclination, surface heat flux and Reynolds number on the onset of instability were investigated. Results related to the buoyancy affected secondary flow and the onset of instability have been discussed. Some of the results obtained from the experimental measurements are also compared with the literature, and a good agreement was observed. The onset of instability was found to move upstream for increasing Grashof number and increasing aspect ratio. On the other hand, onset of instability was delayed for increasing Reynolds number and increasing inclination angle. Received on 19 March 2001 / Published online: 29 November 2001     

Heat transfer by laminar flow in a vertical pipe with twisted-tape inserts

Heat transfer for laminar flow of water in an air-cooled vertical copper pipe with four twisted-tape inserts was determined experimentally. The tests were executed for laminar flow within 110 ≤ Re ≤ 1500, 8.1 ≤ Gz ≤  82.0 and 1.62 ≤ y ≤ 5.29. The correlation equation for heat transfer was defined for the tested range. The obtained results were compared to the results of other authors. Received on 28 April 1998     

Rheological and mechanical properties of silica colloids: from Newtonian liquid to brittle behaviour

Rheological and mechanical properties of aqueous mono-disperse silica suspensions (Ludox? HS40) are investigated as a function of particle volume fraction (ϕ _p ranging from 0.22 to 0.51) and water content, using shear rate tests, oscillatory methods, indentation and an ultrasonic technique. As the samples are progressively dried, four regimes are identified; they are related to the increasing particle content and the existence and behaviour of the electrical double layer (EDL) around each particle. For 0.22 ≤ ϕ _p ≤ 0.30), the suspensions are stable due to the strong electrostatic repulsion between particles and show Newtonian behaviour (I). As water is removed, the solution pH decreases and the ionic strength increases. The EDL thickness therefore slowly decreases, and screening of the electrostatic repulsion increases. For 0.31 ≤ ϕ _p ≤ 0.35, the suspensions become turbid and exhibit viscoelastic (VE) shear thinning behaviour (II), as they progressively flocculate. For 0.35 ≤ ϕ _p ≤ 0.47, the suspensions turn transparent again and paste-like, with VE shear thinning behaviour and high elastic modulus (III). At higher particle concentration, the suspensions undergo a glass transition and behave as an elastic brittle solid (IV, ϕ _p = 0.51).     

Shear-strain-rate effects in a high-strength aluminum alloy

Methods of investigating the plastic flow of materials at high rates are reviewed, and experiments are described in which thin-walled tubular specimens were subjected to dynamic torsional loading. These experiments were performed using a modified version of a torsional Hopkinsonbar apparatus used in earlier work. The results show that, at strain rates of order 10³ sec⁻¹, the stress-strain curve of the alloy tested does not differ significantly from that found at 10⁻³ sec⁻¹. In tests involving the propagation of a stress increment, however, it was found that the speed of propagation was that of elastic shear waves, indicating that the initial response of the material is essentially rate dependent.     

Effect of an Exit-Wedge Angle on Pinch-off and Mass Entrainment of Vortex Rings in Air

The effects of exit-wedge angle on evolution, formation, pinch-off, propagation and diffusive mass entrainment of vortex rings in air were studied using digital particle image velocimetry. Vortex rings were generated by passing a solenoid-valve-controlled air jet through a cylindrical nozzle. Experiments were performed over a wide range of exit-wedge angles (10° ≤ α ≤ 90°) of the cylindrical nozzle, initial Reynolds numbers (450 ≤ Re ≤ 4,580) and length-to-diameter ratios (0.9 ≤ L/D ≤ 11) of the air jet. For sharp edges (α ≤ 10°), a secondary ring may emerge at high Reynolds numbers, which tended to distort the vortex ring if ingested into it. For blunt edges (α ≥ 45°), by contrast, stable vortex rings were produced. The formation phase of a vortex ring was found to be closely related to its evolution pattern. An exit-wedge angle of 45° was found to be optimal for rapid pinch-off and faster propagation and better stability of a vortex ring. Diffusive mass entrainment was found to be between 35% and 40% in the early stages of a vortex ring propagation and it gradually increased throughout the course of vortex ring propagation. Entrainment fraction was found to be sensitive to the L/D ratio of the initial jet and decreases when the L/D ratio is increased.     

Rheology of SiO2/(acrylic polymer/epoxy) suspensions. I. Linear viscoelasticity

Linear viscoelastic properties of SiO₂/(AP/EP) suspension with various SiO₂ volume fractions (ϕ) in a blend of acrylic polymer (AP) and epoxy (EP) were investigated at various temperatures (T). The AP/EP contained 70 vol.% of EP. The SiO₂ particles were treated with epoxy silane coupling agent. The effects of the SiO₂ particles are more pronounced in the terminal zone: a transition from viscoelastic liquid (ϕ ≤ 30 vol.%) to viscoelastic solid (ϕ ≥ 40 vol.%) was observed which can be interpreted as a critical gelation occurring at a critical particle content and critical gel temperature. The SiO₂/(AP/EP) systems exhibited a critical gel behavior at ϕ ≅ 35 vol.% and T ≅ 100°C characterized with a power–law relationship between the storage and loss moduli (G ^′ and G ^″) and frequency (ω); G ^′ = G ^″/tan(nπ/2) ∝ ω ⁿ . The critical gel exponent (n) was estimated to be about 0.45. The gelation occurred with increasing T.     

Effect of local forcing on a turbulent boundary layer

An experimental study is performed to analyze flow structures behind local suction and blowing in a flat-plate turbulent boundary layer. The local forcing is given to the boundary layer flow by means of a sinusoidally oscillating jet issuing from a thin spanwise slot at the wall. The Reynolds number based on the momentum thickness is about Re _θ =1700. The effects of local forcing are scrutinized by altering the forcing frequency (0.011 ≤ f ⁺≤ 0.044). The forcing amplitude is fixed at A ₀=0.4. It is found that a small local forcing reduces the skin friction and the skin friction reduction increases with the forcing frequency. A phase-averaging technique is employed to capture the large-scale vortex evolution. An organized spanwise vortical structure is generated by the local forcing. The cross-sectional area of vortex and the time fraction of vortex are examined by changing the forcing frequency. An investigation of the random fluctuation components reveals that turbulent energy is concentrated near the center of vortical structures. Received: 17 March 2000/Accepted: 3 April 2001     

Measurements of heat transfer coefficients between turbulent liquid jets and a radially finned heated disk placed at the bottom of an open vertical circular cavity

 Experiments have been performed to assess the impact of an extended surface on the heat transfer enhancement for axisymmetric, turbulent liquid jet impingement on a heated round disk. The disk, with an array of integral radial fins mounted on its surface, is placed at the bottom of an open vertical circular cavity. Hydrodynamic and heat transfer data were obtained for a dielectric fluorocarbon liquid FC-77. For a fixed circular heater of diameter D=22.23 mm, several geometric parameters were tested: the nozzle diameter (4.42≤d≤9.27 mm), the confining wall diameter of the vertical cavity (22.23≤D _c≤30.16 mm), and the nozzle-to-heater spacing (0.5≤S/d≤5.0). The FC-77 flow rates varied from =0.2 to 11.0 l/min producing Reynolds numbers in the wide interval 700≤Re _d ≤44,000. For d=4.42 mm, the heat transfer response to the separation distance S/d was small but increased gradually with increasing nozzle diameter up to d=9.27 mm. The thermal resistance R _th increased with the confining wall diameter D _c and also with the nozzle diameter d. A minimum value of the thermal resistance of R _th,min=0.4 cm² K/W was attained for a combination of d=4.42 mm, D _c=22.23 mm, S/d=1, and =7.5 l/min. Based on a simplified heat transfer model, reasonable agreement was obtained between measured values of the thermal resistance and the R _th-predictions. The total fin effectiveness ɛ_f was shown to increase with increasing nozzle diameter, but was invariant with the flow rate (or the jet exit velocity). More than a three-fold heat transfer enhancement was realized through the addition of the array of integral radial fins on the heated round disk. Received on 30 August 2000 / Published online: 29 November 2001     

Global stability of a virus dynamics model with Beddington–DeAngelis incidence rate and CTL immune response

In this paper, the global stability of virus dynamics model with Beddington–DeAngelis infection rate and CTL immune response is studied by constructing Lyapunov functions. We derive the basic reproduction number R ₀ and the immune response reproduction number R ⁰ for the virus infection model, and establish that the global dynamics are completely determined by the values of R ₀. We obtain the global stabilities of the disease-free equilibrium E ₀, immune-free equilibrium E ₁ and endemic equilibrium E ^∗ when R ₀≤1, R ₀>1, R ⁰>1, respectively.     

Effect of initial conditions on decaying grid turbulence at low R λ

The transport equations for the second-order velocity structure functions 〈(δu)²〉 and 〈(δq)²〉 are used as a scale-by-scale budget to quantify the effect of initial conditions at low Reynolds numbers, typical of grid turbulence. The validity of these equations is first investigated via hot-wire measurements of velocity and transverse vorticity fluctuations. The transport equation for 〈(δq)²〉 is shown to be balanced at all scales, while anisotropy of the large scales leads to a significant imbalance in the equation for 〈(δu)²〉. The effect of using similarity to evaluate the transport equation is rigorously tested. This approach has the desirable benefit of requiring less extensive measurements to calculate the inhomogeneous term of the transport equation. The similarity form of the 〈(δq)²〉 equation produces nearly identical results as those obtained without the similarity assumption. In the case of the 〈(δu)²〉 equation, the similarity method forces a balance at large separation, although the imbalance due to large scale anisotropy remains. The initial conditions of the turbulence at constant R _M ≃ 10,400 (28≤ R _λ≤ 55) are changed by using three grids of different geometries. Initial conditions affect the shape and magnitude of the second- and third-order structure functions, as well as the anisotropy of the large scales. The effect of initial conditions on the scale-by-scale budget is restricted to the inhomogeneous term of the transport equations, while the dissipation term remains unaffected despite the low R _λ. Scales as small as λ are affected by the changes in initial conditions.     

Miniaturized Compression Test at Very High Strain Rates by Direct Impact

A modified miniaturized version of the Direct Impact Compression Test (DICT) technique is described in this paper. The method permits determination of the rate-sensitive plastic properties of materials up to strain rate ∼10⁵ s⁻¹. Miniaturization of the experimental setup with specimen dimensions: diameter d _S = 2.0 mm and thickness l _S = 1.0 mm, Hopkinson bar diameter 5.2 mm, with application of a novel optical arrangement in measurement of specimen strain, makes possible compression tests at strain rates from ∼10³ s⁻¹ to ∼10⁵ s⁻¹. In order to estimate the rate sensitivity of a low-alloy construction steel, quasi-static, Split Hopkinson Pressure Bar (SHPB) and DICT tests have been performed at room temperature within the rate spectrum ranging from 5*10⁻⁴ s⁻¹ to 5*10⁴ s⁻¹. Adiabatic heating and friction effects are analyzed and the final true stress versus true strain curves at different strain rates are corrected to a constant temperature and zero friction. The results have been analyzed in the form of true stress versus the logarithm of strain rate and they show two regions of a constant rate sensitivity : relatively low up to the strain rate threshold ∼50 s⁻¹, and relatively high above the threshold, up to strain rate ∼4.5*10⁴ s⁻¹.     

Forced convection of turbulent flow in triangular ducts with different angles and surface roughnesses

The experimental investigations were consisting of two parts. The first part was carried out to study the effect of corner geometry on the steady-state forced convection inside horizontal isosceles triangular ducts with sharp corners. The electrically-heated triangular duct was used to simulate the triangular passage of a plate-fin compact heat exchanger. The isosceles triangular ducts were manufactured with duralumin, and fabricated with the same length of 2.4 m and hydraulic diameter of 0.44 m, but five different apex angles (i.e. θ _a =15^∘,30^∘, 40^∘,60^∘, and 90^∘) respectively. The investigation was performed under turbulent flow condition covering a wide range of Reynolds number (i.e. 7000≤Re _D ≤20000). It was found that the best thermal performance is achieved with the apex angle of 60^∘. The second part was performed to investigate the effect of surface roughness on the forced convection of the same system. Horizontal equilateral triangular ducts with an apex angle of 60^∘ were fabricated with the same length and hydraulic diameter, but different average surface roughnesses of 1.2 m,3.0 m and 11.5 m respectively. It was concluded that the duct with a higher surface roughness will have a better heat transfer performance. Non-dimensional expressions for the determination of the heat transfer coefficient of the triangular ducts with different apex angles and surface roughnesses were also developed. Received on 15 December 1997     

Non-Darcy buoyancy driven flows in a fluid saturated porous medium: the use of asymptotic computational fluid dynamics (ACFD) approach

Natural convection in a fluid saturated porous medium has been numerically investigated using a generalized non-Darcy approach. The governing equations are solved by using Finite Volume approach. First order upwind scheme is employed for convective formulation and SIMPLE algorithm for pressure velocity coupling. Numerical results are presented to study the influence of parameters such as Rayleigh number (10⁶ ≤Ra ≤10⁸), Darcy number (10⁻⁵ ≤ Da ≤ 10⁻²), porosity (0.4 ≤ ɛ ≤ 0.9) and Prandtl number (0.01 ≤ Pr ≤ 10) on the flow behavior and heat transfer. By combining the method of matched asymptotic expansions with computational fluid dynamics (CFD), so called asymptotic computational fluid dynamics (ACFD) technique has been employed to generate correlation for average Nusselt number. The technique is found to be an attractive option for generating correlation and also in the analysis of natural convection in porous medium over a fairly wide range of parameters with fewer simulations for numerical solutions.     

Finite Element Model Tuning with Spatially-Dense 3D Modes

An investigation into varying the spatial density of three-dimensional (3D) modes in order to accurately tune a finite element (FE) model to larger numbers of modes is conducted. This FE model tuning approach is evaluated on a 60,000+ FE degree of freedom (DOF) model of the United States Air Force Academy’s fifth small satellite, FalconSAT-5 (FS-5), structural engineering model (SEM). The evaluation of this FE model tuning approach starts with the collection of natural frequencies and 3D modes extracted from scanning laser Doppler velocimeter (LDV) frequency response function measurements at 2,165 closely-spaced points on the surface of the test article. The measured modes and their associated natural frequencies serve as target values in a gradient-based tuning approach. The FE model is tuned to have differences with the measured natural frequencies less than 2% and in many cases cross-orthogonality values greater than 0.90. Using both a QR-decomposition and the Triaxial Effective Independence (EfI3⁺) sensor selection strategies, the effect of adding more tuning points on tuning accuracy is studied. This study demonstrates that increasing the number of experimental data points using the EfI3⁺ sensor selection strategy for FE model optimization results in increased tuning accuracy.     

Some bounded results of θ(t)-type singular integral operators

Let B be a Banach space in UMD with an unconditional basis. The boundedness of the θ(t)-type singular integral operators in L _B ^p (R ⁿ), (1≤p<+∞) and H _B ¹ (R ⁿ) spaces are discussed. Foundation item: the Education Commission of Shandong Province (J98P51) Biography: Zhao Kai (1960-)     

Effect of natural ventilation on the boundary layer separation and near-wake vortex shedding characteristics of a sphere

Experiments were conducted in water and wind tunnels on spheres in the Reynolds number range 6 × 10³ to 6.5 × 10⁵ to study the effect of natural ventilation on the boundary layer separation and near-wake vortex shedding characteristics. In the subcritical range of Re (<2 × 10⁵), ventilation caused a marginal downstream shift in the location of laminar boundary layer separation; there was only a small change in the vortex shedding frequency. In the supercritical range (Re > 4 × 10⁵), ventilation caused a downstream shift in the mean locations of boundary layer separation and reattachment; these lines showed significant axisymmetry in the presence of venting. No distinct vortex shedding frequency was found. Instead, a dramatic reduction occurred in the wake unsteadiness at all frequencies. The reduction of wake unsteadiness is consistent with the reduction in total drag already reported. Based on the present results and those reported earlier, the effects of natural ventilation on the flow past a sphere can be categorized in two broad regimes, viz., weak and strong interaction regimes. In the weak interaction regime (subcritical Re), the broad features of the basic sphere are largely unaltered despite the large addition of mass in the near wake. Strong interaction is promoted by the closer proximity of the inner and outer shear layers at supercritical Re. This results in a modified and steady near-wake flow, characterized by reduced unsteadiness and small drag. Received: 8 September 1998 / Accepted: 1 January 2000     

On near-wall hot-wire measurements

A specially constructed hot-wire probe was used to obtain very near-wall velocity measurements in both a fully developed turbulent channel flow and flat plate boundary layer flow. The near-wall hot-wire probe, having been calibrated in a specially constructed laminar flow calibration rig, was used to measure the mean streamwise velocity profile, distributions of streamwise and spanwise intensities of turbulence and turbulence kinetic energy k in the viscous sublayer and beyond; these distributions compare very favorably with available DNS results obtained for channel flow. While low Reynolds number effects were clearly evident for the channel flow, these effects are much less distinct for the boundary layer flow. By assuming the dissipating range of eddy sizes to be statistically isotropic and the validity of Taylor's hypothesis, the dissipation rate ɛ _iso in the very near-wall viscous sublayer region and beyond was determined for both the channel and boundary layer flows. It was found that if the convective velocity U _c in Taylor's hypothesis was assumed to be equal to the mean velocity Uˉ at the point of measurement, the value of (ɛ⁺ _iso)₁ thus obtained agrees well with that of (ɛ ⁺)_DNS for y ⁺ ≥ 80 for channel flow; this suggests the validity of assuming U _c=Uˉ and local isotropy for large values of y ⁺. However, if U _c was assumed to be 10.6u _τ , the value of (ɛ⁺ _iso)₂ thus obtained was found to compare reasonably well with the distribution of (ɛ⁺ _iso)_DNS for y ⁺≤ 15. Received: 31 May 1999/Accepted: 20 December 1999