Miscible displacement fronts of shear thinning fluids inside rough fractures

The miscible displacement of a shear-thinning fluid by another of same rheological properties is studied experimentally in a transparent fracture by an optical technique imaging relative concentration distributions. The fracture walls have complementary self-affine geometries and are shifted laterally in the direction perpendicular to the mean flow velocity {\bf U} : the flow field is strongly channelized and macro dispersion controls the front structure for P\'{e}clet numbers above a few units. The global front width increases then linearly with time and reflects the velocity distribution between the different channels. In contrast, at the local scale, front spreading is similar to Taylor dispersion between plane parallel surfaces. Both dispersion mechanisms depend strongly on the fluid rheology which shifts from Newtonian to shear-thinning when the flow rate increases. In the latter domain, increasing the concentration enhances the global front width but reduces both Taylor dispersion (due to the flattening of the velocity profile in the gap of the fracture) and the size of medium scale front structures

Soil remediation using an electrokinetic method

Abstract Remediation of soils contaminated by heavy metals can be accomplished by subjecting the soil to an electric DC field. The heavy metals that are dissolved will move to either the cathode (mostly) or the anode, depending on their charges. The electrical field also induces electro-osmotic flow that can carry uncharged species. Precipitated and sorbed species will dissolve as the solution is depleted. The acid that is generated at the anode can considerably enhance the dissolution. The acid moves toward the cathode and generates a region with high concentration and mobility of the cationic metals. We present the theoretical basis of the method, the results of some remediation experiments, and the results of some numerical simulations that were used to interpret the experimental results. The remediation experiments are supple mented by determinations of the acid neutralization capacity as well as the sorption capacity and sorption distribution constants for the metal ions. These are essential data for the modelling and can be used to assess whether remediation could be accomplished within a reasonable length of time

Chemical Transport in a Fissured Bock: Verification of a Numerical Model

Numerical models for simulating chemical transport in fissured rocks constitute powerful tools for evaluating the acceptability of geological nuclear waste repositories. Due to the very long-term, high toxicity of some nuclear waste products, the models are required to predict, in certain cases, the spatial and temporal distribution of chemical concentration less than 0.001% of the concentration released from the repository. Whether numerical models can provide such accuracies is a major question addressed in the present work. To this end, we have verified a numerical model, TRUMP, which solves the advective diffusion equation in general three dimensions with or without decay and source terms. The method is based on an integrated finite-difference approach. The model was verified against known analytic solution of the one-dimensional advection-diffusion problem as well as the problem of advection-diffusion in a system of parallel fractures separated by spherical particles. The studies show that as long as the magnitude of advectance is equal to or less than that of conductance for the closed surface bounding any volume element in the region (that is, numerical Peclet number <2), the numerical method can indeed match the analytic solution within errors of ±10{sup -3} % or less. The realistic input parameters used in the sample calculations suggest that such a range of Peclet numbers is indeed likely to characterize deep groundwater systems in granitic and ancient argillaceous systems. Thus TRUMP in its present form does provide a viable tool for use in nuclear waste evaluation studies. A sensitivity analysis based on the analytic solution suggests that the errors in prediction introduced due to uncertainties in input parameters is likely to be larger than the computational inaccuracies introduced by the numerical model. Currently, a disadvantage in the TRUMP model is that the iterative method of solving the set of simultaneous equations is rather slow when time constants vary widely over the flow region. Although the iterative solution may be very desirable for large three-dimensional problems in order to minimize computer storage, it seems desirable to use a direct solver technique in conjunction with the mixed explicit-implicit approach whenever possible. work in this direction is in progress

A numerical study on differences in using Navier-Stokes and Reynolds equations for modeling the fluid flow and particle transport in single rock fractures with shear

The study on fluid flow and transport processes of rock fractures in most practical applications involves two fundamental issues: the validity of Reynolds equation for fluid flow (as most often assumed) and the effects of shear displacements on the magnitudes and anisotropy of the fluid flow velocity field. The reason for such concerns is that the impact of the surface roughness of rock fractures is still an unresolved challenging issue. The later has been systematically investigated with results showing that shear displacement plays a dominant role on evolutions of fluid velocity fields, for both magnitudes and anisotropy, but the former has not received examinations in details due to the numerical complexities involving solution of the Navier–Stokes (NS) equations and the representations of fracture geometry during shear. The objective of this paper aims to solve this problem through a FEM modeling effort. Applying the COMSOL Multiphysics code (FEM) and assuming a 2D problem, we consider the coupled hydromechanical effect of fracture geometry change due to shear on fluid flow (velocity patterns) and particle transport (streamline/velocity dispersion), using measured topographical data of natural rock fracture surfaces. The fluid flow in the vertical 2D cross-sections of single rock fractures was simulated by solving both the Navier–Stokes and the Reynolds equation, and the particle transport was predicted by the streamline particle tracking method with calculated flow velocity fields (vectors) from the flow simulations, obtaining results such as flow velocity profiles, total flow rates, particle travel time, breakthrough curves and the Péclet number, Pe, respectively. The results obtained using NS and Reynolds equations were compared to illustrate the degree of the validity of the Reynolds equation for general applications in practice since the later is mush more computationally efficient for solving large-scale problems. The flow simulation results show that both the total flow rate and the flow velocity fields in a rough rock fracture predicted by the NS equation were quite different from those predicted by the Reynolds equation. The results show that a roughly 5–10% overestimation on the flow rate is produced when the Reynolds equation is used, and the ideal parabolic velocity profiles defined by the local cubic law, when Reynolds equation is used, is no longer valid, especially when the roughness feature of the fracture surfaces changes with shear. These deviations of flow rate and flow velocity profiles across the fracture aperture have a significant impact on the particle transport behavior and the associated properties, such as the travel time and Péclet number. The deviations increase with increasing flow velocity and become more significant when fracture aperture geometry changes with shear