Search NASA⌕ Search

SEARCH · Search NASA

Results for “model determination”

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 433 records · Page 24

Pyrolysis Modeling of PMMA decomposition studied by TGA

Data from four TGA experiments conducted at Sandia National Laboratories was used for determination of a pyrolysis model using a commercial thermokinetics program developed by Netzsch Instruments (Kinetics NEO, version 2.1). The data measured at 1 K/min and the average of three measurements at 50 K/min were used as input into Kinetics NEO. The model was developed using data in the range 373 to 773 K. An initial estimate of the energy of activation (E) and pre-exponential constant (A) were determined from the model-free Friedman approach.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

CHARACTERIZATION ON ANISOTROPIC THERMAL CONDUCTIVITY FOR BIG AREA ADDITIVE MANUFACTURING WITH POLYMERS

Additive manufacturing with polymers has been used mainly for prototyping. A recent development of Big Area Additive Manufacturing (BAAM) at Oak Ridge National Laboratory has opened its applications in the mold and die industry. A numerical simulation and prediction for a mold heating performance requires accurate anisotropic thermal properties of the printed material, which are challenging to obtain, and often requires the use of multiple techniques. The transient plane source (TPS) technique has been widely used due to its ability to measure the thermal properties of an extensive range of materials (solids, liquids, powder). Despite the capability to characterize thermal conductivity k of isotropic and anisotropic materials, the measurements of latter materials are limited to the cases, where the samples have the same thermal conductivity k along x- and y-axis that form the radial plane. In this work, the method for a characterization of k in all three dimensions is developed, and the application of TPS is extended to the determination of thermal properties along the x-, y-, and z-axis individually. The materials are represented by additively manufactured polymers including polylactic acid (PLA) and styrene maleic anhydride (SMA). The developed method consists of (1) a determination of the heat capacity of the polymers by means of TPS in combination with the developed in this work data analysis procedure, (2) a machining three types of cylindrical samples from the same material, with the height corresponding either to x-, y-, or z-direction of printing, and (3) a determination of axial thermal conductivity employing anisotropic model and using previously determined heat capacity

Trofimov, Artem↗

Empirical evidence that mean soil contaminant concentration is not a conservative indicator of external exposure to wildlife

Adequate information on external exposure is often a limiting factor in ecological risk assessments because of a lack of exposure data for free-ranging animals as they interact with a heterogeneous distribution of contaminants. Instead, most wildlife exposure data are estimated from model simulations. However, empirical exposure data are rarely available to verify model results, or to test fundamental assumptions inherent in exposure assessments. We used novel, GPS-coupled radiation dosimeters on reindeer (Rangifer tarandus) to quantify their external exposure for five months in a Chernobyl fallout affected mountainous area of Norway. The exposure data were then used to test a fundamental assumption in screening-level risk assessments: Mean soil contaminant concentrations can be used to conservatively estimate chronic external exposures to representative animals. External exposure and GPS coordinates were obtained hourly and compared to model simulations of exposure based on airborne surveys of 137 Cs soil contamination. Our objectives were to determine 1) if modelled dose estimates based on mean 137 Cs soil concentrations were conservative, and 2) if modelled dose estimates could be improved by adjusting for spatial–temporal heterogeneity of contaminants using 137 Cs concentrations within reindeers’ GPS-defined home ranges and core areas. We found that external exposure modelled from mean soil contaminant concentrations underestimated exposure by 70% compared to empirical measurements from GPS-dosimeter collars worn on individual animals. Using area-weighted mean contaminant levels within the animals’ home range improved the model estimates, but still underestimated the field dosimeter data by 53%. The best estimate of external exposure was obtained when high-use areas within their home range, along with spatial heterogeneity of the contaminant, was accounted for by using area-weighted mean soil contaminant levels within GPS-derived core areas. Accounting for temporal use of core areas within animals’ home range provided a 24% improvement. Here, a comparison of our analyses to similar work performed on Chernobyl wolves found the two studies concur that risk assessors should be cautious in relying on the paradigm of using averaged soil contaminant levels to provide conservative estimates of external exposure to representative individuals. Results of modelled external exposure based on mean soil contaminant concentrations could misguide assessors and cause inappropriate management decisions.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗

Integrated Hourly Meteorological Database of 20 Meteorological Stations (1981-2022) for Watershed Function SFA Hydrological Modeling

This dataset contains (a) a script “R_met_integrated_for_modeling.R”, and (b) associated input CSV files: 3 CSV files per location to create a 5-variable integrated meteorological dataset file (air temperature, precipitation, wind speed, relative humidity, and solar radiation) for 19 meteorological stations and 1 location within Trail Creek from the modeling team within the East River Community Observatory as part of the Watershed Function Scientific Focus Area (SFA). As meteorological forcings varied across the watershed, a high-frequency database is needed to ensure consistency in the data analysis and modeling. We evaluated several data sources, including gridded meteorological products and field data from meteorological stations. We determined that our modeling efforts required multiple data sources to meet all their needs. As output, this dataset contains (c) a single CSV data file (*_1981-2022.csv) for each location (20 CSV output files total) containing hourly time series data for 1981 to 2022 and (d) five PNG files of time series and density plots for each variable per location (100 PNG files). Detailed location metadata is contained within the Integrated_Met_Database_Locations.csv file for each point location included within this dataset, obtained from Varadharajan et al., 2023 doi:10.15485/1660962. This dataset also includes (e) a file-level metadata (flmd.csv) file that lists each file contained in the dataset with associated metadata and (f) a data dictionary (dd.csv) file that contains column/row headers used throughout the files along with a definition, units, and data type. Review the (g) ReadMe_Integrated_Met_Database.pdf file for additional details on the script, methods, and structure of the dataset.The script integrates Northwest Alliance for Computational Science and Engineering’s PRISM gridded data product, National Oceanic and Atmospheric Administration’s NCEP-NCAR Reanalysis 1 gridded data product (through the `RCNEP` R package, Kemp et al., doi:10.32614/CRAN.package.RNCEP), and analytical-based calculations. Further, this script downscales the input data into hourly frequency, which is necessary for the modeling efforts.

54 ENVIRONMENTAL SCIENCES↗

River Dissolved Oxygen Prediction Using Machine Learning Models and Wireless Sensor Measurements

Simultaneous flooding&heat and droughts&heat events can potentially destabilize hydro-meteorological conditions to deteriorate the water quality of Neches River. Machine learning (ML) models utilizing wireless sensor measurements have been applied to predict water quality and optimize various water management strategies. This study aims to develop ML models to predict dissolved oxygen (DO) prediction under various hydro-meteorological conditions and enhance water management decision-making. Wireless sensor measurements of DO, water temperature, sample depth, conductivity, turbidity, and pH, along with discharge from the United States Geological Survey stations, are collected for model inputs at the Pine Island Bayou C749 station (PIB-C749) and Neches River Saltwater Barrier (SWB). Multilayer perceptron neural networks, recurrent neural networks, long short-term memory (LSTM), and bidirectional LSTM (BiLSTM) with and without attention mechanism (AT) are tested to determine the best model, which is applied the rolling forecast method to predict 14-day DO. Traditional and recurrent transfer learning (TL and RTL) methods are adopted to overcome insufficient data at the SWB. The input feature importance analysis using the integrated gradients (IG) algorithm is applied to determine dominant inputs. The results show LSTM-based models are capable handling long sequential data. AT-BiLSTM and RTL-LSTM demonstrate the best performance at the PIB-C749 (RMSE=0.054) and the SWB (RMSE=0.028), respectively. TL and RTL methods significantly improve model performance at the SWB. DO, temperature, and pH show higher importance, consistent with hydrodynamics and water chemistry. Both best models are applied to predict 14-day DO and demonstrate reasonable performance for decision-making. Hydro-meteorological conditions of 2017 flood and 2012 drought events are simulated and reveal that possible hypoxia occurs after flooding due to increasing temperature and turbidity, and DO concentration decreases significantly under heat and drought conditions. In conclusion, LSTM-based models utilizing wireless sensor data can be a timely and effective approach to make appropriate decisions on water resource management.

54 ENVIRONMENTAL SCIENCES↗

Mechanistic understanding of pH effects on the oxygen evolution reaction

The oxygen-evolution reaction (OER) is pivotal in many energy-conversion technologies as it is an important counter reaction to others that convert stable chemicals to higher-value products using electrochemistry. The local microenvironment and pH for the anode OER can vary from acidic to neutral to alkaline depending on the system being explored, making definitive mechanistic insights difficult. In this paper, we couple experiments, first-principles calculations based on density functional theory, microkinetics, and transport modeling to explore the entire pH range of the OER. At low current densities, neutral pH values unexpectedly perform better than the acidic and alkaline conditions, and this trend is reversed at higher current densities (> 20 mA cm -2 ). Using multiscale modeling, this switch is rationalized by a change from a dual-reaction mechanism to a single rate-determining step. The model also shows how the alkaline reaction rates dominate in the middle to high pH range. Furthermore, we explore that the local pH for near-neutral conditions is much different (e.g., 2.4 at the reaction surface vs. 9 in the bulk) than the pH extremes, demonstrating the criticality that transport phenomena plays in kinetic activity.

36 MATERIALS SCIENCE↗

Shell-model study of calcium isotopes toward their drip line

Here, we report in this paper a study in terms of the nuclear shell model about the location of the calcium isotopes drip line. The starting point is considering the realistic two-body potential derived by Entem and Machleidt within chiral perturbation theory at next-to-next-to-next-to-leading order (N 3 LO), as well as a chiral three-body force at next-to-next-to-leading order (N 2 LO) whose structure and low-energy constants are consistent with the two-body potential. Then we construct the effective single-particle energies and residual interaction needed to diagonalize the shell-model Hamiltonian. The calculated two-neutron separation energies agree nicely with experiment until 56 Ca, which is the heaviest isotope whose mass has been measured, and do not show any sign of two-neutron emission until 70 Ca . We discuss the role of the choice of the model space in determining the neutron drip line, and also the dependence of the results on the parameters of the shell-model Hamiltonian.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Affinity propagation clustering of full-field, high-spatial-dimensional measurements for robust output-only modal identification: A proof-of-concept study

Determination of the model order is a challenging problem in system identification, especially in output-only or operational modal identification where some modes are weakly excited. While existing methods such as the stabilization diagram method (spectral information) are effective, they do not scale to high-dimensional data, which is usually needed for high-fidelity characterization of structural dynamics and has been made available in the emerging full-field measurement techniques using optical methods such as photogrammetry and laser vibrometers. In this proof-of-concept study we present a new non-parametric, data-driven approach for robust output-only identification of high-spatial-dimensional modal parameters of basic structures by efficiently processing and interactively exploiting the full-field measurement (i.e., very dense spatial measurement points). Specifically, we first over-estimate the system model once, producing a pool of candidate modes associated with their modal frequencies and full-field, high-spatial-dimensional mode shapes. This is accomplished by a data-driven method termed affinity propagation clustering (APC), where the active clusters, which are the active modes in our formulations, emerge from the "message-passing" procedure and does not require a pre-determination of the cluster number (mode or model order). Next, rather than using the spectral information to distinguish the physical and spurious modes in the stabilization diagram method, we exploit and visualize the spatial, full-field mode shape associated with each candidate mode to do so. We conduct extensive experiments on basic structural models with comparisons to a few existing methods. The results indicate that the new method is computationally efficient for identifying high-spatial-dimensional modal parameters, and robust to identify weak modes by exploiting the full-field measurement. We also discuss its applicability and limitations for structures with complex geometry (shapes).

42 ENGINEERING↗

Uncertainty Quantification Enabled by Automatic Differentiation for Hydrodynamic Simulation of Shock‐to‐Detonation Transition in High Explosives

Quantifying the effects of uncertainty in a reactive burn model on the run-to-detonation time in high explosives (HEs) provides a robust methodology for assessing the probability of an HE failing the IHE qualification standard. Moreover, uncertainty quantification helps evaluate whether the model calibration accurately represents data outside the calibration set. This study uses a specialized hydrodynamic simulation code for modeling detonation to determine the run-to-detonation time of the HE PBX 9502 for various impact velocities. To quickly approximate uncertainties in the model, a surrogate was constructed using a Taylor series expansion centered at the mean of the input parameters. To obtain the sensitivities required for constructing the Taylor series, HYP-percomplex Automatic Differentiation (HYPAD) was implemented. HYPAD is a methodology for infusing existing codes with automatic differentiation capabilities by augmenting variables with one or more imaginary units to compute step-size independent partial derivatives. These derivatives are accurate to machine precision with respect to the implemented numerical algorithm, meaning their accuracy reflects that of the underlying method (e.g., integration or discretization schemes). Using reduced order modeling techniques, the mean and standard deviation of the run-to-detonation time of a shock within PBX 9502 were computed for a number of initial impact velocities. A weighted least squares regression was then performed to obtain a best fit curve and prediction interval for the computed statistics. Historical data points from explosively driven wedge tests were utilized to validate the prediction interval, ensuring its reliability in predicting future outcomes. With this prediction interval and a known safety constraint curve, the most probable point of failure and the probability of failure for the HE PBX 9502 were determined.

97 MATHEMATICS AND COMPUTING↗

Cascadia Subduction Zone Fault Heterogeneities From Newly Detected Small Magnitude Earthquakes

The Cascadia subduction zone (CSZ) is known to host M9 megathrust ruptures; however, no such event has occurred in historical observation. The distribution and characteristics of small- to moderate-sized earthquakes can be used to determine the behavior of the megathrust fault but are notably absent offshore the CSZ due to the distance from onshore seismometers. We use automated subspace detection coupled with an onshore-offshore seismic deployment to find small-magnitude earthquakes in the offshore seismogenic zone and analyze their locations in the context of interseismic locking and seismogenic zone extent. Here we detected and located 5,282 earthquakes, 4,096 of which had been previously undetected. We find that the downdip extent of the seismogenic zone as defined by interplate seismicity agrees with the 20% locking contour of the Schmalzle et al. (2014, https://doi.org/10.1002/2013GC005172) geodetic model and extends deeper than predicted by previous thermal models. We cannot determine the updip extent of the seismogenic zone; this may be due to a lack of templates for detection in the updip source area, stress shadows updip of asperity loading, and/or strong locking to the trench. We present a map of possible asperities determined by the small earthquakes in this study. Our asperity locations and extents show some, but not complete, agreement with the asperities modeled from the 1700 M9 rupture and geodetic locking models, and good agreement with the paleo-rupture extents determined from offshore turbidites and forearc basin-based asperity estimates. This highlights the need of continued offshore observations over time, and to elucidate fine-scale variation in locking.

58 GEOSCIENCES↗

An Open‐Source, Physics‐Based, Tropical Cyclone Downscaling Model With Intensity‐Dependent Steering

Abstract An open‐source, physics‐based tropical cyclone (TC) downscaling model is developed, in order to generate a large climatology of TCs. The model is composed of three primary components: (a) a random seeding process that determines genesis, (b) an intensity‐dependent beta‐advection model that determines the track, and (c) a non‐linear differential equation set that determines the intensification rate. The model is entirely forced by the large‐scale environment. Downscaling ERA5 reanalysis data shows that the model is generally able to reproduce observed TC climatology, such as the global seasonal cycle, genesis locations, track density, and lifetime maximum intensity distributions. Inter‐annual variability in TC count and power‐dissipation is also well captured, on both basin‐wide and global scales. Regional TC hazard estimated by this model is also analyzed using return period maps and curves. In particular, the model is able to reasonably capture the observed return period curves of landfall intensity in various sub‐basins around the globe. The incorporation of an intensity‐dependent steering flow is shown to lead to regionally dependent changes in power dissipation and return periods. Advantages and disadvantages of this model, compared to other downscaling models, are also discussed.

Meteorology & Atmospheric Sciences↗

Photodiode-based machine learning for optimization of laser powder bed fusion parameters in complex geometries

We report the quality of parts produced through laser powder bed fusion additive manufacturing can be irregular, with complex geometries sometimes exhibiting dimensional inaccuracies and defects. For optimal part quality, laser process parameters should be selected carefully prior to printing and adjusted during the print if necessary. This is challenging since approaches to control and optimize the build parameters need to take into account the part geometry, the material, and the complex physics of laser powder bed fusion. This work describes a data-driven approach using experimental diagnostics for the optimization of laser process parameters prior to printing. A training dataset is generated by collecting high speed photodiode signal data while printing simple parts containing key geometry features with various process parameter strategies. Supervised learning approaches are employed to train both a forward model and an inverse model. The forward model takes as inputs track-wise geometry features and laser parameters and outputs the photodiode signal along the scan path. The inverse model takes as inputs the geometry features and photodiode signal and predicts the laser parameters. Given the part geometry and a desired photodiode signal, the inverse model can thus determine the required laser parameters. Two test parts which contain defect-prone features are used to assess the validity of the inverse model. The use of the model leads to improved part quality (higher dimensional accuracy, reduced dross, reduced distortion) for both test geometries.

36 MATERIALS SCIENCE↗

Unexpected Thermomechanical Behavior of Off-Stoichiometry Epoxy/Amine Materials

Recent studies on off-stoichiometric thermosets reveal unique viscoelastic behavior derived from increased free volume and physical interactions between chain ends. To understand structural characteristics arising from cure and its effect on properties, we developed a Monte Carlo model based on step-growth polymerization. Our model accurately predicted structure–property trends for a two-component system of EPON 828 (EPON) and ethylenediamine. A second epoxy monomer, D.E.R. 732 (DER), was investigated to modulate T g . Binary mixtures of EPON and DER in off-stoichiometric, amine-rich formulations resulted in nonlinear evolution of thermomechanical properties with respect to initial formulation stoichiometry. Modifying our model with kinetic parameters allowing for differential epoxide/amine reaction kinetics only partially accounted for trends in T g , suggesting that spatiotemporal contributions-not captured by our model-were significant determinants of material properties compared to polymer architecture for three-component systems. Here these findings underpin the importance of spatial awareness in modeling to inform the development of dynamic thermosets.

36 MATERIALS SCIENCE↗

A high-resolution pseudo-polygon discrete element model for regional sea ice

Here, this work presents a pseudo-polygon discrete element model for high-resolution sea ice simulations. A scale-invariant bonded particle contact model is proposed to model joints between sea ice floes based on the smeared fracture model and a lattice spring beam model, and the Mohr–Coulomb failure criterion is implemented to represent the shearing failure mechanism of sea ice packings under complex loadings. All mechanical parameters of the bond model can be directly determined from laboratory tests. Validations of the proposed model are made by investigations of mechanical response and failure criteria of field sea ice sheets. Compared with the field observations of sea ice from satellite radar and in situ stress sensors, the proposed model is capable of reproducing the typical constitutive behavior and the Coulomb friction envelope of field sea ice. Finally, the proposed discrete element sea ice model is used to study the effect of loading rates on mechanical behavior including failure strength of regional sea ice.

54 ENVIRONMENTAL SCIENCES↗

Development and Evaluation of Ensemble Learning-based Environmental Methane Detection and Intensity Prediction Models

The environmental impacts of global warming driven by methane (CH 4 ) emissions have catalyzed significant research initiatives in developing novel technologies that enable proactive and rapid detection of CH 4 . Several data-driven machine learning (ML) models were tested to determine how well they identified fugitive CH 4 and its related intensity in the affected areas. Various meteorological characteristics, including wind speed, temperature, pressure, relative humidity, water vapor, and heat flux, were included in the simulation. We used the ensemble learning method to determine the best-performing weighted ensemble ML models built upon several weaker lower-layer ML models to (i) detect the presence of CH 4 as a classification problem and (ii) predict the intensity of CH 4 as a regression problem. The classification model performance for CH 4 detection was evaluated using accuracy, F1 score, Matthew’s Correlation Coefficient (MCC), and the area under the receiver operating characteristic curve (AUC ROC), with the top-performing model being 97.2%, 0.972, 0.945 and 0.995, respectively. The R 2 score was used to evaluate the regression model performance for CH 4 intensity prediction, with the R 2 score of the best-performing model being 0.858. The ML models developed in this study for fugitive CH 4 detection and intensity prediction can be used with fixed environmental sensors deployed on the ground or with sensors mounted on unmanned aerial vehicles (UAVs) for mobile detection.

Majumder, Reek↗

Optimization of Geometric Perturbations on a Rod Moving Through a High Explosive Target

After completing a study to ensure the simulation results were converged, several high resolution 3D Smoothed Particle Hydrodynamic (SPH) simulations of copper rods impacting a high explosive (LX14) target were performed. This was then formulated into an optimization problem: I wanted to find the optimum shape and location of a perturbation on the rod that would maximize its erosion after it left the target. The shape of the perturbation was modeled as a 2D Gaussian bump and parameterized by its location along the rod axis (z 0 ) and amplitude (A). The final mass of the coherent part of the rod as it leaves the target was used as a metric to represent the erosion of the rod, and the optimization was formulated to maximize this metric with respect to the aforementioned design variables. Due to the expensive nature of the high-fidelity 3D SPH simulations, a surrogate model needed to be chosen so that many function calls to the optimizer would be feasible. Thus, a strategic full factorial sampling plan was chosen to build a dataset, which consisted of 24 high-fidelity simulations. Two surrogate models, a third order polynomial regression model and a Gaussian Process Model, were analyzed using a 14%/86% test/train holdout technique. The root mean square and R2 score of the test set was used to determine the best model, and the third order polynomial regression model was chosen as the surrogate model. Finally, the Nelder-Mead Simplex and Basin-hopping optimization algorithms were implemented, and it was found that the two algorithms gave slightly different optimum values. Nelder-Mead gave an optimum point of [z* 0 ;A*] = [9:9;0:4] and Basin-Hopping gave an optimum value of x* = [z* 0 ;A*] = [9:2;0:1].

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Spatiotemporal Route to Understanding Metal Halide Perovskitoid Crystallization

Here, a spatiotemporal experimental route is reported for the antisolvent vapor diffusion crystal growth of metal halide perovskitoids. A computational analysis combining automated image capture and diffusion modeling enables the determination of the critical concentrations required for nucleation and crystal growth from a single experiment. Five different solvent systems and ten distinct organic ammonium iodide salts were investigated with lead iodide, from which nine previously unreported compounds were discovered. Automated image capture of the mother liquor and antisolvent vials was used to determine changes in solution meniscus positions and detect the nucleation event location. Matching the observations to a numerical solution of Fick's second law diffusion model enables the calculation of reactant, solvent, and antisolvent concentrations at both the time and position of the first stable nucleation and crystal growth. A machine learning model was trained on the resulting data, and it reveals solvent- and amine-specific crystallization tendencies. Solvent systems that interact more weakly with dissolved lead species promote crystallization, while those with stronger interactions can prevent crystallization through increased solubilities. Organic amines that interact more strongly with inorganic components and exhibit greater rigidity are more likely to be incorporated into crystalline products.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Verification of Reactor Activation Modelling with Physical Characterization Data and the Impact on Long Term Assessments - 20303

Whiteshell Laboratories, located in Manitoba, Canada, provided research facilities for the Canadian nuclear industry since the early 1960's. The facility was centered on Whiteshell Reactor 1 (WR-1), an organically cooled, heavy water moderated nuclear reactor. WR-1 was safely shut down in 1985, defueled and drained, and has been in a safe Storage with Surveillance (Deferred Decommissioning) state until present day. In support of an Environmental Assessment for proposed in-situ disposal of WR-1, detailed characterization information regarding the reactor core components (calandria, calandria tubes, and pressure tubes) was necessary to provide confidence that the long-term impacts of the facility can be assessed adequately. A simplified estimate of the total activation of core components was completed in 1992 using conservative assumptions, expected reactor configuration and computer models built with WIMS-CRNL, ONEDANT and ORIGEN-S. Due to the simplification, and age of the work, there was considerable doubt that the modelled inventory estimate of the core components was sufficient to support long-term performance modeling of the planned in-situ decommissioning. Additional data would be required to provide this confidence. Several options were considered for characterizing the WR-1 core components, but it was decided to collect only a limited number of samples in an effort to validate the 1992 simplified estimate as bounding. The limited number of samples reduced risks to workers, both radiological and non-radiological, and significantly reduced the costs required to confidently define and bound the radiological inventory of activation products in the core. The results of the sampling were compared to the simplified activation estimation. There was good agreement between the modelled and sampled results, with all major contributors to total activity showing similar relative quantities. The samples were generally lower than predicted by the models, with few exceptions. Notable exceptions included Nb-94, which were higher in the sample results than the model. This was determined to be caused by the use of Niobium containing steel alloys in the stainless steel pressure tube that were not accounted for in the models. The comparison of results provided confidence that the original inventory estimates produced through the simplified model are conservative and bounding, and therefore suitable for use in the long-term assessment of in situ disposal of the WR-1 reactor. This paper provides details of the simplified computer model results, their comparison to results of the sampling, and the changes adopted in the assessment of the in-situ decommissioning of WR-1 as a result. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗