Search NASA⌕ Search

SEARCH · Search NASA

Results for “expanded uncertainty”

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

A Component-Level Model of Automatic Dependent Surveillance - Broadcast (ADS-B)

Automatic Dependent Surveillance – Broadcast (ADS-B) is being employed in numerous peer-to-peer initiatives attempting to expand the capacity of the National Airspace System (NAS) or enable mixed operations of manned and unmanned vehicles. Safety assessments of these initiatives rely, in part, on modeling the accuracy of ADS-B in reporting the position and direction of an ownship and surrounding traffic. Frequently, these initiatives utilize a position uncertainty model that applies a reported ADS-B estimation position uncertainty (EPU) value to a Rayleigh distribution and uses a Gauss-Markov random walk to add error to the ADS-B output of a vehicle. This model of ADS-B state error is easy to implement and apply to numerous problems. However, it has a couple of draw-backs. First, the ADS-B state errors are equally probable in all directions. This is a good assumption in situations where aircraft maneuvering is not constrained. However, in situations where the aircraft maneuvering is constrained such as landing, the error distribution is likely to exhibit directionality and the non-directional model may skew results especially when assessing very low probabilities (e.g., 10(exp -9)) of catastrophic encounters. Second, the model does not account for processing latency in the receiving aircraft. NASA Langley Research Center (LaRC) recently examined the feasibility of decreasing the spacing of aircraft on parallel approaches to runways separated by as little as 700 feet. For Monte-Carlo analysis using a high-fidelity simulation of a large transport, LaRC started with a Gauss-Markov model of ADS-B error but then developed a component level model of ADS-B error to increase the fidelity of results.

Madden, Michael M.↗

The Four-Dimensional Carbon Cycle of the Southern Ocean

The Southern Ocean plays a fundamental role in the global carbon cycle, dominating the oceanic uptake of heat and carbon added by anthropogenic activities and modulating atmospheric carbon concentrations in past, present, and future climates. However, the remote and extreme conditions found there make the Southern Ocean perpetually one of the most difficult places on the planet to observe and to model, resulting in significant and persistent uncertainties in our knowledge of the oceanic carbon cycle there. The flow of carbon in the Southern Ocean is traditionally understood using a zonal mean framework, in which the meridional overturning circulation drives the latitudinal variability observed in both air–sea flux and interior ocean carbon concentration. However, recent advances, based largely on expanded observation and modeling capabilities in the region, reveal the importance of processes acting at smaller scales, including basin-scale zonal asymmetries in mixed-layer depth, mesoscale eddies, and high-frequency atmospheric variability. Assessing the current state of knowledge and remaining gaps emphasizes the need to move beyond the zonal mean picture and embrace a four-dimensional understanding of the carbon cycle in the Southern Ocean.

54 ENVIRONMENTAL SCIENCES↗

Silicon and carbon abundances in the Orion nebula

We assess the gas-phase abundances of Si and C from our recent measurements of Si(2+), C(2+), and C(+) in the Orion Nebula by expanding on our earlier 'blister' models. The C and Si abundances are derived from new IUE high-dispersion spectra of the C(2+) 1907, 1909 A and Si(2+) 1883, 1892 A lines and archival IUE data. Gas-phase Si/C = 0.016 in the Orion ionized volume and is particularly insensitive to uncertainties in extinction and temperature structure. The solar value is 0.098. Gas-phase C/H = 2.8 x 10 exp -4 and Si/H = 4.5 x 10 exp -6. Compared to solar, Si is depleted by a factor of about 8 in the ionized region, while C is much less depleted (factor of 1.3), if depleted at all. This suggests that most Si resides in dust grains even in the ionized volume. Thus, most of the observed forbidden Si II 34.8-micron emission in Orion does not arise in the H II region.

Rubin, Robert H.↗

Studies of Error Sources in Geodetic VLBI

Achieving the goal of millimeter uncertainty in three dimensional geodetic positioning on a global scale requires significant improvement in the precision and accuracy of both random and systematic error sources. For this investigation we proposed to study errors due to instrumentation in Very Long Base Interferometry (VLBI) and due to the atmosphere. After the inception of this work we expanded the scope to include assessment of error sources in GPS measurements, especially as they affect the vertical component of site position and the measurement of water vapor in the atmosphere. The atmosphere correction 'improvements described below are of benefit to both GPS and VLBI.

Rogers, A. E. E.↗

An Economics-by-Design Approach Applied to a Heat Pipe Microreactor Concept

Microreactors present a potential paradigm shift in the nuclear industry. Emphasis thus far has been on large-scale multi-billion-dollar projects that cater solely to grid electricity market. These projects can be challenging to finance and execute. On the other hand, microreactors are intended to target a wide variety of smaller niche markets and are expected to be factory-fabricated and more readily deployable. While diseconomies of scale for microreactors may tend to raise their costs per energy output (MWh) relative to large nuclear plants, offsetting gains can be expected from standardization, simplification, passive safety, lower radionuclide inventories, factory fabrication, fast installation, and low financing costs. To adequately assess these contributions, designers should have a different perspective on cost drivers than for large nuclear plants and can utilize novel approaches for systematic cost reduction. To account for these important aspects of microreactors, this report proposes an economics-by-design approach that places economic considerations at the center of the design process. The methodology builds on existing frameworks such as design-to-cost and value engineering, expanding them to new markets (beyond the grid), new attributes (beyond costs alone), and introducing the approach at earlier points in the design cycle. Design parameters and technical specifications are systematically evaluated until costs meet market entry points, while also providing the high-priority performance attributes of the particular use case. Determining first-order estimates for different components early in the process enables designers to focus R&D efforts on the biggest overall cost contributors and components with the most cost uncertainty. The analysis is always guided by market needs and threshold prices. In addition to microreactors, the approach is expected to be useful for other classes of nuclear reactors as well. The analysis was applied to a concept found in the open literature (the Design A heat-pipe reactor). A comprehensive bottom-up estimate was generated by leveraging a new microreactor-specific code of accounts and a range of cost equations. The initial estimate for levelized cost of electricity (LCOE) unsurprisingly exceeded market ranges since the use case had prioritized technological readiness over economic considerations in design choices. An alternate concept was then proposed, with various assumptions/targets made to reduce the largest cost contributors. Changes in the neutron spectrum, the power output, and building structures were found to make even the first-of-a-kind of this modified concept competitive with diesel generation in some remote communities. Learning rate (LR) assumptions indicated cost reductions achieved from sequential unit deployments could expand the range of competitiveness to include additional markets as deployments proceed.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Understanding Crew Decision-Making in the Presence of Complexity: A Flight Simulation Experiment

Crew decision making and response have long been leading causal and contributing factors associated with aircraft accidents. Further, it is anticipated that future aircraft and operational environments will increase exposure to risks related to these factors if proactive steps are not taken to account for ever-increasing complexity. A flight simulation study was designed to collect data to help in understanding how complexity can, or may, be manifest. More specifically, an experimental apparatus was constructed that allowed for manipulation of information complexity and uncertainty, while also manipulating operational complexity and uncertainty. Through these manipulations, and the aid of experienced airline pilots, several issues have been discovered, related most prominently to the influence of information content, quality, and management. Flight crews were immersed in an environment that included new operational complexities suggested for the future air transportation system as well as new technological complexities (e.g. electronic flight bags, expanded data link services, synthetic and enhanced vision systems, and interval management automation). In addition, a set of off-nominal situations were emulated. These included, for example, adverse weather conditions, traffic deviations, equipment failures, poor data quality, communication errors, and unexpected clearances, or changes to flight plans. Each situation was based on one or more reference events from past accidents or incidents, or on a similar case that had been used in previous developmental tests or studies. Over the course of the study, 10 twopilot airline crews participated, completing over 230 flights. Each flight consisted of an approach beginning at 10,000 ft. Based on the recorded data and pilot and research observations, preliminary results are presented regarding decision-making issues in the presence of the operational and technological complexities encountered during the flights.

Young, Steven D.↗

Computational Modeling of Mars Retropropulsion Concepts in the Langley Unitary Plan Wind Tunnel

Future human Mars missions will require powered descent starting at supersonic conditions, something which has never been done before at Mars. Computational powered descent flowfield simulations have been completed at full-scale Mars conditions, but the available ground test data are not appropriate for calibrating computational uncertainties for aerodynamic interference on proposed Mars descent vehicles. A test will be conducted in the NASA Langley Unitary Plan Wind Tunnel to begin addressing powered descent aerodynamics risks for large-scale human Mars entry concepts and to identify gaps in computational predictive capabilities. This paper covers pre-test computational flowfield predictions of two different models derived from full-scale reference vehicles: a blunt low lift-to-drag vehicle and a more slender geometry. Calculations of the blunt model include variations in nozzle configuration: nozzle location, size, area ratio, and pointing direction. There exist some significant differences between solvers, but some general trends are observed from simulations of the blunt model. First, aerodynamic axial force from the heatshield decreases with increasing thrust due to expanding plume blockage. Second, nozzles that point along the model axis result in lower aerodynamic axial force compared to nozzles that have a radial thrust component. Finally, placing the nozzles further away from the model nose preserves more heatshield axial force with increasing thrust compared to nozzles that are closer to the nose. For the slender model, the axial force from the heatshield is similar to the non-blowing axial force regardless of thrust magnitude, due to the nozzle arrangement on the heatshield. Once the test is completed, direct comparisons between the computations and test data will be made to determine computational uncertainties in a wind tunnel environment, to identify gaps in predictive capabilities, and to inform planning for future ground and flight test programs for Mars powered descent vehicles.

Supersonic Retropropulsion↗

Computational Modeling of Mars Retropropulsion Concepts in the Langley Unitary Plan Wind Tunnel

Future human Mars missions will require powered descent starting at supersonic conditions, something which has never been done before at Mars. Computational powered descent flowfield simulations have been completed at full-scale Mars conditions, but the available ground test data are not appropriate for calibrating computational uncertainties for aerodynamic interference on proposed Mars descent vehicles. A test will be conducted in the NASA Langley Unitary Plan Wind Tunnel to begin addressing powered descent aerodynamics risks for large-scale human Mars entry concepts and to identify gaps in computational predictive capabilities. This paper covers pre-test computational flowfield predictions of two different models derived from full-scale reference vehicles: a blunt low lift-to-drag vehicle and a more slender geometry. Calculations of the blunt model include variations in nozzle configuration: nozzle location, size, area ratio, and pointing direction. There exist some significant differences between solvers, but some general trends are observed from simulations of the blunt model. First, aerodynamic axial force from the heatshield decreases with increasing thrust due to expanding plume blockage. Second, nozzles that point along the model axis result in lower aerodynamic axial force compared to nozzles that have a radial thrust component. Finally, placing the nozzles further away from the model nose preserves more heatshield axial force with increasing thrust com- pared to nozzles that are closer to the nose. For the slender model, the axial force from the heatshield is similar to the non-blowing axial force regardless of thrust magnitude, due to the nozzle arrangement on the heatshield. Once the test is completed, direct comparisons between the computations and test data will be made to determine computational uncertainties in a wind tunnel environment, to identify gaps in predictive capabilities, and to inform planning for future ground and flight test programs for Mars powered descent vehicles.

Mars↗

Planck 2018 results XII. Galactic astrophysics using polarized dust emission

Observations of the submillimetre emission from Galactic dust, in both total intensity I and polarization, have received tremendous interest thanks to the Planck full-sky maps. In this paper we make use of such full-sky maps of dust polarized emission produced from the third public release of Planck data. As the basis for expanding on astrophysical studies of the polarized thermal emission from Galactic dust, we present full-sky maps of the dust polarization fraction p, polarization angle ψ, and dispersion function of polarization angles S . The joint distribution (one-point statistics) of p and NH confirms that the mean and maximum polarization fractions decrease with increasing NH. The uncertainty on the maximum observed polarization fraction, pmax = 22.0 -1.4 +3.5% at 353 GHz and 80' resolution, is dominated by the uncertainty on the Galactic emission zero level in total intensity, in particular towards diffuse lines of sight at high Galactic latitudes. Furthermore, the inverse behaviour between p and S found earlier is seen to be present at high latitudes. This follows the S ∝ p -1 relationship expected from models of the polarized sky (including numerical simulations of magnetohydrodynamical turbulence) that include effects from only the topology of the turbulent magnetic field, but otherwise have uniform alignment and dust properties. Thus, the statistical properties of p, ψ, and S for the most part reflect the structure of the Galactic magnetic field. Nevertheless, we search for potential signatures of varying grain alignment and dust properties. First, we analyse the product map S × p, looking for residual trends. While the polarization fraction p decreases by a factor of 3-4 between NH = 10 20 cm -2 and NH = 2 x 10 22 cm -2 , out of the Galactic plane, this product S × p only decreases by about 25%. Because S is independent of the grain alignment efficiency, this demonstrates that the systematic decrease in p with N H is determined mostly by the magnetic-field structure and not by a drop in grain alignment. This systematic trend is observed both in the diffuse interstellar medium (ISM) and in molecular clouds of the Gould Belt. Second, we look for a dependence of polarization properties on the dust temperature, as we would expect from the radiative alignment torque (RAT) theory. We find no systematic trend of S × p with the dust temperature T d , whether in the diffuse ISM or in the molecular clouds of the Gould Belt. In the diffuse ISM, lines of sight with high polarization fraction p and low polarization angle dispersion S tend, on the contrary, to have colder dust than lines of sight with low p and high S . We also compare the Planck thermal dust polarization with starlight polarization data in the visible at high Galactic latitudes. The agreement in polarization angles is remarkable, and is consistent with what we expect from the noise and the observed dispersion of polarization angles in the visible on the scale of the Planck beam. The two polarization emission-to-extinction ratios, R P/p and R S/V , which primarily characterize dust optical properties, have only a weak dependence on the column density, and converge towards the values previously determined for translucent lines of sight. We also determine an upper limit for the polarization fraction in extinction, p V /E(B - V), of 13% at high Galactic latitude, compatible with the polarization fraction p ≈ 20% observed at 353 GHz. Taken together, these results provide strong constraints for models of Galactic dust in diffuse gas.

79 ASTRONOMY AND ASTROPHYSICS↗

Use of GPS TEC Maps for Calibrating Single Band VLBI Sessions

GPS TEC ionosphere maps were first applied to a series of K and Q band VLBA astrometry sessions to try to eliminate a declination bias in estimated source positions. Their usage has been expanded to calibrate X-band only VLBI observations as well. At K-band, approx.60% of the declination bias appears to be removed with the application of GPS ionosphere calibrations. At X-band however, it appears that up to 90% or more of the declination bias is removed, with a corresponding increase in RA and declination uncertainties of approx.0.5 mas. GPS ionosphere calibrations may be very useful for improving the estimated positions of the X-only and S-only sources in the VCS and RDV sessions.

Gordon, David↗

New Tools for Automating Arcjet Sample Recession Tracking and Analysis

Arcjet Computer Vision (arcjetCV) has been significantly upgraded to enhance accuracy and performance in tracking material recession and shock-material standoff in test videos. These improvements include integrating new machine learning models, developing a specialized edge detection class, and incorporating a more comprehensive training dataset. These upgrades have refined the software’s ability to automate time-resolved recession tracking, making it more precise and reliable for analyzing complex physical processes. In parallel, a new tool called STARscan (Spatial Targeting and Alignment Rig for Scanning) is being developed to capture detailed 3D surface data before and after testing. By comparing these pre- and post-test scans with arcjetCV’s automated video analysis results, users can achieve a more comprehensive assessment of material recession. This method enables cross-validation of results, improving confidence in the analysis of tested materials. The expanded capabilities of arcjetCV have been successfully demonstrated on videos from various facilities, including the NASA Ames arcjets, UIUC’s PlasmatronX, and the VKI Plasmatron. It has been adopted as a new standard for in-situ recession tracking by the Mars Sample Return Project and Orion. ArcjetCV’s improved efficiency and accuracy are critical for reducing testing uncertainties and validating heatshield material performance under extreme conditions. The software’s user-friendly graphical interface ensures ease of use, enabling seamless processing and precise analysis of arcjet videos, providing deeper insights into material behavior in hypersonic environments. ArcjetCV is now available on both PyPI and Conda, allowing easy installation via "pip install arcjetCV" or through the Conda package manager, ensuring broad accessibility and streamlined deployment for users across various platforms.

Ablation↗

Inelastic dark matter scattering off Thallium cannot save DAMA

We study the compatibility of the observed DAMA modulation signal with inelastic scattering of dark matter (DM) off of the 0.1% Thallium (Tl) dopant in DAMA. In this work we test whether there exist regions of parameter space where the Tl interpretation gives a good fit to the most recent data from DAMA, and whether these regions are compatible with the latest constraints from other direct detection experiments. Previously, Chang et al. in 2010, had proposed the Tl interpretation of the DAMA data, and more recently (in 2019) the DAMA/LIBRA collaboration found regions in parameter space of Tl inelastic scattering that differ by more than 10σ from a no modulation hypothesis. We have expanded upon their work by testing whether the regions of parameter space where inelastic DM-Tl scattering gives a good fit to the most recent DAMA data survive the constraints placed by the lack of a DM signal in XENON1T and CRESST-II. In addition, we have tested how these regions change with the main sources of uncertainty: the Tl quenching factor, which has never been measured directly, and the astrophysical uncertainties in the DM distribution. We conclude that inelastic DM scattering off Tl cannot explain the DAMA data in light of null results from other experiments.

79 ASTRONOMY AND ASTROPHYSICS↗

Evaluation of Energy Storage Potential of Unconventional Shale Reservoirs Using Numerical Simulation of Cyclic Gas Injection

Compressed air energy storage (CAES) stores energy as compressed air in underground formations, typically salt dome caverns. When electricity demand grows, the compressed air is released through a turbine to produce electricity. CAES in the US is limited to one plant built in 1991, due in part to the inherent risk and uncertainty of developing subsurface storage reservoirs. As an alternative to CAES, we propose using some of the hundreds of thousands of hydraulically fractured horizontal wells to store energy as compressed natural gas in unconventional shale reservoirs. To store energy, produced or “sales” natural gas is injected back into the formation using excess electricity and is later produced through an expander to generate electricity. To evaluate this concept, we performed numerical simulations of cyclic natural gas injection into unconventional shale reservoirs using cmg-gem commercial reservoir modeling software. We tested short-term (diurnal) and long-term (seasonal) energy storage potential by modeling well injection and production gas flowrates as a function of bottom-hole pressure. First, we developed a conceptual model of a single fracture stage in an unconventional shale reservoir to characterize reservoir behavior during cyclic injection and production. Next, we modeled cyclic injection in the Marcellus shale gas play using published data. Results indicate that Marcellus unconventional shale reservoirs could support both short- and long-term energy storage at capacities of 100–1000 kWe per well. The results indicate that energy storage in unconventional shale gas wells may be feasible and warrants further investigation.

25 ENERGY STORAGE↗

Precipitation and Latent Heating Distributions from Satellite Passive Microwave Radiometry: Method and Uncertainties - Part 1

A revised Bayesian algorithm for estimating surface rain rate, convective rain proportion, and latent heating/drying profiles from satellite-borne passive microwave radiometer observations over ocean backgrounds is described. The algorithm searches a large database of cloud-radiative model simulations to find cloud profiles that are radiatively consistent with a given set of microwave radiance measurements. The properties of these radiatively consistent profiles are then composited to obtain best estimates of the observed properties. The revised algorithm is supported by an expanded and more physically consistent database of cloud-radiative model simulations. The algorithm also features a better quantification of the convective and non-convective contributions to total rainfall, a new geographic database, and an improved representation of background radiances in rain-free regions. Bias and random error estimates are derived from applications of the algorithm to synthetic radiance data, based upon a subset of cloud resolving model simulations, and from the Bayesian formulation itself. Synthetic rain rate and latent heating estimates exhibit a trend of high (low) bias for low (high) retrieved values. The Bayesian estimates of random error are propagated to represent errors at coarser time and space resolutions, based upon applications of the algorithm to TRMM Microwave Imager (TMI) data. Errors in instantaneous rain rate estimates at 0.5 deg resolution range from approximately 50% at 1 mm/h to 20% at 14 mm/h. These errors represent about 70-90% of the mean random deviation between collocated passive microwave and spaceborne radar rain rate estimates. The cumulative algorithm error in TMI estimates at monthly, 2.5 deg resolution is relatively small (less than 6% at 5 mm/day) compared to the random error due to infrequent satellite temporal sampling (8-35% at the same rain rate).

Olson, William S.↗

Precipitation and Latent Heating Distributions from Satellite Passive Microwave Radiometry: Improved Method and Uncertainties - Part 1

A revised Bayesian algorithm for estimating surface rain rate, convective rain proportion, and latent heating profiles from satellite-borne passive microwave radiometer observations over ocean backgrounds is described. The algorithm searches a large database of cloud-radiative model simulations to find cloud profiles that are radiatively consistent with a given set of microwave radiance measurements. The properties of these radiatively consistent profiles are then composited to obtain best estimates of the observed properties. The revised algorithm is supported by an expanded and more physically consistent database of cloud-radiative model simulations. The algorithm also features a better quantification of the convective and nonconvective contributions to total rainfall, a new geographic database, and an improved representation of background radiances in rain-free regions. Bias and random error estimates are derived from applications of the algorithm to synthetic radiance data, based upon a subset of cloud-resolving model simulations, and from the Bayesian formulation itself. Synthetic rain-rate and latent heating estimates exhibit a trend of high (low) bias for low (high) retrieved values. The Bayesian estimates of random error are propagated to represent errors at coarser time and space resolutions, based upon applications of the algorithm to TRMM Microwave Imager (TMI) data. Errors in TMI instantaneous rain-rate estimates at 0.5 -resolution range from approximately 50% at 1 mm/h to 20% at 14 mm/h. Errors in collocated spaceborne radar rain-rate estimates are roughly 50%-80% of the TMI errors at this resolution. The estimated algorithm random error in TMI rain rates at monthly, 2.5deg resolution is relatively small (less than 6% at 5 mm day.1) in comparison with the random error resulting from infrequent satellite temporal sampling (8%-35% at the same rain rate). Percentage errors resulting from sampling decrease with increasing rain rate, and sampling errors in latent heating rates follow the same trend. Averaging over 3 months reduces sampling errors in rain rates to 6%-15% at 5 mm day.1, with proportionate reductions in latent heating sampling errors.

Olson, William S.↗

Application of On-Line Monitoring and Real-Time Characterization of Low Level Waste Samples from Hanford Tanks

Waste and treatment stream characterization is an essential part of safe, efficient, and cost-effective processing of Hanford materials. The demand for materials characterization, including chemical, radionuclide, and physical attributes during waste management, transfer, and staging operations in the Hanford tank farms, is expected to significantly increase as the U.S. Department of Energy Office of River Protection (ORP) progresses toward Hanford waste treatment and immobilization. The application of on-line monitoring can greatly improve operation timeframes and reduce costs by providing real-time data and supporting the reduction of necessary grab sample collection. Recently, our team applied Raman online monitoring and chemometric modeling to the supernate of Hanford tank 241-AP-105, and demonstrated the quantitative measurement of nine analytes within this waste. Application of this method has been expanded to additional Hanford tank wastes, AW-102 and AP-107. The performance of the process monitoring instruments and the accuracy of the chemometric models will be assessed. On-line monitoring results will be compared to anion results determined from the current standard, ion-chromatographic and inductively coupled plasma mass spectrometry methods. Comparisons of both determined values and uncertainties of results will be presented. Results will be shared with personnel working on flowsheet interfaces for the Direct-Feed Low Activity Waste mission at Hanford.

Lines, Amanda M.↗

Application of On-Line Monitoring and Real-Time Characterization of Low Level Waste Samples from Hanford Tanks - 20380

Waste and treatment stream characterization is an essential part of safe, efficient, and cost-effective processing of Hanford materials. The demand for materials characterization, including chemical, radionuclide, and physical attributes during waste management, transfer, and staging operations in the Hanford tank farms, is expected to significantly increase as the U.S. Department of Energy Office of River Protection (ORP) progresses toward Hanford waste treatment and immobilization. The application of on-line monitoring can greatly improve operation time frames and reduce costs by providing real-time data and supporting the reduction of necessary grab sample collection. Recently, our team applied Raman online monitoring and chemometric modeling to the supernate of Hanford tank 241-AP-105, and demonstrated the quantitative measurement of nine analytes within this waste. Application of this method has been expanded to additional Hanford tank wastes, AW-102 and AP-107. The performance of the process monitoring instruments and the accuracy of the chemometric models will be assessed. On-line monitoring results will be compared to anion results determined from the current standard, ion-chromatographic and inductively coupled plasma mass spectrometry methods. Comparisons of both determined values and uncertainties of results will be presented. Results will be shared with personnel working on flowsheet interfaces for the Direct-Feed Low Activity Waste mission at Hanford. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Probabilistic flux limiters

The stable numerical integration of shocks in compressible flow simulations relies on the reduction or elimination of Gibbs phenomena (unstable, spurious oscillations). A popular method to virtually eliminate Gibbs oscillations caused by numerical discretization in under-resolved simulations is to use a flux limiter. A wide range of flux limiters have been studied in the literature, with recent interest in their optimization via machine learning methods trained on high-resolution datasets. The common use of flux limiters in numerical codes as plug-and-play blackbox components makes them key targets for design improvement. Even for deterministic dynamical models, numerical uncertainty is introduced via coarse-graining required by insufficient computational power to solve all scales of motion. Conventional flux limiters are deterministic and lack the capacity to address uncertainties, both aleatoric (inherent randomness) and epistemic (modeling uncertainty due to limited knowledge), which arise in coarse-grained numerical simulations. Here, we introduce a conceptually distinct type of flux limiter that is designed to handle the effects of randomness in the model and uncertainty in model parameters. Unlike traditional single-function flux limiters, these new probabilistic flux limiters incorporate multiple flux limiting functions, each applied with a learned probability drawn from high-resolution data to mitigate the effects of uncertainty in numerical simulations. This approach departs from traditional single-function limiters by explicitly modeling and incorporating uncertainty into the shock capturing process. Using the example of Burgers' equation as a testbed, we show that a machine learned, probabilistic flux limiter may be used in a shock capturing code to more accurately capture shock profiles. In particular, we show that our probabilistic flux limiter outperforms standard limiters and can be successively improved upon (up to a point) by expanding the set of probabilistically chosen flux limiting functions.

97 MATHEMATICS AND COMPUTING↗