2018 journal article
A SCALAR CONSERVATION LAW FOR PLUME MIGRATION IN CARBON SEQUESTRATION
SIAM JOURNAL ON APPLIED MATHEMATICS, 78(3), 1823–1841.

A quasi-linear hyperbolic partial differential equation with a discontinuous flux models geologic carbon dioxide (CO$_2$) migration and storage [M. Hesse, F. Orr, and H. Tchelepi, Fluid Mech., 611 (2008), pp. 35--60]. Two flux functions characterize the model, giving rise to flux discontinuities. One convex flux describes the invasion of the plume into pore space, and the other captures the flow as the plume leaves CO$_2$ bubbles behind, which are then trapped in the pore space. We investigate the method of characteristics, the structure of shock and rarefaction waves, and the result of binary wave interactions. The dual flux property introduces unexpected differences between the structure of these solutions and those of a scalar conservation law with a convex flux. During our analysis, we introduce a new construction of cross-hatch characteristics in regions of the space-time plane where the solution is constant, and there are two characteristic speeds. This construction is used to generalize the notion of the Lax entropy condition for admissible shocks and is crucial to continuing the propagation of a shock wave if its speed becomes characteristic.