Comparison of Different Numerical Interface Capturing Methods for the Simulation of Faraday Waves

Faraday instability is a classic problem that occurs due to the relative displacement of the interface that separates two immiscible fluids placed in a closed container under oscillating acceleration parallel to gravity. The interface deformation and the induced flow patterns of this two-phase flow are very complex and numerical simulations could allow a deeper understanding of the dynamics of these systems. Some tests have been performed to establish a reference solution, but further validation is needed in order to ensure the validity of these solutions. In this work, we compare some numerical solutions for the linear and nonlinear regimes using the phase field scheme with predictions obtained using different numerical schemes such as Front Tracking, Volume of Fluid, and Element-based Finite Volume Method. The results show that, in both linear and nonlinear regimes, some important differences in the prediction of the interface dynamics between the methods are observed, and the need to provide a reference numerical solution for future benchmarks is highlighted.

Faraday waves in soft elastic solids

Recent experiments have observed the emergence of standing waves at the free surface of elastic bodies attached to a rigid oscillating substrate and subjected to critical values of forcing frequency and amplitude. This phenomenon, known as Faraday instability, is now well understood for viscous fluids but surprisingly eluded any theoretical explanation for soft solids. Here, we characterize Faraday waves in soft incompressible slabs using the Floquet theory to study the onset of harmonic and subharmonic resonance eigenmodes. We consider a ground state corresponding to a finite homogeneous deformation of the elastic slab. We transform the incremental boundary value problem into an algebraic eigenvalue problem characterized by the three dimensionless parameters, that characterize the interplay of gravity, capillary and elastic waves. Remarkably, we found that Faraday instability in soft solids is characterized by a harmonic resonance in the physical range of the material parameters. This seminal result is in contrast to the subharmonic resonance that is known to characterize viscous fluids, and opens the path for using Faraday waves for a precise and robust experimental method that is able to distinguish solid-like from fluid-like responses of soft matter at different scales.

Effect of porous layer on the Faraday instability in viscous liquid

A linear stability analysis of a viscous liquid on a vertically oscillating porous plane is performed for infinitesimal disturbances of arbitrary wavenumbers. A time-dependent boundary value problem is derived and solved based on the Floquet theory along with the complex Fourier series expansion. Numerical results show that the Faraday instability is dominated by the subharmonic solution at high forcing frequency, but it responds harmonically at low forcing frequency. The unstable regions corresponding to both subharmonic and harmonic solutions enhance with the increasing value of permeability and yields a destabilizing effect on the Faraday instability. Further, the presence of porous layer makes faster the transition process from subharmonic instability to harmonic instability in the wavenumber regime. In addition, the first harmonic solution shrinks gradually and becomes an unstable island, and ultimately disappears from the neutral curve if the porous layer thickness is increased. In contrast, the first and second subharmonic solutions coalesce, and the onset of Faraday instability is dominated by the subharmonic solution. In a special case, the study of Faraday instability of a viscous liquid on a porous substrate can be replaced by a study of Faraday instability of a viscous liquid on a slippery substrate when the permeability of the porous substrate is very low. Further, the Faraday instability can be destabilized by introducing a slip effect at the bottom plane.