skip to main content


Title: Two-dimensional stability analysis of finite-amplitude interfacial gravity waves in a two-layer fluid
We explore basic mechanisms for the instability of finite-amplitude interfacial gravity waves through a two-dimensional linear stability analysis of the periodic and irrotational plane motion of the interface between two unbounded homogeneous fluids of different density in the absence of background currents. The flow domains are conformally mapped into two half-planes, where the time-varying interface is always mapped onto the real axis. This unsteady conformal mapping technique with a suitable representation of the interface reduces the linear stability problem to a generalized eigenvalue problem, and allows us to accurately compute the growth rates of unstable disturbances superimposed on steady waves for a wide range of parameters. Numerical results show that the wave-induced Kelvin–Helmholtz (KH) instability due to the tangential velocity jump across the interface produced by the steady waves is the major instability mechanism. Any disturbances whose dominant wavenumbers are greater than a critical value grow exponentially. This critical wavenumber that depends on the steady wave steepness and the density ratio can be approximated by a local KH stability analysis, where the spatial variation of the wave-induced currents is neglected. It is shown, however, that the growth rates need to be found numerically with care and the successive collisions of eigenvalues result in the wave-induced KH instability. In addition, the present study extends the previous results for the small-wavenumber instability, such as modulational instability, of relatively small-amplitude steady waves to finite-amplitude ones.  more » « less
Award ID(s):
2108524
NSF-PAR ID:
10345098
Author(s) / Creator(s):
;
Date Published:
Journal Name:
Journal of Fluid Mechanics
Volume:
938
ISSN:
0022-1120
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. We examine the linear stability of a homogeneous gas–solid suspension of small Stokes number particles, with a moderate mass loading, subject to a simple shear flow. The modulation of the gravitational force exerted on the suspension, due to preferential concentration of particles in regions of low vorticity, in response to an imposed velocity perturbation, can lead to an algebraic instability. Since the fastest growing modes have wavelengths small compared with the characteristic length scale ( $U_{g}/{\it\Gamma}$ ) and oscillate with frequencies large compared with ${\it\Gamma}$ , $U_{g}$ being the settling velocity and ${\it\Gamma}$ the shear rate, we apply the WKB method, a multiple scale technique. This analysis reveals the existence of a number density mode which travels due to the settling of the particles and a momentum mode which travels due to the cross-streamline momentum transport caused by settling. These modes are coupled at a turning point which occurs when the wavevector is nearly horizontal and the most amplified perturbations are those in which a momentum wave upstream of the turning point creates a downstream number density wave. The particle number density perturbations reach a finite, but large amplitude that persists after the wave becomes aligned with the velocity gradient. The growth of the amplitude of particle concentration and fluid velocity disturbances is characterised as a function of the wavenumber and Reynolds number ( $\mathit{Re}=U_{g}^{2}/{\it\Gamma}{\it\nu}$ ) using both asymptotic theory and a numerical solution of the linearised equations. 
    more » « less
  2. ABSTRACT

    We present numerical simulation results for the propagation of Alfvén waves in the charge starvation regime. This is the regime where the plasma density is below the critical value required to supply the current for the wave. We analyse a conservative scenario where Alfvén waves pick up charges from the region where the charge density exceeds the critical value and advect them along at a high Lorentz factor. The system consisting of the Alfvén wave and charges being carried with it, which we call charge-carrying Alfvén wave (CC-AW), moves through a medium with small, but non-zero, plasma density. We find that the interaction between CC-AW and the stationary medium has a two-stream like instability which leads to the emergence of a strong electric field along the direction of the unperturbed magnetic field. The growth rate of this instability is of the order of the plasma frequency of the medium encountered by the CC-AW. Our numerical code follows the system for hundreds of wave periods. The numerical calculations suggest that the final strength of the electric field is of the order of a few per cent of the AW amplitude. Little radiation is produced by the sinusoidally oscillating currents associated with the instability during the linear growth phase. However, in the non-linear phase, the fluctuating current density produces strong EM radiation near the plasma frequency and limits the growth of the instability.

     
    more » « less
  3. We investigate the effect of inertial particles dispersed in a circular patch of finite radius on the stability of a two-dimensional Rankine vortex in semi-dilute dusty flows. Unlike the particle-free case where no unstable modes exist, we show that the feedback force from the particles triggers a novel instability. The mechanisms driving the instability are characterized using linear stability analysis for weakly inertial particles and further validated against Eulerian–Lagrangian simulations. We show that the particle-laden vortex is always unstable if the mass loading $M>0$ . Surprisingly, even non-inertial particles destabilize the vortex by a mechanism analogous to the centrifugal Rayleigh–Taylor instability in radially stratified vortex with density jump. We identify a critical mass loading above which an eigenmode $m$ becomes unstable. This critical mass loading drops to zero as $m$ increases. When particles are inertial, modes that fall below the critical mass loading become unstable, whereas modes above it remain unstable but with lower growth rates compared with the non-inertial case. Comparison with Eulerian–Lagrangian simulations shows that growth rates computed from simulations match well the theoretical predictions. Past the linear stage, we observe the emergence of high-wavenumber modes that turn into spiralling arms of concentrated particles emanating out of the core, while regions of particle-free flow are sucked inward. The vorticity field displays a similar pattern which leads to the breakdown of the initial Rankine structure. This novel instability for a dusty vortex highlights how the feedback force from the disperse phase can induce the breakdown of an otherwise resilient vortical structure. 
    more » « less
  4. null (Ed.)
    The Rayleigh–Taylor instability of a linear viscoelastic fluid overlying a passive gas is considered, where, under neutral conditions, the key dimensionless groups are the Bond number and the Weissenberg number. The branching behaviour upon instability to sinusoidal disturbances is determined by weak nonlinear analysis with the Bond number advanced from its critical value at neutral stability. It is shown that the solutions emanating from the critical state either branch out supercritically to steady waves at predictable wavelengths or break up subcritically with a wavelength having a single node. The nonlinear analysis leads to the counterintuitive observation that Rayleigh–Taylor instability of a viscoelastic fluid in a laterally unbounded layer must always result in saturated steady waves. The analysis also shows that the subcritical breakup in a viscoelastic fluid can only occur if the layer is laterally bounded below a critical horizontal width. If the special case of an infinitely deep viscoelastic layer is considered, a simple expression is obtained from which the transition between steady saturated waves and subcritical behaviour can be determined in terms of the leading dimensionless groups. This expression reveals that the supercritical saturation of the free surface is due to the influence of the normal elastic stresses, while the subcritical rupture of the free surface is attributed to the influence of capillary effects. In short, depending upon the magnitude of the scaled shear modulus, there exists a wavenumber at which a transition from saturated waves to subcritical breakup occurs. 
    more » « less
  5. - (Ed.)

    We study a thin, laterally confined heated liquid layer subjected to mechanical parametric forcing without gravity. In the absence of parametric forcing, the liquid layer exhibits the Marangoni instability, provided the temperature difference across the layer exceeds a threshold. This threshold varies with the perturbation wavenumber, according to a curve with two minima, which correspond to long- and short-wave instability modes. The most unstable mode depends on the lateral confinement of the liquid layer. In wide containers, the long-wave mode is typically observed, and this can lead to the formation of dry spots. We focus on this mode, as the short-wave mode is found to be unaffected by parametric forcing. We use linear stability analysis and nonlinear computations based on a reduced-order model to investigate how parametric forcing can prevent the formation of dry spots. At low forcing frequencies, the liquid film can be rendered linearly stable within a finite range of forcing amplitudes, which decreases with increasing frequency and ultimately disappears at a cutoff frequency. Outside this range, the flow becomes unstable to either the Marangoni instability (for small amplitudes) or the Faraday instability (for large amplitudes). At high frequencies, beyond the cutoff frequency, linear stabilization through parametric forcing is not possible. However, a nonlinear saturation mechanism, occurring at forcing amplitudes below the Faraday instability threshold, can greatly reduce the film surface deformation and therefore prevent dry spots. Although dry spots can also be avoided at larger forcing amplitudes, this comes at the expense of generating large-amplitude Faraday waves.

     
    more » « less