Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “representative volume element”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5

Critical mineral source potential from oil & gas produced waters in the United States

The volume of produced water, a by-product of oil & gas operations and other energy processes, has been growing across the United States (U.S.) along with the need to manage or recycle this wastewater. Produced water contains many naturally occurring elements of varying concentrations, including critical minerals which are essential to the clean energy transition. However, the current understanding of critical mineral concentrations in produced water and the associated volumes across the U.S. is limited. This study has assessed available databases and literature to gain insight into the presence and concentration of five high priority critical minerals, namely cobalt, lithium, magnesium, manganese, and nickel. The U.S. Geological Survey's National Produced Waters Geochemical Database was the main data source used for determining average critical mineral concentrations in produced water from the major oil and gas reservoirs in the U.S. The volumes of produced water for these major reservoirs were coupled with these concentrations to provide insights into where critical minerals are likely to have high abundance and therefore more recovery options. The analysis indicated the highest recovery potential for lithium and magnesium from produced water in the Permian basin and the Marcellus shale region. However, these assessments should be considered conservative due to the limited availability of reliable concentration data. Finally, it is expected more critical mineral recovery options could emerge with comprehensive characterization data from more recent and representative sources of produced water.

54 ENVIRONMENTAL SCIENCES↗

Finite-element boundary-integral simulation of thin wires and inhomogeneous penetrable bodies in subsurface multilayered anisotropic media

With the prevailing presence of drilling wells near the subsurface in mature oil and gas fields, the application of electromagnetic methods can be particularly challenging where the electromagnetic field is affected by the steel casing. In the past decades, borehole-to-surface and crosswell electromagnetic methods have been utilized for monitoring of reservoir and underground CO 2 storage. This paper presents a unified finite-element boundary-integral (FEBI) method capable of simultaneously modeling the complex electromagnetic interactions between thin metallic wires (representing steel casings) with 3D trajectory and arbitrary 3D inhomogeneous penetrable bodies (such as CO 2 plumes or hydrocarbon reservoirs) within anisotropic multilayered subsurface environments. Unlike existing approaches that treat these components separately or require dense discretization, or are limited to vertical wells, our unified formulation preserves flexible electromagnetic coupling while delivering improved computational efficiency. Assuming the background formation is multilayered anisotropic media, the surface integral equation method is applied to model the thin wires and boundaries of the inhomogeneous bodies. Meanwhile, the finite element method is applied to model the volume of inhomogeneous bodies. Here, the performance of the proposed FEBI method is assessed through comparison with reference numerical results and its practical significance is demonstrated through CO 2 plume monitoring scenarios.

97 MATHEMATICS AND COMPUTING↗

Comparison of Fission Product Release Predictions using PARFUME and BISON with Results from the AGR-3/4 Irradiation Experiment

The PARFUME (PARticle Fuel ModEl) fuel performance modeling code and the BISON nuclear fuel performance application built on the Multiphysics Object-Oriented Simulation Environment (MOOSE) finite element library were used to predict the fission product release from tristructural isotropic (TRISO) coated fuel particles and compacts during the third and fourth irradiation experiment of the Advanced Gas Reactor (AGR-3/4) Fuel Development and Qualification Program. The fuel performance modeling codes PARFUME and BISON modeled the AGR-3/4 irradiation experiment using the fuel compact time-averaged volume averaged (TAVA) daily temperatures for a total irradiation duration of 369.1 effective full power days (EFPD) to predict the release fraction of the fission product silver (Ag-110m) from a representative TRISO-coated fuel particle from AGR-3/4 compacts. Post-irradiation examination (PIE) measurements provided data on the release of these fission products in the compacts outside of the silicon carbide (SIC) layer. The PARFUME and BISON results were then compared to the silver release measured from compact gamma scanning. The results showed good agreement between PARFUME and BISON but both codes under-predicted the silver release fraction for all the compacts. In addition, BISON was used to model and predict the fission product concentration radial profile outside of the compacts in capsules’ inner and outer rings. These rings were either comprised of matrix and/or structural graphite. To obtain the concentration profiles of silver, cesium, and strontium, a sorption isotherm model was developed in BISON to capture the effects of fission product transport across the gaps between the concentric rings. The general shape of the concentration radial profiles as calculated by BISON were similar in the inner ring (IR) but varied in the outer ring (OR) depending on the fission product of interest or capsule temperature. Using this methodology and model, BISON now has the capability to aid in developing new fission product diffusion coefficients for matrix or structural graphite materials.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Level-set topology optimization with PDE generated conformal meshes

This paper presents a level-set topology optimization approach that uses conformal meshes for the analysis of the displacement field. The structure’s boundary is represented by the iso-contour of a level-set field discretized on a fixed background design mesh. The conformal mesh is updated for each design iteration via a PDE based mesh morphing process that identifies the set of facets in the background mesh that are homeomorphic to the boundary and relaxes the homeomorphic mesh to conform to the structure’s boundary and ensure high element quality. The conformal mesh allows for a more accurate computation of the response versus density and some level-set based methods which interpolate material properties using the volume fraction. Numerical examples illustrate the proposed approach by optimizing linear-elastic two- and three-dimensional structures, wherein insight into the performance of the mesh morphing process is provided. The examples also highlight the scalability of the approach.

42 ENGINEERING↗

Microstructure Generation via Generative Adversarial Network for Heterogeneous, Topologically Complex 3D Materials

Using a large-scale, experimentally captured 3D microstructure dataset, we implement the generative adversarial network (GAN) framework to learn and generate 3D microstructures of solid oxide fuel cell electrodes. The generated microstructures are visually, statistically, and topologically realistic, with distributions of microstructural parameters, including volume fraction, particle size, surface area, tortuosity, and triple phase boundary density, being highly similar to those of the original microstructure. These results are compared and contrasted with those from an established, grain-based generation algorithm (DREAM.3D). Importantly, simulations of electrochemical performance, using a locally resolved finite element model, demonstrate that the GAN generated microstructures closely match the performance distribution of the original, while DREAM.3D leads to significant differences. Finally, the ability of the generative machine learning model to recreate microstructures with high fidelity suggests that the essence of complex microstructures may be captured and represented in a compact and manipulatable form.

36 MATERIALS SCIENCE↗

Science Capsule: Towards Sharing and Reproducibility of Scientific Workflows

Workflows are increasingly processing large volumes of data from scientific instruments, experiments and sensors. These workflows often consist of complex data processing and analysis steps that might include a diverse ecosystem of tools and also often involve human-in-the-loop steps. Sharing and reproducing these workflows with collaborators and the larger community is critical but hard to do without the entire context of the workflow including user notes and execution environment. In this paper, we describe Science Capsule, which is a framework to capture, share, and reproduce scientific workflows. Science Capsule captures, manages and represents both computational and human elements of a workflow. It automatically captures and processes events associated with the execution and data life cycle of workflows, and lets users add other types and forms of scientific artifacts. Science Capsule also allows users to create `workflow snapshots' that keep track of the different versions of a workflow and their lineage, allowing scientists to incrementally share and extend workflows between users. Our results show that Science Capsule is capable of processing and organizing events in near real-time for high-throughput experimental and data analysis workflows without incurring any significant performance overheads.

Ghoshal, Devarshi↗

Permanent Magnets Featuring Heavy Main Group Elements for Magnetic Anisotropy

Permanent magnets are the functional component of electric motors and generators found in numerous renewable energy applications. To improve energy conversion in such applications, we require fundamentally new magnets that generate higher magnetic flux per volume while retaining the properties conferred by rare-earth elements incorporated into current technologies. We hypothesize that by engendering a covalent interaction between two elements, we can access a new regime of magnetic materials where the two components of a magnetic moment—spin and orbital angular momentum—come from two separate atoms to form a complete magnetic moment. Our previous research utilized high-pressure conditions to discover two new candidate materials ideal for assessing this hypothesis. The first material, FeBi 2 , enables the study of an unprecedented solid state metal-metal bonding interaction. The second, MnBi 2 represents the second member of the promising Mn–Bi family known for its magnetic properties. Importantly, MnBi 2 is isostructural to FeBi 2 . Together these chemically simple but magnetically rich materials provide an elegant platform for elucidating fundamental design principals of magnetic anisotropy while inspiring the synthesis of new magnetic materials. More generally, solid-state chemistry remains a synthetic black box. To fully harness the potential of such new materials, it is vital to create a window into both the structure and the properties of our high-pressure materials. One missing area of importance to energy science is high-pressure magnetometry. To this end, we assessed the magnetic structure of these materials through two objectives: (1) we performed magnetometry measurements at high-pressures, including performing high-pressure X-ray Magnetic Circular Dichroism (XMCD) experiments on these materials that are only synthesized at high-pressure and cannot be recovered to ambient conditions and (2) we targeted the recovery of these materials using dynamic compression approaches, starting with laser shock techniques. Future directions include expanding our focus to ternary phase space, employing high-throughput computational approaches to direct our search for magnetic materials across the Periodic Table. This award was originated at Northwestern, and reopened as a new grant at MIT. This is the close out for Northwestern.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Critical Design Elements for an On-Site Mercury Concentration and Holding Facility - 20175

Administrative changes and funding issues have delayed the Department of Energy (DOE) from establishing and commissioning a National Repository for mercury as required under the Mercury Export Ban Act (MEBA) of 2008. In the interim, a few hazardous waste Treatment, Storage, and Disposal Facilities (TSDFs) have been authorized to accept mercury classified as a hazardous waste as under Subtitle C of the Resource Conservation and Recovery Act (RCRA). Waste generators are charged (either by weight or by volume) for the waste they manifest to these TSDFs. In many cases, mercury and mercury compounds have chemically sorbed to (or physically lodged into) the interstitial spaces of a substrate material. This almost always causes the amount of waste generated to be unnecessarily large, due to the extra volume (or weight) of the accompanying, non-mercury materials. Depending on a multitude of factors, it could be financially beneficial for waste generators to take greater control of their mercury affairs by separating these substrate materials from the mercury and mercury compounds that have become associated with them. Building and operating an on-site mercury Concentration and Holding Facility (CHF) could help relieve a waste generator's overall cost burden. The purpose of such CHFs is not to treat or dispose of mercury, but to recover and concentrate mercury from bulky mercury-laden materials to the extent that the host substrate to which the mercury may have previously been sorbed could pass a Toxicity Characteristic Leaching Procedure (TCLP) test and either be disposed of at a conventional landfill, or be disposed of at a hazardous waste landfill, but absent additional restrictions specifically for mercury. A CHF represents a way for waste generators to save money by separating mercury from common host materials, concentrating it, and then safely storing it until transportation to an authorized TSDF takes place This paper focuses on the critical design elements waste generators should consider when designing and operating a mercury CHF. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Study of Additive Manufacturing Application to Geothermal Technologies

Geothermal reservoir characterization, field construction, and reservoir operations are very technology intensive activities that contribute significantly to the cost of delivering electricity produced from geothermal resources. Many geothermal technologies, such as downhole tools and drilling equipment, have unusual material, design, and manufacturing considerations dictated by the harsh geothermal environment and extreme aspect ratios required for deployment in a borehole. An additional challenge that faces geothermal applications is the low tool production volume needed to support the industry. Whereas tens of thousands of Oil & Gas wells are drilled and completed in the U.S. annually, there are typically only tens of geothermal wells that are drilled and completed. If a tool typically used in Oil & Gas applications cannot be directly used for geothermal, then the cost associated with making the tool suitable for geothermal is often prohibitive. There is therefore a much smaller inventory of technologies available to geothermal as compared to Oil & Gas and the level of efficiency and sophistication associated with field practice suffers accordingly. A number of advanced manufacturing methods, such as additive manufacturing, have received increased R&D as well as commercial attention in recent years because of their ability to rapidly prototype complex parts. Additive manufacturing in particular provides an opportunity to increase the technology available to the geothermal industry by either reducing fabrication costs associated with complex components or enabling economic production of low volume parts where specialized tooling is often required. Additional potential benefits of additive manufacturing include increased design freedom to make higher performing parts that cannot be made conventionally, the ability to integrate components into assemblies without joining operations, and the ability to economically fabricate variations on design in cases, such as casting molds, where there are large up-front costs associated with tooling. We have recently completed a study that investigated technology needs, representative use cases, manufacturability, and a techno-economic framework for comparing conventional to additive manufacturing methods for geothermal applications. This paper will provide an overview of this recent effort, describe the different elements of the assessment, and summarize the key takeaways related to both the feasibility of using additive manufacturing for geothermal technology applications as well as the potential benefits and impacts.

Polsky, Yarom↗

Impacts of Lanthanide Ion Complexation on the Radiation Robustness of Diglycolamide Extractants

Nuclear power reactors represent a clean and reliable source of baseload energy. However, nuclear reactors produce used nuclear fuel (UNF), which, if not recycled/reprocessed, must be disposed of as high-level radioactive waste, necessitating long-term storage options. UNF contains approximately a third of the periodic table, including the radioactive minor actinides (MA) americium and curium. Removal of these MA from UNF would greatly reduce the radiotoxic burden, volume, and cost of a storage facility or geological repository. Thus, a significant global effort has been devoted to the design of reprocessing strategies for the clean separation of the MA from UNF. As MA separations are difficult, due to the presence of lanthanide elements in UNF and their remarkably similar chemical properties, many extractants have been studied. The diglycolamide (DGA) class of extractants are promising for this separation as they exhibit high distribution coefficients for the MAs, high metal loading capacities, and are robust in highly acidic radiation environments.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Characterization of Incremental Markings in the Sagittal Otolith of the Pacific Sardine ( Sardinops sagax ) Using Different Imaging Modalities

Teleost fish possess calcium carbonate otoliths located in separate chambers (utriculus, sacculus, and lagena) of their membranous labyrinth. This study analyzed the surface topography of the sagittal otolith of the Pacific sardine (Sardinops sagax) and the daily and annual increments in these otoliths. The otolith surface, characterized by laser scanning confocal microscopy for the first time, consisted of a system of prominent ridges and valleys (grooves), but it is unclear whether these structures are functional or represent time-resolving markings reflecting growth periodicity. Within the first-year volume, daily increments, each consisting of an incremental (more mineralized) and a discontinuous (less mineralized) zone, were resolved by optical microscopy and backscattered electron (BSE) imaging in the scanning electron microscope (SEM). Daily growth increments could, however, not be resolved in volumes formed after the first year, presumably because otolith growth markedly slows down and spacing of incremental markings narrows in older fish. Throughout otolith growth, the crystalline network continues across the discontinuous zones. Fluorochrome labeling provides additional information on growth after the first year. Compared with optical and BSE imaging, synchrotron microComputed Tomography of intact otoliths (with 0.69 µm volume elements) was less able to resolve daily increments; X-ray phase contrast reconstructions provided more detail than reconstructions with absorption contrast. Future research directions are proposed.

59 BASIC BIOLOGICAL SCIENCES↗

Experimental and Computational Studies on Superhard Material Rhenium Diboride under Ultrahigh Pressures

An emerging class of superhard materials for extreme environment applications are compounds formed by heavy transition metals with light elements. In this work, ultrahigh pressure experiments on transition metal rhenium diboride (ReB 2 ) were carried out in a diamond anvil cell under isothermal and non-hydrostatic compression. Two independent high-pressure experiments were carried out on ReB 2 for the first time up to a pressure of 241 GPa (volume compression V/V 0 = 0.731 ± 0.004), with platinum as an internal pressure standard in X-ray diffraction studies. The hexagonal phase of ReB 2 was stable under highest pressure, and the anisotropy between the a-axis and c-axis compression increases with pressure to 241 GPa. The measured equation of state (EOS) above the yield stress of ReB 2 is well represented by the bulk modulus K 0 = 364 GPa and its first pressure derivative K 0 ´ = 3.53. Corresponding density-functional-theory (DFT) simulations of the EOS and elastic constants agreed well with the experimental data. DFT results indicated that ReB 2 becomes more ductile with enhanced tendency towards metallic bonding under compression. The DFT results also showed strong crystal anisotropy up to the maximum pressure under study. The pressure-enhanced electron density distribution along the Re and B bond direction renders the material highly incompressible along the c-axis. Our study helps to establish the fundamental basis for anisotropic compression of ReB 2 under ultrahigh pressures

36 MATERIALS SCIENCE↗

Influence of cementation on the yield surface of rocks numerically determined from digital microstructures

Digital Rock Physics has reached a level of maturity on the characterisation of primary properties that depend on the microstructure – such as porosity, permeability or elastic moduli – by numerically solving field equations on μCT scan images of rock. After small deformations or at depth though, most rocks eventually reach their limit of elasticity and the complementary plastic properties are needed to describe the full mechanical behaviour. Currently, determination of a rock’s yield surface from its microstructure is often restricted to semi-analytical criteria derived by limit analysis or numerical simulations performed on idealised geometries. Such simplification lacks representativeness, particularly for processes that affect directly the pore-grain interface such as the cementation phenomenon, happening during diagenesis. Eventually, only direct numerical simulation of elasto-plasticity performed on digitalised microstructures can be used to assess the strength of different cemented materials and its evolution with the alteration of the microstructure. In this study, we provide a comprehensive parametric study on the impact of cementation on rock strength for real microstructures of cemented granular materials. Compared to most previous studies, the whole yield surface is determined numerically (using Finite Element Method) in order to assess the influence of cementation for different stress-paths. The previously known tendency of rock to strengthen with increasing cementation volume is verified. New results on the influence of cement property namely Young’s modulus, friction and cohesion on the rock’s yield surface are explored. The envelopes obtained are compared to the ones obtained by experimental data and existing models. Finally, the framework presented in this study showcases the wider possibility of determining any rock’s or porous material’s yield surface from its microstructure.

36 MATERIALS SCIENCE↗

Simultaneous bright- and dark-field X-ray microscopy at X-ray free electron lasers

Abstract The structures, strain fields, and defect distributions in solid materials underlie the mechanical and physical properties across numerous applications. Many modern microstructural microscopy tools characterize crystal grains, domains and defects required to map lattice distortions or deformation, but are limited to studies of the (near) surface. Generally speaking, such tools cannot probe the structural dynamics in a way that is representative of bulk behavior. Synchrotron X-ray diffraction based imaging has long mapped the deeply embedded structural elements, and with enhanced resolution, dark field X-ray microscopy (DFXM) can now map those features with the requisite nm-resolution. However, these techniques still suffer from the required integration times due to limitations from the source and optics. This work extends DFXM to X-ray free electron lasers, showing how the $$10^{12}$$ 10 12 photons per pulse available at these sources offer structural characterization down to 100 fs resolution (orders of magnitude faster than current synchrotron images). We introduce the XFEL DFXM setup with simultaneous bright field microscopy to probe density changes within the same volume. This work presents a comprehensive guide to the multi-modal ultrafast high-resolution X-ray microscope that we constructed and tested at two XFELs, and shows initial data demonstrating two timing strategies to study associated reversible or irreversible lattice dynamics.

47 OTHER INSTRUMENTATION↗

Thermomechanical Microstructural Predictions of Fracture Nucleation of Zircaloy-4 Alloys With δ and ε Hydride Distributions

Abstract A crystalline dislocation-density formulation that was incorporated with a nonlinear finite-element (FE) method was utilized to understand and to predict the thermomechanical behavior of an hexagonal closest packed (h.c.p.) zircaloy system with hydrides with either face-centered cubic (f.c.c.) or body-centered cubic (b.c.c.) hydrides. This formulation was then used with a recently developed fracture methodology that is adapted for finite inelastic strains and multiphase crystalline systems to understand how different microstructurally based fracture modes nucleate and propagate. The interrelated microstructural characteristics of the different crystalline hydride and matrix phases with the necessary orientation relationships (ORs) have been represented, such that a detailed physical understanding of fracture nucleation and propagation can be predicted for the simultaneous thermomechanical failure modes of hydride populations and the matrix. The effects of volume fraction, morphology, crystalline structure, and orientation and distribution of the hydrides on simultaneous and multiple fracture modes were investigated for radial, circumferential, and mixed distributions. Another key aspect was accounting for temperatures changes due to the effects of thermal conduction and dissipated plastic work and their collective effects on fracture. For hydrided aggregates subjected to high temperatures, thermal softening resulted in higher ductility due to increased dislocation-density activity, which led to higher shear strain accumulation and inhibited crack nucleation and growth. The predictions provide validated insights into why circumferential hydrides are more fracture-resistant than radial hydrides for different volume fractions and thermomechanical loading conditions.

Engineering↗

In situ measurement of three-dimensional intergranular stress localizations and grain yielding under elastoplastic axial-torsional loading

The three-dimensional grain-averaged response of solid bar samples under non-proportional (NP) elastoplastic axial-torsional loading was investigated using in situ high energy diffraction microscopy (HEDM) and companion crystal plasticity finite element (CPFE) modeling. Important stress metrics including applied shear (σ θZ ) and axial (σ ZZ ) stress tensor components, stress and stress deviator tensor invariants (I 1 , J 2 , and J 3 ), von Mises equivalent stress (σ$^{grain}_{VM}$), maximum resolved shear stress (mRSS), stress triaxiality (η), and lode angle parameter ($\barθ$) values were tracked for ~300 grains under two different loading conditions: (1) Torsion-dominated loading (low NP) and (2) Tension-torsion loading (high NP) in equiatomic NiCoCr, a representative multicomponent face-centered cubic (FCC) superalloy. Overall, significant stress localizations existed within both samples as evidenced by the radial dependence of grain-resolved σ θZ , σ$^{grain}_{VM}$, and J 2 ; by comparison, I 1 , J 3 , η, and $\barθ$ metrics did not show discernible trends within the volume. These stress localizations reveal a complex interplay between axial and shear stress components (e.g., stress coupling) resulting in grain yielding near the sample surface largely driven by shear stress, whereas internal grain yielding was largely accommodated by axial stress. Grain-resolved stress localization trends were described well by the CPFE model, although some discrepancies in magnitude occurred, particularly for volumetric stress metrics (I 1 and η) due to initial type II residual stress distributions. The superposition of initial residual stress states onto CPFE grain-resolved data significantly improved model accuracy for η. This suggests that residual stresses more strongly influence the simulation of volumetric rather than deviatoric (yield) stress metrics.

36 MATERIALS SCIENCE↗

Comparing and contrasting In-Vial and full-scale systems for sparging volatile analytes

Here, in-vial sparging was demonstrated as an effective, practical alternative to a full-scale sparging system for supporting the analysis of volatile constituents. Using elemental mercury (Hg 0 ) and toluene as representative purgeable analytes, the mass removal for various sparge configurations was measured and a reduced order model was developed and validated. In the primary experiments, Hg 0 in the sparge gas was trapped on activated carbon or gold, thermally desorbed, and quantified using atomic absorption or atomic fluorescence spectroscopy. Toluene experiments using the same in-vial sparge apparatus and sparge parameters were performed to demonstrate the applicability of the reduced order model to a broad range of compounds. Toluene removal was tracked by measuring the remaining toluene in sparged aliquots using Ultraviolet-visible (UV–Vis) spectroscopy. For the sparging, flow rates varied from 25 to 75 mL/min for periods from 0 to 30 min. Sparge performance, mass removal as a function of time, and sparge gas volume were measured for both in-vial and full-scale systems. A model based on dimensionless Henry's Law coefficient, normalized sparge gas volume, and fractional extent of equilibrium matched the experimental data for both compounds and provides a practical tool for future applications. For the conditions tested in this study, the calibrated model indicated that the sparge gas in the in-vial system reached approximately 33% of its equilibrium value before exiting the water surface, while a full-scale system reached approximately 100%. The tests validated the quality, reproducibility, and predictability of sparging performance for both full scale and in-vial sparge systems. Related factors such as waste generation, worker risk, and labor were also assessed. Full scale sparge systems provide the advantage of lower detection levels due to larger sample volume, while the in-vial sparge systems provide advantages for most other factors; including automatability, reducing secondary wastes, lessening the need to clean and check the sparge apparatus, and lowering labor and costs. The data and associated reduced order model support continued development and deployment of in-vial sparge platforms as a practical option for analysis of purgeable analytes such as volatile organic compounds and volatile metals/organometallics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Scaling hyporheic nitrogen cycling in large river alluvial aquifers

Research Element 1: Denitrification is a difficult flux to measure when fluxes are low and when the ecosystem is difficult to sample. These attributes describe Nyack floodplain, a large gravel-bedded aquifer adjacent to the Middle Fork Flathead River. Given that denitrification is simply respiration with nitrate as a terminal electron acceptor, we reasoned that N 2 production should negatively covary with oxygen respiration so long that oxygen is available. We sampled 16 wells 8 times throughout a year and analyzed these water samples on a membrane-inlet mass spectrometer for N 2 :Ar and O 2 . We fit multilevel models to N 2 :Ar as a function of dissolved oxygen. Within wells, there was no relationship between O 2 and N 2 :Ar. Among wells N 2 :Ar increased as O 2 declined. Based on the slope of this relationship, ~ -0.07, and estimated oxygen respiration in Nyack derived from prior studies, a possible range for denitrification is 0-8 mmol m -2 d -1 corresponding to a turnover rate for nitrate of 0-0.2/d. This estimated turnover rate is in the lower 20% of rates measured in the LINX 2 study, yet shows that even in a low nitrate (0.005 mmol/L) aquifer denitrification can be a substantial sink for dissolved N. Research Element 3 and 4: We developed a new theoretical approach to describe the relationship between water inflow, water outflow, and water age in porous media control volumes under dynamic flow conditions. Our approach had three distinct steps. First we described the application of well-established equations borrowed from Chemical Engineering to describe water age distributions in the hyporheic zone under steady state conditions. We published the associated results in the scientific journal PLOS ONE. Second, we developed a new theory ("flow spans") to describe how predicted steady-state water age distributions will respond to time varying imbalances in hyporheic discharge and recharge. Finally, we developed a proof-of-concept simulation model based on the new theory and used the model to examine the effects of flow regulation on a simulated river with characteristics representative of western U.S. montane rivers with expansive hyporheic zones.

54 ENVIRONMENTAL SCIENCES↗