Search NASA⌕ Search

SEARCH · Search NASA

Results for “manage simulation input data”

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 217 records · Page 12

Experimental and Modeling of Low Pressure Plasma as a Complementary Process for Thermal Waste Treatment - 19270

Polymeric waste contaminated with radionuclides is a challenging wasteform in its raw state. This waste is characterized by its high volume to weight ratio, together with the hazard of its radiolysis and consequently, the generation of flammable and explosive gases (e.g. H2 and light hydrocarbons). Thus, minimizing waste volumes together with stabilizing the wasteform are of high significance in the treatment of contaminated polymers. Thermal decomposition of polymeric waste results in gas products, mainly hydrocarbon molecules (C2-C6), which cannot be released to the atmosphere for environmental reasons. Here, we study a secondary gas treatment reactor, using plasma formed by inductive radio-frequency (13.56 MHz RF), in order to decompose the hydrocarbon gas products into smaller molecules. As a case study, we compare the plasma treatment of gases from two sources: a mixed gaseous product collected from the pyrolytic decomposition of high-density polyethylene (HDPE), and industrial butane gas. We also examine the effects of systematic addition of Ar, N{sub 2} and Ar+N{sub 2} mixtures to the treated gases. In this work, we utilize optical emission spectroscopy (OES) to obtain in situ data from the region in the center of the inductive RF coil. The main emissions from discharges in the gases or their mixtures occur in the 200-700 nm spectral region. Their relative intensities depend strongly on the plasma excitation mode, N{sub 2} flow rate, and gas feed composition. The addition of N{sub 2} to the input gas enhanced hydrocarbon decomposition as compared to the addition of pure Ar. This can be attributed to decomposition by free-radical reactions rather than by particles collisions processes. We performed preliminary modeling of the plasma behavior from ignition to steady state, which takes less than 0.01 s under these conditions. These simulations will help us to better understand and optimize a variety of decomposition reactions. The results of this study suggest that coupling low pressure plasma to conventional thermal decomposition techniques can provide an effective treatment process for contaminated solid plastic waste. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Exfoliation Propensity of Oxide Scale in Heat Exchangers Used for Supercritical CO2 Power Cycles

Supercritical CO2 (sCO2) Brayton cycle systems offer the possibility of improved efficiency in future fossil energy power generation plants operating at temperatures of 650 C and above. As there are few data on the oxidation/corrosion behavior of structural alloys in sCO2 at these temperatures, modeling to predict the propensity for oxide exfoliation is not well developed, thus hindering materials selection for these novel cycles. The ultimate goal of this effort is to provide needed data on scale exfoliation behavior in sCO2 for confident alloy selection. To date, a model developed by ORNL and EPRI for the exfoliation of oxide scales formed on boiler tubes in high-temperature, high-pressure steam has proven useful for managing exfoliation in conventional steam plants. A major input provided by the model is the ability to predict the likelihood of scale failure and loss based on understanding of the evolution of the oxide morphologies and the conditions that result in susceptibility to exfoliation. This paper describes initial steps taken to extend the existing model for exfoliation of steam-side oxide scales to sCO2 conditions. The main differences between high-temperature, high-pressure steam and sCO2 that impact the model involve (i) significant geometrical differences in the heat exchangers, ranging from standard pressurized tubes seen typically in steam-producing boilers to designs for sCO2 that employ variously-curved thin walls to create shaped flow paths for extended heat transfer area and small channel cross-sections to promote thermal convection and support pressure loads; (ii) changed operating characteristics with sCO2 due to the differences in physical and thermal properties compared to steam; and (iii) possible modification of the scale morphologies, hence properties that influence exfoliation behavior, due to reaction with carbon species from sCO2. The numerical simulations conducted were based on an assumed sCO2 operating schedule and several generic heat exchanger channel shapes and cross-sectional areas. Implications for the evolution of stresses in the oxide scales formed on sCO2 heat exchangers, and ensuing critical oxide thicknesses for exfoliation, were derived and compared with expectations for an equivalent conventional tubular heat exchanger in a steam cycle (for a given alloy).

Sabau, Adrian S↗

Designs for Risk Evaluation and Management

The Designs for Risk Evaluation and Management (DREAM) tool was developed as part of the effort to quantify the risk of geologic storage of carbon dioxide (CO 2 ) under the U.S. Department of Energy's National Risk Assessment Partnership (NRAP). DREAM is an optimization tool created to identify optimal monitoring schemes that minimize the time to first detection of CO 2 leakage from a subsurface storage formation. DREAM acts as a post-processer on user-provided output from subsurface leakage simulations. While DREAM was developed for CO 2 leakage scenarios, it is applicable to any subsurface leakage simulation of the same output format. The DREAM tool is comprised of three main components: (1) a Java wizard used to configure and execute the simulations, (2) a visualization tool to view the domain space and optimization results, and (3) a plotting tool used to analyze the results. A secondary Java application is provided to aid users in converting common American Standard Code for Information Interchange (ASCII) output data to the standard DREAM hierarchical data format (HDF5). DREAM employs a simulated annealing approach that searches the solution space by iteratively mutating potential monitoring schemes built of various configurations of monitoring locations and leak detection parameters. This approach has proven to be orders of magnitude faster than an exhaustive search of the entire solution space. The user's manual illustrates the program graphical user interface (GUI), describes the tool inputs, and includes an example application.

(Ruprecht) Yonkofski, Catherine↗

Effects of Media Nutrient Variation on Microalgae Productivity and Economics During Semi-Continuous Cultivation

The development of large-scale microalgae growth for biofuel production is currently limited by the cost of biomass production. However, new approaches to infrastructure and cultivation practices are bringing the field closer to realization. Macronutrients in the cultivation media contribute significant costs, especially since their concentrations have not been optimized for specific strains and conditions. Environmental photobioreactors (ePBRs) were used to simulate cultivation under outdoor conditions, during which the nitrogen and phosphorus levels in the media were varied. The growth of two potential biofuel production strains, Picochlorum celeri and Tetraselmis striata, with varying nutrient inputs during summer and winter scripts, respectively, was studied. This study demonstrated that nitrogen and phosphorus in f/2 media could be reduced by more than 60% from the standard formulation, while maintaining growth rates in a semi-continuous harvesting approach. Experiments comparing the standard and reduced nutrient input concentrations were also conducted for both species in 820 L outdoor raceway ponds, in Mesa, AZ. P. celeri grown in these ponds in October had a growth rate of 10.6 ± 0.7 g/m 2 /day and 10.6 ± 0.3 g/m 2 /day for the standard and low-nutrient P. celeri ponds, respectively. T. striata grown in April–May had a growth rate of 16.6 ± 1.4 g/m 2 /day for the standard nutrient input ponds and 17.4 ± 1.1 g/m 2 /day for the low-nutrient input ponds, and in October 14.5 ± 0.6 g/m 2 /day for standard nutrient ponds and 14.4 ± 0.6 g/m 2 /day for low-nutrient ponds. These outdoor data therefore confirmed the indoor ePBR data. Techno-economic analysis shows that, if high growth rates can be attained at lower nutrient concentrations, a reduction of at least 60% in nutrient costs can be achieved. Such results highlight the importance of managing macronutrient media inputs, as these have a considerable contribution to biomass production costs in large-scale facilities. The analysis also points to the importance of maintaining high spent medium recycling rates in an industrial deployment, so as to minimize the losses of nitrogen and phosphorus compounds.

09 BIOMASS FUELS↗

Scaling Relationships Based on Scaled Tank Mixing and Transfer Test Results

This report documents the statistical analyses performed (by Pacific Northwest National Laboratory for Washington River Protection Solutions) on data from 26 tests conducted using two scaled tanks (43 and 120 inches) in the Small Scale Mixing Demonstration platform. The 26 tests varied several test parameters, including mixer-jet nozzle velocity, base simulant, supernatant viscosity, and capture velocity. For each test, samples were taken pre-transfer and during five batch transfers. The samples were analyzed for the concentrations (lbs/gal slurry) of four primary components in the base simulants (gibbsite, stainless steel, sand, and ZrO2). The statistical analyses including modeling the component concentrations as functions of test parameters using stepwise regression with two different model forms. The resulting models were used in an equivalent performance approach to calculate values of scaling exponents (for a simple geometric scaling relationship) as functions of the parameters in the component concentration models. The resulting models and scaling exponents are displayed in tables and graphically. The sensitivities of component concentrations and scaling exponents to the test parameters are presented graphically. These results will serve as inputs to subsequent work by other researchers to develop scaling relationships that are applicable to full-scale tanks.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Technology Development for Dry Storage of Aluminum-Clad Spent Nuclear Fuel - 20490

A candidate disposition pathway for the > 13 MTHM of aluminum-clad spent nuclear fuel (ASNF), owned and managed by the U.S. Department of Energy, is the drying and placement of the SNF into sealed-canister dry storage, with the ASNF-in- canisters 'road-ready' for transportation to and final direct disposal in a repository waste package. Technical information gaps in fuel drying, and fuel dry storage behavior, have previously challenged the declaration of technology readiness for drying and placement of this fuel into the DOE Standard Canister design for > 50 years of safe dry storage. The principal technical information gaps included: i) characterization and thermal dehydration behavior of aluminum (oxy)hydroxide films attendant on the cladding due to film formation during reactor operation and during post-discharge up to long-term wet storage histories; and ii) G-value data to enable estimation of the radiolytic gas generation from the cladding with its (oxy)hydroxide films. Thus, the oxide films on the ASNF challenged the safety of a sealed storage canister with thermal and radiolytic decomposition of the waters on the films that can lead to corrosion, pressurization, and flammability issues. These gaps in the technical information base have largely been closed. This paper discusses the investigations at the Idaho National Laboratory (INL) and the Savannah River National laboratory (SRNL), and outlines the pending technology development work for input to an engineering design to enable a road-ready dry storage system for ASNF. The ASNF inventory considered for road-ready dry storage is stored at the Savannah River Site (SRS) and at the Idaho Nuclear Technology and Engineering Center (INTEC) at the INL. The ASNF inventory in the SRS L Basin is from foreign and domestic research reactors (FRR and DRR), and is diverse in terms of design, irradiation, and post-reactor-discharge storage conditions; these factors yield a range of characteristics of cladding oxide films on ASNF. Mixed aluminum (oxy)hydroxide (boehmite and bayerite/gibbsite) films, non-uniform in thicknesses up to a maximum local thickness bounded by 25 μm, were observed on ASNF materials removed from wet storage in the L Basin and in non-sealed dry storage at the SRS following reactor service and a long-term (up to 40+ years) interim storage history. The ASNF inventory at INTEC, in both wet and vented dry storage, is predominantly from the Advanced Test Reactor (ATR), but it also includes DRR and FRR fuel. To address a profound behavior of these films, radiolysis testing of aluminum specimens with mixed type boehmite/bayerite oxide films was performed to develop basic data on G-values for production of radiolytic hydrogen under dry storage conditions with nominal relative humidity, temperature, and cover gas. Modeling and simulation of canister internal environments with postulated inventories of oxide films provides estimation of the evolution of the conditions of the canister loaded with ASNF. Simultaneous Thermal Analysis (TGA/DSC) of hydrated oxide powders, and laboratory-scale aluminum specimens with a bayerite film (∼10 μm) using TGA methods, inform drying time/temperature conditions to be used for the ASNF. A demonstration project is recommended for Verification and Validation of the drying and storage of the ASNF. Remaining major tasks leading up to the hot demonstration include scale-up radiolysis testing and scale-up drying testing. Engineering design with the information from the technology program will establish the safety basis and enable long-term (> 50 years) dry storage compatibility with ASNF in the DOE Standard Canister pending its transportation to and disposal in a repository. This full capability would show the ASNF-in-canister storage system to be road-ready. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Influence of Agricultural Managed Aquifer Recharge (AgMAR) and Stratigraphic Heterogeneities on Nitrate Reduction in the Deep Subsurface

This data package contains input files for TOUGHREACT for a modeling study examining the effects of managed aquifer recharge on agricultural lands on nitrate cycling and transport in the Central Valley of CA near Modesto. The files contain all the geochemical species, reactions, and hydrological parameters for the model. The files are text files used for the TOUGH family of code created by LBNL. To use the files a license is required. https://tough.lbl.gov/licensing-download/toughreact-licensing-download/Accompanying Paper Abstract: Agricultural managed aquifer recharge (AgMAR) is a proposed management strategy whereby surface water flows are used to intentionally flood croplands with the purpose of recharging underlying aquifers. However, legacy nitrate (NO3-) contamination in agriculturally-intensive regions poses a threat to groundwater resources under AgMAR. To address these concerns, we use a reactive transport modeling framework to better understand the effects of AgMAR management strategies (i.e., by varying the frequency, duration between flooding events, and amount of water) on N leaching to groundwater under different stratigraphic configurations and antecedent moisture conditions. In particular, we examine the potential of denitrification and nitrogen retention in deep vadose zone sediments (~15 m) using variable AgMAR application rates on two-dimensional representations of differently textured soils, soils with discontinuous bands/channels, and soils with preferential flow paths characteristic of typical agricultural field sites. Our results indicate that finer textured sediments, such as silt loams, alone or embedded within high flow regions, are important reducing zones providing conditions needed for denitrification. Simulation results further suggest that applying water all-at-once rather than in increments for a fixed volume of recharge transports higher concentrations of NO3- deeper into the profile, which has the potential to exacerbate groundwater quality. This transport into deeper depths can be aggravated by wetter antecedent soil moisture conditions. However, applying water all-at-once also increases denitrification within the vadose zone by promoting anoxic conditions. We conclude that AgMAR management strategies can be designed to enhance denitrification in the subsurface and reduce N leaching to groundwater, while specifically accounting for lithologic heterogeneity, antecedent soil moisture conditions, and depth to the water table. Our findings are potentially relevant to other systems that experience flooding inundation such as riparian corridors, floodplains, wetlands, and other managed landscapes like dedicated recharge basins.

54 ENVIRONMENTAL SCIENCES↗

Data-model files associated with the manuscript "The Effects of Spatial and Temporal Resolution of Gridded Meteorological Forcing on Watershed Hydrological Responses" (Shuai et al., 2022 HESS)

This data package contains the model inputs and outputs used in "The Effects of Spatial and Temporal Resolution of Gridded Meteorological Forcing on Watershed Hydrological Responses" (Shuai et al., 2022 HESS). The data.zip file contains the data used to drive the model simulations. The model.zip file contains the XML input file for ATS. The notebook.zip file contains the Jupyter notebooks for pre- and post- processing model results. The figures.zip file contains the raw figures associated with the manuscript.Meteorological forcing plays a critical role in accurately simulating the watershed hydrological cycle. With the advancement of high-performance computing and the development of integrated watershed models, simulating the watershed hydrological cycle at high temporal (hourly to daily) and spatial resolution (10s of meters) has become efficient and computationally affordable. These hyperresolution watershed models require high resolution of meteorological forcing as model input to ensure the fidelity and accuracy of simulated responses. In this study, we utilized the Advanced Terrestrial Simulator (ATS), an integrated watershed model, to simulate surface and subsurface flow and land surface processes using unstructured meshes at the Coal Creek Watershed near Crested Butte (Colorado). We compared simulated watershed hydrologic responses including streamflow, and distributed variables such as evapotranspiration, snow water equivalent (SWE), and groundwater table driven by three publicly available, gridded meteorological forcings (GMFs) -- Daily Surface Weather and Climatological Summaries (Daymet), Parameter-elevation Regressions on Independent Slopes Model (PRISM), and North American Land Data Assimilation System (NLDAS). By comparing various spatial resolutions (ranging from 400 m to 4 km) of PRISM, the simulated streamflow only becomes marginally worse when spatial resolution of meteorological forcing is coarsened to 4 km (or 30% of the watershed area). However, the 4 km resolution has much worse performance than finer resolution in spatially distributed variables such as SWE. Using temporally disaggregated PRISM, we compared models forced by different temporal resolutions (hourly to daily), sub-daily resolution preserves the dynamic watershed responses (e.g., diurnal fluctuation of streamflow) that are absent in results forced by daily resolution. Conversely, the simulated streamflow shows better performance using daily resolution compared to that using sub-daily resolution. Our findings suggest that the choice of GMF and its spatiotemporal resolution depends on the quantity of interest and its spatial and temporal scale, which may have important implications on model calibration and watershed management decisions.

54 ENVIRONMENTAL SCIENCES↗

Estimates of Acetonitrile Generation from Scale Melter Testing of LAW Simulants (VSL-19S4573-1, Rev A)

With high nitrate feeds, such as Hanford low activity waste (LAW) feeds, the addition of reductants is necessary in order to control foaming of the molten glass pool in the melter because such foaming can lead to extensive downtime. Sugar, which was used for this purpose at the West Valley Demonstration Project vitrification facility, has been selected as the baseline reductant for the Hanford Waste Treatment and Immobilization Plant (WTP). The amount of sugar required increases with the amount of nitrates present in the feed and decreases with the amount of waste organics present in the feed, which themselves act as reductants. Excessive additions of reductants can be deleterious, leading to over-reduction of the melt and formation of sulfides and molten metals. Consequently, the oxidants and reductants in the feed must be suitably balanced. The basis for achieving this balance was developed by VSL and Atkins for the vitrification of high-sodium-nitrate feeds at Savannah River's M-Area and has been successfully applied to the processing of a wide variety of simulated WTP feeds over many years and was ultimately incorporated into the WTP process control models. The reaction of sugar with nitrates and nitrites in the cold cap controls the redox state of the underlying melt, which prevents foaming, and also significantly decreases the amount of NOx produced. In addition, however, small amount of organics reaction products are formed. Extensive testing at VSL over the past more than 20 years has shown that acetonitrile is the most prevalent such organic species. That testing has also shown that a significant amount of acetonitrile is captured in the liquid scrubbers in the WTP off-gas treatment system and will therefore be present in the secondary liquid effluents. In particular, acetonitrile is prevalent in the submerged bed scrubber (SBS) and wet electrostatic precipitator (WESP) liquid effluents from WTP LAW off-gas treatment. When these liquids are concentrated in the WTP Effluent Management Facility (EMF) evaporator in the direct feed LAW (DFLAW) flow-sheet, VSL testing has shown that the majority of the acetonitrile partitions to the evaporator condensate. Since the evaporator condensate is directed to the Hanford Effluent Treatment Facility (ETF), this creates a potential issue with the ETF waste acceptance criteria. Consequently, there is a need to assess the available test data in order to develop projections of the likely concentrations of acetonitrile in the streams to ETF. The objective of the present report is to provide such an assessment in order to provide input to flow-sheet model projections.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

COMSOL MULTIPHYSICS MODEL FOR DWPF CANISTER FILLING, REVISION 1

This revision is an extension of the COMSOL Multiphysics model previously developed and documented to simulate the temperatures of the glass during pouring a Defense Waste Processing Facility (DWPF) canister. In that report the COMSOL Multiphysics model used a lumped heat loss term derived from experimental thermocouple data based on a nominal pour rate of 228 lbs./hr. As such, the model developed using the lumped heat loss term had limited application without additional experimental data. Therefore, the COMSOL Multiphysics model was modified to simulate glass pouring and subsequent heat input which, replaced the heat loss term in the initial model. This new model allowed for changes in flow geometry based on pour rate as well as the ability to increase and decrease flow and stop and restart flow to simulate varying process conditions. A revised COMSOL Multiphysics model was developed to predict temperatures of the glass within DWPF canisters during filling and cooldown. The model simulations and experimental data were in good agreement. The largest temperature deviations were {approx} 40 C for the 87 inch thermocouple location at 3000 minutes and during the initial cool down at the 51 inch location occurring at approximately 600 minutes. Additionally, the model described in this report predicts the general temperature trends during filling and cooling as observed experimentally. The revised model incorporates a heat flow region corresponding to the glass pouring down the centerline of the canister. The geometry of this region is dependent on the flow rate of the glass and can therefore be used to see temperature variations for various pour rates. The equations used for this model were developed by comparing simulation output to experimental data from a single pour rate. Use of the model will predict temperature profiles for other pour rates but the accuracy of the simulations is unknown due to only a single flow rate comparison.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

The effects of spatial and temporal resolution of gridded meteorological forcing on watershed hydrological responses

Abstract. Meteorological forcing plays a critical role in accurately simulating the watershed hydrological cycle. With the advancement of high-performance computing and the development of integrated watershed models, simulating the watershed hydrological cycle at high temporal (hourly to daily) and spatial resolution (tens of meters) has become efficient and computationally affordable. These hyperresolution watershed models require high resolution of meteorological forcing as model input to ensure the fidelity and accuracy of simulated responses. In this study, we utilized the Advanced Terrestrial Simulator (ATS), an integrated watershed model, to simulate surface and subsurface flow and land surface processes using unstructured meshes at the Coal Creek Watershed near Crested Butte (Colorado). We compared simulated watershed hydrologic responses including streamflow and distributed variables such as evapotranspiration, snow water equivalent (SWE), and groundwater table driven by three publicly available, gridded meteorological forcings (GMFs) – Daily Surface Weather and Climatological Summaries (Daymet), the Parameter-elevation Regressions on Independent Slopes Model (PRISM), and the North American Land Data Assimilation System (NLDAS). By comparing various spatial resolutions (ranging from 400 m to 4 km) of PRISM, the simulated streamflow only becomes marginally worse when spatial resolution of meteorological forcing is coarsened to 4 km (or 30 % of the watershed area). However, the 4 km-resolution has much worse performance than finer resolution in spatially distributed variables such as SWE. Using the temporally disaggregated PRISM, we compared models forced by different temporal resolutions (hourly to daily), and sub-daily resolution preserves the dynamic watershed responses (e.g., diurnal fluctuation of streamflow) that are absent in results forced by daily resolution. Conversely, the simulated streamflow shows better performance using daily resolution compared to that using sub-daily resolution. Our findings suggest that the choice of GMF and its spatiotemporal resolution depends on the quantity of interest and its spatial and temporal scale, which may have important implications for model calibration and watershed management decisions.

54 ENVIRONMENTAL SCIENCES↗

Multi-machine analysis of termination scenarios with comparison to simulations of controlled shutdown of ITER discharges

To improve our understanding of the dynamics and control of ITER terminations, a study has been carried out on data from existing tokamaks. The aim of this joint analysis is to compare the assumptions for ITER terminations with the present experience basis. The study examined the parameter ranges in which present day devices operated during their terminations, as well as the dynamics of these parameters. The analysis of a database, built using a selected set of experimental termination cases, showed that, the H-mode density decays slower than the plasma current ramp-down. The consequential increase in f<sub>GW</sub> limits the duration of the H-mode phase or result in disruptions. The lower temperatures after the drop out of H-mode will allow the plasma internal inductance to increase. But vertical stability control remains manageable in ITER at high internal inductance when accompanied by a strong elongation reduction. This will result in ITER terminations remaining longer at low q (q<sub>95</sub>~3) than most present-day devices during the current ramp-down. A fast power ramp-down leads to a larger change in β<sub>p</sub> at the H-L transition, but the experimental data showed that these are manageable for the ITER radial position control. The analysis of JET data shows that radiation and impurity levels significantly alter the H-L transition dynamics. Self-consistent calculations of the impurity content and resulting radiation should be taken into account when modelling ITER termination scenarios. Here, the results from this analysis can be used to better prescribe the inputs for the detailed modelling and preparation of ITER termination scenarios.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Assessment of Methane Generation from the Thermolysis of Organics in Caustic Waste

Methane generation has been detected during thermolytic Hydrogen Generation Rate (HGR) testing of SRS tank waste and in simulant samples. Simulant testing has shown methane generation in the presence of trimethylsilanol (TMS) and radioactive sample testing has shown methane generation in presence of methylmercury, both of which are found within SRS waste tanks. The 2019 report Evaluation of the Current State of Knowledge for Thermolysis of Organics within SRS Waste Forming Volatile Organic Compounds (VOCs) reviewed historical data and assessed literature pertaining to thermolytic methane generation. That review recommended undertaking thermolysis studies on simulated waste to detect and quantify VOCs. The purpose of this assessment is to review chemistry and current knowledge of organics in SRS waste to recommend which organics in SRS waste tanks, if any, should be tested for thermolytic methane generation. An evaluation of observed methane generation rates has been performed to assess the implication for current CSTF controls. Based on the inputs used in this report, it was concluded that no methane generation data observed to date has challenged the existing 5% organic contribution to the Composite Lower Flammability Limit (CLFL) assumption in the CSTF tank farm at temperatures below 100 °C. Evidence shows that thermolytic pathways to methane generation become significant at temperatures >100 °C in caustic waste. Currently, the influences of thermolytic methane production on the flammability of vessels at >100 °C (e.g., evaporators) is unclear. It is recommended that testing should be undertaken using the HGR experimental approach for methane generation. The initial simulant testing should obtain quantitative methane generation rates at elevated temperatures for methylated siloxanes/silanes and methylmercury with the goal of determining the current mechanism for thermolytic methane generation seen in SRS tank waste. The results of this testing should then be compiled and analyzed to ascertain the risk of challenging the current projected organic contribution of 5% to the CLFL at elevated temperatures. Following this testing, an evaluation of possible flammable gas generation in the 2H and 3H Evaporators should be performed to assess the risk of methane generation at the high temperatures employed by the evaporators.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

O’Hare Airport Short-Term Ground Transportation Modal Demand Forecast Using Gaussian Processes

Here, the principal objective of this study is to analyze the spatial and temporal variation of ground transportation airport demand and provide demand forecast to inform planning capability and explore alternatives for investments to accommodate airport growth. Because of its good adaptability and strong generalization ability for dealing with high-dimensional input, small-sample, and nonlinear spatial data, Gaussian process (GP) regression is used to provide forecast estimates using data from transportation network company (TNC) trips and urban rail passengers at Chicago's O'Hare International Airport. TNC airport trips differ significantly, with three times more distance, more than twice the travel time, and half of the share requests compared with nonairport trips. This highlights the need for separate demand models. Hourly analysis of the rail service indicates that this is likely heavily used by airport workers, whereas TNC services focus on travelers because of variations in the peak demand hours. Heteroscedastic GP regression is implemented because of differences in trip variance between night and day hours. Estimates are given for weekdays and weekend trips, and the 95% confidence intervals are calculated. The introduction of flight schedule information into the models shows marginal improvements in their performance. However, fitting a GP regression becomes computationally expensive with increased sample size and the introduction of spatial components. Transportation planners and policymakers can use the results and methods implemented in this study to optimize transportation assets and provide long-range simulations of the current and future conditions in the area.

42 ENGINEERING↗

Numerical and Experimental Investigations of Polyurethane Foam for Use as Cask Impact Limiter in Accidental Drop Scenarios - 12099

Rigid, closed-cell polyurethane foams are frequently used as cask impact limiters in nuclear materials and hazardous waste transport due to their high energy-absorption potential. When assessing the cask integrity in accidental scenarios based on numerical simulations, a description of the foam damping properties is required for different strain rates and for a wide temperature range with respect to waste heat generation in conjunction with critical operating and environmental conditions. Implementation and adaption of a respective finite element material model strongly relies on an appropriate experimental data base. Even though extensive impact experiments were conducted e.g. in Sandia National Laboratories, Savannah River National Laboratory and by Rolls Royce plc, not all relevant factors were taken into account. Hence, BAM who is in charge of the mechanical evaluation of such packages within the approval procedure in Germany, incorporated systematic test series into a comprehensive research project aimed to develop numerical methods for a couple of damping materials. In a first step, displacement driven compression tests have been performed on confined, cubic specimens at five loading rates ranging from 0.02 mm/s to 3 m/s at temperatures between +90 deg. C and -40 deg. C. Materials include two different polyurethane foam types called FR3718 and FR3730 having densities of 280 kg/m{sup 3} and 488 kg/m{sup 3} from the product line-up of General Plastics Manufacturing Company. Their data was used to adapt an advanced plasticity model allowing for reliably simulating cellular materials under multi-axial compression states. Therefore, an automated parameter identification procedure had been established by combining an artificial neural network with local optimization techniques. Currently, the selected numerical material input values are validated and optimized by means of more complex loading configurations with the prospect of establishing methods applicable to impact limiters under severe accidental conditions. The reference data base is provided by experiments, where weights between 212 kg and 1200 kg have been dropped from heights between 1.25 m and 7 m on confined 10 cm cubic foam specimens. By presenting the deviations between experimental values and the corresponding output of finite element simulations, the potentials and restrictions of the resulting models are highlighted. Systematic compression tests on polyurethane foams had been performed at BAM test site within the framework of a research project on impact limiters for handling casks for radioactive waste. The experimental results had been used to adapt numerical models for simulating the behaviour of different foam types at different temperatures. The loading speed, however, turned out to have a major influence on their flow curves that can not be captured by simple strain-rate dependent multipliers. Especially for guided drop tests that come close to real accidental scenarios there is a significant gap between experimental and numerical results even when applying such advanced material models. Hence, the extensive data base is currently deployed for expanding the standard algorithms to include adequate dynamic hardening factors. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Modeling Approach for the Aluminum-clad Dry Storage Pilot using HFIR Fuel

To confirm that the dry storage of aluminum-clad research reactor spent nuclear fuel (ASNF) will remain within the safety envelope after applied drying schemes and that the resulting evolution of the gas space composition, temperature, and pressure conditions are understood, a dry storage pilot project is being established. The pilot will incorporate an instrumented lid for discrete interval or for on-demand gas composition and temperature monitoring of two DOE Standard Canisters (DSCs) loaded with three High Flux Isotope Reactor (HFIR) inner cores per DSC. Each DSC would be subjected to a separate alternative candidate drying scheme. Canisters will undergo 1 to 5 years of monitoring, including internal temperature and gas sampling to track pressure and composition changes. This report outlines the approach for modeling the ASNF-in-canister behavior in terms of evolving gas space conditions for the ASNF dry storage pilot using HFIR fuel. The ASNF has an adherent surface oxyhydroxide layer comprised of boehmite/bayerite that generates hydrogen when subjected to irradiation. Three-dimensional multi-physics computational fluid dynamics simulations will be executed to compute the thermal field within the DSC and provide inputs to a chemical model employed to compute pressure buildup as hydrogen is generated in the system. Implemented in Cantera, the chemical model solves gas phase and aluminum oxyhydroxide surface-mediated radiolysis reactions. Gas phase reactions are sourced from Wittman and Hanson (2015), whereas surface-mediated reactions are incorporated by fitting experimental data using an optimization algorithm (Abboud, 2023). Water radiolysis reactions from Wren and Ball (2001) are adopted with modifications as described in Abboud (2023c). Understanding the effect of the hydrogen buildup over time is important for long-term storage safety considerations. Modeling results will include the canister pressure, temperature, and composition evolution from the initial helium backfill with the addition of radiolytically-evolved chemical species (e.g., hydrogen and oxygen). The specific HFIR cores for the pilot program have not yet been selected, and the overall design is still in development. The CFD-chemical model used for this work will be based on prior models with necessary updates to allow for improved accuracy and efficiency. The experimental data obtained from the HFIR demonstration will be used to improve and validate the computational models to predict the ASNF-in-canister behavior.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Level-2 Milestone 4468: Lorenz Simulation Interface Beta Release

The goal of Lorenz is to develop a sophisticated web-based framework as a front end to Livermore Computing Center resources. The primary features include: a personalized dashboard for each user, giving an at-a-glace overview of all computing center resources; a front-end to job submission and management; and a software framework to extend the convenience of web interfaces into simulation codes, including support for setting up, launching, steering, and monitoring simulations. Because the breadth of desktop operating systems in use both inside and outside LLNL spans all flavors of Linux, Mac, and Windows, the decision was made to use exclusively web-based technologies in Lorenz. Recent advancements in AJAX and Javascript programming, along with coming advances in the form of HTML5 and CSS3, made the browser the preferred development platform. The goal is to provide the user with a robust application interface that requires nothing more than a modern, standards-based web browser. This approach not only provides convenient access to Lorenz for internal and external users, but also has significant advantages over conventional applications in terms of rapid deployment and updates. The first phase of the Lorenz project is a system “dashboard”, known as MyLC, that gives LC users a rich set of portal applications, or “portlets”, from which to choose. Information that previously required knowledge of dozens of different commands, files, and web pages is now available in one convenient location. This includes account information, bank and disk usage, quotas, cluster utilization, news items, active jobs on LC systems, and more. The extensible software design of the dashboard allows users and staff members to develop their own portlets, and for them to be shared with all users or just within work groups. The second major phase of Lorenz was built on the idea that HPC centers are managed around batch processing, and “jobs” (or “runs”) progress through several discrete stages: setup, submission, queuing/waiting, running, completion, analysis, and archiving. Lorenz aims to ultimately tie all of these pieces together, including integration with the dashboard and application portals. In addition to helping guide users through the process of defining batch parameters, monitoring job progress, and interacting with jobs and job output, Lorenz also aims to promote access to additional tools and capabilities which to this point have been limited to use by power users. Ultimately Lorenz will assist users in choosing the appropriate resources, such as cluster and file system, given the requirements of a job. The third major phase of Lorenz, the subject of this milestone, is the development of an application portal. This involves teaming with application codes to provide a simple graphical interface to the simulations, which in turn seamlessly ties into job submission and monitoring. This will allow users to directly define simulation parameters and launch applications using a simple but powerful web front-end. Because many complex multi-physics codes can have thousands of inputs, providing a general GUI for building an input deck from scratch is impractical. However, most users of large applications start with an existing input file, and make some number of small changes to it. The application portal will feature a process that combines trusted and validated input files with a GUI that guides the user toward making the relatively minor changes needed for their particular problem. The input files will be in the form of templates, containing most of the data needed for a particular problem along with parameterization of changeable inputs. For expert users, there ultimately will be a way to edit the final input file directly - in the event that there are input changes needed that were not captured in the GUI.

97 MATHEMATICS AND COMPUTING↗

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↗