Search NASA⌕ Search

SEARCH · Search NASA

Results for “Code assessment”

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 775 records · Page 43

Validation of Ullage Collapse Due to Lateral Violent Slosh

Understanding the coupling between thermodynamics and propellant sloshing is important for cryogenic upper stages, especially during the lift-off phase and in missions including multiple restarts, as the fluid dynamic condition of the propellant is mainly influenced by various flight maneuvers. As a matter of fact, one of the major technical challenges associated with cryogenic fluid management (CFM) is the phenomenon of ullage gas collapse. This collapse is mainly caused by heat transfer from ullage gas to tank walls and interfacing propellants, which are both at temperatures well below those of this gas. The understanding of this phenomenon is of major importance concerning the next generation of cryogenic propellants. This section presents our assessment of a multi-phase CFD code developed at MSFC, Loci/STREAM-VOF, in predicting ullage thermodynamics and collapse due to violent slosh. The experimental data of Himeno et al. will be used for validation.

H. Q Yang↗

Navigating Team Dynamics: Automated Detection of Micro-Behaviors Between Team Members Through Longitudinal Interaction Data

The success in future long term space exploration missions will depend on the cooperation, coordination, and mutual understanding among the crew members. Micro-behaviors are momentary, subtle linguistic and paralinguistic indicators of thinking and feeling toward another member of the team (Cortina et al., 2001; Smith & Griffiths, 2022) that can significantly impact team dynamics and influence the overall team performance (Paromita & Chaspari, 2024). Due to their interactive nature, micro-behaviors have a sender (i.e., the team member expressing the micro-behavior) and a target (the team member impacted by the micro-behavior). Detection of these behaviors can assist in avoiding possible conflict among crew members and promoting the overall team success. Our prior research focused on an initial proof of concept of machine learning (ML) models and natural language processing (NLP) techniques that were used for automatically detect micro-behaviors among crew members of the US National Aeronautics and Space Administration’s (NASA) Human Exploration Research Analog (HERA) Campaigns 4 and 5 missions (Paromita et al., 2023). Results underscored the importance of incorporating contextual information in the ML models in the form of sentiment analysis, type of task, and dyadic interaction among team members. Here, we expand the scope of our prior work in two ways. First, we assess ML/NLP methods on new behavioral annotations coded using an adapted version of Smith & Griffins (2022) theoretical framework in terms of Violation (i.e., presence of valenced behavior, uplifting/positive or discouraging/negative), Intensity (i.e., force of behavior in terms of how uplifting or discouraging is the behavior), and Intent (i.e., motive of the behavior in terms of whether it was deliberate or unintentional). Second, we expand the design of the ML model to preserve information about the role of each team member within the occurrence of the micro-behavior (in contrast to the previous model that only considered the sender and the target without determining the team member role). This allows to consider all team members' contributions in the conversation and model long-term dependencies in the dialogue. Our experiments for this study are conducted on data from 5 teams of the NASA HERA C4 (NASA grant NNX16AQ48G (PI: Bell)). Conversations were extracted from the 1.5 hour Team Interaction Battery (TIB) task that occurred 5 times in-mission per crew. This resulted in a total of 13,058 conversational turns (i.e., 17.8% uplifting, 3.3% discouraging, 75.76% neutral, 3.14% nulls). Our findings with the revised behavioral coding and ML/NLP models indicate a 43.66% macro F1-score (i.e., 38.29% precision (P), 50.8% recall (R)) for a dialog state-tracking model that includes information from the sender only, and a 40.9% F1-score (i.e., 38.7% P, 43.36% R) for the same model that includes information from both the sender and the target of the micro-behavior. These are significantly higher compared to simple random forest models that classify behaviors strictly based on speech content and do not consider iterative team dynamics, achieving a 36.07% F1-score (i.e., 39.04% R, 33.53% P). Our findings demonstrate potential ways to leverage large conversational datasets to better capture complex team dynamics. We will discuss future directions including proposed models that can incorporate additional mission days and tasks beyond the TIB for objectively quantifying team behavior at high temporal resolution in space exploration missions.

Projna Paromita↗

SOARCA Peach Bottom Atomic Power Station Long-Term Station Blackout Uncertainty Analysis: Knowledge Advancement.

This paper describes the knowledge advancements from the uncertainty analysis for the State-of- the-Art Reactor Consequence Analyses (SOARCA) unmitigated long-term station blackout accident scenario at the Peach Bottom Atomic Power Station. This work assessed key MELCOR and MELCOR Accident Consequence Code System, Version 2 (MACCS2) modeling uncertainties in an integrated fashion to quantify the relative importance of each uncertain input on potential accident progression, radiological releases, and off-site consequences. This quantitative uncertainty analysis provides measures of the effects on consequences, of each of the selected uncertain parameters both individually and in interaction with other parameters. The results measure the model response (e.g., variance in the output) to uncertainty in the selected input. Investigation into the important uncertain parameters in turn yields insights into important phenomena for accident progression and off-site consequences. This uncertainty analysis confirmed the known importance of some parameters, such as failure rate of the Safety Relief Valve in accident progression modeling and the dry deposition velocity in off-site consequence modeling. The analysis also revealed some new insights, such as dependent effect of cesium chemical form for different accident progressions. (auth)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Measurements of impurity concentrations and transport in the Lithium Tokamak Experiment

This thesis presents new measurements of core impurity concentrations and transport in plasmas with lithium coatings on all-metal plasma facing components (PFCs) in the Lithium Tokamak Experiment (LTX). LTX is a modest-sized spherical tokamak uniquely capable of operating with large area solid and/or liquid lithium coatings essentially surrounding the entire plasma (as opposed to just the divertor or limiter region in other devices). Lithium (Li) wall-coatings have improved plasma performance and confinement in several tokamaks with carbon (C) PFCs, including the National Spherical Torus Experiment (NSTX). In NSTX, contamination of the core plasma with Li impurities was very low (<0.1%) despite extensive divertor coatings. Low Li levels in NSTX were found to be largely due to neoclassical forces from the high level of C impurities. Studying impurity levels and transport with Li coatings on stainless steel surfaces in LTX is relevant to future devices (including future enhancements to NSTX-Upgrade) with all-metal PFCs. The new measurements in this thesis were enabled by a refurbished Thomson scattering system and improved impurity spectroscopy, primarily using a novel visible spectrometer monitoring several Li, C, and oxygen (O) emission lines. A simple model was used to account for impurities in unmeasured charge states, assuming constant density in the plasma core and constant concentration in the edge. In discharges with solid Li coatings, volume averaged impurity concentrations were low but non-negligible, with~2-4% Li, ~0.6-2% C, ~0.4-0.7% O, and Z_eff<1.2. Transport was assessed using the TRANSP, NCLASS, and MIST codes. Collisions with the main H ions dominated the neoclassical impurity transport, unlike in NSTX, where collisions with C dominated. Furthermore, neoclassical transport coefficients calculated with NCLASS were similar across all impurity species and differed no more than a factor of two, in contrast to NSTX where they differed by an order of magnitude. However, time-independent simulations with MIST indicated that unlike NSTX, neoclassical theory did not fully capture the impurity transport and anomalous transport likely played a significant role in determining impurity profiles.

Lithium plasma facing components↗

Assessment of ALEGRA Computation for Magnetostatic Configurations

Here, a closed-form solution is described here for the equilibrium configurations of the magnetic field in a simple heterogeneous domain. This problem and its solution are used for rigorous assessment of the accuracy of the ALEGRA code in the quasistatic limit. By the equilibrium configuration we understand the static condition, or the stationary states without macroscopic current. The analysis includes quite a general class of 2D solutions for which a linear isotropic metallic matrix is placed inside a stationary magnetic field approaching a constant value H i ° at infinity. The process of evolution of the magnetic fields inside and outside the inclusion and the parameters for which the quasi-static approach provides for self-consistent results is also explored. Lastly, it is demonstrated that under spatial mesh refinement, ALEGRA converges to the analytic solution for the interior of the inclusion at the expected rate, for both body-fitted and regular rectangular meshes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

VALIDATION OF A RELAP5 POINT KINETICS MODEL OF TREAT

Development of advanced nuclear fuels requires testing of fuel specimens under various power-cooling mismatch conditions. The Transient Reactor Test (TREAT) facility at the Idaho National Laboratory will be used to provide transient nuclear heating of fuel specimens. Experiments will be conducted to support assessment of accident tolerant fuel designs being proposed for light water reactor applications. The experiments to be conducted at TREAT will provide data for the assessment of advanced, multi-physics computer codes. However, before the experiments can be performed, calculations are required to demonstrate that they will meet program objectives and that they can be performed safely in TREAT in the event of an unplanned reactivity insertion accident that occurs during the experiment. The safety calculations for the first experiments need to be performed in the next year and therefore will be performed with an existing computer code, RELAP5-3D. A point kinetics model of TREAT was developed for RELAP5-3D. The purpose of the model is to calculate the reactor power during experiments and to support the accident analyses of those experiments. The RELAP5-3D model was validated using data from nine TREAT experiments conducted previously. The RELAP5-3D model generates results that are in reasonable agreement with measured values of maximum core power, energy deposition, and maximum fuel temperature for a wide range of reactivity insertions.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

A variable-order time-dependent neutron transport method for nuclear reactor kinetics using analytically-integrated space-time characteristics

A new time-dependent neutron transport method based on the method of characteristics (MOC) has been developed. Whereas most spatial kinetics methods treat time dependence through temporal discretization, this new method treats time dependence by defining the characteristics to span space and time. In this implementation regions are defined in space-time where the thickness of the region in time fulfills an analogous role to the time step in discretized methods. The time dependence of the local source is approximated using a truncated Taylor series expansion with high order derivatives approximated using backward differences, permitting the solution of the resulting space-time characteristic equation. To avoid a drastic increase in computational expense and memory requirements due to solving many discrete characteristics in the space-time planes, the temporal variation of the boundary source is similarly approximated. This allows the characteristics in the space-time plane to be represented analytically rather than discretely, resulting in an algorithm comparable in implementation and expense to one that arises from conventional time integration techniques. Furthermore, by defining the boundary flux time derivative in terms of the preceding local source time derivative and boundary flux time derivative, the need to store angularly-dependent data is avoided without approximating the angular dependence of the angular flux time derivative. The accuracy of this method is assessed through implementation in the neutron transport code DeCART. The method is employed with variable-order local source representation to model a TWIGL transient. The results demonstrate that this method is accurate and more efficient than the discretized method. (authors)

97 MATHEMATICS AND COMPUTING↗

Opportunities for Artificial Intelligence in Emergency Management

This article for the Domestic Preparedness Journal highlights how PNNL's in-depth landscape assessment dug deep into research articles, preprints, code repositories, and surveys at the intersection of applied artificial intelligence and emergency management.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Palo Verde Water Cycle Model

SAND2022-12901 O The Palo Verde Water Cycle Model (PVWCM) is a dynamic multi-physics and economics code that allows managers and engineers to assess the economic benefit of potential water-saving technologies and water sources before money is invested in them. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Middleton, Bobby↗

A Benchmark Example for Delamination Propagation Predictions Based on the Single Leg Bending Specimen Under Quasi-Static and Fatigue Loading

Benchmark examples based on Single Leg Bending (SLB) specimens with equal and unequal bending arm thicknesses were used to assess the performance of delamination prediction capabilities in finite element codes. First, the development of the quasi-static benchmark cases using the Virtual Crack Closure Technique (VCCT) is discussed in detail. Second, based on the quasi-static benchmark results, additional benchmark cases to assess delamination propagation under fatigue loading are created. Third, the application is demonstrated for the commercial finite element code Abaqus Standard 2018. The benchmark cases are compared to results obtained from VCCT-based, automated quasi-static propagation analysis. A comparison with results from automated fatigue propagation analysis was not performed at this point since the current version of Abaqus does not include this capability under variable mixed-mode conditions. In general, good agreement between the results obtained from the quasi-static propagation analysis and the benchmark results were achieved. Overall, the benchmarking procedure proved valuable for analysis verification.

Krueger, Ronald↗

A Benchmark Example for Delamination Propagation Predictions Based on the Single Leg Bending Specimen under Quasi-static and Fatigue Loading

Benchmark examples based on Single Leg Bending (SLB) specimens with equal and unequal bending arm thicknesses were used to assess the performance of delamination prediction capabilities in finite element codes. First, the development of the quasi-static benchmark cases using the Virtual Crack Closure Technique (VCCT) is discussed in detail. Second, based on the quasi-static benchmark results, additional benchmark cases to assess delamination propagation under fatigue loading are created. Third, the application is demonstrated for the commercial finite element code Abaqus Standard 2018. The benchmark cases are compared to results obtained from VCCT-based, automated quasi-static propagation analysis. A comparison with results from automated fatigue propagation analysis was not performed at this point since the current version of Abaqus does not include this capability under variable mixed-mode conditions. In general, good agreement between the results obtained from the quasi-static propagation analysis and the benchmark results were achieved. Overall, the benchmarking procedure proved valuable for analysis verification.

Ronald Krueger↗

Wind tunnel requirements for computational fluid dynamics code verification

The role of experiment in the development of Computational Fluid Dynamics (CFD) for aerodynamic flow field prediction is discussed. Requirements for code verification from two sources that pace the development of CFD are described for: (1) development of adequate flow modeling, and (2) establishment of confidence in the use of CFD to predict complex flows. The types of data needed and their accuracy differs in detail and scope and leads to definite wind tunnel requirements. Examples of testing to assess and develop turbulence models, and to verify code development, are used to establish future wind tunnel testing requirements. Versatility, appropriate scale and speed range, accessibility for nonintrusive instrumentation, computerized data systems, and dedicated use for verification were among the more important requirements identified.

Marvin, Joseph G.↗

Wind tunnel requirements for computational fluid dynamics code verification

The role of experiment in the development of computational fluid dynamics (CFD) for aerodynamic flow field prediction is discussed. Requirements for code verification from two sources that pace the development of CFD are described for: (1) development of adequate flow modeling, and (2) establishment of confidence in the use of CFD to predict complex flows. The types of data needed and their accuracy differs in detail and scope and leads to definite wind tunnel requirements. Examples of testing to assess and develop turbulence models, and to verify code development, are used to establish future wind tunnel testing requirements. Versatility, appropriate scale and speed range, accessibility for nonintrusive instrumentation, computerized data systems, and dedicated use for verification were among the more important requirements identified.

Marvin, Joseph G.↗

Approaching the investigation of plasma turbulence through a rigorous verification and validation procedure: A practical example

In the present work, a Verification and Validation procedure is presented and applied showing, through a practical example, how it can contribute to advancing our physics understanding of plasma turbulence. Bridging the gap between plasma physics and other scientific domains, in particular, the computational fluid dynamics community, a rigorous methodology for the verification of a plasma simulation code is presented, based on the method of manufactured solutions. This methodology assesses that the model equations are correctly solved, within the order of accuracy of the numerical scheme. The technique to carry out a solution verification is described to provide a rigorous estimate of the uncertainty affecting the numerical results. A methodology for plasma turbulence code validation is also discussed, focusing on quantitative assessment of the agreement between experiments and simulations. The Verification and Validation methodology is then applied to the study of plasma turbulence in the basic plasma physics experiment TORPEX [Fasoli et al., Phys. Plasmas 13, 055902 (2006)], considering both two-dimensional and three-dimensional simulations carried out with the GBS code [Ricci et al., Plasma Phys. Controlled Fusion 54, 124047 (2012)]. The validation procedure allows progress in the understanding of the turbulent dynamics in TORPEX, by pinpointing the presence of a turbulent regime transition, due to the competition between the resistive and ideal interchange instabilities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Preparatory Verification of a New Monte Carlo Code cosRMC based on VERA Core Physics Benchmark

To promote nuclear power plant system design and simulation software self-reliance development, DOE of China approved a R and D key project to develop an integrated engine for pressurized water reactor core and system design and analysis, referred as COSINE hereafter. The scope of this code package includes reactor physics, thermal hydraulics, safety analysis, severe accidents, probability safety assessment and etc. Many universities, research institutes and nuclear related corporations in China are involved in the code development, among which the State Nuclear Power Software Development Center (SNPSDC) plays the leading role. With the requirement of accurate three-dimensional modeling of the advanced, special or new complex core in reactor neutronics and shielding analysis, Monte Carlo method is becoming a more powerful tool for reactor analysis. Thus, a new Monte Carlo transport code named cosRMC (abbreviation for COSINE Reactor Monte Carlo) has been developed jointly by SNPSDC and Department of Engineering Physics, Tsinghua University, the main developer of RMC code. The cosRMC now has many functions as criticality calculation, fixed-source calculation, burnup calculation, and etc. Some techniques for geometry treatment, new burnup algorithm, source convergence acceleration, massive tally and parallel calculation have been implemented in cosRMC to improve the efficiency. In order to ensure its reliability, several types of V and V work have been and are being carried out for cosRMC. In this paper, some preparatory verification results based on VERA Core Physics Benchmark are presented. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Validation of radiative transfer computation with Monte Carlo method for ultra-relativistic background flow

We developed a three-dimensional radiative transfer code for an ultra-relativistic background flow-field by using the Monte Carlo (MC) method in the context of gamma-ray burst (GRB) emission. For obtaining reliable simulation results in the coupled computation of MC radiation transport with relativistic hydrodynamics which can reproduce GRB emission, we validated radiative transfer computation in the ultra-relativistic regime and assessed the appropriate simulation conditions. The radiative transfer code was validated through two test calculations: (1) computing in different inertial frames and (2) computing in flow-fields with discontinuous and smeared shock fronts. The simulation results of the angular distribution and spectrum were compared among three different inertial frames and in good agreement with each other. If the time duration for updating the flow-field was sufficiently small to resolve a mean free path of a photon into ten steps, the results were thoroughly converged. The spectrum computed in the flow-field with a discontinuous shock front obeyed a power-law in frequency whose index was positive in the range from 1 to 10 MeV. The number of photons in the high-energy side decreased with the smeared shock front because the photons were less scattered immediately behind the shock wave due to the small electron number density. The large optical depth near the shock front was needed for obtaining high-energy photons through bulk Compton scattering. Even one-dimensional structure of the shock wave could affect the results of radiation transport computation. Although we examined the effect of the shock structure on the emitted spectrum with a large number of cells, it is hard to employ so many computational cells per dimension in multi-dimensional simulations. Therefore, a further investigation with a smaller number of cells is required for obtaining realistic high-energy photons with multi-dimensional computations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Studies of the charging of a thin dust layer in a plasma

The present study is designed to extend Wilson's model to include refinements in the actual operational design of the simulation code itself and to include additional physics in order to assess the resulting electric field structure and grain charges. Most of the initial phases of the work have focussed on the fact that the code, though in principle reliable, suffers from severe time constraints, requiring as must as 10-50 hours of CPU time for typical runs on a VAX computer. It is likely that by staging the code into steps in which early steps are carried out crudely and later steps successively fine tuned the total CPU time can be cut significantly. Preliminary results indicate the possibility to cut CPU time to as much as 1/3 or 1/5 the former time. The original Wilson code has been carefully examined, tested and documented. Shortcuts and staging systems have been introduced and tested which should make the code produce more accurate results with a finer grid system than would have been feasible heretofor. Continued testing and implementation of these refinements will be carried out at the author's home institution in the near future. In addition to the effort at reducing CPU time to a manageable amount, work is in progress to include a dust-size distribution instead of assuming all dust grains have the same size.

Peterson, Lennart R.↗

Scatter factors assessment in microbeam radiation therapy

Purpose: The success of the preclinical studies in Microbeam Radiation Therapy (MRT) paved the way to the clinical trials under preparation at the Biomedical Beamline of the European Synchrotron Radiation Facility. Within this framework, an accurate determination of the deposited dose is crucial. With that aim, the scatter factors, which translate the absolute dose measured in reference conditions (2 x 2 cm{sup 2} field size at 2 cm-depth in water) to peak doses, were assessed. Methods: Monte Carlo (MC) simulations were performed with two different widely used codes, PENELOPE and GEANT4, for the sake of safety. The scatter factors were obtained as the ratio of the doses that are deposited by a microbeam and by a field of reference size, at the reference depth. The calculated values were compared with the experimental data obtained by radiochromic (ISP HD-810) films and a PTW 34070 large area chamber. Results: The scatter factors for different microbeam field sizes assessed by the two MC codes were in agreement and reproduced the experimental data within uncertainty bars. Those correction factors were shown to be non-negligible for the future MRT clinical settings: an average 30% lower dose was deposited by a 50 {mu}m microbeam with respect to the reference conditions. Conclusions: For the first time, the scatter factors in MRT were systematically studied. They constitute an essential key to deposit accurate doses in the forthcoming clinical trials in MRT. The good agreement between the different calculations and the experimental data confirms the reliability of this challenging micrometric dose estimation.

62 RADIOLOGY AND NUCLEAR MEDICINE↗