Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffusion processes”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Adsorption Studies of Dimetridazole and Metronidazole onto Biochar Derived from Sugarcane Bagasse: Kinetic, Equilibrium, and Mechanisms

Biochar derived from sugarcane bagasse was prepared at different pyrolysis temperatures and holding times for adsorption of dimetridazole and metronidazole from aqueous solution. The optimal pyrolysis temperature and holding time were 500 °C and 120 min, respectively. The removal efficiencies of dimetridazole and metronidazole were 98.0% and 88.8%, respectively, at biochar dosage of 4.0 g L{sup −1} and solution temperature of 30 °C. Kinetics data were fitted by four kinetics models, and results indicated that adsorption was governed by pseudo-second-order kinetics. Adsorption is a multi-step process that involves film diffusion and pore filling. In addition, five isotherm models were employed to describe the adsorption equilibrium. Brunauer–Emmett–Teller model presented better fitting for the adsorption equilibrium, and the maximum adsorption capacities were 72.17 and 23.61 mg g{sup − 1} for dimetridazole and metronidazole, respectively. The calculated values of ΔG{sup 0} and ΔH{sup 0} indicated the non-spontaneous and exothermic nature of the adsorption process at the range of temperature studied. Thermodynamic studies also revealed that physical and chemical adsorptions were co-action. The adsorption mechanisms of dimetridazole and metronidazole onto biochar were mainly hydrogen bonding and π–π interaction. All the results revealed that sugarcane bagasse biochar can be used as alternative to costly adsorbents for the removal of dimetridazole and metronidazole from aqueous solution.

36 MATERIALS SCIENCE↗
DOE OSTI2024

2D-EFICACY: Control of Metastable 2D Carbide[1]Chalcogenide Heterolayers: Strain and Moire Engineering

The experimental isolation of graphene led to the discovery of an entirely new world of two-dimensional (2D) materials in which the 2D nature often leads to emergent behaviors not seen in bulk systems. 2D transition metal dichalcogenides (TMDs) exhibit physico-chemical properties that depend on the transition metal, polymorph, thickness, and presence and type of defects. Recently, a group of thin (10-100nm) transition metal carbides (TMCs), such as Mo2C, has been synthesized that exhibit a thickness-dependent superconducting critical temperature (Tc). These thin TMCs are different from MXenes, another class of 2D materials consisting of few layers of nitrides or carbides (<5nm) produced by chemical etching and delamination. The goal of this renewal proposal is to combine experiment and computation to synthesize and elucidate the guiding principles that control the growth, orientation and strain of heterostacks of thin TMCs and TMDs composed with Nb, Ti and W. We expect to stabilize metastable hybrid phases of TMCs sandwiched between TMDs (H-TMD/Cs) with unprecedented physico-chemical properties. As part of previous DOE-funded work by the Terrones/Sinnott groups, thin (10-100 nm thick) Mo2C flakes were successfully synthesized by chemical vapor deposition (CVD). By subsequently exposing Mo2C to H2S, partial chalcogenization was achieved, resulting in heterostacks of MoCx phases and MoS2. The formation of MoS2 led to a deficiency of Mo atoms in the underlying Mo2C, resulting in an inhomogeneous phase change from α-Mo2C to γ’-MoCx and then to γ-MoC. The γ’-MoCx is a strained metastable phase and the heterostack of all three phases demonstrated an increased Tc relative to that of α-Mo2C, from 4 to 6K; its interleaved layered structure consisting of superconducting and semiconducting phases is ideal for future studies of Josephson junction series arrays. Moiré patterns in these heterostacked systems could result in new phenomena, as moiré patterns in bilayer graphene showed unconventional superconductivity and moiré excitons have been observed in twisted TMD heterobilayers. The scientific hypothesis of the proposed synergistic computational and experimental research is that orientation and strain control within confined thin metastable TMCs, sandwiched by stable phases of TMCs and layered TMDs, will depend on kinetic and thermodynamic “knobs” that include fast temperature changes, chalcogen diffusion through preferred crystallographic planes, reaction times, pressure, reactive atmosphere, precursors, and surfactants, which will also tailor properties such as superconductivity, magnetism, ferroelectricity, piezoelectricity, and catalytic performance. We will develop the guiding principles for the synthesis and stabilization of metastable H-TMD/Cs based on Nb, Ti and W. In order to validate the hypothesis, four tasks are proposed: The first task will synthesize ultra-thin TMCs based on Nb, W and Ti, by: 1) adapting the CVD method used for Mo2C, 2) plasma assisted CVD, 3) defect-mediated CVD processes, and 4) cryo-milling of carbide powders. The second task will accomplish the synthesis and basic physico-chemical characterizations of H-TMD/Cs by chalcogenization of the materials synthesized in task one, and by carbonization of TMDs. H-TMD/Cs will also be investigated for their suitability in energy conversion applications such as supercapacitors, Li and multivalent ion batteries, and electrocatalysts, topics of interest to DoE. These tasks will be carried out in close conjunction with density functional theory (DFT) calculations with insights into energetics, lattice parameters, stability, phase diagrams, band structures, and density of states of H-TMD/Cs. The third task will characterize and evaluate strain and moiré patterns at the interfaces of different H-TMD/Cs by high-resolution scanning transmission electron microscopy (HR-STEM), scanning tunneling microscopy (STM), and conductive tip atomic force microscopy. Nudged elastic band calculations with DFT will be performed to understand the chalcogen diffusion process, which will provide insights into the interfaces between different phases of TMCs and TMDs. The fourth task aims at quantifying the stability and dynamics of H-TMD/Cs by in-situ TEM and Raman studies under heating, strain, and electrical biasing. Phonon calculations using DFT will provide a basis for interpreting Raman spectra. This coherent framework involving synthesis, characterization, and computation will result in a broad scientific impact for energy related applications. The ability to develop new H-TMD/Cs will enhance a range of applications that include batteries, catalysts, switches, sensors, quantum computing components and smart coatings.

2-Dimensional materials↗
DOE OSTI2024

2D-EFICACY: Control of Metastable 2D Carbide Chalcogenide Heterolayers: Strain and Moire Engineering

The experimental isolation of graphene led to the discovery of an entirely new world of two-dimensional (2D) materials in which the 2D nature often leads to emergent behaviors not seen in bulk systems. 2D transition metal dichalcogenides (TMDs) exhibit physico-chemical properties that depend on the transition metal, polymorph, thickness, and presence and type of defects. Recently, a group of thin (10-100nm) transition metal carbides (TMCs), such as Mo 2 C, has been synthesized that exhibit a thickness-dependent superconducting critical temperature (Tc). These thin TMCs are different from MXenes, another class of 2D materials consisting of few layers of nitrides or carbides (<5nm) produced by chemical etching and delamination. The goal of this renewal proposal is to combine experiment and computation to synthesize and elucidate the guiding principles that control the growth, orientation and strain of heterostacks of thin TMCs and TMDs composed with Nb, Ti and W. We expect to stabilize metastable hybrid phases of TMCs sandwiched between TMDs (H-TMD/Cs) with unprecedented physico-chemical properties. As part of previous DOE-funded work by the Terrones/Sinnott groups, thin (10-100 nm thick) Mo 2 C flakes were successfully synthesized by chemical vapor deposition (CVD). By subsequently exposing Mo 2 C to H 2 S, partial chalcogenization was achieved, resulting in heterostacks of MoCx phases and MoS 2 . The formation of MoS 2 led to a deficiency of Mo atoms in the underlying Mo 2 C, resulting in an inhomogeneous phase change from α-Mo 2 C to γ’-MoCx and then to γ-MoC. The γ’-MoCx is a strained metastable phase and the heterostack of all three phases demonstrated an increased Tc relative to that of α-Mo 2 C, from 4 to 6K; its interleaved layered structure consisting of superconducting and semiconducting phases is ideal for future studies of Josephson junction series arrays. Moiré patterns in these heterostacked systems could result in new phenomena, as moiré patterns in bilayer graphene showed unconventional superconductivity and moiré excitons have been observed in twisted TMD heterobilayers. The scientific hypothesis of the proposed synergistic computational and experimental research is that orientation and strain control within confined thin metastable TMCs, sandwiched by stable phases of TMCs and layered TMDs, will depend on kinetic and thermodynamic “knobs” that include fast temperature changes, chalcogen diffusion through preferred crystallographic planes, reaction times, pressure, reactive atmosphere, precursors, and surfactants, which will also tailor properties such as superconductivity, magnetism, ferroelectricity, piezoelectricity, and catalytic performance. We will develop the guiding principles for the synthesis and stabilization of metastable H-TMD/Cs based on Nb, Ti and W. In order to validate the hypothesis, four tasks are proposed: The first task will synthesize ultra-thin TMCs based on Nb, W and Ti, by: 1) adapting the CVD method used for Mo 2 C, 2) plasma assisted CVD, 3) defect-mediated CVD processes, and 4) cryo-milling of carbide powders. The second task will accomplish the synthesis and basic physico-chemical characterizations of H-TMD/Cs by chalcogenization of the materials synthesized in task one, and by carbonization of TMDs. H-TMD/Cs will also be investigated for their suitability in energy conversion applications such as supercapacitors, Li and multivalent ion batteries, and electrocatalysts, topics of interest to DoE. These tasks will be carried out in close conjunction with density functional theory (DFT) calculations with insights into energetics, lattice parameters, stability, phase diagrams, band structures, and density of states of H-TMD/Cs. The third task will characterize and evaluate strain and moiré patterns at the interfaces of different H-TMD/Cs by high-resolution scanning transmission electron microscopy (HR-STEM), scanning tunneling microscopy (STM), and conductive tip atomic force microscopy. Nudged elastic band calculations with DFT will be performed to understand the chalcogen diffusion process, which will provide insights into the interfaces between different phases of TMCs and TMDs. The fourth task aims at quantifying the stability and dynamics of H-TMD/Cs by in-situ TEM and Raman studies under heating, strain, and electrical biasing. Phonon calculations using DFT will provide a basis for interpreting Raman spectra. This coherent framework involving synthesis, characterization, and computation will result in a broad scientific impact for energy related applications. The ability to develop new H-TMD/Cs will enhance a range of applications that include batteries, catalysts, switches, sensors, quantum computing components and smart coatings.

2-Dimensional Materials↗
NASA NTRS2021

Groundwater Recharge Estimated by Land Surface Models: An Evaluation in the Conterminous U.S.

Estimating diffuse recharge of precipitation is fundamental to assessing groundwater sustainability. Diffuse recharge is also the process through which climate and climate change directly affect groundwater. In this study, we evaluated diffuse recharge over the conterminous U.S. simulated by a suite of land surface models (LSMs) that were forced using a common set of meteorological input data. Simulated annual recharge exhibited spatial patterns that were similar among the LSMs, with the highest values in the eastern U.S. and Pacific Northwest. However, the magnitudes of annual recharge varied significantly among the models and were associated with differences in simulated ET, runoff and snow. Evaluation against two independent datasets did not answer the question of whether the ensemble mean performs the best, due to inconsistency between those datasets. The amplitude and timing of seasonal maximum recharge differed among the models, influenced strongly by model physics governing deep soil moisture drainage rates and, in cold regions, snowmelt. Evaluation using in situ soil moisture observations suggested that true recharge peaks 1-3 months later than simulated recharge, indicating systematic biases in simulating deep soil moisture. However, recharge from lateral flows and through preferential flows cannot be inferred from soil moisture data, and the seasonal cycle of simulated groundwater storage actually compared well with in situ groundwater observations. Long-term trends in recharge were not consistently correlated with either precipitation trends or temperature trends. This study highlights the need to employ dynamic flow models in LSMs, among other improvements, to enable more accurate simulation of recharge.

snowmelt↗
DOE OSTI2023

Review—Fundamental Uranium Electrochemistry and Spectroscopy in Molten Salt Systems

Uranium is a key element used for nuclear energy production. Some advanced reactor designs, specifically molten salt reactors, will continue to use uranium as the fissile material for energy production. These new technologies require an intimate understanding of uranium chemistry during and after energy production. This review covers contemporary research on the coordination chemistry and behavior of uranium with the coolant and pyroprocessing salts as proposed for use in future reactor designs. Discussed topics include the nature of U redox reactions involving the reduction of U(III) to U metal and oxidation of U(III) to U(IV). These systems have been interrogated using cyclic voltammetry, chronopotentiometry, and optical and X-ray absorption spectroscopies. Insights obtained into the electrode potentials, the uranium species, and their diffusion coefficients in alkali halide melts from decades of research are summarized selectively. Further, perspectives are provided on the importance of unifying studies for comparison across multiple institutions. The application of synchrotron radiation research and multimodal approaches involving two (or more) probes, such as the widespread combination of UV–visible spectroscopy and electroanalysis known as spectroelectrochemistry, can provide new knowledge about the main process of uranium electrorefining—diffusion, as will be demonstrated in this review through the lack of comparable results.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2001

Combustion of Metals in Reduced-Gravity and Extraterrestrial Environments

As a result of the ongoing exploration of Mars and the several unmanned and possibly manned missions planned for the near future, increased attention has been given to the use of the natural resources of the planet for rocket propellant production and energy generation. Since the atmosphere of Mars consists of approximately 95% carbon dioxide (CO2), this gas is the resource of choice to be employed for these purposes. Since many metals burn vigorously with CO2, these may be used as an energy source or as propellants for a research vehicle on the surface of Mars. Shafirovich and Goldshleger conducted experiments with spherical particles up to 2.5 mm in diameter and found that the burning process was controlled by diffusion and that the particles exhibited pulsating combustion due to superheating of the Mg vapor trapped inside a protective oxide shell. They also proposed a reaction mechanism based on the gas-phase reaction, Mg + CO2 yields MgO + CO and the heterogeneous reaction Mg + CO yields MgO + C occurring on the sample surface. In all the above studies with large Mg particles, the burning process is invariably influenced by strong convective currents that accelerate the combustion reaction and shorten the burning times. Although these currents are nearly absent in the burning of small particles, the high emissivity of the flames, rapid reaction, and small length scales make the gathering of any useful information on burning rates and flame structure very difficult. The goal of this investigation is to provide a detailed study of flame structure by taking advantage of large, free-floating spherical metal samples and their corresponding long burning times available in a weightless environment. The use of reduced gravity is essential to eliminate the intrusive buoyant flows that plague high temperature metal reactions, to remove the destructive effect of gravity on the shape of molten metal samples, and to study the combustion behavior of metals in the presence of solid oxides undisturbed by natural convection. This work presents the most complete modeling of metal particle burning to date for Mg with CO2 and O2.

Branch, M. C.↗
NASA NTRS1993

Immiscible phase incorporation during directional solidification of hypermonotectics

Solidification processes in immiscible samples were investigated by directly observing the events taking place at the solid-liquid interface during directional solidification. Visualization of these events was made possible through the use of a transparent metal analog system and a temperature gradient stage assembly fitted to an optical microscope. The immiscible transparent analog system utilized was the succinonitrile-glycerol system. This system has been shown to exhibit the same morphological transitions as observed in metallic alloys of monotectic composition. Both monotectic and hypermonotectic composition samples were directionally solidified in order to gain an improved understanding of the manner in which the excess hypermonotectic liquid is incorporated into the solidifying structure. The processing conditions utilized prevented sedimentation of the excess hypermonotectic liquid by directionally solidifying the samples in very thin (13 microns), horizontally oriented cells. High thermal gradient to growth rate ratios (G/R) were used in an effort to prevent constitutional supercooling and the subsequent formation of L(sub 2) droplets in advance of the solidification front during the growth of fibrous composite structures. Results demonstrated that hypermonotectic composites could be produced in samples up to two weight percent off of the monotectic composition by using a G/R ratio greater than or equal to 4.6 x 10(exp 4) C(s)/mm(sup 2) to avoid constitutional supercooling. For hypermonotectic samples processed with G/R ratios below 4.6 x 10(exp 4) C(s)/mm(sup 2), constitutional supercooling occurred and resulted in slight interfacial instability. For these samples, two methods of incorporation of the hypermonotectic liquid were observed and are reported. The correlation between the phase spacing, lambda, and the growth rate, R, was examined and was found to obey a relationship generally associated with a diffusion controlled coupled growth process. For samples with compositions ranging from the monotectic composition up to 2 percent off of the monotectic composition, data indicated that the square of the phase spacing (lambda) varied linearly with the inverse of the growth rate (R).

Andrews, J. Barry↗
DOE OSTI2018

The Kinetics of the Liquid Phase Nucleation in a Stretched FCC Crystal: A Molecular Dynamics Simulation

The kinetics of a spontaneous formation of liquid phase in a stretched (superheated) Lennard-Jones fcc crystal is studied. Molecular dynamics experiments are used to determine the main parameters of the nucleation process: nucleation frequency J, diffusion coefficient of nuclei D{sub *}, nonequilibrium Zel’dovich factor Z, and critical nucleus size R{sub *}. The calculations are performed at negative pressures from the endpoint of the melting line and at positive pressures that are higher by a factor of eight than the critical pressure. The simulation results are compared to the classical homogeneous nucleation theory. It is found that the theory qualitatively correctly reproduces the dynamics of developing the process. The theory and the simulation demonstrate good quantitative agreement for the transition rate of the liquid phase nucleus through the critical size, but there is large difference in the numbers of critical nuclei in the unit volume of the metastable phase. In the case of significant superheatings and negative pressures, the contribution of the energy of elastic stresses to the moving force of the phase transformation is small and it can be neglected in a first approximation. The mismatch between the theory and the simulation results can be eliminated taking that the surface free energy of a curved “crystal–liquid droplet” interface is smaller than that of a plane interface by 30–35%.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Anomalous Oxidative Diffusion in Titanium Pyrotechnic Powders

It has long been observed that oxidation processes in metals tend to follow a parabolic rate law associated with the growth of a surface oxide layer. Here we observe that for certain titanium powders, the expected parabolic law (∝t 1/2 ) is recovered, yet for others, the exponent differs significantly. One explanation for this non-parabolic, anomalous diffusion arises from fractal geometry. Theoretical considerations indicate that the time response of diffusion-limited processes in an object closely follow a power-law in time (t n ) with n=(E–D)/2, where E is the object's Euclidean dimension and D is its boundary's Hausdorff dimension. Non-integer, (fractal) values of D will result in n≠1/2. Finite element simulations of several canonical fractal objects were performed to verify the application of this theory; the results matched the theory well. Two different types of titanium powder were tested in isothermal thermogravimetric tests under dilute oxygen. Time-dependent mass uptake data were fit with power-law forms and the associated exponents were used to determine an equivalent fractal dimension. One Ti powder type has an implied surface dimension of ca. 2.3 to 2.5, suggesting fractal geometry may be operative. Finally, the other has a dimension near 2.0, indicating it behaves like traditional material.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Perpendicular diffusion in the transport of solar energetic particles from unconnected sources: The counter-streaming particle beams revisited

In some solar energetic particle (SEP) events, a counter-streaming particle beam with a deep depression of flux at ∼90{sup ∘} pitch angle during the beginning phase is observed. Two different interpretations exist within the community to explain this interesting phenomenon. One explanation invokes the hypothesis of an outer reflecting boundary or a magnetic mirror beyond the observer. The other one considers the effect of perpendicular diffusion on the transport process of SEPs in interplanetary space. In this work, we revisit the problem of counter-streaming particle beams observed in SEP events and discuss the possible mechanisms responsible for the formation of this phenomenon. We clarify some results in previous works.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Volume nanograting formation in laser-silica interaction as a result of the 1D plasma-resonance ionization instability

The initial stage of the small-scale ionization-induced instability developing inside the fused silica volume exposed to the femtosecond laser pulse is studied as a possible initial cause of the self-organized nanograting formation. We have calculated the spatial spectra of the instability with the electron-hole diffusion taken into account for the first time and have found that it results in the formation of some hybrid (diffusion-wave) 1D structure with the spatial period determined as the geometrical mean of the laser wavelength and characteristic diffusion length of the process considered. Near the threshold of the instability, this period occurs to be approximately equal to the laser half-wavelength in the silica, close to the one experimentally observed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Innovative compact braze-free accelerating cavity

We report that high energy physics experiments using particle accelerators as well as industrial and medical applications are continuously seeking more compact, robust and cheaper accelerating structures. As of today, stable operating gradients, exceeding 100 MV/m, have been demonstrated by the SLAC group in the X-Band (11.424 GHz). These experiments show that hard structures, fabricated without high-temperature processes, achieve a better high gradient performance in terms of accelerating gradients. Therefore, we present an innovative and compact type of accelerating cavity that avoids any high-temperature processes like brazing or diffusion bonding. All cells are joined together by means of specifically designed and proprietary screws which ensure good vacuum and RF contacts. Two three-cell standing-wave accelerating structures, designed to operate in the pi-mode at 11.424 GHz, have been successfully built and cold tested. In order to guarantee a vacuum envelope and mechanically robust assembly, we used the Electron Beam Welding (EBW) and the Tungsten Inert Gas (TIG) processes. Lastly, this work has been carried out in the framework of a funded project by the Vth Committee of the INFN for the Laboratori Nazionali di Frascati (LNF), within a large international collaboration between LNF, SLAC and KEK for the development of X-Band accelerating cavities using "hard structure" technology.

43 PARTICLE ACCELERATORS↗
NASA NTRS1971

Investigation of the feasibility of developing low permeability polymeric films

The feasibility of reducing the gas permeability rate of Mylar and Kapton films without drastically effecting their flexibility characteristics at cryogenic temperatures was considered. This feasibility was established using a concept of diffusion bonding two layers of metallized films together forming a film-metal-film sandwich laminate. The permeability of kapton film to gaseous helium was reduced from a nominal ten = to the minus 9 power cc-mm/sq cm sec. cm Hg to ten to the minus 13 power cc-mm/ sq cm - sec. cm Hg with some values as low as ten to the minus 15 power cc - mm/sq cm m-sec - cm Hg being obtained. Similar reductions occurred in the liquid hydrogen permeability at -252 C. In the course of the program the permeability, flexibility and bond strength of plain, metalized and diffusion bond film were determined at +25 C, -195 C and -252 C. The cryogenic flexibility of Kapton film was reduced slightly due to the metallization process but no additional loss in flexibility resulted from the diffusion bonding process.

Hoggatt, J. T.↗
DOE OSTI2026

Achieving high rate performance in hybrid pristine-recycled cathodes using model-informed electrode designs

Direct recycling lithium-ion battery cathodes, a process that retains the engineered oxide structures from end-of-life materials, presents a cost-effective and energy-efficient alternative to other battery recycling methods. However, while direct-recycled cathodes have demonstrated performance comparable to that of pristine materials at low cycling rates, their high-rate performance remains uncertain. Morphology changes in cathode particles, a main mode of degradation, directly impact rate performance by limiting surface kinetics and solid-phase diffusion. If direct recycling processes do not sufficiently restore pristine-like morphologies, the recycled materials may retain structural defects that hinder high-rate performance. The present work uses a physics-based pseudo-2D model to simulate hybrid electrodes with pristine and artificially “aged/recycled” NMC materials to investigate potential impacts of incorporating performance-limited aged cathode materials into cells. The study highlights how differences in transport and kinetic properties can influence rate capabilities in mixed electrodes — particularly in high-loading cells in high-demand applications. However, model results also reveal a possible mitigation strategy via dual-layer electrode architectures with lower-performing materials positioned near the current collector. Simulations of 4.0 mAh cm −2 cells cycled at 4C using a dual-layer architecture provided approximately 5%–30% more capacity in constant-current protocols compared to homogeneously blended electrode architectures with the same loadings and mixed-material compositions. These findings highlight the importance of strategic electrode design in minimizing potential performance losses and facilitating the integration of recycled materials into high-performance batteries, advancing sustainable and cost-effective battery manufacturing.

25 ENERGY STORAGE↗
DOE OSTI2024

DECOVALEX-2023: Task F1 Final Report

DECOVALEX-2023 Task F is a comparison of models and methods for post-closure performance assessment (PA) of a deep geologic repository for radioactive waste. The general aims of Task F are to build confidence in the models, methods, and software used for PA and to stimulate additional research and development in PA methodologies. The task objectives are to motivate development of PA modelling skills and capabilities, to examine the influence of model choices on calculated repository performance, and to compare the uncertainties introduced by model choices to other sources of uncertainty. Task F involves no actual experiment or site. It is a PA modelling exercise that requires the conceptual development of hypothetical repository designs and geologic settings. Because three of the teams were interested in salt and the rest of the teams were interested in crystalline rock, Task F was split into two branches: Task F1 for crystalline rock and Task F2 for salt. This report is for Task F1, crystalline rock. Teams from seven countries (Canada, Czech Republic, Germany, Korea, Sweden, Taiwan, and United States) participated in Task F1. The teams worked together to define the features, events, and processes of the reference case repository and established a set of performance measures. In addition, they defined a set of benchmark problems designed to test and compare modelling capabilities for fracture flow and transport at different scales. The repository design and benchmark problems are documented in a Task Specification that evolved over time as the group honed the specifications. The benchmark problems verified that each team can aptly model flow and transport in fractured media in 1-, 2-, and 3-dimensions. Two general approaches were used for the 3-dimensional benchmarks: discrete fracture network (DFN) and equivalent continuous porous medium (ECPM). DFN modelling involves explicit meshing of each fracture while ECPM modelling aims to capture the effective porosity and directional permeability of each cell in a space-filling mesh as affected by intersecting fractures. In some models, a combination of the two is used, i.e., DFN for large known fractures and ECPM for the rest of the domain. Transport is solved by using either the advection-dispersion equation or particle tracking. Although some variation is observed among model breakthrough curves in the benchmark problems, there is strong agreement in breakthrough behaviour up to at least the 75 th percentile for all benchmarks. At the 90 th percentile, breakthrough results show larger differences, suggesting several models retain substantially higher fractions of tracer in regions of slower moving water. In addition to the flow and transport benchmarks, several teams completed the source term benchmark, verifying capabilities for modelling radionuclide decay and ingrowth, waste package breach, instant release fractions, fuel matrix degradation rates, and radionuclide solubility limitations. The reference case is conceptualized as a generic spent fuel repository at a depth of 450 m in fractured crystalline rock. The repository has 50 parallel backfilled drifts, each with 50 deposition holes 6 m apart. Each deposition hole contains a 4-PWR waste package and bentonite buffer. The rock domain is 5 km in length, 2 km in width, and 1 km in depth. It has 6 deterministic fractured deformation zones and a multitude of stochastic fractures. Teams generally used the ECPM approach for the entire rock or a hybrid approach in which the deterministic fracture zones are modelled with a DFN and the rest of the rock is modelled by ECPM. Of the reference case problems specified, only the results of the initial reference case problem are compared in this report. The initial problem focuses on transport from the deposition holes to the surface, i.e., it neglects waste package performance. Tracers are released at all waste package locations at time zero and tracked for their releases to the near field and ground surface. The water fluxes calculated at the ground surface entry and exit regions of the domain are similar for all models except for two that have considerably lower fluxes. For tracer transport, large differences are observed among models in the magnitude of tracer transported. Much of the difference appears to be due to how the repository is implemented and hence the different degrees of repository simplification. Models that exclude the drifts, buffer, and backfill from the domain tend to show greater release of tracers and radionuclides from the repository. The initial study presented here indicates that major differences in modelling important processes within the repository (e.g., diffusion through buffer and backfill) can produce broadly different release and transport results, especially when those processes are excluded. Even for the models that included all specified features, events, and processes, the results show significant differences and demonstrate the importance of examining multiple modelling approaches in performance assessment. The differences in results observed in this study are expected to motivate teams to either increase complexity in future versions of the reference case models or to improve methods to account for the effects of simplified features and processes. Either way, future improvements in these models are expected to produce results that more closely agree.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS1992

Space Station Freedom combustion research

Extended operations in microgravity, on board spacecraft like Space Station Freedom, provide both unusual opportunities and unusual challenges for combustion science. On the one hand, eliminating the intrusion of buoyancy provides a valuable new perspective for fundamental studies of combustion phenomena. On the other hand, however, the absence of buoyancy creates new hazards of fires and explosions that must be understood to assure safe manned space activities. These considerations - and the relevance of combustion science to problems of pollutants, energy utilization, waste incineration, power and propulsion systems, and fire and explosion hazards, among others - provide strong motivation for microgravity combustion research. The intrusion of buoyancy is a greater impediment to fundamental combustion studies than to most other areas of science. Combustion intrinsically heats gases with the resulting buoyant motion at normal gravity either preventing or vastly complicating measurements. Perversely, this limitation is most evident for fundamental laboratory experiments; few practical combustion phenomena are significantly affected by buoyancy. Thus, we have never observed the most fundamental combustion phenomena - laminar premixed and diffusion flames, heterogeneous flames of particles and surfaces, low-speed turbulent flames, etc. - without substantial buoyant disturbances. This precludes rational merging of theory, where buoyancy is of little interest, and experiments, that always are contaminated by buoyancy, which is the traditional path for developing most areas of science. The current microgravity combustion program seeks to rectify this deficiency using both ground-based and space-based facilities, with experiments involving space-based facilities including: laminar premixed flames, soot processes in laminar jet diffusion flames, structure of laminar and turbulent jet diffusion flames, solid surface combustion, one-dimensional smoldering, ignition and flame spread of liquids, drop combustion, and quenching of panicle-air flames. Unfortunately, the same features that make microgravity attractive for fundamental combustion experiments, introduce new fire and explosion hazards that have no counterpart on earth. For example, microgravity can cause broader flammability limits, novel regimes of flame spread, enhanced effects of flame radiation, slower fire detector response, and enhanced combustion upon injecting fire extinguishing agents, among others. On the other hand, spacecraft provide an opportunity to use 'fire-safe' atmospheres due to their controlled environment. Investigation of these problems is just beginning, with specific fire safety experiments supplementing the space based fundamental experiments listed earlier; thus, much remains to be done to develop an adequate technology base for fire and explosion safety considerations for spacecraft.

Faeth, G. M.↗
DOE OSTI2020

A Comparison of RESRAD and GoldSim Models for Assessing Radiological Dose for a RCRA Landfill - 20418

US Ecology Idaho (USEI) operates a Resource Conservation and Recovery Act (RCRA) Subtitle C, Hazardous Waste Treatment, Storage, and Disposal (TSD) Facility on a 640-acre property in Owyhee County, Idaho. USEI accepts a wide variety of RCRA-exempt low activity radioactive wastes, including naturally occurring radioactive material (NORM) and technologically enhanced NORM (TENORM) (USEI 2009). The performance of the landfill with regard to radiological operating permit requirements was previously assessed using a model constructed with the residual radioactivity (RESRAD) computer program, developed by Argonne National Laboratory for the U.S. Department of Energy (DOE). RESRAD supports the evaluation of several environmental transport pathways related to a radionuclide-contaminated soil source term, but it was not conceived as a model for the evaluation of landfill radiological performance. A more comprehensive and realistic model is desired in order to better support submissions to both the State of Idaho's Department of Environmental Quality as well as the U.S. Nuclear Regulatory Commission. A Performance Assessment (PA) computer model was developed using GoldSim software for the USEI RCRA Subtitle C landfill. The RESRAD computer model was specifically developed for calculating soil cleanup criteria and radiological dose and cancer risk from residual radioactive material in soil. To evaluate site-specific conditions, users may select from among a number of available environmental transport and exposure pathways and modify 'default' parameter values. In GoldSim, a user must construct the model 'from scratch.' However, this modeling is facilitated by a number of specialized elements available in GoldSim to support a radiological mass transport model. These include elements for defining radiological decay and ingrowth, container failure and radiological release, advective and diffusive transport, and other processes. Advantages afforded by the use of RESRAD for modeling radiological dose for a RCRA Subtitle C landfill include ease of use and reasonable flexibility in specifying site-specific conditions. GoldSim allows for considerably more flexibility and site-specificity than RESRAD, including evaluation of potentially relevant environmental transport processes not supported in RESRAD. The probabilistic modeling capabilities of GoldSim also far exceed those of RESRAD. The pros and cons of RESRAD and GoldSim for this modeling problem are explored with a focus on identifying approaches and critical factors in identifying the appropriate platform. The capabilities of RESRAD and GoldSim for mathematically modeling the disposal system will be contrasted and compared. For example, the existing RESRAD radiological safety assessment evaluated potentially complete exposure pathways related to infiltration to groundwater and upwards diffusion of radon. Additional transport pathways identified in the USEI Idaho facility Conceptual Site Model include deposition of radon decay products in cover material, root uptake of radionuclides in disposed wastes by native plants, and mixing of cover material by burrowing animals. Results of the two models and how the different transport and dose pathways affect results will be discussed. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗