Studies on the mixing of two non-axial jets in a confined passage — turbulence characteristics

The mean characteristics of the flow due to mixing of two non-axial plane jets in a confined passage, were reported previously by Manjunath et al. (Exp. Fluids 11 (1991), 17–24), as part I. In this paper, the turbulence characteristics of the flow are discussed. The various components of the Reynolds stress tensor are reported for the four inlet angles considered, viz., 15°, 30°, 45° and 60°. As expected from the mean flow field, the inlet angle influences the distribution of the various turbulence components through the relative size of the recirculation zones and the turbulence field is different for the smaller angles of 15° and 30° in comparison to those for 45° and 60°.     

Bagley correction: the effect of contraction angle and its prediction

The excess pressure losses due to end effects (mainly entrance) in the capillary flow of a branched polypropylene melt were studied both experimentally and theoretically. These losses were first determined experimentally as a function of the contraction angle ranging from 10° to 150°. It was found that the excess pressure loss function decreases for the same apparent shear rate with increasing contraction angle from 10° to about 45°, and consequently slightly increases from 45° up to contraction angles of 150°. Numerical simulations using a multimode K-BKZ viscoelastic and a purely viscous (Cross) model were used to predict the end pressures. It was found that the numerical predictions do agree well with the experimental results for small contraction angles up to 30°. However, the numerical simulations under-predict the end pressure for larger contraction angles. The effects of viscoelasticity, shear, and elongation on the numerical predictions are also assessed in detail. Shear is the dominant factor controlling the overall pressure drop in flows through small contraction angles. Elongation becomes important at higher contraction angles (greater than 45°). It is demonstrated in abrupt contractions (angle of 180°) that both the entrance pressure loss and the vortex size are strongly dependent on the extensional viscosity for this branched polymer. It is suggested that such an experiment (visualisation of entrance flow) can be useful in evaluating the validity of constitutive equations and it can also be used to fitting parameters of rheological models that control the elongational viscosity.     

Assessment of 2D/3D numerical modeling for deep dynamic stall experiments

The results of computational fluid dynamics (CFD) simulations in two and three spatial dimensions are compared to pressure measurements and particle image velocimetry (PIV) flow surveys to assess the suitability of numerical models for the simulation of deep dynamic stall experiments carried out on a pitching NACA 23012 airfoil. A sinusoidal pitching motion with a 10° amplitude and a reduced frequency of 0.1 is imposed around two different mean angles of attack of 10° and 15°. The comparison of the airloads curves and of the pressure distribution over the airfoil surface shows that a three-dimensional numerical model can better reproduce the flow structures and the airfoil performance for the deep dynamic stall regime. Also, the vortical structures observed by PIV in the flow field are better captured by the three-dimensional model. This feature highlighted the relevance of three-dimensional effects on the flow field in deep dynamic stall.     

Natural convection in inclined two dimensional rectangular cavities

Steady two-dimensional natural convection in fluid filled cavities is numerically investigated. The channel is heated from below and cooled from the top with insulated side walls and the inclination angle is varied. The field equations for a Newtonian Boussinesq fluid are solved numerically for three cavity height based Rayleigh numbers, Ra = 10⁴, 10⁵ and 10⁶, and several aspect ratios. The calculations are in excellent agreement with previously published benchmark results. The effect of the inclination of the cavity to the horizontal with the angle varying from 0° to 180° and the effect of the startup conditions on the flow pattern, temperature distribution and the heat transfer rates have been investigated. Flow admits different configurations at different angles as the angle of inclination is increased depending on the initial conditions. Regardless of the initial conditions Nusselt number Nu exhibits discontinuities triggered by gradual transition from multiple cell to a single cell configuration. The critical angle of inclination at which the discontinuity occurs is strongly influenced by the assumed startup field. The hysteresis effect previously reported is not always present when the calculations are reversed from 90° to 0°. A comprehensive study of the flow structure, the Nu variation with varying angle of inclination, the effect of the initial conditions and the hysteresis effect are presented.     

Investigation of size-dependent quasistatic response of electrically actuated nonlinear viscoelastic microcantilevers and microbridges

This study investigates the size-dependent quasistatic response of a nonlinear viscoelastic microelectromechanical system (MEMS) under an electric actuation. To have this problem in view, the deformable electrode of the MEMS is modelled using cantilever and doubly-clamped viscoelastic microbeams. The modified couple stress theory in conjunction with Bernoulli–Euler beam theory are used for mathematical modeling of the size-dependent instability of microsystems in the framework of linear viscoelastic theory. Simultaneous effect of electrostatic actuation including fringing field, residual stress, mid-plane stretching and Casimir and van der Waals intermolecular forces are considered in the theoretical model. A single element of the standard linear solid element is used to simulate the viscoelastic behavior. Based on the extended Hamilton’s variational principle, the nonlinear governing integro-differential equation and boundary conditions are derived. Thereafter, a new generalized differential-integral quadrature solution for the nonlinear quasistatic response of electrically actuated viscoelastic micro/nanobeams under two different boundary conditions; doubly-clamped microbridge and clamped-free microcantilever. The developed model is verified and a good agreement is obtained. Finally, a comprehensive study is conducted to investigate the effects of various parameters such as material relaxation time, durable modulus, material length scale parameter, Casimir force, van der Waals force, initial gap and beam length on the pull-in response of viscoelastic microbridges and microcantilevers in the framework of viscoelasticity.     

横摆角下汽车尾部湍流特征量的PIV试验分析

汽车尾部湍流场是汽车压差阻力的主要来源,在HD-2汽车模型风洞中,首先使用测力天平和测压系统,对横摆角工况下汽车模型的气动六分力和纵对称截面48个测点的表面压力进行了测量,然后利用PIV测量技术对模型在横摆角分别为0°、15°的尾部湍流场进行了测量,获得该模型尾流场的速度场、涡量场和雷诺应力流场信息,通过计算得出尾流场区域空间相关系数和湍流积分尺度。结果表明:在横摆角工况下,汽车模型尾部涡流的结构呈现向上发展的趋势;尾流场拖拽涡的范围和强度的增大导致了模型气动力出现较大的增加;湍流积分尺度的变化表明,尾部涡流区的分离噪声与涡流分离位置有关,在汽车尾部造型设计中,要尽量推迟尾部涡流的分离。     

Two-dimensional planar flow of a viscoelastic plastic medium

The present study is concerned with finite element simulation of the planar entry flow of a viscoelastic plastic medium exhibiting yield stress. The numerical scheme is based on the Galerkin formulation. Flow experiments are carried out on a carbon black filled rubber compound. Steady-state pressure drops are measured on two sets of contraction or expansion dies having different lengths and a constant contraction or expansion ratio of 4:1 with entrance angles of 90, 45 and 15 degrees. The predicted and measured pressure drops are compared. The predicted results indicate that expansion flow has always a higher pressure drop than contraction flow. This prediction is in agreement with experimental data only at low flow rates, but not at high flow rates. The latter disagreement is possibly an indication that the assumption of fully-developed flow in the upstream and downstream regions is not realistic at high flow rates, even for the large length-to-thickness ratio channels employed. The evolution of the velocity, shear stress, and normal stress fields in the contraction or expansion flow and the location of pseudo-yield surfaces are also calculated.     

Experimental study on turbulent mixing of spray droplets in crossflow

Centrifugal spray injected at various angles in gas crossflow has been studied experimentally using PIV visualization system and image-processing techniques. Experiments were carried out inside a rectangular duct (95 mm × 95 mm in cross-section) at ambient temperature and pressure, with different gas Reynolds numbers (vary from 12,900 to 45,000) and three injection angles (60°, 90° and 120°). The spray angle of the centrifugal nozzle is 80°, with D₃₂ of 80 μm. The instantaneous images of droplets distribution and the values of the line-averaged D₃₂ at different positions on the cross-sections along the flow field for each condition were obtained, and their flow field configurations were achieved. Quantitative assessments of mixing degree between two phases for different injection angles were determined using a spatial unmixedness parameter. It is found that the addition of droplets into the gas crossflow enhanced the turbulence intensity of the gas crossflow and caused different-scale vortices. The flow field structure, to a great extent, is dependent on the injection angle. The entrainment and centrifugal force of large vortex lead to uneven droplet distribution and moreover influence the mixing of droplets and gas crossflow. A better mixing result can be obtained with the injection angle of nozzles of 60°.     

Yielding on inclined planes at the tip of a crack loaded in uniform tension

Plane-strain yielding from a crack in an infinite elastic body is represented here by a distribution of edge dislocations on two planes inclined at angles ±ga to the crack plane, and the equilibrium condition is solved numerically. Approximate analytical expressions are obtained for the plastic-zone length, the crack opening displacement, and the J-integral, as functions of the applied stress and α. A comparison with a co-planar model of the plastic zone gives very similar results for α ≈ 65°. It is shown that fracture criteria based either on a critical crack opening displacement (COD) or on a critical value of J are always different, and the use of the former may lead to critical defect-sizes which are twice as large as those given by the latter. Furthermore, COD appears not to be a well-defined material property. The critical J criterion gives a fracture stress which is α-dependent : this may be responsible for deviations towards results of linear elastic fracture mechanics when post-yield fracture mechanics is used to analyse extensive yielding. The changes in the stress field of the crack due to the existence of the plastic zone are also discussed.     

Two-dimensional viscoelastic flows: experimentation and modeling

Extensive experimental data on the birefringence in converging and diverging flows of a polymeric melt have been obtained. The birefringence and pressure drop measurements were carried out in working cells of planar geometry having different contraction angles and contraction ratios. For investigation of diverging or abrupt expansion flow, the direction of flow in the cells was reversed. The theoretical predictions are based upon the Leonov constitutive equation and a finite element scheme with streamwise integration.In contrast to Newtonian and second-order fluids, viscoelastic fluids at high shear rates show significant differences in pressure drop and birefringence (i.e. stresses) in converging and diverging flows. For a constant flow rate, the pressure drop is higher and the birefringence smaller in diverging flows than in converging flows. This difference increases with increasing flow rate. Further, for the same contraction ratio but different contraction angles, the birefringence maximum increases considerably with contraction angle. In addition, an increase in contraction ratio has the same effect.The viscoelastic constitutive equation of Leonov has been shown to describe all the above viscoelastic effects observed in the experiments. In general, a reasonable agreement between theory and experiment has been obtained, which shows the usefulness of the Leonov model in describing actual flows.     

Cycle resolved multi-planar flow measurements in a four-valve combustion engine

The non-reacting flow in a one-cylinder four-valve combustion engine is measured via cycle resolved two-component/two-dimensional (2C/2D) particle-image velocimetry (PIV). The three-dimensional structure of the velocity field is analyzed based on the flow field measured in eight planar planes within the cylinder for several crank angles during the intake and compression phase. Using the mean and statistical values of the single planes quasi three-dimensional flow fields are reconstructed for crank angles of 80°, 160°, and 240° atdc. This enables the detailed analysis of the spatial distribution of the large and small scale flow structures, e.g., by visualizing large vortical structures and the distribution of the turbulent kinetic energy. It was found that two ring vortices evolving beneath the inlet valves are the dominant large scale structures that seem to be of major concern for the mixing process in the cylinder of a four-valve combustion engine operated at 1500 rpm. Furthermore, the temporal evolution of the flow field within the symmetry plane of the cylinder, measured for crank angles between 40° and 320° atdc in steps of 20°, is discussed. The results give new insight into the complex three-dimensional flow in the combustion chamber of a one-cylinder four-valve combustion engine. That is, the tumble vortex only seems to be of secondary importance for the flow concerning the mixing process at 1500 rpm. This is an essential result for future work considering the fluid mechanics of fuel-air-interaction processes and mixing principles in combustion engines.     

Experimental study of compressible boundary layer on a cone at angles of attack

Experimental study was conducted for boundarylayers on a sharp 5° half-angle cone of 400mm length at angles of attack. The model was tested in the T-326 hypersonic wind tunnel （ITAM） at freestream Mach number M = 5.95. Mean and fluctuation wall characteristics of the boundary layer are measured at 0°, 2°, 3° and 4° angles of attack for different stagnation pressures. Pulsation measurements are carried out by means of ALTP sensor. Pressure and temperature distributions along the model are obtained, and transition beginning and end locations have been found. Boundary layer stabilization with the increase of angle of attack and the decrease of stagnation pressure is observed. High frequency pulsations inherent to hypersonic boundary layer （second mode） have been detected.     

Experimental and numerical study of the response of flexible laminates to impact loading

Fabrics and flexible laminates comprising highly oriented polymers possess high impact resistance and are often used in flexible armour applications. This study presents an idealized computational model for flexible [0°/90°] polyethylene fiber-reinforced laminates and show how the model can simulate actual impact tests conducted on the laminates. As these materials are viscoelastic, accurate modeling of their impact and perforation responses requires the formulation of constitutive equations representing such behavior. The material is idealized as networks of one-dimensional pin-jointed fiber elements defined by viscoelastic constitutive relations. Three-element viscoelastic Zener models are used as they are simple yet sufficient to account for the effects of intermolecular and intramolecular bonds, as well as viscous slippage between molecular chains, on the mechanical properties of oriented polymeric fibers. The effects of delamination in laminates are also taken into consideration by modeling flexible laminates as two fiber network layers bonded to each other at corresponding element cross-over points in the adjacent layers. Inter-ply bonding is represented by infinitesimal rigid links which break when the inter-ply bond strength is exceeded resulting in delamination between the plies. Predictions of residual velocity, development of deformation and delamination correlate well with experimental results.     

Time-resolved and volumetric PIV measurements of a transitional separation bubble on an SD7003 airfoil

To comprehensively understand the effects of Kelvin–Helmholtz instabilities on a transitional separation bubble on the suction side of an airfoil regarding as to flapping of the bubble and its impact on the airfoil performance, the temporal and spatial structure of the vortices occurring at the downstream end of the separation bubble is investigated. Since the bubble variation leads to a change of the pressure distribution, the investigation of the instantaneous velocity field is essential to understand the details of the overall airfoil performance. This vortex formation in the reattachment region on the upper surface of an SD7003 airfoil is analyzed in detail at different angles of attack. At a Reynolds number Re _c < 100,000 the laminar boundary layer separates at angles of attack >4°. Due to transition processes, turbulent reattachment of the separated shear layer occurs enclosing a locally confined recirculation region. To identify the location of the separation bubble and to describe the dynamics of the reattachment, a time-resolved PIV measurement in a single light-sheet is performed. To elucidate the spatial structure of the flow patterns in the reattachment region in time and space, a stereo scanning PIV set-up is applied. The flow field is recorded in at least ten successive light-sheet planes with two high-speed cameras enclosing a viewing angle of 65° to detect all three velocity components within a light-sheet leading to a time-resolved volumetric measurement due to a high scanning speed. The measurements evidence the development of quasi-periodic vortex structures. The temporal dynamics of the vortex roll-up, initialized by the Kelvin–Helmholtz (KH) instability, is shown as well as the spatial development of the vortex roll-up process. Based on these measurements a model for the evolving vortex structure consisting of the formation of c-shape vortices and their transformation into screwdriver vortices is introduced.     

Effect of dimples on glancing shock wave turbulent boundary layer interactions

An experimental study has been conducted to examine the control effectiveness of dimples on the glancing shock wave turbulent boundary layer interaction produced by a series of hemi-cylindrically blunted fins at Mach numbers 0.8 and 1.4, and at angles of sweep 0°, 15°, 30° and 45°. Schlieren photography, oil flow, pressure sensitive paints, and pressure tappings were employed to examine the characteristics of the induced flow field. The passive control technique used a series of 2 mm diameter, 1 mm deep indents drilled across the hemi-cylindrical leading edge at angles 0°, 45° and 90°. The effects of dimples were highly dependent on their orientation relative to the leading edge apex, and the local boundary layer properties.      

Velocity overshoots in gradual contraction flows

This study reports the results of a systematic numerical investigation, using the upper-convected Maxwell (UCM) and Phan-Thien–Tanner (PTT) models, of viscoelastic fluid flow through three-dimensional gradual planar contractions of various contraction ratios with the aim of investigating experimental observations of extremely large near-wall velocity overshoots in similar geometries [R.J. Poole, M.P. Escudier, P.J. Oliveira, Laminar flow of a viscoelastic shear-thinning liquid through a plane sudden expansion preceded by a gradual contraction, Proc. Roy. Soc. Lond. Ser. A 461 (2005) 3827]. We are able to obtain good qualitative agreement with the experiments, even using the UCM model in creeping-flow conditions, showing that neither inertia, second normal-stress difference nor shear-thinning effects are required for the phenomenon to be observed. Guided by the numerical results we propose a simple explanation for the occurrence of the velocity overshoots and the conditions under which they arise.     

拉伸荷载下分支裂隙破坏机理研究

为研究拉伸荷载下分支裂隙对破坏模式的影响,保持主裂隙参数不变,改变分支裂隙倾角和长度,利用扩展有限元方法模拟了弯折裂隙的动态扩展,总结了分支裂隙参数变化对破坏模式的影响.利用ABAQUS中的轮廓积分计算了分支裂隙尖端应力强度因子,并根据最大周向应力准则计算起裂角.结果表明:拉伸荷载下分支裂隙出现三种破坏模式;分支裂隙倾...     

压缩-剪切复合应力下PMMA材料的静态力学特性研究

通过压缩具有一定倾斜角(0°,10°,15°,20°和25°)试件和双剪切模型试件,实现了单轴压缩、压缩-剪切复合应力以及纯剪切三种应力状态,得到PMMA(聚甲基丙烯酸甲酯)在相应应力状态下的应力-应变曲线,同时对不同应力状态下试件的破坏模式进行了分析。结果表明:在不同受力环境中材料的强度和破坏的机理不同;同单轴压缩状态下相比,材料在压缩-剪切复合应力状态下屈服极限、强度极限以及破坏应变均不同程度的增大,呈现明显的"剪切增强"现象。单轴压缩与压缩-剪切应力状态下试件的破坏模式均为在试件短对角面上出现明显的剪切屈服带,由应力分析得出试件剪应力在短对角面上达到最大,引起在此平面上分子链间滑动从而产生应变软化形成剪切屈服带;双剪切试件的破坏模式为与剪切面呈45°的斜面。     

Experimental investigation of jointed rock breaking under a disc cutter with different confining stresses

Extensive and detailed investigations have been made to better understand the rock-breaking mechanism of the tunnel boring machine (TBM) disc cutter, but the crack propagation and failure modes induced by the disc cutter when the confining stresses and joint characteristics vary have not been comprehensively investigated. To address this area of research, a triaxial testing machine (TRW-3000) is modified to investigate the effect of different confining stresses (0, 2.5, 5, 7.5 and 10 MPa) on the rock breaking of different joint angles (0°, 30°, 60°, and 90°) induced by the disc cutter. In this series of tests, the crack propagation and failure modes of the intact and jointed rock with different confining stresses are analysed. During the experiments, four different types of failure modes have been observed. The failure mode is affected by the joint orientation at low confining stress. The existence of joints has no obvious effect on the failure mode when the confining stress increases to a certain extent.