Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Effective diffusivity”

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 1,351 records · Page 75

Morphology and Crystallinity Effects of Nanochanneled Niobium Oxide Electrodes for Na-Ion Batteries

Niobium pentoxide (Nb 2 O 5 ) is a promising negative electrode for sodium ion batteries (SIBs). By engineering the morphology and crystallinity of nanochanneled niobium oxides (NCNOs) the kinetic behavior and charge storage mechanism of Nb 2 O 5 electrodes were investigated. Amorphous and crystalline NCNO samples were made by modulating anodization conditions (20 - 40V and 140 - 180 °C) to synthesize nanostructures of differing pore sizes and wall thicknesses with identical chemical composition. The electrochemical energy storage properties of the NCNOs were studied with the amorphous samples showing better overall rate performance than the crystalline samples. The enhanced rate performance of the amorphous samples is attributed to the higher capacitive contributions and Na-ion diffusivity analyzed from cyclic voltammetry (CV) and galvanostatic intermittent titration technique (GITT). It was found that the amorphous samples with smaller wall thicknesses facilitated improved kinetics. Among samples with similar pore size and wall thickness, the difference in their power performance stems from crystallinity effect, which plays a more significant role in the resulting kinetics of the materials for Na-ion batteries.

25 ENERGY STORAGE↗

Irreversible Multielement Diffusion and the Resulting Compositional and Processing Flexibility in the Synthesis and Densification of Lithium Aluminum Lanthanum Zirconium Oxide

Despite the broad interest in Al-doped lithium lanthanum zirconium oxide (LLZO), the wide range of reported ionic conductivities suggests that the effect of processing parameters on the resulting phase purity and bulk conductivity is not well-understood. Here, in this work, we synthesized six separate series of LLZO with variations in Al concentration, Li excess, Li addition order, Li source, densification method, and mother powder to determine the effect on the composition, phase purity, density, and bulk conductivity. We found that a wider range of compositions (6.08–7.61 mol of Li and 0.06–0.23 mol of Al) than previously reported can result in cubic phase stability and that nearly all elements (Li, Al, Zr, and La) are lost to the MgO crucible during calcination and sintering. The manner in which different elements are lost is affected by the processing parameters. We observed that Al-doped LLZO shows great compositional flexibility in stabilizing the cubic phase, offering an explanation for the range of electrochemical performance metrics reported in the literature.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Understanding the Origin of Negative Temperature Dependence and Activity of N-Coordinated Cobalt Sites During Ethylene Dimerization

The on-demand production of short-chain linear alpha olefins (LAOs; C4-C8) via C2H4 dimerization and oligomerization is industrially attractive, prompting extensive research on designing active, selective, and stable catalysts for industrial use. Cobalt supported on ammoniated carbon (Co(NH3)x/C) catalysts have shown remarkable activity and selectivity in this process. However, critical aspects such as the active phase, active site structure, the role of the catalyst support, cobalt loading effects, and the inverse correlation of the reaction rate with temperature remain inadequately understood. This study systematically explores these factors using a combination of steady-state differential catalytic tests, in situ molecular characterization including diffuse reflectance UV-Vis (DR-UV-Vis), Infrared, and Raman spectroscopies, and ex situ X-ray diffraction (XRD) and high annular aberration-corrected dark field transmission electron microscopy (HAADF-STEM). Various supports (SiO2, Al2O3, NH4-ZSM-5, g-C3N4, and C) and cobalt loadings (1.0-3.0 Co nm-2) were studied to determine the optimal catalyst composition and identify the active phase and sites. Carbon-supported catalysts uniquely produce C4-8 LAOs during C2H4 dimerization, with site-time-yield remaining constant (~10-3 s-1) for 1.0-4.0 Co nm-2 at prolonged reaction times (24-48?h time-on-stream). At higher loadings of 6.0 Co nm-2, the formation of crystalline CoO and Co3O4 phases reduces catalytic activity and LAO selectivity. Our findings show that active catalysts lack crystalline cobalt oxides and instead feature dispersed Co2+ sites, tetra-coordinated to a mix of N/NH3 and O/H2O ligands, which catalyze C2H4 dimerization via the Cossee-Arlman mechanism, exhibiting 1st order dependence on C2H4 concentration. The observed inverse rate-temperature correlation is attributed to compensation effects (i.e., presence of Cremer-Constable relationship) linked to changes in adsorption enthalpic and entropic factors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Role of Wadsley Defects and Cation Disorder to Enhance MoNb 12 O 33 Diffusion

Wadsley-Roth (WR) niobates have emerged as high-rate anode materials that can combine rapid ionic diffusion with good electronic conductivity. WR compounds have been defect-enhanced by limited annealing, however, such materials often contain multiple types of defects. In particular, both Wadsley defects (variable block size) and transition metal disorder have the potential to modify transport rates, however the corresponding effects are not well understood mechanistically. Here, MoNb 12 O 33 (MNO) was calcined at two different temperatures to compare a defect-rich condition (MNO-800) with a proximal order-rich condition (MNO-900) as assessed through XRD, XANES, EXAFS, and STEM characterizations. Galvanostatically cycled lithium half cells of MNO-800 exhibited additional capacity (307 mAhg −1 at 0.1C, 4.66% higher) and improved high-rate capacity of 200 mAhg −1 at 10C. ICI-based overpotential analysis identified solid state diffusion as the dominant rate limiting process where MNO-800 correspondingly exhibited ∼3X faster capacity-weighted diffusivity. A machine-learning interatomic potential was trained to density functional theory and then applied with molecular dynamics (MLIP-MD) to examine the possible roles of Wadsley defects and transition metal disorder. For both defect-types, Li was found to populate and activate fast diffusion paths from window sites at lower extents of lithiation as compared to the order-rich model.

defect↗

Development of Corrosion Resistant Coatings for Structural Materials for Liquid Fueled Molten Salts Reactors Applications (Final Report)

Current structural alloys, code certified for the MSRs temperature ranges, contain high levels of chromium, making them highly susceptible to molten halide salt corrosion. One potential solution to circumvent the need for code certification of novel alloys, which is expensive and time consuming, is to design claddings that will protect the underlying code certified materials from corrosion damage during operation. In this work, we examined the corrosion of Ni and Cu electroplated, Ni and CuNi weld overlay, Mo-laser clad, and carburized claddings on SS316H for use in molten salt reactor environments. After characterization of the cladded materials (Task 1), static corrosion tests were used to assess corrosion resistance of the claddings in typical molten fluoride salts (Task 2). The corrosion tests were performed in molten FLiNaK at 700°C up to 1000 hours. Pre- and post-corrosion Scanning Electron Microscopy (SEM), energy-dispersive-spectroscopy (EDS), Transmission Electron Microscopy (TEM), X-ray Diffraction (XRD) and glow-discharge-optical-emission-spectroscopy (GDOES) were performed on cladding cross-sections and surfaces to evaluate degradation. The Cu and Ni electroplated samples, as well as the carburized samples, showed excellent corrosion resistance relative to the bare SS316H. To assess high temperature cladding stability, ageing experiments were performed at temperatures up to 900 C in inert atmosphere for the electroplated and carburized samples and a diffusion model was developed to predict long term cladding behavior (Task 3). It was found that the Cu cladding was basically insensitive to the high temperature ageing, except for small secondary phases forming at the interface. On the other hand, the Ni electroplated cladding experienced significant interdiffusion. Nevertheless, the gain in corrosion resistance for a 100μm Ni electroplated cladding is phenomenal, with more than 50% reduction in chromium dissolution from the substrate material for the first 15 years of salt exposure at 700°C. To assess radiation resistance and phase stability, high-temperature 4MeV Ni heavy ion irradiation was performed across the cladding/substrate interface up to 50 displacement-per-atom (DPA) at 500°C and 700°C to assess the phase stability and irradiation behaviors of the electroplated systems (Task 4). The Ni and Cu electroplated systems did not experience void swelling at the contrary to the SS316H substrate due to their nanocrystalline nature. Some level of recrystallization was observed in the cladding, as well as radiation induced segregation and enhanced diffusion. The interface acted as a potent sink for point defects with the presence of a void denuded zone. The mechanical properties of the claddings were assessed using thermal shock, micro-indentation, nano-indentation, and four-point bend testing experiments (Task 5). These experiments were performed on the claddings in as-received, corroded, irradiated, and thermally aged states to determine the effects of typical molten salt reactor environments on cladding integrity. No significant mass loss was observed after repeated thermal shocks. While the electroplated samples softened after high-temperature ageing, irradiation hardening compensate this effect, such that there is little different with the as-received materials. The results of these experiments suggest that the electroplated (copper/nickel) show the greatest promise for application in molten salt reactor development. These claddings prevented any chromium dissolution from occurring during the static corrosion experiments. Additionally, they demonstrated favorable properties for high temperature diffusion, phase stability, interfacial mechanical properties, and irradiation resistance. Weld-overlay cladding are also of great interest and could reach properties similar to the electroplated samples upon optimization (especially with multiple weld passes). The carburized SS316H showed excellent behavior as well, but more work needs to be done to assess their long-term stability. Finally, the Mo-laser clad system was not pursued further due to the manufacturing challenges. It is also worth noting that the Ni cladding systems should behave relatively well in terms of weldability since the Ni-weld overlay microstructure and associated corrosion rate are sound. While the NiCu weld-overlay has even lower corrosion rates than the Ni-weld overlay, more studies on welding of Cu-electroplated systems are necessary since Cu clusters are known to embrittle steels.

36 MATERIALS SCIENCE↗

Polymer-based microfluidic device for on-chip counter-diffusive crystallization and in situ X-ray crystallography at room temperature

Proteins are long chains of amino acid residues that perform a myriad of functions in living organisms, including enzymatic reactions, signalling, and maintaining structural integrity. Protein function is determined directly by the protein structure. X-ray crystallography is the primary technique for determining the 3D structure of proteins, and facilitates understanding the effects of protein structure on function. The first step towards structure determination is crystallizing the protein of interest. We have developed a centrifugally-actuated microfluidic device that incorporates the fluid handling and metering necessary for protein crystallization. Liquid handling takes advantage of surface forces to control fluid flow and enable metering, without the need for any fluidic or pump connections. Our approach requires only the simple steps of pipetting the crystallization reagents into the device followed by either spinning or shaking to set up counter-diffusive protein crystallization trials. The use of thin, UV-curable polymers with a high level of X-ray transparency allows for in situ X-ray crystallography, eliminating the manual handling of fragile protein crystals and streamlining the process of protein structure analysis. We demonstrate the utility of our device using hen egg white lysozyme as a model system, followed by the crystallization and in situ, room temperature structural analysis of the hub domain of calcium–calmodulin dependent kinase II (CaMKIIβ).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Digital quantum magnetism on a trapped-ion quantum computer

Digital quantum matter—realized when discrete quantum gates approximate continuous time evolution—is susceptible to heating into chaotic, structureless states. If digitization errors are adequately suppressed, a long-lived transient regime of approximately energy-conserving dynamics can be observed on gate-based quantum computers. Conservation of energy, in turn, enables the exploration of a wide variety of complex behaviours observed in equilibrium systems, ranging from the non-trivial microscopic origins of thermalization itself to the stabilization of effective models hosting exotic emergent properties. Here we use Quantinuum’s H2 quantum computer to simulate digitized dynamics of the quantum Ising model, suppressing digitization errors well enough to observe thermalization on timescales that severely challenge classical simulation methods. Relaxation of an inhomogeneous state reveals an emergent hydrodynamics owing to approximate energy conservation and we compute the associated diffusion constant. By reprogramming our simulations to take place on a triangular lattice with periodic boundary conditions, we observe thermalization consistent with emergent gauge and topological constraints resulting from lattice frustration. Furthermore, our results were enabled by continued advances in two-qubit gate quality (native partial entangler fidelities of 99.94(1)%) and establish digital quantum computers as powerful tools for studying (effectively) continuous-time dynamics.

Information theory and computation↗

Control of Two Solid Electrolyte Interphases at the Negative Electrode of an Anode‐Free All Solid‐State Battery based on Argyrodite Electrolyte

Abstract Anode‐free all solid‐state batteries (AF‐ASSBs) employ “empty” current collector with three active interfaces that determine electrochemical stability; lithium metal – Solid electrolyte (SE) interphase (SEI‐1), lithium – current collector interface, and collector – SE interphase (SEI‐2). Argyrodite Li 6 PS 5 Cl (LPSCl) solid electrolyte (SE) displays SEI‐2 containing copper sulfides, formed even at open circuit. Bilayer of 140 nm magnesium/30 nm tungsten (Mg/W‐Cu) controls the three interfaces and allows for state‐of‐the‐art electrochemical performance in half‐cells and fullcells. AF‐ASSB with NMC811 cathode achieves 150 cycles with Coulombic efficiency (CE) above 99.8%. With high mass‐loading cathode (8.6 mAh cm −2 ), AF‐ASSB retains 86.5% capacity after 45 cycles at 0.2C. During electrodeposition of Li, gradient Li‐Mg solid solution is formed, which reverses upon electrodissolution. This promotes conformal wetting/dewetting by Li and stabilizes SEI‐1 by lowering thermodynamic driving force for SE reduction. Inert refractory W underlayer is required to prevent ongoing formation of SEI‐2 that also drives electrochemical degradation. Inert Mo and Nb layers likewise protect Cu from corroding, while Li‐alloying layers (Mg, Sn) are less effective due to ongoing volume changes and associated pulverization. Mechanistic explanation for observed Li segregation within alloying Li x Mg layer is provided through mesoscale modelling, considering opposing roles of diffusivity differences and interfacial stresses.

Wang, Yixian [Materials Science and Engineering Pr↗

Grazing-Incidence Texture Tomography and Diffuse Reflectivity Tomography of an Organic Semiconductor Device Array

The use of grazing-incidence scattering methods for the characterization of 2D patterned organic thin films is limited due to the elongated 1D footprint of the X-ray beam on the sample. However, this characteristic feature can be turned into an advantage, when combined with tomographic reconstruction. In this pilot study we show, how the use of a chosen texture reflection and a diffuse reflectivity signal can each provide 2D images of the deposits, simultaneously revealing the organic film's crystal orientation and the location of the metal electrodes in a field-effect transistor structure from a single sequence of diffraction images.

36 MATERIALS SCIENCE↗

Computational Methodology to Simulate Pyrometallurgical Processes in a Secondary Lead Furnace

Pyrometallurgical recovery of nonferrous metals involves a combination of thermally intensive transformations during exothermic gas-phase reactions, endothermic decomposition of solid charge, and melting of simpler solids. In the recovery of secondary lead, simultaneous thermal effects in a reverberatory-style furnace cause a melt pool to accumulate at the bottom, with lighter solids (slag) floating above and gaseous products from decomposition of the charge diffusing through the gas–slag interface. Species from oxy-fuel combustion of natural gas, species profiles from smelting reactions, and the formation of a melt pool consisting primarily of lead are simulated via a time-averaged formulation. Predictions of outflow are compared with preset inflow profiles to ensure conservation of mass. Thermal profiles for solid, liquid, and gas phases are presented by species. A novel method is implemented to model the latent heat of fusion using a heterogeneous chemical reaction. The simulation is conducted in Simcenter STAR-CCM+ v. 16.02.009-R8.

Rao, Vivek↗

Limiting Current in Nanostructured Block Copolymer Electrolytes

The next-generation electrolytes for lithium batteries must be able to conduct ions at sufficiently high current densities; yet this regime is rarely studied directly. The limiting current density of an electrolyte quantifies the highest possible rate of ion transport under an applied dc potential. Herein, we report on the limiting current density in twelve nanostructured polystyrene-block-poly(ethylene oxide) (PS-b-PEO, or SEO) copolymer electrolytes. We find that the limiting current at a given salt concentration increases systematically with increasing volume fraction of the PEO block (φ EO ). In contrast, the effective-medium theory, commonly used to analyze conductivity in block copolymer electrolytes, predicts that limiting current is independent of φ EO . To resolve this conundrum, the ionic conductivity, the mutual diffusion coefficient of the salt, and the steady-state current fraction of the block copolymer electrolytes were measured. These measurements enable predictions of limiting current with no adjustable parameters using the concentrated solution theory. We found quantitative agreement between experimentally measured limiting current densities and predictions based on the concentrated solution theory. This work sheds light on how to reliably measure and predict limiting current density in composite electrolytes.

25 ENERGY STORAGE↗

Diversity-oriented synthesis of polymer membranes with ion solvation cages

Microporous polymers feature shape-persistent free volume elements (FVEs), which are permeated by small molecules and ions when used as membranes for chemical separations, water purification, fuel cells, and batteries. It remains a significant challenge to identify FVEs with analyte specificity, due to difficulties in generating microporous polymer libraries with sufficient diversity for screening their properties. Here, we describe a diversity-oriented synthetic (DOS) strategy for microporous polymer membranes from which we identified those whose FVEs serve as solid solvation cages for lithium ions (Li + ). Furthermore, key elements of our strategy included diversification of bis(catechol)-type monomers via multi-component Mannich reactions to introduce Li+-coordinating functionality within individual FVEs, topology-enforcing polymerizations for generating macromolecular skeletal diversity for networking FVEs into different pore architectures, and several classes of on-polymer reactions for diversifying pore geometries and dielectric properties. Lead candidate polymer membranes featuring explicit ion solvation cages exhibited both higher ionic conductivity and higher cation transference number than control membranes where FVEs were aspecific, which indicates conventional bounds for membrane permeability and selectivity for ion transport can be overcome.4 These advantages are tied to enhanced Li + partitioning from the electrolyte when the cages are present, higher diffusion barriers for anions within the pores, and network-enforced restrictions on the number of solvent molecules bound to Li+ by comparison to the bulk electrolyte, which reduces the effective mass of the working ion. Such membranes show promise as anode-stabilizing interlayers in high-voltage lithium-metal batteries for electric mobility.

42 ENGINEERING↗

Unconventional short-range structural fluctuations in cuprate superconductors

Abstract The interplay between structural and electronic degrees of freedom in complex materials is the subject of extensive debate in physics and materials science. Particularly interesting questions pertain to the nature and extent of pre-transitional short-range order in diverse systems ranging from shape-memory alloys to unconventional superconductors, and how this microstructure affects macroscopic properties. Here we use neutron and X-ray diffuse scattering to uncover universal structural fluctuations in La 2-x Sr x CuO 4 and Tl 2 Ba 2 CuO 6+δ , two cuprate superconductors with distinct point disorder effects and with optimal superconducting transition temperatures that differ by more than a factor of two. The fluctuations are present in wide doping and temperature ranges, including compositions that maintain high average structural symmetry, and they exhibit unusual, yet simple scaling behaviour. The scaling regime is robust and universal, similar to the well-known critical fluctuations close to second-order phase transitions, but with a distinctly different physical origin. We relate this behaviour to pre-transitional phenomena in a broad class of systems with structural and magnetic transitions, and propose an explanation based on rare structural fluctuations caused by intrinsic nanoscale inhomogeneity. We also uncover parallels with superconducting fluctuations, which indicates that the underlying inhomogeneity plays an important role in cuprate physics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

The genetic architecture of leaf stable carbon isotope composition in Zea mays and the effect of transpiration efficiency on leaf elemental accumulation

With increased demand on freshwater resources for agriculture, it is imperative that more water-use efficient crops are developed. Leaf stable carbon isotope composition, δ 13 C, is a proxy for transpiration efficiency and a possible tool for breeders, but the underlying mechanisms effecting δ 13 C in C 4 plants are not known. It has been suggested that differences in specific leaf area (SLA), which potentially reflects variation in internal CO 2 diffusion, can impact leaf δ 13 C. Furthermore, although it is known that water movement is important for elemental uptake, it is not clear how manipulation of transpiration for increased water-use efficiency may impact nutrient accumulation. Here, we characterize the genetic architecture of leaf δ 13 C and test its relationship to SLA and the ionome in five populations of maize. Five significant QTL for leaf δ 13 C were identified, including novel QTL as well as some that were identified previously in maize kernels. One of the QTL regions contains an Erecta-like gene, the ortholog of which has been shown to regulate transpiration efficiency and leaf δ 13 C in Arabidopsis. QTL for δ 13 C were located in the same general chromosome region, but slightly shifted, when comparing data from two different years. Our data does not support a relationship between δ 13 C and SLA, and of the 19 elements analyzed, only a weak correlation between molybdenum and δ 13 C was detected. Together these data add to the genetic understanding of leaf δ 13 C in maize and suggest that improvements to plant water use may be possible without significantly influencing elemental homeostasis.

59 BASIC BIOLOGICAL SCIENCES↗

Role of carbon and hydrogen in limiting n -type doping of monoclinic ( Al x Ga 1 – x ) 2 O 3

In this study, we use hybrid density functional calculations to assess n -type doping in monoclinic ( Al x Ga 1 – x ) 2 O 3 alloys. We focus on silicon, the most promising donor dopant, and study the structural properties, formation energies, and charge-state transition levels of its various configurations. We also explore the impact of carbon and hydrogen, which are common impurities in metal-organic chemical vapor deposition (MOCVD). In Ga 2 O 3 , Si Ga is an effective shallow donor, but in Al 2 O 3 Si Al acts as a DX center with a ( + / – ) transition level in the band gap. Interstitial hydrogen acts as a shallow donor in Ga 2 O 3 but behaves as a compensating acceptor in n -type Al 2 O 3 . Interpolation indicates that Si is an effective donor in ( Al x Ga 1 – x ) 2 O 3 up to 70% Al, but it can be compensated by hydrogen already at 1% Al. We also assess the diffusivity of hydrogen and study complex formation. Si cation – H complexes have relatively low binding energies. Substitutional carbon on a cation site acts as a shallow donor in Ga 2 O 3 , but can be stable in a negative charge state in ( Al x Ga 1 – x ) 2 O 3 when x > 5 % . Substitutional carbon on an oxygen site ( C O ) always acts as an acceptor in n -type ( Al x Ga 1 – x ) 2 O 3 , but will incorporate only under relatively oxygen-poor conditions. C O – H complexes can actually incorporate more easily, explaining observations of carbon-related compensation in Ga 2 O 3 grown by MOCVD. We also investigate C cation – H complexes, finding they have high binding energies and act as compensating acceptors when x > 56 % ; otherwise the hydrogen just passivates the unintentional carbon donors. C-H complex formation explains why MOCVD-grown Ga 2 O 3 can exhibit record-low free-carrier concentrations, in spite of the unavoidable incorporation of carbon. Our study highlights that, while Si is in principle a suitable shallow donor in ( Al x Ga 1 – x ) 2 O 3 alloys up to high Al compositions, control of unintentional impurities is essential to avoid compensation.

36 MATERIALS SCIENCE↗

High-dimensional and permutation invariant anomaly detection

Methods for anomaly detection of new physics processes are often limited to low-dimensional spaces due to the difficulty of learning high-dimensional probability densities. Particularly at the constituent level, incorporating desirable properties such as permutation invariance and variable-length inputs becomes difficult within popular density estimation methods. In this work, we introduce a permutation-invariant density estimator for particle physics data based on diffusion models, specifically designed to handle variable-length inputs. We demonstrate the efficacy of our methodology by utilizing the learned density as a permutation-invariant anomaly detection score, effectively identifying jets with low likelihood under the background-only hypothesis. To validate our density estimation method, we investigate the ratio of learned densities and compare to those obtained by a supervised classification algorithm.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Mixing Cell Model: A One-Dimensional Numerical Model for Assessment of Water Flow and Contaminant Transport in the Unsaturated Zone

This report describes the Mixing Cell Model code, a one dimensional model for water flow and solute transport in the unsaturated zone under steady state or transient flow conditions. The model is based on the principles and assumptions underlying mixing-cell model formulations. The unsaturated zone is discretized into a series of independent mixing cells. Each cell may have unique hydrologic, lithologic, and sorptive properties. Ordinary differential equations describe the material (water and solute) balance within each cell. Water-flow equations are derived from the continuity equation, assuming that unit gradient conditions exist at all times in each cell. Pressure gradients are considered implicitly through model discretization. Unsaturated hydraulic conductivity and moisture contents are determined by the material specific moisture-characteristic curves. Solute-transport processes include explicit treatment of advective processes, first order chain decay, and linear sorption reactions. Dispersion is addressed through implicit and explicit dispersion. Implicit dispersion is an inherent feature of all mixing-cell models and originates from the formulation of the problem in terms of mass balance around fully mixed volume elements. Expressions are provided that relate implicit dispersion to the physical dispersion of the system. Two FORTRAN codes were developed to solve the water flow and solute-transport equations: (1) the Mixing Cell Model for Flow (MCMF) solves transient water-flow problems and (2) the Mixing Cell Model for Transport (MCMT) solves the solute-transport problem. The transient water-flow problem is typically solved first by estimating the water flux through each cell in the model domain as a function of time using the MCMF code. These data are stored in either ASCII or binary files that are later read by the solute transport code (MCMT). Code output includes solute pore water concentrations, water and solute inventories in each cell and at each specified output time, and water and solute fluxes through each cell and specified output time. Computer run times for coupled transient water flow and solute transport were typically several seconds on a 2 GHz Intel Pentium IV desktop computer. The model was benchmarked against analytical solutions and finite element approximations to the partial differential equations (PDE) describing unsaturated flow and transport. Differences between the maximum solute flux estimated by the mixing cell model and the PDE models were typically less than two percent. This revision includes an option for a fixed concentration lower boundary condition for diffusive fluxes for versions 020321 and later.

54 ENVIRONMENTAL SCIENCES↗

Spin Susceptibility in Neutron Matter from Quantum Monte Carlo Calculations

The spin susceptibility in pure neutron matter is computed from auxiliary field diffusion Monte Carlo calculations over a wide range of densities. The calculations are performed for different spin asymmetries, while using twist-averaged boundary conditions to reduce finite-size effects. The employed nuclear interactions include both the phenomenological Argonne AV8′ + UIX potential and local interactions that are derived from chiral effective field theory up to next-to-next-to-leading order.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗