Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Hubbard parameter”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Influence of high magnetic field on access to stationary H-modes and pedestal characteristics in Alcator C-Mod

We report recent Alcator C-Mod experiments have explored access to and characteristics of H-modes at magnetic fields approaching 8 T, the highest field achieved to date in a diverted tokamak. The H-modes originated from L-mode densities ranging from 1.1 x 10 20 m -3 to 2.8 x 10 20 m -3 , allowing insight into the density dependence of the H-mode power threshold at high magnetic field. This dependence is compared to predictions from the ITPA scaling law, finding that the law is approximately accurate at 7.8 T. However, the law underpredicted the high density H-mode threshold at lower magnetic field in previous C-Mod experiments, suggesting that the overall dependence of the threshold on magnetic field is weaker than predicted by the scaling law. The threshold data at 7.8 T also indicates that the onset of a low density branch at this magnetic field on C-Mod occurs below 1.4 x 10 20 m -3 , which is lower than predicted by an existing model for low density branch onset. The H-modes achieved steady-state densities ranging from 2.3 x 10 20 m -3 to 4.4 x 10 20 m -3 , and higher transient densities, and had values of q 95 from 3.3 to 6.0. This parameter range allowed the achievement of all three types of H-mode routinely observed at lower magnetic field on C-Mod: the stationary, ELM-suppressed Enhanced D α (EDA) regime, seen at high densities and high values of q 95 ; the nonstationary ELM-free regime, seen at lower densities and values of q 95 ; and the ELMy regime, seen at low density, moderate q 95 , and specialized plasma shape. The parameter space in which these regimes occur at 7.8 T is consistent with lower magnetic field experience. Pressure pedestal height at 7.8 T is compared to EPED predictions, and a scaling law for EDA density pedestal height developed between 4.5 T and 6.0 T is updated to include fields from 2.7 T to 7.8 T. Overall, this analysis increases confidence in the use of low magnetic field experience to predict some elements of high magnetic field tokamak behavior.

7.8 T↗
DOE OSTI2016

FERROMAGNETISM AND PARTICLE COLLISIONS: APPLICATIONS TO PROTOPLANETARY DISKS AND THE METEORITICAL RECORD

The meteoritical record shows both iron partitioning and tungsten isotopic partitioning between the matrix and chondrules. Tungsten is not abundant enough to have driven its own isotopic partitioning, but if tungsten were correlated with iron then ferromagnetic interaction grains could help explain both observations. We derive a practical parameterization for the increase in particle–particle collision rates caused by mutually attracting particle magnetic dipole moments. While the appropriate magnetic parameters remain uncertain, we show that ambient magnetic fields in protoplanetary disks are expected to be strong enough to magnetize iron metal bearing dust grains sufficiently to drive large increases in their collision rates. Such increased collision rates between iron-metal-rich grains could help preserve primordial iron and W isotopic inhomogeneities, and would help explain why the meteoritical record shows their partitioning in the solar nebula. The importance of magnetic interactions for larger grains whose growth is balanced by fragmentation is less clear and will require future laboratory or numerical studies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

MINERAL PROCESSING BY SHORT CIRCUITS IN PROTOPLANETARY DISKS

Meteoritic chondrules were formed in the early solar system by brief heating of silicate dust to melting temperatures. Some highly refractory grains (Type B calcium-aluminum-rich inclusions, CAIs) also show signs of transient heating. A similar process may occur in other protoplanetary disks, as evidenced by observations of spectra characteristic of crystalline silicates. One possible environment for this process is the turbulent magnetohydrodynamic flow thought to drive accretion in these disks. Such flows generally form thin current sheets, which are sites of magnetic reconnection, and dissipate the magnetic fields amplified by a disk dynamo. We suggest that it is possible to heat precursor grains for chondrules and other high-temperature minerals in current sheets that have been concentrated by our recently described short-circuit instability. We extend our work on this process by including the effects of radiative cooling, taking into account the temperature dependence of the opacity; and by examining current sheet geometry in three-dimensional, global models of magnetorotational instability. We find that temperatures above 1600 K can be reached for favorable parameters that match the ideal global models. This mechanism could provide an efficient means of tapping the gravitational potential energy of the protoplanetary disk to heat grains strongly enough to form high-temperature minerals. The volume-filling nature of turbulent magnetic reconnection is compatible with constraints from chondrule-matrix complementarity, chondrule-chondrule complementarity, the occurrence of igneous rims, and compound chondrules. The same short-circuit mechanism may perform other high-temperature mineral processing in protoplanetary disks such as the production of crystalline silicates and CAIs.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Influence of High Magnetic Field on Access to Stationary H-modes and Pedestal Characteristics in Alcator C-Mod

Recent Alcator C-Mod experiments have explored access to and characteristics of H-modes at magnetic fields approaching 8 T, the highest field achieved to date in a diverted tokamak. The H-modes originated from L-mode densities ranging from $1.1 \times 10^{20} m^{-3}$ to $2.8 \times 10^{20} m^{-3}$, allowing insight into the density dependence of the H-mode power threshold at high magnetic field. This dependence is compared to predictions from the ITPA scaling law (\cite{martin2008power}), finding that the law is approximately accurate at 7.8 T. However, the law underpredicted the high density H-mode threshold at lower magnetic field in previous C-Mod experiments (\cite{ma2012scaling}), \hl{suggesting that the overall dependence of the threshold on magnetic field is weaker than predicted by the scaling law.} The threshold data at 7.8 T also indicates that the onset of a low density branch at this magnetic field on C-Mod occurs below $1.4 \times 10^{20} m^{-3}$, which is lower than predicted by an existing model for low density branch onset. The H-modes achieved steady-state densities ranging from $2.3 \times 10^{20} m^{-3}$ to $4.4 \times 10^{20} m^{-3}$, and higher transient densities, and had values of $q_{95}$ from 3.3 to 6.0. This parameter range allowed the achievement of all three types of H-mode routinely observed at lower magnetic field on C-Mod: the stationary, ELM-suppressed enhanced D-alpha (EDA) regime, seen at high densities and high values of $q_{95}$; the nonstationary ELM-free regime, seen at lower densities and values of $q_{95}$; and the ELMy regime, seen at low density, moderate $q_{95}$, and specialized plasma shape. The parameter space in which these regimes occur at 7.8 T is consistent with lower magnetic field experience. Pressure pedestal height at 7.8 T is compared to EPED \cite{snyder2009development, snyder2011first} predictions, and a scaling law for EDA density pedestal height developed between 4.5 and 6.0 T is updated to include fields from 2.7 T to 7.8 T. Overall, this analysis increases confidence in the use of low magnetic field experience to predict some elements of high magnetic field tokamak behavior.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2018

Fuzzy Modeling and Parallel Distributed Compensation for Aircraft Flight Control from Simulated Flight Data

A method is described that combines fuzzy system identification techniques with Parallel Distributed Compensation (PDC) to develop nonlinear control methods for aircraft using minimal a priori knowledge, as part of NASA’s Learn-to-Fly initiative. A fuzzy model was generated with simulated flight data, and consisted of a weighted average of multiple linear time invariant state-space cells having parameters estimated using the equation-error approach and a least-squares estimator. A compensator was designed for each subsystem using Linear Matrix Inequalities (LMI) to guarantee closed-loop stability and performance requirements. This approach is demonstrated using simulated flight data to automatically develop a fuzzy model and design control laws for a simplified longitudinal approximation of the F-16 nonlinear flight dynamics simulation. Results include a comparison of flight data with the estimated fuzzy models and simulations that illustrate the feasibility and utility of the combined fuzzy modeling and control approach.

Weinstein, Rose↗
NASA NTRS1977

The occultation of Epsilon Geminorum by Mars - Analysis of McDonald data

An analysis of observations of the occultation of Epsilon Gem by Mars on April 8, 1976, is presented. The data were obtained by three neighboring telescopes at McDonald Observatory. Intensity fluctuations on time scales of the order of 100 ms were observed simultaneously at the three telescopes. As the observations compare well with predictions of turbulent scintillation theory, it is concluded that such fluctuations were probably largely the effect of stellar scintillations in the Martian atmosphere. The stellar diameter is included as a parameter in the theory but in a way which differs from previously published interpretations of occultations of extended sources by planetary atmospheres. Scintillations govern the experimental uncertainty in the deduction of the scale height of the high Martian atmosphere. A density scale height of 9.9 + or - 2.5 km is obtained at an altitude of 74 + or - 8 km above the mean surface. For CO 2 gas, this result corresponds to a temperature of 190 + or - 50 K.

Africano, J.↗
DOE OSTI2017

Physics and performance of the I-mode regime over an expanded operating space on Alcator C-Mod

New results on the I-mode regime of operation on the Alcator C-Mod tokamak are reported. This ELM-free regime features high energy confinement and a steep temperature pedestal, while particle confinement remains at L-mode levels, giving stationary density and avoiding impurity accumulation. I-mode has now been obtained over nearly all of the magnetic fields and currents possible in this high field tokamak (I p 0.55–1.7 MA, B T 2.8–8 T) using a configuration with B × ∇B drift away from the X-point. Results at 8 T confirm that the L–I power threshold varies only weakly with B T , and that the power range for I-mode increases with B T ; no 8 T discharges transitioned to H-mode. Parameter dependences of energy confinement are investigated. Core transport simulations are giving insight into the observed turbulence reduction, profile stiffness and confinement improvement. Pedestal models explain the observed stability to ELMs, and can simulate the observed weakly coherent mode. Conditions for I–H transitions have complex dependences on density as well as power. I-modes have now been maintained in near-DN configurations, leading to improved divertor power flux sharing. Overall, prospects for I-mode on future fusion devices such as ITER and ARC are encouraging. Further experiments on other tokamaks are needed to improve confidence in extrapolation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

The emergence of hydrogeophysics for improved understanding of subsurface processes over multiple scales

Geophysics provides a multidimensional suite of investigative methods that are transforming our ability to see into the very fabric of the subsurface environment, and monitor the dynamics of its fluids and the biogeochemical reactions that occur within it. Here we document how geophysical methods have emerged as valuable tools for investigating shallow subsurface processes over the past two decades and offer a vision for future developments relevant to hydrology and also ecosystem science. The field of “hydrogeophysics” arose in the late 1990s, prompted, in part, by the wealth of studies on stochastic subsurface hydrology that argued for better field-based investigative techniques. These new hydrogeophysical approaches benefited from the emergence of practical and robust data inversion techniques, in many cases with a view to quantify shallow subsurface heterogeneity and the associated dynamics of subsurface fluids. Furthermore, the need for quantitative characterization stimulated a wealth of new investigations into petrophysical relationships that link hydrologically relevant properties to measurable geophysical parameters. Development of time-lapse approaches provided a new suite of tools for hydrological investigation, enhanced further with the realization that some geophysical properties may be sensitive to biogeochemical transformations in the subsurface environment, thus opening up the new field of “biogeophysics.” Early hydrogeophysical studies often concentrated on relatively small “plot-scale” experiments. More recently, however, the translation to larger-scale characterization has been the focus of a number of studies. In conclusion, geophysical technologies continue to develop, driven, in part, by the increasing need to understand and quantify key processes controlling sustainable water resources and ecosystem services.

58 GEOSCIENCES↗
DOE OSTI2018

Surface heat flux feedback controlled impurity seeding experiments with Alcator C-Mod’s high-Z vertical target plate divertor: performance, limitations and implications for fusion power reactors

The Alcator C-Mod team has recently developed a feedback system to measure and control surface heat flux in real-time. The system uses real-time measurements of surface heat flux from surface thermocouples and a pulse-width modulated piezo valve to inject low-Z impurities (typically N2) into the private flux region. It has been used in C-Mod to mitigate peak surface heat fluxes >40 MW/m2 down to <10 MW/m2 while maintaining excellent core confinement, H98>1. While the system works quite well under relatively steady conditions, use of it during transients has revealed important limitations on feedback control of impurity seeding in conventional vertical target plate divertors. In some cases, the system is unable to avoid plasma reattachment to the divertor plate or the formation of a confinement-damaging x-point MARFE. This is due to the small operational window for mitigated heat flux in the parameters of incident plasma heat flux, plasma density, and impurity density as well as the relatively slow response of the impurity gas injection system compared to plasma transients. Given the severe consequences for failure of such a system to operate reliably in a reactor, there is substantial risk that the conventional vertical target plate divertor will not provide an adequately controllable system in reactor-class devices. These considerations motivate the need to develop passively stable, highly compliant divertor configurations and experimental facilities that can test such possible solutions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Surface heat flux feedback controlled impurity seeding experiments with Alcator C-Mod's high-Z vertical target plate divertor: performance, limitations and implications for fusion power reactors

The Alcator C-Mod team has recently developed a feedback system to measure and control surface heat flux in real-time. The system uses real-time measurements of surface heat flux from surface thermocouples and a pulse-width modulated piezo valve to inject low-Z impurities (typically N 2 ) into the private flux region. It has been used in C-Mod to mitigate peak surface heat fluxes >40 MW m -2 down to <10 MW m -2 while maintaining excellent core confinement, H 98 > 1. While the system works quite well under relatively steady conditions, use of it during transients has revealed important limitations on feedback control of impurity seeding in conventional vertical target plate divertors. In some cases, the system is unable to avoid plasma reattachment to the divertor plate or the formation of a confinement-damaging x-point MARFE. This is due to the small operational window for mitigated heat flux in the parameters of incident plasma heat flux, plasma density, and impurity density as well as the relatively slow response of the impurity gas injection system compared to plasma transients. Given the severe consequences for failure of such a system to operate reliably in a reactor, there is substantial risk that the conventional vertical target plate divertor will not provide an adequately controllable system in reactor-class devices. Finally, these considerations motivate the need to develop passively stable, highly compliant divertor configurations and experimental facilities that can test such possible solutions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Determination of confinement regime boundaries via separatrix parameters on Alcator C-Mod based on a model for interchange-drift-Alfvén turbulence

The separatrix operational space (SepOS) model (Eich and Manz 2021 Nucl. Fusion 61 086017) is shown to predict the L–H transition, the L-mode density limit, and the ideal magnetohydrodynamic ballooning limit in terms of separatrix parameters for a wide range of Alcator C-Mod plasmas. The model is tested using Thomson scattering measurements across a wide range of operating conditions on C-Mod, spanning $\overline{n}_{\mathrm{e}}$ = 0.3–5.5 $\times 10^{20}\,$m −3 , $B_{\mathrm{t}} = 2.5$–8.0 T, and $B_{\mathrm{p}}$= 0.1–1.2 T. An empirical regression for the electron pressure gradient scale length, $\lambda_{{p}_{\mathrm{e}}}$, against a turbulence control parameter, $\alpha_{\mathrm{t}}$, and the poloidal fluid gyroradius, $\rho_{\mathrm{s,p}}$, is constructed for H-modes and found to require positive exponents for both regression parameters, indicating turbulence widening of near-scrape-off layer widths at high $\alpha_{\mathrm{t}}$ and an inverse scaling with $B_{\mathrm{p}}$, consistent with results on ASDEX Upgrade. The SepOS model is also tested in the unfavorable drift direction and found to apply well to all three boundaries, including the L–H transition as long as a correction to the Reynolds energy transfer term, $\alpha_\mathrm{RS} \lt 1$ is applied. I-modes typically exist in the unfavorable drift direction for values of $\alpha_{\mathrm{t}} \lesssim 0.3$. Finally, an experiment studying the transition between the Type-I ELMy and EDA H-mode is analyzed using the same framework. It is found that a recently identified boundary $\alpha_{\mathrm{t}} = 0.55$ at excludes most EDA H-modes but that the balance of wavenumbers responsible for the L-mode density limit, namely $k_\mathrm{EM} = k_\mathrm{RBM}$, may better describe the transition on C-Mod. The ensemble of boundaries validated and explored is then applied to project regime access and limit avoidance for the SPARC primary reference discharge parameters.

ELM suppression↗
DOE OSTI2015

Many-body localization in disorder-free systems: The importance of finite-size constraints

Recently it has been suggested that many-body localization (MBL) can occur in translation-invariant systems, and candidate 1D models have been proposed. We find that such models, in contrast to MBL systems with quenched disorder, typically exhibit much more severe finite-size effects due to the presence of two or more vastly different energy scales. In a finite system, this can artificially split the density of states (DOS) into bands separated by large gaps. We argue for such models to faithfully represent the thermodynamic limit behavior, the ratio of relevant coupling must exceed a certain system-size depedent cutoff, chosen such that various bands in the DOS overlap one another. Setting the parameters this way to minimize finite-size effects, we study several translation-invariant MBL candidate models using exact diagonalization. Based on diagnostics including entanglement and local observables, we observe thermal (ergodic), rather than MBL-like behavior. Our results suggest that MBL in translation-invariant systems with two or more very different energy scales is less robust than perturbative arguments suggest, possibly pointing to the importance of non-perturbative effects which induce delocalization in the thermodynamic limit.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

Remote sensing of freeze-thaw transitions in Arctic soils using the complex resistivity method

Our ability to monitor freeze - thaw transitions is critical to developing a predictive understanding of biogeochemical transitions and carbon dynamics in high latitude environments. In this study, we conducted laboratory column experiments to explore the potential of the complex resistivity method for monitoring the freeze - thaw transitions of the arctic permafrost soils. Samples for the experiment were collected from the upper active layer of Gelisol soils at the Barrow Environmental Observatory, Barrow Alaska. Freeze - thaw transitions were induced through exposing the soil column to controlled temperature environments at 4 C and -20 C. Complex resistivity and temperature measurements were collected regularly during the freeze - thaw transitions using electrodes and temperature sensors installed along the column. During the experiments, over two orders of magnitude of resistivity variations were observed when the temperature was increased or decreased between -20 C and 0 C. Smaller resistivity variations were also observed during the isothermal thawing or freezing processes that occurred near 0 C. Single frequency electrical phase response and imaginary conductivity at 1 Hz were found to be exclusively related to the unfrozen water in the soil matrix, suggesting that these geophysical 24 attributes can be used as a proxy for the monitoring of the onset and progression of the freeze - thaw transitions. Spectral electrical responses and fitted Cole Cole parameters contained additional information about the freeze - thaw transition affected by the soil grain size distribution. Specifically, a shift of the observed spectral response to lower frequency was observed during isothermal thawing process, which we interpret to be due to sequential thawing, first from fine then to coarse particles within the soil matrix. Our study demonstrates the potential of the complex resistivity method for remote monitoring of freeze - thaw transitions in arctic soils. Although conducted at the laboratory scale, this study provides the foundation for exploring the potential of the complex resistivity signals for monitoring spatiotemporal variations of freeze - thaw transitions over field-relevant scales.

climate↗
DOE OSTI2021

Dynamical vertex approximation for many-electron systems with spontaneously broken SU(2) symmetry

We generalize the formalism of the dynamical vertex approximation (DΓA)—a diagrammatic extension of the dynamical mean-field theory (DMFT)—to treat magnetically ordered phases. To this aim, we start by concisely illustrating the many-electron formalism for performing ladder resummations of Feynman diagrams in systems with broken SU(2) symmetry associated to ferromagnetic (FM) or antiferromagnetic (AF) order. We then analyze the algorithmic simplifications introduced by taking the local approximation of the two-particle irreducible vertex functions in the Bethe-Salpeter equations, which defines the ladder implementation of DΓA for magnetic systems. The relation of this assumption with the DMFT limit of large coordination-number/high dimensions is explicitly discussed. As a last step, we derive the expression for the ladder DΓA self-energy in the FM- and AF-ordered phases of the Hubbard model. The physics emerging in the AF-ordered case is explicitly illustrated by means of approximated calculations based on a static mean-field input for DΓA equations. The results obtained capture fundamental aspects of both metallic and insulating ground states of two-dimensional antiferromagnets, providing a reliable compass for future, more extensive applications of our approach. Furthermore, possible routes to further develop diagrammatic-based treatments of magnetic phases in correlated electron systems are briefly outlined in the Conclusions.

2-dimensional systems↗
DOE OSTI2011

Fibrous filter efficiency and pressure drop in the viscous-inertial transition flow regime.

Fibrous filter pressure drop and aerosol collection efficiency were measured at low air pressures (0.2 to 0.8 atm) and high face velocities (5 to 20 meters per second) to give fiber Reynolds numbers in the viscous-inertial transition flow regime (1 to 16). In this regime, contemporary filtration theory based on Kuwabara's viscous flow through an ensemble of fibers under-predicts single fiber impaction by several orders of magnitude. Streamline curvature increases substantially as inertial forces become dominant. Dimensionless pressure drop measurements followed the viscous-inertial theory of Robinson and Franklin rather than Darcy's linear pressure-velocity relationship (1972). Sodium chloride and iron nano-agglomerate test aerosols were used to evaluate the effects of particle density and shape factor. Total filter efficiency collapsed when plotted against the particle Stokes and fiber Reynolds numbers. Efficiencies were then fitted with an impactor type equation where the cutpoint Stokes number and a steepness parameter described data well in the sharply increasing portion of the curve (20% to 80% efficiency). The cutpoint Stokes number was a linearly decreasing function of fiber Reynolds number. Single fiber efficiencies were calculated from total filter efficiencies and compared to contemporary viscous flow impaction theory (Stechkina et al. 1969), and numerical simulations from the literature. Existing theories under-predicted measured single fiber efficiencies although the assumption of uniform flow conditions for each successive layer of fibers is questionable; the common exponential relationship between single fiber efficiency and total filter efficiency may not be appropriate in this regime.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1983

A comparison with theory of peak to peak sound level for a model helicopter rotor generating blade slap at low tip speeds

Mini-tuft and smoke flow visualization techniques have been developed for the investigation of model helicopter rotor blade vortex interaction noise at low tip speeds. These techniques allow the parameters required for calculation of the blade vortex interaction noise using the Widnall/Wolf model to be determined. The measured acoustics are compared with the predicted acoustics for each test condition. Under the conditions tested it is determined that the dominating acoustic pulse results from the interaction of the blade with a vortex 1-1/4 revolutions old at an interaction angle of less than 8 deg. The Widnall/Wolf model predicts the peak sound pressure level within 3 dB for blade vortex separation distances greater than 1 semichord, but it generally over predicts the peak S.P.L. by over 10 dB for blade vortex separation distances of less than 1/4 semichord.

Fontana, R. R.↗
DOE OSTI2013

HUBBLE PARAMETER MEASUREMENT CONSTRAINTS ON THE COSMOLOGICAL DECELERATION-ACCELERATION TRANSITION REDSHIFT

We compile a list of 28 independent measurements of the Hubble parameter between redshifts 0.07 {<=} z {<=} 2.3 and use this to place constraints on model parameters of constant and time-evolving dark energy cosmologies. These H(z) measurements by themselves require a currently accelerating cosmological expansion at about, or better than, 3{sigma} confidence. The mean and standard deviation of the six best-fit model deceleration-acceleration transition redshifts (for the three cosmological models and two Hubble constant priors we consider) are z{sub da} = 0.74 {+-} 0.05, in good agreement with the recent Busca et al. determination of z{sub da} = 0.82 {+-} 0.08 based on 11 H(z) measurements between redshifts 0.2 {<=} z {<=} 2.3, almost entirely from baryon-acoustic-oscillation-like data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Using ASCEM modeling and visualization to inform stakeholders of contaminant plume evolution and remediation efficacy at F-Basin Savannah river, SC - 15156

Communication with stakeholders, regulatory agencies, and the public is an essential part of implementing different remediation and monitoring activities, and developing site closure strategies at contaminated sites. Modeling of contaminant plume evolution plays a critical role in estimating the benefit, cost, and risk of particular options. At the same time, effective visualization of monitoring data and modeling results are particularly important for conveying the significance of the results and observations. In this paper, we present the results of the Advanced Simulation Capability for Environmental Management (ASCEM) project, including the discussion of the capabilities of newly developed ASCEM software package, along the its application to the F-Area Seepage Basins located in the U.S. Department of Energy Savannah River Site (SRS). ASCEM software includes state-of-the-art numerical methods for simulating complex flow and reactive transport, as well as various tool-sets such as a graphical user interface (GUI), visualization, data management, uncertainty quantification, and parameter estimation. Using this software, we have developed an advanced visualization of tritium plume migration coupled with a data management system, and simulated a three-dimensional model of flow and plume evolution on a high-performance computing platform. We evaluated the effect of engineered flow barriers on a nonreactive tritium plume, through advanced plume visualization and modeling of tritium plume migration. In addition, we developed a geochemical reaction network to describe complex geochemical processes at the site, and evaluated the impact of coupled hydrological and geochemical heterogeneity. These results are expected to support SRS's monitoring activities and operational decisions. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗