We carry out a linear stability analysis of the density driven, miscible flow in porous
media in the context of CO2 sequestration. The physical phenomenon is analogous to
the thermal convective instability in a fluid cooled uniformly from above or heated from
below, in a semi-infinite half-space. CO2 dissolution into the underlying brine leads to a
local density increase that results in a gravitational instability. The problem is cast into
the self-similar form to benefit from localized eigenfunctions of the streamwise diffusion
operator. In contrast to the flow in an unbounded domain, the current problem in the
semi-infinite domain yields a critical time for instability as well as a longwave cutoff. In
order to avoid the quasi-steady-state approximation (QSSA), a dominant mode approach
is used to highlight the early time instability mechanism. Comparison with the initial-
value-problem shows that this method gives exact results for both the critical time and
the longwave cutoff. We also show by comparing with the initial-value-problem that the
QSSA in the self-similar setting provides reasonably accurate solutions for all times.
The linear stability theory is extended to the nonlinear regime by very high accuracy
direct numerical simulations. The nonlinear simulations show good agreement with the
linear stability predictions for small times. Macroscopic fingers display intense nonlinear
interactions at inception. This gravitational instability significantly increases the rate of
CO2 dissolution into the brine in high permeability reservoirs.
for more details please see:
We are currently investigating the long term behavior of the convecting system. In particular we are interested in the statistical time averaged flux in a semi infinite convecting system. Some preliminary results have been presented at the Supri-C Affiliate meeting 2005 [talk]
I hope to provide some links to other interesting research on density driven convection in CO2 storage in saline aquifers: