Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “boundary effects”

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,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Size effect in stainless steel thin wires under tension

Miniaturization of modern devices drives the exploration into how materials behave in small dimensions. However, the mechanical behaviors of micro-scale specimens under varied loading conditions are still challenging to characterize and predict. In this work, specimen size-dependent weakening effect was experimentally revealed in 316L stainless steel thin wires under tension with wire diameter decreasing from 200 μm down to 61 μm. The underlying mechanisms were primarily associated with surface grain softening and strain localization, which were largely dependent on the number of grains across wire diameter. Through an analog composite model with surface and interior grains that follows the Hall-Petch relation, the extent of surface grain softening was found to coincide with the flow stress reduction due to lack of boundary strengthening. Further using elasto-plastic self-consistent polycrystal modeling, the stress-strain behaviors of interior and surface grains were well described in analogy to those of bulk materials with similar-size grains and extremely large grains, respectively. Subsequently, the analog composite model well simulated the specimen size-dependent weakening effect under tension. This article provided more insights into the size effect exploration and proposed a practical model for describing size-dependent mechanical behaviors of small-scale materials to be applied to miniaturized devices.

316L stainless steel↗
DOE OSTI2020

Mechanisms of energetic-particle transport in magnetically confined plasmas

Super-thermal ions and electrons occur in both space and fusion plasmas. Because these energetic particles (EP) have large velocities, EP orbits necessarily deviate substantially from magnetic surfaces. Orbits are described by conserved constants of motion that define topological boundaries for different orbit types. Electric and magnetic field perturbations produced by instabilities can disrupt particle orbits, causing the constants of motion to change. The statistics of the “kicks” associated with these perturbations determines the resulting cross field transport. A unifying theme of this tutorial is the importance of the perturbation’s phase at the particle’s position Θ = k • r – ωt, where k and ω are the wavevector and frequency of the perturbation, r is the EP position, and t is the time. A distinction is made between field perturbations that resonate with an aspect of the orbital motion and those that do not. Resonance occurs when the wave phase returns to its initial value in an integer multiple of an orbital period. Convective transport occurs when resonant particles experience an unvarying wave phase. Alternatively, multiple wave-particle resonances usually decorrelate the phase, resulting in diffusive transport. Large orbits increase the number of important resonances and can cause chaotic orbits even for relatively small amplitude waves. In contrast, in the case of non-resonant perturbations, orbital phase averaging reduces transport. Here, large field perturbations introduce additional effects, including nonlinear resonances at fractional values of the orbital motion. In summary, large orbits are a blessing and a curse: For non-resonant modes, orbit-averaging reduces transport but, for resonant transport, large orbits facilitate jumps across topological boundaries and enhance the number of important resonances.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Half-life determination of heavy ions in a storage ring considering feeding and depleting background processes

Heavy-ion storage rings have relatively large momentum acceptance which allows for multiple ion species to circulate at the same time. This needs to be considered in radioactive decay measurements of highly charged ions, where atomic charge exchange reactions can significantly alter the intensities of parent and daughter ions. In this study, we investigate this effect using the decay curves of ion numbers in the recent 205 Tl 81+ bound-state beta decay experiment conducted using the Experimental Storage Ring at GSI Darmstadt. To understand the intricate dynamics of ion numbers, we present a set of differential equations that account for various atomic and nuclear reaction processes—bound state beta decay, atomic electron recombination and capture, and electron ionization. By incorporating appropriate boundary conditions, we develop a set of differential equations that accurately simulate the decay curves of various simultaneously stored ions in the storage ring: 205 Tl 81+ , 205 Pb 81+ , 205 Pb 82+ , 200 Hg 79+ , and 200 Hg 80+ . Through a quantitative comparison between simulations and experimental data, we provide insights into the detailed reaction mechanisms governing stored heavy ions within the storage ring. Our approach effectively models charge-changing processes, reduces the complexity of the experimental setup, and provides a simpler method for measuring the decay half-lives of highly charged ions in storage rings.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Investigating coastal backwater effects and flooding in the coastal zone using a global river transport model on an unstructured mesh

Abstract. Coastal backwater effects are caused by the downstream water level increase as a result of elevated sea level, high river discharge and their compounding influence. Such effects have crucial impacts on floods in densely populated regions but have not been well represented in large-scale river models used in Earth system models (ESMs), partly due to model mesh deficiency and oversimplifications of river hydrodynamics. Using two mid-Atlantic river basins as a testbed, we perform the first attempt to simulate the backwater effects comprehensively over a coastal region using the MOSART river transport model under an ESM framework, i.e., Energy Exascale Earth System Model (E3SM) configured on a regionally refined unstructured mesh, with a focus on understanding the backwater drivers and their long-term variations. By including sea level variations at the river downstream boundary, the model performance in capturing backwaters is greatly improved. We also propose a new flood event selection scheme to facilitate the decomposition of backwater drivers into different components. Our results show that while storm surge is a key driver, the influence of extreme discharge cannot be neglected, particularly when the river drains to a narrow river-like estuary. Compound flooding, while not necessarily increasing the flood peaks, exacerbates the flood risk by extending the duration of multiple coastal and fluvial processes. Furthermore, our simulations and analysis highlight the increasing strength of backwater effects due to sea level rise and more frequent storm surge during 1990–2019. Thus, backwaters need to be properly represented in ESMs to improve the predictive understanding of coastal flooding.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Combined Effect of the Wegener–Bergeron–Findeisen Mechanism and Large Eddies on Microphysics of Mixed-Phase Stratiform Clouds

Abstract The process of glaciation in mixed-phase stratiform clouds was investigated by a novel Lagrangian–Eulerian model (LEM) in which thousands of adjoining Lagrangian parcels moved within a turbulent-like velocity field with statistical parameters typical of the Arctic boundary layer. We used detailed bin microphysics to describe the condensation/evaporation processes in each parcel, in which droplets, aerosols, and ice particles were described using size distributions of 500 mass bins. The model also calculated aerosol mass inside droplets and ice particles. Gravitational sedimentation of droplets and ice particles was also accounted for. Assuming that droplet freezing is the primary source of ice particles, the Arctic clouds observed in Indirect and Semi-Direct Aerosol Campaign (ISDAC) were successfully simulated. The model showed that at a low ice particle concentration typical of ISDAC, large vortices (eddies) led to a quasi-stationary regime, in which mixed-phase St existed for a long time. The large eddies controlled the water partitioning in the mixed-phase clouds. Droplets formed and grew in updrafts, typically reaching the cloud top, and evaporated in downdrafts. Ice particles grew in updrafts and downdrafts. The Wegener–Bergeron–Findeisen (WBF) mechanism was efficient in downdrafts and some parts of updrafts, depending on ice concentration and vertical velocity. At low ice concentrations, the effect of ice on the phase partitioning was negligible. In this regime, liquid droplets were found near the cloud top, whereas ice particles precipitated through the cloud base. When ice concentration exceeded about 10 L −1 , the WBF mechanism led to glaciation of almost the entire cloud, with the exception of narrow cloud regions associated with strong updrafts. At ice particle concentrations of a few tens per liter, the oscillatory regime took place due to the ice–liquid interaction. The microphysical structure of mixed-phase St forms as a combined effect of cloud dynamics (large eddies) and the WBF mechanism.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Multiscale Understanding of Surface Structural Effects on High-Temperature Operational Resiliency of Layered Oxide Cathodes

The worldwide energy demand in electric vehicles and the increasing global temperature have called for development of high-energy and long-life lithium-ion batteries (LIBs) with improved high temperature operational resiliency. However, current attention has been mostly focused on cycling aging at elevated temperature, leaving considerable gaps of knowledge in the failure mechanism and practical control of abusive calendar aging and thermal runaway that are highly related to the eventual operational lifetime and safety performance of LIBs. Herein, using a combination of various in situ synchrotron X-ray and electron microscopy techniques, we report a multiscale understanding of surface structure effects involved in regulating the high temperature operational tolerance of polycrystalline Ni-rich layered cathodes. Our results collectively show that an ultraconformal poly (3,4-ethylenedioxythiophene) coating can effectively prevent LiNi 0.8 Co 0.1 Mn 0.1 O 2 cathode from undesired phase transformation and transition metal dissolution on the surface, atomic displacement and dislocations within primary particles, intergranular cracking along the grain boundaries within secondary particles, and intensive bulk oxygen release during high state-of-charge and high temperature aging. As a result, the present work highlights the essential role of surface structure controls in overcoming the multiscale degradation pathways of high-energy battery materials at extreme temperature.

25 ENERGY STORAGE↗
DOE OSTI2022

Intergranular ductile failure of materials with plastically heterogeneous grains

In several precipitation hardened alloys that are susceptible to intergranular ductile failure, precipitation does not always occur uniformly throughout the microstructure, and regions close to grain boundaries may remain precipitate-free. These precipitate-free zones (PFZs) in the material microstructures result in plastically heterogeneous grains, since PFZs are expected to have lower yield strength but higher strain-hardenability compared to the precipitate containing grain interior. Experimentally, the presence of PFZs in precipitation-hardened alloys has been associated with both increase and decrease in materials’ ductility, with or without significant change in the strength. Thus, to understand and rationalize the experimental observations, we have carried out extensive microstructure-based finite element calculations of intergranular ductile failure in materials under tensile loading conditions. In the calculations, both the grain boundaries and PFZs are discretely modeled, and a wide range of the values of yield strength, strain-hardenability and width of PFZs in the material microstructures are analyzed. Our results show that the effects of PFZs on the overall mechanical response of the material strongly depend on the values of yield strength and strain-hardenability of PFZs. Here, there exists an optimum combination of the values of these two parameters that can result in the overall ductility and tensile strength of the material microstructures with PFZs being greater than the microstructures without PFZs.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Origin of Disorder Tolerance in Piezoelectric Materials and Design of Polar Systems

Current high-performing piezoelectric materials are dominated by perovskites that rely on soft optical phonon modes stabilized by disorder near a morphotropic phase boundary and a unique resilience of the polar response to that disorder. To identify structural families with similar resilience, in this study we develop a first-principles sensitivity analysis approach to determine the effect of disorder on the piezoelectric response for structures in the Materials Project database. In well-known piezoelectric systems, the lattice dynamics, rather than internal strain or dielectric, control the polar response. Additionally, multiple stable optical phonon modes are found to contribute to the piezoelectric response, providing a fingerprint for disorder tolerance. A multiple-phonon mode criterion is used to evaluate candidate materials for disorder-tolerant piezoelectric prototype systems. Five promising structures are altered through chemical substitution, generating potential MPB end points with large piezoelectric responses beyond perovskites including Akermanite Sr 2 x Ca 2 – 2 x CoSi 2 O 7 , which exhibits a nearly 20% increase in response at the 50% composition.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Toward an Envelope of Design Solutions for Combined/Intensified Reaction/Separation Systems

In this work, we develop an envelope of design solutions for combined intensified reaction/separation systems. The attainable region-based continuous flow stirred tank reactor equivalence principle is applied to characterize design boundaries for a given chemistry independent of process design. The thermodynamic-based Generalized Modular Representation Framework is then employed to generate candidate process alternatives along the design boundaries. The proposed approach is showcased via a case study on olefin metathesis. As a result, we also highlight the need to incorporate the attainable region-based constraints into synthesis strategies to assist effective bounding of design/intensification space.

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

Correlating Crystallographic Orientation and Ferroic Properties of Twin Domains in Metal Halide Perovskites

Metal halide perovskite (MHP) solar cells have attracted worldwide research interest. Although it has been well established that grain, grain boundary, and grain facet affect MHPs optoelectronic properties, less is known about subgrain structures. Recently, MHP twin stripes, a subgrain feature, have stimulated extensive discussion due to the potential for both beneficial and detrimental effects of ferroelectricity on optoelectronic properties. Furthermore, connecting the ferroic behavior of twin stripes in MHPs with crystal orientation will be a vital step to understand the ferroic nature and the effects of twin stripes. In this work, we studied the crystallographic orientation and ferroic properties of CH 3 NH 3 PbI 3 twin stripes, using electron backscatter diffraction (EBSD) and advanced piezoresponse force microscopy (PFM), respectively. Using EBSD, we discovered that the orientation relationship across the twin walls in CH 3 NH 3 PbI 3 is a 90° rotation about ($\bar{11}0$), with the (030) and (111) directions parallel to the direction normal to the surface. By careful inspection of CH 3 NH 3 PbI 3 PFM results including in-plane and out-of-plane PFM measurements, we demonstrate some nonferroelectric contributions to the PFM responses of this CH 3 NH 3 PbI 3 sample, suggesting that the PFM signal in this CH 3 NH 3 PbI 3 sample is affected by nonferroelectric and nonpiezoelectric forces. If there is piezoelectric response, it is below the detection sensitivity of our interferometric displacement sensor PFM (<0.615 pm/V). Overall, this work offers an integrated picture describing the crystallographic orientations and the origin of PFM signal of MHPs twin stripes, which is critical to understanding the ferroicity in MHPs.

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

Tunable magnons of an antiferromagnetic Mott insulator via interfacial metal-insulator transitions

Antiferromagnetic insulators present a promising alternative to ferromagnets due to their ultrafast spin dynamics essential for low-energy terahertz spintronic device applications. Magnons, i.e., quantized spin waves capable of transmitting information through excitations, serve as a key functional element in this paradigm. However, identifying external mechanisms to effectively tune magnon properties has remained a major challenge. Here we demonstrate that interfacial metal-insulator transitions offer an effective method for controlling the magnons of Sr 2 IrO 4 , a strongly spin-orbit coupled antiferromagnetic Mott insulator. Resonant inelastic x-ray scattering experiments reveal a significant softening of zone-boundary magnon energies in Sr 2 IrO 4 films epitaxially interfaced with metallic 4 d transition-metal oxides. Therefore, the magnon dispersion of Sr 2 IrO 4 can be tuned by metal-insulator transitions of the 4 d transition-metal oxides. We tentatively attribute this non-trivial behavior to a long-range phenomenon mediated by magnon-acoustic phonon interactions. Our experimental findings introduce a strategy for controlling magnons and underscore the need for further theoretical studies to better understand the underlying microscopic interactions between magnons and phonons.

magnetic properties and materials↗
DOE OSTI2022

Crystallographic effects on transgranular chloride-induced stress corrosion crack propagation of arc welded austenitic stainless steel

Abstract The effect of crystallography on transgranular chloride-induced stress corrosion cracking (TGCISCC) of arc welded 304L austenitic stainless steel is studied on >300 grains along crack paths. Schmid and Taylor factor mismatches across grain boundaries (GBs) reveal that cracks propagate either from a hard to soft grain, which can be explained merely by mechanical arguments, or soft to hard grain. In the latter case, finite element analysis reveals that TGCISCC will arrest at GBs without sufficient mechanical stress, favorable crystallographic orientations, or crack tip corrosion. GB type does not play a significant role in determining TGCISCC cracking behavior nor susceptibility. TGCISCC crack behaviors at GBs are discussed in the context of the competition between mechanical, crystallographic, and corrosion factors.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Reflection and refraction of directrons at the interface

Reflection and refraction are ubiquitous phenomena with extensive applications, yet minimizing energy loss and information distortion during these processes remains a significant challenge. This study examines the behavior of structurally stable solitons, known as directrons, in nematic liquid crystals interacting with an interface where the director field orientation changes, despite identical physical properties, external potentials, and boundary anchoring in the two regions. During reflection and refraction, the directrons maintain nearly constant structure and velocity, ensuring energy conservation and information integrity. Microscopic analyses of the director field and macroscopic evaluations of effective potential are employed to elucidate the dependence of reflection and refraction probabilities on the directron’s incident angle and the orientation difference across the interface. The findings provide valuable insights into the dynamics of solitary waves in structured liquid crystal systems, offering significant implications for the development of tunable photonic devices, reconfigurable optical systems, and nanoscale material engineering.

Science & Technology - Other Topics↗
DOE OSTI2022

Searching for dark radiation at the LHC

There is an interesting connection between early universe cosmology and searches for long-lived particles (LLPs) at the LHC. Light particles can be produced via freeze-in and act as dark radiation, contributing to the effective number of relativistic species $N_\text{eff}$. The parameter space of interest for future CMB missions points to LLP decay lengths in the mm to cm range. These decay lengths lie at the boundary between prompt and displaced signatures at the LHC and can be comprehensively explored only by combining searches for both. We consider a model where the LLP decays into a charged lepton and a (nearly) massless invisible particle. By reinterpreting searches for promptly decaying sleptons and for displaced leptons at both ATLAS and CMS we can then directly compare LHC exclusions with cosmological observables. Our results show how in this model the target value of CMB-S4 is already excluded by current LHC searches.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Boundary of the Distribution of Solar Wind Proton Beta versus Temperature Anisotropy

Abstract The frequency distribution of solar wind protons, measured in the vicinity of Earth’s orbit, is customarily plotted in ( β ∥ , T ⊥ / T ∥ ) phase space. Here, T ⊥ / T ∥ is the ratio of perpendicular and parallel temperatures, and β ∥ = 8 π nT ∥ / B 2 is the ratio of parallel thermal energy to background magnetic field energy, the so-called “parallel beta,” with ⊥ and ∥ denoting directions with respect to the ambient magnetic field. Such a frequency distribution, plotted as a two-dimensional histogram, forms a peculiar rhombic shape defined with an outer boundary in the said phase space. Past studies reveal that the threshold conditions for temperature anisotropy–driven plasma instability partially account for the boundary on the high- β ∥ side. The low- β ∥ side remains largely unexplained despite some efforts. Work by Vafin et al. recently showed that certain contours of collisional relaxation frequency, ν pp , when parameterized by T ⊥ / T ∥ and β ∥ , could match the overall shape of the left-hand boundary, thus suggesting that the collisional relaxation process might be closely related to the formation of the left-hand boundary. The present paper extends the analysis by Vafin et al. and carries out the dynamical computation of the collisional relaxation process for an ensemble of initial proton states with varying degrees of anisotropic temperatures. The final states of the relaxed protons are shown to closely match the observed boundary to the left of the ( β ∥ , T ⊥ / T ∥ ) phase space. When coupled with a similar set of calculations for the ensemble in the collective instability regime, it is found that the combined collisional/collective effects provide the baseline explanation for the observation.

Astronomy & Astrophysics↗
DOE OSTI2022

Micro-spectroscopic and freezing characterization of ice-nucleating particles collected in the marine boundary layer in the eastern North Atlantic

Abstract. Formation of atmospheric ice plays a crucial role in the microphysical evolution of mixed-phase and cirrus clouds and thus climate. How aerosol particles impact ice crystal formation by acting as ice-nucleating particles (INPs) is a subject of intense research activities. To improve understanding of atmospheric INPs, we examined daytime and nighttime particles collected during the Aerosol and Cloud Experiments in the Eastern North Atlantic (ACE-ENA) field campaign conducted in summer 2017. Collected particles, representative of a remote marine environment, were investigated for their propensity to serve as INPs in the immersion freezing (IMF) and deposition ice nucleation (DIN) modes. The particle population was characterized by chemical imaging techniques such as computer-controlled scanning electron microscopy with energy-dispersive X-ray analysis (CCSEM/EDX) and scanning transmission X-ray microscopy with near-edge X-ray absorption fine-structure spectroscopy (STXM/NEXAFS). Four major particle-type classes were identified where internally mixed inorganic–organic particles make up the majority of the analyzed particles. Following ice nucleation experiments, individual INPs were identified and characterized by SEM/EDX. The identified INP types belong to the major particle-type classes consisting of fresh sea salt with organics or processed sea salt containing dust and sulfur with organics. Ice nucleation experiments show IMF events at temperatures as low as 231 K, including the subsaturated regime. DIN events were observed at lower temperatures of 210 to 231 K. IMF and DIN observations were analyzed with regard to activated INP fraction, ice-nucleation active site (INAS) densities, and a water activity-based immersion freezing model (ABIFM) yielding heterogeneous ice nucleation rate coefficients. Observed IMF and DIN events of ice formation and corresponding derived freezing rates demonstrate that the marine boundary layer aerosol particles can serve as INPs under typical mixed-phase and cirrus cloud conditions. The derived IMF and DIN parameterizations allow for implementation in cloud and climate models to evaluate predictive effects of atmospheric ice crystal formation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

INVESTIGATION OF TURBULENT JETS, SPRAYS, AND SUPERCRITICAL MIXING BY TIME RESOLVED RAINBOW SCHLIEREN DEFLECTOMETRY

As technological advancements have enabled higher efficiency compression ignition (CI) engines (diesel engines) by operating at higher pressures and temperatures, the mechanisms of fuel-air mixing are postulated to change. Thus, there is a greater need for experimental data to develop and verify new phenomena. However, the fast time scales and harsh environments associated with CI engines present many experimental difficulties and requires high quality non-intrusive diagnostics to obtain useful data. The diagnostics itself must be validated prior to its implementation in an engine environment. The purpose of this research is to advance rainbow schlieren deflectometry (RSD) to study turbulent fuel sprays including the effects of supercritical mixing in CI engine environments. This objective is met by developing and applying quantitative RSD, for the first time, to three separate aspects of fuel sprays: 1) turbulent mixing, 2) phase boundaries, and 3) non-ideal gas mixing. First, quantitative RSD is demonstrated for turbulent mixing in a canonical helium jet. Rayleigh scattering data in the literature is used to validate the concentration measurements obtained by RSD. The results also provide explanations for phenomena occurring in the helium jet and provide data not yet present in the literature. Second, RSD is applied to high pressure multi-phase fuel sprays and new methodologies are developed to distinguish the liquid region from the vapor region and to quantify the in-between region on a probabilistic basis. Isolating liquid and vapor regions generally requires two separate diagnostics, e.g., Mie scattering (liquid region) and Rayleigh scattering (vapor region). However, the results of this study demonstrate that RSD can be used as a single diagnostic to measure both liquid and vapor boundaries. Third, the optical-to-thermodynamic relations associated with refrac-tometry are investigated under non-ideal gas mixing conditions to develop a generalized relationship valid for both ideal or non-ideal gas mixing. The last part of this research develops the theoretical tools to understand and experimentally realize supercritical mixing at diesel conditions. Guided by the theoretical analysis, RSD is applied to compare supercritical versus transcritical mixing at CI conditions using a constant pressure test rig. It is found that while fully supercritical mixing ers benefits in diesel engines, i.e, faster mixing times (50% higher velocities), the absence of droplets which cause high levels of particulate emissions, etc., it would be difficult to achieve in practice because of the requirements of fuels with low critical temperature and/or substantial fuel preheating.

02 PETROLEUM↗
DOE OSTI2021

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↗