Vortical structures in a turbulent far-wake: effect of Reynolds number

The effect of the Reynolds number on vortical structures in a turbulent far-wake has been investigated for Re_d (based on the free stream velocity and the cylinder diameter) =2800 and 9750. Velocity data were obtained using two orthogonal arrays of 16 X-wires, eight in the (x,y)-plane and eight in the (x,z)-plane. Structures were detected in both planes using a technique based on vorticity concentration and circulation. Conditional streamlines and contours of vorticity based on spanwise structures, i.e. detections in the (x,y)-plane, reveal that the streamwise size of spanwise structures increases as Re_d increases. The interrelationship is investigated between detections simultaneously identified in the two planes. Transverse structures, i.e. detections in the (x,z)-plane, correspond, with a relatively high probability, to spanwise structures, in conformity with a distortion in the (y,z)-plane of spanwise structures. Those that correspond, with relatively high probability, to the saddle between consecutive spanwise structures are interpreted in terms of ribs, whose signatures are detectable in instantaneous data. The probability is also high for transverse structures to occur between the focus of a spanwise structure and its associated saddle when Re_d=9750, but not when Re_d=2800. This is consistent with an increased vortex pairing frequency at the higher Re_d, as observed in instantaneous sectional streamlines.     

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.     

Uniformly valid approximations for non-linear oscillations with small damping

A uniformly valid zeroth-order approximation is obtained for the general equation y + εH(y)y + M(y)y = 0, where ε is a small parameter. The notion of multiple scaling is utilized to set up a systematic approximation scheme. Examples are given for simple polynomials for H(y) and M(y), which lead to results involving elliptic integrals. Further restrictions allow progress to be made in terms of gamma functions.     

The effect of axial stretch on cavitation in an elastic cylinder

The phenomenon of cavitation in an elastic cylinder subjected to a combined axial stretch λ_z and radial traction P is examined. One possible response of the cylinder to this loading, for all λ_z; and P, is a pure homogeneous deformation. However, for some materials, the homogeneous deformation becomes unstable at certain critical values of (λ_z, P), and a second deformation bifurcates from the homogeneous one; this latter deformation involves a cylindrical cavity co-axial with the elastic cylinder. We determine the values of (λ_z, P) at which this happens. The results are displayed by constructing a curve in the (λ_z, P)-plane which divides this plane into regions where the homogeneous deformation is stable and unstable.     

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.     

Convection regimes and heat transfer characteristics along a continuously moving heated vertical plate

The steady laminar flow and thermal characteristics of a continuously moving vertical sheet of extruded material are studied close to and far downstream from the extrusion slot. The velocity and temperature variations, obtained by a finite volume method, are used to map out the entire forced, mixed and natural convection regimes. The effects of the Prandtl number (Pr) and the buoyancy force parameter (B) on the friction and heat transfer coefficients are investigated. Comparisons with experimental measurements and solutions by others in the pure forced and pure natural convection regions are made. In the mixed convection region, the results are compared with available finite-difference solutions of the boundary layer equations showing excellent agreement. The region close to the extrusion slot is characterized as a non-similar forced-convection dominated region in which Nu_xRe_x^−1/2 drops sharply with increasing Richardson number (Ri_x). This is followed by a self-similar forced-convection dominated region in which Nu_xRe_x^−1/2 levels off with increasing Ri_x until the buoyancy effect sets in. The existence and extent of the latter region depend upon the value of B. A non-similar mixed convection region where increasing buoyancy effect enhances the heat transfer rate follows. Finally, this region is followed downstream by a self-similar natural-convection dominated region in which Nu_xRe_x^−1/2 approaches the pure natural convection asymptote at large Ri_x. Critical values of Ri_x to distinguish the various convection regimes are determined for different Pr and B.     

On rubber disks under rotation or axisymmetric stretching

The governing equations for a class of axisymmetric problems under large elastic deformation, concerned with a circular rubber disk with body force as well as non-uniform initial thickness, are formulated in terms of two coupled first-order ordinary differential equations with explicit derivatives. The following two problems subjected to different boundary conditions are solved: (a) Rotating disks with uniform initial thickness (b) Circular disks with non-uniform initial thickness under axisymmetric stretching at the outer boundary. In problem (b), the rubber disk whose initial thickness contour is h₀ = crⁿ (where c and n are any constants), or whose final thickness is a constant, is considered. Highly elastic materials with a Mooney strain-energy function are used for numerical calculations.     

Large-scale structure of a leading-edge separation bubble with local forcing

The Karhunen-Loeve (K-L) expansion is used to extract coherent structures from a leading-edge separation bubble with local forcing. A leading-edge separation bubble is simulated using the discrete vortex method, where a time-dependent source forcing is perturbed near the separation point. Based on the wealth of numerical data, the K-L procedure is applied in a range of the forcing amplitude (A₀ = 0, 0.5, 1.0 and 1.5) and forcing frequency (0 f_FH/U_∞ 0.3). Application of K-L procedure reveals that the eigenstructures are changed noticeably by local forcings. In an effort to investigate the mechanism of decreasing reattachment length (x_R), dynamic behaviors of the expansion coefficients and contributions of the eigenfunctions are scrutinized. As the forcing amplitude increases, large-scale vortex structures are formed near the separation point. Furthermore, the flow becomes more organized, which results in the reduction of x_R. Two distinctive regimes are classified: the regime of decreaing x_R and the regime of increasing x_R. The K-L global entropy indicates that x_R is closely linked to the organization of the flow structure.     

Can a single-wall carbon nanotube be modeled as a thin shell?

Single-wall carbon nanotubes (SWCNT) have been frequently modeled as thin shells, but the shell thickness and Young's modulus reported in literatures display large scattering. The order of error to approximate SWCNTs as thin shells is studied in this paper via an atomistic-based finite-deformation shell theory, which avoids the shell thickness and Young's modulus, but links the tension and bending rigidities directly to the interatomic potential. The ratio of atomic spacing (Δ≈0.14 nm) to the radius of SWCNT, Δ/R, which ranges from zero (for graphene) to 40% [for a small (5,5) armchair SWCNT (R=0.35 nm)], is used to estimate the order of error. For the order of error O[(Δ/R)³], SWCNTs cannot be represented by a conventional thin shell because their constitutive relation involves the coupling between tension and curvature and between bending and strain. For the order of error O[(Δ/R)²], the tension and bending (shear and torsion) rigidities of SWCNTs can be represented by an elastic orthotropic thin shell, but the thickness and elastic modulus cannot. Only for the order of error O(Δ/R), a universal constant shell thickness can be defined and SWCNTs can be modeled as an elastic isotropic thin shell.     

Conjugate natural convection from an array of discrete heat sources: part 2 — a numerical parametric study

Coupled conduction and natural convection transport within a discretely heated cavity have been investigated numerically. One vertical wall of the cavity is composed of discrete, isoflux heat sources mounted in a substrate of finite thermal conductivity. The opposite vertical wall and the horizontal walls are assumed to be isothermal and adiabatic, respectively. The governing steady-state partial differential equations for the fluid and solid region are solved simultaneously using a control volume formulation, coupled with an additive correction multigrid procedure that increases the convergence rate of the solution. The fluid Prandtl number and heater/fluid thermal conductivity ratio are fixed at 25 and 2350, respectively, corresponding to a dielectric fluid (FC-77) and heaters manufactured from silicon. With increasing modified Rayleigh number (10⁴ < Ra_Lz^* < 109), the cavity flow becomes more boundary layer-like along the vertical walls, and multiple fluid cells develop in the central region. Thermal spreading in the substrate increases with decreasing modified Rayleigh number and with increasing values of the substrate/fluid thermal conductivity ratio (10⁻¹ <- R_s ≤ 10³). For large R_s, the discrete heat sources lose their thermal identity, and the streamlines and isotherms resemble those associated with a differentially heated cavity. Thermal spreading in the substrate also has a significant effect on circulation in the cavity and on maximum surface temperatures.     

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.     

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).     

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.     

Pressure spectra in homogeneous low Reynolds number turbulent shear flow subjected to weak rotation

In this paper, pressure spectra have been derived from the authors’ model (Eur. J. Mech., B/Fluids 12 (1) (1993) 31–42) developed by means of rapid distortion theory (RDT) of homogeneous low Reynolds number turbulent shear flow subjected to weak rotation. The combined effects of uniform shear dU₁/dx₂ and weak rotation Ω₃ on the evolution of pressure spectra have been examined in terms of the rotation number 2Ω₃/(dU₁/dx₂). It is found that the system rotation exhibits the opposite effect on the pressure field as compared with the influence of rotation on the velocity fluctuations.     

二硫化钼/硫化锌纳米粉体的制备及其改性聚酰亚胺复合材料的摩擦学性能研究

采用一步水热法设计制备了二硫化钼/硫化锌(MoS₂/ZnS)纳米杂化体,并利用热压成型技术得到聚酰亚胺/二硫化钼/硫化锌(PI/MoS₂/ZnS)复合材料. 采用扫描电子显微镜、透射电子显微镜、X射线衍射仪以及光电子能谱仪对所制备材料的形貌和化学组成进行表征,结果表明MoS₂纳米薄片均匀致密地包覆在ZnS纳米颗粒表面. 热重分析和差示扫描量热曲线结果表明,MoS₂/ZnS纳米杂化体的引入显著地提升了PI基体的热稳定性能. 摩擦磨损测试结果表明,三种填料(MoS₂,ZnS和MoS₂/ZnS)均能有效改善PI基体的摩擦学性能,其中MoS₂/ZnS纳米杂化体的增强效应最为显著,这主要归因于MoS₂纳米片和ZnS纳米粒子之间的协同增强效应. 当MoS₂/ZnS纳米杂化体的质量分数为1.5%时, PI/MoS₂/ZnS复合材料的摩擦学性能达到最优,相较于纯的PI,复合材料的摩擦系数和磨损率分别下降了15. 9%和34. 3%.      

The X-structure/mechanism approach to beneficial nonlinear design in engineering

Nonlinearity can take an important and critical role in engineering systems, and thus cannot be simply ignored in structural design, dynamic response analysis, and parameter selection. A key issue is how to analyze and design potential nonlinearities introduced to or inherent in a system under study. This is a must-do task in many practical applications involving vibration control, energy harvesting, sensor systems, robotic technology, etc. This paper presents an up-to-date review on a cutting-edge method for nonlinearity manipulation and employment developed in recent several years, named as the X-structure/mechanism approach. The method is inspired from animal leg/limb skeletons, and can provide passive low-cost high-efficiency adjustable and beneficial nonlinear stiffness (high static & ultra-low dynamic), nonlinear damping (dependent on resonant frequency and/or relative vibration displacement), and nonlinear inertia (low static & high dynamic) individually or simultaneously. The X-structure/mechanism is a generic and basic structure/mechanism, representing a class of structures/mechanisms which can achieve beneficial geometric nonlinearity during structural deflection or mechanism motion, can be flexibly realized through commonly-used mechanical components, and have many different forms (with a basic unit taking a shape like X/K/Z/S/V, quadrilateral, diamond, polygon, etc.). Importantly, all variant structures/mechanisms may share similar geometric nonlinearities and thus exhibit similar nonlinear stiffness/damping properties in vibration. Moreover, they are generally flexible in design and easy to implement. This paper systematically reviews the research background, motivation, essential bio-inspired ideas, advantages of this novel method, the beneficial nonlinear properties in stiffness, damping, and inertia, and the potential applications, and ends with some remarks and conclusions.

     

血液-血管耦合特性与脉搏波传播特性的关系

脉搏波既不可简单地理解为可压缩血液流体中的压力纵波，也不可简单地理解为沿固体血管传播的涨缩位移横波，而是超乎普通想象的流-固耦合和纵波-横波耦合的复杂波。从分析耦合本构关系的新途径出发，本文中提出了一个流-固耦合/纵波-横波耦合的串联模型，可为解读“位数形势”中医脉诊提供更丰富的信息。结果表明，脉搏波耦合系统的等效体积压缩模量K_s以及相应的耦合系统脉搏波传播速度c_s主要依赖于两个无量纲参数:血液-血管模量比K_b(p)/E(p)和薄壁血管径厚比D(p)/h₀，它们因人而异、因人的不同脉搏位置而异。文中定量分析了它们对c_s的影响，显示人体的K_b/E值在103数量级，从而c_s值在100～101 m/s数量级，以适应人体生理生化反应。由临床有创测量，证实脉搏体积横波与脉搏压力纵波是相耦合地以相同速度传播；还显示脉搏波是在其波阵面上具有氧合生化反应的“生物波”。此外，还讨论了“脉压放大”现象与非线性本构关系和与血管分叉处加载增强反射之间的关系，并讨论了Lewis关于重搏波形成的假设。     

Effects of polymer additives on the large scale structure of a two-dimensional submerged jet

The large structures on the boundary of a two-dimensional submerged jet were examined. In the first experiment, a jet of a 100 ppm aqueous polyacrylamid solution was fed into pure water (P/W). In the second experiment, the jet as well as the ambient fluid were pure water (W/W), whereas in the third experiment, polymer solution (W/P) was in the ambient only.The large structures were found to be of the same size for all three experiments. The spread-angle was lowest for P/W and highest for W/P. These effects are attributed to a strong flow which stretches the polymer molecules in the mixing layer.     

Coherent fine scale eddies in turbulence transition of spatially-developing mixing layer

To investigate the relationship between characteristics of the coherent fine scale eddy and a laminar–turbulent transition, a direct numerical simulation (DNS) of a spatially-developing turbulent mixing layer with Re_ω,0 = 700 was conducted. On the onset of the transition, strong coherent fine scale eddies appears in the mixing layer. The most expected value of maximum azimuthal velocity of the eddy is 2.0 times Kolmogorov velocity (u_k), and decreases to 1.2u_k, which is an asymptotic value in the fully-developed state, through the transition. The energy dissipation rate around the eddy is twice as high compared with that in the fully-developed state. However, the most expected diameter and eigenvalues ratio of strain rate acting on the coherent fine scale eddy are maintained to be 8 times Kolmogorov length (η) and :β:γ = −5:1:4 in the transition process. In addition to Kelvin–Helmholtz rollers, rib structures do not disappear in the transition process and are composed of lots of coherent fine scale eddies in the fully-developed state instead of a single eddy observed in early stage of the transition or in laminar flow.