High-order LES modeling of turbulent incompressible flow

This article presents the high-order algorithms that we have developed for large-eddy simulation of incompressible flows, and the results that have been obtained for the 3D turbulent wake of a cylinder at a Reynolds number of $\mathit{Re}=3900$. To cite this article: R. Pasquetti, C. R. Mecanique 333 (2005).     

On the slow motion of a cluster of bubbles under the combined action of gravity and thermocapillarity

The slow migration of N spherical bubbles under combined buoyancy and thermocapillarity effects is investigated by appealing solely to $3N+1$ boundary-integral equations. In addition to the theory and the associated implementation strategy, preliminary numerical results are both presented and discussed for a few clusters involving 2, 3, 4 or 5 bubbles with a special attention paid to the case of rigid configurations. To cite this article: A. Sellier, C. R. Mecanique 333 (2005).     

Stabilité de l'écoulement de Hartmann chauffé par le basStability of the Hartmann flow heated from below

A numerical study is conducted in order to determine the influence of a vertical magnetic field, the Reynolds number and a temperature stratification on the instabilities occurring in the Hartmann flow heated from below. For $\mathit{Pr}=0.001$ and $\mathit{Ha}?2.5$, the results show that the vertical magnetic field has a stabilizing effect on both transverse oscillatory travelling waves $(T)$ and longitudinal stationary rolls $(L)$. The temperature stratification is responsible of a destabilization of the transverse $(T)$ modes and the appearance of longitudinal $(L)$ modes non-existent for the isothermal Hartmann flow. Moreover, the extent of the domains of Re where the transverse modes $(T)$ prevail is found to narrow when Ha increases and to widen when Ra increases for a given value of Ha. On the other hand, for the $(L)$ modes, the extent of the domains of Re where they prevail increases when Ha grows. To cite this article: W. Fakhfakh et al., C. R. Mecanique 334 (2006).     

Instability of Navier slip flow of liquidsInstabilité des écoulements glissant de liquides.

We investigate the stability problem related to the basic slip flows of liquids in plane microchannels by using the Navier slip concept. We found that if the Navier slip parameter ($N_s$) equals 0.06, the critical Reynolds number (${\mathit{Re}}_{\mathrm{cr}}$) becomes 213.6. There are short-wave instabilities, however, when we further increase $N_s$ to 0.07 or 0.08. ${\mathit{Re}}_{\mathrm{cr}}$ becomes 132.9 for $N_s=0.08$ if we neglect the short-wave instability. To cite this article: A.K.-H. Chu, C. R. Mecanique 332 (2004).     

Modélisation et simulation numérique du changement de phase liquide–vapeur en cavitéModeling and numerical simulation of liquid–vapor phase change in an enclosed cavity

A model for the simulation of boiling flow with phase change in a closed cavity is presented. A front-tracking method is used to deal with the liquid–vapor interface. The liquid phase is incompressible while the vapor phase is weakly compressible and obeys to the perfect gas law. This model can deal with large density ratio (${\rho }_l/{\rho }_v?1000$) flows while accounting for the saturation curve. Computations are performed on a 1D validation case, idealizing a pressure cooker. Results are compared with a low Mach number approximation. To cite this article: V. Daru et al., C. R. Mecanique 334 (2006).     

Déformation cumulée tensorielle dans le référentiel en rotation logarithmiqueCumulated tensorial strain measure in the logarithmic rotating frame.

The time integration of strain rate tensor $\stackrel{?}{\mathit{D}}$ is a central problem in large transformations even if it is often an underlying one. The cumulated tensorial strains, obtained by the time integration of strain rate tensor $\stackrel{?}{\mathit{D}}$, allow the tackling of this problem from a geometrical point of view, and independently of material behaviour considerations. The time integration here takes place in the local objective frame defined by the logarithmic spin proposed by Lehmann et al. and Xiao et al. The numerical results obtained in a closed deformation path are presented here. The advantages and drawbacks of this novel integration for the development of behaviour laws are described. To cite this article: V. Mora et al., C. R. Mecanique 332 (2004).     

Settling motion of interacting solid particles in the vicinity of a plane solid boundary

The sedimentation of $N?1$ small arbitrarily-shaped solid bodies near a solid plane is addressed by discarding inertial effects and using 6N boundary-integral equations. Numerical results for 2 or 3 identical spheres reveal that combined wall–particle and particle–particle interactions deeply depend on the cluster's geometry and distance to the wall and may even cancel for a sphere which then moves as it were isolated. To cite this article: A. Sellier, C. R. Mecanique 333 (2005).     

Asymptotic analysis of layered elastic beams

We consider an elastic beam formed by three layers, fixed at one end and loaded at the free end. We call adherents the upper and lower layers ${\Omega }_{+}^{?}$ and ${\Omega }_{?}^{?}$ and an adhesive layer ${\Omega }_m^{?}$. We denote by $?h_{\pm ,m}$ the thickness of each layer and we suppose that the stiffness of the adhesive layer is ${?}^2$, with respect to that of the adherents. By an asymptotic analysis we obtain the zeroth order limit problem and the form of the second order displacements. To cite this article: M. Serpilli, C. R. Mecanique 333 (2005).     

Control of three-dimensional wakes using evolution strategies

We investigate three-dimensional cylinder wakes of incompressible fully developed flows at $\mathit{Re}=300$, resulting from control induced by tangential motions of the cylinder surface. The motion of the cylinder surface, in two dimensions, is optimized using evolution strategies, resulting in significant drag reduction and drastic modification of the wake as compared to the uncontrolled flow. The quasi-optimal velocity profile obtained in 2D is modified by spanwise harmonics and applied to 3D flows. The results indicate important differences in the flow physics induced by two and three dimensional control strategies. To cite this article: P. Poncet et al., C. R. Mecanique 333 (2005).     

A symmetry-preserving Cartesian grid method for computing a viscous flow past a circular cylinder

This article deals with a numerical method for solving the unsteady, incompressible Navier–Stokes equations in domains with arbitrarily-shaped boundaries, where the boundary is represented using the Cartesian grid approach. We introduce a novel cut-cell discretization which preserves the spectral properties of convection and diffusion. Here, convection is discretized by a skew-symmetric operator and diffusion is approximated by a symmetric, positive-definite coefficient matrix. Such a symmetry-preserving discretization conserves the kinetic energy (if the dissipation is turned off) and is stable on any grid. The method is successfully tested for an incompressible, unsteady flow around a circular cylinder at $\mathit{Re}=100$. To cite this article: R. Verstappen, M. Dröge, C. R. Mecanique 333 (2005).     

Identification of the convective instability in a multi-component solution by 3D simulations

Three-dimensional calculations have been done to simulate the onset of convective motion in ternary nondilute solution under phase transition conditions. The process is considered for Rayleigh number in the range $[1\times {10}^3,1.4\times {10}^4]$, where subcritical convective motion with hexagonal flow pattern is identified. The results are in good agreement with the linear and finite amplitude theory of hydrodynamics instability. To cite this article: V.V. Kolmychkov et al., C. R. Mecanique 333 (2005).     

A fourth order accurate finite volume scheme for numerical simulations of turbulent Rayleigh–Bénard convection in cylindrical containers

A finite volume scheme, which is based on fourth order accurate central differences in spatial directions and on a hybrid explicit/semi-implicit time stepping scheme, was developed to solve the incompressible Navier–Stokes and energy equations on cylindrical staggered grids. This includes a new fourth order accurate discretization of the velocity and temperature fields at the singularity of the cylindrical coordinate system and a new stability condition [J. Appl. Numer. Anal. Comput. Math. 1 (2004) 315–326]. The method was applied in direct numerical simulations of turbulent Rayleigh–Bénard convection for different Rayleigh numbers $\mathit{Ra}={10}^{\gamma }$, $\gamma =5\text{,}\dots \text{,}8$, in wide cylinders with the aspect ratios $a\equiv H/R=0.2$ and $a=0.4$ (where R denotes the radius and H – the height of the cylinder). To cite this article: O. Shishkina, C. Wagner, C. R. Mecanique 333 (2005).     

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

Contrôle par rotation ou par aspiration de l'écoulement autour d'un cylindre calculé par Simulation des Grandes Échelles

Control of the flow around a circular cylinder is studied using Large Eddy Simulation. The influence of control by rotation and suction on the flow characteristics is considered for several Reynolds numbers. Comparisons with experiments were conducted at $\mathit{Re}={10}^5$ for the flow with and without control. A drag reduction up to 30% is obtained for an usual suction intensity. To cite this article: G. Fournier et al., C. R. Mecanique 333 (2005).     

A method for the calculation of nonsymmetric steady periodic capillary–gravity waves on water of arbitrary uniform depthUne méthode de calcul des ondes de gravité–capillarité périodiques et non-symétriques en profondeur arbitraire

In this Note the method developed by Aider and Debiane (2004) for the calculation of nonsymmetric water waves on infinite depth is extended to finite depth. The water-wave problem is reduced to a system of nonlinear algebraic equations which is solved by using Newton's method. Solutions are computed up to their limiting forms by decrementing the depth from the infinity to a value of the depth-wavelength ratio $h/\lambda$ less than 0.025. It is found that the waves become symmetric when the depth becomes very small. Relations giving some integral properties are derived. To cite this article: R. Aider, M. Debiane, C. R. Mecanique 334 (2006).     

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

Cut-off analysis of coherent vortical structure identification in a three-dimensional external flow

Vortical structure identification has more recently been applied in the study of the transport of vortical structures in low Reynolds number three-dimensional complex geometry flows. An important issue in this identification procedure is to choose an appropriate cut-off value ${\lambda }_2$ which takes into consideration the finite precision vortex interfaces. This cut-off choice is studied in this Note and applied to an external flow around a curved cylinder. The vortex identification technique at different cut-off values is compared to the threshold of the vorticity field showing the efficiency of choosing the optimal tolerance gap. The computations are performed with a fully three-dimensional spectral/hp element method. To cite this article: A. Miliou et al., C. R. Mecanique 333 (2005).     

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

The two-dimensional problem of steady waves on water of finite depth: regimes without waves of small amplitude

The two-dimensional problem of steady waves on water of finite depth is considered without assumptions about periodicity and symmetry of waves. A new form of Bernoulli's equation is derived, and it involves a new bifurcation parameter which is the product of the Froude number μ and the rate of flow ω. The main result obtained from this equation is the absence of waves, having sufficiently small amplitude, provided $|\mu \omega |>1$. To cite this article: V. Kozlov, N. Kuznetsov, C. R. Mecanique 333 (2005).