No Arabic abstract
Flows forced by a precessional motion can exhibit instabilities of crucial importance, whether they concern the fuel of a flying object or the liquid core of a telluric planet. So far, stability analyses of these flows have focused on the special case of a resonant forcing. Here, we address the instability of the flow inside a precessing cylinder in the general case. We first show that the base flow forced by the cylinder precession is a superposition of a vertical or horizontal shear flow and an infinite sum of forced modes. We then perform a linear stability analysis of this base flow by considering its triadic resonance with two free Kelvin modes. Finally, we derive the amplitude equations of the free Kelvin modes and obtain an expression of the instability threshold and growth rate.
A cylinder undergoes precession when it rotates around its axis and this axis itself rotates around another direction. In a precessing cylinder full of fluid, a steady and axisymmetric component of the azimuthal flow is generally present. This component is called a zonal flow. Although zonal flows have been often observed in experiments and numerical simulations, their origin has eluded theoretical approaches so far. Here, we develop an asymptotic analysis to calculate the zonal flow forced in a resonant precessing cylinder, that is when the harmonic response is dominated by a single Kelvin mode. We find that the zonal flow originates from three different sources: (1) the nonlinear interaction of the inviscid Kelvin mode with its viscous correction; (2) the steady and axisymmetric response to the nonlinear interaction of the Kelvin mode with itself; and (3) the nonlinear interactions in the end boundary layers. In a precessing cylinder, two additional sources arise due to the equatorial Coriolis force and the forced shear flow. However, they cancel exactly. The study thus generalises to any Kelvin mode, forced by precession or any other mechanism. The present theoretical predictions of the zonal flow are confirmed by comparison with numerical simulations and experimental results. We also show numerically that the zonal flow is always retrograde in a resonant precessing cylinder (m=1) or when it results from resonant Kelvin modes of azimuthal wavenumbers m=2, 3, and presumably higher.
We analyze theoretically and experimentally the triadic resonance instability (TRI) of a plane inertial wave in a rotating fluid. Building on the classical triadic interaction equations between helical modes, we show by numerical integration that the maximum growth rate of the TRI is found for secondary waves that do not propagate in the same vertical plane as the primary wave (the rotation axis is parallel to the vertical). In the inviscid limit, we prove this result analytically, in which case the change in the horizontal propagation direction induced by the TRI evolves from $60^circ$ to $90^circ$ depending on the frequency of the primary wave. Thanks to a wave generator with a large spatial extension in the horizontal direction of invariance of the forced wave, we are able to report experimental evidence that the TRI of a plane inertial wave is three-dimensional. The wavevectors of the secondary waves produced by the TRI are shown to match the theoretical predictions based on the maximum growth rate criterion. These results reveal that the triadic resonant interactions between inertial waves are very efficient at redistributing energy in the horizontal plane, normal to the rotation axis.
We report on a new class of electromagnetically-driven fluid interface instability. Using the optical radiation pressure of a cw laser to bend a very soft near-critical liquid-liquid interface, we show that it becomes unstable for sufficiently large beam power P, leading to the formation of a stationary beam-centered liquid micro-jet. We explore the behavior of the instability onset by tuning the interface softness with temperature and varying the size of the exciting beam. The instability mechanism is experimentally demonstrated. It simply relies on total reflection of light at the deformed interface whose condition provides the universal scaling relation for the onset Ps of the instability.
We study the Richtmyer--Meshkov (RM) instability of a relativistic perfect fluid by means of high order numerical simulations with adaptive mesh refinement (AMR). The numerical scheme adopts a finite volume Weighted Essentially Non-Oscillatory (WENO) reconstruction to increase accuracy in space, a local space-time discontinuous Galerkin predictor method to obtain high order of accuracy in time and a high order one-step time update scheme together with a cell-by-cell space-time AMR strategy with time-accurate local time stepping. In this way, third order accurate (both in space and in time) numerical simulations of the RM instability are performed, spanning a wide parameter space. We present results both for the case in which a light fluid penetrates into a higher density one (Atwood number $A>0$), and for the case in which a heavy fluid penetrates into a lower density one (Atwood number $A<0$). We find that, for large Lorentz factors gamma_s of the incident shock wave, the relativistic RM instability is substantially weakened and ultimately suppressed. More specifically, the growth rate of the RM instability in the linear phase has a local maximum which occurs at a critical value of gamma_s ~ [1.2,2]. Moreover, we have also revealed a genuine relativistic effect, absent in Newtonian hydrodynamics, which arises in three dimensional configurations with a non-zero velocity component tangent to the incident shock front. In this case, the RM instability is strongly affected, typically resulting in less efficient mixing of the fluid.
We study fluid-structure interactions (FSIs) in a long and shallow microchannel, conveying a non-Newtonian fluid, at steady state. The microchannel has a linearly elastic and compliant top wall, while its three other walls are rigid. The fluid flowing inside the microchannel has a shear-dependent viscosity described by the power-law rheological model. We employ lubrication theory to solve for the flow problem inside the long and shallow microchannel. For the structural problem, we employ two plate theories, namely Kirchhoff-Love theory of thin plates and Reissner-Mindlin first-order shear deformation theory. The hydrodynamic pressure couples the flow and deformation problem by acting as a distributed load onto the soft top wall. Within our perturbative (lubrication theory) approach, we determine the relationship between flow rate and the pressure gradient, which is a nonlinear first-order ordinary differential equation for the pressure. From the solution of this differential equation, all other quantities of interest in non-Newtonian microchannel FSIs follow. Through illustrative examples, we show the effect of FSI coupling strength and the plate thickness on the pressure drop across the microchannel. Through direct numerical simulation of non-Newtonian microchannel FSIs using commercial computational engineering tools, we benchmark the prediction from our mathematical prediction for the flow rate-pressure drop relation and the structural deformation profile of the top wall. In doing so, we also establish the limits of applicability of our perturbative theory.