Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Fractional operators”

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 937 records · Page 52

Titanium Optics for Ion Thrusters

Ion thruster total impulse capability is limited, in part, by accelerator grid sputter erosion. A development effort was initiated to identify a material with a lower accelerator grid volumetric sputter erosion rate than molybdenum, but that could utilize the present NSTAR thruster grid design and fabrication techniques to keep development costs low, and perform as well as molybdenum optics. After comparing the sputter erosion rates of several atomic materials to that of molybdenum at accelerator voltages, titanium was found to offer a 45% reduction in volumetric erosion rates. To ensure that screen grid sputter erosion rates are not higher at discharge chamber potentials, titanium and molybdenum sputter erosion rates were measured at these potentials. Preliminary results showed only a slightly higher volumetric erosion rate for titanium, so that screen grid erosion is insignificant. A number of material, thermal, and mechanical properties were also examined to identify any fabrication, launch environment, and thruster operation issues. Several titanium grid sets were successfully fabricated. A titanium grid set was mounted onto an NSTAR 30 cm engineering model ion thruster and tested to determine optics performance. The titanium optics operated successfully over the entire NSTAR power range of 0.5 to 2.3 kW. Differences in impingement-limited perveances and electron backstreaming limits were found to be due to a larger cold gap for the titanium optics. Discharge losses for titanium grids were lower than those for molybdenum, likely due to a slightly larger titanium screen grid open area fraction. Radial distributions of beam current density with titanium optics were very similar to those with molybdenum optics at all power levels. Temporal electron backstreaming limit measurements showed that titanium optics achieved thermal equilibrium faster than molybdenum optics.

Soulas, George C.↗

Fast Image Texture Classification Using Decision Trees

Texture analysis would permit improved autonomous, onboard science data interpretation for adaptive navigation, sampling, and downlink decisions. These analyses would assist with terrain analysis and instrument placement in both macroscopic and microscopic image data products. Unfortunately, most state-of-the-art texture analysis demands computationally expensive convolutions of filters involving many floating-point operations. This makes them infeasible for radiation- hardened computers and spaceflight hardware. A new method approximates traditional texture classification of each image pixel with a fast decision-tree classifier. The classifier uses image features derived from simple filtering operations involving integer arithmetic. The texture analysis method is therefore amenable to implementation on FPGA (field-programmable gate array) hardware. Image features based on the "integral image" transform produce descriptive and efficient texture descriptors. Training the decision tree on a set of training data yields a classification scheme that produces reasonable approximations of optimal "texton" analysis at a fraction of the computational cost. A decision-tree learning algorithm employing the traditional k-means criterion of inter-cluster variance is used to learn tree structure from training data. The result is an efficient and accurate summary of surface morphology in images. This work is an evolutionary advance that unites several previous algorithms (k-means clustering, integral images, decision trees) and applies them to a new problem domain (morphology analysis for autonomous science during remote exploration). Advantages include order-of-magnitude improvements in runtime, feasibility for FPGA hardware, and significant improvements in texture classification accuracy.

Thompson, David R.↗

Liquid Oxygen Thermodynamic Vent System Testing with Helium Pressurization

This report presents the results of several thermodynamic vent system (TVS) tests with liquid oxygen plus a test with liquid nitrogen. In all tests, the liquid was heated above its normal boiling point to 111 K for oxygen and 100 K for nitrogen. The elevated temperature was representative of tank conditions for a candidate lunar lander ascent stage. An initial test series was conducted with saturated oxygen liquid and vapor at 0.6 MPa. The initial series was followed by tests where the test tank was pressurized with gaseous helium to 1.4 to 1.6 MPa. For these tests, the helium mole fraction in the ullage was quite high, about 0.57 to 0.62. TVS behavior is different when helium is present than when helium is absent. The tank pressure becomes the sum of the vapor pressure and the partial pressure of helium. Therefore, tank pressure depends not only on temperature, as is the case for a pure liquid-vapor system, but also on helium density (i.e., the mass of helium divided by the ullage volume). Thus, properly controlling TVS operation is more challenging with helium pressurization than without helium pressurization. When helium was present, the liquid temperature would rise with each successive TVS cycle if tank pressure was kept within a constant control band. Alternatively, if the liquid temperature was maintained within a constant TVS control band, the tank pressure would drop with each TVS cycle. The final test series, which was conducted with liquid nitrogen pressurized with helium, demonstrated simultaneous pressure and temperature control during TVS operation. The simultaneous control was achieved by systematic injection of additional helium during each TVS cycle. Adding helium maintained the helium partial pressure as the liquid volume decreased because of TVS operation. The TVS demonstrations with liquid oxygen pressurized with helium were conducted with three different fluid-mixer configurations-a submerged axial jet mixer, a pair of spray hoops in the tank ullage, and combined use of the axial jet and spray hoops. A submerged liquid pump and compact heat exchanger located inside the test tank were used with all the mixer configurations. The initial series without helium and the final series with liquid nitrogen both used the axial jet mixer. The axial jet configuration successfully demonstrated the ability to control tank pressure; but in the normal-gravity environment, the temperature in the upper tank region (ullage and unwetted wall) was not controlled. The spray hoops and axial jet combination also successfully demonstrated pressure control as well as temperature control of the entire tank and contents. The spray-hoops-only configuration was not expected to be a reliable means of tank mixing because there was no direct means to produce liquid circulation. However, surprisingly good results also were obtained with the sprayhoops- only configuration (i.e., performance metrics such as cycle-averaged vent flowrate were similar to those obtained with the other configurations). A simple thermodynamic model was developed that correctly predicted the TVS behavior (temperature rise or pressure drop per TVS cycle) when helium was present in the ullage. The model predictions were correlated over a range of input parameters. The correlations show that temperature rise or pressure drop per cycle was proportional to both helium mole fraction and tank heat input. The response also depended on the tank fill fraction: the temperature rise or pressure drop (per TVS cycle) increased as the ullage volume decreased.

Thermodynamic vent system↗

Evolution of parton distribution functions in the short-distance factorization scheme

Lattice QCD offers the possibility of computing parton distributions from first principles, although not in the usual $\overline{MS}$ factorization scheme. Calculations are therefore matched to $\overline{MS}$ using a perturbative procedure which is the source of significant uncertainty within the currently accessible kinematics. We present the possibility of computing the z 2 evolution of non-singlet pseudo-parton distribution functions within the short factorization scheme in a numerically improvable way. The goal is to have tools to evolve a calculation to a scale where perturbative uncertainties are less pronounced. We compare a numerical extraction of the evolution operator from lattice data to the computation of z 2 dependence in perturbation theory. Finally, we discuss how this numerical work may be extended to address the two-scale problem that arises when the Ioffe time range must be made large to extend the reach of the calculation of the pseudo-PDF to smaller values of the momentum fraction.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Search for charged-lepton flavour violation in top quark interactions with an up-type quark, a muon, and a $τ$ lepton in proton-proton collisions at $\sqrt{s}$ = 13 TeV

A search for charged-lepton flavour violation (CLFV) in top quark (t) production and decay is presented. The search uses proton-proton collision data corresponding to 138 fb$^{-1}$ collected with the CMS experiment at $\sqrt{s}$ = 13 TeV. The signal consists of the production of a single top quark via a CLFV interaction or top quark pair production followed by a CLFV decay. The analysis selects events containing a hadronically decaying $τ$ lepton and a muon of opposite electric charge, as well as at least three jets, one of which is identified as originating from the fragmentation of a bottom quark. Machine learning classification techniques are used to distinguish signal from standard model background events. The results of this search are consistent with the standard model expectations. The upper limits at 95% confidence level on the branching fraction $\mathcal{B}$ for CLFV top quark decays to a muon, a $τ$ lepton, and an up or a charm quark are set at $\mathcal{B}$(t $\to$ $μτ$u) $\lt$ (0.04, 0.08, and 0.12) $\times$ 10$^{-6}$, and $\mathcal{B}$(t $\to$ $μτ$c) $\lt$ (0.81, 1.71, and 2.05) $\times$ 10$^{-6}$ for scalar, vector, and tensor-like operators, respectively.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Deep Learning Coordinate-Free Quantum Chemistry

Computing quantum chemical properties of small molecules and polymers can provide insights valuable to physicists, chemists and biologists when designing new materials, catalysts, biological probes and drugs. Deep learning can compute quantum chemical properties accurately in a fraction of the time required by commonly used methods such as density functional theory (DFT). However, many of these deep learning architectures require energy minimized molecular geometries as input, which is also computationally expensive, and decreasing the reproducibility and throughput of these methods. In this study, we demonstrate that accurate quantum chemical computations can be performed without optimized geometries by operating in the coordinate-free domain using deep learning on graph encodings. Furthermore, we also find that the choice of graph-encoding architecture substantially affects the performance of these methods. The Wave architecture outperforms graph convolution architectures, particularly on complex molecules. Furthermore, the structures of these graph encoding architectures provide an opportunity to probe an important, outstanding question in quantum mechanics: What types of quantum chemical properties can be represented by local-variable models? We find that Wave, a local-variable model, is more accurately calculates quantum chemical properties. Graph convolutional architectures require global variables, and are not as effective as as Wave. We anticipate that coordinate-free, deep-learning models of quantum chemistry will become valuable tools in chemistry and biology, enabling researchers to rapidly screen chemical databases or identify new molecules using automated, de-novo design algorithms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Topological Data Analysis for Particulate Gels

Soft gels, formed via the self-assembly of particulate materials, exhibit intricate multiscale structures that provide them with flexibility and resilience when subjected to external stresses. Here, this work combines particle simulations and topological data analysis (TDA) to characterize the complex multiscale structure of soft gels. Our TDA analysis focuses on the use of the Euler characteristic, which is an interpretable and computationally scalable topological descriptor that is combined with filtration operations to obtain information on the geometric (local) and topological (global) structure of soft gels. We reduce the topological information obtained with TDA using principal component analysis (PCA) and show that this provides an informative low-dimensional representation of the gel structure. We use the proposed computational framework to investigate the influence of gel preparation (e.g., quench rate, volume fraction) on soft gel structure and to explore dynamic deformations that emerge under oscillatory shear in various response regimes (linear, nonlinear, and flow). Our analysis provides evidence of the existence of hierarchical structures in soft gels, which are not easily identifiable otherwise. Moreover, our analysis reveals direct correlations between topological changes of the gel structure under deformation and mechanical phenomena distinctive of gel materials, such as stiffening and yielding. In summary, we show that TDA facilitates the mathematical representation, quantification, and analysis of soft gel structures, extending traditional network analysis methods to capture both local and global organization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Overview of ASDEX upgrade results in view of ITER and DEMO

Experiments on ASDEX Upgrade (AUG) in 2021 and 2022 have addressed a number of critical issues for ITER and EU DEMO. A major objective of the AUG programme is to shed light on the underlying physics of confinement, stability, and plasma exhaust in order to allow reliable extrapolation of results obtained on present day machines to these reactor-grade devices. Concerning pedestal physics, the mitigation of edge localised modes (ELMs) using resonant magnetic perturbations (RMPs) was found to be consistent with a reduction of the linear peeling-ballooning stability threshold due to the helical deformation of the plasma. Conversely, ELM suppression by RMPs is ascribed to an increased pedestal transport that keeps the plasma away from this boundary. Candidates for this increased transport are locally enhanced turbulence and a locked magnetic island in the pedestal. The enhanced D-alpha (EDA) and quasi-continuous exhaust (QCE) regimes have been established as promising ELM-free scenarios. Here, the pressure gradient at the foot of the H-mode pedestal is reduced by a quasi-coherent mode, consistent with violation of the high-n ballooning mode stability limit there. This is suggestive that the EDA and QCE regimes have a common underlying physics origin. In the area of transport physics, full radius models for both L- and H-modes have been developed. These models predict energy confinement in AUG better than the commonly used global scaling laws, representing a large step towards the goal of predictive capability. A new momentum transport analysis framework has been developed that provides access to the intrinsic torque in the plasma core. In the field of exhaust, the X-Point Radiator (XPR), a cold and dense plasma region on closed flux surfaces close to the X-point, was described by an analytical model that provides an understanding of its formation as well as its stability, i.e., the conditions under which it transitions into a deleterious MARFE with the potential to result in a disruptive termination. With the XPR close to the divertor target, a new detached divertor concept, the compact radiative divertor, was developed. Here, the exhaust power is radiated before reaching the target, allowing close proximity of the X-point to the target. No limitations by the shallow field line angle due to the large flux expansion were observed, and sufficient compression of neutral density was demonstrated. With respect to the pumping of non-recycling impurities, the divertor enrichment was found to mainly depend on the ionisation energy of the impurity under consideration. In the area of MHD physics, analysis of the hot plasma core motion in sawtooth crashes showed good agreement with nonlinear 2-fluid simulations. This indicates that the fast reconnection observed in these events is adequately described including the pressure gradient and the electron inertia in the parallel Ohm’s law. Concerning disruption physics, a shattered pellet injection system was installed in collaboration with the ITER International Organisation. Thanks to the ability to vary the shard size distribution independently of the injection velocity, as well as its impurity admixture, it was possible to tailor the current quench rate, which is an important requirement for future large devices such as ITER. Progress was also made modelling the force reduction of VDEs induced by massive gas injection on AUG. The H-mode density limit was characterised in terms of safe operational space with a newly developed active feedback control method that allowed the stability boundary to be probed several times within a single discharge without inducing a disruptive termination. Regarding integrated operation scenarios, the role of density peaking in the confinement of the ITER baseline scenario (high plasma current) was clarified. The usual energy confinement scaling ITER98(p,y) does not capture this effect, but the more recent H20 scaling does, highlighting again the importance of developing adequate physics based models. Advanced tokamak scenarios, aiming at large non-inductive current fraction due to non-standard profiles of the safety factor in combination with high normalised plasma pressure were studied with a focus on their access conditions. A method to guide the approach of the targeted safety factor profiles was developed, and the conditions for achieving good confinement were clarified. Based on this, two types of advanced scenarios (‘hybrid’ and ‘elevated’

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Evaluation of the Radioactive Material Released in the Harborview Research and Training Building and Some Implications for Emergency Response

On 2 May 2019, during the 137 Cs source recovery operation, a source capsule in a research irradiator containing approximately 77.1 TBq was breached. Based on a geometric reconstruction analysis of the damage to the capsule, approximately 46.3 GBq (0.04%) was impacted by the chop saw (grinder) inside a mobile hot cell on the loading dock at the University of Washington Harborview Research and Training (HRT) Building. A very small fraction of the material impacted, less than 1%, was released from the mobile hot cell and then to the rest of the HRT Building. The objectives of this project were to assess the accidental release of 137 CsCl and its implications related to emergency response methods and the ramifications of 137 CsCl transport. The phenomenology of this event was also compared with past alkali halide dispersal events. The vast number of measurements and samples collected by the remediation contractors, the Department of Energy’s Nuclear Emergency Support Team, and the small number of retrospective samples collected by the authors informed the analysis. The techniques included (1) autoradiography and electron microscopy of samples collected from the HRT Building and the irradiator, (2) 3D visualization of deposition on surfaces and within the ventilation system, and (3) a study of the damage to the source capsule to evaluate the Cs particle size and particle composition due to the grinding accident. Subsequently, the cesium contaminant transport through the numerous pathways in the building was reconstructed to assess the deposition on surfaces as a function of particle size. Furthermore, the implications for emergency response are relevant to data quality and management. A Data Quality Objective guides data collection methods so that they have appropriate accuracy and precision for the intended application. Recommendations were made with respect to the sample collection protocols and archiving of samples.

46 - INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AN↗

Coherent Lidar Turbulence Measurement for Gust Load Alleviation

Atmospheric turbulence adversely affects operation of commercial and military aircraft and is a design constraint. The airplane structure must be designed to survive the loads imposed by turbulence. Reducing these loads allows the airplane structure to be lighter, a substantial advantage for a commercial airplane. Gust alleviation systems based on accelerometers mounted in the airplane can reduce the maximum gust loads by a small fraction. These systems still represent an economic advantage. The ability to reduce the gust load increases tremendously if the turbulent gust can be measured before the airplane encounters it. A lidar system can make measurements of turbulent gusts ahead of the airplane, and the NASA Airborne Coherent Lidar for Advanced In-Flight Measurements (ACLAIM) program is developing such a lidar. The ACLAIM program is intended to develop a prototype lidar system for use in feasibility testing of gust load alleviation systems and other airborne lidar applications, to define applications of lidar with the potential for improving airplane performance, and to determine the feasibility and benefits of these applications. This paper gives an overview of the ACLAIM program, describes the lidar architecture for a gust alleviation system, and describes the prototype ACLAIM lidar system.

Bogue, Rodney K.↗

Algorithm for Estimating the Plume Centerline Temperature and Ceiling Jet Temperature in the Presence of a Hot Upper Layer

The experiments were designed to provide insight into the behavior of jet fuel fires in aircraft hangars and to study the impact of these fires on the design and operation of a variety of fire protection systems. As a result, the test series included small fires designed to investigate the operation of UV/IR detectors and smoke detectors as well as large fires which were used to investigate the operation of ceiling mounted heat detectors and sprinklers. The impact of the presence or absence of draft curtains was also studied in the 15 m hangar. It is shown that in order to predict the plume centerline temperature within experimental uncertainty, the entrainment of the upper layer gas must be modeled. For large fires, the impact of a changing radiation fraction must also be included in the calculation. The dependence of the radial temperature profile of the ceiling jet as a function of layer development is demonstrated and a ceiling jet temperature algorithm which includes the impact of a growing layer is developed.

ALGORITHMS↗

Normal Gravity Testing of a Microchannel Phase Separator for In Situ Resource Utilization

A microchannel separator, with 2.7 millimeters as the smallest dimension, was tested, and a pore throat structure captured and removed liquid from a gas-liquid stream. The microchannel device was tested over a of gas and liquid flow rates ranging from 0.0005 up to 0. 14 volume fraction of liquid. Four liquids were tested with air. The biggest factor affecting the throughput is the capacity of liquid flow through the pore throat, which is dictated by permeability, liquid viscosity, flow area, pore throat thickness, and pressure difference across the pore throat. Typically, complete separation of gas and liquid fractions was lost when the liquid flow rate reached about 40 to 60% of the pore throat capacity. However, this could occur over a range of 10 to 90% utilization of pore throat capacity. Breakthrough occurs in the microchannel phase separator at conditions similar to the annular to plug flow transition of two-phase microgravity pipe flow implying that operating in the proper flow regime is crucial. Analysis indicates that the Bond number did not affect performance, supporting the premise that hydrodynamic, interfacial, and capillary forces are more important than gravity. However, the relative importance of gravity is better discerned through testing under reduced gravity conditions.

TeGrotenhuis, Ward E.↗

Progress in extending high poloidal beta scenarios on DIII-D towards a steady-state fusion reactor and impact of energetic particles

To prepare for steady-state operation of future fusion reactors (e.g. the International Thermonuclear Experimental Reactor and China Fusion Engineering Test Reactor (CFETR)), experiments on DIII-D have extended the high poloidal beta ( β P ) scenario to reactor-relevant edge safety factor q 95 ~ 6.0, while maintaining a large-radius internal transport barrier (ITB) using negative magnetic shear. Excellent energy confinement quality ( H 98y2 > 1.5) is sustained at high normalized beta ( β N ~ 3.5). This high-performance ITB state with Greenwald density fraction near 100% and q min ≥ 3 is achieved with toroidal plasma rotation V tor ~ 0 at ρ ≥ 0.6. This is a key result for reactors expected to have low V tor . At high β P (≥1.9), large Shafranov shift can stabilize turbulence leading to a high confinement state with a low pedestal and an ITB. At lower β P (<1.9), negative magnetic shear in the plasma core contributes to turbulence suppression and can compensate for reduced Shafranov shift to continue to access a large-radius ITB and excellent confinement with low V tor , consistent with the results of gyrofluid transport simulations. These high- β P cases are characterized by weak/no Alfvén eigenmodes (a.e.) and classical fast-ion transport. At high density, the fast-ion deceleration time decreases and Δ β fast is lower; these reduce a.e. drive. The reverse-shear Alfvén eigenmodes are weaker or stable because the negative magnetic shear region is located at higher radius, away from the peaked fast-ion profile. Resistive wall modes can be a limitation at simultaneous high β N , low internal inductance, and low rotation. Analysis suggests that additional off-axis external current drive could provide a more stable path at reduced q 95 . Based on a DIII-D high- β P plasma with large-radius ITB, two scenarios are proposed for CFETR Q = 5 steady-state operation with ~1 GW fusion power: a lower-$l_i$($l_i$ ~ 0.66) and a higher-$l_i$($l_i$ ~ 0.75) case. Using a Landau closure model, multiple energetic particle (EP) effects on the a.e. stability are analyzed modifying the growth rate of the a.e.s triggered by the neutral-beam-injection EPs and alpha particles, although the stabilizing/destabilizing effect is weak for the cases analyzed. The stabilizing effects of the combined EP species β , energy, and density profile in CFETR need further investigation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA Low-Power Stirling Convertor for Small Landers, Probes, and Rovers Operating in Darkness

A variety of mission concepts have been studied by NASA and the U. S. Department of Energy that would utilize low power Radioisotope Power Systems (RPS) for landers, probes, and rovers. These missions would contain science measuring instruments distributed across planetary surfaces or near objects of interest where solar flux is insufficient for using solar cells. Landers could be used to provide data like wind, temperature, pressure, seismic activity, and other planetary measurements to mission planners or future explorers. The studies proposed using fractional versions of the General Purpose Heat Source (GPHS) or multiple Light Weight Radioisotope Heater Units (LWRHU) to heat power conversion technologies for science instruments and communication. Dynamic power systems are capable of higher conversion efficiencies, which could enable equal power using less fuel or more power using equal fuel, when compared to less efficient static power conversion technologies. Providing spacecraft with more power would decrease duty cycling of basic functions and, therefore, increase the quality and abundance of science data. Low power Stirling convertors are being developed at NASA Glenn Research Center (GRC) to provide future micro spacecraft with electrical power by converting heat from one or more LWRHUs. An initial design converts multiple Watts of heat to one Watt of electrical power output using a Stirling convertor. Development of the concept includes maturation of convertor and controller designs, performance evaluation of an evacuated metal foil insulation, and development of system interfaces.

Wilson, Scott D.↗

NASA Low Power Stirling Convertor for Small Landers, Probes, and Rovers Operating in Darkness

A variety of mission concepts have been studied by NASA and the U. S. Department of Energy that would utilize low power Radioisotope Power Systems (RPS) for landers, probes, and rovers. These missions would contain science measuring instruments distributed across planetary surfaces or near objects of interest where solar flux is insufficient for using solar cells. Landers could be used to provide data like wind, temperature, pressure, seismic activity, and other planetary measurements to mission planners or future explorers. The studies proposed using fractional versions of the General Purpose Heat Source (GPHS) or multiple Light Weight Radioisotope Heater Units (LWRHU) to heat power conversion technologies for science instruments and communication. Dynamic power systems are capable of higher conversion efficiencies, which could enable equal power using less fuel or more power using equal fuel, when compared to less efficient static power conversion technologies. Providing spacecraft with more power would decrease duty cycling of basic functions and, therefore, increase the quality and abundance of science data. Low power Stirling convertors are being developed at NASA Glenn Research Center (GRC) to provide future micro spacecraft with electrical power by converting heat from one or more LWRHUs. An initial design converts multiple Watts of heat to one Watt of electrical power output using a Stirling convertor. Development of the concept includes maturation of convertor and controller designs, performance evaluation of an evacuated metal foil insulation, and development of system interfaces.

Wilson, Scott↗

Modification in Applying Appendix D of 40 CFR Part 61 to Heated Solid Radionuclide Materials With High Melting and Boiling Points

Appendix D of Title 40 Part 61 of the US Code of Federal Regulations (CFR) provides a procedure that US Department of Energy (US DOE) facility owners and operators can use to estimate radionuclide emissions to the atmosphere for dose calculations instead of measuring emissions for minor sources under the 40 CFR Part 61, Subpart H, National Emission Standards for Emissions of Radionuclides Other Than Radon From Department of Energy Facilities, regulation. The procedure assumes that any radioactive material heated above 100 °C is completely vaporized and emitted to the atmosphere. In 1991, the US DOE Oak Ridge Reservation (ORR) requested approval to use different release fractions (RFs) for uranium because of its high melting and boiling points. In response to the request, the US Environmental Protection Agency (US EPA) Region IV approved the use of modified RFs for elemental uranium provided no reaction had taken place to alter its chemical form. In 2015, the ORR requested approval to use different RFs for tungsten, again because of its high melting and boiling points. EPA Region IV approved the use of modified RFs for heated radioactive tungsten metal. In accordance with the two precedents set for heating uranium and radioactive tungsten metals, in 2016, the ORR requested approval to use modified RFs in a similar fashion for other radioactive solid metals and compounds with melting and boiling points above 500 °C that might be heated above 100 °C in future research projects and experiments. EPA Region IV again granted approval to use modified RFs for the list of compounds. This note discusses the proposed modified RFs and their development.

36 MATERIALS SCIENCE↗

Technical Gaps in Hanford High-Level Waste Solids Settling Behavior and Settling Time Evaluation for Direct Feed High-Level Waste (DFHLW) Operations

Settling of high-level waste (HLW) solids in process vessels is a key conceptual process step in providing HLW feed directly to the Hanford Waste Treatment and Immobilization Plant (WTP) HLW Vitrification Facility. Direct Feed High-Level Waste (DFHLW) is a potential approach to initiating HLW vitrification prior to completing of the WTP Pretreatment Facility. Settling would be used with subsequent supernatant decant to concentrate HLW feed. To support planning for DFHLW, Washington River Protection Solutions (WRPS) requested support from the Pacific Northwest National Laboratory to evaluate the current data set available to predict the time needed for HLW solids to settle, to identify gaps in the understanding and predictive capability of HLW solids waste settling times, and to provide scoping estimates of the potential settling time. Eight technical gaps were identified for predicting settling times and characteristics of the formed sediment layers including: Gap 1: In-Tank Settling Rates Faster than Settling of Laboratory Samples, Gap 2: Effect of Sludge Leaching/Washing on Predicted Settling Times, Gap 3: Predicting Waste Settling from Waste Chemistry (Waste Type), Gap 4: Predicting Waste Settling from Particle Size and Density Distributions (PSDDs), Gap 5: Insufficient Laboratory and In-Tank Settling Data to Represent Hanford Waste, Gap 6: Methods for Real-Time, In-Tank Tracking of Settling, Gap 7: Prediction of Sediment Erosion Resistance as a Function of Settling Time, and Gap 8: Prediction of Sediment Solids Content as a Function of Settling Time. In addition to the data gaps, an overarching observation of the settling rate and settled layer data is the significant variation in behavior. At similar solids concentrations, settling rates can vary by as much as 3 orders of magnitude depending on the source waste tank, and significantly different settling rates are noted between laboratory and in situ tests for the same waste tank. The range of average solids concentration in existing HLW sediment, which may have been quiescent for decades, can vary from less than 7 wt% to greater than 74 wt% solids. The shear strengths (or yield stresses) measured on laboratory samples range from less than 27 Pa to greater than 6400 Pa. These variations can challenge process planning for the application of a settle/decant process for DFHLW. This report describes the significance of the gaps to the settle/decant process and presents uncertainties by way of examples. Potential technical approaches for resolving these gaps are described and the estimated difficulty in resolving these gaps is evaluated. Based on the significance of the gap and the difficulty of resolution, recommendations are made to address specific gaps. Scoping estimates of the potential settling times for DFHLW solids have been made based on the existing data set with its associated gaps. Depending on the process vessel depth and final sediment concentration, substantial fractions of the scoping estimate results for settling times for characterized HLW exceed the 2-week period that has been previously assumed for process planning. There is also significant disparity, potentially greater than a factor of 5000 difference, in the estimated settling times depending on process vessel depth and final sediment solid concentration. This variation in results underscores the significance of the identified gaps and uncertainties with respect to process planning for utilizing settle/decant operations for DFHLW.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Evaluating Probability of Containment Effectiveness at a GCS Sites using integrated assessment modeling approach with Bayesian decision Networks

Improved scientific and engineering understanding of the behavior of geologic CO2 storage together with established regulatory framework and incentive structures raise the prospects for accelerated, large-scale deployment of this greenhouse gas emissions reduction approach. Incentive structures call for the establishment of appropriate verification and accounting approaches to support claims of the integrity of a geologic storage complex and to justify taking credit for long-term storage. In this study, we present a framework for assessing the probability of containment effectiveness over the lifetime of a geologic carbon storage site (e.g., after 70 years of injection and post-injection site performance) using forward stochastic model realizations based on site characterization data and using a monitoring-informed Bayesian network based on hypothetical detectability from surface seismic surveys over the site injection and post-injection phases. The National Risk Assessment Partnership’s open-source Integrated Assessment Model (NRAP-Open-IAM) was utilized to develop an ensemble of 10,000 a priori stochastic forecasts of CO2 containment. Those simulations were used to train the Bayesian network model to estimate the prior probabilities of the CO2 leakage mass into overlying, monitorable aquifers considering the uncertainties in the reservoir properties, permeability of potentially leaky wells and the overlying aquifers. The conditional probabilities in the Bayesian network were either learned from the NRAP-Open-IAM simulations or derived from the predefined detection thresholds for the monitoring method. Observations obtained from monitoring, over time during the site operation phases were then used to generate updated posterior probabilities of containment (and any loss from containment) in the Bayesian network by propagating the prior probabilities through the conditional probabilities. We demonstrate how to construct and use the Bayesian network for verifying the long-term storage complex effectiveness informed by monitoring based on the NRAP-Open-IAM simulations previously developed for the FutureGen 2.0 site. This approach may have relevance for stake holders to demonstrate secure geologic storage, provide a defensible, probabilistic approach to claim credit for geologic storage, and to estimate the likelihood that any fraction of the claimed credit may need to be refunded to the creditor based on available monitoring information.

Bayesian network, Risk assessment, Monitoring, car↗