National Library of Energy BETA

Sample records for reactive transport models

  1. Three-Dimensional Modeling of the Reactive Transport of CO2 and...

    Office of Scientific and Technical Information (OSTI)

    ... Subject: CO2 reservoir, reactive transport, geochemistry, geomechanics, geomechanical modeling, mineral composition, homogenization, fault, hydraulic fracture, elastic modulus Word ...

  2. Uranium transport in a crushed granodiorite: Experiments and reactive transport modeling

    DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)

    Dittrich, T. M.; Reimus, P. W.

    2015-02-12

    The primary objective of this study was to develop and demonstrate an experimental method to refine and better parameterize process models for reactive contaminant transport in aqueous subsurface environments and to reduce conservatism in such models without attempting to fully describe the geochemical system.

  3. Maximum likelihood Bayesian model averaging and its predictive analysis for groundwater reactive transport models

    DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)

    Lu, Dan; Ye, Ming; Curtis, Gary P.

    2015-08-01

    While Bayesian model averaging (BMA) has been widely used in groundwater modeling, it is infrequently applied to groundwater reactive transport modeling because of multiple sources of uncertainty in the coupled hydrogeochemical processes and because of the long execution time of each model run. To resolve these problems, this study analyzed different levels of uncertainty in a hierarchical way, and used the maximum likelihood version of BMA, i.e., MLBMA, to improve the computational efficiency. Our study demonstrates the applicability of MLBMA to groundwater reactive transport modeling in a synthetic case in which twenty-seven reactive transport models were designed to predict themore » reactive transport of hexavalent uranium (U(VI)) based on observations at a former uranium mill site near Naturita, CO. Moreover, these reactive transport models contain three uncertain model components, i.e., parameterization of hydraulic conductivity, configuration of model boundary, and surface complexation reactions that simulate U(VI) adsorption. These uncertain model components were aggregated into the alternative models by integrating a hierarchical structure into MLBMA. The modeling results of the individual models and MLBMA were analyzed to investigate their predictive performance. The predictive logscore results show that MLBMA generally outperforms the best model, suggesting that using MLBMA is a sound strategy to achieve more robust model predictions relative to a single model. MLBMA works best when the alternative models are structurally distinct and have diverse model predictions. When correlation in model structure exists, two strategies were used to improve predictive performance by retaining structurally distinct models or assigning smaller prior model probabilities to correlated models. Since the synthetic models were designed using data from the Naturita site, the results of this study are expected to provide guidance for real-world modeling. Finally

  4. Maximum likelihood Bayesian model averaging and its predictive analysis for groundwater reactive transport models

    SciTech Connect (OSTI)

    Lu, Dan; Ye, Ming; Curtis, Gary P.

    2015-08-01

    While Bayesian model averaging (BMA) has been widely used in groundwater modeling, it is infrequently applied to groundwater reactive transport modeling because of multiple sources of uncertainty in the coupled hydrogeochemical processes and because of the long execution time of each model run. To resolve these problems, this study analyzed different levels of uncertainty in a hierarchical way, and used the maximum likelihood version of BMA, i.e., MLBMA, to improve the computational efficiency. Our study demonstrates the applicability of MLBMA to groundwater reactive transport modeling in a synthetic case in which twenty-seven reactive transport models were designed to predict the reactive transport of hexavalent uranium (U(VI)) based on observations at a former uranium mill site near Naturita, CO. Moreover, these reactive transport models contain three uncertain model components, i.e., parameterization of hydraulic conductivity, configuration of model boundary, and surface complexation reactions that simulate U(VI) adsorption. These uncertain model components were aggregated into the alternative models by integrating a hierarchical structure into MLBMA. The modeling results of the individual models and MLBMA were analyzed to investigate their predictive performance. The predictive logscore results show that MLBMA generally outperforms the best model, suggesting that using MLBMA is a sound strategy to achieve more robust model predictions relative to a single model. MLBMA works best when the alternative models are structurally distinct and have diverse model predictions. When correlation in model structure exists, two strategies were used to improve predictive performance by retaining structurally distinct models or assigning smaller prior model probabilities to correlated models. Since the synthetic models were designed using data from the Naturita site, the results of this study are expected to provide guidance for real-world modeling. Finally, limitations of

  5. Pore scale modeling of reactive transport involved in geologic CO2 sequestration

    SciTech Connect (OSTI)

    Kang, Qinjin; Lichtner, Peter C; Viswanathan, Hari S; Abdel-fattah, Amr I

    2009-01-01

    We apply a multi-component reactive transport lattice Boltzmann model developed in previolls studies to modeling the injection of a C02 saturated brine into various porous media structures at temperature T=25 and 80 C. The porous media are originally consisted of calcite. A chemical system consisting of Na+, Ca2+, Mg2+, H+, CO2(aq), and CI-is considered. The fluid flow, advection and diHusion of aqueous species, homogeneous reactions occurring in the bulk fluid, as weB as the dissolution of calcite and precipitation of dolomite are simulated at the pore scale. The effects of porous media structure on reactive transport are investigated. The results are compared with continuum scale modeling and the agreement and discrepancy are discussed. This work may shed some light on the fundamental physics occurring at the pore scale for reactive transport involved in geologic C02 sequestration.

  6. Reaction-based reactive transport modeling of Fe(III)

    SciTech Connect (OSTI)

    Kemner, K.M.; Kelly, S.D.; Burgos, Bill; Roden, Eric

    2006-06-01

    This research project (started Fall 2004) was funded by a grant to Argonne National Laboratory, The Pennsylvania State University, and The University of Alabama in the Integrative Studies Element of the NABIR Program (DE-FG04-ER63914/63915/63196). Dr. Eric Roden, formerly at The University of Alabama, is now at the University of Wisconsin, Madison. Our project focuses on the development of a mechanistic understanding and quantitative models of coupled Fe(III)/U(VI) reduction in FRC Area 2 sediments. This work builds on our previous studies of microbial Fe(III) and U(VI) reduction, and is directly aligned with the Scheibe et al. NABIR FRC Field Project at Area 2.

  7. Coupled reactive mass transport and fluid flow: Issues in model verification

    SciTech Connect (OSTI)

    Freedman, Vicky L.; Ibaraki, Motomu

    2003-01-03

    Model verification and validation are both important steps in the development of reactive transport models. In this paper, a distinction is made between verification and validation, and the focus is on codifying the issues of verification for a numerical, reactive transport flow model. First, the conceptual basis of model verification is reviewed, which shows that verification should be understood as a first step in model development, and be followed by a protocol that assures that the model accurately represents system behavior. Second, commonly used procedures and methods of model verification are presented. In the third part of this paper, an intercomparison of models is used to demonstrate that model verification can be performed despite differences in hydrogeochemical transport code formulations. Results of an example simulation of transport are presented in which the numerical model is tested against other hydrogeochemical codes. Different kinetic formulations between solid and aqueous phases used among numerical models complicates model verification. This test problem involves uranium transport under conditions of varying pH and oxidation potential, with reversible precipitation of calcium uranate and coffinite. Results between the different hydrogeochemical transport codes show differences in oxidation potentials, but similarities in mineral assemblages and aqueous transport patterns. Because model verification can be further complicated by differences in the approach for solving redox problems, a comparison of a fugacity approach to both the external approach (based on hypothetical electron activity) and effective internal approach (based on conservation of electrons) is performed. The comparison demonstrates that the oxygen fugacity approach produces different redox potentials and mineral assemblages than both the effective internal and external approaches.

  8. Application of a data assimilation method via an ensemble Kalman filter to reactive urea hydrolysis transport modeling

    SciTech Connect (OSTI)

    Juxiu Tong; Bill X. Hu; Hai Huang; Luanjin Guo; Jinzhong Yang

    2014-03-01

    With growing importance of water resources in the world, remediations of anthropogenic contaminations due to reactive solute transport become even more important. A good understanding of reactive rate parameters such as kinetic parameters is the key to accurately predicting reactive solute transport processes and designing corresponding remediation schemes. For modeling reactive solute transport, it is very difficult to estimate chemical reaction rate parameters due to complex processes of chemical reactions and limited available data. To find a method to get the reactive rate parameters for the reactive urea hydrolysis transport modeling and obtain more accurate prediction for the chemical concentrations, we developed a data assimilation method based on an ensemble Kalman filter (EnKF) method to calibrate reactive rate parameters for modeling urea hydrolysis transport in a synthetic one-dimensional column at laboratory scale and to update modeling prediction. We applied a constrained EnKF method to pose constraints to the updated reactive rate parameters and the predicted solute concentrations based on their physical meanings after the data assimilation calibration. From the study results we concluded that we could efficiently improve the chemical reactive rate parameters with the data assimilation method via the EnKF, and at the same time we could improve solute concentration prediction. The more data we assimilated, the more accurate the reactive rate parameters and concentration prediction. The filter divergence problem was also solved in this study.

  9. Chombo-Crunch: Modeling Pore-Scale Reactive Transport in Carbon

    Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)

    Sequestration | Argonne Leadership Computing Facility Transient pH on calcite grains in capillary tube experiment. 1 billion grid points computed on 48K cores. 1 micron resolution. Transient pH on calcite grains in capillary tube experiment. 1 billion grid points computed on 48K cores. 1 micron resolution. David Trebotich, Lawrence Berkeley National Laboratory Chombo-Crunch: Modeling Pore-Scale Reactive Transport in Carbon Sequestration PI Name: David Trebotich PI Email: treb@lbl.gov

  10. A Dual Regime Reactive Transport Model for Simulation of High Level Waste Tank Closure Scenarios - 13375

    SciTech Connect (OSTI)

    Sarkar, Sohini; Kosson, David S.; Brown, Kevin; Garrabrants, Andrew C.; Meeussen, Hans; Van der Sloot, Hans

    2013-07-01

    A numerical simulation framework is presented in this paper for estimating evolution of pH and release of major species from grout within high-level waste tanks after closure. This model was developed as part of the Cementitious Barriers Partnership. The reactive transport model consists of two parts - (1) transport of species, and (2) chemical reactions. The closure grout can be assumed to have varying extents of cracking and composition for performance assessment purposes. The partially or completely degraded grouted tank is idealized as a dual regime system comprising of a mobile region having solid materials with cracks and macro-pores, and an immobile/stagnant region having solid matrix with micropores. The transport profiles of the species are calculated by incorporating advection of species through the mobile region, diffusion of species through the immobile/stagnant region, and exchange of species between the mobile and immobile regions. A geochemical speciation code in conjunction with the pH dependent test data for a grout material is used to obtain a mineral set that best describes the trends in the test data of the major species. The dual regime reactive transport model predictions are compared with the release data from an up-flow column percolation test. The coupled model is then used to assess effects of crack state of the structure, rate and composition of the infiltrating water on the pH evolution at the grout-waste interface. The coupled reactive transport model developed in this work can be used as part of the performance assessment process for evaluating potential risks from leaching of a cracked tank containing elements of human health and environmental concern. (authors)

  11. A reaction-based paradigm to model reactive chemical transport in groundwater with general kinetic and equilibrium reactions

    SciTech Connect (OSTI)

    Zhang, Fan; Yeh, Gour-Tsyh; Parker, Jack C; Brooks, Scott C; Pace, Molly; Kim, Young Jin; Jardine, Philip M; Watson, David B

    2007-01-01

    This paper presents a reaction-based water quality transport model in subsurface flow systems. Transport of chemical species with a variety of chemical and physical processes is mathematically described by M partial differential equations (PDEs). Decomposition via Gauss-Jordan column reduction of the reaction network transforms M species reactive transport equations into two sets of equations: a set of thermodynamic equilibrium equations representing NE equilibrium reactions and a set of reactive transport equations of M-NE kinetic-variables involving no equilibrium reactions (a kinetic-variable is a linear combination of species). The elimination of equilibrium reactions from reactive transport equations allows robust and efficient numerical integration. The model solves the PDEs of kinetic-variables rather than individual chemical species, which reduces the number of reactive transport equations and simplifies the reaction terms in the equations. A variety of numerical methods are investigated for solving the coupled transport and reaction equations. Simulation comparisons with exact solutions were performed to verify numerical accuracy and assess the effectiveness of various numerical strategies to deal with different application circumstances. Two validation examples involving simulations of uranium transport in soil columns are presented to evaluate the ability of the model to simulate reactive transport with complex reaction networks involving both kinetic and equilibrium reactions.

  12. Coupling a Genome-Scale Metabolic Model with a Reactive Transport Model to Describe In Situ Uranium Bioremediation

    SciTech Connect (OSTI)

    Scheibe, Timothy D.; Mahadevan, Radhakrishnan; Fang, Yilin; Garg, Srinath; Long, Philip E.; Lovley, Derek R.

    2009-03-01

    Quantitative numerical simulation codes known as reactive transport models are widely used for simulating the hydrologic transport and geochemical speciation of dissolved constituents in the subsurface (Steefel et al., 2005). Because the activity of microorganisms strongly influences the fate of many constituents, both organic and inorganic, such models often include microbially-mediated reactions in their reaction networks (Hunter et al., 1998; Burgos et al., 2002; Fang et al., 2006; Scheibe et al., 2006; Yabusaki et al., 2007). However, the canonical form and stoichiometry of microbial reactions, reaction rate formulations and parameters, and biomass growth yield coefficients are prescribed a priori and applied over the entire range of simulated conditions. This approach does not account for the fact that fundamental microbial functions vary in response to local variations in environmental conditions(Stewart and Franklin, 2008). Multiple alternative reaction pathways are encoded in microbial genomes; specific pathways become active or inactive in response to, for example, nutrient limitation. Recent advances in genomic analysis allow us to define cellular metabolic networks, and accurate predictions of active pathways and reaction fluxes have been made using constraint-based metabolic models (Mahadevan et al., 2002; Price et al., 2003; Reed and Palsson, 2003; Mahadevan et al., 2006). Here, we demonstrate for the first time a methodology of coupling constraint-based metabolic models with reactive transport models. Our approach integrates advanced microbiological characterization, hydrology, and geochemistry in a powerful manner that will significantly improve subsurface reactive transport models.

  13. Multiphase flow and multicomponent reactive transport model of the ventilation experiment in Opalinus clay

    SciTech Connect (OSTI)

    Zheng, L.; Samper, J.; Montenegro, L.; Major, J.C.

    2008-10-15

    During the construction and operational phases of a high-level radioactive waste (HLW) repository constructed in a clay formation, ventilation of underground drifts will cause desaturation and oxidation of the rock. The Ventilation Experiment (VE) was performed in a 1.3 m diameter unlined horizontal microtunnel on Opalinus clay at Mont Terri underground research laboratory in Switzerland to evaluate the impact of desaturation on rock properties. A multiphase flow and reactive transport model of VE is presented here. The model accounts for liquid, vapor and air flow, evaporation/condensation and multicomponent reactive solute transport with kinetic dissolution of pyrite and siderite and local-equilibrium dissolution/precipitation of calcite, ferrihydrite, dolomite, gypsum and quartz. Model results reproduce measured vapor flow, liquid pressure and hydrochemical data and capture the trends of measured relative humidities, although such data are slightly overestimated near the rock interface due to uncertainties in the turbulence factor. Rock desaturation allows oxygen to diffuse into the rock and triggers pyrite oxidation, dissolution of calcite and siderite, precipitation of ferrihydrite, dolomite and gypsum and cation exchange. pH in the unsaturated rock varies from 7.8 to 8 and is buffered by calcite. Computed changes in the porosity and the permeability of Opalinus clay in the unsaturated zone caused by oxidation and mineral dissolution/precipitation are smaller than 5%. Therefore, rock properties are not expected to be affected significantly by ventilation of underground drifts during construction and operational phases of a HLW repository in clay.

  14. Reactive chemical transport in ground-water hydrology: Challenges to mathematical modeling

    SciTech Connect (OSTI)

    Narasimhan, T.N.; Apps, J.A.

    1990-07-01

    For a long time, earth scientists have qualitatively recognized that mineral assemblages in soils and rocks conform to established principles of chemistry. In the early 1960's geochemists began systematizing this knowledge by developing quantitative thermodynamic models based on equilibrium considerations. These models have since been coupled with advective-dispersive-diffusive transport models, already developed by ground-water hydrologists. Spurred by a need for handling difficult environmental issues related to ground-water contamination, these models are being improved, refined and applied to realistic problems of interest. There is little doubt that these models will play an important role in solving important problems of engineering as well as science over the coming years. Even as these models are being used practically, there is scope for their improvement and many challenges lie ahead. In addition to improving the conceptual basis of the governing equations, much remains to be done to incorporate kinetic processes and biological mediation into extant chemical equilibrium models. Much also remains to be learned about the limits to which model predictability can be reasonably taken. The purpose of this paper is to broadly assess the current status of knowledge in modeling reactive chemical transport and to identify the challenges that lie ahead.

  15. Geophysical monitoring and reactive transport modeling of ureolytically-driven calcium carbonate precipitation

    SciTech Connect (OSTI)

    Wu, Y.; Ajo-Franklin, J.B.; Spycher, N.; Hubbard, S.S.; Zhang, G.; Williams, K.H.; Taylor, J.; Fujita, Y.; Smith, R.

    2011-07-15

    Ureolytically-driven calcium carbonate precipitation is the basis for a promising in-situ remediation method for sequestration of divalent radionuclide and trace metal ions. It has also been proposed for use in geotechnical engineering for soil strengthening applications. Monitoring the occurrence, spatial distribution, and temporal evolution of calcium carbonate precipitation in the subsurface is critical for evaluating the performance of this technology and for developing the predictive models needed for engineering application. In this study, we conducted laboratory column experiments using natural sediment and groundwater to evaluate the utility of geophysical (complex resistivity and seismic) sensing methods, dynamic synchrotron x-ray computed tomography (micro-CT), and reactive transport modeling for tracking ureolytically-driven calcium carbonate precipitation processes under site relevant conditions. Reactive transport modeling with TOUGHREACT successfully simulated the changes of the major chemical components during urea hydrolysis. Even at the relatively low level of urea hydrolysis observed in the experiments, the simulations predicted an enhanced calcium carbonate precipitation rate that was 3-4 times greater than the baseline level. Reactive transport modeling results, geophysical monitoring data and micro-CT imaging correlated well with reaction processes validated by geochemical data. In particular, increases in ionic strength of the pore fluid during urea hydrolysis predicted by geochemical modeling were successfully captured by electrical conductivity measurements and confirmed by geochemical data. The low level of urea hydrolysis and calcium carbonate precipitation suggested by the model and geochemical data was corroborated by minor changes in seismic P-wave velocity measurements and micro-CT imaging; the latter provided direct evidence of sparsely distributed calcium carbonate precipitation. Ion exchange processes promoted through NH{sub 4}{sup

  16. Geophysical Monitoring and Reactive Transport Modeling of Ureolytically-Driven Calcium Carbonate Precipitation

    SciTech Connect (OSTI)

    Yuxin Wu; Jonathan B. Ajo-Franklin; Nicolas Spycher; Susan S. Hubbard; Guoxiang Zhang; Kenneth H. Williams; Joanna Taylor; Yoshiko Fujita; Robert Smith

    2011-09-01

    Ureolytically-driven calcium carbonate precipitation is the basis for a promising in-situ remediation method for sequestration of divalent radionuclide and trace metal ions. It has also been proposed for use in geotechnical engineering for soil strengthening applications. Monitoring the occurrence, spatial distribution, and temporal evolution of calcium carbonate precipitation in the subsurface is critical for evaluating the performance of this technology and for developing the predictive models needed for engineering application. In this study, we conducted laboratory column experiments using natural sediment and groundwater to evaluate the utility of geophysical (complex resistivity and seismic) sensing methods, dynamic synchrotron x-ray computed tomography (micro-CT), and reactive transport modeling for tracking ureolytically-driven calcium carbonate precipitation processes under site relevant conditions. Reactive transport modeling with TOUGHREACT successfully simulated the changes of the major chemical components during urea hydrolysis. Even at the relatively low level of urea hydrolysis observed in the experiments, the simulations predicted an enhanced calcium carbonate precipitation rate that was 3-4 times greater than the baseline level. Reactive transport modeling results, geophysical monitoring data and micro-CT imaging correlated well with reaction processes validated by geochemical data. In particular, increases in ionic strength of the pore fluid during urea hydrolysis predicted by geochemical modeling were successfully captured by electrical conductivity measurements and confirmed by geochemical data. The low level of urea hydrolysis and calcium carbonate precipitation suggested by the model and geochemical data was corroborated by minor changes in seismic P-wave velocity measurements and micro-CT imaging; the latter provided direct evidence of sparsely distributed calcium carbonate precipitation. Ion exchange processes promoted through NH4

  17. Subsurface Multiphase Flow and Multicomponent Reactive Transport Modeling using High-Performance Computing

    SciTech Connect (OSTI)

    Hammond, Glenn E.; Lichtner, Peter C.; Lu, Chuan

    2007-07-16

    Numerical modeling has become a critical tool to the U.S. Department of Energy for evaluating the environmental impact of alternative energy sources and remediation strategies for legacy waste sites. Unfortunately, the physical and chemical complexity of many sites overwhelms the capabilities of even most state of the art groundwater models. Of particular concern are the representation of highly-heterogeneous stratified rock/soil layers in the subsurface and the biological and geochemical interactions of chemical species within multiple fluid phases. Clearly, there is a need for higher-resolution modeling (i.e. more spatial, temporal, and chemical degrees of freedom) and increasingly mechanistic descriptions of subsurface physicochemical processes. We present SciDAC-funded research being performed in the development of PFLOTRAN, a parallel multiphase flow and multicomponent reactive transport model. Written in Fortran90, PFLOTRAN is founded upon PETSc data structures and solvers. We are employing PFLOTRAN in the simulation of uranium transport at the Hanford 300 Area, a contaminated site of major concern to the Department of Energy, the State of Washington, and other government agencies. By leveraging the billions of degrees of freedom available through high-performance computation using tens of thousands of processors, we can better characterize the release of uranium into groundwater and its subsequent transport to the Columbia River, and thereby better understand and evaluate the effectiveness of various proposed remediation strategies.

  18. Subsurface Multiphase Flow and Multicomponent Reactive Transport Modeling using High-Performance Computing

    SciTech Connect (OSTI)

    Hammond, Glenn E.; Lichtner, Peter C.; Lu, Chuan

    2007-08-01

    Numerical modeling has become a critical tool to the Department of Energy for evaluating the environmental impact of alternative energy sources and remediation strategies for legacy waste sites. Unfortunately, the physical and chemical complexity of many sites overwhelms the capabilities of even most state of the art groundwater models. Of particular concern are the representation of highly-heterogeneous stratified rock/soil layers in the subsurface and the biological and geochemical interactions of chemical species within multiple fluid phases. Clearly, there is a need for higher-resolution modeling (i.e. more spatial, temporal, and chemical degrees of freedom) and increasingly mechanistic descriptions of subsurface physicochemical processes. We present research being performed in the development of PFLOTRAN, a parallel multiphase flow and multicomponent reactive transport model. Written in Fortran90, PFLOTRAN is founded upon PETSc data structures and solvers and has exhibited impressive strong scalability on up to 4000 processors on the ORNL Cray XT3. We are employing PFLOTRAN in the simulation of uranium transport at the Hanford 300 Area, a contaminated site of major concern to the Department of Energy, the State of Washington, and other government agencies where overly-simplistic historical modeling erroneously predicted decade removal times for uranium by ambient groundwater flow. By leveraging the billions of degrees of freedom available through high-performance computation using tens of thousands of processors, we can better characterize the release of uranium into groundwater and its subsequent transport to the Columbia River, and thereby better understand and evaluate the effectiveness of various proposed remediation strategies.

  19. Reaction-Based Reactive Transport Modeling of Fe(III) and U(V) Reduction

    SciTech Connect (OSTI)

    Burgos, William D.; Roden, Eric E.; Yeh, Gour-Tsyh

    2005-06-01

    Our new research project (started Fall 2004) was funded by a grant to The Pennsylvania State University, University of Central Florida, and The University of Alabama in the Integrative Studies Element of the NABIR Program (DE-FG04-ER63914/63915/63196). Our previous NABIR project (DE-FG02-01ER63180/63181/63182, funded within the Biotransformation Element) focused on (1) microbial reduction of Fe(III) and U(VI) individually, and concomitantly in natural sediments, (2) Fe(III) oxide surface chemistry, specifically with respect to reactions with Fe(II) and U(VI), (3) the influence of humic substances on Fe(III) and U(VI) bioreduction, and on U(VI) complexation, and (4) the development of reaction-based reactive transport biogeochemical models to numerically simulate our experimental results. The new project focuses on the development of a mechanistic understanding and quantitative models of coupled Fe(III)/U(VI) reduction in FRC Area 2 sediments. This work builds on our previous studies of microbial Fe(III) and U(VI) reduction, and is directly aligned with the Scheibe et al. NABIR FRC Field Project at Area 2.

  20. The Development and Application of Reactive Transport Modeling Techniques to Study Radionuclide Migration at Yucca Mountain, NV

    SciTech Connect (OSTI)

    Viswanathan, Hari Selvi

    1999-09-01

    Yucca Mountain, Nevada has been chosen as a possible site for the first high level radioactive waste repository in the United States. As part of the site investigation studies, we need to make scientifically rigorous estimations of radionuclide migration in the event of a repository breach. Performance assessment models used to make these estimations are computationally intensive. We have developed two reactive transport modeling techniques to simulate radionuclide transport at Yucca Mountain: (1) the selective coupling approach applied to the convection-dispersion-reaction (CDR) model and (2) a reactive stream tube approach (RST). These models were designed to capture the important processes that influence radionuclide migration while being computationally efficient. The conventional method of modeling reactive transport models is to solve a coupled set of multi-dimensional partial differential equations for the relevant chemical components in the system. We have developed an iterative solution technique, denoted the selective coupling method, that represents a versatile alternative to traditional uncoupled iterative techniques and the filly coupled global implicit method. We show that selective coupling results in computational and memory savings relative to these approaches. We develop RST as an alternative to the CDR method for solving large two- or three-dimensional reactive transport simulations for cases in which one is interested in predicting the flux across a specific control plane. In the RST method, the multidimensional problem is reduced to a series of one-dimensional transport simulations along streamlines. The key assumption with RST is that mixing at the control plane approximates the transverse dispersion between streamlines. We compare the CDR and RST approaches for several scenarios that are relevant to the Yucca Mountain Project. For example, we apply the CDR and RST approaches to model an ongoing field experiment called the Unsaturated Zone

  1. Development of Modeling Methods and Tools for Predicting Coupled Reactive Transport Processes in Porous Media at Multiple Scales

    SciTech Connect (OSTI)

    Clement, T Prabhakar; Barnett, Mark O; Zheng, Chunmiao; Jones, Norman L

    2010-05-05

    DE-FG02-06ER64213: Development of Modeling Methods and Tools for Predicting Coupled Reactive Transport Processes in Porous Media at Multiple Scales Investigators: T. Prabhakar Clement (PD/PI) and Mark O. Barnett (Auburn), Chunmiao Zheng (Univ. of Alabama), and Norman L. Jones (BYU). The objective of this project was to develop scalable modeling approaches for predicting the reactive transport of metal contaminants. We studied two contaminants, a radioactive cation [U(VI)] and a metal(loid) oxyanion system [As(III/V)], and investigated their interactions with two types of subsurface materials, iron and manganese oxyhydroxides. We also developed modeling methods for describing the experimental results. Overall, the project supported 25 researchers at three universities. Produced 15 journal articles, 3 book chapters, 6 PhD dissertations and 6 MS theses. Three key journal articles are: 1) Jeppu et al., A scalable surface complexation modeling framework for predicting arsenate adsorption on goethite-coated sands, Environ. Eng. Sci., 27(2): 147-158, 2010. 2) Loganathan et al., Scaling of adsorption reactions: U(VI) experiments and modeling, Applied Geochemistry, 24 (11), 2051-2060, 2009. 3) Phillippi, et al., Theoretical solid/solution ratio effects on adsorption and transport: uranium (VI) and carbonate, Soil Sci. Soci. of America, 71:329-335, 2007

  2. Reactive transport modeling of stable carbon isotope fractionation in a multi-phase multi-component system during carbon sequestration

    SciTech Connect (OSTI)

    Zhang, Shuo [Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States); DePaolo, Donald J. [Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States); Zheng, Liange [Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States); Mayer, Bernhard [Univ. of Calgary (Canada). Dept. of Geosciences

    2014-12-31

    Carbon stable isotopes can be used in characterization and monitoring of CO2 sequestration sites to track the migration of the CO2 plume and identify leakage sources, and to evaluate the chemical reactions that take place in the CO2-water-rock system. However, there are few tools available to incorporate stable isotope information into flow and transport codes used for CO2 sequestration problems. We present a numerical tool for modeling the transport of stable carbon isotopes in multiphase reactive systems relevant to geologic carbon sequestration. The code is an extension of the reactive transport code TOUGHREACT. The transport module of TOUGHREACT was modified to include separate isotopic species of CO2 gas and dissolved inorganic carbon (CO2, CO32-, HCO3-,). Any process of transport or reaction influencing a given carbon species also influences its isotopic ratio. Isotopic fractionation is thus fully integrated within the dynamic system. The chemical module and database have been expanded to include isotopic exchange and fractionation between the carbon species in both gas and aqueous phases. The performance of the code is verified by modeling ideal systems and comparing with theoretical results. Efforts are also made to fit field data from the Pembina CO2 injection project in Canada. We show that the exchange of carbon isotopes between dissolved and gaseous carbon species combined with fluid flow and transport, produce isotopic effects that are significantly different from simple two-component mixing. These effects are important for understanding the isotopic variations observed in field demonstrations.

  3. Reactive transport modeling of stable carbon isotope fractionation in a multi-phase multi-component system during carbon sequestration

    DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)

    Zhang, Shuo; DePaolo, Donald J.; Zheng, Liange; Mayer, Bernhard

    2014-12-31

    Carbon stable isotopes can be used in characterization and monitoring of CO2 sequestration sites to track the migration of the CO2 plume and identify leakage sources, and to evaluate the chemical reactions that take place in the CO2-water-rock system. However, there are few tools available to incorporate stable isotope information into flow and transport codes used for CO2 sequestration problems. We present a numerical tool for modeling the transport of stable carbon isotopes in multiphase reactive systems relevant to geologic carbon sequestration. The code is an extension of the reactive transport code TOUGHREACT. The transport module of TOUGHREACT was modifiedmore » to include separate isotopic species of CO2 gas and dissolved inorganic carbon (CO2, CO32-, HCO3-,…). Any process of transport or reaction influencing a given carbon species also influences its isotopic ratio. Isotopic fractionation is thus fully integrated within the dynamic system. The chemical module and database have been expanded to include isotopic exchange and fractionation between the carbon species in both gas and aqueous phases. The performance of the code is verified by modeling ideal systems and comparing with theoretical results. Efforts are also made to fit field data from the Pembina CO2 injection project in Canada. We show that the exchange of carbon isotopes between dissolved and gaseous carbon species combined with fluid flow and transport, produce isotopic effects that are significantly different from simple two-component mixing. These effects are important for understanding the isotopic variations observed in field demonstrations.« less

  4. Pore Scale Modeling of the Reactive Transport of Chromium in the Cathode of a Solid Oxide Fuel Cell

    SciTech Connect (OSTI)

    Ryan, Emily M.; Tartakovsky, Alexandre M.; Recknagle, Kurtis P.; Khaleel, Mohammad A.; Amon, Cristina

    2011-01-01

    We present a pore scale model of a solid oxide fuel cell (SOFC) cathode. Volatile chromium species are known to migrate from the current collector of the SOFC into the cathode where over time they decrease the voltage output of the fuel cell. A pore scale model is used to investigate the reactive transport of chromium species in the cathode and to study the driving forces of chromium poisoning. A multi-scale modeling approach is proposed which uses a cell level model of the cathode, air channel and current collector to determine the boundary conditions for a pore scale model of a section of the cathode. The pore scale model uses a discrete representation of the cathode to explicitly model the surface reactions of oxygen and chromium with a cathode material. The pore scale model is used to study the reaction mechanisms of chromium by considering the effects of reaction rates, diffusion coefficients, chromium vaporization, and oxygen consumption on chromiums deposition in the cathode. The study shows that chromium poisoning is most significantly affected by the chromium reaction rates in the cathode and that the reaction rates are a function of the local current density in the cathode.

  5. Building Conceptual Models of Field-Scale Uranium Reactive Transport in a Dynamic Vadose Zone-Aquifer-River System

    SciTech Connect (OSTI)

    Yabusaki, Steven B.; Fang, Yilin; Waichler, Scott R.

    2008-12-04

    Subsurface simulation is being used to build, test, and couple conceptual process models to better understand controls on a 0.4 km by 1.0 km uranium plume that has persisted above the drinking water standard in the groundwater of the Hanford 300 Area over the last 15 years. At this site, uranium-contaminated sediments in the vadose zone and aquifer are subject to significant variations in water levels and velocities driven by the diurnal, weekly, seasonal, and episodic Columbia River stage dynamics. Groundwater flow reversals typically occur twice a day with significant exchange of river water and groundwater in the near-river aquifer. Mixing of the dilute solution chemistry of the river with the groundwater complicates the uranium sorption behavior as the mobility of U(VI) has been shown experimentally to be a function of pH, carbonate, calcium, and uranium. Furthermore, uranium mass transfer between solid and aqueous phases has been observed to be rate-limited in the context of the high groundwater velocities resulting from the river stage fluctuations and the highly transmissive sediments (hydraulic conductivities ~1500 m/d). One- and two-dimensional vertical cross-sectional simulations of variably-saturated flow and reactive transport, based on laboratory-derived models of distributed rate mass transfer and equilibrium multicomponent surface complexation, are used to assess uranium transport at the dynamic vadose zone aquifer interface as well as changes to uranium mobility due to incursions of river water into the aquifer.

  6. Review of existing reactive transport software

    SciTech Connect (OSTI)

    Glassley, W., LLNL

    1998-02-03

    Simulations of thermal and hydrological evolution following the potential emplacement of a subterranean nuclear waste repository at Yucca Mountain, NV provide data that suggest the inevitability of dependent, simultaneous chemical evolution in this system. These chemical changes will modify significantly both the magnitude and structure of local porosity and permeability; hence, they will have a dynamic feedback effect on the evolving thermal and hydrological regime. Yet, despite this intimate interdependence of transport and chemical processes, a rigorous quantitative analysis of the post- emplacement environment that incorporates this critical feedback mechanism has not been completed to date. As an initial step in this direction, the present document outlines the fundamental chemical and transport processes that must be accounted for in such an analysis, and reviews the inventory of existing software that encodes these processed in explicitly coupled form. A companion report describes the prioritization of specific capabilities that are needed for modeling post-emplacement reactive transport at Yucca Mountain.

  7. quantifying and Predicting Reactive Transport

    SciTech Connect (OSTI)

    Peter C. Burns, Department of Civil Engineering and Geological Sciences, University of Notre Dame

    2009-12-04

    This project was led by Dr. Jiamin Wan at Lawrence Berkeley National Laboratory. Peter Burns provided expertise in uranium mineralogy and in identification of uranium minerals in test materials. Dr. Wan conducted column tests regarding uranium transport at LBNL, and samples of the resulting columns were sent to Dr. Burns for analysis. Samples were analyzed for uranium mineralogy by X-ray powder diffraction and by scanning electron microscopy, and results were provided to Dr. Wan for inclusion in the modeling effort. Full details of the project can be found in Dr. Wan's final reports for the associated effort at LBNL.

  8. Modeling reactive transport in deformable porous media using the theory of interacting continua.

    SciTech Connect (OSTI)

    Turner, Daniel Zack

    2012-01-01

    This report gives an overview of the work done as part of an Early Career LDRD aimed at modeling flow induced damage of materials involving chemical reactions, deformation of the porous matrix, and complex flow phenomena. The numerical formulation is motivated by a mixture theory or theory of interacting continua type approach to coupling the behavior of the fluid and the porous matrix. Results for the proposed method are presented for several engineering problems of interest including carbon dioxide sequestration, hydraulic fracturing, and energetic materials applications. This work is intended to create a general framework for flow induced damage that can be further developed in each of the particular areas addressed below. The results show both convincing proof of the methodologies potential and the need for further validation of the models developed.

  9. A Uranium Bioremediation Reactive Transport Benchmark

    SciTech Connect (OSTI)

    Yabusaki, Steven B.; Sengor, Sevinc; Fang, Yilin

    2015-06-01

    A reactive transport benchmark problem set has been developed based on in situ uranium bio-immobilization experiments that have been performed at a former uranium mill tailings site in Rifle, Colorado, USA. Acetate-amended groundwater stimulates indigenous microorganisms to catalyze the reduction of U(VI) to a sparingly soluble U(IV) mineral. The interplay between the flow, acetate loading periods and rates, microbially-mediated and geochemical reactions leads to dynamic behavior in metal- and sulfate-reducing bacteria, pH, alkalinity, and reactive mineral surfaces. The benchmark is based on an 8.5 m long one-dimensional model domain with constant saturated flow and uniform porosity. The 159-day simulation introduces acetate and bromide through the upgradient boundary in 14-day and 85-day pulses separated by a 10 day interruption. Acetate loading is tripled during the second pulse, which is followed by a 50 day recovery period. Terminal electron accepting processes for goethite, phyllosilicate Fe(III), U(VI), and sulfate are modeled using Monod-type rate laws. Major ion geochemistry modeled includes mineral reactions, as well as aqueous and surface complexation reactions for UO2++, Fe++, and H+. In addition to the dynamics imparted by the transport of the acetate pulses, U(VI) behavior involves the interplay between bioreduction, which is dependent on acetate availability, and speciation-controlled surface complexation, which is dependent on pH, alkalinity and available surface complexation sites. The general difficulty of this benchmark is the large number of reactions (74), multiple rate law formulations, a multisite uranium surface complexation model, and the strong interdependency and sensitivity of the reaction processes. Results are presented for three simulators: HYDROGEOCHEM, PHT3D, and PHREEQC.

  10. Using reactive transport codes to provide mechanistic biogeochemistry representations in global land surface models: CLM-PFLOTRAN 1.0

    DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)

    Tang, G.; Yuan, F.; Bisht, G.; Hammond, G. E.; Lichtner, P. C.; Kumar, J.; Mills, R. T.; Xu, X.; Andre, B.; Hoffman, F. M.; et al

    2015-12-17

    We explore coupling to a configurable subsurface reactive transport code as a flexible and extensible approach to biogeochemistry in land surface models; our goal is to facilitate testing of alternative models and incorporation of new understanding. A reaction network with the CLM-CN decomposition, nitrification, denitrification, and plant uptake is used as an example. We implement the reactions in the open-source PFLOTRAN code, coupled with the Community Land Model (CLM), and test at Arctic, temperate, and tropical sites. To make the reaction network designed for use in explicit time stepping in CLM compatible with the implicit time stepping used in PFLOTRAN,more » the Monod substrate rate-limiting function with a residual concentration is used to represent the limitation of nitrogen availability on plant uptake and immobilization. To achieve accurate, efficient, and robust numerical solutions, care needs to be taken to use scaling, clipping, or log transformation to avoid negative concentrations during the Newton iterations. With a tight relative update tolerance to avoid false convergence, an accurate solution can be achieved with about 50 % more computing time than CLM in point mode site simulations using either the scaling or clipping methods. The log transformation method takes 60–100 % more computing time than CLM. The computing time increases slightly for clipping and scaling; it increases substantially for log transformation for half saturation decrease from 10−3 to 10−9 mol m−3, which normally results in decreasing nitrogen concentrations. The frequent occurrence of very low concentrations (e.g. below nanomolar) can increase the computing time for clipping or scaling by about 20 %; computing time can be doubled for log transformation. Caution needs to be taken in choosing the appropriate scaling factor because a small value caused by a negative update to a small concentration may diminish the update and result in false convergence even with very

  11. REDUCING UNCERTAINTIES IN MODEL PREDICTIONS VIA HISTORY MATCHING OF CO2 MIGRATION AND REACTIVE TRANSPORT MODELING OF CO2 FATE AT THE SLEIPNER PROJECT

    SciTech Connect (OSTI)

    Zhu, Chen

    2015-03-31

    An important question for the Carbon Capture, Storage, and Utility program is “can we adequately predict the CO2 plume migration?” For tracking CO2 plume development, the Sleipner project in the Norwegian North Sea provides more time-lapse seismic monitoring data than any other sites, but significant uncertainties still exist for some of the reservoir parameters. In Part I, we assessed model uncertainties by applying two multi-phase compositional simulators to the Sleipner Benchmark model for the uppermost layer (Layer 9) of the Utsira Sand and calibrated our model against the time-lapsed seismic monitoring data for the site from 1999 to 2010. Approximate match with the observed plume was achieved by introducing lateral permeability anisotropy, adding CH4 into the CO2 stream, and adjusting the reservoir temperatures. Model-predicted gas saturation, CO2 accumulation thickness, and CO2 solubility in brine—none were used as calibration metrics—were all comparable with the interpretations of the seismic data in the literature. In Part II & III, we evaluated the uncertainties of predicted long-term CO2 fate up to 10,000 years, due to uncertain reaction kinetics. Under four scenarios of the kinetic rate laws, the temporal and spatial evolution of CO2 partitioning into the four trapping mechanisms (hydrodynamic/structural, solubility, residual/capillary, and mineral) was simulated with ToughReact, taking into account the CO2-brine-rock reactions and the multi-phase reactive flow and mass transport. Modeling results show that different rate laws for mineral dissolution and precipitation reactions resulted in different predicted amounts of trapped CO2 by carbonate minerals, with scenarios of the conventional linear rate law for feldspar dissolution having twice as much mineral trapping (21% of the injected CO2) as scenarios with a Burch-type or Alekseyev et al.–type rate law for feldspar dissolution (11%). So far, most reactive transport modeling (RTM) studies for

  12. Addressing numerical challenges in introducing a reactive transport code into a land surface model: a biogeochemical modeling proof-of-concept with CLM–PFLOTRAN 1.0

    DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)

    Tang, Guoping; Yuan, Fengming; Bisht, Gautam; Hammond, Glenn E.; Lichtner, Peter C.; Kumar, Jitendra; Mills, Richard T.; Xu, Xiaofeng; Andre, Ben; Hoffman, Forrest M.; et al

    2016-03-04

    We explore coupling to a configurable subsurface reactive transport code as a flexible and extensible approach to biogeochemistry in land surface models. A reaction network with the Community Land Model carbon–nitrogen (CLM-CN) decomposition, nitrification, denitrification, and plant uptake is used as an example. We implement the reactions in the open-source PFLOTRAN (massively parallel subsurface flow and reactive transport) code and couple it with the CLM. To make the rate formulae designed for use in explicit time stepping in CLMs compatible with the implicit time stepping used in PFLOTRAN, the Monod substrate rate-limiting function with a residual concentration is used to represent the limitation ofmore » nitrogen availability on plant uptake and immobilization. We demonstrate that CLM–PFLOTRAN predictions (without invoking PFLOTRAN transport) are consistent with CLM4.5 for Arctic, temperate, and tropical sites.Switching from explicit to implicit method increases rigor but introduces numerical challenges. Care needs to be taken to use scaling, clipping, or log transformation to avoid negative concentrations during the Newton iterations. With a tight relative update tolerance (STOL) to avoid false convergence, an accurate solution can be achieved with about 50 % more computing time than CLM in point mode site simulations using either the scaling or clipping methods. The log transformation method takes 60–100 % more computing time than CLM. The computing time increases slightly for clipping and scaling; it increases substantially for log transformation for half saturation decrease from 10−3 to 10−9 mol m−3, which normally results in decreasing nitrogen concentrations. The frequent occurrence of very low concentrations (e.g. below nanomolar) can increase the computing time for clipping or scaling by about 20 %, double for log transformation. Overall, the log transformation method is accurate and robust, and the clipping and scaling

  13. Addressing numerical challenges in introducing a reactive transport code into a land surface model: A biogeochemical modeling proof-of-concept with CLM PFLOTRAN 1.0

    DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)

    Tang, Guoping; Yuan, Fengming; Bisht, Gautam; Hammond, Glenn E.; Lichtner, Peter C.; Kumar, Jitendra; Mills, Richard T.; Xu, Xiaofeng; Andre, Ben; Hoffman, Forrest M.; et al

    2016-03-04

    Here, we explore coupling to a configurable subsurface reactive transport code as a flexible and extensible approach to biogeochemistry in land surface models. A reaction network with the Community Land Model carbon nitrogen (CLM-CN) decomposition, nitrification, denitrification, and plant uptake is used as an example. We implement the reactions in the open-source PFLOTRAN (massively parallel subsurface flow and reactive transport) code and couple it with the CLM. To make the rate formulae designed for use in explicit time stepping in CLMs compatible with the implicit time stepping used in PFLOTRAN, the Monod substrate rate-limiting function with a residual concentration ismore » used to represent the limitation of nitrogen availability on plant uptake and immobilization. We demonstrate that CLM PFLOTRAN predictions (without invoking PFLOTRAN transport) are consistent with CLM4.5 for Arctic, temperate, and tropical sites. Switching from explicit to implicit method increases rigor but introduces numerical challenges. Care needs to be taken to use scaling, clipping, or log transformation to avoid negative concentrations during the Newton iterations. With a tight relative update tolerance (STOL) to avoid false convergence, an accurate solution can be achieved with about 50 % more computing time than CLM in point mode site simulations using either the scaling or clipping methods. The log transformation method takes 60–100 % more computing time than CLM. The computing time increases slightly for clipping and scaling; it increases substantially for log transformation for half saturation decrease from 10–3 to 10–9 mol m–3, which normally results in decreasing nitrogen concentrations. The frequent occurrence of very low concentrations (e.g. below nanomolar) can increase the computing time for clipping or scaling by about 20 %, double for log transformation. Overall, the log transformation method is accurate and robust, and the clipping and scaling methods are

  14. Upscaling of Long-Term U9VI) Desorption from Pore Scale Kinetics to Field-Scale Reactive Transport Models

    SciTech Connect (OSTI)

    Andy Miller

    2009-01-25

    Environmental systems exhibit a range of complexities which exist at a range of length and mass scales. Within the realm of radionuclide fate and transport, much work has been focused on understanding pore scale processes where complexity can be reduced to a simplified system. In describing larger scale behavior, the results from these simplified systems must be combined to create a theory of the whole. This process can be quite complex, and lead to models which lack transparency. The underlying assumption of this approach is that complex systems will exhibit complex behavior, requiring a complex system of equations to describe behavior. This assumption has never been tested. The goal of the experiments presented is to ask the question: Do increasingly complex systems show increasingly complex behavior? Three experimental tanks at the intermediate scale (Tank 1: 2.4m x 1.2m x 7.6cm, Tank 2: 2.4m x 0.61m x 7.6cm, Tank 3: 2.4m x 0.61m x 0.61m (LxHxW)) have been completed. These tanks were packed with various physical orientations of different particle sizes of a uranium contaminated sediment from a former uranium mill near Naturita, Colorado. Steady state water flow was induced across the tanks using constant head boundaries. Pore water was removed from within the flow domain through sampling ports/wells; effluent samples were also taken. Each sample was analyzed for a variety of analytes relating to the solubility and transport of uranium. Flow fields were characterized using inert tracers and direct measurements of pressure head. The results show that although there is a wide range of chemical variability within the flow domain of the tank, the effluent uranium behavior is simple enough to be described using a variety of conceptual models. Thus, although there is a wide range in variability caused by pore scale behaviors, these behaviors appear to be smoothed out as uranium is transported through the tank. This smoothing of uranium transport behavior transcends

  15. Fluid-rock interaction: A reactive transport approach

    SciTech Connect (OSTI)

    Steefel, C.; Maher, K.

    2009-04-01

    Fluid-rock interaction (or water-rock interaction, as it was more commonly known) is a subject that has evolved considerably in its scope over the years. Initially its focus was primarily on interactions between subsurface fluids of various temperatures and mostly crystalline rocks, but the scope has broadened now to include fluid interaction with all forms of subsurface materials, whether they are unconsolidated or crystalline ('fluid-solid interaction' is perhaps less euphonious). Disciplines that previously carried their own distinct names, for example, basin diagenesis, early diagenesis, metamorphic petrology, reactive contaminant transport, chemical weathering, are now considered to fall under the broader rubric of fluid-rock interaction, although certainly some of the key research questions differ depending on the environment considered. Beyond the broadening of the environments considered in the study of fluid-rock interaction, the discipline has evolved in perhaps an even more important way. The study of water-rock interaction began by focusing on geochemical interactions in the absence of transport processes, although a few notable exceptions exist (Thompson 1959; Weare et al. 1976). Moreover, these analyses began by adopting a primarily thermodynamic approach, with the implicit or explicit assumption of equilibrium between the fluid and rock. As a result, these early models were fundamentally static rather than dynamic in nature. This all changed with the seminal papers by Helgeson and his co-workers (Helgeson 1968; Helgeson et al. 1969) wherein the concept of an irreversible reaction path was formally introduced into the geochemical literature. In addition to treating the reaction network as a dynamically evolving system, the Helgeson studies introduced an approach that allowed for the consideration of a multicomponent geochemical system, with multiple minerals and species appearing as both reactants and products, at least one of which could be

  16. Using laboratory flow experiments and reactive chemical transport modeling for designing waterflooding of the Agua Fria Reservoir, Poza Rica-Altamira Field, Mexico

    SciTech Connect (OSTI)

    Birkle, P.; Pruess, K.; Xu, T.; Figueroa, R.A. Hernandez; Lopez, M. Diaz; Lopez, E. Contreras

    2008-10-01

    Waterflooding for enhanced oil recovery requires that injected waters must be chemically compatible with connate reservoir waters, in order to avoid mineral dissolution-and-precipitation cycles that could seriously degrade formation permeability and injectivity. Formation plugging is a concern especially in reservoirs with a large content of carbonates, such as calcite and dolomite, as such minerals typically react rapidly with an aqueous phase, and have strongly temperature-dependent solubility. Clay swelling can also pose problems. During a preliminary waterflooding pilot project, the Poza Rica-Altamira oil field, bordering the Gulf coast in the eastern part of Mexico, experienced injectivity loss after five months of reinjection of formation waters into well AF-847 in 1999. Acidizing with HCl restored injectivity. We report on laboratory experiments and reactive chemistry modeling studies that were undertaken in preparation for long-term waterflooding at Agua Frma. Using analogous core plugs obtained from the same reservoir interval, laboratory coreflood experiments were conducted to examine sensitivity of mineral dissolution and precipitation effects to water composition. Native reservoir water, chemically altered waters, and distilled water were used, and temporal changes in core permeability, mineral abundances and aqueous concentrations of solutes were monitored. The experiments were simulated with the multi-phase, nonisothermal reactive transport code TOUGHREACT, and reasonable to good agreement was obtained for changes in solute concentrations. Clay swelling caused an additional impact on permeability behavior during coreflood experiments, whereas the modeled permeability depends exclusively on chemical processes. TOUGHREACT was then used for reservoir-scale simulation of injecting ambient-temperature water (30 C, 86 F) into a reservoir with initial temperature of 80 C (176 F). Untreated native reservoir water was found to cause serious porosity and

  17. System-Scale Model of Aquifer, Vadose Zone, and River Interactions for the Hanford 300 Area - Application to Uranium Reactive Transport

    SciTech Connect (OSTI)

    Rockhold, Mark L.; Bacon, Diana H.; Freedman, Vicky L.; Parker, Kyle R.; Waichler, Scott R.; Williams, Mark D.

    2013-10-01

    This report represents a synthesis and integration of basic and applied research into a system-scale model of the Hanford 300 Area groundwater uranium plume, supported by the U.S. Department of Energy’s Richland Operations (DOE-RL) office. The report integrates research findings and data from DOE Office of Science (DOE-SC), Office of Environmental Management (DOE-EM), and DOE-RL projects, and from the site remediation and closure contractor, Washington Closure Hanford, LLC (WCH). The three-dimensional, system-scale model addresses water flow and reactive transport of uranium for the coupled vadose zone, unconfined aquifer, and Columbia River shoreline of the Hanford 300 Area. The system-scale model of the 300 Area was developed to be a decision-support tool to evaluate processes of the total system affecting the groundwater uranium plume. The model can also be used to address “what if” questions regarding different remediation endpoints, and to assist in design and evaluation of field remediation efforts. For example, the proposed cleanup plan for the Hanford 300 Area includes removal, treatment, and disposal of contaminated sediments from known waste sites, enhanced attenuation of uranium hot spots in the vadose and periodically rewetted zone, and continued monitoring of groundwater with institutional controls. Illustrative simulations of polyphosphate infiltration were performed to demonstrate the ability of the system-scale model to address these types of questions. The use of this model in conjunction with continued field monitoring is expected to provide a rigorous basis for developing operational strategies for field remediation and for defining defensible remediation endpoints.

  18. Reaction-Based Reactive Transport Modeling of Iron Reduction and Uranium Immobilization at Area 2 of the NABIR Field Research Center

    SciTech Connect (OSTI)

    Burgos, W.D.

    2009-09-02

    This report summarizes research conducted in conjunction with a project entitled “Reaction-Based Reactive Transport Modeling of Iron Reduction and Uranium Immobilization at Area 2 of the NABIR Field Research Center”, which was funded through the Integrative Studies Element of the former NABIR Program (now the Environmental Remediation Sciences Program) within the Office of Biological and Environmental Research. Dr. William Burgos (The Pennsylvania State University) was the overall PI/PD for the project, which included Brian Dempsey (Penn State), Gour-Tsyh (George) Yeh (Central Florida University), and Eric Roden (formerly at The University of Alabama, now at the University of Wisconsin) as separately-funded co-PIs. The project focused on development of a mechanistic understanding and quantitative models of coupled Fe(III)/U(VI) reduction in FRC Area 2 sediments. The work builds on our previous studies of microbial Fe(III) and U(VI) reduction, and was directly aligned with the Scheibe et al. ORNL FRC Field Project at Area 2.

  19. Nonlocal reactive transport with physical and chemical heterogeneity: Localization errors

    SciTech Connect (OSTI)

    Cushman, J.H.; Hu, B.X.; Deng, F.W.

    1995-09-01

    The origin of nonlocality in {open_quotes}macroscale{close_quotes} models for subsurface chemical transport is illustrated. It is argued that media that are either nonperiodic (e.g., media with evolving heterogeneity) or periodic viewed on a scale wherein a unit cell is discernible must display some nonlocality in the mean. A metaphysical argument suggests that owing to the scarcity of information on natural scales of heterogeneity and on scales of observation associated with an instrument window, constitutive theories for the mean concentration should at the outset of any modeling effort always be considered nonlocal. The intuitive appeal to nonlocality is reinforced with an analytical derivation of the constitutive theory for a conservative tracer without appeal to any mathematical approximations. Comparisons are made between the fully nonlocal (FNL), nonlocal in time (NLT), and fully localized (FL) theories. For conservative transport, there is little difference between the first-order FL and FNL models for spatial moments up to and including the third. However, for conservative transport the first-order NLT model differs significantly from the FNL model in the third spatial moments. For reactive transport, all spatial moments differ between the FNL and FL models. The second transverse-horizontal and third longitudinal-horizontal moments for the NLT model differ from the FNL model. These results suggest that localized first-order transport models for conservative tracers are reasonable if only lower-order moments are desired. However, when the chemical reacts with its environment, the localization approximation can lead to significant error in all moments, and a FNL model will in general be required for accurate simulation. 18 refs., 9 figs., 1 tab.

  20. Reactive transport modeling to study changes in water chemistry induced by CO2 injection at the Frio-I brine pilot

    SciTech Connect (OSTI)

    Kharaka, Y.K; Doughty, C.; Freifeld, B.M.; Daley, T.M.; Xu, T.

    2009-11-01

    To demonstrate the potential for geologic storage of CO{sub 2} in saline aquifers, the Frio-I Brine Pilot was conducted, during which 1600 tons of CO{sub 2} were injected into a high-permeability sandstone and the resulting subsurface plume of CO{sub 2} was monitored using a variety of hydrogeological, geophysical, and geochemical techniques. Fluid samples were obtained before CO{sub 2} injection for baseline geochemical characterization, during the CO{sub 2} injection to track its breakthrough at a nearby observation well, and after injection to investigate changes in fluid composition and potential leakage into an overlying zone. Following CO{sub 2} breakthrough at the observation well, brine samples showed sharp drops in pH, pronounced increases in HCO{sub 3}{sup -} and aqueous Fe, and significant shifts in the isotopic compositions of H{sub 2}O and dissolved inorganic carbon. Based on a calibrated 1-D radial flow model, reactive transport modeling was performed for the Frio-I Brine Pilot. A simple kinetic model of Fe release from the solid to aqueous phase was developed, which can reproduce the observed increases in aqueous Fe concentration. Brine samples collected after half a year had lower Fe concentrations due to carbonate precipitation, and this trend can be also captured by our modeling. The paper provides a method for estimating potential mobile Fe inventory, and its bounding concentration in the storage formation from limited observation data. Long-term simulations show that the CO{sub 2} plume gradually spreads outward due to capillary forces, and the gas saturation gradually decreases due to its dissolution and precipitation of carbonates. The gas phase is predicted to disappear after 500 years. Elevated aqueous CO{sub 2} concentrations remain for a longer time, but eventually decrease due to carbonate precipitation. For the Frio-I Brine Pilot, all injected CO{sub 2} could ultimately be sequestered as carbonate minerals.

  1. End-Member Formulation of Solid Solutions and Reactive Transport

    SciTech Connect (OSTI)

    Lichtner, Peter C.

    2015-09-01

    A model for incorporating solid solutions into reactive transport equations is presented based on an end-member representation. Reactive transport equations are solved directly for the composition and bulk concentration of the solid solution. Reactions of a solid solution with an aqueous solution are formulated in terms of an overall stoichiometric reaction corresponding to a time-varying composition and exchange reactions, equivalent to reaction end-members. Reaction rates are treated kinetically using a transition state rate law for the overall reaction and a pseudo-kinetic rate law for exchange reactions. The composition of the solid solution at the onset of precipitation is assumed to correspond to the least soluble composition, equivalent to the composition at equilibrium. The stoichiometric saturation determines if the solid solution is super-saturated with respect to the aqueous solution. The method is implemented for a simple prototype batch reactor using Mathematica for a binary solid solution. Finally, the sensitivity of the results on the kinetic rate constant for a binary solid solution is investigated for reaction of an initially stoichiometric solid phase with an undersaturated aqueous solution.

  2. Transportation Systems Modeling

    Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)

    TRACC RESEARCH Computational Fluid Dynamics Computational Structural Mechanics Transportation Systems Modeling TRANSPORTATION SYSTEMS MODELING Overview of TSM Transportation systems modeling research at TRACC uses the TRANSIMS (Transportation Analysis SIMulation System) traffic micro simulation code developed by the U.S. Department of Transportation (USDOT). The TRANSIMS code represents the latest generation of traffic simulation codes developed jointly under multiyear programs by USDOT, the

  3. Reactive facies: An approach for parameterizing field-scale reactive...

    Office of Scientific and Technical Information (OSTI)

    transport models using geophysical methods Citation Details In-Document Search Title: Reactive facies: An approach for parameterizing field-scale reactive transport models ...

  4. Influence of calcite on uranium(VI) reactive transport in the groundwaterriver mixing zone

    SciTech Connect (OSTI)

    Ma, Rui; Liu, Chongxuan; Greskowiak, Janek; Prommer, Henning; Zachara, John M.; Zheng, Chunmiao

    2014-01-23

    Calcite is an important mineral that can affect uranyl reactive transport in subsurface sediments. This study investigated the distribution of calcite and its influence on uranyl adsorption and reactive transport in the groundwater-river mixing zone at US Hanford 300A, Washington State. Simulations using a 2D reactive transport model under field-relevant hydrogeochemical conditions revealed a complex distribution of calcite concentration as a result of dynamic groundwater-river interactions. The calcite concentration distribution in turn affected the spatial and temporal changes in aqueous carbonate, calcium, and pH, which subsequently influenced U(VI) mobility and discharge rates into the river. The results implied that calcite distribution and its concentration dynamics is an important consideration for field characterization, monitoring, and reactive transport prediction.

  5. Coupling a Reactive Transport Code with a Global Land Surface Model for Mechanistic Biogeochemistry Representation: 1. Addressing the Challenge of Nonnegativity

    SciTech Connect (OSTI)

    Tang, Guoping; Yuan, Fengming; Bisht, Gautam; Hammond, Glenn E.; Lichtner, Peter C.; Collier, Nathaniel O.; Kumar, Jitendra; Mills, Richard T.; Xu, Xiaofeng; Andre, Ben; Hoffman, Forrest M.; Painter, Scott L.; Thornton, Peter E.

    2016-01-01

    Reactive transport codes (e.g., PFLOTRAN) are increasingly used to improve the representation of biogeochemical processes in terrestrial ecosystem models (e.g., the Community Land Model, CLM). As CLM and PFLOTRAN use explicit and implicit time stepping, implementation of CLM biogeochemical reactions in PFLOTRAN can result in negative concentration, which is not physical and can cause numerical instability and errors. The objective of this work is to address the nonnegativity challenge to obtain accurate, efficient, and robust solutions. We illustrate the implementation of a reaction network with the CLM-CN decomposition, nitrification, denitrification, and plant nitrogen uptake reactions and test the implementation at arctic, temperate, and tropical sites. We examine use of scaling back the update during each iteration (SU), log transformation (LT), and downregulating the reaction rate to account for reactant availability limitation to enforce nonnegativity. Both SU and LT guarantee nonnegativity but with implications. When a very small scaling factor occurs due to either consumption or numerical overshoot, and the iterations are deemed converged because of too small an update, SU can introduce excessive numerical error. LT involves multiplication of the Jacobian matrix by the concentration vector, which increases the condition number, decreases the time step size, and increases the computational cost. Neither SU nor SE prevents zero concentration. When the concentration is close to machine precision or 0, a small positive update stops all reactions for SU, and LT can fail due to a singular Jacobian matrix. The consumption rate has to be downregulated such that the solution to the mathematical representation is positive. A first-order rate downregulates consumption and is nonnegative, and adding a residual concentration makes it positive. For zero-order rate or when the reaction rate is not a function of a reactant, representing the availability limitation of each

  6. Coupling a Reactive Transport Code with a Global Land Surface Model for Mechanistic Biogeochemistry Representation: 1. Addressing the Challenge of Nonnegativity

    DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)

    Tang, Guoping; Yuan, Fengming; Bisht, Gautam; Hammond, Glenn E.; Lichtner, Peter C.; Collier, Nathaniel O.; Kumar, Jitendra; Mills, Richard T.; Xu, Xiaofeng; Andre, Ben; et al

    2016-01-01

    Reactive transport codes (e.g., PFLOTRAN) are increasingly used to improve the representation of biogeochemical processes in terrestrial ecosystem models (e.g., the Community Land Model, CLM). As CLM and PFLOTRAN use explicit and implicit time stepping, implementation of CLM biogeochemical reactions in PFLOTRAN can result in negative concentration, which is not physical and can cause numerical instability and errors. The objective of this work is to address the nonnegativity challenge to obtain accurate, efficient, and robust solutions. We illustrate the implementation of a reaction network with the CLM-CN decomposition, nitrification, denitrification, and plant nitrogen uptake reactions and test the implementation atmore » arctic, temperate, and tropical sites. We examine use of scaling back the update during each iteration (SU), log transformation (LT), and downregulating the reaction rate to account for reactant availability limitation to enforce nonnegativity. Both SU and LT guarantee nonnegativity but with implications. When a very small scaling factor occurs due to either consumption or numerical overshoot, and the iterations are deemed converged because of too small an update, SU can introduce excessive numerical error. LT involves multiplication of the Jacobian matrix by the concentration vector, which increases the condition number, decreases the time step size, and increases the computational cost. Neither SU nor SE prevents zero concentration. When the concentration is close to machine precision or 0, a small positive update stops all reactions for SU, and LT can fail due to a singular Jacobian matrix. The consumption rate has to be downregulated such that the solution to the mathematical representation is positive. A first-order rate downregulates consumption and is nonnegative, and adding a residual concentration makes it positive. For zero-order rate or when the reaction rate is not a function of a reactant, representing the availability limitation

  7. Shock Desensitization Experiments and Reactive Flow Modeling...

    Office of Scientific and Technical Information (OSTI)

    Shock Desensitization Experiments and Reactive Flow Modeling on Self-Sustaining LX-17 Detonation Waves Citation Details In-Document Search Title: Shock Desensitization Experiments ...

  8. UZ Colloid Transport Model

    SciTech Connect (OSTI)

    M. McGraw

    2000-04-13

    The UZ Colloid Transport model development plan states that the objective of this Analysis/Model Report (AMR) is to document the development of a model for simulating unsaturated colloid transport. This objective includes the following: (1) use of a process level model to evaluate the potential mechanisms for colloid transport at Yucca Mountain; (2) Provide ranges of parameters for significant colloid transport processes to Performance Assessment (PA) for the unsaturated zone (UZ); (3) Provide a basis for development of an abstracted model for use in PA calculations.

  9. integrated-transportation-models

    Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)

    support a wider application of integrated transportation models, especially focusing on travel demand and network ... irrevocable worldwide license in said article to ...

  10. STOMP-ECKEChem: An Engineering Perspective on Reactive Transport in Geologic Media

    SciTech Connect (OSTI)

    White, Mark D.; Fang, Yilin

    2012-04-04

    ECKEChem (Equilibrium, Conservation, Kinetic Equation Chemistry) is a reactive transport module for the STOMP suite of multifluid subsurface flow and transport simulators that was developed from an engineering perspective. STOMP comprises a suite of operational modes that are distinguished by the solved coupled conservation equations with capabilities for a variety of subsurface applications (e.g., environmental remediation and stewardship, geologic sequestration of greenhouse gases, gas hydrate production, and oil shale production). The ECKEChem module was designed to provide integrated reactive transport capabilities across the suite of STOMP simulator operational modes. The initial application for the ECKEChem module was in the simulation of the mineralization reactions that occurred with the injection of supercritical carbon dioxide into deep Columbia River basalt formations, where it was implemented in the STOMP-CO2 simulator. The STOMP-ECKEChem solution approach to modeling reactive transport in multifluid geologic media is founded on an engineering perspective: (1) sequential non-iterative coupling between the flow and reactive transport is sufficient, (2) reactive transport can be modeled by operator splitting with local geochemistry and global transport, (3) geochemistry can be expressed as a system of coupled nonlinear equilibrium, conservation and kinetic equations, (4) a limited number of kinetic equation forms are used in geochemical practice. This chapter describes the conceptual approach to converting a geochemical reaction network into a series of equilibrium, conservation and kinetic equations, the implementation of ECKEChem in STOMP, the numerical solution approach, and a demonstration of the simulator on a complex application involving desorption of uranium from contaminated field-textured sediments.

  11. Comparison of parameter sensitivities between a laboratory and field scale model of uranium transport in a dual domain, distributed-rate reactive system

    SciTech Connect (OSTI)

    Greskowiak, Janek; Prommer, Henning; Liu, Chongxuan; Post, Vincent; Ma, Rui; Zheng, Chunmiao; Zachara, John M.

    2010-09-16

    A laboratory-derived conceptual and numerical model for U(VI) transport at the Hanford 300A site, Washington, USA, was applied to a range of field-scale scenarios of different complexity to systematically evaluate model parameter sensitivities. The model, originally developed from column experiment data, included distributed-rate surface complexation kinetics of U(VI), aqueous speciation, and physical non-equilibrium transport processes. A rigorous parameter sensitivity analysis was carried out with respect to different state variables: concentrations, mass fluxes, total mass and spatial moments of dissolved U(VI) for laboratory systems, and various simulation scenarios that represented the field-scale characteristics at the Hanford 300A site. The field-scenarios accounted for transient groundwater flow and variable geochemical conditions driven by frequent water level changes of the nearby Columbia River. Simulations indicated that the transient conditions significantly affected U(VI) plume migration at the site. The parameter sensitivities were largely similar between the laboratory and field scale systems. Where differences existed, they were shown to result from differing degrees of U(VI) adsorption disequilibrium caused by hydraulic or hydrogeochemical conditions. Adorption disequilibrium was found to differ (i) between short duration peak flow events at the field scale and much longer flow events in the laboratory, (ii) for changing groundwater chemical compositions due to river water intrusion, and (iii) for different sampling locations at the field scale. Parameter sensitivities were also found to vary with respect to the different investigated state variables. An approach is demonstrated that elucidates the most important parameters of a laboratory-scale model that must constrained in both the laboratory and field for meaningful field application.

  12. Modeling of fault reactivation and induced seismicity during...

    Office of Scientific and Technical Information (OSTI)

    Modeling of fault reactivation and induced seismicity during hydraulic fracturing of shale-gas reservoirs Citation Details In-Document Search Title: Modeling of fault reactivation ...

  13. Developing a robust geochemical and reactive transport model to evaluate possible sources of arsenic at the CO[subscript 2] sequestration natural analog site in Chimayo, New Mexico

    SciTech Connect (OSTI)

    Viswanathan, Hari; Dai, Zhenxue; Lopano, Christina; Keating, Elizabeth; Hakala, J. Alexandra; Scheckel, Kirk G.; Zheng, Liange; Gutherie, George D.; Pawar, Rajesh

    2012-10-24

    simulating the drop in As concentration as a function of time in the batch experiments. Based on modeling, kaolinite precipitation is anticipated to occur during the experiment, which allows for additional sorption sites to form with time resulting in the slow decrease in As concentration. This mechanism can be viewed as trace metal 'scavenging' due to sorption caused secondary mineral precipitation. Since deep geologic transport of these trace metals to the shallow subsurface by brine or CO{sub 2} intrusion is critical to assessing environmental impacts, the effective retardation of trace metal transport is an important parameter to estimate and it is dependent on multiple coupled reactions. At the field scale, As mobility is retarded due to the influence of sorption reactions, which can affect environmental performance assessment studies of a sequestration site.

  14. Reactive transport modeling of the enhancement of density-driven CO2 convective mixing in carbonate aquifers and its potential implication on geological carbon sequestration

    DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)

    Islam, Akand; Sun, Alexander Y.; Yang, Changbing

    2016-04-20

    We study the convection and mixing of CO2 in a brine aquifer, where the spread of dissolved CO2 is enhanced because of geochemical reactions with the host formations (calcite and dolomite), in addition to the extensively studied, buoyancy-driven mixing. The nonlinear convection is investigated under the assumptions of instantaneous chemical equilibrium, and that the dissipation of carbonate rocks solely depends on flow and transport and chemical speciation depends only on the equilibrium thermodynamics of the chemical system. The extent of convection is quantified in term of the CO2 saturation volume of the storage formation. Our results suggest that the densitymore » increase of resident species causes significant enhancement in CO2 dissolution, although no significant porosity and permeability alterations are observed. Furthermore, early saturation of the reservoir can have negative impact on CO2 sequestration.« less

  15. Final Report Coupled In Silico Microbial and Geochemical Reactive Transport Models: Extension to Multi-Organism Communities, Upscaling, and Experimental Validation

    SciTech Connect (OSTI)

    Lovley, Derek R.

    2014-03-05

    The project was highly successful in improving the understanding of physiological and ecological factors controlling the growth and interaction of subsurface microorganisms and in developing better strategies for in silico modeling of the interactions of subsurface microorganisms with other species and their environment.

  16. Chemically Reactive Working Fluids for the Capture and Transport...

    Office of Energy Efficiency and Renewable Energy (EERE) Indexed Site

    Planar Optical Waveguide Coupler Transformers for High-Power Solar Enegy Collection and Transmission Chemically Reactive Working Fluids Low-Cost Light Weigh Thin Film Solar ...

  17. Assessment of Controlling Processes for Field-Scale Uranium Reactive Transport under Highly Transient Flow Conditions

    SciTech Connect (OSTI)

    Ma, Rui; Zheng, Chunmiao; Liu, Chongxuan; Greskowiak, Janek; Prommer, Henning; Zachara, John M.

    2014-02-13

    This paper presents the results of a comprehensive model-based analysis of a uranium tracer test conducted at the U.S Department of Energy Hanford 300 Area (300A) IFRC site. A three-dimensional multi-component reactive transport model was employed to assess the key factors and processes that control the field-scale uranium reactive transport. Taking into consideration of relevant physical and chemical processes, the selected conceptual/numerical model replicates the spatial and temporal variations of the observed U(VI) concentrations reasonably well in spite of the highly complex field conditions. A sensitivity analysis was performed to interrogate the relative importance of various processes and factors for reactive transport of U(VI) at the field-scale. The results indicate that multi-rate U(VI) sorption/desorption, U(VI) surface complexation reactions, and initial U(VI) concentrations were the most important processes and factors controlling U(VI) migration. On the other hand, cation exchange reactions, the choice of the surface complexation model, and dual-domain mass transfer processes, which were previously identified to be important in laboratory experiments, played less important roles under the field-scale experimental condition at the 300A site. However, the model simulations also revealed that the groundwater chemistry was relatively stable during the uranium tracer experiment and therefore presumably not dynamic enough to appropriately assess the effects of ion exchange reaction and the choice of surface complexation models on U(VI) sorption and desorption. Furthermore, it also showed that the field experimental duration (16 days) was not sufficiently long to precisely assess the role of a majority of the sorption sites that were accessed by slow kinetic processes within the dual domain model. The sensitivity analysis revealed the crucial role of the intraborehole flow that occurred within the long-screened monitoring wells and thus significantly

  18. A sequential partly iterative approach for multicomponent reactive transport with CORE2D

    SciTech Connect (OSTI)

    Samper, J.; Xu, T.; Yang, C.

    2008-11-01

    CORE{sup 2D} V4 is a finite element code for modeling partly or fully saturated water flow, heat transport and multicomponent reactive solute transport under both local chemical equilibrium and kinetic conditions. It can handle coupled microbial processes and geochemical reactions such as acid-base, aqueous complexation, redox, mineral dissolution/precipitation, gas dissolution/exsolution, ion exchange, sorption via linear and nonlinear isotherms, sorption via surface complexation. Hydraulic parameters may change due to mineral precipitation/dissolution reactions. Coupled transport and chemical equations are solved by using sequential iterative approaches. A sequential partly-iterative approach (SPIA) is presented which improves the accuracy of the traditional sequential noniterative approach (SNIA) and is more efficient than the general sequential iterative approach (SIA). While SNIA leads to a substantial saving of computing time, it introduces numerical errors which are especially large for cation exchange reactions. SPIA improves the efficiency of SIA because the iteration between transport and chemical equations is only performed in nodes with a large mass transfer between solid and liquid phases. The efficiency and accuracy of SPIA are compared to those of SIA and SNIA using synthetic examples and a case study of reactive transport through the Llobregat Delta aquitard in Spain. SPIA is found to be as accurate as SIA while requiring significantly less CPU time. In addition, SPIA is much more accurate than SNIA with only a minor increase in computing time. A further enhancement of the efficiency of SPIA is achieved by improving the efficiency of the Newton-Raphson method used for solving chemical equations. Such an improvement is obtained by working with increments of log-concentrations and ignoring the terms of the Jacobian matrix containing derivatives of activity coefficients. A proof is given for the symmetry and non-singularity of the Jacobian matrix

  19. Nitrate Biogeochemistry and Reactive Transport in California Groundwater: LDRD Final Report

    SciTech Connect (OSTI)

    Esser, B K; Beller, H; Carle, S; Cey, B; Hudson, G B; Leif, R; LeTain, T; Moody-Bartel, C; Moore, K; McNab, W; Moran, J; Tompson, A

    2006-02-24

    Nitrate is the number one drinking water contaminant in the United States. It is pervasive in surface and groundwater systems,and its principal anthropogenic sources have increased dramatically in the last 50 years. In California alone, one third of the public drinking-water wells has been lost since 1988 and nitrate contamination is the most common reason for abandonment. Effective nitrate management in groundwater is complicated by uncertainties related to multiple point and non-point sources, hydrogeologic complexity, geochemical reactivity, and quantification of denitrification processes. In this paper, we review an integrated experimental and simulation-based framework being developed to study the fate of nitrate in a 25 km-long groundwater subbasin south of San Jose, California, a historically agricultural area now undergoing rapid urbanization with increasing demands for groundwater. The modeling approach is driven by a need to integrate new and archival data that support the hypothesis that nitrate fate and transport at the basin scale is intricately related to hydrostratigraphic complexity, variability of flow paths and groundwater residence times, microbial activity, and multiple geochemical reaction mechanisms. This study synthesizes these disparate and multi-scale data into a three-dimensional and highly resolved reactive transport modeling framework.

  20. Reactive Transport and Coupled THM Processes in Engineering Barrier Systems (EBS)

    SciTech Connect (OSTI)

    Steefel, Carl; Rutqvist, Jonny; Tsang, Chin-Fu; Liu, Hui-Hai; Sonnenthal, Eric; Houseworth, Jim; Birkholzer, Jens

    2010-08-31

    Geological repositories for disposal of high-level nuclear wastes generally rely on a multi-barrier system to isolate radioactive wastes from the biosphere. The multi-barrier system typically consists of a natural barrier system, including repository host rock and its surrounding subsurface environment, and an engineering barrier system (EBS). EBS represents the man-made, engineered materials placed within a repository, including the waste form, waste canisters, buffer materials, backfill and seals (OECD, 2003). EBS plays a significant role in the containment and long-term retardation of radionuclide release. EBS is involved in complex thermal, hydrogeological, mechanical, chemical and biological processes, such as heat release due to radionuclide decay, multiphase flow (including gas release due to canister corrosion), swelling of buffer materials, radionuclide diffusive transport, waste dissolution and chemical reactions. All these processes are related to each other. An in-depth understanding of these coupled processes is critical for the performance assessment (PA) for EBS and the entire repository. Within the EBS group of Used Fuel Disposition (UFD) Campaign, LBNL is currently focused on (1) thermal-hydraulic-mechanical-chemical (THMC) processes in buffer materials (bentonite) and (2) diffusive transport in EBS associated with clay host rock, with a long-term goal to develop a full understanding of (and needed modeling capabilities to simulate) impacts of coupled processes on radionuclide transport in different components of EBS, as well as the interaction between near-field host rock (e.g., clay) and EBS and how they effect radionuclide release. This final report documents the progress that LBNL has made in its focus areas. Specifically, Section 2 summarizes progress on literature review for THMC processes and reactive-diffusive radionuclide transport in bentonite. The literature review provides a picture of the state-of-the-art of the relevant research areas

  1. Double Shock Experiments and Reactive Flow Modeling of High Pressure...

    Office of Scientific and Technical Information (OSTI)

    Double Shock Experiments and Reactive Flow Modeling of High Pressure LX-17 Detonation Reaction Product States Citation Details In-Document Search Title: Double Shock Experiments ...

  2. Minority Transportation Expenditure Allocation Model

    Energy Science and Technology Software Center (OSTI)

    1993-04-12

    MITRAM (Minority TRansportation expenditure Allocation Model) can project various transportation related attributes of minority (Black and Hispanic) and majority (white) populations. The model projects vehicle ownership, vehicle miles of travel, workers, new car and on-road fleet fuel economy, amount and share of household income spent on gasoline, and household expenditures on public transportation and taxis. MITRAM predicts reactions to sustained fuel price changes for up to 10 years after the change.

  3. Reactive Transport and Coupled THM Processes in Engineering Barrier...

    Office of Scientific and Technical Information (OSTI)

    This may be important for identifying future modeling activities within the EBS group with ... Section 4 presents the results of geomechanical modeling using the Barcelona Basic Model ...

  4. Addressing numerical challenges in introducing a reactive transport code into a land surface model: A biogeochemical modeling proof-of-concept with CLM-PFLOTRAN 1.0: Modeling Archive

    DOE Data Explorer [Office of Scientific and Technical Information (OSTI)]

    R. T. Mills; G. Bisht; G. E. Hammond; J. Kumar; P. C. Lichtner; F. M. Hoffman; X. Xu; F. Yuan; S. L. Painter; G. Tang; and P.E. Thornton; B. Andre

    2016-04-19

    This Modeling Archive is in support of an NGEE Arctic discussion paper under review and available at doi:10.5194/gmd-9-927-2016. The purpose is to document the simulations to allow verification, reproducibility, and follow-up studies. This dataset contains shell scripts to create the CLM-PFLOTRAN cases, specific input files for PFLOTRAN and CLM, outputs, and python scripts to make the figures using the outputs in the publication. Through these results, we demonstrate that CLM-PFLOTRAN can approximately reproduce CLM results in selected cases for the Arctic, temperate and tropic sites. In addition, the new framework facilitates mechanistic representations of soil biogeochemistry processes in the land surface model.

  5. Uranium Transport Modeling

    SciTech Connect (OSTI)

    Bostick, William D.

    2008-01-15

    Uranium contamination is prevalent at many of the U.S. DOE facilities and at several civilian sites that have supported the nuclear fuel cycle. The potential off-site mobility of uranium depends on the partitioning of uranium between aqueous and solid (soil and sediment) phases. Hexavalent U (as uranyl, UO{sub 2}{sup 2+}) is relatively mobile, forming strong complexes with ubiquitous carbonate ion which renders it appreciably soluble even under mild reducing conditions. In the presence of carbonate, partition of uranyl to ferri-hydrate and select other mineral phases is usually maximum in the near-neutral pH range {approx} 5-8. The surface complexation reaction of uranyl with iron-containing minerals has been used as one means to model subsurface migration, used in conjunction with information on the site water chemistry and hydrology. Partitioning of uranium is often studied by short-term batch 'equilibrium' or long-term soil column testing ; MCLinc has performed both of these methodologies, with selection of method depending upon the requirements of the client or regulatory authority. Speciation of uranium in soil may be determined directly by instrumental techniques (e.g., x-ray photoelectron spectroscopy, XPS; x-ray diffraction, XRD; etc.) or by inference drawn from operational estimates. Often, the technique of choice for evaluating low-level radionuclide partitioning in soils and sediments is the sequential extraction approach. This methodology applies operationally-defined chemical treatments to selectively dissolve specific classes of macro-scale soil or sediment components. These methods recognize that total soil metal inventory is of limited use in understanding bioavailability or metal mobility, and that it is useful to estimate the amount of metal present in different solid-phase forms. Despite some drawbacks, the sequential extraction method can provide a valuable tool to distinguish among trace element fractions of different solubility related to

  6. Air Transport Optimization Model | NISAC

    Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)

    NISACAir Transport Optimization Model content top Network Optimization Models (RNAS and ATOM) Posted by Admin on Mar 1, 2012 in | Comments 0 comments Many critical infrastructures can be represented by a network of interconnected nodes and links. Mathematically sound nonlinear optimization techniques can then be applied to these networks to understand their behavior under normal and disrupted situations. Network optimization models are particularly useful for evaluating transportation system

  7. transportation-system-modeling-webinar

    Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)

    Webinar Announcement Webinar for the Intelligent Transportation Society of the Midwest (ITS Midwest) May 16, 2011 1:00 PM(CST) Hubert Ley Director, TRACC Argonne National Laboratory Argonne, Illinois High Performance Computing in Transportation Research - High Fidelity Transportation Models and More The Role of High-Performance Computing Because ITS relies on a very diverse collection of technologies, including communication and control technologies, advanced computing, information management

  8. Transport Modeling Working Group Meeting Reports | Department...

    Office of Energy Efficiency and Renewable Energy (EERE) Indexed Site

    Modeling Working Group Meeting Reports Transport Modeling Working Group Meeting Reports View reports from meetings of the Transport Modeling Working Group, which meets twice per ...

  9. Release of aged contaminants from weathered sediments: Effects of sorbate speciation on scaling of reactive transport

    SciTech Connect (OSTI)

    Chorover, Jon; Perdrial, Nico; Mueller, Karl; Strepka, Caleb; O???¢????????Day, Peggy; Rivera, Nelson; Um, Wooyong; Chang, Hyun-Shik; Steefel, Carl; Thompson, Aaron

    2012-11-05

    Hanford sediments impacted by hyperalkaline high level radioactive waste have undergone incongruent silicate mineral weathering concurrent with contaminant uptake. In this project, we studied the impact of background pore water (BPW) on strontium, cesium and iodine desorption and transport in Hanford sediments that were experimentally weathered by contact with simulated hyperalkaline tank waste leachate (STWL) solutions. Using those lab-weathered Hanford sediments (HS) and model precipitates formed during nucleation from homogeneous STWL solutions (HN), we (i) provided thorough characterization of reaction products over a matrix of field-relevant gradients in contaminant concentration, partial pressure of carbon dioxide, and reaction time; (ii) improved molecular-scale understanding of how sorbate speciation controls contaminant desorption from weathered sediments upon removal of caustic sources; and (iii) developed a mechanistic, predictive model of meso- to field-scale contaminant reactive transport under these conditions. In this final report, we provide detailed descriptions of our results from this three-year study, completed in 2012 following a one-year no cost extension.

  10. A Reactive Transport Simulator for Biogeochemical Processes in Subsurface System

    Energy Science and Technology Software Center (OSTI)

    2003-04-01

    BIOGEOCHEM is a Fortran code that mumerically simulates the coupled processes of solute transport, microbial population dynamics, microbial metabolism, and geochemical reactions. The potential applications of the code include, but not limited to, (a) sensitivity and uncertainty analyses for assessing the impact of microbial activity on subsurface geochemical systems; (b) extraction of biogeochemical parameter values from field observations or laboratory measurements, (c) helping to design and optimize laboratory biogeochemical experiments, and (d) data integration. Methodmore » of Solution: A finite difference method and a Newton-Raphson technique are used to solve a set of coupled nonlinear partial differential equations and algebraic equations. Practical Application: Environmental analysis, bioremediation performance assessments of radioactive or non-radioactive wase disposal, and academic research.« less

  11. COMBINED MODELING OF ACCELERATION, TRANSPORT, AND HYDRODYNAMIC...

    Office of Scientific and Technical Information (OSTI)

    Title: COMBINED MODELING OF ACCELERATION, TRANSPORT, AND HYDRODYNAMIC RESPONSE IN SOLAR FLARES. I. THE NUMERICAL MODEL Acceleration and transport of high-energy particles and fluid ...

  12. Reactive transport of uranium in fractured crystalline rock: Upscaling in time and distance

    SciTech Connect (OSTI)

    Dittrich, Timothy M.; Reimus, Paul W.

    2015-09-29

    In this study, batch adsorption and breakthrough column experiments were conducted to evaluate uranium transport through altered material that fills fractures in a granite rock system at the Grimsel Test Site in Switzerland at pH 6.9 and 7.9. The role of adsorption and desorption kinetics was evaluated with reactive transport modeling by comparing one-, two-, and three-site models. Emphasis was placed on describing long desorption tails that are important for upscaling in time and distance. The effect of increasing pH in injection solutions was also evaluated. For pH 6.9, a three-site model with forward rate constants between 0.07 and 0.8 ml g–1 h–1, reverse rate constants between 0.001 and 0.06 h–1, and site densities of 1.3, 0.104, and 0.026 μmol g–1 for ‘weak/fast’, ‘strong/slow’, and ‘very strong/very slow’ sites provided the best fits. For pH 7.9, a three-site model with forward rate constants between 0.05 and 0.8 mL g–1 h–1, reverse rate constants between 0.001 and 0.6 h–1, and site densities of 1.3, 0.039, and 0.013 μmol g–1 for a ‘weak/fast’, ‘strong/slow’, and ‘very strong/very slow’ sites provided the best fits. Column retardation coefficients (Rd) were 80 for pH 6.9 and 10.3 for pH 7.9. Model parameters determined from the batch and column experiments were used in 50 year large-scale simulations for continuous and pulse injections and indicated that a three-site model is necessary at pH 6.9, although a Kd-type equilibrium partition model with one-site was adequate for large scale predictions at pH 7.9. Batch experiments were useful for predicting early breakthrough times in the columns while column experiments helped differentiate the relative importance of sorption sites and desorption rate constants on transport.

  13. Reactive transport of uranium in fractured crystalline rock: Upscaling in time and distance

    DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)

    Dittrich, Timothy M.; Reimus, Paul W.

    2015-09-29

    In this study, batch adsorption and breakthrough column experiments were conducted to evaluate uranium transport through altered material that fills fractures in a granite rock system at the Grimsel Test Site in Switzerland at pH 6.9 and 7.9. The role of adsorption and desorption kinetics was evaluated with reactive transport modeling by comparing one-, two-, and three-site models. Emphasis was placed on describing long desorption tails that are important for upscaling in time and distance. The effect of increasing pH in injection solutions was also evaluated. For pH 6.9, a three-site model with forward rate constants between 0.07 and 0.8more » ml g–1 h–1, reverse rate constants between 0.001 and 0.06 h–1, and site densities of 1.3, 0.104, and 0.026 μmol g–1 for ‘weak/fast’, ‘strong/slow’, and ‘very strong/very slow’ sites provided the best fits. For pH 7.9, a three-site model with forward rate constants between 0.05 and 0.8 mL g–1 h–1, reverse rate constants between 0.001 and 0.6 h–1, and site densities of 1.3, 0.039, and 0.013 μmol g–1 for a ‘weak/fast’, ‘strong/slow’, and ‘very strong/very slow’ sites provided the best fits. Column retardation coefficients (Rd) were 80 for pH 6.9 and 10.3 for pH 7.9. Model parameters determined from the batch and column experiments were used in 50 year large-scale simulations for continuous and pulse injections and indicated that a three-site model is necessary at pH 6.9, although a Kd-type equilibrium partition model with one-site was adequate for large scale predictions at pH 7.9. Batch experiments were useful for predicting early breakthrough times in the columns while column experiments helped differentiate the relative importance of sorption sites and desorption rate constants on transport.« less

  14. transportation-systems-modeling-training

    Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)

    Training Table of Contents Date Location Integrated Transportation Models Workshop at ITM 2012 April 29, 2012 Hyatt Regency Tampa, FL TRANSIMS Training Course April 14-15, 2011 James E. Clyburn University Transportation Center Orangeburg, SC TRANSIMS RTSTEP Guest Lecturer March 29, 2011 Argonne TRACC Argonne, IL TRANSIMS Training Course January 19-21 2011 Argonne TRACC Argonne, IL TRANSIMS Training Course September 7-8, 2010 Turner Fairbank Highway Research Center Washington D.C. Network

  15. Predictive modeling of reactive wetting and metal joining.

    SciTech Connect (OSTI)

    van Swol, Frank B.

    2013-09-01

    The performance, reproducibility and reliability of metal joints are complex functions of the detailed history of physical processes involved in their creation. Prediction and control of these processes constitutes an intrinsically challenging multi-physics problem involving heating and melting a metal alloy and reactive wetting. Understanding this process requires coupling strong molecularscale chemistry at the interface with microscopic (diffusion) and macroscopic mass transport (flow) inside the liquid followed by subsequent cooling and solidification of the new metal mixture. The final joint displays compositional heterogeneity and its resulting microstructure largely determines the success or failure of the entire component. At present there exists no computational tool at Sandia that can predict the formation and success of a braze joint, as current capabilities lack the ability to capture surface/interface reactions and their effect on interface properties. This situation precludes us from implementing a proactive strategy to deal with joining problems. Here, we describe what is needed to arrive at a predictive modeling and simulation capability for multicomponent metals with complicated phase diagrams for melting and solidification, incorporating dissolutive and composition-dependent wetting.

  16. Final Project Report: Release of aged contaminants from weathered sediments: Effects of sorbate speciation on scaling of reactive transport

    SciTech Connect (OSTI)

    Jon Chorover, University of Arizona; Peggy O'????Day, University of California, Merced; Karl Mueller, Penn State University; Wooyong Um, Pacific Northwest National Laboratory; Carl Steefel, Lawrence Berkeley National Laboratory

    2012-10-01

    Hanford sediments impacted by hyperalkaline high level radioactive waste have undergone incongruent silicate mineral weathering concurrent with contaminant uptake. In this project, we studied the impact of background pore water (BPW) on strontium, cesium and iodine desorption and transport in Hanford sediments that were experimentally weathered by contact with simulated hyperalkaline tank waste leachate (STWL) solutions. Using those lab-weathered Hanford sediments (HS) and model precipitates formed during nucleation from homogeneous STWL solutions (HN), we (i) provided detailed characterization of reaction products over a matrix of field-relevant gradients in contaminant concentration, PCO2, and reaction time; (ii) improved molecular-scale understanding of how sorbate speciation controls contaminant desorption from weathered sediments upon removal of caustic sources; and (iii) developed a mechanistic, predictive model of meso- to field-scale contaminant reactive transport under these conditions.

  17. Semianalytical Solutions of Radioactive or Reactive Tracer Transport in Layered Fractured Media

    SciTech Connect (OSTI)

    G.J. Moridis; G. S. Bodvarsson

    2001-10-01

    In this paper, semianalytical solutions are developed for the problem of transport of radioactive or reactive tracers (solutes or colloids) through a layered system of heterogeneous fractured media with misaligned fractures. The tracer transport equations in the matrix account for (a) diffusion, (b) surface diffusion (for solutes only), (c) mass transfer between the mobile and immobile water fractions, (d) linear kinetic or equilibrium physical, chemical, or combined solute sorption or colloid filtration, and (e) radioactive decay or first order chemical reactions. Any number of radioactive decay daughter products (or products of a linear, first-order reaction chain) can be tracked. The tracer-transport equations in the fractures account for the same processes, in addition to advection and hydrodynamic dispersion. Additionally, the colloid transport equations account for straining and velocity adjustments related to the colloidal size. The solutions, which are analytical in the Laplace space, are numerically inverted to provide the solution in time and can accommodate any number of fractured and/or porous layers. The solutions are verified using analytical solutions for limiting cases of solute and colloid transport through fractured and porous media. The effect of important parameters on the transport of {sup 3}H, {sup 237}Np and {sup 239}Pu (and its daughters) is investigated in several test problems involving layered geological systems of varying complexity. {sup 239}Pu colloid transport problems in multilayered systems indicate significant colloid accumulations at straining interfaces but much faster transport of the colloid than the corresponding strongly sorbing solute species.

  18. A Computer Program for the Calculation of Reactivity and Kinetic Parameters by Two-Dimensional Neutron Transport Perturbation Theory.

    Energy Science and Technology Software Center (OSTI)

    1985-02-01

    Version 00 TP2 is a transport theory code, developed to determine reactivity effects and kinetic parameters such as effective delayed neutron fractions and mean generation time by applying the usual perturbation formalism for two-dimensional geometry.

  19. A Computer Code System for the Calculation of Reactivity and Kinetic Parameters by One-Dimensional Neutron Transport Perturbation Theory.

    Energy Science and Technology Software Center (OSTI)

    1985-02-01

    Version 00 TP1 is a transport theory code, developed to determine reactivity effects and kinetic parameters such as effective delayed neutron fractions and mean generation time by applying the usual perturbation formalism for one-dimensional geometry.

  20. Model Recovery Procedure for Response to a Radiological Transportation...

    Office of Environmental Management (EM)

    for Response to a Radiological Transportation Incident Model Recovery Procedure for Response to a Radiological Transportation Incident This Transportation Emergency...

  1. A Computer Program for the Reactivity and Kinetic Parameters for Two-Dimensional Triangular Geometry by Transport Perturbation Theory.

    Energy Science and Technology Software Center (OSTI)

    1990-04-25

    Version 00 TPTRIA calculates reactivity, effective delayed neutron fractions and mean generation time for two-dimensional triangular geometry on the basis of neutron transport perturbation theory. DIAMANT2 (also designated as CCC-414), is a multigroup two-dimensional discrete ordinates transport code system for triangular and hexagonal geometry which calculates direct and adjoint angular fluxes.

  2. RADIONUCLIDE TRANSPORT MODELS UNDER AMBIENT CONDITIONS

    SciTech Connect (OSTI)

    S. Magnuson

    2004-11-01

    The purpose of this model report is to document the unsaturated zone (UZ) radionuclide transport model, which evaluates, by means of three-dimensional numerical models, the transport of radioactive solutes and colloids in the UZ, under ambient conditions, from the repository horizon to the water table at Yucca Mountain, Nevada.

  3. New NIR Calibration Models Speed Biomass Composition and Reactivity Characterization

    SciTech Connect (OSTI)

    2015-09-01

    Obtaining accurate chemical composition and reactivity (measures of carbohydrate release and yield) information for biomass feedstocks in a timely manner is necessary for the commercialization of biofuels. This highlight describes NREL's work to use near-infrared (NIR) spectroscopy and partial least squares multivariate analysis to develop calibration models to predict the feedstock composition and the release and yield of soluble carbohydrates generated by a bench-scale dilute acid pretreatment and enzymatic hydrolysis assay. This highlight is being developed for the September 2015 Alliance S&T Board meeting.

  4. Technical Basis for Peak Reactivity Burnup Credit for BWR Spent Nuclear Fuel in Storage and Transportation Systems

    SciTech Connect (OSTI)

    Marshall, William BJ J; Ade, Brian J; Bowman, Stephen M; Gauld, Ian C; Ilas, Germina; Mertyurek, Ugur; Radulescu, Georgeta

    2015-01-01

    Oak Ridge National Laboratory and the United States Nuclear Regulatory Commission have initiated a multiyear project to investigate application of burnup credit for boiling-water reactor (BWR) fuel in storage and transportation casks. This project includes two phases. The first phase (1) investigates applicability of peak reactivity methods currently used in spent fuel pools (SFPs) to storage and transportation systems and (2) evaluates validation of both reactivity (keff) calculations and burnup credit nuclide concentrations within these methods. The second phase will focus on extending burnup credit beyond peak reactivity. This paper documents the first phase, including an analysis of lattice design parameters and depletion effects, as well as both validation components. Initial efforts related to extended burnup credit are discussed in a companion paper. Peak reactivity analyses have been used in criticality analyses for licensing of BWR fuel in SFPs over the last 20 years. These analyses typically combine credit for the gadolinium burnable absorber present in the fuel with a modest amount of burnup credit. Gadolinium burnable absorbers are used in BWR assemblies to control core reactivity. The burnable absorber significantly reduces assembly reactivity at beginning of life, potentially leading to significant increases in assembly reactivity for burnups less than 15–20 GWd/MTU. The reactivity of each fuel lattice is dependent on gadolinium loading. The number of gadolinium-bearing fuel pins lowers initial lattice reactivity, but it has a small impact on the burnup and reactivity of the peak. The gadolinium concentration in each pin has a small impact on initial lattice reactivity but a significant effect on the reactivity of the peak and the burnup at which the peak occurs. The importance of the lattice parameters and depletion conditions are primarily determined by their impact on the gadolinium depletion. Criticality code validation for BWR burnup

  5. Air pollution transport modeling. Master's thesis

    SciTech Connect (OSTI)

    Paal, D.M.

    1993-12-01

    This research effort addresses modeling of the transportation of air pollution in the atmosphere and the numerical analysis of the partial differential equations used in such modeling. Three Gaussian models are examined and compared using example problems. Several finite difference schemes are developed to solve the partial differential equations used in air pollution transport modeling. This study examines three Gaussian models: SCREEN, AFTOX, and the program GAUSPLUM. The model GAUSPLUM is developed in this study and uses the Ada programming language and the analytic solution to the advection-diffusion equation. Numerical analysis of the partial differential equations (PDE) used in air pollution modeling is also examined. The equations are generally parabolic or hyperbolic PDE's. The following are examined in this research: the advection equation; the one-, two-, and three-dimensional advection-diffusion equations; and the two-dimensional steady-state equation. Air Pollution Transport, Modeling, Finite Difference Scheme, Stability, Consistency, Convergence, Advection-Diffusion Equations.

  6. Developing a robust geochemical and reactive transport model...

    Office of Scientific and Technical Information (OSTI)

    Migration of carbon dioxide (COsub 2) from deep storage formations into shallow drinking ... LEAKS; LEVELS; MATERIALS; METALS; MIGRATION; MINERALS; MOBILITY; OXIDES; ...

  7. Modeling of coupled heat transfer and reactive transport processesin...

    Office of Scientific and Technical Information (OSTI)

    heating and boiling, and through local convection. In cooler regions, the vapor condenses on fracture walls, where it drains through the fracture network. Slow imbibition of water ...

  8. Modeling of coupled heat transfer and reactive transport processesin...

    Office of Scientific and Technical Information (OSTI)

    Subsequently, vapor migrates out of the matrix pore space, moving away from the tunnel through the permeable fracture network. This migration is propelled by buoyancy, by the ...

  9. Vertical transport and sources in flux models

    SciTech Connect (OSTI)

    Canavan, G.H.

    1997-01-01

    Vertical transport in flux models in examined and shown to reproduce expected limits for densities and fluxes. Disparities with catalog distributions are derived and inverted to find the sources required to rectify them.

  10. The Transportation Sector Model of the National Energy Modeling...

    Gasoline and Diesel Fuel Update (EIA)

    Information Administration NEMS Transportation Demand Model Documentation Report 2005 25 manufacturing, and design advances. Manufacturing advances can generally be thought of as...

  11. Pore-scale simulation of coupled reactive transport and dissolution in fractures and porous media using the level set interface tracking method

    SciTech Connect (OSTI)

    Hai Huang; Xiaoyi Li

    2011-01-01

    A level set simulation methodology developed for modeling coupled reactive transport and structure evolution has been applied to dissolution in fracture apertures and porous media. The coupled processes such as fluid flow, reactant transport and dissolution at the solid-liquid interfaces are handled simultaneously. The reaction-induced evolution of solid-liquid interfaces is captured using the level set method, with the advantage of representing the interface with sub-grid scale resolution. The coupled processes are simulated for several geometric models of fractures and porous media under various flow conditions and reaction rates. Quantitative relationships between permeability and porosity are obtained from some of the simulation results and compared with analytical constitutive relations (i.e., the conventional cubic law and the Carman-Kozeny law) based on simplified pore space geometries and reaction induced geometric evolutions. The drastic deviation of the simulation results from these analytical theories is explained by the development of large local concentration gradients of reactants within fracture apertures and individual pores observed in the simulation results and consequently the complex geometric evolution patterns of fracture apertures and pores due to mineral dissolution. The simulation results support the argument that traditional constitutive relations based on simplified geometries and conditions have limited applicability in predicting field scale reactive transport and that incorporation of micro-scale physics is necessary.

  12. Reactive sputter deposition of pyrite structure transition metal disulfide thin films: Microstructure, transport, and magnetism

    SciTech Connect (OSTI)

    Baruth, A.; Manno, M.; Narasimhan, D.; Shankar, A.; Zhang, X.; Johnson, M.; Aydil, E. S.; Leighton, C. [Department of Chemical Engineering and Materials Science, University of Minnesota, Minneapolis, Minnesota 55455 (United States)

    2012-09-01

    Transition metal disulfides crystallizing in the pyrite structure (e.g., TMS{sub 2}, with TM = Fe, Co, Ni, and Cu) are a class of materials that display a remarkably diverse array of functional properties. These properties include highly spin-polarized ferromagnetism (in Co{sub 1-x}Fe{sub x}S{sub 2}), superconductivity (in CuS{sub 2}), an antiferromagnetic Mott insulating ground state (in NiS{sub 2}), and semiconduction with close to optimal parameters for solar absorber applications (in FeS{sub 2}). Exploitation of these properties in heterostructured devices requires the development of reliable and reproducible methods for the deposition of high quality pyrite structure thin films. In this manuscript, we report on the suitability of reactive sputter deposition from metallic targets in an Ar/H{sub 2}S environment as a method to achieve exactly this. Optimization of deposition temperature, Ar/H{sub 2}S pressure ratio, and total working gas pressure, assisted by plasma optical emission spectroscopy, reveals significant windows over which deposition of single-phase, polycrystalline, low roughness pyrite films can be achieved. This is illustrated for the test cases of the ferromagnetic metal CoS{sub 2} and the diamagnetic semiconductor FeS{sub 2}, for which detailed magnetic and transport characterization are provided. The results indicate significant improvements over alternative deposition techniques such as ex situ sulfidation of metal films, opening up exciting possibilities for all-sulfide heterostructured devices. In particular, in the FeS{sub 2} case it is suggested that fine-tuning of the sputtering conditions provides a potential means to manipulate doping levels and conduction mechanisms, critical issues in solar cell applications. Parenthetically, we note that conditions for synthesis of phase-pure monosulfides and thiospinels are also identified.

  13. World Energy Projection System Plus Model Documentation: Transportation Model

    Reports and Publications (EIA)

    2011-01-01

    This report documents the objectives, analytical approach and development of the World Energy Projection System Plus (WEPS ) International Transportation model. It also catalogues and describes critical assumptions, computational methodology, parameter estimation techniques, and model source code.

  14. Radionuclide Transport Models Under Ambient Conditions

    SciTech Connect (OSTI)

    G. Moridis; Q. Hu

    2001-12-20

    The purpose of Revision 00 of this Analysis/Model Report (AMR) is to evaluate (by means of 2-D semianalytical and 3-D numerical models) the transport of radioactive solutes and colloids in the unsaturated zone (UZ) under ambient conditions from the potential repository horizon to the water table at Yucca Mountain (YM), Nevada.

  15. Radionuclide Transport Models Under Ambient Conditions

    SciTech Connect (OSTI)

    G. Moridis; Q. Hu

    2000-03-12

    The purpose of this Analysis/Model Report (AMR) is to evaluate (by means of 2-D semianalytical and 3-D numerical models) the transport of radioactive solutes and colloids in the unsaturated zone (UZ) under ambient conditions from the potential repository horizon to the water table at Yucca Mountain (YM), Nevada. This is in accordance with the ''AMR Development Plan U0060, Radionuclide Transport Models Under Ambient Conditions'' (CRWMS M and O 1999a). This AMR supports the UZ Flow and Transport Process Model Report (PMR). This AMR documents the UZ Radionuclide Transport Model (RTM). This model considers: the transport of radionuclides through fractured tuffs; the effects of changes in the intensity and configuration of fracturing from hydrogeologic unit to unit; colloid transport; physical and retardation processes and the effects of perched water. In this AMR they document the capabilities of the UZ RTM, which can describe flow (saturated and/or unsaturated) and transport, and accounts for (a) advection, (b) molecular diffusion, (c) hydrodynamic dispersion (with full 3-D tensorial representation), (d) kinetic or equilibrium physical and/or chemical sorption (linear, Langmuir, Freundlich or combined), (e) first-order linear chemical reaction, (f) radioactive decay and tracking of daughters, (g) colloid filtration (equilibrium, kinetic or combined), and (h) colloid-assisted solute transport. Simulations of transport of radioactive solutes and colloids (incorporating the processes described above) from the repository horizon to the water table are performed to support model development and support studies for Performance Assessment (PA). The input files for these simulations include transport parameters obtained from other AMRs (i.e., CRWMS M and O 1999d, e, f, g, h; 2000a, b, c, d). When not available, the parameter values used are obtained from the literature. The results of the simulations are used to evaluate the transport of radioactive solutes and colloids, and

  16. PFLOTRAN User Manual: A Massively Parallel Reactive Flow and...

    Office of Scientific and Technical Information (OSTI)

    Technical Report: PFLOTRAN User Manual: A Massively Parallel Reactive Flow and Transport Model for Describing Surface and Subsurface Processes Citation Details In-Document Search...

  17. PFLOTRAN User Manual: A Massively Parallel Reactive Flow and...

    Office of Scientific and Technical Information (OSTI)

    PFLOTRAN User Manual: A Massively Parallel Reactive Flow and Transport Model for Describing Surface and Subsurface Processes Lichtner, Peter OFM Research; Karra, Satish Los...

  18. Current Transportation Models Used in the Vehicle Technologies Program

    SciTech Connect (OSTI)

    2009-04-06

    A summary of various transportation models (VISION, TRUCK, GREET, Oil Peaking Model, Feebate Model, Oil Security Metrics Model, ORNL PHEV Choice Model: Version 1, PSAT, PSAT-PRO,

  19. Performance Assessment Transport Modeling of Uranium at the Area 5 Radioactive Waste Management Site at the Nevada National Security Site

    SciTech Connect (OSTI)

    NSTec Radioactive Waste

    2010-10-12

    Following is a brief summary of the assumptions that are pertinent to the radioactive isotope transport in the GoldSim Performance Assessment model of the Area 5 Radioactive Waste Management Site, with special emphasis on the water-phase reactive transport of uranium, which includes depleted uranium products.

  20. A Many-Task Parallel Approach for Multiscale Simulations of Subsurface Flow and Reactive Transport

    SciTech Connect (OSTI)

    Scheibe, Timothy D.; Yang, Xiaofan; Schuchardt, Karen L.; Agarwal, Khushbu; Chase, Jared M.; Palmer, Bruce J.; Tartakovsky, Alexandre M.

    2014-12-16

    Continuum-scale models have long been used to study subsurface flow, transport, and reactions but lack the ability to resolve processes that are governed by pore-scale mixing. Recently, pore-scale models, which explicitly resolve individual pores and soil grains, have been developed to more accurately model pore-scale phenomena, particularly reaction processes that are controlled by local mixing. However, pore-scale models are prohibitively expensive for modeling application-scale domains. This motivates the use of a hybrid multiscale approach in which continuum- and pore-scale codes are coupled either hierarchically or concurrently within an overall simulation domain (time and space). This approach is naturally suited to an adaptive, loosely-coupled many-task methodology with three potential levels of concurrency. Each individual code (pore- and continuum-scale) can be implemented in parallel; multiple semi-independent instances of the pore-scale code are required at each time step providing a second level of concurrency; and Monte Carlo simulations of the overall system to represent uncertainty in material property distributions provide a third level of concurrency. We have developed a hybrid multiscale model of a mixing-controlled reaction in a porous medium wherein the reaction occurs only over a limited portion of the domain. Loose, minimally-invasive coupling of pre-existing parallel continuum- and pore-scale codes has been accomplished by an adaptive script-based workflow implemented in the Swift workflow system. We describe here the methods used to create the model system, adaptively control multiple coupled instances of pore- and continuum-scale simulations, and maximize the scalability of the overall system. We present results of numerical experiments conducted on NERSC supercomputing systems; our results demonstrate that loose many-task coupling provides a scalable solution for multiscale subsurface simulations with minimal overhead.

  1. Climate Impact of Transportation A Model Comparison

    SciTech Connect (OSTI)

    Girod, Bastien; Van Vuuren, Detlef; Grahn, Maria; Kitous, Alban; Kim, Son H.; Kyle, G. Page

    2013-06-01

    Transportation contributes to a significant and rising share of global energy use and GHG emissions. Therefore modeling future travel demand, its fuel use, and resulting CO2 emission is highly relevant for climate change mitigation. In this study we compare the baseline projections for global service demand (passenger-kilometers, ton-kilometers), fuel use, and CO2 emissions of five different global transport models using harmonized input assumptions on income and population. For four models we also evaluate the impact of a carbon tax. All models project a steep increase in service demand over the century. Technology is important for limiting energy consumption and CO2 emissions, but quite radical changes in the technology mix are required to stabilize or reverse the trend. While all models project liquid fossil fuels dominating up to 2050, they differ regarding the use of alternative fuels (natural gas, hydrogen, biofuels, and electricity), because of different fuel price projections. The carbon tax of US$200/tCO2 in 2050 stabilizes or reverses global emission growth in all models. Besides common findings many differences in the model assumptions and projections indicate room for improvement in modeling and empirical description of the transport system.

  2. Model for assessing bronchial mucus transport

    SciTech Connect (OSTI)

    Agnew, J.E.; Bateman, J.R.M.; Pavia, D.; Clarke, S.W.

    1984-02-01

    The authors propose a scheme for the assessment of regional mucus transport using inhaled Tc-99m aerosol particles and quantitative analysis of serial gamma-camera images. The model treats input to inner and intermediate lung regions as the total of initial deposition there plus subsequent transport into these regions from more peripheral airways. It allows for interregional differences in the proportion of particles deposited on the mucus-bearing conducting airways, and does not require a gamma image 24 hr after particle inhalation. Instead, distribution of particles reaching the respiratory bronchioles or alveoli is determined from a Kr-81m ventilation image, while the total amount of such deposition is obtained from 24-hr Tc-99m retention measured with a sensitive counter system. The model is applicable to transport by mucociliary action or by cough, and has been tested in ten normal and ten asthmatic subjects.

  3. Modeling heat conduction and radiation transport with the diffusion...

    Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)

    heat conduction and radiation transport with the diffusion equation in NIF ALE-AMR This ... IOPscience Modeling Heat Conduction and Radiation Transport with the Diffusion Equation in ...

  4. Model Recovery Procedure for Response to a Radiological Transportation Incident

    Broader source: Energy.gov [DOE]

    This Transportation Emergency Preparedness Program (TEPP) Model Recovery Procedure contains the recommended elements for developing and conducting recovery planning at transportation incident scene...

  5. Fire and materials modeling for transportation systems

    SciTech Connect (OSTI)

    Skocypec, R.D.; Gritzo, L.A.; Moya, J.L.; Nicolette, V.F.; Tieszen, S.R.; Thomas, R.

    1994-10-01

    Fire is an important threat to the safety of transportation systems. Therefore, understanding the effects of fire (and its interaction with materials) on transportation systems is crucial to quantifying and mitigating the impact of fire on the safety of those systems. Research and development directed toward improving the fire safety of transportation systems must address a broad range of phenomena and technologies, including: crash dynamics, fuel dispersion, fire environment characterization, material characterization, and system/cargo thermal response modeling. In addition, if the goal of the work is an assessment and/or reduction of risk due to fires, probabilistic risk assessment technology is also required. The research currently underway at Sandia National Laboratories in each of these areas is summarized in this paper.

  6. Final Scientific/Technical Report – DE-FG02-06ER64172 – Reaction-Based Reactive Transport Modeling of Iron Reduction and Uranium Immobilization at Area 2 of the NABIR Field Research Center – Subproject to Co-PI Eric E. Roden

    SciTech Connect (OSTI)

    Eric E. Roden

    2009-03-17

    This report summarizes research conducted in conjunction with a project entitled “Reaction-Based Reactive Transport Modeling of Iron Reduction and Uranium Immobilization at Area 2 of the NABIR Field Research Center”, which was funded through the Integrative Studies Element of the former NABIR Program (now the Environmental Remediation Sciences Program) within the Office of Biological and Environmental Research. Dr. William Burgos (The Pennsylvania State University) was the overall PI/PD for the project, which included Brian Dempsey (Penn State), Gour-Tsyh (George) Yeh (Central Florida University), and Eric Roden (formerly at The University of Alabama, now at the University of Wisconsin) as separately-funded co-PIs. The project focused on development of a mechanistic understanding and quantitative models of coupled Fe(III)/U(VI) reduction in FRC Area 2 sediments. The work builds on our previous studies of microbial Fe(III) and U(VI) reduction, and was directly aligned with the Scheibe et al. ORNL FRC Field Project at Area 2. Area 2 is a shallow pathway for migration of contaminated groundwater to seeps in the upper reach of Bear Creek at ORNL, mainly through a ca. 1 m thick layer of gravel located 4-5 m below the ground surface. The gravel layer is sandwiched between an overlying layer of disturbed fill material, and 2-3 m of undisturbed shale saprolite derived from the underlying Nolichucky Shale bedrock. The fill was put in place when contaminated soils were excavated and replaced by native saprolite from an uncontaminated area within Bear Creek Valley; the gravel layer was presumably installed prior to addition of the fill in order to provide a stable surface for the operation of heavy machinery. The undisturbed saprolite is highly weathered bedrock that has unconsolidated character but retains much of the bedding and fracture structure of the parent rock (shale with interbedded limestone). Hydrological tracer studies conducted during the Scheibe et al. field

  7. Modeling uranium transport in acidic contaminated groundwater with base addition

    SciTech Connect (OSTI)

    Zhang, Fan; Luo, Wensui; Parker, Jack C.; Brooks, Scott C; Watson, David B; Jardine, Philip; Gu, Baohua

    2011-01-01

    This study investigates reactive transport modeling in a column of uranium(VI)-contaminated sediments with base additions in the circulating influent. The groundwater and sediment exhibit oxic conditions with low pH, high concentrations of NO{sub 3}{sup -}, SO{sub 4}{sup 2-}, U and various metal cations. Preliminary batch experiments indicate that additions of strong base induce rapid immobilization of U for this material. In the column experiment that is the focus of the present study, effluent groundwater was titrated with NaOH solution in an inflow reservoir before reinjection to gradually increase the solution pH in the column. An equilibrium hydrolysis, precipitation and ion exchange reaction model developed through simulation of the preliminary batch titration experiments predicted faster reduction of aqueous Al than observed in the column experiment. The model was therefore modified to consider reaction kinetics for the precipitation and dissolution processes which are the major mechanism for Al immobilization. The combined kinetic and equilibrium reaction model adequately described variations in pH, aqueous concentrations of metal cations (Al, Ca, Mg, Sr, Mn, Ni, Co), sulfate and U(VI). The experimental and modeling results indicate that U(VI) can be effectively sequestered with controlled base addition due to sorption by slowly precipitated Al with pH-dependent surface charge. The model may prove useful to predict field-scale U(VI) sequestration and remediation effectiveness.

  8. ENVIRONMENTAL REACTIVITY OF SOLID STATE HYDRIDE MATERIALS: MODELING AND TESTING FOR AIR AND WATER EXPOSURE

    SciTech Connect (OSTI)

    Anton, D.; James, C.; Cortes-Concepcion, J.; Tamburello, D.; Brinkman, K.; Gray, J.

    2010-05-18

    To make commercially acceptable condensed phase hydrogen storage systems, it is important to understand quantitatively the risks involved in using these materials. A rigorous set of environmental reactivity tests have been developed based on modified testing procedures codified by the United Nations for the transportation of dangerous goods. Potential hydrogen storage material, 2LiBH4{center_dot}MgH2 and NH3BH3, have been tested using these modified procedures to evaluate the relative risks of these materials coming in contact with the environment in hypothetical accident scenarios. It is apparent that an ignition event will only occur if both a flammable concentration of hydrogen and sufficient thermal energy were available to ignite the hydrogen gas mixture. In order to predict hydride behavior for hypothesized accident scenarios, an idealized finite element model was developed for dispersed hydride from a breached system. Empirical thermodynamic calculations based on precise calorimetric experiments were performed in order to quantify the energy and hydrogen release rates and to quantify the reaction products resulting from water and air exposure. Both thermal and compositional predictions were made with identification of potential ignition event scenarios.

  9. Modelling of multi-ion-beam reactive cosputtering for metal oxide thin films

    SciTech Connect (OSTI)

    Xiao, D.Q.; Zhu, J.G.; Qian, Z.H.; Peng, W.B.; Wei, L.F.; Li, Z.S.

    1995-12-31

    Very recently a new technique named multi-ion-beam reactive cosputtering (MIBRECS) was developed for preparing multi-component metal oxide thin films. Epitaxial or highly oriented (Pb, La) TiO{sub 3} thin films sputtered from pure metals of lead, titanium and lathanium were deposited by using this technique. In order to consummate the technique and to study the mechanism of reactive cosputtering, a general model of multi-ion-beam reactive cosputtering was proposed for the first time based on the well-known gas kinetics under stable sputtering circumstances, and a computer numerical simulation of the model was carried out with the parameters adopted in the experiments. The relationships among the sputtering ratios of the targets, and the coverage ratios of simple substances and oxides of the target metals on substrate surface with the total reactive gas flux and the densities of the sputtering ion beam were obtained respectively, and the hysteresis effect of the characteristic of reactive sputtering and the interactions during multi-ion-beam reactive cosputtering processes were also obtained. The numerical simulation results are at least qualitatively in agreement with the experiments.

  10. Documentation of TRU biological transport model (BIOTRAN)

    SciTech Connect (OSTI)

    Gallegos, A.F.; Garcia, B.J.; Sutton, C.M.

    1980-01-01

    Inclusive of Appendices, this document describes the purpose, rationale, construction, and operation of a biological transport model (BIOTRAN). This model is used to predict the flow of transuranic elements (TRU) through specified plant and animal environments using biomass as a vector. The appendices are: (A) Flows of moisture, biomass, and TRU; (B) Intermediate variables affecting flows; (C) Mnemonic equivalents (code) for variables; (D) Variable library (code); (E) BIOTRAN code (Fortran); (F) Plants simulated; (G) BIOTRAN code documentation; (H) Operating instructions for BIOTRAN code. The main text is presented with a specific format which uses a minimum of space, yet is adequate for tracking most relationships from their first appearance to their formulation in the code. Because relationships are treated individually in this manner, and rely heavily on Appendix material for understanding, it is advised that the reader familiarize himself with these materials before proceeding with the main text.

  11. Neutral gas transport modeling with DEGAS 2

    SciTech Connect (OSTI)

    Stotler, D.; Karney, C.

    1993-10-01

    We are currently rewriting the neutral gas transport code, DEGAS with a view to not only making it faster, but also easing the process of including new physics. The goal is to make adding new species and reactions relatively simple so that the code can be rapidly adapted to new divertor physics regimes. DEGAS 2 will also be optimized for coupling to fluid plasma codes, incorporating many of the techniques utilized in B2-EIRENE. Finally, it is our intention that DEGAS 2, like DEGAS, be well-documented and easy to use. We ill present model calculations including ionization and charge exchange which will illustrate the way reactions are included into DEGAS 2 and will demonstrate operation of the code on a distributed network of workstations.

  12. Diffusion Dominant Solute Transport Modelling in Fractured Media Under Deep Geological Environment - 12211

    SciTech Connect (OSTI)

    Kwong, S.; Jivkov, A.P.

    2012-07-01

    Deep geologic disposal of high activity and long-lived radioactive waste is gaining increasing support in many countries, where suitable low permeability geological formation in combination with engineered barriers are used to provide long term waste contaminant and minimise the impacts to the environment and risk to the biosphere. This modelling study examines the solute transport in fractured media under low flow velocities that are relevant to a deep geological environment. In particular, reactive solute transport through fractured media is studied using a 2-D model, that considers advection and diffusion, to explore the coupled effects of kinetic and equilibrium chemical processes. The effects of water velocity in the fracture, matrix porosity and diffusion on solute transport are investigated and discussed. Some illustrative modelled results are presented to demonstrate the use of the model to examine the effects of media degradation on solute transport, under the influences of hydrogeological (diffusion dominant) and microbially mediated chemical processes. The challenges facing the prediction of long term degradation such as cracks evolution, interaction and coalescence are highlighted. The potential of a novel microstructure informed modelling approach to account for these effects is discussed, particularly with respect to investigating multiple phenomena impact on material performance. The GRM code is used to examine the effects of media degradation for a geological waste disposal package, under the combined hydrogeological (diffusion dominant) and chemical effects in low groundwater flow conditions that are typical of deep geological disposal systems. An illustrative reactive transport modelling application demonstrates the use of the code to examine the interplay of kinetic controlled biogeochemical reactive processes with advective and diffusive transport, under the influence of media degradation. The initial model results are encouraging which show the

  13. On Numerical Considerations for Modeling Reactive Astrophysical Shocks

    SciTech Connect (OSTI)

    Papatheodore, Thomas L; Messer, Bronson

    2014-01-01

    Simulating detonations in astrophysical environments is often complicated by numerical approximations to shock structure. A common prescription to ensure correct detonation speeds (and associated quantities) is to prohibit burning inside the numerically broadened shock (Fryxell et al. 1989). We have performed a series of simulations to verify the efficacy of this approximation and to understand how resolution and dimensionality might affect its use. Our results show that, in one dimension, prohibiting burning in the shock is important wherever the carbon burning length is not resolved, in keeping with the results of Fryxell et al. (1989). In two dimensions, we find that the prohibition of shock burning effectively inhibits the development of cellular structure for all but the most highly-resolved cases. We discuss the possible impacts this outcome may have on sub-grid models and detonation propagation in Type Ia supernovae.

  14. An Updated Site Scale Saturated Zone Ground Water Transport Model...

    Office of Scientific and Technical Information (OSTI)

    An Updated Site Scale Saturated Zone Ground Water Transport Model for Yucca Mountain. Citation Details In-Document Search Title: An Updated Site Scale Saturated Zone Ground Water ...

  15. STOMP Subsurface Transport Over Multiple Phases Version 1.0 Addendum: ECKEChem Equilibrium-Conservation-Kinetic Equation Chemistry and Reactive Transport

    SciTech Connect (OSTI)

    White, Mark D.; McGrail, B. Peter

    2005-12-01

    flow and transport simulator, STOMP (Subsurface Transport Over Multiple Phases). Prior to these code development activities, the STOMP simulator included sequential and scalable implementations for numerically simulating the injection of supercritical CO2 into deep saline aquifers. Additionally, the sequential implementations included operational modes that considered nonisothermal conditions and kinetic dissolution of CO2 into the saline aqueous phase. This addendum documents the advancement of these numerical simulation capabilities to include reactive transport in the STOMP simulator through the inclusion of the recently PNNL developed batch geochemistry solution module ECKEChem (Equilibrium-Conservation-Kinetic Equation Chemistry). Potential geologic reservoirs for sequestering CO2 include deep saline aquifers, hydrate-bearing formations, depleted or partially depleted natural gas and petroleum reservoirs, and coal beds. The mechanisms for sequestering carbon dioxide in geologic reservoirs include physical trapping, dissolution in the reservoir fluids, hydraulic trapping (hysteretic entrapment of nonwetting fluids), and chemical reaction. This document and the associated code development and verification work are concerned with the chemistry of injecting CO2 into geologic reservoirs. As geologic sequestration of CO2 via chemical reaction, namely precipitation reactions, are most dominate in deep saline aquifers, the principal focus of this document is the numerical simulation of CO2 injection, migration, and geochemical reaction in deep saline aquifers. The ECKEChem batch chemistry module was developed in a fashion that would allow its implementation into all operational modes of the STOMP simulator, making it a more versatile chemistry component. Additionally, this approach allows for verification of the ECKEChem module against more classical reactive transport problems involving aqueous systems.

  16. Transport Modeling Working Group | Department of Energy

    Energy Savers [EERE]

    for and understanding of transport phenomena in polymer electrolyte fuel cell stacks. ... to the efficient and cost-effective operation of polymer electrolyte fuel cells. ...

  17. Understanding the Subsurface Reactive Transport of Transuranic Contaminants at DOE Sites

    SciTech Connect (OSTI)

    Barnett, Mark O.; Albrecht-Schmitt, Thomas E.; Saiers, James E.; Shuh, David K.

    2013-12-20

    Our primary hypothesis is that actinides can interact with surfaces in fundamentally different ways than other metals, metalloids, and oxyanions and that this fundamental difference requires new approaches to studying and modeling transuranic sorption to minerals and geomedia. This project supports a key mission of the SBR program to develop sufficient scientific understanding such that DOE sites will be able to incorporate coupled physical, chemical, and biological processes into decision making for environmental management and long-term stewardship, while also supporting DOE’s commitment to education, training, and collaboration with DOE user facilities.

  18. Transportation and Vehicle Energy Modeling | Argonne National Laboratory

    Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)

    Transportation and Vehicle Energy Modeling Transportation and Vehicle Energy Modeling For the past 15 years, Argonne has been at the forefront of research in energy-efficient transportation. In recent years, the vehicle technologies have become increasingly complex with the introduction of new powertrain configurations (such as electrified vehicles), new component technologies (such as advanced transmissions and engines) and control strategies (such eco-routing). In addition, with increased

  19. AEO2017 Modeling updates in the transportation sector

    U.S. Energy Information Administration (EIA) Indexed Site

    7 For AEO2017 Transportation Working Group August 31, 2016 | Washington, DC By Melissa Lynes, John Maples, Mark Schipper, and David Stone Office of Energy Consumption and Efficiency Analysis Modeling updates in the transportation sector Updates to the Annual Energy Outlook 2017 * Transportation demand model highlights - 10-year extension of last-year projection, AEO2016 is 2040 and AEO2017 is 2050 - Battery costs for electric vehicles - Phase 2 greenhouse gas and fuel efficiency standards for

  20. A quantum energy transport model for semiconductor device simulation

    SciTech Connect (OSTI)

    Sho, Shohiro; Odanaka, Shinji

    2013-02-15

    This paper describes numerical methods for a quantum energy transport (QET) model in semiconductors, which is derived by using a diffusion scaling in the quantum hydrodynamic (QHD) model. We newly drive a four-moments QET model similar with a classical ET model. Space discretization is performed by a new set of unknown variables. Numerical stability and convergence are obtained by developing numerical schemes and an iterative solution method with a relaxation method. Numerical simulations of electron transport in a scaled MOSFET device are discussed. The QET model allows simulations of quantum confinement transport, and nonlocal and hot-carrier effects in scaled MOSFETs.

  1. Water Transport in PEM Fuel Cells: Advanced Modeling, Material Selection,

    Office of Energy Efficiency and Renewable Energy (EERE) Indexed Site

    Testing and Design Optimization | Department of Energy in PEM Fuel Cells: Advanced Modeling, Material Selection, Testing and Design Optimization Water Transport in PEM Fuel Cells: Advanced Modeling, Material Selection, Testing and Design Optimization Part of a $100 million fuel cell award announced by DOE Secretary Bodman on Oct. 25, 2006. 2_cfd.pdf (21.58 KB) More Documents & Publications Fuel Cell Kickoff Meeting Agenda Water Transport in PEM Fuel Cells: Advanced Modeling, Material

  2. Large-Scale Transport Model Uncertainty and Sensitivity Analysis: Distributed Sources in Complex Hydrogeologic Systems

    SciTech Connect (OSTI)

    Sig Drellack, Lance Prothro

    2007-12-01

    The Underground Test Area (UGTA) Project of the U.S. Department of Energy, National Nuclear Security Administration Nevada Site Office is in the process of assessing and developing regulatory decision options based on modeling predictions of contaminant transport from underground testing of nuclear weapons at the Nevada Test Site (NTS). The UGTA Project is attempting to develop an effective modeling strategy that addresses and quantifies multiple components of uncertainty including natural variability, parameter uncertainty, conceptual/model uncertainty, and decision uncertainty in translating model results into regulatory requirements. The modeling task presents multiple unique challenges to the hydrological sciences as a result of the complex fractured and faulted hydrostratigraphy, the distributed locations of sources, the suite of reactive and non-reactive radionuclides, and uncertainty in conceptual models. Characterization of the hydrogeologic system is difficult and expensive because of deep groundwater in the arid desert setting and the large spatial setting of the NTS. Therefore, conceptual model uncertainty is partially addressed through the development of multiple alternative conceptual models of the hydrostratigraphic framework and multiple alternative models of recharge and discharge. Uncertainty in boundary conditions is assessed through development of alternative groundwater fluxes through multiple simulations using the regional groundwater flow model. Calibration of alternative models to heads and measured or inferred fluxes has not proven to provide clear measures of model quality. Therefore, model screening by comparison to independently-derived natural geochemical mixing targets through cluster analysis has also been invoked to evaluate differences between alternative conceptual models. Advancing multiple alternative flow models, sensitivity of transport predictions to parameter uncertainty is assessed through Monte Carlo simulations. The

  3. Simple Model Representations of Transport in a Complex Fracture...

    Office of Scientific and Technical Information (OSTI)

    Effects on Long-Term Predictions Citation Details In-Document Search Title: Simple Model Representations of Transport in a Complex Fracture and Their Effects on Long-Term ...

  4. CO2 leakage impacts on shallow groundwater. Field-scale reactive-transport simulations informed by observations at a natural analog site

    DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)

    Keating, Elizabeth H.; Hakala, J. Alexandra; Viswanathan, Hari; Carey, J. William; Pawar, Rajesh; Guthrie, George D.; Fessenden-Rahn, Julianna

    2013-03-01

    It is challenging to predict the degree to which shallow groundwater might be affected by leaks from a CO2 sequestration reservoir, particularly over long time scales and large spatial scales. In this study observations at a CO2 enriched shallow aquifer natural analog were used to develop a predictive model which is then used to simulate leakage scenarios. This natural analog provides the opportunity to make direct field observations of groundwater chemistry in the presence of elevated CO2, to collect aquifer samples and expose them to CO2 under controlled conditions in the laboratory, and to test the ability of multiphase reactivemore » transport models to reproduce measured geochemical trends at the field-scale. The field observations suggest that brackish water entrained with the upwelling CO2 are a more significant source of trace metals than in situ mobilization of metals due to exposure to CO2. The study focuses on a single trace metal of concern at this site: U. Experimental results indicate that cation exchange/adsorption and dissolution/precipitation of calcite containing trace amounts of U are important reactions controlling U in groundwater at this site, and that the amount of U associated with calcite is fairly well constrained. Simulations incorporating these results into a 3-D multi-phase reactive transport model are able to reproduce the measured ranges and trends between pH, pCO2, Ca, total C, U and Cl-at the field site. Although the true fluxes at the natural analog site are unknown, the cumulative CO2 flux inferred from these simulations are approximately equivalent to 37.8E-3 MT, approximately corresponding to a .001% leak rate for injection at a large (750 MW) power plant. The leakage scenario simulations suggest that if the leak only persists for a short time the volume of aquifer contaminated by CO2-induced mobilization of U will be relatively small, yet persistent over 100 a.« less

  5. Elucidating reactivity regimes in cyclopentane oxidation: Jet stirred reactor experiments, computational chemistry, and kinetic modeling

    DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)

    Al Rashidi, Mariam J.; Thion, Sebastien; Togbe, Casimir; Dayma, Guillaume; Mehl, Marco; Dagaut, Philippe; Pitz, William J.; Zador, Judit; Sarathy, S. Mani

    2016-06-22

    This study is concerned with the identification and quantification of species generated during the combustion of cyclopentane in a jet stirred reactor (JSR). Experiments were carried out for temperatures between 740 and 1250 K, equivalence ratios from 0.5 to 3.0, and at an operating pressure of 10 atm. The fuel concentration was kept at 0.1% and the residence time of the fuel/O2/N2 mixture was maintained at 0.7 s. The reactant, product, and intermediate species concentration profiles were measured using gas chromatography and Fourier transform infrared spectroscopy. The concentration profiles of cyclopentane indicate inhibition of reactivity between 850-1000 K for φ=2.0more » and φ=3.0. This behavior is interesting, as it has not been observed previously for other fuel molecules, cyclic or non-cyclic. A kinetic model including both low- and high-temperature reaction pathways was developed and used to simulate the JSR experiments. The pressure-dependent rate coefficients of all relevant reactions lying on the PES of cyclopentyl + O2, as well as the C-C and C-H scission reactions of the cyclopentyl radical were calculated at the UCCSD(T)-F12b/cc-pVTZ-F12//M06-2X/6-311++G(d,p) level of theory. The simulations reproduced the unique reactivity trend of cyclopentane and the measured concentration profiles of intermediate and product species. Furthermore, sensitivity and reaction path analyses indicate that this reactivity trend may be attributed to differences in the reactivity of allyl radical at different conditions, and it is highly sensitive to the C-C/C-H scission branching ratio of the cyclopentyl radical decomposition.« less

  6. Chemical Kinetic Modeling of Advanced Transportation Fuels

    SciTech Connect (OSTI)

    PItz, W J; Westbrook, C K; Herbinet, O

    2009-01-20

    Development of detailed chemical kinetic models for advanced petroleum-based and nonpetroleum based fuels is a difficult challenge because of the hundreds to thousands of different components in these fuels and because some of these fuels contain components that have not been considered in the past. It is important to develop detailed chemical kinetic models for these fuels since the models can be put into engine simulation codes used for optimizing engine design for maximum efficiency and minimal pollutant emissions. For example, these chemistry-enabled engine codes can be used to optimize combustion chamber shape and fuel injection timing. They also allow insight into how the composition of advanced petroleum-based and non-petroleum based fuels affect engine performance characteristics. Additionally, chemical kinetic models can be used separately to interpret important in-cylinder experimental data and gain insight into advanced engine combustion processes such as HCCI and lean burn engines. The objectives are: (1) Develop detailed chemical kinetic reaction models for components of advanced petroleum-based and non-petroleum based fuels. These fuels models include components from vegetable-oil-derived biodiesel, oil-sand derived fuel, alcohol fuels and other advanced bio-based and alternative fuels. (2) Develop detailed chemical kinetic reaction models for mixtures of non-petroleum and petroleum-based components to represent real fuels and lead to efficient reduced combustion models needed for engine modeling codes. (3) Characterize the role of fuel composition on efficiency and pollutant emissions from practical automotive engines.

  7. Transportation Sector Model of the National Energy Modeling System. Volume 2 -- Appendices: Part 1

    SciTech Connect (OSTI)

    1998-01-01

    This volume contains input data and parameters used in the model of the transportation sector of the National Energy Modeling System. The list of Transportation Sector Model variables includes parameters for the following: Light duty vehicle modules (fuel economy, regional sales, alternative fuel vehicles); Light duty vehicle stock modules; Light duty vehicle fleet module; Air travel module (demand model and fleet efficiency model); Freight transport module; Miscellaneous energy demand module; and Transportation emissions module. Also included in these appendices are: Light duty vehicle market classes; Maximum light duty vehicle market penetration parameters; Aircraft fleet efficiency model adjustment factors; and List of expected aircraft technology improvements.

  8. A model for radionuclide transport in the Cooling Water System

    SciTech Connect (OSTI)

    Kahook, S.D.

    1992-08-01

    A radionuclide transport model developed to assess radiological levels in the K-reactor Cooling Water System (CWS) in the event of an inadvertent process water (PW) leakage to the cooling water (CW) in the heat exchangers (HX) is described. During and following a process water leak, the radionuclide transport model determines the time-dependent release rates of radionuclide from the cooling water system to the environment via evaporation to the atmosphere and blow-down to the Savannah River. The developed model allows for delay times associated with the transport of the cooling water radioactivity through cooling water system components. Additionally, this model simulates the time-dependent behavior of radionuclides levels in various CWS components. The developed model is incorporated into the K-reactor Cooling Tower Activity (KCTA) code. KCTA allows the accident (heat exchanger leak rate) and the cooling tower blow-down and evaporation rates to be described as time-dependent functions. Thus, the postulated leak and the consequence of the assumed leak can be modelled realistically. This model is the first of three models to be ultimately assembled to form a comprehensive Liquid Pathway Activity System (LPAS). LPAS will offer integrated formation, transport, deposition, and release estimates for radionuclides formed in a SRS facility. Process water and river water modules are forthcoming as input and downstream components, respectively, for KCTA.

  9. Modeling carbon nanotube growth on the catalyst-substrate surface subjected to reactive plasma [

    SciTech Connect (OSTI)

    Tewari, Aarti; Sharma, Suresh C.

    2014-06-15

    The paper presents a theoretical model to study the growth of the carbon nanotube (CNT) on the catalyst substrate surface subjected to reactive plasma. The charging rate of the CNT, kinetics of electron, ions and neutral atoms, the growth rate of the CNT because of diffusion and accretion of ions on the catalyst nanoparticle inclusion of the issue of the plasma sheath is undertaken in the present model. Numerical calculations on the effect of ion density and temperature and the substrate bias on the growth of the CNT have been carried out for typical glow discharge plasma parameters. It is found that the height of CNT increases with the ion density of carbon ions and radius of CNT decreases with hydrogen ion density. The substrate bias also affects the growth rate of the CNT. The field emission characteristics from the CNTs can be analyzed from the results obtained.

  10. Multiscale modeling for fluid transport in nanosystems.

    SciTech Connect (OSTI)

    Lee, Jonathan W.; Jones, Reese E.; Mandadapu, Kranthi Kiran; Templeton, Jeremy Alan; Zimmerman, Jonathan A.

    2013-09-01

    Atomistic-scale behavior drives performance in many micro- and nano-fluidic systems, such as mircrofludic mixers and electrical energy storage devices. Bringing this information into the traditionally continuum models used for engineering analysis has proved challenging. This work describes one such approach to address this issue by developing atomistic-to-continuum multi scale and multi physics methods to enable molecular dynamics (MD) representations of atoms to incorporated into continuum simulations. Coupling is achieved by imposing constraints based on fluxes of conserved quantities between the two regions described by one of these models. The impact of electric fields and surface charges are also critical, hence, methodologies to extend finite-element (FE) MD electric field solvers have been derived to account for these effects. Finally, the continuum description can have inconsistencies with the coarse-grained MD dynamics, so FE equations based on MD statistics were derived to facilitate the multi scale coupling. Examples are shown relevant to nanofluidic systems, such as pore flow, Couette flow, and electric double layer.

  11. Lessons Learned from Alternative Transportation Fuels: Modeling Transition Dynamics

    Alternative Fuels and Advanced Vehicles Data Center [Office of Energy Efficiency and Renewable Energy (EERE)]

    Lessons Learned from Alternative Transportation Fuels: Modeling Transition Dynamics C. Welch Technical Report NREL/TP-540-39446 February 2006 Lessons Learned from Alternative Transportation Fuels: Modeling Transition Dynamics C. Welch Prepared under Task Nos. HS04.2000 and HS06.1002 Technical Report NREL/TP-540-39446 February 2006 National Renewable Energy Laboratory 1617 Cole Boulevard, Golden, Colorado 80401-3393 303-275-3000 * www.nrel.gov Operated for the U.S. Department of Energy Office of

  12. Upscaling of U(VI) Desorption and Transport from Decimeter-Scale Heterogeneity to Plume-Scale Modeling

    SciTech Connect (OSTI)

    Curtis, Gary P; Kohler, Matthias; Kannappan, Ramakrishnan; Briggs, Martin; Day-Lewis, Fred

    2015-02-24

    Reactive solute transport in aquifers is commonly affected by rate limited mass transfer. This slow mass transfer can exhibit significant control on the times required to restore contaminated aquifers to near-pristine conditions under both ambient and forced-gradient flow systems and is therefore important to understand. Both nonreactive and reactive tracer experiments provide valuable insight into the exchange of solute between mobile and immobile porosity. At the grain scale and column scale, mass transfer limitations were manifested as a concentration rebound when contaminated sediments were contacted with pristine groundwater. This behavior was successfully modeled using the multirate mass transfer model. Mass transfer observed in a 2 m long intermediate laboratory scale experiment showed significant concentration rebound in the first half meter along a flowpath through the tank and negligible rebound near the exit of the tank. Experimental observations and model simulations show that although concentration rebound was small at the end of the tank, the overall elution of uranium from of the tank was still controlled by mass transfer which was manifested by a long tail. At the field scale, mass transfer parameters inferred from geo-electrical measurements of bulk conductivity and traditional conductivity measurements of fluid samples showed significant spatial variability. Overall the improved understanding of mass transfer across multiple scales should lead to more robust reactive transport simulations and site management.

  13. Unsaturated Groundwater and Heat Transport Model

    Energy Science and Technology Software Center (OSTI)

    2008-05-15

    TOUGH2-MP is a massive parallel version of the TOUGH2 Code, designed for computationally efficient parallel simulation of isothermal and nonisothermal flows of multicomponent, multiphase fluids in one, two, and three-dimensional porous and fractured media. The code runs on computers with parallel architecture or clusters and can be used for applications, such as radioactive waste disposal, CO2 geological sequestration, environmental assessment and remediation, reservoir engineering, and groundwater hydrology. The parallel simulator has achieved orders-of-the-magnitude improvement inmore » computational time and/or modeling problem size. The parallel simulator uses fully implicit time differencing and solves large, sparse linear systems arising from discretization of the partial differential equations for mass and energy balance. A domain decomposition approach is adopted for multiphase flow simulations with coarse- granularity parallel computation. Current version of TOUGH2-MP includes following modules: EOS1, EOS2, EOS3, EOS4, EOS5, EOS7, EOS7R, E0S8, EOS9, ECO2N, EWASG, and T2R3D.« less

  14. Methods for modeling impact-induced reactivity changes in small reactors.

    SciTech Connect (OSTI)

    Tallman, Tyler N.; Radel, Tracy E.; Smith, Jeffrey A.; Villa, Daniel L.; Smith, Brandon M.; Radel, Ross F.; Lipinski, Ronald J.; Wilson, Paul Philip Hood

    2010-10-01

    This paper describes techniques for determining impact deformation and the subsequent reactivity change for a space reactor impacting the ground following a potential launch accident or for large fuel bundles in a shipping container following an accident. This technique could be used to determine the margin of subcriticality for such potential accidents. Specifically, the approach couples a finite element continuum mechanics model (Pronto3D or Presto) with a neutronics code (MCNP). DAGMC, developed at the University of Wisconsin-Madison, is used to enable MCNP geometric queries to be performed using Pronto3D output. This paper summarizes what has been done historically for reactor launch analysis, describes the impact criticality analysis methodology, and presents preliminary results using representative reactor designs.

  15. Model documentation report: Transportation sector model of the National Energy Modeling System

    SciTech Connect (OSTI)

    Not Available

    1994-03-01

    This report documents the objectives, analytical approach and development of the National Energy Modeling System (NEMS) Transportation Model (TRAN). The report catalogues and describes the model assumptions, computational methodology, parameter estimation techniques, model source code, and forecast results generated by the model. This document serves three purposes. First, it is a reference document providing a detailed description of TRAN for model analysts, users, and the public. Second, this report meets the legal requirements of the Energy Information Administration (EIA) to provide adequate documentation in support of its statistical and forecast reports (Public Law 93-275, 57(b)(1)). Third, it permits continuity in model development by providing documentation from which energy analysts can undertake model enhancements, data updates, and parameter refinements.

  16. Test of models for electron transport in laser produced plasmas

    SciTech Connect (OSTI)

    Colombant, D.G.; Manheimer, W.M.; Busquet, M.

    2005-07-15

    This paper examines five different models of electron thermal transport in laser produced spherical implosions. These are classical, classical with a flux limit f, delocalization, beam deposition model, and Fokker-Planck solutions. In small targets, the results are strongly dependent on f for flux limit models, with small f's generating very steep temperature gradients. Delocalization models are characterized by large preheat in the center of the target. The beam deposition model agrees reasonably well with the Fokker-Planck simulation results. For large, high gain fusion targets, the delocalization model shows the gain substantially reduced by the preheat. However, flux limitation models show gain largely independent of f, with the beam deposition model also showing the same high gain.

  17. Edge turbulence and transport: Text and ATF modeling

    SciTech Connect (OSTI)

    Ritz, C.P.; Rhodes, T.L.; Lin, H.; Rowan, W.L.; Bengtson, R.; Wootton, A.J. . Fusion Research Center); Carreras, B.A.; Leboeuf, J.N.; Lee, D.K.; Harris, J.; Hidalgo, C.; Bell, J.D.; Holmes, J.A.; Isler, R.; Lynch, V.E.; Uckan, T. ); Diamond, P.H.; Ware, A.S. ); Thayer, D.R. (Science Applications Inter

    1990-01-01

    We present experimental results on edge turbulence and transport from the tokamak TEXT and the torsatron ATF. The measured electrostatic fluctuations can explain the edge transport of particles and energy. Certain drive (radiation) and stabilizing (velocity shear) terms are suggested by the results. The experimental fluctuation levels and spectral widths can be reproduced by considering the nonlinear evolution of the reduced MHD equations, incorporating a thermal drive from line radiation. In the tokamak limit (with toroidal electric field) the model corresponds to the resistivity gradient mode, while in the currentless torsatron or stellarator limit it corresponds to a thermally driven drift wave.

  18. Modelling On Photogeneration Of Hydroxyl Radical In Surface Waters And Its Reactivity Towards Pharmaceutical Wastes

    SciTech Connect (OSTI)

    Das, Radha; Vione, Davide; Rubertelli, Francesca; Maurino, Valter; Minero, Claudio; Barbati, Stephane; Chiron, Serge

    2010-10-26

    This paper reports a simple model to describe the formation and reactivity of hydroxyl radicals in the whole column of freshwater lakes. It is based on empirical irradiation data and is a function of the water chemical composition (the photochemically significant parameters NPOC, nitrate, nitrite, carbonate and bicarbonate), the lake conformation best expressed as the average depth, and the water absorption spectrum in a simplified Lambert-Beer approach. The purpose is to derive the lifetime of dissolved molecules, due to reaction with OH, on the basis of their second-order rate constants with the hydroxyl radical. The model was applied to two compounds of pharmaceutical wastes ibuprofen and carbamazepine, for which the second-order rate constants for reaction with the hydroxyl radical were measured by means of the competition kinetics with 2-propanol. The measured values of the rate constants are 1.0x10{sup 10} and 1.6x10{sup 10} M{sup -1} s{sup -1} for ibuprofen and carbamazepine, respectively. The model suggests that the lifetime of a given compound can be very variable in different lakes, even more than the lifetime of different compounds in the same lake. It can be concluded that as far as the reaction with OH, is concerned the concepts of photolability and photostability, traditionally attached to definite compounds, are ecosystem-dependent at least as much as they depend on the molecule under consideration.

  19. Sediment and radionuclide transport in rivers: radionuclide transport modeling for Cattaraugus and Buttermilk Creeks, New York

    SciTech Connect (OSTI)

    Onishi, Y.; Yabusaki, S.B.; Kincaid, C.T.; Skaggs, R.L.; Walters, W.H.

    1982-12-01

    SERATRA, a transient, two-dimensional (laterally-averaged) computer model of sediment-contaminant transport in rivers, satisfactorily resolved the distribution of sediment and radionuclide concentrations in the Cattaraugus Creek stream system in New York. By modeling the physical processes of advection, diffusion, erosion, deposition, and bed armoring, SERATRA routed three sediment size fractions, including cohesive soils, to simulate three dynamic flow events. In conjunction with the sediment transport, SERATRA computed radionuclide levels in dissolved, suspended sediment, and bed sediment forms for four radionuclides (/sup 137/Cs, /sup 90/Sr, /sup 239/ /sup 240/Pu, and /sup 3/H). By accounting for time-dependent sediment-radionuclide interaction in the water column and bed, SERATA is a physically explicit model of radionuclide fate and migration. Sediment and radionuclide concentrations calculated by SERATA in the Cattaraugus Creek stream system are in reasonable agreement with measured values. SERATRA is in the field performance phase of an extensive testing program designed to establish the utility of the model as a site assessment tool. The model handles not only radionuclides but other contaminants such as pesticides, heavy metals and other toxic chemicals. Now that the model has been applied to four field sites, including the latest study of the Cattaraugus Creek stream system, it is recommended that a final model be validated through comparison of predicted results with field data from a carefully controlled tracer test at a field site. It is also recommended that a detailed laboratory flume be tested to study cohesive sediment transport, deposition, and erosion characteristics. The lack of current understanding of these characteristics is one of the weakest areas hindering the accurate assessment of the migration of radionuclides sorbed by fine sediments of silt and clay.

  20. Comparison of a radial fractional transport model with tokamak experiments

    SciTech Connect (OSTI)

    Kullberg, A. Morales, G. J.; Maggs, J. E.

    2014-03-15

    A radial fractional transport model [Kullberg et al., Phys. Rev. E 87, 052115 (2013)], that correctly incorporates the geometric effects of the domain near the origin and removes the singular behavior at the outer boundary, is compared to results of off-axis heating experiments performed in the Rijnhuizen Tokamak Project (RTP), ASDEX Upgrade, JET, and DIII-D tokamak devices. This comparative study provides an initial assessment of the presence of fractional transport phenomena in magnetic confinement experiments. It is found that the nonlocal radial model is robust in describing the steady-state temperature profiles from RTP, but for the propagation of heat waves in ASDEX Upgrade, JET, and DIII-D the model is not clearly superior to predictions based on Fick's law. However, this comparative study does indicate that the order of the fractional derivative, ?, is likely a function of radial position in the devices surveyed.

  1. Anisotropic Elastic Resonance Scattering model for the Neutron Transport equation

    SciTech Connect (OSTI)

    Mohamed Ouisloumen; Abderrafi M. Ougouag; Shadi Z. Ghrayeb

    2014-11-24

    The resonance scattering transfer cross-section has been reformulated to account for anisotropic scattering in the center-of-mass of the neutron-nucleus system. The main innovation over previous implementations is the relaxation of the ubiquitous assumption of isotropic scattering in the center-of-mass and the actual effective use of scattering angle distributions from evaluated nuclear data files in the computation of the angular moments of the resonant scattering kernels. The formulas for the high order anisotropic moments in the laboratory system are also derived. A multi-group numerical formulation is derived and implemented into a module incorporated within the NJOY nuclear data processing code. An ultra-fine energy mesh cross section library was generated using these new theoretical models and then was used for fuel assembly calculations with the PARAGON lattice physics code. The results obtained indicate a strong effect of this new model on reactivity, multi-group fluxes and isotopic inventory during depletion.

  2. Understanding composite explosive energetics: 3, Reactive flow modeling of aluminum reaction kinetics in PETN and TNT

    SciTech Connect (OSTI)

    Tao, W.C.; Tarver, C.M.; Ornellas, D.L.

    1991-12-06

    Using Fabry-Perot interferometry techniques, we have determined that early time rate of energy release from detonating PETN and TNT explosives filled with 5 and 10 wt % of either 5 {mu}m of 18 {mu}m spherical aluminum (Al) particles. From the measured particle velocity data, we are able to infer the reaction rate of aluminum with the detonation products, and calculate the extent of reaction 1--3 {mu}s after the detonation. We observed that a substantional portion of the aluminum metal in all of the PETN and TNE formulations reacted within the timeframe of the one-dimensional experiment. In the PETN formulation filed with 5 wt % of 5 {mu}m aluminum, all of the metal reacted within 1.5 {mu}s, resulting in an increase of 22% in energy compared to pure PETN. A reactive-flow hydrodynamic model based on the Zeldovich-von Neumann-Doring (ZND) description of the reaction zone and subsequent reaction produce expansion (Taylor wave) is used to interpret the reaction rate of the aluminum particles with detonation product gases. The diffusion-controlled reaction mechanism for aluminum and the global kinetic parameters used in the model have been found to be consistent for all the PETN and TNT formulations.

  3. Data Assimilation in the ADAPT Photospheric Flux Transport Model

    SciTech Connect (OSTI)

    Hickmann, Kyle S.; Godinez, Humberto C.; Henney, Carl J.; Arge, C. Nick

    2015-03-17

    Global maps of the solar photospheric magnetic flux are fundamental drivers for simulations of the corona and solar wind and therefore are important predictors of geoeffective events. However, observations of the solar photosphere are only made intermittently over approximately half of the solar surface. The Air Force Data Assimilative Photospheric Flux Transport (ADAPT) model uses localized ensemble Kalman filtering techniques to adjust a set of photospheric simulations to agree with the available observations. At the same time, this information is propagated to areas of the simulation that have not been observed. ADAPT implements a local ensemble transform Kalman filter (LETKF) to accomplish data assimilation, allowing the covariance structure of the flux-transport model to influence assimilation of photosphere observations while eliminating spurious correlations between ensemble members arising from a limited ensemble size. We give a detailed account of the implementation of the LETKF into ADAPT. Advantages of the LETKF scheme over previously implemented assimilation methods are highlighted.

  4. Modeling of transport phenomena in tokamak plasmas with neural networks

    SciTech Connect (OSTI)

    Meneghini, O.; Luna, C. J.; Smith, S. P.; Lao, L. L.

    2014-06-15

    A new transport model that uses neural networks (NNs) to yield electron and ion heat flux profiles has been developed. Given a set of local dimensionless plasma parameters similar to the ones that the highest fidelity models use, the NN model is able to efficiently and accurately predict the ion and electron heat transport profiles. As a benchmark, a NN was built, trained, and tested on data from the 2012 and 2013 DIII-D experimental campaigns. It is found that NN can capture the experimental behavior over the majority of the plasma radius and across a broad range of plasma regimes. Although each radial location is calculated independently from the others, the heat flux profiles are smooth, suggesting that the solution found by the NN is a smooth function of the local input parameters. This result supports the evidence of a well-defined, non-stochastic relationship between the input parameters and the experimentally measured transport fluxes. The numerical efficiency of this method, requiring only a few CPU-?s per data point, makes it ideal for scenario development simulations and real-time plasma control.

  5. Engineered Barrier System Degradation, Flow, and Transport Process Model Report

    SciTech Connect (OSTI)

    E.L. Hardin

    2000-07-17

    The Engineered Barrier System Degradation, Flow, and Transport Process Model Report (EBS PMR) is one of nine PMRs supporting the Total System Performance Assessment (TSPA) being developed by the Yucca Mountain Project for the Site Recommendation Report (SRR). The EBS PMR summarizes the development and abstraction of models for processes that govern the evolution of conditions within the emplacement drifts of a potential high-level nuclear waste repository at Yucca Mountain, Nye County, Nevada. Details of these individual models are documented in 23 supporting Analysis/Model Reports (AMRs). Nineteen of these AMRs are for process models, and the remaining 4 describe the abstraction of results for application in TSPA. The process models themselves cluster around four major topics: ''Water Distribution and Removal Model, Physical and Chemical Environment Model, Radionuclide Transport Model, and Multiscale Thermohydrologic Model''. One AMR (Engineered Barrier System-Features, Events, and Processes/Degradation Modes Analysis) summarizes the formal screening analysis used to select the Features, Events, and Processes (FEPs) included in TSPA and those excluded from further consideration. Performance of a potential Yucca Mountain high-level radioactive waste repository depends on both the natural barrier system (NBS) and the engineered barrier system (EBS) and on their interactions. Although the waste packages are generally considered as components of the EBS, the EBS as defined in the EBS PMR includes all engineered components outside the waste packages. The principal function of the EBS is to complement the geologic system in limiting the amount of water contacting nuclear waste. A number of alternatives were considered by the Project for different EBS designs that could provide better performance than the design analyzed for the Viability Assessment. The design concept selected was Enhanced Design Alternative II (EDA II).

  6. Transport Studies and Modeling in PEM Fuel Cells

    SciTech Connect (OSTI)

    Mittelsteadt, Cortney K.; Xu, Hui; Brawn, Shelly

    2014-07-30

    This project’s aim was to develop fuel cell components (i.e. membranes, gas-diffusion media (GDM), bipolar plates and flow fields) that possess specific properties (i.e. water transport and conductivity). A computational fluid dynamics model was developed to elucidate the effect of certain parameters on these specific properties. Ultimately, the model will be used to determine sensitivity of fuel cell performance to component properties to determine limiting components and to guide research. We have successfully reached our objectives and achieved most of the milestones of this project. We have designed and synthesized a variety of hydrocarbon block polymer membranes with lower equivalent weight, structure, chemistry, phase separation and process conditions. These membranes provide a broad selection with optimized water transport properties. We have also designed and constructed a variety of devices that are capable of accurately measuring the water transport properties (water uptake, water diffusivity and electro-osmatic drag) of these membranes. These transport properties are correlated to the membranes’ structures derived from X-ray and microscopy techniques to determine the structure-property relationship. We successfully integrated hydrocarbon membrane MEAs with a current distribution board (CBD) to study the impact of hydrocarbon membrane on water transport in fuel cells. We have designed and fabricated various GDM with varying substrate, diffusivity and micro-porous layers (MPL) and characterized their pore structure, tortuosity and hydrophobicity. We have derived a universal chart (MacMullin number as function of wet proofing and porosity) that can be used to characterize various GDM. The abovementioned GDMs have been evaluated in operating fuel cells; their performance is correlated to various pore structure, tortuosity and hydrophobicity of the GDM. Unfortunately, determining a universal relationship between the MacMullin number and these properties

  7. Building 235-F Goldsim Fate And Transport Model

    SciTech Connect (OSTI)

    Taylor, G. A.; Phifer, M. A.

    2012-09-14

    Savannah River National Laboratory (SRNL) personnel, at the request of Area Completion Projects (ACP), evaluated In-Situ Disposal (ISD) alternatives that are under consideration for deactivation and decommissioning (D&D) of Building 235-F and the Building 294-2F Sand Filter. SRNL personnel developed and used a GoldSim fate and transport model, which is consistent with Musall 2012, to evaluate relative to groundwater protection, ISD alternatives that involve either source removal and/or the grouting of portions or all of 235-F. This evaluation was conducted through the development and use of a Building 235-F GoldSim fate and transport model. The model simulates contaminant release from four 235-F process areas and the 294-2F Sand Filter. In addition, it simulates the fate and transport through the vadose zone, the Upper Three Runs (UTR) aquifer, and the Upper Three Runs (UTR) creek. The model is designed as a stochastic model, and as such it can provide both deterministic and stochastic (probabilistic) results. The results show that the median radium activity concentrations exceed the 5 ?Ci/L radium MCL at the edge of the building for all ISD alternatives after 10,000 years, except those with a sufficient amount of inventory removed. A very interesting result was that grouting was shown to basically have minimal effect on the radium activity concentration. During the first 1,000 years grouting may have some small positive benefit relative to radium, however after that it may have a slightly deleterious effect. The Pb-210 results, relative to its 0.06 ?Ci/L PRG, are essentially identical to the radium results, but the Pb-210 results exhibit a lesser degree of exceedance. In summary, some level of inventory removal will be required to ensure that groundwater standards are met.

  8. Compendium of Material Composition Data for Radiation Transport Modeling

    SciTech Connect (OSTI)

    Williams, Ralph G.; Gesh, Christopher J.; Pagh, Richard T.

    2006-10-31

    Computational modeling of radiation transport problems including homeland security, radiation shielding and protection, and criticality safety all depend upon material definitions. This document has been created to serve two purposes: 1) to provide a quick reference of material compositions for analysts and 2) a standardized reference to reduce the differences between results from two independent analysts. Analysts are always encountering a variety of materials for which elemental definitions are not readily available or densities are not defined. This document provides a location where unique or hard to define materials will be located to reduce duplication in research for modeling purposes. Additionally, having a common set of material definitions helps to standardize modeling across PNNL and provide two separate researchers the ability to compare different modeling results from a common materials basis.

  9. Model documentation report: Transportation sector model of the National Energy Modeling System

    SciTech Connect (OSTI)

    1997-02-01

    Over the past year, several modifications have been made to the NEMS Transportation Model, incorporating greater levels of detail and analysis in modules previously represented in the aggregate or under a profusion of simplifying assumptions. This document is intended to amend those sections of the Model Documentation Report (MDR) which describe these superseded modules. Significant changes have been implemented in the LDV Fuel Economy Model, the Alternative Fuel Vehicle Model, the LDV Fleet Module, and the Highway Freight Model. The relevant sections of the MDR have been extracted from the original document, amended, and are presented in the following pages. A brief summary of the modifications follows: In the Fuel Economy Model, modifications have been made which permit the user to employ more optimistic assumptions about the commercial viability and impact of selected technological improvements. This model also explicitly calculates the fuel economy of an array of alternative fuel vehicles (AFV`s) which are subsequently used in the estimation of vehicle sales. In the Alternative Fuel Vehicle Model, the results of the Fuel Economy Model have been incorporated, and the program flows have been modified to reflect that fact. In the Light Duty Vehicle Fleet Module, the sales of vehicles to fleets of various size are endogenously calculated in order to provide a more detailed estimate of the impacts of EPACT legislation on the sales of AFV`s to fleets. In the Highway Freight Model, the previous aggregate estimation has been replaced by a detailed Freight Truck Stock Model, where travel patterns, efficiencies, and energy intensities are estimated by industrial grouping. Several appendices are provided at the end of this document, containing data tables and supplementary descriptions of the model development process which are not integral to an understanding of the overall model structure.

  10. Turbulence transport modeling of the temporal outer heliosphere

    SciTech Connect (OSTI)

    Adhikari, L.; Zank, G. P.; Hu, Q.; Dosch, A.

    2014-09-20

    The solar wind can be regarded as a turbulent magnetofluid, evolving in an expanding solar wind and subject to turbulent driving by a variety of in situ sources. Furthermore, the solar wind and the drivers of turbulence are highly time-dependent and change with solar cycle. Turbulence transport models describing low-frequency magnetic and velocity fluctuations in the solar wind have so far neglected solar cycle effects. Here we consider the effects of solar cycle variability on a turbulence transport model developed by Zank et al. This model is appropriate for the solar wind beyond about 1 AU, and extensions have described the steady-state dependence of the magnetic energy density fluctuations, correlation length, and solar wind temperature throughout the outer heliosphere. We find that the temporal solar wind introduces a periodic variability, particularly beyond ?10 AU, in the magnetic energy density fluctuations, correlation length, and solar wind temperature. The variability is insufficient to account for the full observed variability in these quantities, but we find that the time-dependent solutions trace the steady-state solutions quite well, suggesting that the steady-state models are reasonable first approximations.

  11. Transportation Sector Model of the National Energy Modeling System. Volume 1

    SciTech Connect (OSTI)

    1998-01-01

    This report documents the objectives, analytical approach and development of the National Energy Modeling System (NEMS) Transportation Model (TRAN). The report catalogues and describes the model assumptions, computational methodology, parameter estimation techniques, model source code, and forecast results generated by the model. The NEMS Transportation Model comprises a series of semi-independent models which address different aspects of the transportation sector. The primary purpose of this model is to provide mid-term forecasts of transportation energy demand by fuel type including, but not limited to, motor gasoline, distillate, jet fuel, and alternative fuels (such as CNG) not commonly associated with transportation. The current NEMS forecast horizon extends to the year 2010 and uses 1990 as the base year. Forecasts are generated through the separate consideration of energy consumption within the various modes of transport, including: private and fleet light-duty vehicles; aircraft; marine, rail, and truck freight; and various modes with minor overall impacts, such as mass transit and recreational boating. This approach is useful in assessing the impacts of policy initiatives, legislative mandates which affect individual modes of travel, and technological developments. The model also provides forecasts of selected intermediate values which are generated in order to determine energy consumption. These elements include estimates of passenger travel demand by automobile, air, or mass transit; estimates of the efficiency with which that demand is met; projections of vehicle stocks and the penetration of new technologies; and estimates of the demand for freight transport which are linked to forecasts of industrial output. Following the estimation of energy demand, TRAN produces forecasts of vehicular emissions of the following pollutants by source: oxides of sulfur, oxides of nitrogen, total carbon, carbon dioxide, carbon monoxide, and volatile organic compounds.

  12. THE LOS ALAMOS NATIONAL LABORATORY ATMOSPHERIC TRANSPORT AND DIFFUSION MODELS

    SciTech Connect (OSTI)

    M. WILLIAMS

    1999-08-01

    The LANL atmospheric transport and diffusion models are composed of two state-of-the-art computer codes. The first is an atmospheric wind model called HOThlAC, Higher Order Turbulence Model for Atmospheric circulations. HOTMAC generates wind and turbulence fields by solving a set of atmospheric dynamic equations. The second is an atmospheric diffusion model called RAPTAD, Random Particle Transport And Diffusion. RAPTAD uses the wind and turbulence output from HOTMAC to compute particle trajectories and concentration at any location downwind from a source. Both of these models, originally developed as research codes on supercomputers, have been modified to run on microcomputers. Because the capability of microcomputers is advancing so rapidly, the expectation is that they will eventually become as good as today's supercomputers. Now both models are run on desktop or deskside computers, such as an IBM PC/AT with an Opus Pm 350-32 bit coprocessor board and a SUN workstation. Codes have also been modified so that high level graphics, NCAR Graphics, of the output from both models are displayed on the desktop computer monitors and plotted on a laser printer. Two programs, HOTPLT and RAPLOT, produce wind vector plots of the output from HOTMAC and particle trajectory plots of the output from RAPTAD, respectively. A third CONPLT provides concentration contour plots. Section II describes step-by-step operational procedures, specifically for a SUN-4 desk side computer, on how to run main programs HOTMAC and RAPTAD, and graphics programs to display the results. Governing equations, boundary conditions and initial values of HOTMAC and RAPTAD are discussed in Section III. Finite-difference representations of the governing equations, numerical solution procedures, and a grid system are given in Section IV.

  13. MODELING TRITIUM TRANSPORT, DEPOSITION AND RE-EMISSION

    SciTech Connect (OSTI)

    Viner, B.

    2012-04-03

    The atmospheric release of tritium oxide (HTO) potentially impacts human health, typically through inhalation or absorption. Due to HTO's similarity to water, vegetation will absorb HTO by solution in the leaf water and then re-emit it, creating a number of secondary sources of HTO. Currently, models used for emergency response at Savannah River Site incorporate the transport and deposition of HTO but do not provide estimates for its potential re-emission from vegetation or soil surface though re-emission could result in prolonged exposure and greater than predicted dose for an individual downwind. A simple model of HTO transport, deposition and re-emission has been developed to examine the potential increase in exposure and dose. The model simulates an initial release of HTO that moves with a mean wind and expands through diffusion as a Gaussian puff. Deposition is modeled using previous estimates of deposition velocity for HTO and re-emission is modeled using a time constant that describes how quickly HTO is transferred between the surface and atmosphere. Additional puffs are created to simulate re-emission of HTO as well as horizontal diffusion across model grid cells. An evaluation of field data indicates that the use of a re-emission module tends to improve model predictions through improved prediction of peak concentration magnitude and location. When considering dose, nearly all of the released material is included in the dose calculation when re-emission is included. Although exposure to HTO through re-emission occurs over a few hours, the incremental increase in dose is relatively small because the atmospheric concentration of re-emitted HTO is much lower than the initial release.

  14. Effect of a detailed radial core expansion reactivity feedback model on ATWS calculations using SASSYS/SAS4A

    SciTech Connect (OSTI)

    Wigeland, R.A.

    1986-01-01

    The present emphasis on inherent safety and inherently safe designs for liquid-metal reactors has resulted in a need to represent the various reactivity feedback mechanisms as accurately as possible. In particular, the reactivity feedback from radial core expansion has been found to provide the dominant negative feedback contribution in postulated anticipated transient without scram (ATWS) events. Review of the existing modeling in the SASSYS/SAS4A computer code system revealed that while the modeling may be adequate for the early phases of various unprotected transients, the accuracy would be less than desirable for the extended transients which typically occur for inherently safe designs. The existing model for calculating the reactivity feedback from radial core expansion uses a feedback from radial core expansion uses a feedback coefficient in conjunction with changes in the temperatures of the grid support plate and the above-core load pad. The accuracy of this approach is determined partly by the conditions used in deriving the feedback coefficient, and their relevance to the transient being investigated. Accuracy is also affected by the need to include effects other than those that could be directly related to changes in the grid plate and above-core load pad temperatures, such as subassembly bowing and the potential for clearances to occur between subassemblies in the above-core load pad region. As a result, a detailed model was developed in an attempt to account for these and other effects in a more mechanistic form.

  15. Modeling interfacial area transport in multi-fluid systems

    SciTech Connect (OSTI)

    Yarbro, S.L.

    1996-11-01

    Many typical chemical engineering operations are multi-fluid systems. They are carried out in distillation columns (vapor/liquid), liquid-liquid contactors (liquid/liquid) and other similar devices. An important parameter is interfacial area concentration, which determines the rate of interfluid heat, mass and momentum transfer and ultimately, the overall performance of the equipment. In many cases, the models for determining interfacial area concentration are empirical and can only describe the cases for which there is experimental data. In an effort to understand multiphase reactors and the mixing process better, a multi-fluid model has been developed as part of a research effort to calculate interfacial area transport in several different types of in-line static mixers. For this work, the ensemble-averaged property conservation equations have been derived for each fluid and for the mixture. These equations were then combined to derive a transport equation for the interfacial area concentration. The final, one-dimensional model was compared to interfacial area concentration data from two sizes of Kenics in-line mixer, two sizes of concurrent jet and a Tee mixer. In all cases, the calculated and experimental data compared well with the highest scatter being with the Tee mixer comparison.

  16. A computerized coal-water slurry transportation model

    SciTech Connect (OSTI)

    Ljubicic, B.R.; Trostad, B.; Bukurov, Z.; Cvijanovic, P.

    1995-12-01

    Coal-water fuel (CWF) technology has been developed to the point where full-scale commercialization is just a matter of gaining sufficient market confidence in the price stability of alternate fossil fuels. In order to generalize alternative fuel cost estimates for the desired combinations of processing and/or transportation, a great deal of flexibility is required owing to the understood lack of precision in many of the newly emerging coal technologies. Previously, decisions regarding the sequential and spatial arrangement of the various process steps were made strictly on the basis of experience, simplified analysis, and intuition. Over the last decade, computer modeling has progressed from empirically based correlation to that of intricate mechanistic analysis. Nomograms, charts, tables, and many simple rules of thumb have been made obsolete by the availability of complex computer models. Given the ability to view results graphically in real or near real time, the engineer can immediately verify, from a practical standpoint, whether the initial assumptions and inputs were indeed valid. If the feasibility of a project is being determined in the context of a lack of specific data, the ability to provide a dynamic software-based solution is crucial. Furthermore, the resulting model can be used to establish preliminary operating procedures, test control logic, and train plant/process operators. Presented in this paper is a computerized model capable of estimating the delivered cost of CWF. The model uses coal-specific values, process and transport requirements, terrain factors, and input costs to determine the final operating configuration, bill of materials, and, ultimately, the capital, operating, and unit costs.

  17. Dynamic modeling of injection-induced fault reactivation and ground motion and impact on surface structures and human perception

    SciTech Connect (OSTI)

    Rutqvist, Jonny; Cappa, Frederic; Rinaldi, Antonio P.; Godano, Maxime

    2014-12-31

    We summarize recent modeling studies of injection-induced fault reactivation, seismicity, and its potential impact on surface structures and nuisance to the local human population. We used coupled multiphase fluid flow and geomechanical numerical modeling, dynamic wave propagation modeling, seismology theories, and empirical vibration criteria from mining and construction industries. We first simulated injection-induced fault reactivation, including dynamic fault slip, seismic source, wave propagation, and ground vibrations. From co-seismic average shear displacement and rupture area, we determined the moment magnitude to about Mw = 3 for an injection-induced fault reactivation at a depth of about 1000 m. We then analyzed the ground vibration results in terms of peak ground acceleration (PGA), peak ground velocity (PGV), and frequency content, with comparison to the U.S. Bureau of Mines vibration criteria for cosmetic damage to buildings, as well as human-perception vibration limits. For the considered synthetic Mw = 3 event, our analysis showed that the short duration, high frequency ground motion may not cause any significant damage to surface structures, and would not cause, in this particular case, upward CO2 leakage, but would certainly be felt by the local population.

  18. Dynamic modeling of injection-induced fault reactivation and ground motion and impact on surface structures and human perception

    DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)

    Rutqvist, Jonny; Cappa, Frederic; Rinaldi, Antonio P.; Godano, Maxime

    2014-12-31

    We summarize recent modeling studies of injection-induced fault reactivation, seismicity, and its potential impact on surface structures and nuisance to the local human population. We used coupled multiphase fluid flow and geomechanical numerical modeling, dynamic wave propagation modeling, seismology theories, and empirical vibration criteria from mining and construction industries. We first simulated injection-induced fault reactivation, including dynamic fault slip, seismic source, wave propagation, and ground vibrations. From co-seismic average shear displacement and rupture area, we determined the moment magnitude to about Mw = 3 for an injection-induced fault reactivation at a depth of about 1000 m. We then analyzed themore » ground vibration results in terms of peak ground acceleration (PGA), peak ground velocity (PGV), and frequency content, with comparison to the U.S. Bureau of Mines’ vibration criteria for cosmetic damage to buildings, as well as human-perception vibration limits. For the considered synthetic Mw = 3 event, our analysis showed that the short duration, high frequency ground motion may not cause any significant damage to surface structures, and would not cause, in this particular case, upward CO2 leakage, but would certainly be felt by the local population.« less

  19. Modeling two-fluid-phase flow and species transport in porous...

    Office of Scientific and Technical Information (OSTI)

    Publisher's Accepted Manuscript: Modeling two-fluid-phase flow and species transport in porous media This content will become publicly available on March 3, 2017 Title: Modeling ...

  20. Data Assimilation in the ADAPT Photospheric Flux Transport Model

    DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)

    Hickmann, Kyle S.; Godinez, Humberto C.; Henney, Carl J.; Arge, C. Nick

    2015-03-17

    Global maps of the solar photospheric magnetic flux are fundamental drivers for simulations of the corona and solar wind and therefore are important predictors of geoeffective events. However, observations of the solar photosphere are only made intermittently over approximately half of the solar surface. The Air Force Data Assimilative Photospheric Flux Transport (ADAPT) model uses localized ensemble Kalman filtering techniques to adjust a set of photospheric simulations to agree with the available observations. At the same time, this information is propagated to areas of the simulation that have not been observed. ADAPT implements a local ensemble transform Kalman filter (LETKF)more » to accomplish data assimilation, allowing the covariance structure of the flux-transport model to influence assimilation of photosphere observations while eliminating spurious correlations between ensemble members arising from a limited ensemble size. We give a detailed account of the implementation of the LETKF into ADAPT. Advantages of the LETKF scheme over previously implemented assimilation methods are highlighted.« less

  1. Modeling Studies on the Transport of Benzene and H2S in CO2-Water Systems

    SciTech Connect (OSTI)

    Zheng, L.; Spycher, N.; Xu, T.; Apps, J.; Kharaka, Y.; Birkholzer, J.T.

    2010-11-05

    In this study, reactive transport simulations were used to assess the mobilization and transport of organics with supercritical CO{sub 2} (SCC), and the co-injection and transport of H{sub 2}S with SCC. These processes were evaluated at conditions of typical storage reservoirs, and for cases of hypothetical leakage from a reservoir to an overlying shallower fresh water aquifer. Modeling capabilities were developed to allow the simulation of multiphase flow and transport of H{sub 2}O, CO{sub 2}, H{sub 2}S, as well as specific organic compounds (benzene), coupled with multicomponent geochemical reaction and transport. This included the development of a new simulator, TMVOC-REACT, starting from existing modules of the TOUGH2 family of codes. This work also included an extensive literature review, calculation, and testing of phase-partitioning properties for mixtures of the phases considered. The reactive transport simulations presented in this report are primarily intended to illustrate the capabilities of the new simulator. They are also intended to help evaluate and understand various processes at play, in a more qualitative than quantitative manner, and only for hypothetical scenarios. Therefore, model results are not intended as realistic assessments of groundwater quality changes for specific locations, and they certainly do not provide an exhaustive evaluation of all possible site conditions, especially given the large variability and uncertainty in hydrogeologic and geochemical parameter input into simulations. The first step in evaluating the potential mobilization and transport of organics was the identification of compounds likely to be present in deep storage formations, and likely to negatively impact freshwater aquifers if mobilized by SCC. On the basis of a literature review related to the occurrence of these organic compounds, their solubility in water and SCC, and their toxicity (as reflected by their maximum contaminant levels MCL), benzene was

  2. Modelling transport phenomena in a multi-physics context

    SciTech Connect (OSTI)

    Marra, Francesco

    2015-01-22

    Innovative heating research on cooking, pasteurization/sterilization, defrosting, thawing and drying, often focuses on areas which include the assessment of processing time, evaluation of heating uniformity, studying the impact on quality attributes of the final product as well as considering the energy efficiency of these heating processes. During the last twenty years, so-called electro-heating-processes (radio-frequency - RF, microwaves - MW and ohmic - OH) gained a wide interest in industrial food processing and many applications using the above mentioned technologies have been developed with the aim of reducing processing time, improving process efficiency and, in many cases, the heating uniformity. In the area of innovative heating, electro-heating accounts for a considerable portion of both the scientific literature and commercial applications, which can be subdivided into either direct electro-heating (as in the case of OH heating) where electrical current is applied directly to the food or indirect electro-heating (e.g. MW and RF heating) where the electrical energy is firstly converted to electromagnetic radiation which subsequently generates heat within a product. New software packages, which make easier solution of PDEs based mathematical models, and new computers, capable of larger RAM and more efficient CPU performances, allowed an increasing interest about modelling transport phenomena in systems and processes - as the ones encountered in food processing - that can be complex in terms of geometry, composition, boundary conditions but also - as in the case of electro-heating assisted applications - in terms of interaction with other physical phenomena such as displacement of electric or magnetic field. This paper deals with the description of approaches used in modelling transport phenomena in a multi-physics context such as RF, MW and OH assisted heating.

  3. TOUGH2. Unsaturated Groundwater and Heat Transport Model

    SciTech Connect (OSTI)

    Pruess, K.

    1991-05-01

    TOUGH2 is a new and improved version of TOUGH. TOUGH2 offers added capabilities and user features, including the flexibility to handle different fluid mixtures (water, water with tracer; water, CO2; water, air; water, air, with vapor pressure lowering and water, hydrogen), facilities for processing of geometric data (computational grids), and an internal version control system to ensure referenceability of code applications. TOUGH2 is a multi-dimensional numerical model for simulating the coupled transport of water, vapor, air, and heat in porous and fractured media. The program provides options for specifying injection or withdrawal of heat and fluids. Although primarily designed for studies of high-level nuclear waste isolation in partially saturated geological media, it should also be useful for a wider range of problems in heat and moisture transfer, and in the drying of porous materials. For example, geothermal reservoir simulation problems can be handled simply by setting the air mass function equal to zero on input. The TOUGH2 simulator was developed for problems involving strongly heat-driven flow. To describe these phenomena a multi-phase approach to fluid and heat flow is used, which fully accounts for the movement of gaseous and liquid phases, their transport of latent and sensible heat, and phase transitions between liquid and vapor. TOUGH2 takes account of fluid flow in both liquid and gaseous phases occurring under pressure, viscous, and gravity forces according to Darcy`s law. Interference between the phases is represented by means of relative permeability functions. The code handles binary, but not Knudsen, diffusion in the gas phase and capillary and phase adsorption effects for the liquid phase. Heat transport occurs by means of conduction with thermal conductivity dependent on water saturation, convection, and binary diffusion, which includes both sensible and latent heat.

  4. Compendium of Material Composition Data for Radiation Transport Modeling

    SciTech Connect (OSTI)

    McConn, Ronald J.; Gesh, Christopher J.; Pagh, Richard T.; Rucker, Robert A.; Williams III, Robert

    2011-03-04

    Introduction Meaningful simulations of radiation transport applications require realistic definitions of material composition and densities. When seeking that information for applications in fields such as homeland security, radiation shielding and protection, and criticality safety, researchers usually encounter a variety of materials for which elemental compositions are not readily available or densities are not defined. Publication of the Compendium of Material Composition Data for Radiation Transport Modeling, Revision 0, in 2006 was the first step toward mitigating this problem. Revision 0 of this document listed 121 materials, selected mostly from the combined personal libraries of staff at the Pacific Northwest National Laboratory (PNNL), and thus had a scope that was recognized at the time to be limited. Nevertheless, its creation did provide a well-referenced source of some unique or hard-to-define material data in a format that could be used directly in radiation transport calculations being performed at PNNL. Moreover, having a single common set of material definitions also helped to standardize at least one aspect of the various modeling efforts across the laboratory by providing separate researchers the ability to compare different model results using a common basis of materials. The authors of the 2006 compendium understood that, depending on its use and feedback, the compendium would need to be revised to correct errors or inconsistencies in the data for the original 121 materials, as well as to increase (per users suggestions) the number of materials listed. This 2010 revision of the compendium has accomplished both of those objectives. The most obvious change is the increased number of materials from 121 to 372. The not-so-obvious change is the mechanism used to produce the data listed here. The data listed in the 2006 document were compiled, evaluated, entered, and error-checked by a group of individuals essentially by hand, providing no library

  5. transportation

    National Nuclear Security Administration (NNSA)

    security missions undertaken by the U.S. government.

    Pantex Plant's Calvin Nelson honored as Analyst of the Year for Transportation Security http:nnsa.energy.gov...

  6. Active patterning and asymmetric transport in a model actomyosin network

    SciTech Connect (OSTI)

    Wang, Shenshen; Wolynes, Peter G.

    2013-12-21

    Cytoskeletal networks, which are essentially motor-filament assemblies, play a major role in many developmental processes involving structural remodeling and shape changes. These are achieved by nonequilibrium self-organization processes that generate functional patterns and drive intracellular transport. We construct a minimal physical model that incorporates the coupling between nonlinear elastic responses of individual filaments and force-dependent motor action. By performing stochastic simulations we show that the interplay of motor processes, described as driving anti-correlated motion of the network vertices, and the network connectivity, which determines the percolation character of the structure, can indeed capture the dynamical and structural cooperativity which gives rise to diverse patterns observed experimentally. The buckling instability of individual filaments is found to play a key role in localizing collapse events due to local force imbalance. Motor-driven buckling-induced node aggregation provides a dynamic mechanism that stabilizes the two-dimensional patterns below the apparent static percolation limit. Coordinated motor action is also shown to suppress random thermal noise on large time scales, the two-dimensional configuration that the system starts with thus remaining planar during the structural development. By carrying out similar simulations on a three-dimensional anchored network, we find that the myosin-driven isotropic contraction of a well-connected actin network, when combined with mechanical anchoring that confers directionality to the collective motion, may represent a novel mechanism of intracellular transport, as revealed by chromosome translocation in the starfish oocyte.

  7. A Two-Dimensional Model for the Analysis of Contaminant Transport in a Fractured Porous Medium.

    Energy Science and Technology Software Center (OSTI)

    1991-03-05

    CHAINT-MC is a two-dimensional finite element model applicable to the transport of a dissolved radionuclide in a fractured porous medium along with radioactive chain decay and subsequent transport of the dissolved daughters.

  8. A generalized Fokker-Planck model applied to beam transport

    SciTech Connect (OSTI)

    Prinja, A.K.; Pomraning, G.C.

    1998-12-31

    In this paper the authors consider a generalized Fokker-Planck (GFP) model for the transport of charged-particle beams, when scattering is not sufficiently forward peaked for the well-known Fokker-Planck (FP) approximation to hold. This approximation is strictly valid in the limit that the total scattering cross section {sigma}{sub s} {r_arrow} {infinity} and the mean cosine of scattering {bar {mu}} {r_arrow} 1. When scattering is not sufficiently forward peaked, higher order Fokker-Planck expansions, restricted to differential cross sections that fall off sufficiently rapidly, can be developed. For instance, the Henyey-Greenstein differential cross section does not possess an FP limit at all, while the screened Rutherford cross section satisfies the FP limit only marginally, and its higher order FP limits do not exist. For cross sections that fall off faster than screened Rutherford, including the exponential phase function, higher order truncations exist.

  9. Reduced Lorenz models for anomalous transport and profile resilience

    SciTech Connect (OSTI)

    Rypdal, K.; Garcia, O. E.

    2007-02-15

    The physical basis for the Lorenz equations for convective cells in stratified fluids, and for magnetized plasmas imbedded in curved magnetic fields, are reexamined with emphasis on anomalous transport. It is shown that the Galerkin truncation leading to the Lorenz equations for the closed boundary problem is incompatible with finite fluxes through the system in the limit of vanishing diffusion. An alternative formulation leading to the Lorenz equations is proposed, invoking open boundaries and the notion of convective streamers and their back-reaction on the profile gradient, giving rise to resilience of the profile. Particular emphasis is put on the diffusionless limit, where these equations reduce to a simple dynamical system depending only on one single forcing parameter. This model is studied numerically, stressing experimentally observable signatures, and some of the perils of dimension-reducing approximations are discussed.

  10. Critical review: Radionuclide transport, sediment transport, and water quality mathematical modeling; and radionuclide adsorption/desorption mechanisms

    SciTech Connect (OSTI)

    Onishi, Y.; Serne, R.J.; Arnold, E.M.; Cowan, C.E.; Thompson, F.L.

    1981-01-01

    This report describes the results of a detailed literature review of radionuclide transport models applicable to rivers, estuaries, coastal waters, the Great Lakes, and impoundments. Some representatives sediment transport and water quality models were also reviewed to evaluate if they can be readily adapted to radionuclide transport modeling. The review showed that most available transport models were developed for dissolved radionuclide in rivers. These models include the mechanisms of advection, dispersion, and radionuclide decay. Since the models do not include sediment and radionuclide interactions, they are best suited for simulating short-term radionuclide migration where: (1) radionuclides have small distribution coefficients; (2) sediment concentrations in receiving water bodies are very low. Only 5 of the reviewed models include full sediment and radionuclide interactions: CHMSED developed by Fields; FETRA SERATRA, and TODAM developed by Onishi et al, and a model developed by Shull and Gloyna. The 5 models are applicable to cases where: (1) the distribution coefficient is large; (2) sediment concentrations are high; or (3) long-term migration and accumulation are under consideration. The report also discusses radionuclide absorption/desorption distribution ratios and addresses adsorption/desorption mechanisms and their controlling processes for 25 elements under surface water conditions. These elements are: Am, Sb, C, Ce, Cm, Co, Cr, Cs, Eu, I, Fe, Mn, Np, P, Pu, Pm, Ra, Ru, Sr, Tc, Th, {sup 3}H, U, Zn and Zr.

  11. Modeling the Energy Use of a Connected and Automated Transportation...

    Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)

    ... from policies, including legality, insurance, liability, incentives, and treatment ... on Transportation: Energy Use, Greenhouse Gas Emissions, and Other Factors: http:...

  12. Reaction-based Transport Modeling of Iron Reduction and Uranium Immobilization at Area 2 of the NABIR Field Research Center

    SciTech Connect (OSTI)

    Tsyh Yeh, Gour

    2007-12-21

    This research sought to examine biogeochemical processes likely to take place in the less conductive materials above and below the gravel during the in situ ethanol biostimulation experiment conducted at Area 2 during 2005-2006. The in situ experiment in turn examined the hypothesis that injection of electron donor into this layer would induce formation of a redox barrier in the less conductive materials, resulting in decreased mass transfer of uranium out these materials and attendant declines in groundwater U(VI) concentration. Our project focuses on the development of a mechanistic understanding and quantitative models of coupled Fe(III)/U(VI) reduction in FRC Area 2 sediments. This report summarizes research activities conducted at The University of Central Florida (2004-2007), the development of biogeochemical and reactive transport models and the conduction of numerical simulations at laboratory, column, and field scales.

  13. Transportation Sector Model of the National Energy Modeling System. Volume 2 -- Appendices: Part 2

    SciTech Connect (OSTI)

    1998-01-01

    The attachments contained within this appendix provide additional details about the model development and estimation process which do not easily lend themselves to incorporation in the main body of the model documentation report. The information provided in these attachments is not integral to the understanding of the model`s operation, but provides the reader with opportunity to gain a deeper understanding of some of the model`s underlying assumptions. There will be a slight degree of replication of materials found elsewhere in the documentation, made unavoidable by the dictates of internal consistency. Each attachment is associated with a specific component of the transportation model; the presentation follows the same sequence of modules employed in Volume 1. The following attachments are contained in Appendix F: Fuel Economy Model (FEM)--provides a discussion of the FEM vehicle demand and performance by size class models; Alternative Fuel Vehicle (AFV) Model--describes data input sources and extrapolation methodologies; Light-Duty Vehicle (LDV) Stock Model--discusses the fuel economy gap estimation methodology; Light Duty Vehicle Fleet Model--presents the data development for business, utility, and government fleet vehicles; Light Commercial Truck Model--describes the stratification methodology and data sources employed in estimating the stock and performance of LCT`s; Air Travel Demand Model--presents the derivation of the demographic index, used to modify estimates of personal travel demand; and Airborne Emissions Model--describes the derivation of emissions factors used to associate transportation measures to levels of airborne emissions of several pollutants.

  14. Interfacial Structure and Reactivity | Argonne National Laboratory

    Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)

    a robust, molecular-scale understanding of its structure and reactivity? Research Context The transport of ions across the electrodeelectrolyte interface can lead to kinetic...

  15. A generic model for transport in turbulent shear flows

    SciTech Connect (OSTI)

    Newton, Andrew P. L.; Kim, Eun-Jin [Department of Applied Mathematics, University of Sheffield, Sheffield, S3 7RH (United Kingdom)

    2011-05-15

    Turbulence regulation by large-scale shear flows is crucial for a predictive modeling of transport in plasma. In this paper the suppression of turbulent transport by large-scale flows is studied numerically by measuring the turbulent diffusion D{sub t} and scalar amplitude of decaying passive scalar fields n{sup '} advected by various turbulent flows. Both uniform flows and shear flows are shown to suppress turbulence causing the quenching in transport and turbulence amplitude. The uniform flows U{sub 0}={Lambda}y with the advection rate {Lambda} in the case of a finite correlated forcing with {tau}{sub F}=1 gives rise to the advection/sweeping effect which suppresses D{sub t}, and as {proportional_to}{Lambda}{sup -2} for {Lambda}>>{tau}{sub F}{sup -1}. In contrast, no influence of the uniform flow is found in the case of a short correlated forcing {tau}{sub F}{yields}0 due to Galilean invariance. For the shear flow U{sub 0}={Omega}sinxy ({Omega}= constant shearing rate) with the appropriate choice of the forcing ({tau}{sub F}{yields}0) the nature of transport suppression is shown to crucially depend on the properties of the turbulence. Specifically, for prescribed turbulence with a short correlation time {tau}{sub c}={tau}{sub F}<<{Omega}{sup -1}, the turbulence statistics scale as D{sub t{proportional_to}{Omega}}{sup -0.02}, {proportional_to}{Omega}{sup -0.62} and cross-phase cos{theta}{proportional_to}{Omega}{sup 0.29}. For consistently evolved turbulence with a finite correlation time {tau}{sub c{>=}{Omega}}{sup -1}, turbulence statistics are suppressed more strongly as D{sub t{proportional_to}{Omega}}{sup -1.75}, {proportional_to}{Omega}{sup -2.41}, {proportional_to}{Omega}{sup -0.65} and <{omega}{sup '2}>{proportional_to}{Omega}{sup -0.50}. A novel renormalization scheme is then introduced to rescale our results into the regime within which the kinetic energy and enstrophy are unchanged by

  16. RADTRAD: A simplified model for RADionuclide Transport and Removal...

    Office of Scientific and Technical Information (OSTI)

    radioactive material due to containment sprays, natural deposition, filters, and other ... FILTERS; QUALITY ASSURANCE; SOURCE TERMS; SPRAYS; TRANSPORT; RADIATION DOSES; R CODES; ...

  17. Predicting tropospheric ozone and hydroxyl radical in a global, three-dimensional, chemistry, transport, and deposition model

    SciTech Connect (OSTI)

    Atherton, C.S.

    1995-01-05

    Two of the most important chemically reactive tropospheric gases are ozone (O{sub 3}) and the hydroxyl radical (OH). Although ozone in the stratosphere is a necessary protector against the sun`s radiation, tropospheric ozone is actually a pollutant which damages materials and vegetation, acts as a respiratory irritant, and is a greenhouse gas. One of the two main sources of ozone in the troposphere is photochemical production. The photochemistry is initiated when hydrocarbons and carbon monoxide (CO) react with nitrogen oxides (NO{sub x} = NO + NO{sub 2}) in the presence of sunlight. Reaction with the hydroxyl radical, OH, is the main sink for many tropospheric gases. The hydroxyl radical is highly reactive and has a lifetime on the order of seconds. Its formation is initiated by the photolysis of tropospheric ozone. Tropospheric chemistry involves a complex, non-linear set of chemical reactions between atmospheric species that vary substantially in time and space. To model these and other species on a global scale requires the use of a global, three-dimensional chemistry, transport, and deposition (CTD) model. In this work, I developed two such three dimensional CTD models. The first model incorporated the chemistry necessary to model tropospheric ozone production from the reactions of nitrogen oxides with carbon monoxide (CO) and methane (CH{sub 4}). The second also included longer-lived alkane species and the biogenic hydrocarbon isoprene, which is emitted by growing plants and trees. The models` ability to predict a number of key variables (including the concentration of O{sub 3}, OH, and other species) were evaluated. Then, several scenarios were simulated to understand the change in the chemistry of the troposphere since preindustrial times and the role of anthropogenic NO{sub x} on present day conditions.

  18. 2D numerical simulation of the MEP energy-transport model with a finite difference scheme

    SciTech Connect (OSTI)

    Romano, V. . E-mail: romano@dmi.unict.it

    2007-02-10

    A finite difference scheme of Scharfetter-Gummel type is used to simulate a consistent energy-transport model for electron transport in semiconductors devices, free of any fitting parameters, formulated on the basis of the maximum entropy principle. Simulations of silicon n{sup +}-n-n{sup +} diodes, 2D-MESFET and 2D-MOSFET and comparisons with the results obtained by a direct simulation of the Boltzmann transport equation and with other energy-transport models, known in the literature, show the validity of the model and the robustness of the numerical scheme.

  19. WIPP Documents - Transportation

    Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)

    Transportation

  20. Modeling of Diesel Exhaust Systems: A methodology to better simulate soot reactivity

    Broader source: Energy.gov [DOE]

    Discussed development of a methodology for creating accurate soot models for soot samples from various origins with minimal characterization

  1. Modeling two-fluid-phase flow and species transport in porous...

    Office of Scientific and Technical Information (OSTI)

    Modeling two-fluid-phase flow and species transport in porous media Citation Details In-Document Search This content will become publicly available on March 3, 2017 Title: Modeling ...

  2. Model Annex for Preparedness and Response to Radiological Transportation Incidents

    Broader source: Energy.gov [DOE]

    This part should contain a general statement of the intent of this Annex. To provide for the planning, preparedness and coordination of emergency service efforts to respond to a transportation...

  3. Water Transport in PEM Fuel Cells: Advanced Modeling, Material...

    Broader source: Energy.gov (indexed) [DOE]

    in PEM fuel cells, was given by CFDRC's J. Vernon Cole at a DOE fuel cell meeting in February 2007. newfccolecfd.pdf (976.38 KB) More Documents & Publications Water Transport in ...

  4. Boundary layer modeling of reactive flow over a porous surface with angled injection

    SciTech Connect (OSTI)

    Liu, Shiling; Fotache, Catalin G.; Hautman, Donald J.; Ochs, Stuart S. [United Technologies Research Center, MS 129-29, 411 Silver Lane, East Hartford, CT 06108 (United States); Chao, Beei-Huan [Department of Mechanical Engineering, University of Hawaii at Manoa, Honolulu, HI 96822 (United States)

    2008-08-15

    An analytical model was developed to investigate the dynamics of nonpremixed flames in a shear layer established between a mainstream flow of fuel-rich combustion products and a porous surface with an angled injection of air. In the model, a one-step overall chemical reaction was employed, together with boundary layer conservation equations solved using similarity solutions. Parametric studies were performed to understand the effects of equivalence ratio, temperature, and mass flow rate of the fuel and air streams on the flame standoff distance, surface temperature, and heat flux at the surface. The analytical model predictions were compared with computational fluid dynamics results obtained using the FLUENT commercial code for both the laminar and the turbulent flow models. Qualitative agreement in surface temperature was observed. Finally, the flame stability limits predicted by the model were compared with available experimental data and found to agree qualitatively, as well. (author)

  5. A mixture model for shock compression of porous multi-component reactive materials

    SciTech Connect (OSTI)

    Baer, M.R.

    1993-07-01

    A multiphase mixture model is presented to describe shock compression of highly porous, multi-component powders. Volume fractions are represented as independent kinematic variables and thermodynamically-admissible phase interaction constitutive models are formulated in the context of a three phase system. Numerical solutions of the multiphase flow equations simulated impact on a porous layer of mixed powders of Al and Fe{sub 2}O{sub 3}. The multiphase model predicts dispersive compaction waves which have features similar to observed time-resolved pressure measurements.

  6. A turbulent transport network model in MULTIFLUX coupled with TOUGH2

    SciTech Connect (OSTI)

    Danko, G.; Bahrami, D.; Birkholzer, J.T.

    2011-02-15

    A new numerical method is described for the fully iterated, conjugate solution of two discrete submodels, involving (a) a transport network model for heat, moisture, and airflows in a high-permeability, air-filled cavity; and (b) a variably saturated fractured porous medium. The transport network submodel is an integrated-parameter, computational fluid dynamics solver, describing the thermal-hydrologic transport processes in the flow channel system of the cavity with laminar or turbulent flow and convective heat and mass transport, using MULTIFLUX. The porous medium submodel, using TOUGH2, is a solver for the heat and mass transport in the fractured rock mass. The new model solution extends the application fields of TOUGH2 by integrating it with turbulent flow and transport in a discrete flow network system. We present demonstrational results for a nuclear waste repository application at Yucca Mountain with the most realistic model assumptions and input parameters including the geometrical layout of the nuclear spent fuel and waste with variable heat load for the individual containers. The MULTIFLUX and TOUGH2 model elements are fully iterated, applying a programmed reprocessing of the Numerical Transport Code Functionalization model-element in an automated Outside Balance Iteration loop. The natural, convective airflow field and the heat and mass transport in a representative emplacement drift during postclosure are explicitly solved in the new model. The results demonstrate that the direction and magnitude of the air circulation patterns and all transport modes are strongly affected by the heat and moisture transport processes in the surrounding rock, justifying the need for a coupled, fully iterated model solution such as the one presented in the paper.

  7. dfnWorks: A discrete fracture network framework for modeling subsurface flow and transport

    DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)

    Hyman, Jeffrey D.; Karra, Satish; Makedonska, Nataliia; Gable, Carl W.; Painter, Scott L.; Viswanathan, Hari S.

    2015-11-01

    DFNWORKS is a parallelized computational suite to generate three-dimensional discrete fracture networks (DFN) and simulate flow and transport. Developed at Los Alamos National Laboratory over the past five years, it has been used to study flow and transport in fractured media at scales ranging from millimeters to kilometers. The networks are created and meshed using DFNGEN, which combines FRAM (the feature rejection algorithm for meshing) methodology to stochastically generate three-dimensional DFNs with the LaGriT meshing toolbox to create a high-quality computational mesh representation. The representation produces a conforming Delaunay triangulation suitable for high performance computing finite volume solvers in anmore » intrinsically parallel fashion. Flow through the network is simulated in dfnFlow, which utilizes the massively parallel subsurface flow and reactive transport finite volume code PFLOTRAN. A Lagrangian approach to simulating transport through the DFN is adopted within DFNTRANS to determine pathlines and solute transport through the DFN. Example applications of this suite in the areas of nuclear waste repository science, hydraulic fracturing and CO2 sequestration are also included.« less

  8. dfnWorks: A discrete fracture network framework for modeling subsurface flow and transport

    SciTech Connect (OSTI)

    Hyman, Jeffrey D.; Karra, Satish; Makedonska, Nataliia; Gable, Carl W.; Painter, Scott L.; Viswanathan, Hari S.

    2015-11-01

    DFNWORKS is a parallelized computational suite to generate three-dimensional discrete fracture networks (DFN) and simulate flow and transport. Developed at Los Alamos National Laboratory over the past five years, it has been used to study flow and transport in fractured media at scales ranging from millimeters to kilometers. The networks are created and meshed using DFNGEN, which combines FRAM (the feature rejection algorithm for meshing) methodology to stochastically generate three-dimensional DFNs with the LaGriT meshing toolbox to create a high-quality computational mesh representation. The representation produces a conforming Delaunay triangulation suitable for high performance computing finite volume solvers in an intrinsically parallel fashion. Flow through the network is simulated in dfnFlow, which utilizes the massively parallel subsurface flow and reactive transport finite volume code PFLOTRAN. A Lagrangian approach to simulating transport through the DFN is adopted within DFNTRANS to determine pathlines and solute transport through the DFN. Example applications of this suite in the areas of nuclear waste repository science, hydraulic fracturing and CO2 sequestration are also included.

  9. Report of the 7th Meeting of the Transport Modeling Working Group

    Office of Energy Efficiency and Renewable Energy (EERE) Indexed Site

    7 th Meeting of the Transport Modeling Working Group May 22 and 23, 2014 The Transport Modeling Working Group (TMWG) met over 1.5 days at the end of May at Lawrence Berkeley National Laboratory. The chosen topics were to discuss multiphase flow with a focus on the diffusion media and the nexus between durability and transport modeling. The format was similar to previous meetings with ample time for discussion led off by a few presentations pertinent to the topic at hand. Both core TMWG

  10. NREL: Transportation Research - Vehicle Thermal Management Models and Tools

    Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)

    Models and Tools image of three models of semi truck cabs. Truck cab models drawn from CAD geometry using CoolCalc (left and center), and a model with overlay of computational fluid dynamics flow (right) indicate areas of heat absorption and loss. Illustrations by Jason Lustbader, Matt Jeffers, and Larry Chaney, NREL The National Renewable Energy Laboratory's (NREL's) vehicle thermal management modeling tools allow researchers to assess the trade-offs and calculate the potential benefits of

  11. Underground Test Area Subproject Phase I Data Analysis Task. Volume VII - Tritium Transport Model Documentation Package

    SciTech Connect (OSTI)

    1996-12-01

    Volume VII of the documentation for the Phase I Data Analysis Task performed in support of the current Regional Flow Model, Transport Model, and Risk Assessment for the Nevada Test Site Underground Test Area Subproject contains the tritium transport model documentation. Because of the size and complexity of the model area, a considerable quantity of data was collected and analyzed in support of the modeling efforts. The data analysis task was consequently broken into eight subtasks, and descriptions of each subtask's activities are contained in one of the eight volumes that comprise the Phase I Data Analysis Documentation.

  12. Simple Model Representations of Transport in a Complex Fracture...

    Office of Scientific and Technical Information (OSTI)

    It is common, however, to represent the complex fracture by much simpler models consisting ... Simple-model properties are often inferred from the analysis of short-term (one to a few ...

  13. Remarks on a model of thermal transport in nanofluids (Journal...

    Office of Scientific and Technical Information (OSTI)

    of Publication: United States Language: English Subject: 75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND SUPERFLUIDITY; HEAT TRANSFER; FLOW MODELS; FLUIDS; NANOSTRUCTURES...

  14. Modeling of the electronic transport in multijunction solar cells

    SciTech Connect (OSTI)

    Rau, U.; Goldbach, M.

    1994-12-31

    Simulations of the electrical transport in multijunction thin-film solar cells made from polycrystalline silicon are presented. The authors investigate the effect of the grain size on the efficiency of the multijunction solar cell. Here, they concentrate on micro crystalline material with a high recombination velocity at the grain boundaries of 10{sup 4}cm/s. Typical results of their calculations demonstrate that based on the multijunction design structure consisting of 8 or more layers efficiencies of 14% may be obtained from 12--20 {micro}m thick solar cells.

  15. Streamtube Fate and Transport Modeling of the Source Term for the Old Radioactive Waste

    SciTech Connect (OSTI)

    Brewer, K.

    2000-11-16

    The modeling described in this report is an extension of previous fate and transport modeling for the Old Radioactive Waste Burial Ground Corrective Measures Study/Feasibility Study. The purpose of this and the previous modeling is to provide quantitative input to the screening of remedial alternatives for the CMS/FS for this site.

  16. DFNWorks. A discrete fracture network framework for modeling subsurface flow and transport

    DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)

    Hyman, Jeffrey D.; Karra, Satish; Makedonska, Nataliia; Gable, Carl W.; Painter, Scott L.; Viswanathan, Hari S.

    2015-08-10

    DFNWorks is a parallalized computational suite to generate three-dimensional discrete fracture networks (DFN) and simulate flow and transport. Developed at Los Alamos National Laboratory over the past five years, it has been used to study flow and transport in fractured media at scales ranging from millimeters to kilometers. The networks are created and meshed using dfnGen, which combines fram (the feature rejection algorithm for meshing) methodology to stochastically generate three-dimensional DFNs on the basis of site specific data with the LaGriT meshing toolbox to create a high-quality computational mesh representation, specifically a conforming Delaunay triangulation suitable for high performance computingmore » finite volume solvers, of the DFN in an intrinsically parallel fashion. Flow through the network is simulated in dfnFlow, which utilizes the massively parallel subsurface flow and reactive transport finite volume code pflotran. A Lagrangian approach to simulating transport through the DFN is adopted within dfnTrans, which is an extension of the walkabout particle tracking method to determine pathlines through the DFN. Example applications of this suite in the areas of nuclear waste repository science, hydraulic fracturing and CO2 sequestration are also included.« less

  17. DFNWorks. A discrete fracture network framework for modeling subsurface flow and transport

    SciTech Connect (OSTI)

    Hyman, Jeffrey D.; Karra, Satish; Makedonska, Nataliia; Gable, Carl W.; Painter, Scott L.; Viswanathan, Hari S.

    2015-08-10

    DFNWorks is a parallalized computational suite to generate three-dimensional discrete fracture networks (DFN) and simulate flow and transport. Developed at Los Alamos National Laboratory over the past five years, it has been used to study flow and transport in fractured media at scales ranging from millimeters to kilometers. The networks are created and meshed using dfnGen, which combines fram (the feature rejection algorithm for meshing) methodology to stochastically generate three-dimensional DFNs on the basis of site specific data with the LaGriT meshing toolbox to create a high-quality computational mesh representation, specifically a conforming Delaunay triangulation suitable for high performance computing finite volume solvers, of the DFN in an intrinsically parallel fashion. Flow through the network is simulated in dfnFlow, which utilizes the massively parallel subsurface flow and reactive transport finite volume code pflotran. A Lagrangian approach to simulating transport through the DFN is adopted within dfnTrans, which is an extension of the walkabout particle tracking method to determine pathlines through the DFN. Example applications of this suite in the areas of nuclear waste repository science, hydraulic fracturing and CO2 sequestration are also included.

  18. A review of wind field models for atmospheric transport

    SciTech Connect (OSTI)

    Ramsdell, J.V. Jr.; Skyllingstad, E.D.

    1993-06-01

    The primary objective of the Hanford Environmental Dose Reconstruction (HEDR) Project is to estimate the radiation dose that individuals could have received as a result of emissions since 1944 from the US Department of Energy`s (DOE) Hanford Site near Richland, Washington. The HEDR Project is developing a computer code to estimate these doses and their uncertainties. The code, known as the HEDR integrated Code (HEDRIC), consists of four separate component codes. One of the component codes, called the Regional Atmospheric Transport Code for Hanford Emission Tracking (RATCHET) combines meteorological and release data to estimate time-integrated air concentrations and surface contamination at specific locations in the vicinity of the Hanford Site. The RATCHET domain covers approximately 75,000 square miles, extending from the crest of the Cascade Mountains on the west to the eastern edge of the Idaho panhandle and from central Oregon on the south to the Canadian border. This letter report explains the procedures in RATCHET that transform observed wind data into the wind fields used in atmospheric transport calculations. It also describes and evaluates alternative procedures not selected for use in RATCHET.

  19. Centrifuge modeling of LNAPL transport in partially saturated sand

    SciTech Connect (OSTI)

    Esposito, G.; Allersma, H.G.B.; Selvadurai, A.P.S.

    1999-12-01

    Model tests were performed at the Geotechnical Centrifuge Facility of Delft University of Technology, The Netherlands, to examine the mechanics of light nonaqueous phase liquid (LNAPL) movement in a partially saturated porous granular medium. The experiment simulated a 2D spill of LNAPL in an unsaturated sand prepared at two values of porosity. The duration of the centrifuge model tests corresponded to a prototype equivalent of 110 days. The choice of modeling a 2D flow together with the use of a transparent container enabled direct visual observation of the experiments. Scaling laws developed in connection with other centrifuge modeling studies were used to support the test results. Tests were conducted at two different centrifuge accelerations to verify, by means of the modeling of models technique, the similitude between the different experiments. The paper presents details of the experimental methodologies and the measuring techniques used to evaluate the final distribution of water and LNAPL content in the soils.

  20. Water Transport in PEM Fuel Cells: Advanced Modeling, Material...

    Office of Energy Efficiency and Renewable Energy (EERE) Indexed Site

    in PEM Fuel Cells: Advanced Modeling, Material Selection, Testing and Design Optimization ... Optimization Part of a 100 million fuel cell award announced by DOE Secretary Bodman on ...

  1. Lessons Learned from Alternative Transportation Fuels: Modeling Transition Dynamics

    SciTech Connect (OSTI)

    Welch, C.

    2006-02-01

    Report focuses on understanding how analytical system modeling and data from AFV experiences could improve our understanding of the dynamic forces governing the transition to a hydrogen future.

  2. Mass-transport models to predict toxicity of inhaled gases in the upper respiratory tract

    SciTech Connect (OSTI)

    Hubal, E.A.C.; Fedkiw, P.S.; Kimbell, J.S. [North Carolina State Univ., Raleigh, NC (United States)

    1996-04-01

    Mass-transport (the movement of a chemical species) plays an important role in determining toxic responses of the upper respiratory tract (URT) to inhaled chemicals. Mathematical dosimetry models incorporate physical characteristics of mass transport and are used to predict quantitative uptake (absorption rate) and distribution of inhaled gases and vapors in the respiratory tract. Because knowledge of dose is an essential component of quantitative risk assessment, dosimetry modeling plays an important role in extrapolation of animal study results to humans. A survey of existing mathematical dosimetry models for the URT is presented, limitations of current models are discussed, and adaptations of existing models to produce a generally applicable model are suggested. Reviewed URT dosimetry models are categorized as early, lumped-parameter, and distributed-parameter models. Specific examples of other relevant modeling work are also presented. 35 refs., 11 figs., 1 tab.

  3. Transportation Research

    Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)

    transportation-research TRACC RESEARCH Computational Fluid Dynamics Computational Structural Mechanics Transportation Systems Modeling Transportation Research Current Research Overview The U.S. Department of Transportation (USDOT) has established its only high-performance computing and engineering analysis research facility at Argonne National Laboratory to provide applications support in key areas of applied research and development for the USDOT community. The Transportation Research and

  4. Modeling the effect of gas transport on the formation of defects during thermolysis of powder moldings

    SciTech Connect (OSTI)

    Song, J.H.; Edirisinghe, M.J.; Evans, J.R.; Twizell, E.H.

    1996-04-01

    The removal of binder from ceramic or metal moldings by thermolysis involves the transport of degradation products through the parent organic phase and the vacated porous body. A numerical model has been developed to combine an equation which takes into account different gas-flow regimes with an equation for the transport of organic molecules in molten polymers. Computer modeling reveals the critical heating rate above which defects occur due to boiling of the polymer-monomer solution at the center of the molding. The situation in which a porous outer layer of the molding develops, offering resistance to flow of the evolved monomer gas, is then treated. This gives rise to a moving boundary with a variable concentration of diffusant which is dependent on the surface flux, gas transport coefficient and thickness of the porous layer. The contributions of diffusion and viscous flow to gas transport are considered. {copyright} {ital 1996 Materials Research Society.}

  5. A SCOPING STUDY: Development of Probabilistic Risk Assessment Models for Reactivity Insertion Accidents During Shutdown In U.S. Commercial Light Water Reactors

    SciTech Connect (OSTI)

    S. Khericha

    2011-06-01

    This report documents the scoping study of developing generic simplified fuel damage risk models for quantitative analysis from inadvertent reactivity insertion events during shutdown (SD) in light water pressurized and boiling water reactors. In the past, nuclear fuel reactivity accidents have been analyzed both mainly deterministically and probabilistically for at-power and SD operations of nuclear power plants (NPPs). Since then, many NPPs had power up-rates and longer refueling intervals, which resulted in fuel configurations that may potentially respond differently (in an undesirable way) to reactivity accidents. Also, as shown in a recent event, several inadvertent operator actions caused potential nuclear fuel reactivity insertion accident during SD operations. The set inadvertent operator actions are likely to be plant- and operation-state specific and could lead to accident sequences. This study is an outcome of the concern which arose after the inadvertent withdrawal of control rods at Dresden Unit 3 in 2008 due to operator actions in the plant inadvertently three control rods were withdrawn from the reactor without knowledge of the main control room operator. The purpose of this Standardized Plant Analysis Risk (SPAR) Model development project is to develop simplified SPAR Models that can be used by staff analysts to perform risk analyses of operating events and/or conditions occurring during SD operation. These types of accident scenarios are dominated by the operator actions, (e.g., misalignment of valves, failure to follow procedures and errors of commissions). Human error probabilities specific to this model were assessed using the methodology developed for SPAR model human error evaluations. The event trees, fault trees, basic event data and data sources for the model are provided in the report. The end state is defined as the reactor becomes critical. The scoping study includes a brief literature search/review of historical events, developments of

  6. Modeling electron transport in the presence of electric and magnetic fields.

    SciTech Connect (OSTI)

    Fan, Wesley C.; Drumm, Clifton Russell; Pautz, Shawn D.; Turner, C. David

    2013-09-01

    This report describes the theoretical background on modeling electron transport in the presence of electric and magnetic fields by incorporating the effects of the Lorentz force on electron motion into the Boltzmann transport equation. Electromagnetic fields alter the electron energy and trajectory continuously, and these effects can be characterized mathematically by differential operators in terms of electron energy and direction. Numerical solution techniques, based on the discrete-ordinates and finite-element methods, are developed and implemented in an existing radiation transport code, SCEPTRE.

  7. Semiclassical matrix model for quantum chaotic transport with time-reversal symmetry

    SciTech Connect (OSTI)

    Novaes, Marcel

    2015-10-15

    We show that the semiclassical approach to chaotic quantum transport in the presence of time-reversal symmetry can be described by a matrix model. In other words, we construct a matrix integral whose perturbative expansion satisfies the semiclassical diagrammatic rules for the calculation of transport statistics. One of the virtues of this approach is that it leads very naturally to the semiclassical derivation of universal predictions from random matrix theory.

  8. Baseline requirements of the proposed action for the Transportation Management Division routing models

    SciTech Connect (OSTI)

    Johnson, P.E.; Joy, D.S.

    1995-02-01

    The potential impacts associated with the transportation of hazardous materials are important to shippers, carriers, and the general public. This is particularly true for shipments of radioactive material. The shippers are primarily concerned with safety, security, efficiency, and equipment requirements. The carriers are concerned with the potential impact that radioactive shipments may have on their operations--particularly if such materials are involved in an accident. The general public has also expressed concerns regarding the safety of transporting radioactive and other hazardous materials through their communities. Because transportation routes are a central concern in hazardous material transport, the prediction of likely routes is the first step toward resolution of these issues. In response to these routing needs, several models have been developed over the past fifteen years at Oak Ridge National Laboratory (ORNL). The HIGHWAY routing model is used to predict routes for truck transportation, the INTERLINE routing model is used to predict both rail and barge routes, and the AIRPORT locator model is used to determine airports with specified criteria near a specific location. As part of the ongoing improvement of the US Department of Energy`s (DOE) Environmental Management Transportation Management Division`s (EM-261) computer systems and development efforts, a Baseline Requirements Assessment Session on the HIGHWAY, INTERLINE, and AIRPORT models was held at ORNL on April 27, 1994. The purpose of this meeting was to discuss the existing capabilities of the models and data bases and to review enhancements of the models and data bases to expand their usefulness. The results of the Baseline Requirements Assessment Section will be discussed in this report. The discussions pertaining to the different models are contained in separate sections.

  9. Temperature, humidity and air flow in the emplacement drifts using convection and dispersion transport models

    SciTech Connect (OSTI)

    Danko, G.; Birkholzer, J.T.; Bahrami, D.; Halecky, N.

    2009-10-01

    A coupled thermal-hydrologic-airflow model is developed, solving for the transport processes within a waste emplacement drift and the surrounding rockmass together at the proposed nuclear waste repository at Yucca Mountain. Natural, convective air flow as well as heat and mass transport in a representative emplacement drift during post-closure are explicitly simulated, using the MULTIFLUX model. The conjugate, thermal-hydrologic transport processes in the rockmass are solved with the TOUGH2 porous-media simulator in a coupled way to the in-drift processes. The new simulation results show that large-eddy turbulent flow, as opposed to small-eddy flow, dominate the drift air space for at least 5000 years following waste emplacement. The size of the largest, longitudinal eddy is equal to half of the drift length, providing a strong axial heat and moisture transport mechanism from the hot to the cold drift sections. The in-drift results are compared to those from simplified models using a surrogate, dispersive model with an equivalent dispersion coefficient for heat and moisture transport. Results from the explicit, convective velocity simulation model provide higher axial heat and moisture fluxes than those estimated from the previously published, simpler, equivalent-dispersion models, in addition to showing differences in temperature, humidity and condensation rate distributions along the drift length. A new dispersive model is also formulated, giving a time- and location-variable function that runs generally about ten times higher in value than the highest dispersion coefficient currently used in the Yucca Mountain Project as an estimate for the equivalent dispersion coefficient in the emplacement drift. The new dispersion coefficient variation, back-calculated from the convective model, can adequately describe the heat and mass transport processes in the emplacement drift example.

  10. Simulation model for oil slick transport in lakes

    SciTech Connect (OSTI)

    Shen, H.T.; Yapa, P.D.; Petroski, M.E.

    1987-10-01

    A computer model for simulating oil slick movement in lakes by a Lagrangian discrete parcel algorithm is presented. In this model the transformation of an oil slick due to advection, spreading, evaporation, and dissolution is considered. For open water conditions the movement of the oil slick by water current and wind is considered using the drifting factor formulation. For ice-covered conditions the drift velocity is determined according to the ice roughness and current velocity. The current distribution in the lake is determined by a rigid lid circulation model. In the spreading process the mechanical spreading of the oil slick due to the balance in inertia, gravity, viscous, and surface tension forces is considered, in addition to the dispersion of the surface oil layer. Boundary conditions along the shore are formulated according to the storage capacity of the shoreline. The model can be used for simulating either instantaneous or continuous oil spills. Sample simulations for oil spills in Lake St. Clair are presented.

  11. Phase II Contaminant Transport Parameters for the Groundwater Flow and Contaminant Transport Model of Corrective Action Unit 98: Frenchman Flat, Nye County, Nevada, Rev. No.: 0

    SciTech Connect (OSTI)

    DeNovio, Nicole M.; Bryant, Nathan; King, Chrissi B.; Bhark, Eric; Drellack, Sigmund L.; Pickens, John F.; Farnham, Irene; Brooks, Keely M.; Reimus, Paul; Aly, Alaa

    2005-04-01

    This report documents pertinent transport data and data analyses as part of the Phase II Corrective Action Investigation (CAI) for Frenchman Flat (FF) Corrective Action Unit (CAU) 98. The purpose of this data compilation and related analyses is to provide the primary reference to support parameterization of the Phase II FF CAU transport model.

  12. NREL: Transportation Research - NREL's Battery Life Predictive Model Helps

    Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)

    Companies Take Charge NREL's Battery Life Predictive Model Helps Companies Take Charge October 26, 2015 A series of batteries hooked together next to a monitor. An example of a stationary, grid-connected battery is the NREL project from Erigo/EaglePicher Technologies, LLC Technologies. Inverters and nickel cadmium batteries inside of a utility scale 300 kW battery storage system will support Department of Defense micro-grids. Photo by Dennis Schroeder / NREL 32696 Companies that rely on

  13. BLT-EC (Breach, Leach and Transport-Equilibrium Chemistry) data input guide. A computer model for simulating release and coupled geochemical transport of contaminants from a subsurface disposal facility

    SciTech Connect (OSTI)

    MacKinnon, R.J.; Sullivan, T.M.; Kinsey, R.R.

    1997-05-01

    The BLT-EC computer code has been developed, implemented, and tested. BLT-EC is a two-dimensional finite element computer code capable of simulating the time-dependent release and reactive transport of aqueous phase species in a subsurface soil system. BLT-EC contains models to simulate the processes (container degradation, waste-form performance, transport, chemical reactions, and radioactive production and decay) most relevant to estimating the release and transport of contaminants from a subsurface disposal system. Water flow is provided through tabular input or auxiliary files. Container degradation considers localized failure due to pitting corrosion and general failure due to uniform surface degradation processes. Waste-form performance considers release to be limited by one of four mechanisms: rinse with partitioning, diffusion, uniform surface degradation, and solubility. Transport considers the processes of advection, dispersion, diffusion, chemical reaction, radioactive production and decay, and sources (waste form releases). Chemical reactions accounted for include complexation, sorption, dissolution-precipitation, oxidation-reduction, and ion exchange. Radioactive production and decay in the waste form is simulated. To improve the usefulness of BLT-EC, a pre-processor, ECIN, which assists in the creation of chemistry input files, and a post-processor, BLTPLOT, which provides a visual display of the data have been developed. BLT-EC also includes an extensive database of thermodynamic data that is also accessible to ECIN. This document reviews the models implemented in BLT-EC and serves as a guide to creating input files and applying BLT-EC.

  14. Conceptual model for regional radionuclide transport from a salt dome repository: a technical memorandum

    SciTech Connect (OSTI)

    Kier, R.S.; Showalter, P.A.; Dettinger, M.D.

    1980-05-30

    Disposal of high-level radioactive wastes is a major environmental problem influencing further development of nuclear energy in this country. Salt domes in the Gulf Coast Basin are being investigated as repository sites. A major concern is geologic and hydrologic stability of candidate domes and potential transport of radionuclides by groundwater to the biosphere prior to their degradation to harmless levels of activity. This report conceptualizes a regional geohydrologic model for transport of radionuclides from a salt dome repository. The model considers transport pathways and the physical and chemical changes that would occur through time prior to the radionuclides reaching the biosphere. Necessary, but unknown inputs to the regional model involve entry and movement of fluids through the repository dome and across the dome-country rock interface and the effect on the dome and surrounding strata of heat generated by the radioactive wastes.

  15. Collisional particle-in-cell modeling for energy transport accompanied by atomic processes in dense plasmas

    SciTech Connect (OSTI)

    Mishra, R.; Beg, F. N.; Leblanc, P.; Sentoku, Y.; Wei, M. S.

    2013-07-15

    Fully relativistic collisional Particle-in-Cell (PIC) code, PICLS, has been developed to study extreme energy density conditions produced in intense laser-solid interaction. Recent extensions to PICLS, such as the implementation of dynamic ionization, binary collisions in a partially ionized plasma, and radiative losses, enhance the efficacy of simulating intense laser plasma interaction and subsequent energy transport in resistive media. Different ionization models are introduced and benchmarked against each other to check the suitability of the model. The atomic physics models are critical to determine the energy deposition and transport in dense plasmas, especially when they consist of high Z (atomic number) materials. Finally we demonstrate the electron transport simulations to show the importance of target material on fast electron dynamics.

  16. Contaminant Transport Parameters for the Groundwater Flow and Contaminant Transport Model of Corrective Action Units 101 and 102: Central and Western Pahute Mesa, Nye County, Nevada, Revision 0

    SciTech Connect (OSTI)

    Drici, Warda

    2003-08-01

    This report documents the analysis of the available transport parameter data conducted in support of the development of a Corrective Action Unit (CAU) groundwater flow model for Central and Western Pahute Mesa: CAUs 101 and 102.

  17. Carlsbad Area Office unveils full-scale model of new WIPP waste transportation cask

    Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)

    Carlsbad Area Office Unveils Full-Scale Model Of New WIPP Waste Transportation Cask CARLSBAD, N.M., February 23, 2000 - The U.S. Department of Energy's (DOE) Carlsbad Area Office today unveiled a full-scale model of its newest waste transportation cask, the RH-72B, during a ceremony at the local DOE offices. "This is another milestone for the Department of Energy," said Dr. Inés Triay, Manager of the Carlsbad Area Office, describing the importance of the new container for those

  18. Transported PDF Modeling of Nonpremixed Turbulent CO/H-2/N-2 Jet Flames

    SciTech Connect (OSTI)

    Zhao, xinyu; Haworth, D. C.; Huckaby, E. David

    2012-01-01

    Turbulent CO/H{sub 2}/N{sub 2} (syngas) flames are simulated using a transported composition probability density function (PDF) method. A consistent hybrid Lagrangian particle/Eulerian mesh algorithm is used to solve the modeled PDF transport equation. The model includes standard k? turbulence, gradient transport for scalars, and Euclidean minimum spanning tree (EMST) mixing. Sensitivities of model results to variations in the turbulence model, the treatment of radiation heat transfer, the choice of chemical mechanism, and the PDF mixing model are explored. A baseline model reproduces the measured mean and rms temperature, major species, and minor species profiles reasonably well, and captures the scaling that is observed in the experiments. Both our results and the literature suggest that further improvements can be realized with adjustments in the turbulence model, the radiation heat transfer model, and the chemical mechanism. Although radiation effects are relatively small in these flames, consideration of radiation is important for accurate NO prediction. Chemical mechanisms that have been developed specifically for fuels with high concentrations of CO and H{sub 2} perform better than a methane mechanism that was not designed for this purpose. It is important to account explicitly for turbulencechemistry interactions, although the details of the mixing model do not make a large difference in the results, within reasonable limits.

  19. Water-Level Data Analysis for the Saturated Zone Site-Scale Flow and Transport Model

    SciTech Connect (OSTI)

    P. Tucci

    2001-12-20

    This Analysis/Model Report (AMR) documents an updated analysis of water-level data performed to provide the saturated-zone, site-scale flow and transport model (CRWMS M&O 2000) with the configuration of the potentiometric surface, target water-level data, and hydraulic gradients for model calibration. The previous analysis was presented in ANL-NBS-HS-000034, Rev 00 ICN 01, Water-Level Data Analysis for the Saturated Zone Site-Scale Flow and Transport Model (USGS 2001). This analysis is designed to use updated water-level data as the basis for estimating water-level altitudes and the potentiometric surface in the SZ site-scale flow and transport model domain. The objectives of this revision are to develop computer files containing (1) water-level data within the model area (DTN: GS010908312332.002), (2) a table of known vertical head differences (DTN: GS0109083 12332.003), and (3) a potentiometric-surface map (DTN: GS010608312332.001) using an alternate concept from that presented in ANL-NBS-HS-000034, Rev 00 ICN 01 for the area north of Yucca Mountain. The updated water-level data include data obtained from the Nye County Early Warning Drilling Program (EWDP) and data from borehole USW WT-24. In addition to being utilized by the SZ site-scale flow and transport model, the water-level data and potentiometric-surface map contained within this report will be available to other government agencies and water users for ground-water management purposes. The potentiometric surface defines an upper boundary of the site-scale flow model, as well as provides information useful to estimation of the magnitude and direction of lateral ground-water flow within the flow system. Therefore, the analysis documented in this revision is important to SZ flow and transport calculations in support of total system performance assessment.

  20. Gray and multigroup radiation transport models for two-dimensional binary stochastic media using effective opacities

    DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)

    Olson, Gordon L.

    2015-09-24

    One-dimensional models for the transport of radiation through binary stochastic media do not work in multi-dimensions. In addition, authors have attempted to modify or extend the 1D models to work in multidimensions without success. Analytic one-dimensional models are successful in 1D only when assuming greatly simplified physics. State of the art theories for stochastic media radiation transport do not address multi-dimensions and temperature-dependent physics coefficients. Here, the concept of effective opacities and effective heat capacities is found to well represent the ensemble averaged transport solutions in cases with gray or multigroup temperature-dependent opacities and constant or temperature-dependent heat capacities. Inmore » every case analyzed here, effective physics coefficients fit the transport solutions over a useful range of parameter space. The transport equation is solved with the spherical harmonics method with angle orders of n=1 and 5. Although the details depend on what order of solution is used, the general results are similar, independent of angular order.« less

  1. Gray and multigroup radiation transport models for two-dimensional binary stochastic media using effective opacities

    SciTech Connect (OSTI)

    Olson, Gordon L.

    2015-09-24

    One-dimensional models for the transport of radiation through binary stochastic media do not work in multi-dimensions. In addition, authors have attempted to modify or extend the 1D models to work in multidimensions without success. Analytic one-dimensional models are successful in 1D only when assuming greatly simplified physics. State of the art theories for stochastic media radiation transport do not address multi-dimensions and temperature-dependent physics coefficients. Here, the concept of effective opacities and effective heat capacities is found to well represent the ensemble averaged transport solutions in cases with gray or multigroup temperature-dependent opacities and constant or temperature-dependent heat capacities. In every case analyzed here, effective physics coefficients fit the transport solutions over a useful range of parameter space. The transport equation is solved with the spherical harmonics method with angle orders of n=1 and 5. Although the details depend on what order of solution is used, the general results are similar, independent of angular order.

  2. Nondestructive determination of plutonium mass in spent fuel: prelliminary modeling results using the passive neutron Albedo reactivity technique

    SciTech Connect (OSTI)

    Evans, Louise G; Tobin, Stephen J; Schear, Melissa A; Menlove, Howard O; Lee, Sang Y; Swinhoe, Martyn T

    2009-01-01

    There are a variety of motivations for quantifying plutonium (Pu) in spent fuel assemblies by means of nondestructive assay (NDA) including the following: strengthening the capability of the International Atomic Energy Agency (LAEA) to safeguard nuclear facilities, quantifying shipper/receiver difference, determining the input accountability value at pyrochemical processing facilities, providing quantitative input to burnup credit and final safeguards measurements at a long-term repository. In order to determine Pu mass in spent fuel assemblies, thirteen NDA techniques were identified that provide information about the composition of an assembly. A key motivation of the present research is the realization that none of these techniques, in isolation, is capable of both (1) quantifying the Pu mass of an assembly and (2) detecting the diversion of a significant number of rods. It is therefore anticipated that a combination of techniques will be required. A 5 year effort funded by the Next Generation Safeguards Initiative (NGSI) of the U.S. DOE was recently started in pursuit of these goals. The first two years involves researching all thirteen techniques using Monte Carlo modeling while the final three years involves fabricating hardware and measuring spent fuel. Here, we present the work in two main parts: (1) an overview of this NGSI effort describing the motivations and approach being taken; (2) The preliminary results for one of the NDA techniques - Passive Neutron Albedo Reactivity (PNAR). The PNAR technique functions by using the intrinsic neutron emission of the fuel (primarily from the spontaneous fission of curium) to self-interrogate any fissile material present. Two separate measurements of the spent fuel are made, both with and without cadmium (Cd) present. The ratios of the Singles, Doubles and Triples count rates obtained in each case are analyzed; known as the Cd ratio. The primary differences between the two measurements are the neutron energy spectrum

  3. A POD reduced order model for resolving angular direction in neutron/photon transport problems

    SciTech Connect (OSTI)

    Buchan, A.G.; Calloo, A.A.; Goffin, M.G.; Dargaville, S.; Fang, F.; Pain, C.C.; Navon, I.M.

    2015-09-01

    This article presents the first Reduced Order Model (ROM) that efficiently resolves the angular dimension of the time independent, mono-energetic Boltzmann Transport Equation (BTE). It is based on Proper Orthogonal Decomposition (POD) and uses the method of snapshots to form optimal basis functions for resolving the direction of particle travel in neutron/photon transport problems. A unique element of this work is that the snapshots are formed from the vector of angular coefficients relating to a high resolution expansion of the BTE's angular dimension. In addition, the individual snapshots are not recorded through time, as in standard POD, but instead they are recorded through space. In essence this work swaps the roles of the dimensions space and time in standard POD methods, with angle and space respectively. It is shown here how the POD model can be formed from the POD basis functions in a highly efficient manner. The model is then applied to two radiation problems; one involving the transport of radiation through a shield and the other through an infinite array of pins. Both problems are selected for their complex angular flux solutions in order to provide an appropriate demonstration of the model's capabilities. It is shown that the POD model can resolve these fluxes efficiently and accurately. In comparison to high resolution models this POD model can reduce the size of a problem by up to two orders of magnitude without compromising accuracy. Solving times are also reduced by similar factors.

  4. User's manual for the Sandia Waste-Isolation Flow and Transport model (SWIFT).

    SciTech Connect (OSTI)

    Reeves, Mark; Cranwell, Robert M.

    1981-11-01

    This report describes a three-dimensional finite-difference model (SWIFT) which is used to simulate flow and transport processes in geologic media. The model was developed for use by the Nuclear Regulatory Commission in the analysis of deep geologic nuclear waste-disposal facilities. This document, as indicated by the title, is a user's manual and is intended to facilitate the use of the SWIFT simulator. Mathematical equations, submodels, application notes, and a description of the program itself are given herein. In addition, a complete input data guide is given along with several appendices which are helpful in setting up a data-input deck. Computer code SWIFT (Sandia Waste Isolation, Flow and Transport Model) is a fully transient, three-dimensional model which solves the coupled equations for transport in geologic media. The processes considered are: (1) fluid flow; (2) heat transport; (3) dominant-species miscible displacement; and (4) trace-species miscible displacement. The first three processes are coupled via fluid density and viscosity. Together they provide the velocity field on which the fourth process depends.

  5. Uncertainty Analysis Framework - Hanford Site-Wide Groundwater Flow and Transport Model

    SciTech Connect (OSTI)

    Cole, Charles R.; Bergeron, Marcel P.; Murray, Christopher J.; Thorne, Paul D.; Wurstner, Signe K.; Rogers, Phillip M.

    2001-11-09

    Pacific Northwest National Laboratory (PNNL) embarked on a new initiative to strengthen the technical defensibility of the predictions being made with a site-wide groundwater flow and transport model at the U.S. Department of Energy Hanford Site in southeastern Washington State. In FY 2000, the focus of the initiative was on the characterization of major uncertainties in the current conceptual model that would affect model predictions. The long-term goals of the initiative are the development and implementation of an uncertainty estimation methodology in future assessments and analyses using the site-wide model. This report focuses on the development and implementation of an uncertainty analysis framework.

  6. A DOUBLE-RING ALGORITHM FOR MODELING SOLAR ACTIVE REGIONS: UNIFYING KINEMATIC DYNAMO MODELS AND SURFACE FLUX-TRANSPORT SIMULATIONS

    SciTech Connect (OSTI)

    Munoz-Jaramillo, Andres; Martens, Petrus C. H.; Nandy, Dibyendu; Yeates, Anthony R. E-mail: dnandi@iiserkol.ac.i E-mail: anthony@maths.dundee.ac.u

    2010-09-01

    The emergence of tilted bipolar active regions (ARs) and the dispersal of their flux, mediated via processes such as diffusion, differential rotation, and meridional circulation, is believed to be responsible for the reversal of the Sun's polar field. This process (commonly known as the Babcock-Leighton mechanism) is usually modeled as a near-surface, spatially distributed {alpha}-effect in kinematic mean-field dynamo models. However, this formulation leads to a relationship between polar field strength and meridional flow speed which is opposite to that suggested by physical insight and predicted by surface flux-transport simulations. With this in mind, we present an improved double-ring algorithm for modeling the Babcock-Leighton mechanism based on AR eruption, within the framework of an axisymmetric dynamo model. Using surface flux-transport simulations, we first show that an axisymmetric formulation-which is usually invoked in kinematic dynamo models-can reasonably approximate the surface flux dynamics. Finally, we demonstrate that our treatment of the Babcock-Leighton mechanism through double-ring eruption leads to an inverse relationship between polar field strength and meridional flow speed as expected, reconciling the discrepancy between surface flux-transport simulations and kinematic dynamo models.

  7. GREET 1.0 -- Transportation fuel cycles model: Methodology and use

    SciTech Connect (OSTI)

    Wang, M.Q.

    1996-06-01

    This report documents the development and use of the Greenhouse Gases, Regulated Emissions, and Energy Use in Transportation (GREET) model. The model, developed in a spreadsheet format, estimates the full fuel-cycle emissions and energy use associated with various transportation fuels for light-duty vehicles. The model calculates fuel-cycle emissions of five criteria pollutants (volatile organic compounds, Co, NOx, SOx, and particulate matter measuring 10 microns or less) and three greenhouse gases (carbon dioxide, methane, and nitrous oxide). The model also calculates the total fuel-cycle energy consumption, fossil fuel consumption, and petroleum consumption using various transportation fuels. The GREET model includes 17 fuel cycles: petroleum to conventional gasoline, reformulated gasoline, clean diesel, liquefied petroleum gas, and electricity via residual oil; natural gas to compressed natural gas, liquefied petroleum gas, methanol, hydrogen, and electricity; coal to electricity; uranium to electricity; renewable energy (hydropower, solar energy, and wind) to electricity; corn, woody biomass, and herbaceous biomass to ethanol; and landfill gases to methanol. This report presents fuel-cycle energy use and emissions for a 2000 model-year car powered by each of the fuels that are produced from the primary energy sources considered in the study.

  8. Basic knowledge on radiative and transport properties to begin in thermal plasmas modelling

    SciTech Connect (OSTI)

    Cressault, Y.

    2015-05-15

    This paper has for objectives to present the radiative and the transport properties for people beginning in thermal plasmas. The first section will briefly recall the equations defined in numerical models applied to thermal plasmas; the second section will particularly deal with the estimation of radiative losses; the third part will quickly present the thermodynamics properties; and the last part will concern the transport coefficients (thermal conductivity, viscosity and electrical conductivity of the gas or mixtures of gases). We shall conclude the paper with a discussion about the validity of these results the lack of data for some specific applications, and some perspectives concerning these properties for non-equilibrium thermal plasmas.

  9. Conditional stochastic modeling of transport of contaminant in the vadose zone. Final project report

    SciTech Connect (OSTI)

    Yeh, T.C.J.; Harter, T.

    1995-06-01

    Spatial heterogeneity media leads to uncertainty in predicting both flow and transport in the vadose zone. In this work an efficient and flexible, combined analytical-numerical Monte Carlo approach is developed for the analysis of steady-state flow and transient transport processes in highly heterogeneous, variably saturated porous media. The approach is also used for the investigation of the validity of linear, first order analytical stochastic models. A combined analytical-numerical conditional simulation algorithm is developed to estimate the impact of in-situ soil hydraulic measurements on reducing the uncertainty of concentration and solute flux predictions.

  10. Link between hopping models and percolation scaling laws for charge transport in mixtures of small molecules

    DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)

    Ha, Dong -Gwang; Kim, Jang -Joo; Baldo, Marc A.

    2016-04-29

    Mixed host compositions that combine charge transport materials with luminescent dyes offer superior control over exciton formation and charge transport in organic light emitting devices (OLEDs). Two approaches are typically used to optimize the fraction of charge transport materials in a mixed host composition: either an empirical percolative model, or a hopping transport model. We show that these two commonly-employed models are linked by an analytic expression which relates the localization length to the percolation threshold and critical exponent. The relation is confirmed both numerically and experimentally through measurements of the relative conductivity of Tris(4-carbazoyl-9-ylphenyl) amine (TCTA) :1,3-bis(3,5-dipyrid-3-yl-phenyl) benzene (BmPyPb)more » mixtures with different concentrations, where the TCTA plays a role as hole conductor and the BmPyPb as hole insulator. Furthermore, the analytic relation may allow the rational design of mixed layers of small molecules for high-performance OLEDs.« less

  11. Surface chemical reactivity of ultrathin Pd(111) films on Ru(0001): Importance of orbital symmetry in the application of the d-band model

    SciTech Connect (OSTI)

    Yin, Xiangshi; Cooper, Valentino R.; Weitering, Hanno H.; Snijders, Paul C.

    2015-09-22

    The chemical bonding of adsorbate molecules on transition-metal surfaces is strongly influenced by the hybridization between the molecular orbitals and the metal d-band. The strength of this interaction is often correlated with the location of the metal d-band center relative to the Fermi level. Here, we exploit finite size effects in the electronic structure of ultrathin Pd(111) films grown on Ru(0001) to tune their reactivity by changing the film thickness one atom layer at a time, while keeping all other variables unchanged. Interestingly, while bulk Pd(111) is reactive toward oxygen, Pd(111) films below five monolayers are surprisingly inert. This observation is fully in line with the d-band model prediction when applied to the orbitals involved in the bonding. The shift of the d-band center with film thickness is primarily attributed to shifts in the partial density of states associated with the 4dxz and 4dyz orbitals. This study provides an in-depth look into the orbital specific contributions to the surface chemical reactivity, providing new insights that could be useful in surface catalysis.

  12. Surface chemical reactivity of ultrathin Pd(111) films on Ru(0001): Importance of orbital symmetry in the application of the d-band model

    DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)

    Yin, Xiangshi; Cooper, Valentino R.; Weitering, Hanno H.; Snijders, Paul C.

    2015-09-22

    The chemical bonding of adsorbate molecules on transition-metal surfaces is strongly influenced by the hybridization between the molecular orbitals and the metal d-band. The strength of this interaction is often correlated with the location of the metal d-band center relative to the Fermi level. Here, we exploit finite size effects in the electronic structure of ultrathin Pd(111) films grown on Ru(0001) to tune their reactivity by changing the film thickness one atom layer at a time, while keeping all other variables unchanged. Interestingly, while bulk Pd(111) is reactive toward oxygen, Pd(111) films below five monolayers are surprisingly inert. This observationmore » is fully in line with the d-band model prediction when applied to the orbitals involved in the bonding. The shift of the d-band center with film thickness is primarily attributed to shifts in the partial density of states associated with the 4dxz and 4dyz orbitals. This study provides an in-depth look into the orbital specific contributions to the surface chemical reactivity, providing new insights that could be useful in surface catalysis.« less

  13. The Role of Electron Transport and Trapping in MOS Total-Dose Modeling

    SciTech Connect (OSTI)

    Fleetwood, D.M.; Winokur, P.S.; Riewe, L.C.; Flament, O.; Paillet, P.; Leray, J.L.

    1999-07-19

    Radiation-induced hole and electron transport and trapping are fundamental to MOS total-dose models. Here we separate the effects of electron-hole annihilation and electron trapping on the neutralization of radiation-induced charge during switched-bias irradiation for hard and soft oxides, via combined thermally stimulated current (TSC) and capacitance-voltage measurements. We also show that present total-dose models cannot account for the thermal stability of deeply trapped electrons near the Si/SiO{sub 2} interface, or the inability of electrons in deep or shallow traps to contribute to TSC at positive bias following (1) room-temperature, (2) high-temperature, or (3) switched-bias irradiation. These results require revisions of modeling parameters and boundary conditions for hole and electron transport in SiO{sub 2}. The nature of deep and shallow electron traps in the near-interfacial SiO{sub 2} is discussed.

  14. Characterizing Indoor Airflow and Pollutant Transport using Simulation Modeling for Prototypical Buildings. I. Office Buildings

    SciTech Connect (OSTI)

    Sohn, M.D.; Daisey, J.M.; Feustel, H.E.

    1999-06-01

    This paper describes the first efforts at developing a set of prototypical buildings defined to capture the key features affecting airflow and pollutant transport in buildings. These buildings will be used to model airflow and pollutant transport for emergency response scenarios when limited site-specific information is available and immediate decisions must be made, and to better understand key features of buildings controlling occupant exposures to indoor pollutant sources. This paper presents an example of this approach for a prototypical intermediate-sized, open style, commercial building. Interzonal transport due to a short-term source release, e.g., accidental chemical spill, in the bottom and the upper floors is predicted and corresponding HVAC system operation effects and potential responses are considered. Three-hour average exposure estimates are used to compare effects of source location and HVAC operation.

  15. FINAL REPORT:Observation and Simulations of Transport of Molecules and Ions Across Model Membranes

    SciTech Connect (OSTI)

    MURAD, SOHAIL; JAMESON, CYNTHIA J

    2013-10-22

    During the this new grant we developed a robust methodology for investigating a wide range of properties of phospho-lipid bilayers. The approach developed is unique because despite using periodic boundary conditions, we can simulate an entire experiment or process in detail. For example, we can follow the entire permeation process in a lipid-membrane. This includes transport from the bulk aqueous phase to the lipid surface; permeation into the lipid; transport inside the lipid; and transport out of the lipid to the bulk aqueous phase again. We studied the transport of small gases in both the lipid itself and in model protein channels. In addition, we have examined the transport of nanocrystals through the lipid membrane, with the main goal of understanding the mechanical behavior of lipids under stress including water and ion leakage and lipid flip flop. Finally we have also examined in detail the deformation of lipids when under the influence of external fields, both mechanical and electrostatic (currently in progress). The important observations and conclusions from our studies are described in the main text of the report

  16. MODELING OF THE GROUNDWATER TRANSPORT AROUND A DEEP BOREHOLE NUCLEAR WASTE REPOSITORY

    SciTech Connect (OSTI)

    N. Lubchenko; M. Rodríguez-Buño; E.A. Bates; R. Podgorney; E. Baglietto; J. Buongiorno; M.J. Driscoll

    2015-04-01

    The concept of disposal of high-level nuclear waste in deep boreholes drilled into crystalline bedrock is gaining renewed interest and consideration as a viable mined repository alternative. A large amount of work on conceptual borehole design and preliminary performance assessment has been performed by researchers at MIT, Sandia National Laboratories, SKB (Sweden), and others. Much of this work relied on analytical derivations or, in a few cases, on weakly coupled models of heat, water, and radionuclide transport in the rock. Detailed numerical models are necessary to account for the large heterogeneity of properties (e.g., permeability and salinity vs. depth, diffusion coefficients, etc.) that would be observed at potential borehole disposal sites. A derivation of the FALCON code (Fracturing And Liquid CONvection) was used for the thermal-hydrologic modeling. This code solves the transport equations in porous media in a fully coupled way. The application leverages the flexibility and strengths of the MOOSE framework, developed by Idaho National Laboratory. The current version simulates heat, fluid, and chemical species transport in a fully coupled way allowing the rigorous evaluation of candidate repository site performance. This paper mostly focuses on the modeling of a deep borehole repository under realistic conditions, including modeling of a finite array of boreholes surrounded by undisturbed rock. The decay heat generated by the canisters diffuses into the host rock. Water heating can potentially lead to convection on the scale of thousands of years after the emplacement of the fuel. This convection is tightly coupled to the transport of the dissolved salt, which can suppress convection and reduce the release of the radioactive materials to the aquifer. The purpose of this work has been to evaluate the importance of the borehole array spacing and find the conditions under which convective transport can be ruled out as a radionuclide transport mechanism

  17. PLUVIUS: a generalized one-dimensional model of reactive pollutant behavior, including dry deposition, precipitation formation, and wet removal. Second edition

    SciTech Connect (OSTI)

    Easter, R.C.; Hales, J.M.

    1984-11-01

    This report is a second-edition user's manual for the PLUVIUS reactive-storm model. The PLUVIUS code simulates the formation of storm systems of a variety of types, and characterizes the behavior of air pollutants as they flow through, react within, and are scavenged by the storms. The computer code supplied with this report is known as PLUVIUS MOD 5.0, and is a substantial improvement over the MOD 3.1 version given in the original user's manual. Example applications of MOD 5.0 are given in the report to facilitate rapid application of the code for a variety of specific uses. 22 references, 7 figures, 48 tables.

  18. A charge carrier transport model for donor-acceptor blend layers

    SciTech Connect (OSTI)

    Fischer, Janine Widmer, Johannes; Koerner, Christian; Vandewal, Koen; Leo, Karl; Kleemann, Hans; Tress, Wolfgang; Riede, Moritz

    2015-01-28

    Highly efficient organic solar cells typically comprise donor-acceptor blend layers facilitating effective splitting of excitons. However, the charge carrier mobility in the blends can be substantially smaller than in neat materials, hampering the device performance. Currently, available mobility models do not describe the transport in blend layers entirely. Here, we investigate hole transport in a model blend system consisting of the small molecule donor zinc phthalocyanine (ZnPc) and the acceptor fullerene C{sub 60} in different mixing ratios. The blend layer is sandwiched between p-doped organic injection layers, which prevent minority charge carrier injection and enable exploiting diffusion currents for the characterization of exponential tail states from a thickness variation of the blend layer using numerical drift-diffusion simulations. Trap-assisted recombination must be considered to correctly model the conductivity behavior of the devices, which are influenced by local electron currents in the active layer, even though the active layer is sandwiched in between p-doped contacts. We find that the density of deep tail states is largest in the devices with 1:1 mixing ratio (E{sub t} = 0.14 eV, N{sub t} = 1.2 × 10{sup 18 }cm{sup −3}) directing towards lattice disorder as the transport limiting process. A combined field and charge carrier density dependent mobility model are developed for this blend layer.

  19. COMIS -- an international multizone air-flow and contaminant transport model

    SciTech Connect (OSTI)

    Feustel, H.E.

    1998-08-01

    A number of interzonal models have been developed to calculate air flows and pollutant transport mechanisms in both single and multizone buildings. A recent development in multizone air-flow modeling, the COMIS model, has a number of capabilities that go beyond previous models, much as COMIS can be used as either a stand-alone air-flow model with input and output features or as an infiltration module for thermal building simulation programs. COMIS was designed during a 12 month workshop at Lawrence Berkeley National Laboratory (LBNL) in 1988-89. In 1990, the Executive Committee of the International Energy Agency`s Energy Conservation in Buildings and Community Systems program created a working group on multizone air-flow modeling, which continued work on COMIS. The group`s objectives were to study physical phenomena causing air flow and pollutant (e.g., moisture) transport in multizone buildings, develop numerical modules to be integrated in the previously designed multizone air flow modeling system, and evaluate the computer code. The working group supported by nine nations, officially finished in late 1997 with the release of IISiBat/COMIS 3.0, which contains the documented simulation program COMIS, the user interface IISiBat, and reports describing the evaluation exercise.

  20. cDF Theory Software for mesoscopic modeling of equilibrium and transport phenomena

    SciTech Connect (OSTI)

    2015-12-01

    The approach is based on classical Density Functional Theory ((cDFT) coupled with the Poisson-Nernst-Planck (PNP) transport kinetics model and quantum mechanical description of short-range interaction and elementary transport processes. The model we proposed and implemented is fully atomistic, taking into account pairwise short-range and manybody long-range interactions. But in contrast to standard molecular dynamics (MD) simulations, where long-range manybody interactions are evaluated as a sum of pair-wise atom-atom contributions, we include them analytically based on wellestablished theories of electrostatic and excluded volume interactions in multicomponent systems. This feature of the PNP/cDFT approach allows us to reach well beyond the length-scales accessible to MD simulations, while retaining the essential physics of interatomic interactions from first principles and in a parameter-free fashion.

  1. Flow and Containment Transport Code for Modeling Variably Saturated Porous Media

    Energy Science and Technology Software Center (OSTI)

    1998-05-14

    FACT is a finite element based code designed to model subsurface flow and contaminant transport. It was designed to perform transient three-dimensional calculations that simulate isothermal groundwater flow, moisture movement, and solute transport in variably saturated and fully saturated subsurface porous media. The code is designed specifically to handle complex multi-layer and/or heterogenous aquifer systems in an efficient manner and accommodates a wide range of boundary conditions. Additionally 1-D and 2-D (in Cartesian coordinates) problemsmore » are handled in FACT by simply limiting the number of elements in a particular direction(s) to one. The governing equations in FACT are formulated only in Cartesian coordinates. FACT writes out both ascii and graphical binary files that are TECPLOT-ready. Special features are also available within FACT for handling the typical groundwater modeling needs for remediation efforts at the Savannah River Site.« less

  2. cDF Theory Software for mesoscopic modeling of equilibrium and transport phenomena

    Energy Science and Technology Software Center (OSTI)

    2015-12-01

    The approach is based on classical Density Functional Theory ((cDFT) coupled with the Poisson-Nernst-Planck (PNP) transport kinetics model and quantum mechanical description of short-range interaction and elementary transport processes. The model we proposed and implemented is fully atomistic, taking into account pairwise short-range and manybody long-range interactions. But in contrast to standard molecular dynamics (MD) simulations, where long-range manybody interactions are evaluated as a sum of pair-wise atom-atom contributions, we include them analytically based onmore » wellestablished theories of electrostatic and excluded volume interactions in multicomponent systems. This feature of the PNP/cDFT approach allows us to reach well beyond the length-scales accessible to MD simulations, while retaining the essential physics of interatomic interactions from first principles and in a parameter-free fashion.« less

  3. Flow and Containment Transport Code for Modeling Variably Saturated Porous Media

    Energy Science and Technology Software Center (OSTI)

    1998-05-14

    FACT is a finite element based code designed to model subsurface flow and contaminant transport. It was designed to perform transient three-dimensional calculations that simulate isothermal groundwater flow, moisture movement, and solute transport in variably saturated and fully saturated subsurface porous media. The code is designed specifically to handle complex multi-layer and/or heterogenous aquifer systems in an efficient manner and accommodates a wide range of boundary conditions. Additionally 1-D and 2-D (in Cartesian coordinates) problemsmore »are handled in FACT by simply limiting the number of elements in a particular direction(s) to one. The governing equations in FACT are formulated only in Cartesian coordinates. FACT writes out both ascii and graphical binary files that are TECPLOT-ready. Special features are also available within FACT for handling the typical groundwater modeling needs for remediation efforts at the Savannah River Site.« less

  4. A multi-model assessment of pollution transport to the Arctic

    SciTech Connect (OSTI)

    Shindell, D T; Chin, M; Dentener, F; Doherty, R M; Faluvegi, G; Fiore, A M; Hess, P; Koch, D M; MacKenzie, I A; Sanderson, M G; Schultz, M G; Schulz, M; Stevenson, D S; Teich, H; Textor, C; Wild, O; Bergmann, D J; Bey, I; Bian, H; Cuvelier, C; Duncan, B N; Folberth, G; Horowitz, L W; Jonson, J; Kaminski, J W; Marmer, E; Park, R; Pringle, K J; Schroeder, S; Szopa, S; Takemura, T; Zeng, G; Keating, T J; Zuber, A

    2008-03-13

    We examine the response of Arctic gas and aerosol concentrations to perturbations in pollutant emissions from Europe, East and South Asia, and North America using results from a coordinated model intercomparison. These sensitivities to regional emissions (mixing ratio change per unit emission) vary widely across models and species. Intermodel differences are systematic, however, so that the relative importance of different regions is robust. North America contributes the most to Arctic ozone pollution. For aerosols and CO, European emissions dominate at the Arctic surface but East Asian emissions become progressively more important with altitude, and are dominant in the upper troposphere. Sensitivities show strong seasonality: surface sensitivities typically maximize during boreal winter for European and during spring for East Asian and North American emissions. Mid-tropospheric sensitivities, however, nearly always maximize during spring or summer for all regions. Deposition of black carbon (BC) onto Greenland is most sensitive to North American emissions. North America and Europe each contribute {approx}40% of total BC deposition to Greenland, with {approx}20% from East Asia. Elsewhere in the Arctic, both sensitivity and total BC deposition are dominated by European emissions. Model diversity for aerosols is especially large, resulting primarily from differences in aerosol physical and chemical processing (including removal). Comparison of modeled aerosol concentrations with observations indicates problems in the models, and perhaps, interpretation of the measurements. For gas phase pollutants such as CO and O{sub 3}, which are relatively well-simulated, the processes contributing most to uncertainties depend on the source region and altitude examined. Uncertainties in the Arctic surface CO response to emissions perturbations are dominated by emissions for East Asian sources, while uncertainties in transport, emissions, and oxidation are comparable for European

  5. Mathematical model for oil slick transport and mixing in rivers. Special report

    SciTech Connect (OSTI)

    Shen, H.T.; Yapa, P.D.; Wang, D.S.; Yang, X.Q.

    1993-08-01

    The growing concern over the impacts of oil spills on aquatic environments has led to the development of many computer models for simulating the transport and spreading of oil slicks in surface waters. Almost all of these models were developed for coastal environments. A few river models exist. These models only considered the movement of surface oil slicks. In this study a two-layer model, ROSS2, is developed for simulating oil spills in rivers. This model considers the oil in the river to consist of a surface slick and suspended oil droplets entrained over the depth of the flow. The oil transformation processes considered in the model include advection, mechanical spreading, turbulent diffusion and mixing, evaporation, dissolution, emulsification, shoreline deposition and sinking. The model can be used for simulating instantaneous or continuous spills either on or under the water surface in rivers with or without an ice cover. The model has been implemented for the Ohio-Monongahela-Allegheny river system and the upper St. Lawrence River. This report describes the model formulation and implementation. A case study is presented along with detailed explanations of the program structure and its input and output. Although it is developed for simulating oil spills, the model can be applied to spills of other hazardous materials. Computer models, Oil spills, Oil slicks, Rivers.

  6. Modeling Gas-Phase Transport in Polymer-Electrolyte FuelCells

    SciTech Connect (OSTI)

    Weber, A.Z.; Newman, J.

    2006-08-17

    In this transaction, the equations and methodology for modeling convection and ordinary, Knudsen, and pressure diffusion of gases in a fuel-cell gas-diffusion layer are described. Some results examining the magnitudes of the various terms are also made. This derivation results in a self-consistent description of the various transport mechanisms and is robust for numerical solutions, especially for conditions involving different flow regimes or where the regime is not known a priori.

  7. A model for thermally driven heat and air transport in passive solar buildings

    SciTech Connect (OSTI)

    Jones, G.F.; Balcomb, J.D.; Otis, D.R.

    1985-01-01

    A model for transient interzone heat and air flow transport in passive solar buildings is presented incorporating wall boundary layers in stratified zones, and with interzone transport via apertures (doors and windows). The model includes features that have been observed in measurements taken in more than a dozen passive solar buildings. The model includes integral formulations of the laminar and turbulent boundary layer equations for the vertical walls which are then coupled to a one-dimensional core model for each zone. The cores in each zone exchange mass and energy through apertures that are modeled by an orifice type equation. The procedure is transient in that time dependence is retained only in the core equations which are solved by an explicit method. The model predicts room stratification of about 2/sup 0/C/m (1.1/sup 0/F/ft) for a room-to-room temperature difference of 0.56/sup 0/C(1/sup 0/F) which is in general agreement with the data.

  8. Model for thermally driven heat and air transport in passive solar buildings

    SciTech Connect (OSTI)

    Jones, G.F.; Balcomb, J.D.; Otis, D.R.

    1985-01-01

    A model for transient interzone heat and air flow transport in passive solar buildings is presented incorporating wall boundary layers in stratified zones, and with interzone transport via apertures (doors and windows). The model includes features that have been observed in measurements taken in more than a dozen passive solar buildings. The model includes integral formulations of the laminar and turbulent boundary layer equations for the vertical walls which are then coupled to a one-dimensional core model for each zone. The cores in each zone exchange mass and energy through apertures that are modeled by an orifice type equation. The procedure is transient in that time dependence is retained only in the core equations which are solved by an explicit method. The model predicts room stratification of about 2/sup 0/C/m (1.1/sup 0/F/ft) for a room-to-room temperature difference of 0.56/sup 0/C(1/sup 0/F) which is in general agreement with the data. 38 references, 10 figures, 1 table.

  9. Mark 22 Reactivity

    SciTech Connect (OSTI)

    Buckner, M.R.

    2001-07-02

    Calculations for reactivity held in control rods have underpredicted the observed Mark 22 reactivity. Reactivity predictions by charge designers have accounted for this by including large biases which change with exposure and reactor region. The purpose of this study was to thoroughly investigate the methods and data used in the reactivity calculations. The goal was to identify errors and improvements and make necessary corrections.

  10. Development and applications of GREET 2.7 -- The Transportation Vehicle-CycleModel.

    SciTech Connect (OSTI)

    Burnham, A.; Wang, M. Q.; Wu, Y.

    2006-12-20

    Argonne National Laboratory has developed a vehicle-cycle module for the Greenhouse gases, Regulated Emissions, and Energy use in Transportation (GREET) model. The fuel-cycle GREET model has been cited extensively and contains data on fuel cycles and vehicle operations. The vehicle-cycle model evaluates the energy and emission effects associated with vehicle material recovery and production, vehicle component fabrication, vehicle assembly, and vehicle disposal/recycling. With the addition of the vehicle-cycle module, the GREET model now provides a comprehensive, lifecycle-based approach to compare the energy use and emissions of conventional and advanced vehicle technologies (e.g., hybrid electric vehicles and fuel cell vehicles). This report details the development and application of the GREET 2.7 model. The current model includes six vehicles--a conventional material and a lightweight material version of a mid-size passenger car with the following powertrain systems: internal combustion engine, internal combustion engine with hybrid configuration, and fuel cell with hybrid configuration. The model calculates the energy use and emissions that are required for vehicle component production; battery production; fluid production and use; and vehicle assembly, disposal, and recycling. This report also presents vehicle-cycle modeling results. In order to put these results in a broad perspective, the fuel-cycle model (GREET 1.7) was used in conjunction with the vehicle-cycle model (GREET 2.7) to estimate total energy-cycle results.