Thin-film Rayleigh–Taylor instability in the presence of a deep periodic corrugated wall

Rayleigh–Taylor instability of a thin liquid film overlying a passive fluid is examined when the film is attached to a periodic wavy deep corrugated wall. A reduced-order long-wave model shows that the wavy wall enhances the instability toward rupture when the interface pattern is sub-harmonic to the wall pattern. An expression that approximates the growth constant of instability is obtained for any value of wall amplitude for the special case when the wall consists of two full waves and the interface consists of a full wave. Nonlinear computations of the interface evolution show that sliding is arrested by the wavy wall if a single liquid film residing over a passive fluid is considered but not necessarily when a bilayer sandwiched by a top wavy wall and bottom flat wall is considered. In the latter case interface tracking shows that primary and secondary troughs will evolve and subsequently slide along the flat wall due to symmetry-breaking. It is further shown that this sliding motion of the interface can ultimately be arrested by the top wavy wall, depending on the holdup of the fluids. In other words, there exists a critical value of the interface position beyond which the onset of the sliding motion is observed and below which the sliding is always arrested.

A review of one-phase Hele-Shaw flows and a level-set method for nonstandard configurations

The classical model for studying one-phase Hele-Shaw flows is based on a highly nonlinear moving boundary problem with the fluid velocity related to pressure gradients via a Darcy-type law. In a standard configuration with the Hele-Shaw cell made up of two flat stationary plates, the pressure is harmonic. Therefore, conformal mapping techniques and boundary integral methods can be readily applied to study the key interfacial dynamics, including the Saffman–Taylor instability and viscous fingering patterns. As well as providing a brief review of these key issues, we present a flexible numerical scheme for studying both the standard and nonstandard Hele-Shaw flows. Our method consists of using a modified finite-difference stencil in conjunction with the level-set method to solve the governing equation for pressure on complicated domains and track the location of the moving boundary. Simulations show that our method is capable of reproducing the distinctive morphological features of the Saffman–Taylor instability on a uniform computational grid. By making straightforward adjustments, we show how our scheme can easily be adapted to solve for a wide variety of nonstandard configurations, including cases where the gap between the plates is linearly tapered, the plates are separated in time, and the entire Hele-Shaw cell is rotated at a given angular velocity. doi:10.1017/S144618112100033X