TY - JOUR

T1 - Rayleigh-Taylor Instabilities in Miscible Fluids with Initially Piecewise Linear Density Profiles

AU - Cowell, Scott

AU - Kent, James

AU - Trevelyan, P.M.J.

N1 - Publisher Copyright:
© 2020, Springer Nature B.V.
Copyright:
Copyright 2020 Elsevier B.V., All rights reserved.

PY - 2020/3/19

Y1 - 2020/3/19

N2 - The stability of some simple density profiles in a vertically orientated two-dimensional porous medium is considered. The quasi-steady-state approximation is made so that the stability of the system can be approximated. As the profiles diffuse in time, the instantaneous growth rates evolve in time. For an initial step function density profile, the instantaneous growth rate was numerically found to decay like T- 1 / 2 for large times T, and the corresponding eigenfunctions scale with eωT where ω is a constant. For density profiles initially corresponding to a finite layer, the instantaneous growth rate eventually decayed like T- 1. This corresponds to an instability with algebraic growth, and the eigenfunctions scale with Tp (where p is a constant) for large time. For a species initially linearly distributed in a finite layer, when the concentration has an increasing gradient in the downwards direction, the stability of the system was similar to that found for a uniformly distributed finite layer. However, when the concentration had a decreasing gradient in the downwards direction, the growth rates remained constant for a long period time, but eventually decayed in the same way as found in a uniformly distributed finite layer, for very large times. Numerical simulations were performed to validate the predictions made by the linear stability analysis.

AB - The stability of some simple density profiles in a vertically orientated two-dimensional porous medium is considered. The quasi-steady-state approximation is made so that the stability of the system can be approximated. As the profiles diffuse in time, the instantaneous growth rates evolve in time. For an initial step function density profile, the instantaneous growth rate was numerically found to decay like T- 1 / 2 for large times T, and the corresponding eigenfunctions scale with eωT where ω is a constant. For density profiles initially corresponding to a finite layer, the instantaneous growth rate eventually decayed like T- 1. This corresponds to an instability with algebraic growth, and the eigenfunctions scale with Tp (where p is a constant) for large time. For a species initially linearly distributed in a finite layer, when the concentration has an increasing gradient in the downwards direction, the stability of the system was similar to that found for a uniformly distributed finite layer. However, when the concentration had a decreasing gradient in the downwards direction, the growth rates remained constant for a long period time, but eventually decayed in the same way as found in a uniformly distributed finite layer, for very large times. Numerical simulations were performed to validate the predictions made by the linear stability analysis.

KW - Buoyancy

KW - linear stability

KW - porous media

KW - Darcy's law

KW - non-linear simulations

KW - piecewise linear pro files

U2 - 10.1007/s10665-020-10039-6

DO - 10.1007/s10665-020-10039-6

M3 - Article

VL - 121

SP - 57

EP - 83

JO - Journal of Engineering Mathematics

JF - Journal of Engineering Mathematics

SN - 0022-0833

IS - 1

ER -