Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Material balance time”

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 343 records · Page 19

Extinction of the Stagnation Point Diffusion Flame: Effect of Conductive Heat Loss into Solid Interior

The flammability boundary of a solid is U-shaped with the ordinate as ambient oxygen percentage (or total pressure) and the abscissa represents flow velocity, gravity level, stretch rate, etc. The left side of the flammability boundary is the quenching branch, and the right side is the blow-off branch. Quenching is due to a weak flame whose heat release from combustion cannot overcome the excessive heat losses from, for example, radiation. Blow-off is due to a short gas residence time that is smaller than the chemical reaction time (traditionally referred to as the small Damköhler number effect). For thick solids, the heat conduction into the solid interior has a long transient after ignition before the solid temperature reaches the steady state. During this slow transient, there is excessive (i.e. more than that at the steady state) gas phase heat conduction into the solid interior. From the viewpoint of the quasi-static gaseous flame, this is an additional heat loss that will affect the gas flame's flammability. In this work, a previously developed one-dimensional axisymmetric stagnation-point diffusion flame model is extended to account for the excessive heat loss into the solid beyond the steady burning state. A non-dimensional excessive conductive heat loss parameter Ψ is defined in the paper to account for the transient solid effect. For each Ψ, the gas phase is considered reaching a quasi-steady state. Extensive numerical computations have been performed to determine the extinction boundary and the characteristics of near-limit flames. The U-shaped flammability boundary is now a series of boundaries with Ψ as a parameter. Flame structures under the same environment condition but with different Ψ are compared. For a given oxygen ambient, the extinction boundary Ψ vs stretch rate is an inverted U-shape. The boundary consists of a blowoff branch and a radiative quenching branch. Extensive computed data including flame and pyrolysis temperatures, burning rate, flame standoff distance, species concentration, and reaction rate are presented along the extinction boundaries. Heat balance analysis is performed for both the gas phase and the gas-solid interface. Along the extinction boundary, the relative importance of conduction, convection, radiation, and chemical reaction are analyzed in detail.

Stagnation point diffusion flame↗

Adsorption of rare earth elements in carboxylated mesoporous carbon

The separation of Rare Earth Elements (REEs) from various non-conventional sources is critically important to maintain the supply–demand balance of REEs in the western world. Here, in this research, we have investigated coal Fly Ash (CFA) as the non-conventional source of REEs. In order to recover REEs, we have synthesized carboxylate-functionalized mesoporous carbon (CMC) to harness the chelation-induced adsorption of REEs. These adsorbents were characterized with pore textural properties, FTIR spectroscopy, and SEM-EDX mapping. In the initial study with single-component La(III) Dy(III) and Lu(III), it was revealed that CMC can adsorb those REEs 2–4 times more than that of pristine mesoporous carbon confirming the effective role of the carboxylate group in the REE separation. Furthermore, the affiliation of the carboxylate group towards heavier REEs (Lu>Dy>La, when present is equal amounts) provides an added benefit of this adsorbent owing to the high demand for heavier REEs in the technology sectors. The leachate solution produced from CFA contained sixteen REEs in the range of 50–200 ppb for most of the REEs. It was revealed that CMC can extract 80–90% of all the REEs thereby further confirming the success of the CMC as an effective REE sorbent. Three successive cycles of adsorption and desorption of REEs with the same CMC revealed the consistent adsorption capacity of REEs. In-situ X-ray Absorption Near Edge Spectroscopy (XANES) analysis confirmed the +3 oxidation states of La, Dy, and Lu within CMC. In-situ Extended X-ray Absorption Fine Structure (EXAFS) analysis revealed the shortening of mean La-O, Dy-O, and Lu-O bond distance by 0.03–0.05 Å thereby confirming the coordination of these REEs with carboxylate groups present on CMC.

36 MATERIALS SCIENCE↗

Crystal structure and magnetic properties in semiconducting Eu 3-δ Zn x Sn y As 3 with Eu-Eu dimers

Magnetic structure and crystal symmetry, which primarily determine the time-reversal and inversion symmetry, may give rise to numerous exotic quantum phenomena in magnetic semiconductors and semimetals when arranged in different patterns. Here, a new layered magnetic semiconductor, Eu 3-δ Zn x Sn y As 3 , was discovered and high-quality single crystals were grown using the Sn flux. According to structural characterization by x-ray diffraction and atomic-resolution scanning transmission electron microscopy, Eu 3-δ Zn x Sn y As 3 is found to crystallize in a hexagonal symmetry with the space group P6 3 /mmc (No. 194). After examining different specimens, we conclude that their stoichiometry is fixed at ~Eu 2.6 Zn 0.65 Sn 0.85 As 3 , which meets the chemical charge balance. Eu 3-δ Zn x Sn y As 3 is composed of septuple (Eu 1-δ Sn y As 2 )-Eu-(Zn x As)-Eu sequences. The shortest Eu–Eu distance in the system is between two Eu layers separated by Zn x As along the c-axis. Magnetization measurement shows an antiferromagnetic ordering in Eu 3-δ Zn x Sn y As 3 at T N ~ 12 K, where the magnetic easy-axis is along the c-axis, and Mössbauer spectroscopy observes magnetic hyperfine splitting on Eu and Sn at 6 K. Magnetic anisotropy is significantly different from the ones along the ab-plane in other layered Eu-based magnetic semimetals. Heat capacity measurements confirm the magnetic transition around 12 K. Electrical resistivity measurement indicates semiconductor behavior with a band gap of ~0.86 eV. Finally, various Eu-based magnetic semiconductors could provide a tunable platform to study potential topological and magnetic properties.

36 MATERIALS SCIENCE↗

Modeling Molten Salt Spreading and Heat Transfer using MELTSPREAD – An Uncertainty Analysis

This report summarizes spreading and heat transfer results for molten salt pouring onto a flat stainless steel substrate that were calculated using MELTSPREAD to gain insight into bulk salt behavior during a molten salt spill accident. MELTSPREAD was developed at Argonne to model the one-dimensional flow and freezing behavior of molten corium and is being applied to model molten salts for the first time as part of this work. An uncertainty analysis of thermophysical properties and spill conditions was performed using MELTSPREAD to determine which properties and conditions have the greatest impact on model outcome. The model was run with and without the inclusion of decay heat. Eutectic FLiNaK was used as the salt composition because it is well-characterized and appropriate for model development. Small salt volumes were used for the initial model runs described in this report so that the model results can be compared to the results from ongoing experiments on molten salt spreading and heat transfer currently being conducted at Argonne at a benchtop scale. The spreading of a small volume of FLiNaK (50 mL) was found to be limited by the balance between the molten salt surface tension and gravity and not by freezing for the pour conditions and assumed mechanisms of heat transfer in the model. Heat transfer from the salt was highly inefficient to due to its low thermal conductivity, high heat capacity, and high heat of fusion and varying these properties had no effect on the spreading behavior. Changing the viscosity and density of the salt slightly affected the spreading behavior. Setting the salt surface tension to approximately zero to simulate the salt wetting the substrate had the greatest effect on the spreading behavior and this increased both the spreading velocity and overall spreading area. Including decay heat in the model had a significant effect on the long-term cooling behavior of the salt but had no effect on the overall spreading behavior. Future effort should involve modeling a larger volume of molten salt to represent a reactor case spill scenario because the spreading behavior of large volumes of salt may differ from the small volumes that were the focus of this document.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Mercury Speciation During Vitrification of LAW

n this work, laboratory and engineering-scale tests were conducted that mimic the reactive environment of the Hanford Tank Waste Treatment and Immobilization Plant (WTP) low activity waste (LAW) melter and off-gas system in order to assess the speciation of mercury at select points in the processing system. The experimental protocols, test equipment, feed materials used, and rationale for testing are detailed in the Test Plan for this work. Particular attention was paid to the amount and speciation of mercury in submerged bed scrubber (SBS) solutions, which are intended to be recycled back to the melter feed at the WTP. Results are presented from detailed mercury analysis of these solutions using Environmental Protection Agency (EPA) Method 1630 and modifications to this method for dimethyl mercury and methyl mercury as well as EPA Method 1631 and modifications to this method for total, dissolved, elemental, suspended, and ionic mercury. Testing was conducted in two separate sets of experiments: (1) One set of tests (crucible-scale furnace tests) that involved heating small batches of mercury-spiked melter feed in crucibles to determine the amount of mercury retained in the glass and the species of mercury in the exhaust gases; (2) A second set of tests (DM10 melter tests) was planned that involved creating LAW melter plenum gas compositions using the DM10 melter system, injecting mercury into the off-gas stream, and passing that stream though a reactor that simulates various plenum gas conditions. Operational issues led to the need to use one of the VSL DM100 melters in place of the planned DM10 melter for these tests. The speciation of mercury after exposure to those conditions was monitored. In both sets of tests, the exhaust gases were run through a scrubber that was intended to mimic the LAW SBS in order to determine how much of the mercury exiting the melter would be retained in the primary off-gas system fluids and in what form. After passing through the SBS, the exhaust stream was analyzed to determine particulate, ionic, elemental, and total mercury passing downstream of the SBS. In tests employing gases derived from the DM100 melter that were spiked with elemental mercury, the processing system provided gas temperatures and residence times that are representative of the WTP LAW vitrification system in order to assess the effect of those conditions on mercury speciation. The Decontamination Factor (DF) across the system in the tests with mercury-spiked DM100 melter exhaust was determined using the analytical data from EPA Method 30B (Fluegas Adsorption Mercury Speciation (FAMS TM )) exhaust samples and the amount of mercury detected in the SBS solutions. Mercury species used in melter feed crucible scale tests were divalent (chloride and iodide), monovalent (chloride and fluoride), and elemental mercury. Individual tests included only a single form of mercury in the feed. The results for total mercury mass balances in the crucible tests are also presented.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Hydrogeological assessment of CO2 containment assurance and wellbore integrity at a Gulf Coast storage site

Abstract A large-scale carbon capture and storage (CCS) initiative on the Texas Gulf Coast serves as a premier demonstration of the U.S. Department of Energy’s CarbonSAFE program. Targeting deep saline formations, specifically Oligo-Miocene deltaic sequences, the project aims to establish technical and commercial viability for geologic CO2 storage within a major industrial corridor. This study provides a rigorous hydrogeological assessment to support Class VI permitting by quantifying the high degree of containment security. Utilizing a compositional reservoir simulator, we developed a suite of 27 distinct simulation cases to evaluate vertical plume dynamics near both planned injection wells and proximal legacy infrastructure. To ensure numerical accuracy near wellbores, we implemented a refined mesh strategy, determining that a 5.6 ft × 5.6 ft grid refinement offered the optimal balance between computational efficiency and descriptive precision. The modeling framework utilized a systematic sensitivity-based approach to evaluate the mechanical redundancy of the subsurface system by performing a bounding analysis of wellbore interfaces against hypothetical high-permeability microannuli. By systematically isolating competing physical drivers, including permeability, porosity, gas hysteresis, thermal gradients, salinity, and solubility trapping (quantified via Henry’s law with dynamically adjusted coefficients), this work moves beyond binary assessments to establish a nuanced hierarchy of containment factors. The results confirm that primary trapping mechanisms (e.g., gas hysteresis and solubility), combined with the site's unique geomechanical stratigraphy, significantly restrict vertical mobility and reinforce the robust containment security of the reservoir. Baseline results demonstrate substantial vertical separation between the CO2 plume and the upper confining system, ensuring robust containment. Sensitivity analysis reveals that even under highly conservative bounding scenarios—assuming theoretical 10-Darcy pathways at specific wellbore locations—the 2,900-ft thick multi-layered confining zone remains a reliable barrier. In these hypothetical upper-bound cases, peak upward fluxes of CO2 and saltwater after 15 years of injection remain localized and dissipate rapidly within the lower sections of the confining interval, leaving the integrity of the seal uncompromised. Furthermore, the study identifies that while localized wellbore pathways define theoretical upper bounds of vertical migration, the Area of Review (AoR) is primarily sensitive to regional thermal gradients and hysteresis, which can influence the AoR by over 3,000 acres in pessimistic configurations. Also, primary trapping mechanisms, specifically gas hysteresis and solubility, work in tandem with the Gulf Coast’s unique geomechanical stratigraphy to significantly restrict vertical mobility. Ductile, smectite-rich mudstones facilitate natural borehole convergence and the self-healing of potential conduits, creating a natural geomechanical bridge that effectively mitigates migration potential at both current injection points and legacy-well locations. This comprehensive modeling effort demonstrates that the integration of high-resolution wellbore simulations and regional geomechanical observations confirms the long-term storage security of the studied site, providing a physics-based foundation for industrial-scale CCS deployments. This modeling framework establishes a baseline for future research into coupled geomechanical effects, such as time-dependent borehole convergence, to further refine long-term containment projections. Acknowledgements We thank the Gulf Coast Carbon Center (GCCC) at the Bureau of Economic Geology for foundational research support. We appreciate Alex Bump for technical guidance and David Hoffman for model mesh generation. This work used TACC’s Frontera cluster for simulations and CMG Ltd. software licenses provided to UT-Austin. This material is based upon work supported by the Department of Energy under Award Number DE-FE0032338. Disclaimer This material is based upon work supported by the U.S. Department of Energy’s Fossil Energy and Carbon Management Office under the CarbonSAFE program, award Number DE-FE0032338. The views expressed herein do not necessarily represent the views of the U.S. Department of Energy or the United States Government.

58 GEOSCIENCES↗

Hydrogeological assessment of CO2 containment assurance and wellbore integrity at a Gulf Coast storage site

Abstract A large-scale carbon capture and storage (CCS) initiative on the Texas Gulf Coast serves as a premier demonstration of the U.S. Department of Energy’s CarbonSAFE program. Targeting deep saline formations, specifically Oligo-Miocene deltaic sequences, the project aims to establish technical and commercial viability for geologic CO2 storage within a major industrial corridor. This study provides a rigorous hydrogeological assessment to support Class VI permitting by quantifying the high degree of containment security. Utilizing a compositional reservoir simulator, we developed a suite of 27 distinct simulation cases to evaluate vertical plume dynamics near both planned injection wells and proximal legacy infrastructure. To ensure numerical accuracy near wellbores, we implemented a refined mesh strategy, determining that a 5.6 ft × 5.6 ft grid refinement offered the optimal balance between computational efficiency and descriptive precision. The modeling framework utilized a systematic sensitivity-based approach to evaluate the mechanical redundancy of the subsurface system by performing a bounding analysis of wellbore interfaces against hypothetical high-permeability microannuli. By systematically isolating competing physical drivers, including permeability, porosity, gas hysteresis, thermal gradients, salinity, and solubility trapping (quantified via Henry’s law with dynamically adjusted coefficients), this work moves beyond binary assessments to establish a nuanced hierarchy of containment factors. The results confirm that primary trapping mechanisms (e.g., gas hysteresis and solubility), combined with the site's unique geomechanical stratigraphy, significantly restrict vertical mobility and reinforce the robust containment security of the reservoir. Baseline results demonstrate substantial vertical separation between the CO2 plume and the upper confining system, ensuring robust containment. Sensitivity analysis reveals that even under highly conservative bounding scenarios—assuming theoretical 10-Darcy pathways at specific wellbore locations—the 2,900-ft thick multi-layered confining zone remains a reliable barrier. In these hypothetical upper-bound cases, peak upward fluxes of CO2 and saltwater after 15 years of injection remain localized and dissipate rapidly within the lower sections of the confining interval, leaving the integrity of the seal uncompromised. Furthermore, the study identifies that while localized wellbore pathways define theoretical upper bounds of vertical migration, the Area of Review (AoR) is primarily sensitive to regional thermal gradients and hysteresis, which can influence the AoR by over 3,000 acres in pessimistic configurations. Also, primary trapping mechanisms, specifically gas hysteresis and solubility, work in tandem with the Gulf Coast’s unique geomechanical stratigraphy to significantly restrict vertical mobility. Ductile, smectite-rich mudstones facilitate natural borehole convergence and the self-healing of potential conduits, creating a natural geomechanical bridge that effectively mitigates migration potential at both current injection points and legacy-well locations. This comprehensive modeling effort demonstrates that the integration of high-resolution wellbore simulations and regional geomechanical observations confirms the long-term storage security of the studied site, providing a physics-based foundation for industrial-scale CCS deployments. This modeling framework establishes a baseline for future research into coupled geomechanical effects, such as time-dependent borehole convergence, to further refine long-term containment projections. Acknowledgements We thank the Gulf Coast Carbon Center (GCCC) at the Bureau of Economic Geology for foundational research support. We appreciate Alex Bump for technical guidance and David Hoffman for model mesh generation. This work used TACC’s Frontera cluster for simulations and CMG Ltd. software licenses provided to UT-Austin. This material is based upon work supported by the Department of Energy under Award Number DE-FE0032338. Disclaimer This material is based upon work supported by the U.S. Department of Energy’s Fossil Energy and Carbon Management Office under the CarbonSAFE program, award Number DE-FE0032338. The views expressed herein do not necessarily represent the views of the U.S. Department of Energy or the United States Government.

58 GEOSCIENCES↗

Contaminant Migration Modeling to Support the In Situ Decommissioning of Hardened Nuclear Facilities - 20215

The decommissioning of hardened nuclear facilities provides a unique challenge in balancing current worker dose with exposure to future receptors and media. In situ decommissioning (ISD) is a cost-effective and safe option for the closure of such facilities but requires a detailed understanding of the threat that radioactive and hazardous constituents remaining in the facilities pose to groundwater and surrounding surface water. The Savannah River Site (SRS) extensively employs contaminant migration models within the decommissioning process to develop site-specific removal, grouting, and monitoring strategies in support of safe and effective final end states of nuclear facilities. This work discusses the use of contaminant migration modeling in the ISD process at SRS, including specific examples from the closure of P-Reactor and R-Reactor, Building 235-F, F-Canyon and F/H Laboratory Complex. Migration models developed during the closure of P-Reactor and R-Reactor were divided into four source areas: the reactor vessel, the process area, the disassembly basin, and the purification wing. Each source area was assigned a specific inventory and migration pathway and was then subject to varying hypothetical removal, capping, and grouting schemes. Modeling indicated that groundwater Maximum Concentration Limits (MCLs) might be exceeded if no action was taken for eleven and ten constituents of concern (COCs) at P-Reactor and R-Reactor, respectively. These exceedances could occur in as few as 200 to 500 years. Alternatively, the migration modeling demonstrated that the selected ISD actions reduced contaminant mobility which allowed for significant radioactive decay and resulted in fewer predicted exceedances of groundwater MCLs (five COCs for P-Reactor; eight COCs for R-Reactor). In response to the modeling results for P-Reactor and R-Reactor, effectiveness monitoring programs were developed to target contaminants, specific to each reactor, that may migrate to groundwater. Contaminant migration modeling also revealed that roof collapse was a large factor in the release of COCs to the environment, giving rise to roof improvements, and an inspection and vegetation control program to ensure roof stability over time. Contaminant migration modeling is also aiding in the closure planning for hardened facilities in F Area. Building 235-F housed the Actinide Billet Line, which produced Np-237 billets for irradiation in SRS reactors, and the Plutonium Fuel Form (PuFF) facility that produced Pu-238 heat sources for the space program. As a result of these missions, areas within Building 235-F contain considerable residual amounts of both Pu-238 and Np-237. Contaminant migration modeling of Building 235-F was originally performed in 2012 to identify the feasibility of ISD and the amount of radioactive material removal required to prevent the exceedance of groundwater MCLs. The original model indicated that a 60% reduction in the PuFF facility Pu-238 inventory could keep groundwater concentrations below standards, while Np-237 did not pose a threat to groundwater. However, updates to the model with an emphasis on source impact pathways revealed that, due to the orientation of the source areas relative to groundwater flow, no amount of reasonable removal of Pu-238 would keep groundwater concentrations below MCLs and that Np-237 could be a large contributor to localized MCL exceedances. With this insight, the refined 2019 model is being used to assist in the development of grouting plans specific to each facility source area, where bentonite may be utilized to slow the migration of Pu-238 and its daughter products from the PuFF facility and a reducing grout may decrease Np-237 transport by ensuring the nuclide remains in the less mobile +IV oxidation state. The beginning phases of contaminant migration modeling are underway for F-Canyon and associated facilities using lessons learned from the Reactors and Building 235-F. A contaminant migration pathway, similar to the pathway used for the reactor vessels, is being developed for the hot and warm canyons. The F-Canyon inventories are being spatially refined to identify specific and localized source zones, comparable to Building 235-F, that may impact groundwater. Contaminant migration modeling provides vital input within the ISD process, from facility investigation to post closure effectiveness monitoring, making it a valuable tool in the development of safe and effective end states for hardened nuclear facilities at SRS. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Distribution Development for the RDX Regional Model at Los Alamos National Laboratory - 20385

Representing uncertainty in model inputs often means finding a balance between uncertainty and physical reality. Developing wide distributions may seem conservative in principle, but this approach may lead to unrealistic model results. Characterizing the current state of knowledge of stochastic inputs presents many challenges, especially if the data available are limited or have limited relevance to the site. Relationships among these inputs may also be important to represent but are typically complex or difficult to define. Often special adjustments must be made to account for reduced credibility in particular data. If parameters are strongly related to one another, a correlation structure may be developed for input to the model. Other techniques such as regression models may be used to incorporate relationships between the information available and the desired parameters. The process of developing distributions must consider details of the model in terms of what the distribution is meant to represent. This paper uses the example of a probabilistic fate and transport model for hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX) in the regional aquifer at Los Alamos National Laboratory (LANL). For many parameters, a single draw is applied to all space and time over which the model is run, for a single iteration. This simplification is often made for many reasons, and can often be beneficial, but also adds additional complexity in the distribution development process. Defining the distributional goals as they relate to the modeling process is an important step, which should take place prior to evaluation of the data. Usually, the distributions developed are meant to characterize the average value of the parameter over the spatial and temporal domain of the model. Distribution development requires consideration of many sources of information on the parameter where available, ideally from multiple references. Examples of different sources include data from different references but also from different conditions, measurement methods, or experimental types. Depending on these conditions and the reliability or relevance of particular references, different sources of data may each contribute valuable information but have varying relevance to the site. In these cases, weighting data unequally is a useful way to incorporate this information. As an example, aqueous dispersivity data are available for a variety of rock types. Only a few values are available for the desired rock type, and this is not enough to develop a distribution. Therefore, dispersivity values from other rock types are included in distribution development but are down-weighted such that the best data have the most influence on the distribution developed. In another case, K{sub d} distributions in the model are meant to represent a known composition of multiple soil types. Data from these materials are weighted accordingly to develop a distribution for the weighted average K{sub d} across soil types. Other cases include varying reliability of different sources, and weighting data according to the confidence in these sources. In some cases, input parameters are correlated with one another. An example is advective porosity, which is positively correlated with total porosity and must be less than total porosity. Paired data with both parameters must exist to discern the relationship between the two parameters and if it is necessary to build a correlation structure into the model. In general, correlated parameters may be represented in the model by a multivariate distribution, or perhaps more desirably, capturing the correlation within the developed distributions which may be treated as independent from one another. In the example of porosity, this can be done by transforming advective porosity into a proportion of total porosity which may be drawn independently from the distribution of total porosity. This paper explores how complex data and correlations can be incorporated to meet the distributional goals of the model, using the RDX regional model as a detailed example. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Computer Modeling of Electrostatic Aggregation of Granular Materials in Planetary and Astrophysical Settings

Electrostatic forces strongly influence the behavior of granular materials in both dispersed (cloud) systems and semi-packed systems. These forces can cause aggregation or dispersion of particles and are important in a variety of astrophysical and planetary settings. There are also many industrial and commercial settings where granular matter and electrostatics become partners for both good and bad. This partnership is important for human exploration on Mars where dust adheres to suits, machines, and habitats. Long-range Coulombic (electrostatic) forces, as opposed to contact-induced dipoles and van der Waals attractions, are generally regarded as resulting from net charge. We have proposed that in addition to net charge interactions, randomly distributed charge carriers on grains will result in a dipole moment regardless of any net charge. If grains are unconfined, or fluidized, they will rotate so that the dipole always induces attraction between grains. Aggregates are readily formed, and Coulombic polarity resulting from the dipole produces end-to-end stacking of grains to form filamentary aggregates. This has been demonstrated in USML experiments on Space Shuttle where microgravity facilitated the unmasking of static forces. It has also been demonstrated in a computer model using grains with charge carriers of both sign. Model results very closely resembled micro-g results with actual sand grains. Further computer modeling of the aggregation process has been conducted to improve our understanding of the aggregation process, and to provide a predictive tool for microgravity experiments slated for Space Station. These experiments will attempt to prove the dipole concept as outlined above. We have considerably enhanced the original computer model: refinements to the algorithm have improved the fidelity of grain behavior during grain contact, special attention has been paid to simulation time steps to enable establishment of a meaningful, quantitative time axis, and calibration of rounding accuracies have been conducted to test cumulative numerical influences in the model. The model has been run for larger grain populations, variable initial cloud densities, and we have introduced random net charging to individual grains, as well as a net charge to the cloud as a whole. The model uses 3 positive and 3 negative charges randomly distributed on each grain, with up to 160 grains contained within various size "boxes" that define the initial number densities in the clouds. Each charge represents localized charged region on a grain, but does not necessarily imply single quantized charge carriers. The Coulomb equations are then allowed to interact for each monopole: dipoles and any higher order charge coupling is a natural product of these "free" interactions over which the modeler exerts no influence. The charges are placed on surfaces of grains at random locations. A series of runs was conducted for neutral grains that had a perfect balance of negative and positive char carriers. Runs were also conducted with grains having additional fractional charges ranging between 0 and 1. By adding fractional charges of one sign, the model created grain populations in which all grains had excess charges the same sign, giving the cloud an overall net charge. This simulates clouds subjected to ionizing radiation (e. protoplanetary debris disk around a protosun), or any other process of charge biasing in a grain population (e.g., volcanic plumes). In another run series, random fractional charges of either sign were added to the grains so th some grains had a slight net positive charge while others had a slight net negative charge. This simulates triboelectrically-charged grain populations in which acquisition of an electron by one surface is at the expense creating a hole elsewhere. This dual sign charging was applied in two ways: in one case the cloud remained neutral by ensuring that all grain excess charges added to zero; in the other case, the cloud was permitted slight net char by not imposing a charge-balance condition. Additional information is contained in the original.

Marshall, J.↗

Comparison of Extended and Accelerated VOG Tests

The reprocessing of used nuclear fuel would release volatile radionuclides into the off-gas streams of a processing plant, including ¹²⁹I. The dissolver off-gas and more dilute vessel off-gas streams (VOG) must both be targeted when mitigating ¹²⁹I environmental release because each contains some amount of iodine. Iodine-129 will likely be found as elemental iodine (I₂) and methyl iodide (CH₃I) in these off-gas streams. Reduced silver-exchanged mordenite (AgZ or Ag⁰Z) has been investigated as a potential sorbent for iodine abatement and is studied here under prototypic VOG conditions. Because of the relatively low iodine abundance and high flow rates of the VOG, total iodine concentrations are expected to be in the parts per billion range. Thus, VOG experiments need to run for extended durations at low concentrations to reach sorbent saturation. Because the sorbent will be exposed to oxidizing gas streams for extended periods of time, aging effects and sorbent degradation need to be considered when designing a sorbent-based abatement system. Three sets of experiments were performed with the aim of determining how the adsorption of iodine by AgZ is affected by differing test durations and gas compositions. The first tested the capacity of sorbent aged for 8 months under a humid air stream. The second tested differences in sorbent capacity and mass transfer zone (MTZ) length during high-concentration (1200 ppb v ) CH₃I loading in a humid nitrogen gas (N₂) stream and a humid air gas stream over 28 days. The third tested sorbent capacity and MTZ length during low concentration (< 200 ppb v ) CH₃I and I₂ loading in a humid air stream over 9 months. The results of these tests suggest that AgZ capacity drops by ~50% after 1 month of aging, by ~60% after 2 months of aging, and then does not continue to decrease significantly at up to 8 months of aging. One-month aging tests that compared N₂ and air as the gas stream diluent resulted in similar maximum loading capacities and overall loading curves, indicating that the effects of aging cannot be mitigated by avoiding air as the balance gas. The 9-month extended VOG tests did not result in clear sorbent saturation at the inlet, but it can be inferred using an assumption of a maximum capacity of 45 mg I/g AgZ. Applying this assumption, then the MTZ of CH₃I is 12.8 cm, and the MTZ of I₂ is 12.4 cm. These results are similar to previous tests. Comparisons to tests completed at Idaho National Laboratory suggest that the CH₃I MTZ may be dependent on concentration. Future work should re-evaluate the design of the VOG iodine capture system with this updated data and should seek to understand fundamental characteristics of CH₃I adsorption by AgZ, such as the effect of concentration, gas velocity, and gas composition.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Hybrid HEFA-HDCJ Process for the Production of Jet Fuel Blendstocks

The hydrotreatment of bio-oil derived from the pyrolysis and biocrude from hydrothermal liquefaction of lignocellulosic materials to produce hydrocarbons faces significant technological challenges, mainly due to the high reactivity and poor thermal stability of bio-oil, resulting in the formation of large quantities of coke. This problem has been addressed by existing PNNL patents with a two-step hydrotreatment technology in which the bio-oil is first stabilized with a noble hydrogenation metal (often Pt or Ru). Then, in the second step, the bio-oil is deoxygenated with a Ni-Mo or Co-Mo sulfide catalyst. The main problem with this approach is that the Pt/Ru catalysts deactivate easily in the presence of S or other impurities, which are commonly present in pyrolysis oils. In this project, we explored technological solutions to mitigate coke formation, avoiding the use of Pt/Ru catalysts. Our strategy is based on three actions: (1) Bio-oil stabilization in the presence of alcohols. In this project, we studied the stabilization with butanol. (2) the use of a cosolvent to solubilize the bio-oil. Because coke formation reactions are second-order reactions a reduction in the concentration of reactive bio-oil molecules. In this case, we used yellow greases as a co-solvent. (3) Separation of bio-oil reactive fractions. In this project, we studied the removal of water-soluble fractions. Our batch co-hydrotreatment studies confirmed that the addition of butanol and methanol and the blend with lipids effectively contributed to mitigating coke formation (reducing coke yield to about 1 wt.% %). The removal of sugars did not have a noticeable effect on the overall coke yield, suggesting that coke precursors are present in all bio-oil fractions. Our analytical work suggests that they may be concentrated in the water-insoluble/CH 2 Cl 2 insoluble fractions of pyrolysis oils. Although the technological strategies tested resulted in significant coke reductions, the levels achieved were not sufficiently low to ensure a reliable operation in continuous, fixed-bed trickle-bed reactors. Long runs of more than 100 hours (maximum: 255 h) of co-processing time on stream were achieved in a continuous 40 mL reactor. When the same test was conducted in a larger 400 mL reactor, pressure drop increases associated with coke formation were observed. This increase in coke formation could be due to larger temperature gradients in the bed. Hydrodeoxygenation tests in moving bed reactors and using more active hydrogenation catalysts (for example Ni) could lead to more reliable operations. Unfortunately, our team did not have access to such experimental setups. The technoeconomic analysis suggests that although alcohol use is an effective means to reduce coke formation, the use of alcohol increases production cost. Thus, its use needs to be minimized. A delicate balance needs to be found between the use of technological solutions that allow the reliable operation of the system (stabilization with Ni catalysts, use of small quantities of solvents, processing in moving bed reactors) with a tolerable level of coke formation for the hydrodeoxygenation reactor used and that result in minimum production costs.

09 BIOMASS FUELS↗

Minimizing Cr-Evaporation from Balance of Plant Components by Utilizing Cost-Effective Alumina-Forming Austenitic Steels

A solid oxide fuel cell (SOFC) is a clean and efficient energy conversion device. The development of intermediate-temperature SOFCs has made it preferable to use metallic interconnects (MICs) to greatly reduce the cost and significantly increase the efficiency compared to ceramic interconnect materials. However, gaseous chromium species will evaporate from the chromium-containing layer formed on the surface of commonly used MICs and balance of plant (BoP) components. Volatile chromium species have been shown to form solid deposits which poison the cathodes of SOFCs, causing drastic cell performance degradation and thereby limiting commercialization. In order to alleviate the Cr poisoning and achieve long-term high performance of SOFC stacks, various Al2O3-forming austenitic (AFA) stainless steels applied at different temperatures are evaluated in this project. Based on our Phase I results (500 h operation), it is shown that on the AFAs, an alumina-based protective layer forms under high temperature that is invulnerable to water vapor effects and suppresses the diffusion of chromium and manganese which can prevent the generation of spinels on the alloy surface. Therefore, 310S, OC4 and OC5 at 800 °C and 625, OC11 and OC11LZ at 900 °C are selected to be tested in the Phase II long-term operation. From the results, we find that besides the lower Cr evaporation rates and better oxidation resistance of AFAs than benchmark alloys after short-term (500 h) operation, AFAs also possess the sturdy and compact alumina layer after a long-term operation (5000 h). A protective oxide layer is of great importance for the long-term high-temperature operation of structural materials. The high-temperature oxidation behavior of AFAs and commercial 310S and 625 alloys in Air + 10% H 2 O at 800 °C and 900 °C is systematically investigated. Severe breakaway oxidation and minor spallation are observed for the 310S and 625 after the long-term operation, while the AFAs show high oxidation resistance. It is found that the formation of a continuous alumina layer could greatly prevent the volatilization of chromium vapors. In addition, the corresponding evolving models are discussed. Chromium evaporation from BoP components in high-temperature environment could severely deteriorate the electrochemical performance of SOFC. Several methods were applied to evaluate the Cr evaporation rates of BoP after 500 h exposure at 800 °C to 900 °C in air with 10% H 2 O. An optimal method was designed to exclude the effect of silicon (Si) deposits from quartz tube and sodium (Na) deposits from the sodium carbonate on the oxidation process and the chemical interaction between Cr gaseous species and alumina tube which could provide further quantitative correlation of the evaporated Cr species quantities and degradation rates of SOFC. Based on the great performance of AFAs after long-term operation, AFAs are assembled with Anode-supported cells (ASC) to investigate the the Cr deposition of anode-supported cell under a constant current density of 0.5 A cm -2 at 800 °C with AFA alloys compared with commercial alloys. In addition, the anodic and cathodic processes are deconvoluted by distribution of relaxation times (DRT) method which are comprehensively discussed. It is found out that the voltage exhibited a slight decrease of 5.09 % and 1.54 % in the presence of OC11 and OC11LZA alloy, respectively. However, a considerable decrease of 22.14 % and 12.06 % was determined in the presence of 310S and 625 alloy, respectively. AFAs possessing the low-cost, low Cr evaporation rates and the high oxidation resistance will make it of great potential to replace the existing BoP components.

36 MATERIALS SCIENCE↗

Minimizing Cr-Evaporation from Balance Of Plant Components By Utilizing Cost-Effective Alumina-Forming Austenitic Steels

A solid oxide fuel cell (SOFC) is a clean and efficient energy conversion device. The development of intermediate-temperature SOFCs has made it preferable to use metallic interconnects (MICs) to greatly reduce the cost and significantly increase the efficiency compared to ceramic interconnect materials. However, gaseous chromium species will evaporate from the chromium-containing layer formed on the surface of commonly used MICs and balance of plant (BoP) components. Volatile chromium species have been shown to form solid deposits which poison the cathodes of SOFCs, causing drastic cell performance degradation and thereby limiting commercialization. In order to alleviate the Cr poisoning and achieve long-term high performance of SOFC stacks, various Al 2 O 3 -forming austenitic (AFA) stainless steels applied at different temperatures are evaluated in this project. Based on our Phase I results (500 h operation), it is shown that on the AFAs, an alumina-based protective layer forms under high temperature that is invulnerable to water vapor effects and suppresses the diffusion of chromium and manganese which can prevent the generation of spinels on the alloy surface. Therefore, 310S, OC4 and OC5 at 800 °C and 625, OC11 and OC11LZ at 900 °C are selected to be tested in the Phase II long-term operation. From the results, we find that besides the lower Cr evaporation rates and better oxidation resistance of AFAs than benchmark alloys after short-term (500 h) operation, AFAs also possess the sturdy and compact alumina layer after a long-term operation (5000 h). A protective oxide layer is of great importance for the long-term high-temperature operation of structural materials. The high-temperature oxidation behavior of AFAs and commercial 310S and 625 alloys in Air + 10% H 2 O at 800 °C and 900 °C is systematically investigated. Severe breakaway oxidation and minor spallation are observed for the 310S and 625 after the long-term operation, while the AFAs show high oxidation resistance. It is found that the formation of a continuous alumina layer could greatly prevent the volatilization of chromium vapors. In addition, the corresponding evolving models are discussed. Chromium evaporation from BoP components in high-temperature environment could severely deteriorate the electrochemical performance of SOFC. Several methods were applied to evaluate the Cr evaporation rates of BoP after 500 h exposure at 800 °C to 900 °C in air with 10% H 2 O. An optimal method was designed to exclude the effect of silicon (Si) deposits from quartz tube and sodium (Na) deposits from the sodium carbonate on the oxidation process and the chemical interaction between Cr gaseous species and alumina tube which could provide further quantitative correlation of the evaporated Cr species quantities and degradation rates of SOFC. Based on the great performance of AFAs after long-term operation, AFAs are assembled with Anode-supported cells (ASC) to investigate the the Cr deposition of anode-supported cell under a constant current density of 0.5 A cm -2 at 800 °C with AFA alloys compared with commercial alloys. In addition, the anodic and cathodic processes are deconvoluted by distribution of relaxation times (DRT) method which are comprehensively discussed. It is found out that the voltage exhibited a slight decrease of 5.09 % and 1.54 % in the presence of OC11 and OC11LZA alloy, respectively. However, a considerable decrease of 22.14 % and 12.06 % was determined in the presence of 310S and 625 alloy, respectively. AFAs possessing the low-cost, low Cr evaporation rates and the high oxidation resistance will make it of great potential to replace the existing BoP components.

01 COAL, LIGNITE, AND PEAT↗

Interstellar PAH in the Laboratory and in Space. What have we Learned from the New Generation of Laboratory and Observational Studies?

Polycyclic Aromatic Hydrocarbons (PAHs) are an important and ubiquitous component of carbon-bearing materials in space. PAHs are the best-known candidates to account for the IR emission bands (UIR bands) and PAH spectral features are now being used as new probes of the ISM. PAHs are also thought to be among the carriers of the diffuse interstellar absorption bands (DIBs). In the model dealing with the interstellar spectral features, PAHs are present as a mixture of radicals, ions and neutral species. PAH ionization states reflect the ionization balance of the medium while PAH size, composition, and structure reflect the energetic and chemical history of the medium. A major challenge for laboratory astrophysics is to reproduce (in a realistic way) the physical conditions that exist in the emission and/or absorption interstellar zones. An extensive laboratory program has been developed at NASA Ames to assess the physical and chemical properties of PAHs in such environments and to describe how they influence the radiation and energy balance in space and the interstellar chemistry. In particular, laboratory experiments provide measurements of the spectral characteristics of interstellar PAH analogs from the ultraviolet and visible range to the infrared range for comparison with astronomical data. This paper will focus on the recent progress made in the laboratory to measure the direct absorption spectra of neutral and ionized PAHs in the gas phase in the near-UV and visible range in astrophysically relevant environments. These measurements provide data on PAHs and nanometer-sized particles that can now be directly compared to astronomical observations. The harsh physical conditions of the IS medium - characterized by a low temperature, an absence of collisions and strong VUV radiation fields - are simulated in the laboratory by associating a molecular beam with an ionizing discharge to generate a cold plasma expansion. PAH ions are formed from the neutral precursors in an isolated environment at low temperature (of the order of 100 K). The spectra of neutral and ionized PAHs are measured using the high sensitivity methods of cavity ring down spectroscopy (CRDS). These experiments provide unique information on the spectra of free, cold large carbon molecules and ions in the gas phase. Intrinsic band profiles and band positions of cold gas-phase PAHs can now be measured with high-sensitivity spectroscopy and directly compared to the astronomical data. The electronic bands measured for ionized PAH are found to be intrinsically broad (about 20 cm(sup -1)) while the bands associated with the neutral precursors are narrower (of the order of 2 - 10 cm(sup -1)). The laboratory data are discussed and compared with recent astronomical spectra of large and narrow DIBs and with the spectra of circumstellar environments of selected carbon stars and the implications for the interstellar PAH population are derived. Preliminary results also show that carbon nanoparticles are formed during the short residence time of the precursors in the plasma.

Salama, Farid↗

Tribal Revegetation Project Final Project Report: 92-Acre Area, Area 5 Radioactive Waste Management Complex, Nevada National Security Site, Nevada

Nuwu (Southern Paiute), Newe (Western Shoshone), and Nuumu (Owens Valley Paiute) are linguistically related, Numic-speaking peoples who are part of the broader Uto-Aztecan language group. Numic peoples view the land as a holistic, living, sentient being with feelings and purpose. The land is personified with human characteristics and it needs to be experienced to be understood through “learning by doing.” Numic peoples do not support ground disturbing activities within their ancestral lands, including activities tied to the storage of low-level radioactive waste or classified materials on the NNSS, which they view as culturally inappropriate. These deep-rooted ancestral connections are the impetus for reinforcing Numic responsibility for healing disturbed areas by integrating respect and patience with consistent Tribal interaction. Tribal Ecological Knowledge (TEK) is the science of describing Tribal approaches for understanding natural resources. Numic TEK is embedded in traditional teachings learned incrementally over time though experience and it evolves through lessons learned and responses to environmental changes over millennia. Therefore, TEK can broaden and enhance Western scientific knowledge associated with revegetation, especially in highly disturbed areas. The project blended TEK with Western scientific ecological methods to create a vegetative cover within test plots on the 92-Acre Area located at the Radioactive Waste Management Complex (RWMC) located in Area 5 on the Nevada National Security Site (NNSS). The vegetated test plots were systematically created for the Department of Energy (DOE) in tandem with the existing Federal Facilities Agreement and Consent Order (FFACO) with the Nevada Department of Environmental Protection (NDEP). Three previous contractor-lead attempts at revegetation, one targeting full cover revegetation and two targeting test plot revegetation, did not achieve the anticipated results at this location. When presented to the 16 American Indian Tribal nations and affiliated groups with cultural and historical ties to the NNSS, the group appointed a Tribal Revegetation Committee (TRC) that included six expert Tribal knowledge holders to collaborate with an ethnoecologist/cultural anthropologist and two biologists. Project design, planning, seed and outplant selection, spiritual land preparation, and methodology were guided by the TRC and an ethnoecologist/cultural anthropologist from Portland State University (PSU) and biologists from Desert Research Institute (DRI). Using TEK, the TRC recommended a specific seed mixture that contained nine native plant species and three species of outplants. The revegetation effort included preparing and planting thirty 10 m × 10 m (32.8 ft × 32.8 ft) seeded plots, twelve of which also included outplants; and eight 10 m × 100 m (32.8 ft × 328 ft) plots that only received outplants, all atop a waste cell cover. The TRC and the project team creatively adapted TEK with Western scientific methods so that the planned revegetation efforts could occur within the safety and security parameters of the RWMC. Test plots were subjected to one of five soil treatments with varying combinations of straw or mulch applications, soil amendment, and/or outplant planting and one of two watering treatments (watered or unwatered). Planting was divided into two events: one in the fall season and another during the subsequent spring season based on TEK and a corresponding recommendation from the TRC. Monitoring and spiritual management was conducted by the TRC to evaluate plant progress on a monthly basis (in conjunction with the biologist and anthropologist) during each respective growing season for a period of three years after planting. This approach allowed Tribal members the opportunity to conduct traditional blessings and other culturally appropriate activities to restore balance to the land in accordance with Tribal protocols. Following TEK-guided methods, successful plant establishment from seed stock and outplants was observed in many plots. Overall, plots planted in the spring, as recommended by the TRC, showed higher rates of outplant survival and native seedling emergence than those planted in the fall. This finding is significant because it is contrary to the original guidance and previous approaches provided for planting in this region. The TRC believes the frequent co-occurrence of native seedlings near surviving outplants indicates an important symbiotic relationship understood by Tribal communities and overlooked by others. Watered outplants displayed much higher survivability than unwatered plants, even after watering was reduced after the plants were established. Soil amendments and mulch created higher densities of native plants from seed. Many native seedlings showed significant delays in germination, which is considered a normal adaptation to desert climates. Some native plants did not germinate until the third growing year, whereas others germinated during the first growing year, which demonstrates the complexity of the desert environment. Evidence of native insects, reptiles, mammals, and birds, as well as native plants that were not part of the planted species, were noted and considered culturally significant. Despite the presence of non-native plants, native outplants continued to thrive and the incidence of native plant germination from seed increased over time. These successful revegetation results where previous efforts were unsuccessful reinforce the importance of integrating regionally appropriate, TEK-guided methodology with Western science to achieve positive results and the necessity of integrating Tribal involvement in all stages of the revegetation effort. Expanded approaches coupled with Tribal knowledge and tools from Western science addressed a complex problem tied to revegetating atop a low-level radioactive waste cell. The level of Tribal participation serves as a progressive model for building collaborative relationships and addressing ecological challenges on the NNSS.

54 ENVIRONMENTAL SCIENCES↗

An Evaluation of Actinide Reactivity with CO 2 , O 2 , and O 2 /He Gases using Inductively Coupled Plasma Tandem Mass Spectrometry: Application to Simultaneous Measurement of 241 Am/ 241 Pu Ratios in Unseparated Complex Matrices

Accurate actinide measurements are critical within the field of nuclear science. Traditional methods for actinide quantification require time-consuming sample processing prior to analysis. There is a need for rapid analytical techniques that still maintain a high degree of accuracy. In this work, actinide reactivity was assessed for multiple oxygen-containing reaction gases using quadrupole inductively coupled plasma tandem mass spectrometry (Q-ICP-MS/MS) to evaluate actinide analysis in complex sample matrices without analyte-matrix separation. A novel method was developed to measure 241 Am/ 241 Pu in complex sample matrices using O 2 /He reaction gas with no matrix removal or analyte pre-concentration. This inline method reduces matrix-derived polyatomic interferences that complicate traditional ICP-MS analyses by mass-shifting to 241 Am 16 O + and 241 Pu 16 O 2 + , allowing Am and Pu to be mass separated for simultaneous analysis. While mass shifting is efficient, a small portion of Am + (<1.3%) and Pu + (<1.4%) react to from AmO 2 + and PuO + , respectively. Therefore, a mass balance approach was used, in combination with reactivity determined from 242 Pu and 243 Am standard solutions, to correct for residual 241 PuO + and 241 AmO 2 + . The method was validated by measuring 241 Am/ 241 Pu in Pu isotope standards CRM-136 and CRM-137 (separated in March/April 1970 and February 2022, respectively) in both neat solutions and complex matrices containing diluted soil (NIST SRM 2711a, >1000 µg·g -1 ). Method detection limits of 15.9 and 9.6 fg·g -1 were determined for 241 Am and 241 Pu, respectively, and 241 Am/ 241 Pu ratios were measured with accuracies within <3.5%. In conclusion, this work presents the first direct analysis of 241 Am/ 241 Pu in unseparated complex matrices, advancing capabilities for rapid actinide measurements.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

RELAP-7 Application and Enhancement for FLEX Strategies and ATF Behavior under Extended Loss of AC Power Conditions

This report summarizes the results of a three-year research project sponsored by the U.S. Department of Energy (DOE) Nuclear Energy University Program (NEUP) to enhance and apply the RELAP-7 code by adding and improving several important components (e.g., a mechanistic Reactor Core Isolation Cooling (RCIC) system model) for thermal hydraulic studies of LWRs under ELAP conditions and evaluating the time available for transition to portable FLEX equipment. The project team included University of Massachusetts–Lowell, The Ohio State University, Texas A&M University, Idaho National Laboratory and Oak Ridge National Laboratory. In the Fukushima accident, it was found that the RCIC system played a crucial role in delaying core meltdown by almost three days in Fukushima Daiichi Unit 2, because of self-regulated operation of the steam driven RCIC turbine-pump injection system. Steam flow in the convergent-divergent nozzles of the RCIC Terry turbine is two-phase non-equilibrium transonic flow with homogenous nucleation condensation. To more accurately predict the dynamic process and behavior of the transonic compressible steam flow, a one-dimensional transient two-phase analytical model is presented. A simplified four-fluid model was employed in the present work with the consideration of four separate fluid fields: vapor, liquid film, entrained droplets and condensed droplets. The mass, momentum and energy interactions between the fluids were considered and modeled. An extended seven-equation non-equilibrium critical flow model was developed to obtain the critical pressure and velocities of each phase at the nozzle throat. To predict the wetness in the divergent section, a mechanistic nucleation condensation model was integrated in the nozzle analysis model, considering the generation and consequent growth of droplets. The governing differential equations on a staggered grid were discretized using the second-order Lax-Wendroff scheme with a flux limiter, and the Semi-Implicit Method for Pressure-Linked Equation (SIMPLE) algorithm was employed to solve the discrete linear system. To demonstrate the predictability and reliability of the physical models and the numerical method proposed in the present work, three representative nozzles were modeled and simulated. The results show good agreement with the available experimental data, even for condensation shock. Then, the 1D nozzle model was employed to obtain nozzle flow tables of the Terry turbine nozzle for different working pressures which can cover the operation pressure range of the RCIC system. A mechanistic RCIC turbine-pump system model was developed and implemented in the system code TRACE to simulate dynamic responses of the RCIC system under Beyond Design Basis Accident (BDBA) conditions. The turbine-pump governing equations are based on the control volume approach of the angular momentum balance. The physics based mechanistic RCIC model was developed using the TRACE control system components (i.e., signal variables, control blocks, and tables), and incorporated into a TRACE boiling water reactor (BWR) model. The TRACE model in this report has a detailed nodalization of the reactor pressure vessel (RPV), and all of the major flow paths and system components, including the safety relief valves (SRVs) and the containment suppression pool and drywell. Based on the nozzle flow tables generated from the 1D nozzle model developed, the turbine drive torque can be calculated from table lookup. Since the detailed specifications of the RCIC pump are unavailable, the homologous curves for a Bingham pump were used in the current pump component. A station black-out (SBO) accident test problem was selected to demonstrate the TRACE RCIC model. The short-term SBO simulations were performed for two cladding materials: Zircaloy and FeCrAl, to demonstrate the effect of the accident tolerant fuel cladding on fuel heat-up under BDBA accident conditions. The wetwell plays a vital safety role in SBO and other BWR accident scenarios in that it can reduce containment pressure and supply additional core make-up water. The suppression pool temperature distribution has a very large impact on both RPV and containment pressure. Thus, another novel contribution of the project comes mainly from an improved, systems-level wetwell model which can capture buoyancy-induced thermal stratification effects due to steam injection and condensation. A two-zone stratified wetwell model has been implemented in RELAP-7 and some results from that model are presented. This wetwell model is capable of simulating thermal stratification due to a low steam mass injection rate. With a low mass flow rate, the model assumes that all the steam condenses within the pipe and the resulting plume can be approximated with a purely buoyant, heat-source driven model. The wetwell model developed with these assumptions is adequate to simulate slow transients such as extended SBO transients.

42 ENGINEERING↗