Large Eddy Simulations of convergent–divergent channel flows at moderate Reynolds numbers

The paper is focused on the study of fully turbulent channel flows, using Large Eddy Simulations (LES), in order to address the effects of adverse pressure gradient regions. Analyses of the effects of streak instabilities, which have been shown to be relevant in such regions, are extended to moderate Reynolds numbers. The work considers two different channel geometries in order to further separate influences from wall curvature, flow separation and adverse pressure gradients. Turbulent kinetic energy and Reynolds stress budgets are investigated at separation and re-attachment points. The numerical approach used in the present work is based on the incompressible Navier–Stokes equations, which are solved by a pseudo-spectral methodology for structured grids. Wall-resolved LES calculations are performed using the WALE subgrid scale model. The study shows that the streak instability mechanism persists at higher Reynolds numbers with and without wall curvature in the adverse pressure gradient regions. Moreover, the observed effects are also present regardless of the existence of flow separation regions. Finally, the study of turbulent kinetic energy budgets indicates that, independently of the flow condition, there are well-defined patterns for such turbulent properties at separation and re-attachment points.     

Shock wave bifurcation in convergent–divergent channels of rectangular cross section

This work addresses two- and three-dimensional turbulent flow in simple channels, modeling the air intakes of rectangular cross section. Flow regimes with a supersonic free stream and supersonic velocities at the throat or immediately downstream of the throat are considered. Bifurcations of the shock wave arising ahead of the cowl are studied numerically. Solutions of the Reynolds-averaged Navier–Stokes equations are obtained with a finite-volume solver of second-order accuracy on fine computational meshes. The solutions reveal jumps of the shock leg position with variations of the free-stream Mach number. The dependence of the shock position on the cowl slope and streamwise location of the throat is examined.     

Numerical modeling and experimental investigation of gas–liquid slug formation in a microchannel T-junction

Gas–liquid two-phase flow in a microfluidic T-junction with nearly square microchannels of 113 μm hydraulic diameter was investigated experimentally and numerically. Air and water superficial velocities were 0.018–0.791 m/s and 0.042–0.757 m/s, respectively. Three-dimensional modeling was performed with computational fluid dynamics (CFD) software FLUENT and the volume of fluid (VOF) model. Slug flow (snapping/breaking/jetting) and stratified flow were observed experimentally. Numerically predicted void fraction followed a linear relationship with the homogeneous void fraction, while experimental values depended on the superficial velocity ratio U_G/U_L. Higher experimental velocity slip caused by gas inlet pressure build-up and oscillation caused deviation from numerical predictions. Velocity slip was found to depend on the cross-sectional area coverage of the gas slug, the formation of a liquid film and the presence of liquid at the channel corners. Numerical modeling was found to require improvement to treat the contact angle and contact line slip, and could benefit from the use of a dynamic boundary condition to simulate the compressible gas phase inlet reservoir.     

A combined experimental–numerical investigation of ductile fracture in bending of a class of ferritic–martensitic steel

We present a combined experimental–numerical study on fracture initiation at the convex surface and its propagation during bending of a class of ferritic–martensitic steel. On the experimental side, so-called free bending experiments are conducted on DP1000 steel sheets until fracture, realizing optical and scanning electron microscopy analyses on the post mortem specimens for fracture characterization. A blended Mode I – Mode II fracture pattern, which is driven by cavitation at non-metallic inclusions as well as martensitic islands and resultant softening-based intense strain localization, is observed. Phenomena like crack zig-zagging and crack alternation at the bend apex along the bending axis are introduced and discussed. On the numerical side, based on this physical motivation, the process is simulated in 2D plane strain and 3D, using Gurson’s dilatant plasticity model with a recent shear modification, strain-based void nucleation, and coalescence effects. The effect of certain material parameters (initial porosity, damage at coalescence and failure, shear modification term, etc.), plane strain constraint and mesh size on the localization and the fracture behavior are investigated in detail.     

An experimental investigation on the characteristics of fluid–structure interactions of a wind turbine model sited in microburst-like winds

An experimental investigation is performed to assess the characteristics of the fluid–structure interactions and microburst-induced wind loads acting on a wind turbine model sited in microburst-liked winds. The experiment study was conducted with a scaled wind turbine model placed in microburst-like winds generated by using an impinging-jet-typed microburst simulator. In addition to quantifying complex flow features of microburst-like winds, the resultant wind loads acting on the turbine model were measured by using a high-sensitive force–moment sensor as the turbine model was mounted at different radial locations and with different orientation angles with respect to the oncoming microburst-like winds. The measurement results reveal clearly that, the microburst-induced wind loads acting on the turbine model were distinctly different from those in a conventional atmospheric boundary layer (ABL) wind. With the scales of the wind turbine model and the microburst-like wind used in the present study, the dynamic wind loadings acting on the turbine model were found to be significantly higher (i.e., up to 4 times higher for the mean loads, and up to 10 times higher for the fluctuation amplitudes) than those with the same turbine model sited in ABL winds. Both the mean values and fluctuation amplitudes of the microburst-induced wind loads were found to vary significantly with the changes of the mounted site of the turbine model, the operating status (i.e., with the turbine blades stationary or freely rotating), and the orientation angle of the turbine model with respect to the oncoming microburst-like wind. The dynamic wind load measurements were correlated to the flow characteristics of the microburst-like winds to elucidate underlying physics. The findings of the present study are helpful to gain further insight into the potential damage caused by the violent microbursts to wind turbines to ensure safer and more efficient operation of the wind turbines in thunderstorm-prone areas.     

Frictional pressure drop during vapour–liquid flow in minichannels: Modelling and experimental evaluation

Condensation in minichannels is widely used in air-cooled condensers for the automotive and air-conditioning industry, in heat pipes and other applications for system thermal control. The knowledge of pressure drops in such small channels is important in order to optimize heat transfer surfaces. This paper presents a model for calculation of the frictional pressure gradient during condensation or adiabatic liquid–gas flow inside minichannels with different surface roughness. In order to account for the effects of surface roughness, new experimental frictional pressure gradient data associated to single-phase flow and adiabatic two-phase flow of R134a inside a single horizontal mini tube with rough wall has been used in the modelling. It is a Friedel (1979) [Friedel, L., 1979. Improved friction pressure drop correlations for horizontal and vertical two-phase pipe flow. In: Proceedings of the European Two-Phase Flow Group Meeting, Ispra, Paper E2] based model and it takes into account mass velocity, vapor quality, fluid properties, reduced pressure, tube diameter, entrainment ratio and surface roughness. With respect to the flow pattern prediction capability, it has been built for shear dominated flow regimes inside pipes, thus, annular, annular-mist and mist flow are here predicted. However, the suggested procedure is extended to the intermittent flow in minichannels and it is also applied with success to horizontal macro tubes.     

An experimental and numerical investigation of the behaviour of AA5083 aluminium alloy in presence of the Portevin–Le Chatelier effect

An experimental investigation of the Portevin–Le Chatelier (PLC) effect in the aluminium alloy AA5083-H116 is undertaken in this study through five different tests involving round, prismatic and flat notched specimen geometries. Measurements based on strain gages and digital image correlation (DIC) are used to capture and characterize the spatio-temporal features of the PLC behaviour. Inhomogeneous deformation with various localization bands caused by the PLC effect is observed in all tests, and the band characteristics are measured. The McCormick elastic–viscoplastic constitutive relation, developed for metals exhibiting this type of dynamic strain aging, is then described in detail, before the various parameters required by the model are determined based on available material tests. The model is finally used in full-scale 3D numerical simulations of the physical tests using the explicit solver of the non-linear finite element code LS-DYNA. It will be shown that the numerical results are able to reproduce most of the experimentally observed phenomena with reasonable accuracy. However, if the model is used to study the micromechanical mechanisms controlling the macromechanical behaviour of materials exhibiting PLC effects (such as the band morphology), more advanced constitutive relations may be required.     

Experimental investigation of dynamics of premixed acetylene–air flames in a micro-combustor

This paper describes a detailed experimental study performed to investigate the flame propagation behaviour of premixed flames in micro-channels. A novel, modular, stackable micro-combustor was developed for this purpose. For a chosen planar channel geometry, the flow condition and the mixture equivalence ratio of premixed acetylene–air were varied to investigate various modes of operation. Three different modes of operation were observed; they were (i) stable periodic operation – consisting of ignition, flame propagation, flame extinction, and re-ignition, (ii) a-periodic operation, and (iii) anchored flame condition. The present work also aims to provide quantitative information on the dynamics of premixed acetylene–air flames propagating inside micro-channels. A novel measurement approach based on OH^* chemiluminescence measurements employing a single photomultiplier unit was developed for this purpose. The data recorded were post processed using an in-house developed MATLAB code to evaluate the mean flame propagation speed measured between three different spatial locations along the length of the micro-channel. The results from the flame propagation speed measurements performed during ‘periodic’ mode of operation indicated that the flame travelled at higher propagation speed in the mid-length region of the channel compared to that at the initial entry point, suggesting flame acceleration. This flame acceleration could be attributed to a situation where the flame experienced different local equivalence ratio conditions at different upstream locations. The results suggest that after completion of a cycle of operation consisting of ignition, flame propagation and flame extinction, the fresh mixture that filled the channel was diluted with the exhaust gas from the previous cycle. This pocket of diluted mixture convected downstream with time, thus enabling the spatial variation in local equivalence ratio along the micro-channel.     

A computational study of the wing–wing and wing–body interactions of a model insect

The aerodynamic interaction between the contralateral wings and between the body and wings of a model insect are studied, by using the method of numerically solving the Navier-Stokes equations over moving overset grids, under typical hovering and forward flight conditions. Both the interaction between the contralateral wings and the interaction between the body and wings are very weak, e.g. at hovering, changes in aerodynamic forces of a wing due to the present of the other wing are less than 3% and changes in aerodynamic forces of the wings due to presence of the body are less than 2%. The reason for this is as following. During each down- or up-stroke, a wing produces a vortex ring, which induces a relatively large jet-like flow inside the ring but very small flow outside the ring. The vortex rings of the left and right wings are on the two sides of the body. Thus one wing is outside vortex ring of the other wing and the body is outside the vortex rings of the left and right wings, resulting in the weak interactions.     

Development and experimental investigation of a Quadrotor’s robust generalized dynamic inversion control system

Nonlinear Dynamics - The development of a two-loops Quadrotor’s robust generalized dynamic inversion (RGDI)-based control system is presented. The outer (position) loop utilizes PD position...     

Numerical and experimental studies of stick–slip oscillations in drill-strings

Nonlinear Dynamics - The cyclic nature of the stick–slip phenomenon may cause catastrophic failures in drill-strings or at the very least could lead to the wear of expensive equipment....     

An experimental investigation of a superplastic constitutive equation in Al–Mg–Si alloy composites reinforced with Si3N4 whiskers

Superplastic properties at 818 K were investigated by tensile tests for Al–Mg–Si alloy composites reinforced with Si₃N₄ whiskers whose volume fractions were 0–30%. The 20 and 30 vol.% composites exhibited large elongation of 615 and 285% at a high strain rate of 2×10⁻¹ s⁻¹, respectively. High strain rate superplasticity of the composites is attributed to the very small grain size of less than 3 μm. The stress–strain rate relation for the composites was almost the same as that of the alloy, taking into consideration the influences of threshold stress and grain size, and the relation was independent of the volume fraction of whisker. This is probably because grain boundary sliding was not hampered by the whiskers due to the presence of liquid phase for the composites.     

The effects of channel diameter on flow pattern,void fraction and pressure drop of two-phase air–water flow in circular micro-channels

Two-phase air–water flow characteristics are experimentally investigated in horizontal circular micro-channels. Test sections are made of fused silica. The experiments are conducted based on three different inner diameters of 0.53, 0.22 and 0.15 mm with the corresponding lengths of 320, 120 and 104 mm, respectively. The test runs are done at superficial velocities of gas and liquid ranging between 0.37–42.36 and 0.005–3.04 m/s, respectively. The flow visualisation is facilitated by systems mainly including stereozoom microscope and high-speed camera. The flow regime maps developed from the observed flow patterns are presented. The void fractions are determined based on image analysis. New correlation for two-phase frictional multiplier is also proposed for practical applications.     

DEM investigation of strain behaviour and force chain evolution of gravel–sand mixtures subjected to cyclic loading

The strain characteristic and load transmission of mixed granular matter are different from those of homogeneous granular matter. Cyclic loading renders the mechanical behaviours of mixed granular matter more complex. To investigate the dynamic responses of gravel–sand mixtures, the discrete element method (DEM) was used to simulate the cyclic loading of gravel–sand mixtures with low fines contents. Macroscopically, the evolution of the axial strain and volumetric strain was investigated. Mesoscopically, the coordination number and contact force anisotropy were studied, and the evolution of strong and weak contacts was explored from two dimensions of loading time and local space. The simulation results show that increasing fines content can accelerate the development of the axial strain and volumetric strain but has little effect on the evolution of contact forces. Strong contacts tend to develop along the loading boundary, presenting the spatial difference. Weak contacts are firstly controlled by confining pressure and then controlled by axial stress, while strong contacts are mainly controlled by axial stress throughout the whole cyclic loading. Once compression failure occurs, the release of axial stress causes the reduction of strong contact proportion in all local regions. These findings are helpful to understand the dynamic responses of gravel–sand mixtures, especially in deformation behaviours and the Spatio-temporal evolution of contact forces.     

Visualization of deflagration-to-detonation transitions in a channel with repeated obstacles using a hydrogen–oxygen mixture

The deflagration-to-detonation transition in a 100 mm square cross-section channel was investigated for a highly reactive stoichiometric hydrogen oxygen mixture at 70 kPa. Obstacles of 5 mm width and 5, 10, and 15 mm heights were equally spaced 60 mm apart at the bottom of the channel. The phenomenon was investigated primarily by time-resolved schlieren visualization from two orthogonal directions using a high-speed video camera. The detonation transition occurred over a remarkably short distance within only three or four repeated obstacles. The global flame speed just before the detonation transition was well below the sound speed of the combustion products and did not reach the sound speed of the initial unreacted gas for tests with an obstacle height of 5 and 10 mm. These results indicate that a detonation transition does not always require global flame acceleration beyond the speed of sound for highly reactive combustible mixtures. A possible mechanism for this detonation initiation was the mixing of the unreacted and reacted gas in the vicinity of the flame front convoluted by the vortex present behind each obstacle, and the formation of a hot spot by the shock wave. The final onset of the detonation originated from the unreacted gas pocket, which was surrounded by the obstacle downstream face and the channel wall.     

Study of laminar–turbulent flow transition under pulsatile conditions in a constricted channel

In this paper, direct numerical simulation is performed to investigate a pulsatile flow in a constricted channel to gain physical insights into laminar–turbulent–laminar flow transitions. An in-house computer code is used to conduct numerical simulations based on available high-performance shared memory parallel computing facilities. The Womersley number tested is fixed to 10.5 and the Reynolds number varies from 500 to 2000. The influences of the degree of stenosis and pulsatile conditions on flow transitions and structures are investigated. In the region upstream of the stenosis, the flow pattern is primarily laminar. Immediately after the stenosis, the flow recirculates under an adverse streamwise pressure gradient, and the flow pattern transitions from laminar to turbulent. In the region far downstream of the stenosis, the flow becomes re-laminarised. The physical characteristics of the flow field have been thoroughly analysed in terms of the mean streamwise velocity, turbulence kinetic energy, viscous wall shear stresses, wall pressure and turbulence kinetic energy spectra.     

Hydrodynamics of gas-solid fluidization of a homogeneous ternary mixture in a conical bed： Prediction of bed expansion and bed fluctuation ratios

Hydrodynamic characteristics of fluidization in a conical or tapered bed differ from those in a columnar bed because the superficial velocity in the bed varies in the axial direction. Fixed and fluidized regions could coexist and sharp variations in pressure drop could occur, thereby giving rise to a noticeable pressure drop-flow rate hysteresis loop under incipient fluidization conditions. To explore these unique properties, several experiments were carried out using homogeneous, well-mixed, ternary mixtures with three dif- ferent particle sizes at varying composition in gas-solid conical fluidized beds with varying cone angles. The hydrodynamic characteristics determined include the minimum fluidization velocity, bed fluctuation, and bed expansion ratios. The dependence of these quantities on average particle diameter, mass fraction of the fines in the mixture, initial static bed height, and cone angle is discussed. Based on dimensional analysis and factorial design, correlations are developed using the system parameters, i.e. geometry of the bed （cone angle）, particle diameter, initial static bed height, density of the solid, and superficial velocity of the fluidizing medium. Experimental values of minimum fluidization velocity, bed fluctuation, and bed expansion ratios were found to agree well with the developed correlations.     

Relaxation oscillations and the mechanism in a periodically excited vector field with pitchfork–Hopf bifurcation

Nonlinear Dynamics - The traditional geometric singular perturbation theory and the slow–fast analysis method cannot be directly used to explore the mechanism of the relaxation oscillations...     

Combined PDA and LIF applied to size–temperature correlations measurements in a heated spray

This paper aims to demonstrate the possibility to achieve droplet temperature measurements per droplet size class by combining two-color laser-induced fluorescence (LIF) and phase Doppler analyzer (PDA). For that purpose, PDA and LIF signal acquisitions are synchronized on the same time base. LIF signal is processed on each of the defined size classes in order to derive the droplet temperature. Since PDA is roughly sensitive to D ² and LIF roughly to D ³, the detection range of the combination of the two techniques in term of droplet size is carefully analyzed. Finally, the technique is demonstrated on a spray of n-decane injected in a turbulent over-heated air flow. The influence of the droplet size and Stokes number on the heating process of the droplets is clearly highlighted.     

Experimental investigation of transients-induced fluid–structure interaction in a pipeline with multiple-axial supports

Fluid–Structure Interaction (FSI) in pipes can significantly affect pressure fluctuations during water hammer event. In transmission pipelines, anchors with axial stops have an important role in the waterhammer-induced FSI as they can suppress or allow the propagation of additional stress waves in the pipe wall. More specifically, a reduction in the number of axial stops and/or their stiffness causes significant oscillations in the observed pressure signal due to the enhancement of Poisson’s coupling. To confirm these physical arguments, this research conducts experimental investigations and then processes the collected pressure signals. The laboratory tests were run on an anchored pipeline with multiple axial supports which some of them removed at some sections to emerge Poisson’s coupling. The collected pressure signals are analyzed in the time and frequency domain in order to decipher fluctuations that stem from Poisson coupling and other anchors effects. The analysis of the laboratory data reveals that the pattern of the time signals of pressure is primarily affected by the stiffness and location of the supports. Likewise, the properties of structural boundaries characterize the frequency spectrum of the transient pressures, which is manifested by altering the amplitudes corresponding to dominant frequencies of the system. The study is of particular importance in practice of transient based defect detections and pipe system design.