The Rayleigh–Taylor instability in a porous medium

The classical Rayleigh–Taylor instability occurs when a heavy fluid overlies a lighter one, and the two fluids are separated by a horizontal interface. The configuration is unstable, and a small perturbation to the interface grows with time. Here, we consider such an arrangement for planar flow, but in a porous medium governed by Darcy’s law. First, the fully saturated situation is considered, where the two horizontal fluids are separated by a sharp interface. A classical linearized theory is reviewed, and the nonlinear model is solved numerically. It is shown that the solution is ultimately limited in time by the formation of a curvature singularity at the interface. A partially saturated Boussinesq theory is then presented, and its linearized approximation predicts a stable interface that merely diffuses. Nonlinear Boussinesq theory, however, allows the growth of drips and bubbles at the interface. These structures develop with no apparent overturning at their heads, unlike the corresponding flow for two free fluids.

Reflection and transmission coefficient approximations for P, S1 and S2 waves in triclinic media

SUMMARY The conventional assumptions, in most published approximations of the reflection and transmission (R/T) coefficients of plane waves at a plane interface between two anisotropic half-spaces, confine their applications to weakly anisotropic and/or weak contrast models. We consider the horizontal interface enclosed by two triclinic half-spaces to approximate the R/T coefficients normalized by the vertical energy flux. The homogeneous background medium can be anisotropic with arbitrary symmetry to better simulate the strongly anisotropic media. The second-order approximations are proposed to accommodate the strong contrast interface. We also consider an isotropic background medium under the weak anisotropy assumption. The obtained approximations can be applied to P, S1 and S2 waves, except for the transmission coefficients between the S1 and S2 waves. The S-wave transmission coefficients are insensitive to the model parameter contrasts and predominately rely on the S-wave polarization directions in the half-spaces above and below the interface. The proposed approximations are tested numerically.

Symmetric double bubbles in the Grushin plane

We address the double bubble problem for the anisotropic Grushin perimeter Pα, α ≥ 0, and the Lebesgue measure in ℝ2, in the case of two equal volumes. We assume that the contact interface between the bubbles lies on either the vertical or the horizontal axis. We first prove existence of minimizers via the direct method by symmetrization arguments and then characterize them in terms of the given area by first variation techniques. Even though no regularity theory is available in this setting, we prove that angles at which minimal boundaries intersect satisfy the standard 120-degree rule up to a suitable change of coordinates. While for α = 0 the Grushin perimeter reduces to the Euclidean one and both minimizers coincide with the symmetric double bubble found in Foisy et al. [Pacific J. Math. 159 (1993) 47–59], for α = 1 vertical interface minimizers have Grushin perimeter strictly greater than horizontal interface minimizers. As the latter ones are obtained by translating and dilating the Grushin isoperimetric set found in Monti and Morbidelli [J. Geom. Anal. 14 (2004) 355–368], we conjecture that they solve the double bubble problem with no assumptions on the contact interface.

A multi-block finite difference method for seismic wave equation in auxiliary coordinate system with irregular fluid–solid interface

Purpose The purpose of this paper is to develop a new finite difference method for solving the seismic wave propagation in fluid-solid media, which can be described by the acoustic and viscoelastic wave equations for the fluid and solid parts, respectively. Design/methodology/approach In this paper, the authors introduced a coordinate transformation method for seismic wave simulation method. In the new method, the irregular fluid–solid interface is transformed into a horizontal interface. Then, a multi-block coordinate transformation method is proposed to mesh every layer to curved grids and transforms every interface to horizontal interface. Meanwhile, a variable grid size is used in different regions according to the shape and the velocity within each region. Finally, a Lebedev-standard staggered coupled grid scheme for curved grids is applied in the multi-block coordinate transformation method to reduce the computational cost. Findings The instability in the auxiliary coordinate system caused by the standard staggered grid scheme is resolved using a curved grid viscoelastic wave field separation strategy. Several numerical examples are solved using this new method. It has been shown that the new method is stable, efficient and highly accurate in solving the seismic wave equation defined on domain with irregular fluid–solid interface. Originality/value First, the irregular fluid–solid interface is transformed into a horizontal interface by using the coordinate transformation method. The conversion between pressures and stresses is easy to implement and adaptive to different irregular fluid–solid interface models, because the normal stress and shear stress vanish when the normal angle is 90° in the interface. Moreover, in the new method, the strong false artificial boundary reflection and instability caused by ladder-shaped grid discretion are resolved as well.

Flow distribution in diverging compound channels using improved independent subsection method

Experiments have been conducted in three diverging compound channels for different flow conditions to study the flow distribution in floodplain, upper and lower main channel. In a compound channel, vertical apparent shear exists on the interface between the upper main channel and the floodplain, which generally accelerates the flow on the floodplain and resists the flow in the upper main channel. In addition, a horizontal apparent shear stress also occurs on the interface between the upper and lower main channels, which generally accelerates the flow in the lower one and resists the flow in the upper one. Therefore, it is essential to consider the exchanges of momentum at both vertical and horizontal shear layer regions. In this paper, an attempt is made to improve the classical independent subsection method (ISM) to determine the magnitudes of flow and velocities in both upper and lower main channels. Four subsections are created in improved ISM according to the vertical and horizontal division lines that correspond to the vertical interface between the main channel and floodplain and the horizontal interface between upper and lower main channels respectively. The improved ISM consists in a set of four coupled 1D momentum equations (instead three equations of classical ISM) for subsections and a mass conservation equation for the total cross-section. The computed results show that the method is well capable of predicting the discharge distributions in the floodplain and main channel (both at upper and lower main channel).