Global stability analysis of a 90°-bend pipe flow

The present work investigates the stability properties of the flow in a 90°-bend pipe with curvature $\delta =R/R_c=1/3$, with R being the radius of the cross-section of the pipe and $R_c$ the radius of curvature at the pipe centreline. Direct numerical simulations (DNS) for values of the bulk Reynolds number ${\mathit{Re}}_b=U_{b}D/\nu$ between 2000 and 3000 are performed. The bulk Reynolds number is based on the bulk velocity $U_b$, the pipe diameter D, and the kinematic viscosity $\nu$. The flow is found to be steady for ${\mathit{Re}}_b⩽2500$, with two main pairs of symmetric, counter-rotating vortices in the section of the pipe downstream of the bend. The presence of two recirculation regions is detected inside the bend: one on the outer wall and the other on the inner side. For ${\mathit{Re}}_b⩾2550$, the flow exhibits a periodic behaviour, oscillating with a fundamental non-dimensional frequency $\mathit{St}=\mathit{fD}/U_b=0.23$. A global stability analysis is performed in order to determine the cause of the transition from the steady to the periodic regime. The spectrum of the linearised Navier-Stokes operator reveals a pair of complex conjugate eigenvalues with positive real part, hence the transition is ascribed to a Hopf bifurcation occurring at ${\mathit{Re}}_{b,\mathit{cr}}\approx 2531$, a value much lower than the critical Reynolds number for the flow in a torus with the same curvature. The velocity components of the unstable direct and adjoint eigenmodes are investigated, and they display a large spatial separation, most likely due to the non-normality of the linearised Navier-Stokes operator. Thus, the core of the instability, also known in the literature as the wavemaker, is sought performing an analysis of the structural sensitivity of the unstable eigenmode to spatially localised feedbacks. The region located 15° downstream of the bend inlet, on the outer wall, is the most receptive to this kind of perturbations, and thus corresponds to where the instability originates. Since this region coincides with the outer-wall separation bubble, it is concluded that the instability is linked to the strong shear by the backflow phenomena. The present results are relevant for technical applications where bent pipes are frequently used, and their stability properties have hitherto not been studied.     

Étude numérique du couplage de la convection naturelle avec le rayonnement de surfaces en cavité carrée remplie d'airNumerical study of natural convection-surface radiation coupling in air-filled square cavities

A numerical tool is developed for coupling natural convection in cavities with surface radiation and computations are performed for an air-filled square cavity whose four walls have the same emissivity. Compared to the adiabatic case without radiation, the top wall is cooled, the bottom wall is heated, air flow along the horizontal walls are reinforced and thermal stratification in cavity core is reduced. Detailed analysis shows that net radiative heat flux is linear with ΔT if $\mathrm{\Delta }T?T_0$, which is the case at low Rayleigh number, and that radiative Nusselt number is a linear function of the cavity height. Surface radiation induces an early transition to time-dependent flows: for $?=0.2$ and a cavity height of 0.335 m the critical Rayleigh number is equal to $9.3\times {10}^6$ and the corresponding Hopf bifurcation is supercritical. Furthermore, multiple periodic solutions are observed between $\mathit{Ra}=1.2\times {10}^7$ and $1.3\times {10}^7$. To cite this article: H. Wang et al., C. R. Mecanique 334 (2006).     

Effect of the kinematic hardening on the plastic anisotropy parameters for metallic sheets

The initial plastic anisotropy parameters are conventionally determined from the Lankford strain ratios defined by $r^{\psi }=\frac{{\epsilon }_{22}^{\mathrm{p}\psi }}{{\epsilon }_{33}^{\mathrm{p}\psi }}$ (ψ being the direction of the loading path). They are usually considered as constant parameters that are determined at a given value of the plastic strain far from the early stage of the plastic flow (i.e. equivalent plastic strain of ${\epsilon }_{\mathrm{eq}}^{\mathrm{p}}=0.2\text{\%}$) and typically at an equivalent plastic strain in between 20% and 50% of plastic strain failure (or material ductility). What prompts to question about the relevance of this determination, considering that this ratio does not remain constant, but changes with plastic strain. Accordingly, when the nonlinear evolution of the kinematic hardening is accounted for, the Lankford strain ratios are expected to evolve significantly during the plastic flow.In this work, a parametric study is performed to investigate the effect of the nonlinear kinematic hardening evolution of the Lankford strain ratios for different values of the kinematic hardening parameters. For the sake of clarity, this nonlinear kinematic hardening is formulated together with nonlinear isotropic hardening in the framework of anisotropic Hill-type (1948) yield criterion. Extension to other quadratic or non-quadratic yield criteria can be made without any difficulty. This parametric study is completed by studying the effect of these parameters on simulations of sheet metal forming by large plastic strains.     

Two-phase flow patterns and size distribution of droplets in a microfluidic T-junction: Experimental observations in the squeezing regime

Generating micrometer sized droplets has been studied in a microfluidic system with T-junction geometry 250 μm in internal diameter and with PTFE capillary tubing. Several experiments were conducted by varying the flow rate of the dispersed phase from $2.78\cdot {10}^{-11}{\text{ m}}^3/\text{s}$ to $5.28\cdot {10}^{-9}{\text{ m}}^3/\text{s}$ and that of the continuous phase from $2.78\cdot {10}^{-10}{\text{ m}}^3/\text{s}$ to $1.94\cdot {10}^{-9}{\text{ m}}^3/\text{s}$. The visualization of different flow regimes (drop, plug, and annular) was carried out for three configurations (not inverted in a horizontal position, inverted in a horizontal position, and inverted in a vertical position) for low capillary numbers. The model of Gauss was also chosen for a droplet size distribution in the dispersed phase, with the flow quality x varying from 0.016 to 0.44. The evolution of the drop size distribution as a function of the flow quality in the dispersed phase shows that the variation coefficient of the droplet's diameter is inversely proportional to the flow quality.     

Non-Darcian Bénard convection in eccentric annuli containing spherical particles

In this paper, a numerical investigation of natural convection in a porous medium confined by two horizontal eccentric cylinders is presented. The cylinders are impermeable to fluid motion and retained at uniform different temperatures. While, the annular porous layer is packed with glass spheres and fully-saturated with air, and the cylindrical packed bed is under the condition of local thermal non-equilibrium. The mathematical model describing the thermal and hydrodynamic phenomena consists of the two-phase energy model coupled by the Brinkman-Forchheimer-extended Darcy model under the Boussinesq approximation. The non-dimensional derived system of formulations is numerically discretised and solved using the spectral-element method. The investigation is conducted for a constant cylinder/particle diameter ratio ($D_i/d$) = 30, porosity ($\epsilon$) = 0.5, and solid/fluid thermal conductivity ratio ($k_r$) = 38.6. The effects of the vertical, horizontal and diagonal heat source eccentricity (−0.8 $⩽e⩽$0.8) and the annulus radius ratio (1.5 $⩽\mathit{RR}⩽$ 5.0) on the temperature and velocity distributions as well as the overall heat dissipation within both the fluid and solid phases, for a broad range of Rayleigh number (10⁴ $⩽$ Ra $⩽$ 8 $\times {10}^7$). The results show that uni-cellular, bi-cellular and tri-cellular flow regimes appear in the vertical eccentric annulus at the higher positive eccentricity e = 0.8 as Rayleigh number increases. However, in the diagonal eccentric annulus, the multi-cellular flow regimes are shown to be deformed and the isotherms are particularly distorted when Rayleigh number increases. In contrast, in the horizontal eccentric annulus, it is found that whatever the Rayleigh number is only an uni-cellular flow regime is seen. In addition, it is shown that the fluid flow is always unstable in the diagonal eccentric geometry at e = 0.8 for moderate and higher Rayleigh numbers. However, it loses its stability in the vertical eccentric geometry only at two particular cases, while it is always stable in the horizontal eccentric geometry, for all eccentricities and Rayleigh numbers.     

Acousto-mechanical behaviour of ex-vivo skin: Nonlinear and viscoelastic properties

The mechanical behaviour of skin is significant for some applications including dermatology, surgery, and impact biomechanics science. In this work, we have investigated the study of the acousto-mechanical viscoelastic properties of skin. For that, both tensile-relaxation and ultrasonic tests were conducted on porcine tissue samples in fibre directions. To understand the complex skin aging phenomena, we used strength tensile test correlated with the Nonlinear Time Reversal signal processing tool extension “TR-NEWS”. Uniaxial tensile tests were carried out at a strain rate of $5?{10}^{?3}\text{ mm}{\text{s}}^{?1}$ on skin using a load-relaxation-discharge load path with increasing amplitude and offset. This work is also under way to extend the frequency range of ultrasounds to 50 MHz. Digital Image Correlation was used for 2D strain measurement of the dermis. From this analysis, we conclude that fresh porcine skin should be modelled as a nonlinear viscoelastic material with strain-rate dependence. The obtained hysteresis loop shall be taken as significant skin damage.     

Freely vibrating circular cylinder in the vicinity of a stationary wall

We present a numerical study on vortex-induced vibration (VIV) of a freely vibrating two degree-of-freedom circular cylinder in close proximity to a stationary plane wall. Fully implicit combined field scheme based on Petrov–Galerkin formulation has been employed to analyze the nonlinear effects of wall proximity on the vibrational amplitudes and hydrodynamic forces. Two-dimensional simulations are performed as a function of decreasing gap to cylinder diameter ratio $e/D\in [0.5,10]$ for reduced velocities $U^{⁎}\in [2,10]$ at Re_D=100 and Re_L=2900, where Re_D and Re_L denote the Reynolds numbers based on the cylinder diameter and the upstream distance, respectively. We investigate the origin of enhanced streamwise oscillation of freely vibrating near-wall cylinder as compared to the isolated cylinder counterpart. For that purpose, detailed analysis of the amplitudes, frequency characteristics and the phase relations has been performed for the isolated and near-wall configurations. Initial and lower branches in the amplitude response are found from the gap ratios of 0.75 to 10, similar in nature to the isolated cylinder laminar VIV. A third response branch has been found between the initial and the lower branch at the gap ratio of $e/D\le 0.60$. For near-wall cases, phase relation between drag force and streamwise displacement varies from close to 0° to 180°. Between $e/D\in [5,7.5]$, the effect of wall proximity on the frequency response tends to disappear. The effect of mass-ratio is further investigated. Finally, we introduce new correlations for characterizing peak amplitudes and forces as a function of the gap ratio for a cylinder vibrating in the vicinity of a stationary plane wall.     

Estimation of aerodynamic noise generated by forced compressible round jetsEstimation du bruit d'origine aérodynamique rayonné par des jets ronds forcés compressibles

An acoustic numerical code based on Ligthill's analogy is combined with large-eddy simulations techniques in order to evaluate the noise emitted by subsonic $(M=0.7)$ and supersonic $(M=1.4)$ round jets. We show first that, for centerline Mach number $M=0.9$ and Reynolds number $\mathit{Re}=3.6\times {10}^3$, acoustic intensities compare satisfactorily with experimental data of the literature in terms of levels and directivity. Afterwards, high Reynolds number $(\mathit{Re}=3.6\times {10}^4)$ free and forced jets at Mach 0.7 and 1.4 are studied. Numerical results show that the jet noise intensity depends on the nature of the upstream mixing layer. Indeed, the subsonic jet is 4 dB quieter than the free jet when acting on this shear layer by superposing inlet varicose and flapping perturbations at preferred and first subharmonic frequency, respectively. The maximal acoustic level of the supersonic jet is, on the other hand, 3 dB lower than the free one with a flapping upstream perturbation at the second subharmonic. The results reported in this paper confirm previous works presented in the literature demonstrating that jet noise may be modified according to the inlet conditions. To cite this article: M. Maidi, C. R. Mecanique 334 (2006).     

Poroacoustic solitary waves under the unidirectional Darcy–Jordan model

Working in the context of poroacoustics, we present new, physically relevant, explicit solutions to the Cauchy problem for the model we term the (1D) damped Riemann equation. The solitary waveforms that evolve from both Lorentzian ($C^{\infty }$-smooth) and symmetric-exponential ($C^0$-smooth) initial conditions are analyzed, the focus being on wave overturning and the evolution/structure of the shocks which develop thereafter. In addition to those for both the multi- and single-valued forms of each solution, expressions for the shock amplitude, velocity, and critical values of the physical parameters are derived/compared. Lastly, links to other areas of continuum physics, and possible follow-on investigations, are noted.     

Large eddy simulation of serration effects on an ultra-high-bypass-ratio engine exhaust jet

Serrated jet nozzles are considered to be an efficient and practical passive control approach for jet noise. However, some fundamental mechanisms of serration effects on jet noise are not fully understood, especially in terms of the sound source. In this paper, a high-fidelity simulation framework using large-eddy simulation (LES) is demonstrated to predict near-field turbulence and far-field acoustics from an ultra-high-bypass-ratio engine with round and serrated nozzles. Far-field sound is predicted using Ffowcs Willams–Hawkings (FWH) integration. The results show that the serrated nozzle increases mixing near the nozzle and hence the turbulence decay rate, reducing the turbulence level downstream. The serrations shift the energy from the low frequencies to the high frequencies and decrease overall sound pressure levels by about 3 dB over the low-frequency range. Sound sources are analysed based on fourth-order space–time correlations. There are six major source components ($R_{1111}$, $R_{2222}$, $R_{3333}$, $R_{1313}$, $R_{1212}$, and $R_{2323}$) inside the jet shear layers. The serrations are able to reduce the amplitude of these source terms, causing them to decay rapidly to a level below the round nozzle jet within 2D downstream of the nozzle.     

Experimental investigations on the dynamic behavior of a 2-DOF airfoil in the transitional Re number regime based on digital-image correlation measurements

The present paper investigates the fluid–structure interaction (FSI) of a wing with two degrees of freedom (DOF), i.e., pitch and heave, in the transitional Reynolds number regime. This 2-DOF setup marks a classic configuration in aeroelasticity to demonstrate flutter stability of wings. In the past, mainly analytic approaches have been developed to investigate this challenging problem under simplifying assumptions such as potential flow. Although the classical theory offers satisfying results for certain cases, modern numerical simulations based on fully coupled approaches, which are more generally applicable and powerful, are still rarely found. Thus, the aim of this paper is to provide appropriate experimental reference data for well-defined configurations under clear operating conditions. In a follow-up contribution these will be used to demonstrate the capability of modern simulation techniques to capture instantaneous physical phenomena such as flutter. The measurements in a wind tunnel are carried out based on digital-image correlation (DIC). The investigated setup consists of a straight wing using a symmetric NACA 0012 airfoil. For the experiments the model is mounted into a frame by means of bending and torsional springs imitating the elastic behavior of the wing. Three different configurations of the wing possessing a fixed elastic axis are considered. For this purpose, the center of gravity is shifted along the chord line of the airfoil influencing the flutter stability of the setup. Still air free-oscillation tests are used to determine characteristic properties of the unloaded system (e.g. mass moment of inertia and damping ratios) for one (pitch or heave) and two degrees (pitch and heave) of freedom. The investigations on the coupled 2-DOF system in the wind tunnel are performed in an overall chord Reynolds number range of $9.66\times 10^3\le \text{Re}\le 8.77\times 10^4$. The effect of the fluid-load induced damping is studied for the three configurations. Furthermore, the cases of limit-cycle oscillation (LCO) as well as diverging flutter motion of the wing are characterized in detail. In addition to the DIC measurements, hot-film measurements of the wake flow for the rigid and the oscillating airfoil are presented in order to distinguish effects originating from the flow and the structure.     

Instabilité non linéaire dans un milieu en expansionNonlinear instability in a expanding medium

We consider nonlinear acoustical phenomena, explosive instabilities and a formation of localized structures in nonstationary environment. An example of such a medium is our Universe in expansion considered as a fluid submissive to a gravitational self-concorded force field and governed by the classical hydrodynamics equations. We show that the taking into account of the nonlinear effects allow us to understand the causes of the appearance of the specific nonlinear instability, which is calling explosive instability. This type of instability is more fast, $～\ln [{(t_0?t)}^{\mathrm{?1}}]$ for density fluctuation, that the habitual instability (exponential, $～{\mathrm{e}}^{\gamma t}$): at the end of a finite time, all spatial inhomogeneity of the initials conditions lead to a formation of singularities in the fields. This phenomena will be appear if certains conditions for the initials amplitudes and wavelengths of the fluctuations are observed. To cite this article: F. Henon, V. Pavlov, C. R. Mecanique 334 (2006).     

2D numerical contributions for the study of non-cohesive sediment transport beneath tidal bores

2D numerical simulations of tidal bores were obtained using the OpenFOAM CFD software to solve the Navier–Stokes equations by means of the Finite Volume Method by applying a LES turbulence model. The trajectories of non-cohesive sediment particles beneath tidal bores were estimated using a tracker method. Using the fourth order Runge–Kutta scheme, the tracker method solves the Maxey and Riley equations, which requires the knowledge of the velocity field at time t. From 2D numerical simulations of tidal bores, we proposed a classification of tidal bores with respect to the Froude number Fr (or r the ratio of water depths). For a Froude number $1<Fr<1.43$ ($1<r<1.57$), the tidal bore is undular. For a Froude number $1.43<Fr<1.57$ ($1.57<r<1.75$), the tidal bore is partially breaking, which is similar to the transitional tidal bore defined by Furgerot (2014). And for a Froude number $Fr>1.57$ ($r>1.75$), the tidal bore is totally breaking. The numerical results of trajectories of non-cohesive sediment particles are similar to the type of trajectories given by the analytical model proposed by Chen et al. (2012) with some modifications to take into account the effects of gravity, elevation, and attenuation. The parameters of modified Chen's model, ${\beta }_1$, ${\beta }_2$ and ${\beta }_3$, are linearly proportional to the Froude number Fr. This is because the level of turbulence for undular tidal bores is low. The flow induced by an undular tidal bore is not complex. This physical phenomenon is quasi linear. The parameter ${\beta }_1$, related to the front celerity of the undular tidal bore, decreases when the Froude number Fr increases. The parameter ${\beta }_2$, related to the elevation, increases when the Froude number Fr increases. And the parameter ${\beta }_3$, related to the attenuation of the secondary waves, increases when the Froude number Fr increases.     

Response and wake patterns of two side-by-side elastically supported circular cylinders in uniform laminar cross-flow

Vortex-induced vibrations (VIV) of two side-by-side elastically supported circular cylinders in a uniform flow with the Reynolds number of 100 are numerically investigated by using the immersed boundary method. The cylinders are constrained to oscillate in the cross-flow direction with a center-to-center spacing ratio $T/D$ ranging from 2 to 5. The structural damping is set to zero to enable large vibration amplitudes in the range of reduced velocity $U_r=3-10$. It is found that the proximity of the cylinders does not have a significant impact to the lock-in region and cylinder responses, except at a small spacing ratio of $T/D=2$. The critical spacing ratio is determined as $T/D=4$ and beyond that the interaction between the cylinders is negligible. The following six near-wake patterns are observed; the irregular pattern, in-phase flip-flopping pattern, out-of-phase flip-flopping pattern, in-phase-synchronized pattern, antiphase-synchronized pattern and the biased antiphase-synchronized pattern. These patterns are plotted in a plane of $U_r$ and $T/D$, together with approximate borderlines to distinguish one region from the others. The time histories, spectral features and wavelet transform contours of drag and lift forces are presented to elucidate the mechanisms of the in-phase and out-of-phase flip-flopping phenomena. It is established that the in-phase flip-flopping stems from the long-short near-wake pattern and its low-frequency flip-over, whereas the out-of-phase pattern originates from the large vortex shedding from the fictitious bluff-body with an augmented characteristic length.