Numerically stable Monte Carlo-burnup-thermal hydraulic coupling schemes

This paper presents stochastic implicit coupling method intended for use in Monte-Carlo (MC) based reactor analysis systems that include burnup and thermal hydraulic (TH) feedbacks. Both feedbacks are essential for accurate modeling of advanced reactor designs and analyses of associated fuel cycles. In particular, we investigate the effect of different burnup-TH coupling schemes on the numerical stability and accuracy of coupled MC calculations. First, we present the beginning of time step method which is the most commonly used. The accuracy of this method depends on the time step length and it is only conditionally stable. This work demonstrates that even for relatively short time steps, this method can be numerically unstable. Namely, the spatial distribution of neutronic and thermal hydraulic parameters, such as nuclide densities and temperatures, exhibit oscillatory behavior. To address the numerical stability issue, new implicit stochastic methods are proposed. The methods solve the depletion and TH problems simultaneously and use under-relaxation to speed up convergence. These methods are numerically stable and accurate even for relatively large time steps and require less computation time than the existing methods. © 2013 Elsevier Ltd. All rights reserved.

A coupled hydro-mechanical analysis for prediction of hydraulic fracture propagation in saturated porous media using EFG mesh-less method

The details of the Element Free Galerkin (EFG) method are presented with the method being applied to a study on hydraulic fracturing initiation and propagation process in a saturated porous medium using coupled hydro-mechanical numerical modelling. In this EFG method, interpolation (approximation) is based on nodes without using elements and hence an arbitrary discrete fracture path can be modelled.The numerical approach is based upon solving two governing partial differential equations of equilibrium and continuity of pore water simultaneously. Displacement increment and pore water pressure increment are discretized using the same EFG shape functions. An incremental constrained Galerkin weak form is used to create the discrete system of equations and a fully implicit scheme is used for discretization in the time domain. Implementation of essential boundary conditions is based on the penalty method. In order to model discrete fractures, the so-called diffraction method is used.Examples are presented and the results are compared to some closed-form solutions and FEM approximations in order to demonstrate the validity of the developed model and its capabilities. The model is able to take the anisotropy and inhomogeneity of the material into account. The applicability of the model is examined by simulating hydraulic fracture initiation and propagation process from a borehole by injection of fluid. The maximum tensile strength criterion and Mohr-Coulomb shear criterion are used for modelling tensile and shear fracture, respectively. The model successfully simulates the leak-off of fluid from the fracture into the surrounding material. The results indicate the importance of pore fluid pressure in the initiation and propagation pattern of fracture in saturated soils. © 2013 Elsevier Ltd.

Methane contamination of drinking water accompanying gas-well drilling and hydraulic fracturing.

Directional drilling and hydraulic-fracturing technologies are dramatically increasing natural-gas extraction. In aquifers overlying the Marcellus and Utica shale formations of northeastern Pennsylvania and upstate New York, we document systematic evidence for methane contamination of drinking water associated with shale-gas extraction. In active gas-extraction areas (one or more gas wells within 1 km), average and maximum methane concentrations in drinking-water wells increased with proximity to the nearest gas well and were 19.2 and 64 mg CH(4) L(-1) (n = 26), a potential explosion hazard; in contrast, dissolved methane samples in neighboring nonextraction sites (no gas wells within 1 km) within similar geologic formations and hydrogeologic regimes averaged only 1.1 mg L(-1) (P < 0.05; n = 34). Average δ(13)C-CH(4) values of dissolved methane in shallow groundwater were significantly less negative for active than for nonactive sites (-37 ± 7‰ and -54 ± 11‰, respectively; P < 0.0001). These δ(13)C-CH(4) data, coupled with the ratios of methane-to-higher-chain hydrocarbons, and δ(2)H-CH(4) values, are consistent with deeper thermogenic methane sources such as the Marcellus and Utica shales at the active sites and matched gas geochemistry from gas wells nearby. In contrast, lower-concentration samples from shallow groundwater at nonactive sites had isotopic signatures reflecting a more biogenic or mixed biogenic/thermogenic methane source. We found no evidence for contamination of drinking-water samples with deep saline brines or fracturing fluids. We conclude that greater stewardship, data, and-possibly-regulation are needed to ensure the sustainable future of shale-gas extraction and to improve public confidence in its use.

Radium and barium removal through blending hydraulic fracturing fluids with acid mine drainage.

Wastewaters generated during hydraulic fracturing of the Marcellus Shale typically contain high concentrations of salts, naturally occurring radioactive material (NORM), and metals, such as barium, that pose environmental and public health risks upon inadequate treatment and disposal. In addition, fresh water scarcity in dry regions or during periods of drought could limit shale gas development. This paper explores the possibility of using alternative water sources and their impact on NORM levels through blending acid mine drainage (AMD) effluent with recycled hydraulic fracturing flowback fluids (HFFFs). We conducted a series of laboratory experiments in which the chemistry and NORM of different mix proportions of AMD and HFFF were examined after reacting for 48 h. The experimental data combined with geochemical modeling and X-ray diffraction analysis suggest that several ions, including sulfate, iron, barium, strontium, and a large portion of radium (60-100%), precipitated into newly formed solids composed mainly of Sr barite within the first ∼ 10 h of mixing. The results imply that blending AMD and HFFF could be an effective management practice for both remediation of the high NORM in the Marcellus HFFF wastewater and beneficial utilization of AMD that is currently contaminating waterways in northeastern U.S.A.

Hydraulic redistribution of soil water by roots affects whole-stand evapotranspiration and net ecosystem carbon exchange.

*Hydraulic redistribution (HR) of water via roots from moist to drier portions of the soil occurs in many ecosystems, potentially influencing both water use and carbon assimilation. *By measuring soil water content, sap flow and eddy covariance, we investigated the temporal variability of HR in a loblolly pine (Pinus taeda) plantation during months of normal and below-normal precipitation, and examined its effects on tree transpiration, ecosystem water use and carbon exchange. *The occurrence of HR was explained by courses of reverse flow through roots. As the drought progressed, HR maintained soil moisture above 0.15 cm(3) cm(-3) and increased transpiration by 30-50%. HR accounted for 15-25% of measured total site water depletion seasonally, peaking at 1.05 mm d(-1). The understory species depended on water redistributed by the deep-rooted overstory pine trees for their early summer water supply. Modeling carbon flux showed that in the absence of HR, gross ecosystem productivity and net ecosystem exchange could be reduced by 750 and 400 g C m(-2) yr(-1), respectively. *Hydraulic redistribution mitigated the effects of soil drying on understory and stand evapotranspiration and had important implications for net primary productivity by maintaining this whole ecosystem as a carbon sink.

Aspects of the in vitro bioactivity of hydraulic calcium (alumino)silicate cement

In response to a burgeoning interest in the prospective clinical applications of hydraulic calcium (alumino)silicate cements, the in vitro bioactivity and dissolution characteristics of a white Portland cement have been investigated. The formation of an apatite layer within 6 h of contact with simulated body fluid was attributed to the rapid dissolution of calcium hydroxide from the cement matrix and to the abundance of pre-existing Si-OH nucleation sites presented by the calcium silicate hydrate phase. A simple kinetic model has been used to describe the rate of apatite formation and an apparent pseudo-second-order rate constant for the removal of HPO42- ions frorn solultion has been calculated (k(2) = 5.8 x 10(-4) g mg(-1)). Aspects of the chemistry of hydraulic cements are also discussed with respect to their potential use in the remedial treatment of living tissue. (C) 2008 Wiley Periodicals, Inc. J Biomed Mater Res 90A: 166-174, 2009