Predicting dissolution patterns in variable aperture fractures: 1. Development and evaluation of an enhanced depth-averaged computational model

Water-rock interactions within variable-aperture fractures can lead to dissolution of fracture surfaces and local alteration of fracture apertures, potentially transforming the transport properties of the fracture over time. Because fractures often provide dominant pathways for subsurface flow and transport, developing models that effectively quantify the role of dissolution on changing transport properties over a range of scales is critical to understanding potential impacts of natural and anthropogenic processes. Dissolution of fracture surfaces is controlled by surface-reaction kinetics and transport of reactants and products to and from the fracture surfaces. We present development and evaluation of a depth-averaged model of fracture flow and reactive transport that explicitly calculates local dissolution-induced alterations in fracture apertures. The model incorporates an effective mass transfer relationship that implicitly represents the transition from reaction-limited dissolution to transport-limited dissolution. We evaluate the model through direct comparison to previously reported physical experiments in transparent analog fractures fabricated by mating an inert, transparent rough surface with a smooth single crystal of potassium dihydrogen phosphate (KDP), which allowed direct measurement of fracture aperture during dissolution experiments using well-established light transmission techniques [Detwiler, et al., 2003]. Comparison of experiments and simulations at different flow rates demonstrate the relative impact of the dimensionless Peclet and Damkohler numbers on fracture dissolution and the ability of the computational model to simulate dissolution. Despite some discrepancies in the small-scale details of dissolution patterns, the simulations predict the evolution of large-scale features quite well for the different experimental conditions. This suggests that our depth-averaged approach to simulating fracture dissolution provides a useful approach for extending laboratory results that are often limited in scale to scales that are more representative of geologic processes of interest

Understanding the interplay of increased pressure, residual heat and geologic and geochemical environment on radionuclide migration from underground nuclear tests

Understanding the processes that govern migration of radionuclides from underground nuclear tests is a problem of interest in many places in the world. Whether to gain understanding in the risks to human health, to guide water resources management decisions or as a surrogate for nuclear repository safety, there are a number of complex, coupled processes governing fate and form of radionuclides in the subsurface that need to be understood. These processes include the increased pore water pressure that is present after some underground nuclear tests, the residual heat signature from the nuclear detonation and the natural (i.e. pre-test) and altered (i.e. post-test) hydrogeologic and geochemical setting. We will focus on tests in saturated and unsaturated environments in Yucca Flat at the US Nevada Test Site. These tests all encompass a range of the aforementioned process interactions and provide insight into the broader context of understanding radionuclide migration from underground nuclear tests over a range of scale and environment

The MGDO software library for data analysis in Ge neutrinoless double-beta decay experiments

The GERDA and Majorana experiments will search for neutrinoless double-beta decay of germanium-76 using isotopically enriched high-purity germanium detectors. Although the experiments differ in conceptual design, they share many aspects in common, and in particular will employ similar data analysis techniques. The collaborations are jointly developing a C++ software library, MGDO, which contains a set of data objects and interfaces to encapsulate, store and manage physical quantities of interest, such as waveforms and high-purity germanium detector geometries. These data objects define a common format for persistent data, whether it is generated by Monte Carlo simulations or an experimental apparatus, to reduce code duplication and to ease the exchange of information between detector systems. MGDO also includes general-purpose analysis tools that can be used for the processing of measured or simulated digital signals. The MGDO design is based on the Object-Oriented programming paradigm and is very flexible, allowing for easy extension and customization of the components. The tools provided by the MGDO libraries are used by both GERDA and Majorana.Comment: 4 pages, 1 figure, proceedings for TAUP201

Fracture Permeability Evolution in Desert Peak Quartz Monzonite

Fracture flow experiments are being conducted on quartz monzonite core from the Desert Peak East EGS site, Churchill County, Nevada. The flow experiments are conducted at temperatures of 167-169 C and 5.5 MPa confining pressure through artificial fractures. Two injection fluids, a saline solution and a silica-bearing solution, have been used to date. Flow rates are typically 0.02 mL/min, but other rates have been used. The fracture surfaces are characterized with a contact profilometer. The profilometry data demonstrate that it is possible to fabricate statistically similar fracture surfaces and enable us to map aperture variations, which we use in numerical simulations. Effluent samples are collected for chemical analysis. The fluid pressure gradient is measured across the specimen and effective hydraulic apertures are calculated. The experiments show a reduction in permeability over time for both injection fluids, but a more rapid loss of permeability was observed for the silica-bearing solution. The calculated hydraulic aperture is observed to decrease by 17% for the saline solution and 75% for the silica-bearing fluid, respectively. Electrical resistivity measurements, which are sensitive to the ionic content of the pore fluid, provide additional evidence of fluid-rock interactions

Ground-State of Charged Bosons Confined in a Harmonic Trap

We study a system composed of N identical charged bosons confined in a harmonic trap. Upper and lower energy bounds are given. It is shown in the large N limit that the ground-state energy is determined within an accuracy of $\pm 8$ and that the mean field theory provides a reasonable result with relative error of less than 16% for the binding energy .Comment: 15 page

A basis for variational calculations in d dimensions

In this paper we derive expressions for matrix elements (\phi_i,H\phi_j) for the Hamiltonian H=-\Delta+\sum_q a(q)r^q in d > 1 dimensions. The basis functions in each angular momentum subspace are of the form phi_i(r)=r^{i+1+(t-d)/2}e^{-r^p/2}, i >= 0, p > 0, t > 0. The matrix elements are given in terms of the Gamma function for all d. The significance of the parameters t and p and scale s are discussed. Applications to a variety of potentials are presented, including potentials with singular repulsive terms of the form b/r^a, a,b > 0, perturbed Coulomb potentials -D/r + B r + Ar^2, and potentials with weak repulsive terms, such as -g r^2 + r^4, g > 0.Comment: 22 page

A Monte Carlo simulation of the Sudbury Neutrino Observatory proportional counters

The third phase of the Sudbury Neutrino Observatory (SNO) experiment added an array of 3He proportional counters to the detector. The purpose of this Neutral Current Detection (NCD) array was to observe neutrons resulting from neutral-current solar neutrino-deuteron interactions. We have developed a detailed simulation of the current pulses from the NCD array proportional counters, from the primary neutron capture on 3He through the NCD array signal-processing electronics. This NCD array Monte Carlo simulation was used to model the alpha-decay background in SNO's third-phase 8B solar-neutrino measurement.Comment: 38 pages; submitted to the New Journal of Physic

Investigation of \u3csup\u3e186\u3c/sup\u3eRe via radiative thermal-neutron capture on \u3csup\u3e185\u3c/sup\u3eRe

Partial -ray production cross sections and the total radiative thermal-neutron capture cross section for the 185Re(n,)186Re reaction were measured using the Prompt Gamma Activation Analysis facility at the Budapest Research Reactor with an enriched 185Re target. The 186Re cross sections were standardized using well-known 35Cl(n,)36Cl cross sections from irradiation of a stoichiometric natReCl3 target. The resulting cross sections for transitions feeding the 186Re ground state from low-lying levels below a cutoff energy of Ec=746keV were combined with a modeled probability of ground-state feeding from levels above Ec to arrive at a total cross section of σ0=111(6)b for radiative thermal-neutron capture on 185Re. A comparison of modeled discrete-level populations with measured transition intensities led to proposed revisions for seven tentative spin-parity assignments in the adopted level scheme for 186Re. Additionally, 102 primary rays were measured, including 50 previously unknown. A neutron-separation energy of Sn=6179.59(5)keV was determined from a global least-squares fit of the measured -ray energies to the known 186Re decay scheme. The total capture cross section and separation energy results are comparable to earlier measurements of these values

First Measurement of Coherent Elastic Neutrino-Nucleus Scattering on Argon

We report the first measurement of coherent elastic neutrino-nucleus scattering (\cevns) on argon using a liquid argon detector at the Oak Ridge National Laboratory Spallation Neutron Source. Two independent analyses prefer \cevns over the background-only null hypothesis with greater than $3\sigma$ significance. The measured cross section, averaged over the incident neutrino flux, is (2.2 $\pm$ 0.7) $\times$10$^{-39}$ cm$^2$ -- consistent with the standard model prediction. The neutron-number dependence of this result, together with that from our previous measurement on CsI, confirms the existence of the \cevns process and provides improved constraints on non-standard neutrino interactions.Comment: 8 pages, 5 figures with 2 pages, 6 figures supplementary material V3: fixes to figs 3,4 V4: fix typo in table 1, V5: replaced missing appendix, V6: fix Eq 1, new fig 3, V7 final version, updated with final revision

PT-symetrically regularized Eckart,Poeschl-Teller and Hulthen potentials

Version 1: The well known Eckart's singular s-wave potential is PT-symmetrically regularized and continued to the whole real line. The new model remains exactly solvable and its bound states remain proportional to Jacobi polynomials. Its real and discrete spectrum exhibits several unusual features. Version 2: Parity times time-reversal symmetry of complex Hamiltonians with real spectra is usually interpreted as a weaker mathematical substitute for Hermiticity. Perhaps an equally important role is played by the related strengthened analyticity assumptions. In a constructive illustration we complexify a few potentials solvable only in s-wave. Then we continue their domain from semi-axis to the whole axis and get the new exactly solvable models. Their energies come out real as expected. The new one-dimensional spectra themselves differ quite significantly from their s-wave predecessors.Comment: Original 10-page letter ``PT-symmetrized exact solution of the singular Eckart oscillator" is extended to a full pape