Long-time evolution of sequestered CO2_2 in porous media

CO$_2$ sequestration in subsurface reservoirs is important for limiting atmospheric CO$_2$ concentrations. However, a complete physical picture able to predict the structure developing within the porous medium is lacking. We investigate theoretically reactive transport in the long-time evolution of carbon in the brine-rock environment. As CO$_2$ is injected into a brine-rock environment, a carbonate-rich region is created amid brine. Within the carbonate-rich region minerals dissolve and migrate from regions of high concentration to low concentration, along with other dissolved carbonate species. This causes mineral precipitation at the interface between the two regions. We argue that precipitation in a small layer reduces diffusivity, and eventually causes mechanical trapping of the CO$_2$. Consequently, only a small fraction of the CO$_2$ is converted to solid mineral; the remainder either dissolves in water or is trapped in its original form. We also study the case of a pure CO$_2$ bubble surrounded by brine and suggest a mechanism that may lead to a carbonate-encrusted bubble due to structural diffusion

Scaling of a slope: the erosion of tilted landscapes

We formulate a stochastic equation to model the erosion of a surface with fixed inclination. Because the inclination imposes a preferred direction for material transport, the problem is intrinsically anisotropic. At zeroth order, the anisotropy manifests itself in a linear equation that predicts that the prefactor of the surface height-height correlations depends on direction. The first higher-order nonlinear contribution from the anisotropy is studied by applying the dynamic renormalization group. Assuming an inhomogeneous distribution of soil substrate that is modeled by a source of static noise, we estimate the scaling exponents at first order in \ep-expansion. These exponents also depend on direction. We compare these predictions with empirical measurements made from real landscapes and find good agreement. We propose that our anisotropic theory applies principally to small scales and that a previously proposed isotropic theory applies principally to larger scales. Lastly, by considering our model as a transport equation for a driven diffusive system, we construct scaling arguments for the size distribution of erosion ``events'' or ``avalanches.'' We derive a relationship between the exponents characterizing the surface anisotropy and the avalanche size distribution, and indicate how this result may be used to interpret previous findings of power-law size distributions in real submarine avalanches.Comment: 19 pages, includes 10 PS figures. J. Stat. Phys. (in press

Earth’s carbon cycle: A mathematical perspective

The carbon cycle represents metabolism at a global scale. When viewed through a mathematical lens, observational data suggest that the cycle exhibits an underlying mathematical structure. This review focuses on two types of emerging results: evidence of global dynamical coupling between life and the environment, and an understanding of the ways in which smaller-scale processes determine the strength of that coupling. Such insights are relevant not only to predicting future climate but also to understanding the long-term co-evolution of life and the environment.NASA Astrobiology Institute (NNA08CN84A)NASA Astrobiology Institute (NNA13AA90A)National Science Foundation (U.S.) (OCE-0930866)National Science Foundation (U.S.) (EAR-1338810

Intertwined Functions of Separase and Caspase in Cell Division and Programmed Cell Death.

Timely sister chromatid separation, promoted by separase, is essential for faithful chromosome segregation. Separase is a member of the CD clan of cysteine proteases, which also includes the pro-apoptotic enzymes known as caspases. We report a role for the C. elegans separase SEP-1, primarily known for its essential activity in cell division and cortical granule exocytosis, in developmentally programmed cell death when the predominant pro-apoptotic caspase CED-3 is compromised. Loss of SEP-1 results in extra surviving cells in a weak ced-3(-) mutant, and suppresses the embryonic lethality of a mutant defective for the apoptotic suppressor ced-9/Bcl-2 implicating SEP-1 in execution of apoptosis. We also report apparent non-apoptotic roles for CED-3 in promoting germ cell proliferation, meiotic chromosome disjunction, egg shell formation, and the normal rate of embryonic development. Moreover, loss of the soma-specific (CSP-3) and germline-specific (CSP-2) caspase inhibitors result in CED-3-dependent suppression of embryonic lethality and meiotic chromosome non-disjunction respectively, when separase function is compromised. Thus, while caspases and separases have evolved different substrate specificities associated with their specialized functions in apoptosis and cell division respectively, they appear to have retained the residual ability to participate in both processes, supporting the view that co-option of components in cell division may have led to the innovation of programmed cell suicide early in metazoan evolution

Two-dimensional hydrodynamic lattice-gas simulations of binary immiscible and ternary amphiphilic fluid flow through porous media

The behaviour of two dimensional binary and ternary amphiphilic fluids under flow conditions is investigated using a hydrodynamic lattice gas model. After the validation of the model in simple cases (Poiseuille flow, Darcy's law for single component fluids), attention is focussed on the properties of binary immiscible fluids in porous media. An extension of Darcy's law which explicitly admits a viscous coupling between the fluids is verified, and evidence of capillary effects are described. The influence of a third component, namely surfactant, is studied in the same context. Invasion simulations have also been performed. The effect of the applied force on the invasion process is reported. As the forcing level increases, the invasion process becomes faster and the residual oil saturation decreases. The introduction of surfactant in the invading phase during imbibition produces new phenomena, including emulsification and micellisation. At very low fluid forcing levels, this leads to the production of a low-resistance gel, which then slows down the progress of the invading fluid. At long times (beyond the water percolation threshold), the concentration of remaining oil within the porous medium is lowered by the action of surfactant, thus enhancing oil recovery. On the other hand, the introduction of surfactant in the invading phase during drainage simulations slows down the invasion process -- the invading fluid takes a more tortuous path to invade the porous medium -- and reduces the oil recovery (the residual oil saturation increases).Comment: 48 pages, 26 figures. Phys. Rev. E (in press

Interface Roughening in a Hydrodynamic Lattice-Gas Model with Surfactant

Using a hydrodynamic lattice-gas model, we study interface growth in a binary fluid with various concentrations of surfactant. We find that the interface is smoothed by small concentrations of surfactant, while microemulsion droplets form for large surfactant concentrations. To assist in determining the stability limits of the interface, we calculate the change in the roughness and growth exponents $\alpha$ and $\beta$ as a function of surfactant concentration along the interface.Comment: 4 pages with 4 embedded ps figures. Requires psfig.tex. Will appear in PRL 14 Oct 199

Simulating Three-Dimensional Hydrodynamics on a Cellular-Automata Machine

We demonstrate how three-dimensional fluid flow simulations can be carried out on the Cellular Automata Machine 8 (CAM-8), a special-purpose computer for cellular-automata computations. The principal algorithmic innovation is the use of a lattice-gas model with a 16-bit collision operator that is specially adapted to the machine architecture. It is shown how the collision rules can be optimized to obtain a low viscosity of the fluid. Predictions of the viscosity based on a Boltzmann approximation agree well with measurements of the viscosity made on CAM-8. Several test simulations of flows in simple geometries -- channels, pipes, and a cubic array of spheres -- are carried out. Measurements of average flux in these geometries compare well with theoretical predictions.Comment: 19 pages, REVTeX and epsf macros require

Lattice-Boltzmann Method for Non-Newtonian Fluid Flows

We study an ad hoc extension of the Lattice-Boltzmann method that allows the simulation of non-Newtonian fluids described by generalized Newtonian models. We extensively test the accuracy of the method for the case of shear-thinning and shear-thickening truncated power-law fluids in the parallel plate geometry, and show that the relative error compared to analytical solutions decays approximately linear with the lattice resolution. Finally, we also tested the method in the reentrant-flow geometry, in which the shear-rate is no-longer a scalar and the presence of two singular points requires high accuracy in order to obtain satisfactory resolution in the local stress near these points. In this geometry, we also found excellent agreement with the solutions obtained by standard finite-element methods, and the agreement improves with higher lattice resolution