Saikiran Rapaka

Non-modal growth of perturbations in density-driven convection in porous media

In the context of geologic sequestration of carbon dioxide in saline aquifers, much interest has been focused on the process of density-driven convection resulting from dissolution of CO 2 in brine in the underlying medium. Recent investigations have studied the time and length scales characteristic of the onset of convection based on the framework of linear stability theory. It is well known that the non-autonomous nature of the resulting matrix does not allow a normal mode analysis and previous researchers have either used a quasi-static approximation or solved the initial-value problem with arbitrary initial conditions. In this manuscript, we describe and use the recently developed non-modal stability theory to compute maximum amplifications possible, optimized over all possible initial perturbations. Non-modal stability theory also provides us with the structure of the most-amplified (or optimal) perturbations. We also present the details of three-dimensional spectral calculations of the governing equations. The results of the amplifications predicted by non-modal theory compare well to those obtained from the spectral calculations.

Flow patterns in the sedimentation of an elliptical particle

We study the dynamics of a single two-dimensional elliptical particle sedimenting in a Newtonian fluid using numerical simulations. The main emphasis in this work is to study the effect of boundaries on the flow patterns observed during sedimentation. The simulations were performed using a multi-block lattice Boltzmann method as well as a finite-element technique and the results are shown to be consistent. We have conducted a detailed study on the effects of density ratio, aspect ratio and the channel blockage ratio on the flow patterns during sedimentation. As the channel blockage ratio is varied, our results show that there are five distinct modes of sedimentation: oscillating, tumbling along the wall, vertical sedimentation, horizontal sedimentation and an inclined mode where the particle sediments with a non-trivial orientation to the vertical. The inclined mode is shown to form a smooth bridge between the vertical and horizontal modes of sedimentation. For narrow channels, the mode of sedimentation is found to be sensitively dependent on the initial orientation of the particle. We present a phase diagram showing the transitions between the various modes of sedimentation with changing blockage ratio of the channel.

Onset of convection over a transient base-state in anisotropic and layered porous media

The topic of density-driven convection in porous media has been the focus of many recent studies due to its relevance as a long-term trapping mechanism during geological sequestration of carbon dioxide. Most of these studies have addressed the problem in homogeneous and anisotropic permeability fields using linear-stability analysis, and relatively little attention has been paid to the analysis for heterogeneous systems. Previous investigators have reduced the governing equations to an initial-value problem and have analysed it either with a quasi-steady-state approximation model or using numerical integration with arbitrary initial perturbations. Recently, Rapaka et al . ( J. Fluid Mech. , vol. 609, 2008, pp. 285–303) used the idea of non-modal stability analysis to compute the maximum amplification of perturbations in this system, optimized over the entire space of initial perturbations. This technique is a mathematically rigorous extension of the traditional normal-mode analysis to non-normal and time-dependent problems. In this work, we extend this analysis to the important cases of anisotropic and layered porous media with a permeability variation in the vertical direction. The governing equations are linearized and reduced to a set of coupled ordinary differential equations of the initial-value type using the Galerkin technique. Non-modal stability analysis is used to compute the maximum growth of perturbations along with the optimal wavenumber leading to this growth. We show that unlike the solution of the initial-value problem, results obtained using non-modal analysis are insensitive to the choice of bottom boundary condition. For the anisotropic problem, the dependence of critical time and wavenumber on the anisotropy ratio was found to be in good agreement with theoretical scalings proposed by Ennis-King et al . ( Phys. Fluids , vol. 17, 2005, paper no. 084107). For heterogeneous systems, we show that uncertainty in the permeability field at low wavenumbers can influence the growth of perturbations. We use a Monte Carlo approach to compute the mean and standard deviation of the critical time for a sample permeability field. The results from theory are also compared with finite-volume simulations of the governing equations using fully heterogeneous porous media with strong layering. We show that the results from non-modal stability analysis match extremely well with those obtained from the simulations as long as the assumption of strong layering remains valid.