Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reference model”

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 1,297 records · Page 72

Comparing Theoretically Scaled Biomechanical Models

BACKGROUND An investigation into incorporating space suit aspects required a 50th-percentile male model which was not part of the existing dataset of models. Previously, theoretical models for 5th-percentile female and 95th-percentile male were created utilizing a scaled test subject as close as possible to the target height and weight. Scaling factors were created by taking the height ratio and applying it uniformly to the model body segments, then fine tuning to ensure the theoretical model’s height is as expected. This study was initiated to examine the existing method of isometrically scaling in OpenSim [1,2] and to create alternative methods which do not rely on an existing subject being close in height and weight to the theoretical model of interest. Generating a theoretical biomechanical model provides additional abilities without the reliance on available OpenSim models or real subjects. Gained abilities include creating different percentile models and attaining representative anthropometry for any target height/weight, such as targeting specific crew populations or gaps within the current dataset. METHODS AND RESULTS This study includes 3 methods for generating the theoretical model; utilizing the modified unscaled OpenSim Full Body Rajagopal Model (FBRM) [3,4], a dataset of scaled OpenSim models, and existing scaling factors from Dumas et al. [5]. All methods use the Anthropometric Survey of US Army Personnel (ANSUR II) [6] as an input of necessary anthropometric measurements. The following measurements are retrieved from the ANSUR II collection: mass, stature, cervicale height, acromial height, axilla height, waist height, trochanterion height, lateral femoral epicondyle height, lateral malleolus height, acromion-radiale length, radiale-stylion length, palm length, ball of foot length, bicristal breadth, bimalleolar breadth. Some of the measurements are direct segment lengths and others are utilized to derive segment lengths. Height, weight, and age ranges are the required inputs to parse the collection for a mean value of the measurements. A simple iterative process may be necessary for the output of mean mass and stature consistent with the theoretical model of interest. In some cases, one may specify exact values instead of a range but that is dependent upon whether those exact values pertain to a single subject from the ANSUR II collection. Equations from Dumas et al. were used to calculate scaling factors for the two methods that utilize OpenSim scaled and unscaled models. In the third method, Dumas’ scaling factors were used directly instead of generating our own. After scaling factors are achieved, they are applied along with the mass and segment lengths retrieved from ANSUR II in order to calculate the segment mass, Center of Mass (CoM), mass Moment of Inertia (MoI), and the joint location in the parent frame. The theoretical models from the different methods were compared to the previous existing method of isometrically scaling in OpenSim, as well as the standards found in NASA STD-3000 [7] and the NASA Human Integration Design Handbook [8]. The whole-body center of mass is the main comparison performed between the models and the standards. The body segment mass properties were also compared when possible. In some cases, the standards do not provide complete data, or the center of mass location is not provided with respect to the appropriate joint coordinate system. For a 50th-percentile male subject, results from isometric scaling in OpenSim and the method utilizing the unscaled FBRM were compared to the standards in NASA STD-3000. According to the standard, the whole-body vertical CoM location was taken from the head vertex to the CoM. The standards provide 80.2 cm for the CoM vertical component whereas isometric scaling in OpenSim and utilizing the unscaled FBRM provide 82.2 and 82.3 cm, respectively. Through definition, the whole-body CoM in the medial/lateral direction agreed across the models. Investigation is on-going to clearly identify the reference point pertaining to dorsal/ventral whole-body CoM position in standards and relate that to the theoretical models. REFERENCES [1] Delp, S.L., et al., “OpenSim: Open-source Software to Create and Analyze Dynamic Simulations of Movement”, IEEE Transactions on Biomedical Engineering, (2007). [2] Seth, A., Hicks J.L., Uchida, T.K., Habib, A., Dembia, C.L., Dunne, J.J., Ong, C.F., DeMers, M.S., Rajagopal, A., Millard, M., Hamner, S.R., Arnold, E.M., Yong, J.R., Lakshmikanth, S.K., Sherman,n M.A., Delp, S.L. OpenSim: Simulating musculoskeletal dynamics and neuromuscular control to study human and animal movement. Plos Computational Biology, 14(7). (2018) [3] R. K. Huffman, W. K. Thompson, C. A. Gallo, L. J. Quiocho, “Improvement of Scaling and Inverse Kinematic Results with Additional Upper Body Joints Added to Opensim Rajagopal Model”, NASA Human Research Program Investigator’s Workshop, (2019). [4] Huffman R.K., Thompson W., Gallo C., “Modified OpenSim Rajagopal Full Body Model”, New Technology Report, MSC-26872-1. (August 2020). [5] Dumas R., Cheze L. and Verriest J., 2007. “Adjustments to McConville et al. and Young et al. body segment inertial parameters”. Journal of Biomechanics 40, pp. 543-553 [6] Gordon C.C., 2012 “Anthropometric Survey of U.S. Army Personnel: Methods and Summary Statistics”, US Army Natick Soldier RD&E Center. [7] NASA STD-3000/REV-B, 1995. “The Man-System Integration Standards” [8] NASA/SP-2010-3407/REV-1, 2014. “The Human Integration Design Handbook (HIDH)”

F N Matari↗

Structure and Dynamics of the Quiet Solar Chromosphere

The grant supported research on the structure of the quiet, nonmagnetic chromosphere and on wave excitation and propagation in both the nonmagnetic chromosphere and the magnetic network. The work on the structure of the chromosphere culminated in the recognition that between two competing views of the solar chromosphere, older models by Avrett and collaborators (referred to as VAL) and the newer, dynamical model by Carlsson & Stein (referred to as CS), the clear decision is in favor of the older models, and this in spite of the evident lack of physics, which does not include wave motion and oscillations. The contrast between the static VAL models and the dynamical CS model can be stated most succinctly by comparing the temperature variation implied by the VAL models and the temperature fluctuations of the CS model, which are, respectively, of the order of 10% for the VAL model (at heights where hydrogen is 50% ionized) and a factor of 10 (at the upper boundary of their chromospheric model). The huge fluctuations of the CS model have never been observed, whereas the smaller temperature variations of the VAL models are consistent with ground-based and space-based observations. While it should be obvious which model describes the Sun and which one fails, the case is far from settled in the minds of solar physicists. Thus, much educational work remains to be done and, of course, more research to develop arguments that make the case more convincing. The research on waves and oscillations has been based on a unified theory of excitation of acoustic waves in the field-free atmosphere and of transverse and longitudinal waves in magnetic flux tubes located in the magnetic network by noting, first, that impulsive excitation of all these waves in gravitationally stratified media leads to oscillations at the respective cutoff frequencies and, second, that the observed oscillation frequencies in the nonmagnetic and magnetic parts of the chromosphere match corresponding cutoff frequencies in the upper solar photosphere. The dynamical simulations by Carlsson and Stein have been most instructive and of fundamental importance for understanding wave propagation in a stratified medium by their "flaws", the most important of which is an intensity excess at the H2v emission peak in the H line of Ca II, which is surprising since the observed intensity should have been an upper limit to the simulated intensity. The only plausible explanation for a predicted intensity that is higher than observed is that energy is spread horizontally in upward propagation in the Sun, but not in the plane-wave modeling of CS (as well as by almost everyone else). Investigations of the horizontal size of the region disturbed by the upward-propagating shock in the acoustic-wave propagation implies that the waves in H2v bright-point oscillations emanate from a point source with a diameter corresponding to the width of an intergranular lane, about 100 km, and reach a size of about 4000 km in the upper layers of the chromosphere. Linear, analytic modeling of waves emanating from a point source in a stratified atmosphere shows that the upward-expanding propagation channel does not have sharp boundaries and that the shape of the wave front depends on the order of the wave behind the initial pulse. Otherwise, the behavior of the linear waves resembles that of the nonlinear shock waves observed in the Sun. Research that needs to be done to firm up the conclusions reached above concerns the numerical simulations of nonlinear waves and oscillations in a three-dimensional stratified atmosphere with impulsive excitation, and observations linking directly the horizontal size of the disturbed area in upward propagation to individual waves.

Kalkofen, Wolfgang↗

A coarse-grained deep neural network model for liquid water

We introduce a coarse-grained deep neural network (CG-DNN) model for liquid water that utilizes 50 rotational and translational invariant coordinates and is trained exclusively against energies of ~30 000 bulk water configurations. Our CG-DNN potential accurately predicts both the energies and the molecular forces of water, within 0.9 meV/molecule and 54 meV/Å of a reference (coarse-grained bond-order potential) model. The CG-DNN water model also provides good prediction of several structural, thermodynamic, and temperature dependent properties of liquid water, with values close to those obtained from the reference model. More importantly, CG-DNN captures the well-known density anomaly of liquid water observed in experiments. Our work lays the groundwork for a scheme where existing empirical water models can be utilized to develop a fully flexible neural network framework that can subsequently be trained against sparse data from high-fidelity albeit expensive beyond-DFT calculations.

42 ENGINEERING↗

FOCAL Campaign I: Advanced Wind Turbine Control Strategies

Campaign I of the Floating Offshore-wind Controls Advanced Laboratory Experimental Program (FOCAL) aims to generate a dataset enabling the validation of aerodynamic performance of a scaled turbine mounted on a rigid tower in a fixed condition. The turbine considered in the FOCAL testing campaigns is the IEA-Wind 15MW Reference Wind Turbine. This scaled model is capable of simulating advanced blade-pitch control strategies in a high-quality wind field. The turbine is fully instrumented to record a variety of parameters in real time such as structural loads and dynamics. The test data considered was generated at the University of Maine's Harold Alfond Wind and Wave (W2) testing facility. The Load Cases (LC) considered in this testing campaign are as follows: LC 1.X - Constant wind and blade-pitch with varying rotor speeds LC 2.X - Constant wind and rotor speed with varying blade-pitch LC 3.X - Varying wind with active closed-loop control Detailed properties on the modeled system are found in the following reference: Lenfest E., Floating Offshore-wind Controls Advanced Laboratory (FOCAL) Experimental Program - Campaign I: 1:70 Model-scale Testing of the IEA-Wind 15MW Reference Turbine. UMaine ASCC Report Number 23-40-1183. Details on the results of the verification and validation are found in the following reference: Mendoza, Nicole et al., "Verification and Validation of Model-Scale Turbine Performance and Control for the IEA Wind 15 MW Reference Wind Turbine," Energies, vol. 15, no. 20, 2022, https://doi.org/10.3390/en15207649.

17 WIND ENERGY↗

Model data for a watershed-scale study in the Portage River Basin (OH) examining the effects of subsurface drainage on the hydrologic response of an agricultural watershed.

This study builds on Rathore et al. (2024, WRR) and investigates the role of artificial tile-drainage on various aspects of watershed hydrological response, with a particular focus on peakflow. The model-data for the original modeling-focused paper (Rathore et al., 2024, WRR) is archived at Rathore et al. (2024, ESS-DIVE). Hence, this model-data archive provides scripts that are specific to this study that includes model updates, processing and analysis scripts. For details and models files of original model, readers are referred to Rathore et al. (2024, ESS-DIVE). The key difference between the model configuration in this study and Rathore et al. (2024, WRR) is that the tile drains are applied to the entire domain, to study the impact of tile-drains on different aspects of hydrological response. Additional scenario considering intensified precipitation after a dry period was also simulated. The Watershed Workflow package is implemented in Python3. The Jupyter notebooks can be executed through multiple open-source tools, for example, Anaconda Jupyter Lab, VS Studio Code, etc. Other data files include CSV and HDF5 files, which can be read through Python scripts.

54 ENVIRONMENTAL SCIENCES↗

Optimal estimation of large structure model errors

In-flight estimation of large structure model errors is usually required as a means of detecting inevitable deficiencies in large structure controller/estimator models. The present paper deals with a least-squares formulation which seeks to minimize a quadratic functional of the model errors. The properties of these error estimates are analyzed. It is shown that an arbitrary model error can be decomposed as the sum of two components that are orthogonal in a suitably defined function space. Relations between true and estimated errors are defined. The estimates are found to be approximations that retain many of the significant dynamics of the true model errors. Current efforts are directed toward application of the analytical results to a reference large structure model.

Rodriguez, G.↗

EUVAC: A solar EUV flux model for aeronomic calculations

This paper presents a new solar Extreme Ultraviolet (EUV) flux model for aeronomic calulations (EUVAC), which is based on the measured F74113 solar EUV reference spectrum. The model provides fluxes in the 37 wavelength bins that are in widespread use. This paper also presents cross sections to be used with the EUVAC flux model to calculate photoionization rates. The flux scaling for solar activity is accomplished using a proxy-based on the F10.7 index and its 81-day average together with the measured solar flux variation from the EUVS instrument on the Atmosphere Explorer E satellite. This new model produces 50-575 A integrated EUV fluxes in good agreement with rocket observations. The solar cycle variation of the chromospheric fluxes agrees well with the measured variation of the Lyman alpha flux between 1982 and 1988. In addition, the theoretical photoelectron fluxes, calculated using the new EUV flux model, are in good agreement with the solar minimum photoelectron fluxes from the Atmosphere Explorer E satellite and also with the solar maximum photoelectron fluxes from the Dynamics Explorer satellite. Its relative simplicity coupled with its ability to reproduce the 50-575 A solar EUV flux as well as the measured photoelectron spectrum makes the model well suited for aeronomic applications. However, EUVAC is not designed to accurately predict the solar flux variability for numerous individual lines.

Richards, P. G.↗

F2 peak electron density at Millstone Hill and Hobart: Comparsion of theory and measurement at solar maximum

This paper compares the observed behavior of the (F2) layer of the ionosphere at Millstone Hill and Hobart with calculations from the field line interhemispheric plasma (FLIP) model for solar maximum, solstice conditions in 1990. During the study period the daily F(sub 10.7) index varied by more than a factor of 2 (123 to 280), but the 81-day mean F(sub 10.7) (F(sub 10.7 A)) was almost constant near 190. Calculations were performed with and without the effects of vibrationally excited N2 (N(sup *)(sub 2) which affects the loss rate of atomic oxygen ions. In the case without N(sup *)(sub 2) there is generally good agreement between the model and measurement for the daytime, peak density of the F region (NmF2). Both the model and the measurement show a strong seasonal anomaly with the winter noon densities a factor of 3 to 4 greater than the summer noon densities at Millstone Hill and a factor of 2 greater at Hobart. The seasonal anomaly in the model is caused by changes in the neutral composition as given by the mass spectrometer and incoherent scatter (MSIS) 86 neutral density model. There is generally little or no increase in the observed noon NmF2 as a function of daily F(sub 10.7) except at Millstone Hill in winter. In contrast to the generally good agreement between model and data at noon, the model badly underestimates the density at night at Millstone Hill at all seasons. At Hobart the model reproduces the nighttime density variations well in both winter and summer. The international reference ionosphere (IRI) model generally provides a good representation of the average behavior of noon NmF2 and hmF2 but because the data show a lot of day-to-day variability, there are often large differences. The FLIP model is able to reproduce this variability when hmF2 is specified. The IRI model peak densities are better than the FLIP densities at night, but the IRI model does not represent the Millstone Hill summer data very well at night in 1990.

Richards, P. G.↗

Are Terrestrial Biosphere Models Fit for Simulating the Global Land Carbon Sink?

The Global Carbon Project estimates that the terrestrial biosphere has absorbed about one-third of anthropogenic CO 2 emissions during the 1959–2019 period. This sink-estimate is produced by an ensemble of terrestrial biosphere models and is consistent with the land uptake inferred from the residual of emissions and ocean uptake. The purpose of our study is to understand how well terrestrial biosphere models reproduce the processes that drive the terrestrial carbon sink. One challenge is to decide what level of agreement between model output and observation-based reference data is adequate considering that reference data are prone to uncertainties. To define such a level of agreement, we compute benchmark scores that quantify the similarity between independently derived reference data sets using multiple statistical metrics. Models are considered to perform well if their model scores reach benchmark scores. Our results show that reference data can differ considerably, causing benchmark scores to be low. Model scores are often of similar magnitude as benchmark scores, implying that model performance is reasonable given how different reference data are. While model performance is encouraging, ample potential for improvements remains, including a reduction in a positive leaf area index bias, improved representations of processes that govern soil organic carbon in high latitudes, and an assessment of causes that drive the inter-model spread of gross primary productivity in boreal regions and humid tropics. The success of future model development will increasingly depend on our capacity to reduce and account for observational uncertainties.

58 GEOSCIENCES↗

Harada–Tsutsui gauge recovery procedure: From Abelian gauge anomalies to the Stueckelberg mechanism

Revisiting a path-integral procedure developed by Harada and Tsutsui for recovering gauge invariance from anomalous effective actions, it is shown that there are two ways to achieve gauge symmetry: one already presented by the authors, which is shown to preserve the anomaly in the sense of standard current conservation law, and another one which is anomaly-free, preserving current conservation. It is also shown that the application of the Harada–Tsutsui technique to other models which are not anomalous but do not exhibit gauge invariance allows the identification of the gauge invariant formulation of the Proca model, also done by the referred authors, with the Stueckelberg model, leading to the interpretation of the gauge invariant map as a generalization of the Stueckelberg mechanism. -- Highlights: • A gauge restoration technique from Abelian anomalous models is discussed. • It is shown that there is another way that leads to gauge symmetry restoration from such technique. • It is shown that the first gauge restoration preserves the anomaly, while the proposed second one is free from anomalies. • It is shown that the proposed gauge symmetry restoration can be identified with the Stueckelberg mechanism.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Subseasonal Tropical Convection Characteristics in the Energy Exascale Earth System Model Version 2

Accurate simulation of subseasonal tropical moist convection remains a key challenge for Earth system models. The difficulties stem from the reliance of cumulus cloud processes on model parameterizations and the need to represent the multiscale nature of interactions among clouds, radiation, moisture, circulation, and surface energy fluxes. Equatorial convection drives circulation anomalies that can affect weather patterns and extremes globally, motivating efforts to better understand and simulate these tropical disturbances. Here, a detailed review of subseasonal tropical convective behavior as simulated in the Energy Exascale Earth System Model version 2 (E3SMv2) is presented, with comparison to its predecessor version 1 (E3SMv1) and reference data sets. Model structural changes to the deep convective trigger function and surface fluxes, along with parametric tuning of the cloud and microphysics schemes, together result in an improved depiction of organized tropical convection across scales. In particular, E3SMv2 exhibits a more realistic Madden‐Julian oscillation (MJO) and low‐frequency Kelvin waves—owing to a sharper time mean equatorial meridional moisture gradient and improved convection‐circulation coupling —as well as a better depiction of MJO Northern Hemisphere teleconnections. Despite these improvements, subseasonal precipitation variance continues to be strongly underestimated in E3SMv2. Use of a cloud plume model also reveals that the coupling between daily averaged tropical precipitation and lower tropospheric instability in E3SM is inconsistent with observations, a bias that could potentially impact the simulation of intraseasonal disturbances.

54 ENVIRONMENTAL SCIENCES↗

Extensions to Dynamic System Simulation of Fissile Solution Systems

Previous reports have documented the results of applying dynamic system simulation (DSS) techniques to model a variety of fissile solution systems. The SUPO (Super Power) aqueous homogeneous reactor (AHR) was chosen as the benchmark for comparison of model results to experimental data for steadystate operation.1 Subsequently, DSS was applied to additional AHR to verify results obtained for SUPO and extend modeling to prompt critical excursions, ramp reactivity insertions of various magnitudes and rate, and boiling operations in SILENE and KEWB (Kinetic Experiment Water Boiler).2 Additional models for pressurized cores (HRE: Homogeneous Reactor Experiment), annular core geometries, and accelerator-driven subcritical systems (ADAHR) were developed and results reported.3 The focus of each of these models is core dynamics; neutron kinetics, thermal hydraulics, radiolytic gas generation and transport are coupled to examine the time-based evolution of these systems from start-up through transition to steady-state. A common characteristic of these models is the assumption that (a) core cooling system inlet temperature and flow and (b) plenum gas inlet pressure and flow are held constant; no external (to core) component operations that may result in dynamic change to these parameters are considered. This report discusses extension of models to include explicit reference to cooling structures and radiolytic gas handling. The accelerator-driven subcritical generic system model described in References 3 and 4 is used as a basis for this extension.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

ASME Code change proposal to implement new universal high temperature constitutive models for Section III, Division 5

This report completes work on a universal high temperature constitutive model suitable for use with the ASME Boiler & Pressure Vessel Code Section III, Division 5 rules for the design by inelastic analysis of Class A nuclear reactor components. The goals of this work are to provide a simple model form that adequately captures the high temperature response of materials and can be applied to any future Code material. Additionally, the report describes an automated process for calibrating a model against test data. The idea is to simplify the effort required to generate a constitutive model for an arbitrary material, provided test data is available. This will accelerate the process of qualifying new Code materials in the future. In addition, the report provides calibrated models and detailed validation comparisons to test data for five currently-qualified or soon-to-be qualified materials: 316H, Grade 91, Alloy 800H, Alloy 617, and Alloy 709. The report surveys the available data for the remaining two ASME Class Materials --- 2.25Cr-1Mo and 304H --- concluding that there is enough data data to generate a model for 2.25Cr-1Mo steel provided some additional sources of non-public data can be included in the test database, but that a dedicated cyclic test campaign would be needed for 304H. Supplemental material includes the full text of an ASME Code change proposal to incorporate the models for the four currently-qualified Class A material, detailed validation comparisons to test data for the five material models, and input files for reference implementations of the constitutive models in the NEML and NEML2 modeling frameworks.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Description of the mean behaviour of ionospheric plasma temperatures

A status report on the empirical modeling of ionospheric electron and ion temperatures is given with special emphasis on the models used in the International Reference Ionosphere (IRI). Electron temperature models have now reached a state where reliable prediction of the mean altitudinal, latitudinal and diurnal variations is possible. These models are largely based on satellite measurements, but comparisons with incoherent scatter radar measurements have shown excellent agreement. Variations with season and magnetic and solar activity seem to be small and are not yet included consistently in these models. Similar to the electron temperature, the ion temperature shows the largest variations with altitude, latitude and local time. But due to the larger mass, these variations are smoother and more steady in the case of the ions and therefore easier to model. Nevertheless, very few ion temperature models exist. The IRI model takes advantage of the observed concurrence of the ion temperature with the neutral temperature at low altitudes and with the electron temperature at high altitudes.

Bilitza, Dieter↗

An adaptive identification and control scheme for large space structures

A unified identification and control scheme capable of achieving space at form performance objectives under nominal or failure conditions is described. Preliminary results are also presented, showing that the methodology offers much promise for effective robust control of large space structures. The control method is a multivariable, adaptive, output predictive controller called Model Predictive Control (MPC). MPC uses a state space model and input reference trajectories of set or tracking points to adaptively generate optimum commands. For a fixed model, MPC processes commands with great efficiency, and is also highly robust. A key feature of MPC is its ability to control either nonminimum phase or open loop unstable systems. As an output controller, MPC does not explicitly require full state feedback, as do most multivariable (e.g., Linear Quadratic) methods. Its features are very useful in LSS operations, as they allow non-collocated actuators and sensors. The identification scheme is based on canonical variate analysis (CVA) of input and output data. The CVA technique is particularly suited for the measurement and identification of structural dynamic processes - that is, unsteady transient or dynamically interacting processes such as between aerodynamics and structural deformation - from short, noisy data. CVA is structured so that the identification can be done in real or near real time, using computationally stable algorithms. Modeling LSS dynamics in 1-g laboratories has always been a major impediment not only to understanding their behavior in orbit, but also to controlling it. In cases where the theoretical model is not confirmed, current methods provide few clues concerning additional dynamical relationships that are not included in the theoretical models. CVA needs no a priori model data, or structure; all statistically significant dynamical states are determined using natural, entropy-based methods. Heretofore, a major limitation in applying adaptive controllers to LSS problems has been the requirement for accurate models. We feel that MPC coupled with the CVA identification scheme will resolve this limitation. The results are based on research performed under an SBIR grant, administered by the Air Force Office of Scientific Research.

Carroll, J. V.↗

TERRA Battery Thermal Control Anomaly - Simulation and Corrective Actions

The TERRA spacecraft was launched in December 1999 from Vandenberg Air Force Base, becoming the flagship of NASA's Earth Observing System program to gather data on how the planet's processes create climate. Originally planned as a 5 year mission, it still provides valuable science data after nearly 10 years on orbit. On October 13th, 2009 at 16:23z following a routine inclination maneuver, TERRA experienced a battery cell failure and a simultaneous failure of several battery heater control circuits used to maintain cell temperatures and gradients within the battery. With several cells nearing the minimum survival temperature, preventing the electrolyte from freezing was the first priority. After several reset attempts and power cycling of the control electronics failed to reestablish control authority on the primary side of the controller, it was switched to the redundant side, but anomalous performance again prevented full heater control of the battery cells. As the investigation into the cause of the anomaly and corrective action continued, a battery thermal model was developed to be used in determining the control ability remaining and to simulate and assess corrective actions. Although no thermal model or detailed reference data of the battery was available, sufficient information was found to allow a simplified model to be constructed, correlated against pre-anomaly telemetry, and used to simulate the thermal behavior at several points after the anomaly. It was then used to simulate subsequent corrective actions to assess their impact on cell temperatures. This paper describes the rapid development of this thermal model, including correlation to flight data before and after the anomaly., along with a comparative assessment of the analysis results used to interpret the telemetry to determine the extent of damage to the thermal control hardware, with near-term corrective actions and long-term operations plan to overcome the anomaly.

Grob, Eric W.↗

Using Open Standards and NASA Open Source Simulation Tools to Model Artemis Base Camp Mission Timelines

The United States’ National Aeronautics and Space Administration (NASA) has announced that the Artemis Program will return humans to the Moon, establishing a persistent presence with the Artemis Base Camp (ABC), and extend human exploration to Mars. The NASA Exploration Systems Simulations (NExSyS) team at NASA’s Johnson Space Center is using internationally developed simulation interoperability standards and NASA open source simulation tools to support Artemis concept, analysis, designs, development, training, and ultimately operations. The NExSyS team has been tasked to support early ABC architecture and mission analysis using mission time lines developed by the crew operations mission planning team. The NExSyS team is developing a distributed simulation framework with initial Artemis element implementations to model the ABC mission timelines using the international simulation interoperability standard High Level Architecture (HLA), the Simulation Interoperability Standards Organization’s Space Reference Federation Object Model (SpaceFOM), the NASA open source Trick Simulation Environment, and another NASA open source interface package called TrickHLA. The ABC architecture is composed of a number of key surface elements and resources. Some examples of modeled elements (also known as entities) are landers, habitats, rovers, logistics carriers, and astronauts. Some examples of modeled transferable and consumable resources are power, water, oxygen, nitrogen, scientific samples, and food. These entities and resources are modeled in a collection of individual simulations called Federates. A coordinated collection of interoperable federates is called a Federation and when these federates are tied together in a coordinated simulation run, it is referred to as a Federation Execution. The federates communicate through HLA using data exchange formats defined by a collection of machine readable files called Federation Object Models (FOMs). These FOM files are based on extensions to the SpaceFOM. This enables the instantiation and sharing of objects and interactions between federates in the federation. These provide for entity and resource tracking, object transfer, and data collection. Federate interactions are used to trigger events and notify federates of entity or resource transfers. For the initial implementation, the constituent federates are Trick-based simulations that use TrickHLA to provide the required HLA-base interoperability. These Trick-based simulations provide the required modeling for the individual Artemis elements along with the associated element resources. These federates provide a means to explore traverses between surface elements and exploration sites as scheduled in a mission timeline and explore the affects traverse times have on the overall mission timeline. The mission time lines are modeled using a Trick input file event handling capabilities. Each timeline operation is handled as individual simulation events, and triggered based on previous event status, time of operation, and simulated task completions. In addition, the ABC Federation can be used to perform Monte Carlo analysis. The Monte Carlo tool can vary the inputs, timings, and malfunctions to show how various contingencies in the mission can affect the mission timeline.

Keaton Craig Dodd↗

Bias Estimates Used in Lieu of Validation of Fission Products and Minor Actinides in MCNP K eff Calculations for PWR Burnup Credit Casks

The U.S. Nuclear Regulatory Commission (NRC) Division of Spent Fuel Storage and Transportation recently issued Interim Staff Guidance (ISG) 8, Revision 3. This ISG provides guidance for burnup credit (BUC) analyses supporting transport and storage of PWR pressurized water reactor (PWR) fuel in casks. Revision 3 includes guidance for addressing validation of criticality (k eff ) calculations crediting the presence of a limited set of fission products and minor actinides (FP&MA). Based on previous work documented in NUREG/CR-7109, recommendation 4 of ISG-8, Rev. 3, includes a recommendation to use 1.5 or 3% of the FP&MA worth to conservatively cover the bias due to the specified FP&MAs. This bias is supplementary to the bias and bias uncertainty resulting from validation of k eff calculations for the major actinides in SNF and does not address extension to actinides and fission products beyond those identified herein. The work described in this report involves comparison of FP&MA worths calculated using SCALE and MCNP with ENDF/B-V, -VI, and -VII based nuclear data and supports use of the 1.5% FP&MA worth bias when either SCALE or MCNP codes are used for criticality calculations, provided the other conditions of the recommendation 4 are met. The method used in this report may also be applied to demonstrate the applicability of the 1.5% FP&MA worth bias to other codes using ENDF/B V, VI or VII based nuclear data. The method involves use of the applicant s computational method to generate FP&MA worths for a reference SNF cask model using specified spent fuel compositions. The applicant s FP&MA worths are then compared to reference values provided in this report. The applicants FP&MA worths should not exceed the reference results by more than 1.5% of the reference FP&MA worths.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗