Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “solution chemistry”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Engineering In Situ Catalytic Cleaning Membrane Via Prebiotic-Chemistry-Inspired Mineralization

Pressure-driven membrane separation promises a sustainable energy-water nexus but is hindered by ubiquitous fouling. Natural systems evolved from prebiotic chemistry offer a glimpse of creative solutions. Herein, a prebiotic-chemistry-inspired aminomalononitrile (AMN)/Mn 2+ -mediated mineralization method is reported for universally engineering a superhydrophilic hierarchical MnO 2 nanocoating to endow hydrophobic polymeric membranes with exceptional catalytic cleaning ability. Green hydrogen peroxide catalytically triggered in-situ cleaning of the mineralized membrane and enabled operando flux recovery to reach 99.8%. The mineralized membrane exhibited a 9-fold higher recovery compared to the unmineralized membrane, which is attributed to active catalytic antifouling coupled with passive hydration antifouling. Electron density differences derived from the precursor interaction during mediated mineralization unveiled an electron-rich bell-like structure with an inner electron-deficient Mn core. This work paves the way to construct multifunctional engineered materials for energy-efficient water treatment as well as for diverse promising applications in catalysis, solar steam generation, biomedicine, and beyond.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Adsorption and exchange reactions of iodine molecules at the alumina surface: modelling alumina-iodine reaction mechanisms

Harnessing aluminum oxidation energy requires navigating the particle's passivation shell composed of alumina. The shell is a barrier to aluminum oxidation but can also exothermically react with halogenated species and therefore contribute to the overall energy generated during aluminum particle combustion. Fluorination reactions with alumina have been studied because fluorine is abundant in binder formulations that commonly surround aluminum particles in an energetic mixture. However, iodine has emerged as an alternative halogenated-based binder or oxidizer because iodine gas provides ancillary benefits such as chemical neutralization of biological agents or sterilization of contaminated environments. Furthermore, this study used density functional theory (DFT) calculations to evaluate potential reaction pathways for aluminum–iodine combustion. Relative to fluorinated fragments such as HF and F – , the adsorption energies associated with HI and I– are nearly triple the exchange reaction energy available from fluorination reactions with alumina (–189 and –278 kJ mol –1 for HI and I – , respectively). However, exchange reactions between iodinated species and the alumina surface are energetically unfavorable. These results explain that through adsorption, alumina surface exothermic reactions with iodine are more energetic than with fluorine fragments. Experiments performed with differential scanning calorimetry (DSC) confirm the higher magnitude of energy generated for iodination compared with fluorination reactions with alumina. Additionally, strong adsorption energies can promote synthesis of new shell chemistries. Adsorption in solution will promote alumina dissolution and iodine precipitation reactions to produce hydroxyl complexes and iodinated species synthesized on the surface of the particle, thereby replacing alumina with alternative passivation shell chemistry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Relationships between Water’s Structure and Solute Affinity at Polypeptoid Brush Surfaces

Excellent antifouling surfaces are generally thought to create a tightly bound layer of water that resists solute adsorption, and highly hydrophilic surfaces such as those with zwitterionic functionalities are of significant current interest as antifoulant strategies. Furthermore, despite significant proofs-of-concept, we still lack a fundamental understanding of how the nanoscopic structure of this hydration layer translates to reduced fouling, how surface chemistry can be tuned to achieve antifouling through hydration water, and why, in particular, zwitterionic surfaces seem so promising. Here, we use molecular dynamics simulations and free energy calculations to investigate the molecular relationships among surface chemistry, hydration water structure, and surface–solute affinity across a variety of surface-decorated chemistries. Specifically, we consider polypeptoid-decorated surfaces that display well-known experimental antifouling capabilities and that can be synthesized sequence specifically, with precise backbone positioning of, e.g., charged groups. Through simulations, we calculate the affinities of a range of small solutes to polypeptoid brush surfaces of varied side-chain chemistries. We then demonstrate that measures of the structure of surface hydration water in response to a particular surface chemistry signal solute–surface affinity; specifically, we find that zwitterionic chemistries produce solute–surface repulsion through highly coordinated hydration water while suppressing tetrahedral structuring around the solute, in contrast to uncharged surfaces that show solute–surface affinity. Based on the relationship of this structural perturbation to the affinity of small-molecule solutes, we propose a molecular mechanism by which zwitterionic surface chemistries enhance solute repulsion, with broader implications for the design of antifouling surfaces.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Slow Strain Rate Testing of A537 Tank Wall Material

At Savannah River Site (SRS), High-Level Waste is stored in below-grade carbon steel tanks. This waste in part consists of sludge, salt cake, and/or supernate. Preparation of this waste for future processing involves dissolution of the salt cake layer. The salt dissolution process can create conditions that leave the carbon steel tanks susceptible to localized corrosion. The salt to be dissolved contains high concentrations of nitrate, that once released, create an environment that may be conducive to pitting corrosion and/or stress corrosion cracking (SCC) of carbon steel. The salt dissolution process also liberates interstitial liquid trapped between the salt crystals. This liquid is initially high in nitrite and hydroxide concentration. High pH and greater ratios of nitrite to nitrate act as inhibitors to minimize corrosion of carbon steel in high nitrate environments. However, as dissolution proceeds, the concentration of nitrate will increase, while the hydroxide and nitrite concentration of the interstitial liquid will deplete and become insufficient to prevent the onset of corrosion attack. Tank blending and the addition of inhibitors are used to ensure adequate concentrations of hydroxide and nitrite. However, this is not desirable during salt dissolution as it can reduce process efficiency and increase the amount of waste that needs processing. This testing program was designed to examine the risk of SCC associated with utilizing the pitting factor (PF) and nitrite/nitrate (NO 2 - /NO 3 - ) ratio limits for handling dissolved salt solutions at an elevated temperature in the carbon steel waste tanks. The previously identified limits are a PF of 1.2 and an NO 2 - /NO 3 - ratio of 0.15. The results indicate that as long as the NO 2 - /NO 3 - ratio exceeds 0.1 and the PF is above approximately 0.8, there is a discernible safety margin between the open circuit potential (OCP) and the critical cracking potential (CCP) observed during applied potential testing. However, this margin, defined by the difference between the OCP and CCP, is relatively narrow, ranging from 0.1 to 0.25 volts. This small margin raises concerns about potential shifts in OCP during waste retrieval operations, which could inadvertently increase the risk of SCC if the OCP approaches or exceeds the CCP. These results confirm that dissolved salt solutions provide a potent chemistry that, under certain conditions, makes carbon steel susceptible to SCC. The next question to consider is the influence these results have on decisions for storage and retrieval of waste from the tanks. For Type III/IIIA waste tanks, the risk of SCC remains very low. First, and most importantly, the post-weld stress relief of the tanks has reduced the residual stress near the welds. Thus, without the stress component, SCC risk is minimized. The material of construction (A537 Carbon steel) for the Type III/IIIA tanks is superior to the steel in its resistance to SCC than the steel that was utilized for the Type I, II, and IV tanks (A285 carbon steel). From a chemistry control standpoint for a Type III/IIIA tank directly involved with handling dissolved salt solutions, the PF and NO 2 - /NO 3 - ratio limits may be utilized wherein chemistry control provides an extra layer of defense against SCC. Chemistry control for a Type III/IIIA tank minimizes the risk for a tank that may receive the dissolved salt solution, particularly if that tank is a Type I, II, or IV waste tank. On the other hand, if the dissolved salt solution is handled by a Type I, II, or IV waste tank the risk of SCC is real. The potent chemistry, absence of stress relief, and inferior material result in a condition that is conducive to cracking. Efforts should be made to either avoid transferring waste that may not meet the PF and NO 2 - /NO 3 - ratio criteria to one of these tanks or if it is unavoidable, take measures to minimize the consequences of a leak. As shown by these tests, even if the PF and NO 2 - /NO 3 - ratio criteria are met, there is a risk that the tank potential may be disturbed in the positive direction and the risk of SCC increase.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Dissolution and reprecipitation of amorphous silica in silica Rich shales induces Non-Monotonic evolution of porosity in acidic reactive environments

Advances in sustainable subsurface energy technologies are crucial for meeting our energy and resource needs for a climate-resilient future. Novel strategies to harness subsurface shale reservoirs for recovering valuable metals and for enabling CO 2 storage are influenced by the morphological and mineralogical heterogeneities of these materials. In this context, delineating the interactions of highly acidic solutions such as wet supercritical CO 2 on shales with varying mineralogy is crucial to inform the stability of caprock seal for CO 2 storage and enhancements in permeability for fluid transport, reactivity, and storage. The feedback chemical effects associated with the interactions of acidic solutions on the morphologies and mineralogies of shales have not been extensively investigated. These insights are crucial for assessing temporal changes in the reactivity and the fate of the fluids in subsurface environments. In this study, we investigate the effect of 1M HCl solution on the chemistry and morphology of three different shale samples with varying carbonate, clay and silica contents. An increase in the amorphous content, from 37% to 41.3%, of silica-rich and carbonate/clay lean shale is noted due to reactions with an acidic solution which is attributed to the dissolution of Si-bearing phases such as clays, accompanied by SiO 2 precipitation. In shales bearing high content of clays and carbonates, significant increase in the pore volumes and surface areas are noted. Non-monotonic changes in the micron-scale porosity of silica rich–carbonate/clay lean (e.g., Mowry shale) are noted using in-situ X-ray microtomography experiments. Due to the initial mobilization of silica and dissolution of carbonate/clay phases, the total porosity slightly increases from 6.7% to 10.7% followed by a decrease to similar to 4% caused by SiO 2 reprecipitation. In conclusion, these findings suggest that even though silica is less reactive in acidic environments, the changes in the amorphous and crystalline content due to dissolution and reprecipitation alter the porosity and fluid flow paths.

04 OIL SHALES AND TAR SANDS↗
NASA NTRS1989

Three dimensional PNS solutions of hypersonic internal flows with equilibrium chemistry

An implicit procedure for solving parabolized Navier-Stokes equations under the assumption of a general equation of state for a gas in chemical equilibrium is given. A general and consistent approach for the evaluation of Jacobian matrices in the implicit operator avoids the use of unnecessary auxiliary quantities and approximations, and leads to a simple expression. Applications to two- and three-dimensional flow problems show efficiency in computer time and economy in storage.

Liou, May-Fun↗
NASA NTRS1973

Survey of application of radiation to preparative chemistry

The use of radiation for preparative chemistry in liquid solutions is investigated. General principles are presented and preparations involving reduction, oxidation, polymerization, and decomposition are given. The use of various solvents, water, other inorganic liquids and organic liquids for this purpose is discussed. Finally, a commentary is made on some specific applications where radiation chemistry as a preparative technique may be useful.

Philipp, W. H.↗
DOE OSTI2024

KCl-UCl 3 molten salts investigated by Ab Initio Molecular Dynamics (AIMD) simulations

Ab Initio Molecular Dynamics (AIMD) simulations are performed on molten KCl-UCl3 salt mixtures to determine energies, heat capacities, and densities. The density-dependent energy correction (DFT-dDsC), Grimme et al.'s DFT-D3, and Langreth & Lundqvist (vdW-cx) models are used for dispersion forces and combined with the Perdew-Burke-Ernzerhof (PBE) exchange-correlation potential with a Hubbard U parameter for the 5f electrons of uranium. After validating predictions for the end-member systems to literature data, KCl-UCl 3 mixtures are studied at select temperatures. Densities and energies both deviate from ideal solution behavior, with the maximum deviation occurring around 36% UCl 3 for mixing energies and slightly lower (29% UCl3) for densities. Compared to the NaCl-UCl 3 system, which was previously investigated using the same simulation methodologies, the KCl-UCl 3 density and mixing energy deviations from ideal solution behavior are larger by almost a factor of two. No deviation from ideal solution behavior for heat capacity was observed. The AIMD predictions for mixing energies and densities agree qualitatively with experimental data, though the spread in data obtained from the various dispersion force models utilized, measurements, and empirical estimates makes strong conclusions difficult. The dependence of thermodynamic and thermophysical properties on composition is correlated with the local chemistry of the solution phase, in particular, the tendency of UCl 3 to form network structures.

36 MATERIALS SCIENCE↗
DOE OSTI2024

KCl-UCl 3 molten salts investigated by Ab Initio Molecular Dynamics (AIMD) simulations: A comparative study with three dispersion models

Ab Initio Molecular Dynamics (AIMD) simulations are performed on molten KCl-UCl 3 salt mixtures to determine energies, heat capacities, and densities. The density-dependent energy correction (DFT-dDsC), Grimme et al.’s DFT-D3, and Langreth & Lundqvist (vdW-cx) models are used for dispersion forces and combined with the Perdew-Burke-Ernzerhof (PBE) exchange-correlation potential with a Hubbard U param eter for the 5$f$ electrons of uranium. After validating predictions for the end-member systems to literature data, KCl-UCl 3 mixtures are studied at select temperatures. Densities and energies both deviate from ideal solution behavior, with the maximum deviation occurring around 36% UCl 3 for mixing energies and slightly lower (29% UCl 3 ) for densities. Compared to the NaCl-UCl 3 system, which was previously investigated using the same simulation methodologies, the KCl-UCl 3 density and mixing energy deviations from ideal solution behavior are larger by almost a factor of two. No deviation from ideal solution behavior for heat capacity was observed. The AIMD predictions for mixing energies and densities agree qualitatively with experimental data, though the spread in data obtained from the various dispersion force models utilized, measurements, and empirical estimates makes strong conclusions difficult. The dependence of thermodynamic and thermophysical properties on composition is correlated with the local chemistry of the solution phase, in particular, the tendency of UCl 3 to form network structures.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Simultaneously improved reversibility and hydrogen production of solid oxide cells through infiltrating air electrode

Among the various fuel cells, solid oxide cells (SOCs) are the unique type that can principally operate reversibly as either fuel cells to produce electricity or as an electrolyser to split water and produce green hydrogen (H 2 ). Nevertheless, the SOCs' reversibility presents enormous challenges that are manifested by the fast degradation through the cycling between fuel cell and electrolysis mode. While the La 0.8 Sr 0.2 MnO 3 /yttria-stabilized zirconia (LSM/YSZ) air electrode possesses significant advantages in terms of high electrical conductivity and high thermal stability under fuel cell mode, the SOCs with the LSM/YSZ air electrode experience rapid performance degradation with catastrophic electrode delamination shortly after switching from fuel cell to electrolysis mode. To prevent such catastrophic delamination and enable the electrolysis H 2 production, a chemical solution with the designed chemistry of SrFe 2 O 4-d was infiltrated into the LSM/YSZ air electrode. The infiltration immediately mitigates the catastrophic delamination, and the infiltrated cells exhibit significantly improved reversibility upon the electrochemical operation. Nanostructure examination reveals nanoscale cracks and second-phase nanograins formed in the air electrode from the baseline cell. By contrast, no delamination was observed at either the micron or the nanoscale for the infiltrated cell, which is attributed to the increased ion conductivity of the Fe-doped LSM mixed conductor induced by the significant interdiffusion between the LSM backbone and infiltrates. Here, this study presents a viable method for preventing electrode delamination while enhancing the durability of H 2 production and power generation for reversible fuel cell/electrolysis cell operation. It further opens new research directions of modifying the electrochemical activity of the electrode of inherently functional cells through the infiltration of the solutions with different chemistry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Probing the self-ionization of liquid water with ab initio deep potential molecular dynamics

The chemical equilibrium between self-ionized and molecular water dictates the acid–base chemistry in aqueous solutions, yet understanding the microscopic mechanisms of water self-ionization remains experimentally and computationally challenging. Herein, Density Functional Theory (DFT)–based deep neural network (DNN) potentials are combined with enhanced sampling techniques and a global acid–base collective variable to perform extensive atomistic simulations of water self-ionization for model systems of increasing size. The explicit inclusion of long-range electrostatic interactions in the DNN potential is found to be crucial to accurately reproduce the DFT free energy profile of solvated water ion pairs in small (64 and 128 H 2 O) cells. The reversible work to separate the hydroxide and hydronium to a distance S is found to converge for simulation cells containing more than 500 H 2 O, and a distance of ∼ 8 Å is the threshold beyond which the work to further separate the two ions becomes approximately zero. The slow convergence of the potential of mean force with system size is related to a restructuring of water and an increase of the local order around the water ions. Calculation of the dissociation equilibrium constant illustrates the key role of long-range electrostatics and entropic effects in the water autoionization process.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Cryogenic Carbon Monoxide Oxidation on Cuprous Oxide

Performing oxidation reactions at low temperatures using earth-abundant materials is crucial for advancing solutions for sustainable chemistry. CO oxidation serves as a benchmark reaction to characterize oxidation and to advance fundamental concepts in surface chemistry. While there are several examples of CO oxidation occurring on metal oxides at low temperatures, from 300 K to ∼200 K, reactivity in the cryogenic temperature regime typically requires a metal nanoparticle on a metal oxide. Here, we show oxygen atoms on the (111) facet of Cu 2 O react with CO to form CO 2 at temperatures below 100 K. Combining spectroscopic experimental evidence with calculations, we propose a low barrier path for CO oxidation at reconstructed surface sites on Cu 2 O(111). This finding is a rare example of an earth-abundant metal oxide, in this case copper, that can provide highly reactive multifunctional sites, enabling both adsorption and reaction fundamental steps toward the efficient heterogeneous oxidation of chemicals.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Revealing Hot and Long-Lived Metastable Spin States in the Photoinduced Switching of Solvated Metallogrid Complexes with Femtosecond Optical and X-ray Spectroscopies

An atomistic understanding of the photoinduced spin-state switching (PSS) within polynuclear systems of d 4 –d 7 transition metal ion complexes is required for their rational integration into light-driven reactions of chemical and biological interest. However, in contrast to mononuclear systems, the multidimensional dynamics of the PSS in solvated molecular arrays have not yet been elucidated due to the expected complications associated with the connectivity between the metal centers and the strong interactions with the surroundings. Herein, the PSS in a solvated triiron(II) metallogrid complex is characterized using transient optical absorption and X-ray emission spectroscopies on the femtosecond time scale. The complementary measurements reveal the photoinduced creation of energy-rich (hot) and long-lived quintet states, whose dynamics differ critically from their mononuclear congeners. We report this finding opens major prospects for developing novel schemes in solution-phase spin chemistry that are driven by the dynamic PSS process in compact oligometallic arrays.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Kinetic Probes of the Origin of Activity in MOF-Based C–H Oxidation Catalysis

The development of homogeneous catalysis is enabled by the availability of a rich toolkit of kinetics experiments, such as the Hgdrop test, that differentiate catalytic activity at ligand-supported metal complexes from potential heterogeneous catalysts derived from decomposition of molecular species. Metal-organic frameworks (MOFs) have garnered significant attention as platforms for catalysis at site-isolated, interstitial catalyst sites. Unlike homogeneous catalysis, relatively few strategies have been advanced to evaluate the origin of catalytic activity in MOF-promoted reactions. Many of the MOFs that have been examined as potential catalysts are comprised of molecular constituents that represent viable catalysts in the absence of the extended MOF lattice, and thus interfacial sites and leached homogeneous species represent potential sources of catalyst activity. Here, we demonstrate that analysis of deuterium kinetic isotope effects (KIEs) and olefin epoxidation diastereoselectivity provides probes of the origin of catalytic activity in MOF-promoted oxidation reactions. These analyses support involvement of lattice-based Fe sites in the turnover-limiting step of C–H activation with Fe-MOF-74-based materials (i.e., the MOF functions as a bona fide catalyst) and evidence that Cu 2 -based MOF MIL-125-Cu 2 O 2 functions as a solid-state initiator for solution-phase oxidation chemistry and is not involved in the turnover limiting step (i.e., the MOF does not function as a catalyst for substrate functionalization). furthermore, we anticipate that the simple experiments described here will provide a valuable tool for clarifying the role of MOFs in C–H oxidation reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

The development and application of the stirred‐reactor coupon analysis (SRCA) test method

A new technique, termed the stirred‐reactor coupon analysis (SRCA) method, has been developed to measure the rate of glass dissolution in forward‐rate conditions. Monolithic glass coupons are partially masked with an inert material before placement in a large volume of well‐mixed solution with known chemistry and temperature for a predetermined duration. After the test, the mask is removed, and the difference in step height between the protected area and the exposed corroded portions of the sample coupon is measured to determine the extent of glass dissolution. The step height is converted to a rate measurement using the test duration and glass density. Test parameters such as sample surface preparation and test duration were evaluated to determine their effects on the measured rates. Additionally, results from an interlaboratory study (ILS) consisting of 12 laboratories from 11 different institutions are presented, where each laboratory performed 12 independent tests. When removing experimental outlier data, the 95% reproducibility limits for the SRCA method has no statistical difference with previously published standardized test methods used to determine the forward rate of glass dissolution. Overall, this paper describes steps necessary to perform the test method and provides the statistical calculations to evaluate test accuracy.

chemical durability↗
DOE OSTI2025

A spline-based method to obtain spatially dependent viscosity in confined flows

Coupling chemical physics to continuum theories is a critical step to understanding multi-scale phenomena. This paper will connect non-equilibrium molecular dynamics simulations to a continuum-based Navier-Stokes equation that has relaxed the assumption of spatial uniformity in viscosity. Using a form for viscosity based on spline interpolation, viscosity as a function of position is obtained from the least squares fit of the velocity profile measured from molecular simulations of flow in a nanochannel. Viscosity can vary widely, particularly near the channel boundaries, indicating that uniform viscosity is no longer appropriate. Variations of the viscosity near the channel surfaces imply that considering solution and surface chemistry could be necessary to rigorously understand molecular-scale flows in nanochannels.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Styrene-Based Elastomer Composites with Functionalized Graphene Oxide and Silica Nanofiber Fillers: Mechanical and Thermal Conductivity Properties

The mechanical and thermal conductivity properties of two composite elastomers were studied. Styrene–butadiene rubber (SBR) filled with functionalized graphene oxide (GO) and silica nanofibers, and styrene–butadiene–styrene (SBS) block copolymers filled with graphene oxide. For the SBR composites, GO fillers with two different surface functionalities were synthesized (cysteamine and dodecylamine) and dispersed in the SBR using mechanical and liquid mixing techniques. The hydrophilic cysteamine-based GO fillers were dispersed in the SBR by mechanical mixing, whereas the hydrophobic dodecylamine-based GO fillers were dispersed in the SBR by liquid mixing. Silica nanofibers (SnFs) were fabricated by electrospinning a sol–gel precursor solution. The surface chemistry of the functionalized fillers was studied in detail. The properties of the composites and the synergistic improvements between the GO and SnFs are presented. For the SBS composites, GO fillers were dispersed in the SBS elastomer at several weight percent loadings using liquid mixing. Characterization of the filler material and the composite elastomers was performed using x-ray photoelectron spectroscopy, x-ray diffraction, transmission electron microscopy, scanning electron microscopy, thermogravimetric analysis, dynamic mechanical analysis, tensile testing, nanoindentation, thermal conductivity and abrasion testing.

36 MATERIALS SCIENCE↗
DOE OSTI2022

BioRT-Flux-PIHM v1.0: a biogeochemical reactive transport model at the watershed scale

Watersheds are the fundamental Earth surface functioning units that connect the land to aquatic systems. Many watershed-scale models represent hydrological processes but not biogeochemical reactive transport processes. This has limited our capability to understand and predict solute export, water chemistry and quality, and Earth system response to changing climate and anthropogenic conditions. Here we present a recently developed BioRT-Flux-PIHM (BioRT hereafter) v1.0, a watershed-scale biogeochemical reactive transport model. The model augments the previously developed RT-Flux-PIHM that integrates land-surface interactions, surface hydrology, and abiotic geochemical reactions. It enables the simulation of (1) shallow and deep-water partitioning to represent surface runoff, shallow soil water, and deeper groundwater and of (2) biotic processes including plant uptake, soil respiration, and nutrient transformation. The reactive transport part of the code has been verified against the widely used reactive transport code CrunchTope. BioRT-Flux-PIHM v1.0 has recently been applied in multiple watersheds under diverse climate, vegetation, and geological conditions. This paper briefly introduces the governing equations and model structure with a focus on new aspects of the model. It also showcases one hydrology example that simulates shallow and deep-water interactions and two biogeochemical examples relevant to nitrate and dissolved organic carbon (DOC). These examples are illustrated in two simulation modes of complexity. One is the spatially lumped mode (i.e., two land cells connected by one river segment) that focuses on processes and average behavior of a watershed. Another is the spatially distributed mode (i.e., hundreds of cells) that includes details of topography, land cover, and soil properties. Whereas the spatially lumped mode represents averaged properties and processes and temporal variations, the spatially distributed mode can be used to understand the impacts of spatial structure and identify hot spots of biogeochemical reactions. The model can be used to mechanistically understand coupled hydrological and biogeochemical processes under gradients of climate, vegetation, geology, and land use conditions.

54 ENVIRONMENTAL SCIENCES↗