Linear Stability Analysis of Double-Diffusive Convection in Porous Media,with Application to Geological Storage of CO<Subscript>2</Subscript>

Onset of double-diffusive buoyancy-driven flow resulted from vertical temperature and concentration gradients in a horizontal layer of a saturated and homogenous porous medium is investigated using amplification factor theory. After injection of CO₂ into a deep saline aquifer, the density of the brine saturated with CO₂ increases slightly. This increase in density induces natural convection. The effect of geothermal gradient is also considered in this work as a second incentive for convection and the double-diffusion convection was studied. Linear stability analysis is used to predict the inception of instabilities and initial wavelength of the convective instabilities. The analysis presented is applied to acid gas injection (as an analogue for CO₂ storage) into saline aquifers in the Alberta basin. It is found that the geothermal gradient does not have significant effect on the onset of convection for these aquifers. It is shown that the geothermal effects on the onset of natural convection are negligible as compared to the solutal effects induced by dissolution and diffusion of CO₂ in deep saline aquifers. Therefore, the linear stability analysis and the long-term numerical simulation of CO₂ sequestration into such saline aquifers may be assumed to be isothermal in terms of natural convection occurrence.     

Numerical Simulation Studies of the Long-term Evolution of a CO<Subscript>2</Subscript> Plume in a Saline Aquifer with a Sloping Caprock

We have used the TOUGH2-MP/ECO2N code to perform numerical simulation studies of the long-term behavior of CO₂ stored in an aquifer with a sloping caprock. This problem is of great practical interest, and is very challenging due to the importance of multi-scale processes. We find that the mechanism of plume advance is different from what is seen in a forced immiscible displacement, such as gas injection into a water-saturated medium. Instead of pushing the water forward, the plume advances because the vertical pressure gradients within the plume are smaller than hydrostatic, causing the groundwater column to collapse ahead of the plume tip. Increased resistance to vertical flow of aqueous phase in anisotropic media leads to reduced speed of up-dip plume advancement. Vertical equilibrium models that ignore effects of vertical flow will overpredict the speed of plume advancement. The CO₂ plume becomes thinner as it advances, but the speed of advancement remains constant over the entire simulation period of up to 400 years, with migration distances of more than 80 km. Our simulations include dissolution of CO₂ into the aqueous phase and associated density increase, and molecular diffusion. However, no convection develops in the aqueous phase because it is suppressed by the relatively coarse (sub-) horizontal gridding required in a regional-scale model. A first crude sub-grid-scale model was developed to represent convective enhancement of CO₂ dissolution. This process is found to greatly reduce the thickness of the CO₂ plume, but, for the parameters used in our simulations, does not affect the speed of plume advancement.     

The Footprint of the CO<Subscript>2</Subscript> Plume during Carbon Dioxide Storage in Saline Aquifers: Storage Efficiency for Capillary Trapping at the Basin Scale

We study a sharp-interface mathematical model of CO₂ migration in deep saline aquifers, which accounts for gravity override, capillary trapping, natural groundwater flow, and the shape of the plume during the injection period. The model leads to a nonlinear advection–diffusion equation, where the diffusive term is due to buoyancy forces, not physical diffusion. For the case of interest in geological CO₂ storage, in which the mobility ratio is very unfavorable, the mathematical model can be simplified to a hyperbolic equation. We present a complete analytical solution to the hyperbolic model. The main outcome is a closed-form expression that predicts the ultimate footprint on the CO₂ plume, and the time scale required for complete trapping. The capillary trapping coefficient and the mobility ratio between CO₂ and brine emerge as the key parameters in the assessment of CO₂ storage in saline aquifers. Despite the many approximations, the model captures the essence of the flow dynamics and therefore reflects proper dependencies on the mobility ratio and the capillary trapping coefficient, which are basin-specific. The expressions derived here have applicability to capacity estimates by capillary trapping at the basin scale.     

Laboratory Flow Experiments for Visualizing Carbon Dioxide-Induced,Density-Driven Brine Convection

Injection of carbon dioxide (CO₂) into saline aquifers confined by low- permeability cap rock will result in a layer of CO₂ overlying the brine. Dissolution of CO₂ into the brine increases the brine density, resulting in an unstable situation in which more-dense brine overlies less-dense brine. This gravitational instability could give rise to density-driven convection of the fluid, which is a favorable process of practical interest for CO₂ storage security because it accelerates the transfer of buoyant CO₂ into the aqueous phase, where it is no longer subject to an upward buoyant drive. Laboratory flow visualization tests in transparent Hele-Shaw cells have been performed to elucidate the processes and rates of this CO₂ solute-driven convection (CSC). Upon introduction of CO₂ into the system, a layer of CO₂-laden brine forms at the CO₂-water interface. Subsequently, small convective fingers form, which coalesce, broaden, and penetrate into the test cell. Images and time-series data of finger lengths and wavelengths are presented. Observed CO₂ uptake of the convection system indicates that the CO₂ dissolution rate is approximately constant for each test and is far greater than expected for a diffusion-only scenario. Numerical simulations of our system show good agreement with the experiments for onset time of convection and advancement of convective fingers. There are differences as well, the most prominent being the absence of cell-scale convection in the numerical simulations. This cell-scale convection observed in the experiments may be an artifact of a small temperature gradient induced by the cell illumination.     

Modeling Gas Transport in the Shallow Subsurface During the ZERT CO<Subscript>2</Subscript> Release Test

We used the multiphase and multicomponent TOUGH2/EOS7CA model to carry out predictive simulations of CO₂ injection into the shallow subsurface of an agricultural field in Bozeman, Montana. The purpose of the simulations was to inform the choice of CO₂ injection rate and design of monitoring and detection activities for a CO₂ release experiment. The release experiment configuration consists of a long horizontal well (70 m) installed at a depth of approximately 2.5 m into which CO₂ is injected to mimic leakage from a geologic carbon sequestration site through a linear feature such as a fault. We estimated the permeability of the soil and cobble layers present at the site by manual inversion of measurements of soil CO₂ flux from a vertical-well CO₂ release. Based on these estimated permeability values, predictive simulations for the horizontal well showed that CO₂ injection just below the water table creates an effective gas-flow pathway through the saturated zone up to the unsaturated zone. Once in the unsaturated zone, CO₂ spreads out laterally within the cobble layer, where liquid saturation is relatively low. CO₂ also migrates upward into the soil layer through the capillary barrier and seeps out at the ground surface. The simulations predicted a breakthrough time of approximately two days for the 100kg d⁻¹ injection rate, which also produced a flux within the range desired for testing detection and monitoring approaches. The seepage area produced by the model was approximately five meters wide above the horizontal well, compatible with the detection and monitoring methods tested. For a given flow rate, gas-phase diffusion of CO₂ tends to dominate over advection near the ground surface, where the CO₂ concentration gradient is large, while advection dominates deeper in the system.     

Effect of Vertical Heterogeneity on Long-Term Migration of CO<Subscript>2</Subscript> in Saline Formations

For deep injection of CO₂ in thick saline formations, the movements of both the free gas phase and dissolved CO₂ are sensitive to variations in vertical permeability. A simple model for vertical heterogeneity was studied, consisting of a random distribution of horizontal impermeable barriers with a given overall volume fraction and distribution of lengths. Analytical results were obtained for the distribution of values for the permeability, and compared to numerical simulations of deep CO₂ injection and convection in heterogeneous formations, using multiple realizations for the permeability distribution. It is shown that for a formation of thickness H, the breakthrough times in two dimensions for deep injection scale as H ² for moderate injection rates. In comparison to heterogeneous shale distributions, a homogeneous medium with equivalent effective vertical permeability has a longer breakthrough time for deep injection, and a longer onset time for convection.     

Analysis on heat transfer correlations of supercritical CO<Subscript>2</Subscript> cooled in horizontal circular tubes

The objective of the present study is to analyze the heat transfer correlations of supercritical CO₂ cooled in horizontal circular tubes. In the paper, heat transfer correlations are first reviewed and compared with the experimental data at different heat fluxes. The results show that most of the previous correlations agree well with the experimental data under lower heat flux, but fail to predict the heat transfer coefficient well when the heat flux is as high as 33 kW/m². The study of buoyancy effect on convective heat transfer shows that buoyancy effect significantly affects the heat transfer with the increase of heat flux, and both free and forced convections operate in the turbulence flow during supercritical CO₂ cooling process. The influencing factors on heat transfer coefficient are summarized and the new correlation can be developed with the four dimensionless numbers.     

Modeling of CO<Subscript>2</Subscript> Leakage up Through an Abandoned Well from Deep Saline Aquifer to Shallow Fresh Groundwaters

This article presents a numerical modeling application using the code TOUGHREACT of a leakage scenario occurring during a CO₂ geological storage performed in the Jurassic Dogger formation in the Paris Basin. This geological formation has been intensively used for geothermal purposes and is now under consideration as a site for the French national program of reducing greenhouse gas emissions and CO₂ geological storage. Albian sandstone, situated above the Dogger limestone is a major strategic potable water aquifer; the impacts of leaking CO₂ due to potential integrity failure have, therefore, to be investigated. The present case–study illustrates both the capacity and the limitations of numerical tools to address such a critical issue. The physical and chemical processes simulated in this study have been restricted to: (i) supercritical CO₂ injection and storage within the Dogger reservoir aquifer, (ii) CO₂ upwards migration through the leakage zone represented as a 1D vertical porous medium to simulate the cement–rock formation interface in the abandoned well, and (iii) impacts on the Albian aquifer water quality in terms of chemical composition and the mineral phases representative of the porous rock by estimating fluid–rock interactions in both aquifers. Because of CPU time and memory constraints, approximation and simplification regarding the geometry of the geological structure, the mineralogical assemblages and the injection period (up to 5 years) have been applied to the system, resulting in limited analysis of the estimated impacts. The CO₂ migration rate and the quantity of CO₂ arriving as free gas and dissolving, firstly in the storage water and secondly in the water of the overlying aquifer, are calculated. CO₂ dissolution into the Dogger aquifer induces a pH drop from about 7.3 to 4.9 limited by calcite dissolution buffering. Glauconite present in the Albian aquifer also dissolves, causing an increase of the silicon and aluminum in solution and triggering the precipitation of kaolinite and quartz around the intrusion point. A sensitivity analysis of the leakage rate according to the location of the leaky well and the variability of the petro-physical properties of the reservoir, the leaky well zone and the Albian aquifers is also provided.     

Pore Scale Modeling of Reactive Transport Involved in Geologic CO<Subscript>2</Subscript> Sequestration

We apply a multi-component reactive transport lattice Boltzmann model developed in previous studies for modeling the injection of a CO₂-saturated brine into various porous media structures at temperatures T = 25 and 80°C. In the various cases considered the porous medium consists initially of calcite with varying grain size and shape. A chemical system consisting of Na⁺, Ca²⁺, Mg²⁺, H⁺, CO₂^°(aq){{\rm CO}_2^{\circ}{\rm (aq)}}, and Cl⁻ is considered. Flow and transport by advection and diffusion of aqueous species, combined with homogeneous reactions occurring in the bulk fluid, as well as the dissolution of calcite and precipitation of dolomite are simulated at the pore scale. The effects of the structure of the porous media on reactive transport are investigated. The results are compared with a continuum-scale model and the discrepancies between the pore- and continuum-scale models are discussed. This study sheds some light on the fundamental physics occurring at the pore scale for reactive transport involved in geologic CO₂ sequestration.     

Dissolution of CO2 From Leaking Fractures in Saline Formations

Leakage of CO₂ through fractures in saline formations will increase the CO₂??brine interface and promote CO₂ dissolution. We use a 2D, finite difference MATLAB model to simulate dissolution rates from a vertical fracture, with CO₂ flowing through it, in a secondary storage formation. The instigation of convection currents increases dissolution rates leading to higher dissolution in higher Rayleigh number systems. Comparison of our results with fracture flow rates shows that for typical fracture apertures dissolution from a fracture is small relative to the amount of CO₂ flowing through the fracture. Temporal and spatial variations in fracture permeability may reduce fracture flow rates and increase the relative amount of CO₂ dissolved from the fracture compared to the CO₂ flowing through the fracture. Further work on CO₂ dissolution in relation to fracture heterogeneity, flow of CO₂ within fractures and the interaction of multiple fractures will improve our ability to predict CO₂ dissolution rates for site characterisation.     

Study on heat transfer and pressure drop characteristics of internal heat exchangers in CO<Subscript>2</Subscript> system under cooling condition

In order to study the heat transfer and pressure drop on four types of internal heat exchangers (IHXs) of a CO₂ system, the experiment and numerical analysis were performed under a cooling condition. The configuration of the IHXs was a coaxial type and a micro-channel type. Two loops on the gas cooler part and the evaporator part were made, for experiment. And the section-by-section method and Hardy-Cross method were used for the numerical analysis. The capacity and pressure drop of the IHX are larger at the micro-channel type than at the coaxial type. When increasing the mass flow rate and the IHX length the capacity and pressure drop increase. The pressure drop of the evaporator loop is much larger than that of the gas cooler loop. The performance of the IHX was affected with operating condition of the gas-cooler and evaporator. The deviations between the experimental result and the numerical result are about ±20% for the micro-channel type and ±10% for the coaxial type. Thus, the new CO₂ heat transfer correlation should be developed to precisely predict a CO₂ heat transfer.     

An Experimental Study of CO<Subscript>2</Subscript> Exsolution and Relative Permeability Measurements During CO<Subscript>2</Subscript> Saturated Water Depressurization

Dissolution of CO₂ into brine is an important and favorable trapping mechanism for geologic storage of CO₂. There are scenarios, however, where dissolved CO₂ may migrate out of the storage reservoir. Under these conditions, CO₂ will exsolve from solution during depressurization of the brine, leading to the formation of separate phase CO₂. For example, a CO₂ sequestration system with a brine-permeable caprock may be favored to allow for pressure relief in the sequestration reservoir. In this case, CO₂-rich brine may be transported upwards along a pressure gradient caused by CO₂ injection. Here we conduct an experimental study of CO₂ exsolution to observe the behavior of exsolved gas under a wide range of depressurization. Exsolution experiments in highly permeable Berea sandstones and low permeability Mount Simon sandstones are presented. Using X-ray CT scanning, the evolution of gas phase CO₂ and its spatial distribution is observed. In addition, we measure relative permeability for exsolved CO₂ and water in sandstone rocks based on mass balances and continuous observation of the pressure drop across the core from 12.41 to 2.76 MPa. The results show that the minimum CO₂ saturation at which the exsolved CO₂ phase mobilization occurs is from 11.7 to 15.5%. Exsolved CO₂ is distributed uniformly in homogeneous rock samples with no statistical correlation between porosity and CO₂ saturation observed. No gravitational redistribution of exsolved CO₂ was observed after depressurization, even in the high permeability core. Significant differences exist between the exsolved CO₂ and water relative permeabilities, compared to relative permeabilities derived from steady-state drainage relative permeability measurements in the same cores. Specifically, very low CO₂ and water relative permeabilities are measured in the exsolution experiments, even when the CO₂ saturation is as high as 40%. The large relative permeability reduction in both the water and CO₂ phases is hypothesized to result from the presence of disconnected gas bubbles in this two-phase flow system. This feature is also thought to be favorable for storage security after CO₂ injection.     

Injection and Storage of CO<Subscript>2</Subscript> in Deep Saline Aquifers: Analytical Solution for CO<Subscript>2</Subscript> Plume Evolution During Injection

Injection of fluids into deep saline aquifers is practiced in several industrial activities, and is being considered as part of a possible mitigation strategy to reduce anthropogenic emissions of carbon dioxide into the atmosphere. Injection of CO₂ into deep saline aquifers involves CO₂ as a supercritical fluid that is less dense and less viscous than the resident formation water. These fluid properties lead to gravity override and possible viscous fingering. With relatively mild assumptions regarding fluid properties and displacement patterns, an analytical solution may be derived to describe the space–time evolution of the CO₂ plume. The solution uses arguments of energy minimization, and reduces to a simple radial form of the Buckley–Leverett solution for conditions of viscous domination. In order to test the applicability of the analytical solution to the CO₂ injection problem, we consider a wide range of subsurface conditions, characteristic of sedimentary basins around the world, that are expected to apply to possible CO₂ injection scenarios. For comparison, we run numerical simulations with an industry standard simulator, and show that the new analytical solution matches a full numerical solution for the entire range of CO₂ injection scenarios considered. The analytical solution provides a tool to estimate practical quantities associated with CO₂ injection, including maximum spatial extent of a plume and the shape of the overriding less-dense CO₂ front.     

Enhanced CO<Subscript>2</Subscript> Storage and Sequestration in Deep Saline Aquifers by Nanoparticles: Commingled Disposal of Depleted Uranium and CO<Subscript>2</Subscript>

Geological storage of anthropogenic CO₂ emissions in deep saline aquifers has recently received tremendous attention in the scientific literature. Injected buoyant CO₂ accumulates at the top part of the aquifer under a sealing cap rock. Potential buoyant movement of CO₂ has caused some concern that the high-pressure CO₂ could breach the seal rock. However, CO₂ will diffuse into the brine underneath and generate a slightly denser fluid that may induce instability and convective mixing. Onset times of instability and convective mixing performance depend on the physical properties of the rock and fluids, such as permeability and density contrast. We present the novel idea of adding nanoparticles (NPs) to injected CO₂ to increase density contrast between the CO₂-rich brine and the underlying resident brine and, consequently, decrease onset time of instability and increase convective mixing. The analyses show that 0.001 volume fraction of NPs added to the CO₂ stream shortens onset time of mixing by approximately 80% and increases convective mixing by 50%. If it thus originally takes 5 years for the overlying CO₂ to start convective mixing, by adding NPs, onset time of mixing reduces to 1 year, and after initiation of convective mixing, mixing improves by 50%. A reduction of the CO₂ leakage risk ensues. In addition to other metallic NPs, use of processed depleted uranium oxide (DU) as the NPs is also proposed. DU-NPs are potentially stable and might be safely commingled with CO₂ to store in saline aquifers.     

Methanol diffusion rates through the anode diffusion layer in direct methanol fuel cells from limiting current measurements

Metal foams may be used in direct methanol fuel cells to feed reactants to the catalyst layer and to collect current from the resulting electrochemical reaction. Although the mass transfer from the metal foam to the underlying gas diffusion layer (GDL) is diffusion-dominated, it is found that at a fixed methanol concentration, the limiting current density increases with increasing methanol flow rates. This unexpected result is attributed to the more efficient removal of product CO₂ from the GDL. A methodology is developed to estimate the effective diffusion coefficient of methanol in the anode diffusion layer from limiting current density measurements, and to extract the fraction of GDL volume occupied by CO₂.     

Modeling of Rocks and Cement Alteration due to CO<Subscript>2</Subscript> Injection in an Exploited Gas Reservoir

The injection of CO₂ in exploited natural gas reservoirs as a means to reduce greenhouse gas (GHG) emissions is highly attractive as it takes place in well-known geological structures of proven integrity with respect to gas leakage. The injection of a reactive gas such as CO₂ puts emphasis on the possible alteration of reservoir and caprock formations and especially of the wells’ cement sheaths induced by the modification of chemical equilibria. Such studies are important for injectivity assurance, wellbore integrity, and risk assessment required for CO₂ sequestration site qualification. Within a R&D project funded by Eni, we set up a numerical model to investigate the rock–cement alterations driven by the injection of CO₂ into a depleted sweet natural gas pool. The simulations are performed with the TOUGHREACT simulator (Xu et al. in Comput Geosci 32:145–165, 2006) coupled to the TMGAS EOS module (Battistelli and Marcolini in Int J Greenh Gas Control 3:481–493, 2009) developed for the TOUGH2 family of reservoir simulators (Pruess et al. in TOUGH2 User’s Guide, Version 2.0, 1999). On the basis of field data, the system is considered in isothermal (50°C) and isobaric (128.5 bar) conditions. The effects of the evolving reservoir gas composition are taken into account before, during, and after CO₂ injection. Fully water-saturated conditions were assumed for the cement sheath and caprock domains. The gas phase does not flow by advection from the reservoir into the interacting domains so that molecular diffusion in the aqueous phase is the most important process controlling the mass transport occurring in the system under study.     

Investigation of CO<Subscript>2</Subscript> Plume Behavior for a Large-Scale Pilot Test of Geologic Carbon Storage in a Saline Formation

The hydrodynamic behavior of carbon dioxide (CO₂) injected into a deep saline formation is investigated, focusing on trapping mechanisms that lead to CO₂ plume stabilization. A numerical model of the subsurface at a proposed power plant with CO₂ capture is developed to simulate a planned pilot test, in which 1,000,000 metric tons of CO₂ is injected over a 4-year period, and the subsequent evolution of the CO₂ plume for hundreds of years. Key measures are plume migration distance and the time evolution of the partitioning of CO₂ between dissolved, immobile free-phase, and mobile free-phase forms. Model results indicate that the injected CO₂ plume is effectively immobilized at 25 years. At that time, 38% of the CO₂ is in dissolved form, 59% is immobile free phase, and 3% is mobile free phase. The plume footprint is roughly elliptical, and extends much farther up-dip of the injection well than down-dip. The pressure increase extends far beyond the plume footprint, but the pressure response decreases rapidly with distance from the injection well, and decays rapidly in time once injection ceases. Sensitivity studies that were carried out to investigate the effect of poorly constrained model parameters permeability, permeability anisotropy, and residual CO₂ saturation indicate that small changes in properties can have a large impact on plume evolution, causing significant trade-offs between different trapping mechanisms.     

Characteristics of Salt-Precipitation and the Associated Pressure Build-Up during CO<Subscript>2</Subscript> Storage in Saline Aquifers

Mitigation and control of borehole pressure at the bottom of an injection well is directly related to the effective management of well injectivity during geologic carbon sequestration activity. Researchers have generally accepted the idea that high rates of CO₂ injection into low permeability strata results in increased bottom-hole pressure in a well. However, the results of this study suggested that this is not always the case, due to the occurrence of localized salt precipitation adjacent to the injection well. A series of numerical simulations indicated that in some cases, a low rate of CO₂ injection into high permeability formation induced greater pressure build-up. This occurred because of the different types of salt precipitation pattern controlled by buoyancy-driven CO₂ plume migration. The first type is non-localized salt precipitation, which is characterized by uniform salt precipitation within the dry-out zone. The second type, localized salt precipitation, is characterized by an abnormally high level of salt precipitation at the dry-out front. This localized salt precipitation acts as a barrier that hampers the propagation of both CO₂ and pressure to the far field as well as counter-flowing brine migration toward the injection well. These dynamic processes caused a drastic pressure build-up in the well, which decreased injectivity. By modeling a series of test cases, it was found that low-rate CO₂ injection into high permeability formation was likely to cause localized salt precipitation. Sensitivity studies revealed that brine salinity linearly affected the level of salt precipitation, and that vertical permeability enhanced the buoyancy effect which increased the growth of the salt barrier. The porosity also affected both the level of localized salt precipitation and dry-out zone extension depending on injection rates. High temperature injected CO₂ promoted the vertical movement of the CO₂ plume, which accelerated localized salt precipitation, but at the same time caused a decrease in the density of the injected CO₂. The combination of these two effects eventually decreased bottomhole pressure. Considering the injectivity degradation, a method is proposed for decreasing the pressure build-up and increasing injectivity by assigning a ‘skin zone’ that represents a local region with a transmissivity different from that of the surrounding aquifer.     

Simplified model for mass diffusivity estimation based on dynamic pendant drop volume analysis

We present a simplified correlation for calculating the dissolved gas moles in a pendant drop during the diffusion time, for several drop shapes. After this correlation is determined, the Yang and Gu (Ind Eng Chem Res 44:4474–4483, 2005) dynamic pendant drop volume analysis (DPDVA) method for calculation of mass diffusivity from the pendant drop volume variation against time can be used. We solved the differential equation in cylindrical coordinates for the mass transfer model of the gas diffusion into the liquid inside the pendant drop, using a different characteristic length (L_C), instead of the outer radius of the syringe needle (r_n) used in Yang and Gu (Ind Eng Chem Res 44:4474–4483, 2005) for defining the dimensionless variables. L_C is the relationship between the pendant drop volume and its mass transfer surface area at the initial conditions. The generalized correlation saves time, simplifies the method application and the deviations in the diffusion coefficient calculation respect to the complete Yang and Gu model are below 6%.     

Matched Boundary Extrapolation Solutions for CO<Subscript>2</Subscript> Well-Injection into a Saline Aquifer

The injection of supercritical CO₂ through wells into deep brine reservoirs is a topic of interest for geologic carbon sequestration. The injected CO₂ is predominantly immiscible with the brine and its low density relative to brine leads to strong buoyancy effects. The displacement of brine by CO₂ in general is a multidimensional, complex nonlinear problem that requires numerical methods to solve. The approximations of vertical equilibrium and complete gravity segregation (sharp interface) have been introduced to reduce the complexity and dimensionality of the problem. Furthermore, for the radial displacement process considered here, the problem can be formulated in terms of a similarity variable that reduces spatial and temporal dependencies to a single variable. However, the resulting ordinary differential equation is still nonlinear and exact solutions are not available. The existing analytical solutions are approximations limited to certain parameter ranges that become inaccurate over a large portion of the parameter space. Here, I use a matched boundary extrapolation method to provide much greater accuracy for analytical/semi-analytical approximations over the full parameter range.