SEARCH · Search NASA

Results for “dynamic parameter estimation”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,549 records · Page 86

All sources NASA NTRS DOE OSTI

NASA NTRS2003

A Solar System Survey of Forced Librations in Longitude

Physical librations in longitude are forced periodic variations of a body's rotation rate. If the torque producing the librations can be calculated, then observations of the phase and amplitude of librations can provide information on mass distribution, and effective strength of the body. In the near future prospects for observing physical librations look quite promising. Radio interferometric observations of Venus and Mercury may yield sufficiently accurate rotational observations that librations there may be visible. Range measurements from Earth to networks of landed instrument packages on Mars are likely to yield librational data there as well. We compute expected libration amplitudes from physical and orbital parameters of a set of planets and satellites partially motivated by a desire to identify candidates for future observations. Solar system bodies occupy one of three general rotation states: non-resonant states, resonant states, and the synchronous resonant state. Analytical treatments of forced librations were initially motivated by the Moon. Lunar librations were predicted by Newton, first detected telescopically by Bessel, and definitively resolved through lunar laser ranging which has led to quite thorough analysis of of librations for the synchronous case. The synchronous resonant state is commonly observed among satellites. The only known body to exist in a non-synchronous resonance is Mercury which exists in a 3:2 resonance, completing three rotations for every two revolutions about the sun. The analysis of Goldreich and Peale has lead to improved understanding of the general case of half integer resonance states. The dynamics of forced librations in non-resonant rotators has received less attention. While there are few cases in which non-resonant forced librations have been observed, Earth is an important exception, and current observing techniques may have the capacity to detect them on Venus. A comprehensive observing program spanning a range of solar system bodies can address an array of geophysical issues involving interior mass distribution of planets, satellites, and asteroids. Calculations of expected librations can supply amplitude estimates helpful in identifying the likelihood of detecting librations observationally. We present a survey of expected libration amplitudes for a subset of solar system bodies, identifying those bodies with amplitudes likely to be detectable, and commenting on spin state, and radial structure implications.

Comstock, R. L.↗

DOE OSTI2022

NPP Simulators for Coupled Thermal and Electric Power Dispatch

The Light Water Reactor Sustainability (LWRS) program within the United States Department of Energy supports extending the operation of the U.S. commercial nuclear power plant (NPP) fleet. Within the LWRS program, the Flexible Plant Operation and Generation (FPOG) Pathway works to diversify the revenue streams of light water reactors (LWRs) by opening opportunities for the co-generation of non-electric products in addition to supplying electrical power to the grid. Recent events have added greater motivation to these efforts. For example, the recent Inflation Reduction Act (IRA) passed by the U.S. federal government offers substantial tax incentives for producing clean hydrogen, the technology readiness level of dispatchable and high-efficiency hydrogen production has dramatically increased in a short time, and societal response to world climate change is driving a transition away from fossil fuels. Producing hydrogen with maximum efficiency using nuclear power requires dispatching both electrical and thermal power from the nuclear plant to the hydrogen plant, so testing concepts of operations for combined electrical and thermal power dispatch (TPD) from an NNP to a hydrogen plant is of interest. This report documents achievement of the Light Water Reactor Sustainability (LWRS) program milestone “Install and demonstrate a vendor-developed simulator on the Human Systems Simulation Laboratory (HSS) for dispatch of LWR electrical power to a close-coupled electrolysis plant” with a due date of Dec. 22, 2022. Several factors provide motivation for this effort. Coupling the power generation deck of a nuclear power plant to a hydrogen production facility introduces new possibilities for operational transients that must be addressed. In particular, the performance of the integrated system during startup and shutdown of the hydrogen production facility, as well as offnormal conditions, need to be evaluated to ensure there are no adverse effects on the operation of the existing NPP. The concept of operations involving the NPP, the hydrogen plant, and the electric power grid must be tested using NPP simulators and operating procedures that have been modified for TPD operations. These tests must also include dynamic simulations of the coupled tertiary thermal and electric loads as well as coordinated activities with NPP operators, tertiary load operators and grid power coordinators. The report summarizes progress in developing and testing full-scope NPP simulators at the HSSL, including a generic BWR simulator from GSE Systems, Inc. and generic PWR simulator from Westinghouse. In the case of the TPD-GBWR Simulator from GSE Systems, Inc., a BWR is thermally coupled to a high temperature electrolysis (HTE) plant that produces hydrogen and oxygen from de-ionized water. The hydrogen plant is not explicitly simulated but only included as a transient heat sink. A thermal power dispatch (TPD) system transfers heat between the steam systems at the BWR and the hydrogen plant. Operational results from two versions of the modified simulator are presented. The first version uses synthetic oil as a heat transfer fluid in a closed delivery heat loop (DHL) that generates steam at the hydrogen plant. The second version uses steam as the heat transfer fluid in a delivery steam line (DSL) to provide steam to the hydrogen plant. For both versions, the estimated thermal power delivery distance is approximately one kilometer. The amount of thermal power dispatched in the simulators is 15% of the total reactor thermal power such that the simulators provide a tool to study the feasibility of coupling a BWR to industrial processes that benefit from a combination electrical and thermal power dispatch. Ongoing work within a CRADA is also developing a full-scope PWR simulator provided by Westinghouse for both thermal and electric power coupling. This simulator is based on a PWR plant with two three-loop Westinghouse reactors. Westinghouse PWRs are sufficiently similar that a simulator of a three-loop reactor is an appropriate representation for two-loop and four-loop PWR reactors. The three-loop simulator will initially be modified for close-coupling to a 100 MW HTE hydrogen production plant that will require approximately 25 MW of thermal power while operating at its maximum rated capacity. The simulator testing will include full coupling to dynamic simulations of a hydrogen production plant and a representative bulk electric grid. The simulator provided by Westinghouse is similar to the GPWR simulator that INL has already obtained from GSE Systems but has a few important added benefits. First, the Westinghouse simulator is based on digital controls and has additional screens that can be called up to show parameter trends to assist operators in decision-making. The Westinghouse simulator also has upgrades to the controls and hardware representations, such as valve actuators, that make it more realistic and flexible in terms of accurately sim

99 GENERAL AND MISCELLANEOUS↗

DOE OSTI2023

Multifidelity Approach to Sensitivity Estimation in Large-Eddy Simulation

A novel approach to compute affordable approximate sensitivities in a large-eddy simulation (LES) is proposed and assessed. The approach is based on solving a Reynolds-averaged Navier–Stokes (RANS) problem that has been linearized around the mean LES solution, with closure modeling required for the linearized changes in turbulent Reynolds stresses. In the present study, the closure modeling is based on a linearized algebraic turbulence model. The method is assessed for the flow over a NACA0012 airfoil at a fixed angle of attack, with the Reynolds number as the variable parameter. The results show that, provided an accurate linearized closure model, the method predicts the correct sensitivity of the skin friction coefficient and the mean velocity field at the cost of a linearized RANS, which provides an important proof-of-concept for this approach. Here, the linearized algebraic turbulence model with standard model coefficients produces a reasonably accurate sensitivity, but the results also suggest that accuracy could be gained from recalibrating the model coefficients for this new use of the model.

42 ENGINEERING↗

DOE OSTI2017

Shock-Driven Hydrodynamic Instability Growth Near Phase Boundaries and Material Property Transitions: Final Report

Activities for this grant included: 1) Development of dynamic impact experiments to probe strength and phase transition influence on dynamic deformation, 2) development of modern strength and phase aware simulation capabilities, 3) and post-processing of experimental data with simulation and closed form analytical techniques. Two different dynamic experiments were developed to probe material strengths in solid metals (largely copper and iron in this effort). In the first experiment a flyer plate impacts a flat target with an opposite rippled surface that is partially supported by a weaker window material. Post mortem analysis of the target sample showed a strong and repeatable residual plastic deformation dependence on grain orientation. Yield strengths for strain rates near 10 5 s -1 and plastic strains near ~50% were estimated to be around 180 to 240 MPa, varying in this range with grain orientation. Unfortunately dynamic real-time measurements were difficult with this setup due to diagnostic laser scattering; hence, an additional experimental setup was developed to complement these results. In the second set of experiments a rippled surface was ablated by a controlled laser pulsed, which launched a rippled shock front to an opposite initially flat diagnostic surface that was monitored in real-time with spatially resolved velocimetry techniques, e.g., line VISAR in addition to Transient Imaging Displacement Interferometry (TIDI) displacement measurements. This setup limited the displacements at the diagnostic surface to a reasonable level for TIDI measurements (~ less than one micrometer). These experiments coupled with analytical and numerical solutions provided evidence that viscous and elastic deviatoric strength affect shock front perturbation evolution in clearly different ways. Particularly, normalized shock front perturbation amplitudes evolve with viscosity (η) and perturbation wavelength (λ) as η/λ, such that increasing viscosity (or decreasing the initial wavelength) delays the perturbation decay. Conversely our experimental data, analysis and simulations show that for materials with elastic yield strength Y the normalized shock perturbation amplitude evolves with Yλ/A 0 , which shows wavelength increases have the opposite effect as in viscous materials and perturbation decay is also dependent on initial amplitude A 0 (viscous materials are independent of this parameter). Materials where strength had clear strain rate dependence, e.g., such as a PTW material law, behaved similarly to materials with only an effective yield stress (elastic-perfectly plastic) in the shock front perturbation studies obeying a Y eff λA 0 relationship where Y eff was a constant (near ~400 MPa for Cu for strain rates around 10 6 s -1 ). Magnitude changes in strain rate would increase Y eff as would be expected from the PTW behavior, but small perturbations (typical of regions behind the shock front) near a mean had little effect. Additional work based on simulations showed that phase transformation kinetics can affect the behavior of the perturbed shock front as well as the evolution of the RM-like instability that develops due to the imprint of the perturbed shock front on the initially flat surface as the shock breaks out.

36 MATERIALS SCIENCE↗

NASA NTRS1991

Modeling the spatiotemporal organization of velocity storage in the vestibuloocular reflex by optokinetic studies

1. A generalized three-dimensional state space model of visual vestibular interaction was developed. Matrix and dynamical system operators associated with inputs from the semicircular canals, otolith velocity estimator, and the visual system have been incorporated into the model, which focus on their relationship to the velocity storage integrator. 2. A relationship was postulated between the eigenvalues and the direction of the eigenvectors of the system matrix and the orientation of the spatial vertical. It was assumed that the system matrix for a tilted position was a composition of two linear transformations of the system matrix for the upright position. One transformation modifies the eigenvalues of the system matrix, whereas another rotates the eigenvectors of the system matrix. The pitch and roll eigenvectors rotate with the head, whereas the yaw axis eigenvector remains approximately spatially invariant. 3. Based on the three-dimensional model, a computational procedure was formulated to identify the eigenvalues and eigenvectors of the system matrix with the use of a modification of the marquardt algorithm. With the use of data obtained from a monkey, it was shown that the three-dimensional behavior of velocity storage cannot be predicted solely in terms of its time constants, i.e., the inverse of its eigenvalues. With the use of the same eigenvalues the data could either be fit or not fit, depending on the eigenvector directions. Therefore, it is necessary to specify eigenvector directions when characterizing velocity storage in three dimensions. 4. Parameters found with the use of the Marquardt algorithm were incorporated into the model. Diagonal matrices in a head coordinate frame were introduced for coupling the visual system to the integrator and to the direct optokinetic pathway. Simulations of optokinetic nystagmus (OKN) and optokinetic after-nystagmus (OKAN) were run. The model predicted the behavior of yaw and pitch OKN and OKAN when the animal is upright. It also predicted the cross-coupling in the side down position. The trajectories in velocity space were also accurately simulated. 5. One of the predictions of the model is that when the stimulus direction is along an eigenvector, the trajectory in velocity space is a straight line. Using the "spectral width" of the residuals from a straight line sequence during OKAN, we developed a methodology to estimate how close the OKAN decay was to an eigenvector trajectory. 6. Thus we have developed a model-based approach for studying and interpreting the response characteristics of velocity storage in three dimensions.(ABSTRACT TRUNCATED AT 400 WORDS).

Non-NASA Center↗

NASA NTRS2002

Hierarchies of Models: Toward Understanding Planetary Nebulae

Stars like our sun (initial masses between 0.8 to 8 solar masses) end their lives as swollen red giants surrounded by cool extended atmospheres. The nuclear reactions in their cores create carbon, nitrogen and oxygen, which are transported by convection to the outer envelope of the stellar atmosphere. As the star finally collapses to become a white dwarf, this envelope is expelled from the star to form a planetary nebula (PN) rich in organic molecules. The physics, dynamics, and chemistry of these nebulae are poorly understood and have implications not only for our understanding of the stellar life cycle but also for organic astrochemistry and the creation of prebiotic molecules in interstellar space. We are working toward generating three-dimensional models of planetary nebulae (PNe), which include the size, orientation, shape, expansion rate and mass distribution of the nebula. Such a reconstruction of a PN is a challenging problem for several reasons. First, the data consist of images obtained over time from the Hubble Space Telescope (HST) and spectra obtained from Kitt Peak National Observatory (KPNO) and Cerro Tololo Inter-American Observatory (CTIO). These images are of course taken from a single viewpoint in space, which amounts to a very challenging tomographic reconstruction. Second, the fact that we have two disparate and orthogonal data types requires that we utilize a method that allows these data to be used together to obtain a solution. To address these first two challenges we employ Bayesian model estimation using a parameterized physical model that incorporates much prior information about the known physics of the PN. In our previous works we have found that the forward problem of the comprehensive model is extremely time consuming. To address this challenge, we explore the use of a set of hierarchical models, which allow us to estimate increasingly more detailed sets of model parameters. These hierarchical models of increasing complexity are akin to scientific theories of increasing sophistication, with each new model/theory being a refinement of a previous one by either incorporating additional prior information or by introducing a new set of parameters to model an entirely new phenomenon. We apply these models to both a simulated and a real ellipsoidal PN to initially estimate the position, angular size, and orientation of the nebula as a two-dimensional object and use these estimates to later examine its three-dimensional properties. The efficiency/accuracy tradeoffs of the techniques are studied to determine the advantages and disadvantages of employing a set of hierarchical models over a single comprehensive model.

Knuth, Kevin H.↗

NASA NTRS2003

Hierarchies of Models: Toward Understanding Planetary Nebulae

Knuth, Kevin H.↗

DOE OSTI2016

The trends high-contrast imaging survey. VI. Discovery of a mass, age, and metallicity benchmark brown dwarf

The mass and age of substellar objects are degenerate parameters leaving the evolutionary state of brown dwarfs ambiguous without additional information. Theoretical models are normally used to help distinguish between old, massive brown dwarfs and young, low-mass brown dwarfs but these models have yet to be properly calibrated. We have carried out an infrared high-contrast imaging program with the goal of detecting substellar objects as companions to nearby stars to help break degeneracies in inferred physical properties such as mass, age, and composition. Rather than using imaging observations alone, our targets are pre-selected based on the existence of dynamical accelerations informed from years of stellar radial velocity (RV) measurements. In this paper, we present the discovery of a rare benchmark brown dwarf orbiting the nearby (d = 18.69 ± 0.19 pc), solar-type (G9V) star HD 4747 ([Fe/H] = −0.22 ± 0.04) with a projected separation of only ρ = 11.3 ± 0.2 au (θ = 0.″6). Precise Doppler measurements taken over 18 years reveal the companion’s orbit and allow us to place strong constraints on its mass using dynamics (msin i=55.3±1.9M{sub Jup}). Relative photometry (ΔK {sub s} = 9.05 ± 0.14, M{sub K{sub s}}=13.00±0.14, K{sub s}−L{sup ′}=1.34±0.46) indicates that HD 4747 B is most likely a late-type L-dwarf and, if near the L/T transition, an intriguing source for studying cloud physics, variability, and polarization. We estimate a model-dependent mass of m=72{sub −13}{sup +3} M{sub Jup} for an age of 3.3{sub −1.9}{sup +2.3} Gyr based on gyrochronology. Combining astrometric measurements with RV data, we calculate the companion dynamical mass (m=60.2±3.3M{sub Jup}) and orbit (e = 0.740 ± 0.002) directly. As a new mass, age, and metallicity benchmark, HD 4747 B will serve as a laboratory for precision astrophysics to test theoretical models that describe the emergent radiation of brown dwarfs.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Bridging the time scale in exascale computing of chemical systems (Final Technical Report)

This report summarizes the work carried out with support of the United States Department of Energy under Award DE-SC0019441. The theme of this project was to develop and apply methods that allowed for the acceleration of atomistic calculations, particularly in challenging areas such as multiphase systems, electrified interfaces, uncertainty estimation, and applications requiring chemical accuracy, which tend to be applications where simulation time is severely bottlenecked by the computational time requirements. Much of the focus was on the application of emerging machine-learning methodologies, although a wide range of methodologies were employed. This report has two major sections. The first focuses on the methodological advances themselves. Within this part, we report a number of major advances, a few examples of which are described here. We report the first machine-learning scheme for the acceleration of electronically grand-canonical calculations (that is, those applicable to electrochemistry). We report new methods of performing transfer learning, in which physics-based priors can be used to provide predictions, often with uncertainty estimates, of images well outside of training sets; we also offer ways to fine-tune these transfer-learning models. We provide a new systematic means to generate and apply minimal training data sets to very large (10,000’s of atoms) systems, with only small training sets appropriate for electronic structure. We developed new methodologies to integrate surface vibrations into surface adsorption calculations. We made advances to the applicability of diffusion Monte Carlo methods to allow (learned) force prediction, finite-size error correction, and force-free means of searching for transition states. We integrated machine-learned atomistic predictions into mechanism generation codes. Additionally, we released new software including AmpTorch, a modernized version of our original atomistic machine-learning code Amp. The second part of this report focuses on the scientific applications that accompanied, and were often enabled by, the methodological advances described earlier. A few examples follow, but full details are in the individual chapters of the report. For example, we developed a general theory of phonon-induced friction on molecular adsorbates. We showed fundamentally how solvent influences the adsorption and desorption process and how it differs from the processes typically involved at the solid–gas interface, making aqueous-phase and electrocatalysis different from traditional thermocatalysis. We examined how metal–insulator and magnetic transitions can be probed, and accelerated exciton dynamics via Frenkel Hamiltonian parameters. We showed that the nearsighted force-training approach, developed within this project, can predict both the stability and reactivity of large nanoparticles, and can also lead to insights on catalyst coverage on binding energies and entropies. These applied studies, which generally integrated with our method development, allowed us to push forward the theoretical understanding of several reaction classes.

08 HYDROGEN↗

NASA NTRS2001

Observations of Ocean Primary Productivity Using MODIS

Measuring the magnitude and variability of oceanic net primary productivity (NPP) represents a key advancement toward our understanding of the dynamics of marine ecosystems and the role of the ocean in the global carbon cycle. MODIS observations make two new contributions in addition to continuing the bio-optical time series begun with Orbview-2's SeaWiFS sensor. First, MODIS provides weekly estimates of global ocean net primary productivity on weekly and annual time periods, and annual empirical estimates of carbon export production. Second, MODIS provides additional insight into the spatial and temporal variations in photosynthetic efficiency through the direct measurements of solar-stimulated chlorophyll fluorescence. The two different weekly productivity indexes (first developed by Behrenfeld & Falkowski and by Yoder, Ryan and Howard) are used to derive daily productivity as a function of chlorophyll biomass, incident daily surface irradiance, temperature, euphotic depth, and mixed layer depth. Comparisons between these two estimates using both SeaWiFS and MODIS data show significant model differences in spatial distribution after allowance for the different integration depths. Both estimates are strongly dependence on the accuracy of the chlorophyll determination. In addition, an empirical approach is taken on annual scales to estimate global NPP and export production. Estimates of solar stimulated fluorescence efficiency from chlorophyll have been shown to be inversely related to photosynthetic efficiency by Abbott and co-workers. MODIS provides the first global estimates of oceanic chlorophyll fluorescence, providing an important proof of concept. MODIS observations are revealing spatial patterns of fluorescence efficiency which show expected variations with phytoplankton photo-physiological parameters as measured during in-situ surveys. This has opened the way for research into utilizing this information to improve our understanding of oceanic NPP variability. Deriving the ocean bio-optical properties places severe demands on instrument performance (especially band to band precision) and atmospheric correction. Improvements in MODIS instrument characterization and calibration over the first 16 mission months have greatly improved the accuracy of the chlorophyll input fields and FLH, and therefore the estimates of NPP and fluorescence efficiency. Annual estimates now show the oceanic NPP accounts for 40-50% of the global total NPP, with significant interannual variations related to large scale ocean processes. Spatial variations in ocean NPP, and exported production, have significant effects on exchange of CO2 between the ocean and atmosphere. Further work is underway to improve both the primary productivity model functions, and to refine our understanding of the relationships between fluorescence efficiency and NPP estimates. We expect that the MODIS instruments will prove extremely useful in assessing the time dependencies of oceanic carbon uptake and effects of iron enrichment, within the global carbon cycle.

Esaias, Wayne E.↗

NASA NTRS1998

Final Report

This grant was made to fund two accepted proposals to observe with the ASCA satellite: (1) Debunking the myth of two-temperature coronae of active stars, and (2) Dynamic coronae of dMe stars. We obtained the requested observations and have now completed two papers that will be published in the Astrophysical Journal. The RS CVn binary star UX Ari was observed for 14 hours with all four detectors onboard the ASCA satellite. The X-ray emission was at a constant, quiescent level during the first 12 hours, after which time a powerful flare with a peak luminosity of 1.4 x 10(exp 32) ergs/s started. The flare was observed until shortly after its peak. We present a spectral and temporal analysis of the UX Ari observations and analyze the data with a two-ribbon flare model including estimates for cooling losses. A time-dependent reconstruction of the emission measure (EM) distribution shows that two separate plasma components evolve during the flare (one being identified with the quiescent EM). Most of the flare EM reaches temperatures between 50 MK and 100 MK or more. Magnetic confinement requires the loop arcade to be geometrically large, with length scales on the order of one stellar radius. The electron densities inferred from the model decrease from initial values around 10(exp 12)/cc early in the flare to about 10(exp 11)/cc at the flare peak. The best-fit models require surface magnetic field strengths of a few hundred G, compatible with the maximum photospheric fields expected from equipartition. The flare parameters imply a (conductive and radiative) cooling loss time of less than one hour at flare peak. The elemental abundances increase significantly during the flare rise, with the abundances of the low-FIP elements Fe, Mg, Si, and Ni typically increasing to higher levels than the high-FIP elements such as S or Ne. The Fe abundance increases from (17 +/- 4)% of the solar photospheric value during quiescence up to (89 +/- 18)% at flare peak. A fractionation process that occurs during the chromospheric evaporation phase may selectively enrich low-FIP elements as in the solar corona; alternatively, the chromospheric evaporation may itself bring metal-rich plasma into the metal-poor corona.

Linsky, Jeffrey L.↗

DOE OSTI2020

Diffusion coefficients in the envelopes of white dwarfs

The diffusion of elements is a key process in understanding the unusual surface composition of white dwarfs and their spectral evolution. The diffusion coefficients of Paquette et al. have been widely used to model diffusion in white dwarfs. In this work, we perform new calculations of the coefficients of interdiffusion and ionic thermal diffusion with (1) a more advanced model that uses a recent modification of the calculation of the collision integrals that is more suitable for the partially ionized, partially degenerate, and moderately coupled plasma and (2) classical molecular dynamics. The coefficients are evaluated for silicon and calcium in white dwarf envelopes of hydrogen and helium. A comparison of our results with Paquette et al. shows that the latter systematically underestimates the coefficient of interdiffusion yet provides reliable estimates for the relatively weakly coupled plasmas found in nearly all types of stars, as well as in white dwarfs with hydrogen envelopes. In white dwarfs with cool helium envelopes (${T}_{\mathrm{eff}}$ < 15,000 K), the difference grows to more than a factor of two. We also explored the effect of the ionization model used to determine the charges of the ions and found that it can be a substantial source of discrepancy between different calculations. Finally, we consider the relative diffusion timescales of Si and Ca in the context of the pollution of white dwarf photospheres by accreted planetesimals and find factor of ≳3 differences between calculations based on Paquette et al. and our model.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2019

Radiological Performance Modeling: Platform Comparison

Zapadnoe Uranium Mill Tailings Site, Pridneprovsky Chemical Plant Ukraine: Two exposure scenarios were utilized to develop performance assessment models in GoldSim and NORMALYSA. The first scenario is an onsite worker exposed to radioactivity through external ground irradiation as well as radon and radioactive dust inhalation. The second being an offsite Resident 300 years after institutional control exposed through the ingestion of garden vegetables irrigated by contaminated ground water, ingestion of fish from Konoplyanka River, and agricultural activities in contaminated garden. NORMALYSA utilizes a modular interface to connect extensive pre-created modules to develop contaminant transport and risk assessment models for safety assessments and performance assessments at NORM and legacy sites. Developed by Facilia AB, the program was tested with the creation of the Zapadnoe model under the International Atomic Energy Agency MODARIA Il Working Group 1 (WG1). Default values are provided for input parameters where site specific data and distributions lack. Probabilistic simulations and inputs are limited to parameters with assigned Gaussian probability functions. Similar results for transient worker inhalation dose were found in comparison with the Zapadnoe GoldSim model. However, the NORMALYSA model has results larger by an order of magnitude for external radiation in the transient worker and for aquatic food ingestion, inhalation, and garden food ingestion to the resident. Performance Assessment Modeling: Radiological performance models promote good management and support in remediation and environmental decision making by optimizing radioactive waste disposal site performance. a performance assessment model of the Zapadnoe uranium mill tailings site and a model for calculating authorized limits at Los Alamos National Laboratory are created to demonstrate and highlight the respective capabilities of three modeling platforms: GoldSim, NORMALYSA, and RESRAD-Onsite. GoldSim is a general-purpose modeling platform with the ability to simulate site-specific FT radionuclide transport and risk assessment using a graphical interface to create and link elements that act as indivisible building blocks of the model. The Zapadnoe performance assessment model created in GoldSim is composed of input parameters collected from site specific parameters from WG1 as well as standard assumptions seen in the Generic PA Model from Neptune and Company. The multivariate Zapadnoe simulation contains both deterministic and probabilistic simulation capabilities with the ability to quantitatively represent the uncertainty and associated risks of the site via the Monte Carlo simulation method. Model parameters, documentation, and assumptions are built directly into the model creating a transparent, reproducible, and more purposeful site specific model. To demonstrate this, an ingestion of water pathway was included in addition to the exposure pathways included in the NORMALYSA model. GoldSim offers a wide variety of probability distributions to choose from as well as many tools to visualize the uncertainty of the models and explore which pathways and radionuclides dominate dose. In contrast to the NORMALYSA model, the GoldSim model uses dynamic simulation to calculate tritium concentrations. Authorized limits support the evaluation of radiological health effects from radionuclides present in bulk contaminated soil in order to assist in the process of clearing land for release under different uses. Calculations of authorized limits in RESRAD-Onsite and GoldSim models support the evaluation for 15 radionuclides using a residential exposure scenario, indoor and outdoor inhalation of dust and vapor phase radionuclides (tritium), incidental soil ingestion, and ingestion grown fruits and Vegetables. RESRAD-Onsite, developed by Argonne National Laboratory is a dedicated computer code for calculating radiation doses and cancer risks to an individual located on top of radioactively contaminated soils. It can also be used to provide guidelines for deriving radionuclide soil levels corresponding to a specific dose criterion. The menu driven user interface facilitates selection of radionuclides and soil concentrations, as well as defines parameter values related to materials, environmental factors, and exposure pathways. The program contains radionuclide-specific data such as dose conversion factors and risk coefficients, as well as soil-water partition coefficients and transfer factors for biotic media. The GoldSim model successfully emulates the RESRAD-Onsite code using identical radiological decay and ingrowth factors, radiation dose conversion factors, and other input parameters. DE EEDE ETS Each model supports dynamic modeling of radionuclide ingrowth and decay over time, and integrate this capability with the transport of radionuclides in the environment via processes such as advection and/or diffusion in fluids, or within porous solid media, The applicability of each program is heavily dependent on project context. In the context of developing Authorized Limits, RESRAD-Onsite may be a more desirable option as it was developed for this particular application. For simple disposal or legacy sites or sites that do not need site-specific model structure, NORMALYSA may be a desirable modeling platform, GoldSim can be used to model any problem for which RESRAD or NORMALYSA is applicable, but its strength lies in its flexibility, analytical, and probabilistic capabilities.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

NASA NTRS2011

Investigating On-Orbit Attitude Determination Anomalies for the Solar Dynamics Observatory Mission

The Solar Dynamics Observatory (SDO) was launched on February 11, 2010 from Kennedy Space Center on an Atlas V launch vehicle into a geosynchronous transfer orbit. SDO carries a suite of three scientific instruments, whose observations are intended to promote a more complete understanding of the Sun and its effects on the Earth's environment. After a successful launch, separation, and initial Sun acquisition, the launch and flight operations teams dove into a commissioning campaign that included, among other things, checkout and calibration of the fine attitude sensors and checkout of the Kalman filter (KF) and the spacecraft s inertial pointing and science control modes. In addition, initial calibration of the science instruments was also accomplished. During that process of KF and controller checkout, several interesting observations were noticed and investigated. The SDO fine attitude sensors consist of one Adcole Digital Sun Sensor (DSS), two Galileo Avionica (GA) quaternion-output Star Trackers (STs), and three Kearfott Two-Axis Rate Assemblies (hereafter called inertial reference units, or IRUs). Initial checkout of the fine attitude sensors indicated that all sensors appeared to be functioning properly. Initial calibration maneuvers were planned and executed to update scale factors, drift rate biases, and alignments of the IRUs. After updating the IRU parameters, the KF was initialized and quickly reached convergence. Over the next few hours, it became apparent that there was an oscillation in the sensor residuals and the KF estimation of the IRU bias. A concentrated investigation ensued to determine the cause of the oscillations, their effect on mission requirements, and how to mitigate them. The ensuing analysis determined that the oscillations seen were, in fact, due to an oscillation in the IRU biases. The low frequencies of the oscillations passed through the KF, were well within the controller bandwidth, and therefore the spacecraft was actually following the oscillating biases, resulting in movement of the spacecraft on the order of plus or minus 20 arcsec. Though this level of error met the ACS attitude knowledge requirement of [35, 70, 70] arcsec, 3 sigma, the desire of the ACS and instrument teams was to remove as much of the oscillation as possible. The Kearfott IRUs have an internal temperature controller, designed to maintain the IRU temperature at a constant temperature of approximately 70 C, thus minimizing the change in the bias drift and scale factors of the mechanical gyros. During ground testing of the observatory, it was discovered that the 83-Hz control cycle of the IRU heaters put a tremendous amount of stress on the spacecraft battery. Analysis by the power systems team indicated that the constant charge/discharge on the battery due to the IRU thermal control cycle could potentially limit the life of the battery. After much analysis, the decision was made not to run the internal IRU heaters. Analysis of on orbit data revealed that the oscillations in the IRU bias had a connection to the temperature of the IRU; changes in IRU temperature resulted in changes in the amplitude and period of the IRU biases. Several mitigating solutions were investigated, the result of which was to tune the KF with larger IRU noise assumptions which allows the KF to follow and correct for the time-varying IRU biases.

Vess, Melissa F.↗

DOE OSTI2018

An evaluation of equilibrium conditions and temperature-dependent speciation in a laser-produced air plasma

A laser-induced air breakdown plasma is a dynamical system and its fundamental parameters change significantly during its lifetime. The temporal evolution of the spectral features from an air plasma shows emission from bremsstrahlung followed by discrete electronic transitions from ions, atoms, and rotational and vibrations bands from molecules. The molecular band emission from air plasma typically appears at times ≳5 μs and persists for 100s of μs. Here, an evaluation is made on the temperature evolution and speciation in a laser-produced air plasma. The air plasmas are generated by focusing 1064 nm, 6 ns pulses from an Nd:YAG laser. High-resolution emission spectra of various molecules such as N 2 , N 2 + , CN, OH, NH, NO are acquired for fits to spectral models. Fitting of atomic and molecular emission features permits tracking of the air spark temperature evolution from 1 to 200 µs. Though the excitation and molecular temperatures show good overlap at times ≲10 μs, the molecular temperatures obtained from different species show a discontinuity at times ~ 30 μs and this is related to shock collapse and subsequent changes in hydrodynamics and chemistry of the plume. The fitting of multiple species in broadband spectra has permitted calculation of the relative concentrations of various molecules as a function of temperature that provides insight into the air spark speciation. Finally, the measured mole fractions from laser-induced air plasma are within a factor of 2 - 6 of the values estimated by speciation model under the assumption of thermal and chemical equilibrium across the temperature ranges studied and probable reasons for the observed deviation are described.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS2014

Mass and Volume Optimization of Space Flight Medical Kits

Resource allocation is a critical aspect of space mission planning. All resources, including medical resources, are subject to a number of mission constraints such a maximum mass and volume. However, unlike many resources, there is often limited understanding in how to optimize medical resources for a mission. The Integrated Medical Model (IMM) is a probabilistic model that estimates medical event occurrences and mission outcomes for different mission profiles. IMM simulates outcomes and describes the impact of medical events in terms of lost crew time, medical resource usage, and the potential for medically required evacuation. Previously published work describes an approach that uses the IMM to generate optimized medical kits that maximize benefit to the crew subject to mass and volume constraints. We improve upon the results obtained previously and extend our approach to minimize mass and volume while meeting some benefit threshold. METHODS We frame the medical kit optimization problem as a modified knapsack problem and implement an algorithm utilizing dynamic programming. Using this algorithm, optimized medical kits were generated for 3 mission scenarios with the goal of minimizing the medical kit mass and volume for a specified likelihood of evacuation or Crew Health Index (CHI) threshold. The algorithm was expanded to generate medical kits that maximize likelihood of evacuation or CHI subject to mass and volume constraints. RESULTS AND CONCLUSIONS In maximizing benefit to crew health subject to certain constraints, our algorithm generates medical kits that more closely resemble the unlimited-resource scenario than previous approaches which leverage medical risk information generated by the IMM. Our work here demonstrates that this algorithm provides an efficient and effective means to objectively allocate medical resources for spaceflight missions and provides an effective means of addressing tradeoffs in medical resource allocations and crew mission success parameters.

Monte Carlo Method↗

DOE OSTI2021

Active microbial biomass decreases, but microbial growth potential remains similar across soil depth profiles under deeply-vs. shallow-rooted plants

Climate-smart land management practices that replace shallow-rooted annual crop systems with deeply-rooted perennial plants can contribute to soil carbon sequestration. However, deep soil carbon accrual may be influenced by active microbial biomass and their capacity to assimilate fresh carbon at depth. Incorporating active microbial biomass, dormancy, and growth in microbially-explicit models can improve our ability to predict soil's capacity to store carbon. But, so far, the microbial parameters that are needed for such modeling are poorly constrained, especially in deep soil layers. Here, we used a lab incubation experiment and growth kinetics model to estimate how microbial parameters vary along 240 cm of soil depth in profiles under shallow- (soy) and deeply-rooted (switchgrass) plants 11 years after plant cover conversion. We also assessed resource origin and availability (total organic carbon, 14 C, extractable organic carbon, specific UV absorbance of K 2 SO 4 extractable organic C, total nitrogen, total dissolved nitrogen) along the soil profiles to examine associations between soil chemical and biological parameters. Even though root biomass was greater and rooting depth was deeper under switchgrass than soy, resource availability and microbial growth parameters were generally similar between vegetation types. Instead, depth significantly influenced soil chemical and biological parameters. For example, resource availability and total and relative active microbial biomass decreased with soil depth. Decreases in the relative active microbial biomass coincided with increased lag time (response time to external carbon inputs) along the soil profiles. Even at a depth of 210–240 cm, microbial communities were activated to grow by added resources within a day. Maximum specific growth rate decreased to a depth of 90 cm and then remained consistent in deeper layers. Our findings show that >10 years of vegetation and rooting depth changes may not be long enough to alter microbial growth parameters, and suggest that at least a portion of the microbial community in deep soils can grow rapidly in response to added resources. Our study determined microbial growth parameters that can be used in microbially-explicit models to simulate carbon dynamics in deep soil layers.

14C↗

DOE OSTI2023

Development of Mechanistic Fission Gas Release and Swelling Models for UN Fuels in BISON

This report describes the work in NEAMS (Nuclear Energy Advanced Modeling and Simulation) to develop a mechanistic fission gas model for uranium nitride fuels in BISON. The existing Sifgrs (Simple integrated fission gas release and swelling) model tracks the average properties of two bubble populations in the bulk and at the grain boundaries. It was recognized that dislocations play a crucial role in the fission gas swelling of UN, and an irradiation-induced dislocation density model was needed, as well as a model describing how fission gas interacts with dislocations creating a third population of bubbles along dislocations. A mechanistic model that tracks dislocation bubbles was implemented in Sifgrs. This model was used to simulate fission gas swelling and release in UN. Lower-length-scale calculations and experimental observations from carbide fuel were leveraged to help populate the model with essential parameters. As a placeholder, an empirical function was formulated for the evolution of the dislocation network. Because of the difference in evolution of this network at different temperatures, the dislocation bubbles are able to capture behavior that the bulk intragranular bubbles cannot. It was found that the bulk bubbles dominate microscopic swelling at low temperatures, and the dislocation bubbles dominate at higher temperatures. Based on this model, the transition between the two bubble types is the main factor behind the breakaway swelling phenomenon in UN. In order to test the model, a couple of assessments were run. For the lower temperature JOYO pins, the model produces reasonable fission gas release and swelling values. For the higher temperature SP1 pin, the model dramatically underestimates the fission gas release. This issue is attributed to the gas being trapped inside dislocation bubbles, unable to escape to grain boundaries to cause fission gas release, and could be remedied by a mechanistic dislocation model that allows the dislocation density to decrease at very high temperatures. A preliminary mechanistic dislocation model was developed supported by first-principle calculations. These calculations provided valuable insight into how interstitial defects cluster in UN in the {110} orientation, and may eventually form dislocation loops, leading to the conclusion that dislocation loops may nucleate from these clusters. An estimate of the dislocation line energy in UN was obtained and will be improved in future work. The free-energy cluster dynamics code Centipede was used to track interstitial and vacancy absorption at dislocation loops and calculate their growth. In addition, improvements to Centipede were made including new convergence criteria for transient simulations and sink driving force updates.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗