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
SN - 1573-2703
VL - 121
SP - 57
EP - 83
JO - Journal of Engineering Mathematics
JF - Journal of Engineering Mathematics
IS - 1
ER -