Search NASA⌕ Search

SEARCH · Search NASA

Results for “code comparison study”

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 415 records · Page 23

An assessment of space shuttle flight software development processes

In early 1991, the National Aeronautics and Space Administration's (NASA's) Office of Space Flight commissioned the Aeronautics and Space Engineering Board (ASEB) of the National Research Council (NRC) to investigate the adequacy of the current process by which NASA develops and verifies changes and updates to the Space Shuttle flight software. The Committee for Review of Oversight Mechanisms for Space Shuttle Flight Software Processes was convened in Jan. 1992 to accomplish the following tasks: (1) review the entire flight software development process from the initial requirements definition phase to final implementation, including object code build and final machine loading; (2) review and critique NASA's independent verification and validation process and mechanisms, including NASA's established software development and testing standards; (3) determine the acceptability and adequacy of the complete flight software development process, including the embedded validation and verification processes through comparison with (1) generally accepted industry practices, and (2) generally accepted Department of Defense and/or other government practices (comparing NASA's program with organizations and projects having similar volumes of software development, software maturity, complexity, criticality, lines of code, and national standards); (4) consider whether independent verification and validation should continue. An overview of the study, independent verification and validation of critical software, and the Space Shuttle flight software development process are addressed. Findings and recommendations are presented.

Source record↗

Design of Experiments for Both Experimental and Analytical Study of Exhaust Plume Effects on Sonic Boom

Computational fluid dynamics (CFD) analysis has been performed to study the plume effects on sonic boom signature for isolated nozzle configurations. The objectives of these analyses were to provide comparison to past work using modern CFD analysis tools, to investigate the differences of high aspect ratio nozzles to circular (axisymmetric) nozzles, and to report the effects of under expanded nozzle operation on boom signature. CFD analysis was used to address the plume effects on sonic boom signature from a baseline exhaust nozzle. Nearfield pressure signatures were collected for nozzle pressure ratios (NPRs) between 6 and 10. A computer code was used to extrapolate these signatures to a ground-observed sonic boom N-wave. Trends show that there is a reduction in sonic boom N-wave signature as NPR is increased from 6 to 10. As low boom designs are developed and improved, there will be a need for understanding the interaction between the aircraft boat tail shocks and the exhaust nozzle plume. These CFD analyses will provide a baseline study for future analysis efforts. For further study, a design of experiments has been conducted to develop a hybrid method where both CFD and small scale wind tunnel testing will validate the observed trends. The CFD and testing will be used to screen a number of factors which are important to low boom propulsion integration, including boat tail angle, nozzle geometry, and the effect of spacing and stagger on nozzle pairs. To design the wind tunnel experiment, CFD was instrumental in developing a model which would provide adequate space to observe the nozzle and boat tail shock structure without interference from the wind tunnel walls.

Castner, Raymond S.↗

How does negative triangularity mitigate ITG turbulence and transport?

Improved confinement in negative triangularity (NT) experiments is attributed to reduced fluxes driven by micro-turbulence. The physical mechanism of why thermal confinement improves in NT relative to PT is unknown. This study employs gyrokinetic flux tube simulations using the GENE code with local Miller equilibrium to elucidate the physical mechanisms behind the beneficial effects of NT flux surface shapes. The focus is on collisionless ion temperature gradient (ITG) driven turbulence with adiabatic electrons. The kinetic profiles are held fixed across a scan of triangularity values, thus enabling comparisons on a level playing field. The reduced linear growth rates for NT is shown to be due to a reduced eigenmode averaged magnetic drift frequency and a wider, stronger negative local magnetic shear region about the outboard mid-plane. The nonlinear heat flux is lower for NT than that for PT, due to reduced radial correlation length and increased correlation time (τ c ) of fluctuations. These, in turn, are due to a comparatively higher level of self-generated zero-frequency E × B zonal shearing rate ω E in NT as compared to PT. Though the linear zonal potential residual is lower for NT, the nonlinearly generated E × B zonal shearing rate is higher for NT than for PT. This outcome is linked to the distinctive features of the radial wavenumber spectra of the zonal potential and the zonal shearing rate. The dimensionless parameter ω E τ c is suggested as a figure of merit. This is higher for NT than for PT. Thus, the reduced heat diffusivity for NT is linked to increased ω E τ c . Self-generated temperature corrugations (i.e. zonal temperature gradients) are much weaker than the background mean temperature gradient. Nevertheless, temperature corrugations are more pronounced in NT than in PT.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Simulations for Planning of Liquid Hydrogen Spill Test

In order to better understand the complex pooling and vaporization of a liquid hydrogen spill, Sandia National Laboratories is conducting a highly instrumented, controlled experiment inside their Shock Tube Facility. Simulations were run before the experiment to help with the planning of experimental conditions, including sensor placement and cross wind velocity. This paper describes the modeling used in this planning process and its main conclusions. Sierra Suite’s Fuego, an in-house computational fluid dynamics code, was used to simulate a RANS model of a liquid hydrogen spill with five crosswind velocities: 0.45, 0.89, 1.34, 1.79, and 2.24 m/s. Two pool sizes were considered: a diameter of 0.85 m and a diameter of 1.7. A grid resolution study was completed on the smaller pool size with a 1.34 m/s crosswind. A comparison of the length and height of the plume of flammable hydrogen vaporizing from the pool shows that the plume becomes longer and remains closer to the ground with increasing wind speed. The plume reaches the top of the facility only in the 0.45 m/s case. From these results, we concluded that it will be best for the spacing and location of the concentration sensors to be reconfigured for each wind speed during the experiment.

08 HYDROGEN↗

Semikinetic modeling and observations of high-latitude plasma outflow

In one fo the most exciting areas of progress, we have now developed a dynamic semikinetic model for examining the synergistic effects of waves and magnetospheric hot plasma populations on the outflowing ionospheric plasma. We have done this by imposing hot biMaxwellian ion and electron distributions at the top of our auroral simulation flux tube (4 R(sub e)), as well as a spectrum of waves with altitude which perpendicularly heats the ionospheric ions. We have also addressed the quasi-statistical properties of out flowing O(+) through bulk parameter analysis of DE-1/RIMS observations when DE-1 was in the midaltitude polar cap magnetosphere. A very exciting paper which was both submitted and appeared during this period in JGR concerned the centrifugal acceleration effect on the polar wind. During this year, we also developed the extension of our GSK code to include the lower ionosphere and associated processes. Finally, perhaps our most exciting project currently is a study of F-region upflows using a modification of the FLIP ionospheric fluid dynamics model of Richards and Torr to allow for the effects of soft electron precipitation ionization and convection-driven ion heating, and performed comparisons with satellite and radar data. A list of reference papers that have been published is listed.

Horwitz, James L.↗

High-Temperature Gas-Cooled Pebble-Bed Reactors Running In And Transient Modeling Capabilities Demonstration

This study presents a comprehensive benchmarking and verification effort of several thermal-hydraulic and multiphysics capabilities for high-temperature gas-cooled reactor (HTGR) applications. The first part of this effort focuses on the running-in verification of Griffin's multiphysics capabilities, specifically for simulating the evolution of Pebble Bed reactor cores from startup to equilibrium. In the absence of validation data, code-to-code comparisons are conducted with Kugelpy, showing good agreement for key quantities like maximum power density and fresh core k-eigenvalue predictions. However, discrepancies in equilibrium core predictions suggest potential issues with cross sections, underscoring the need for further refinement and evaluation. The HTTF system analysis code benchmark involves RELAP5-3D, SAM, and GAMMA+ to assess their predictive capabilities for HTTF behavior under both normal operation and pressurized conduction cooldown (PCC) transient conditions. While there is good agreement in predicting major parameters such as coolant temperature, solid temperature, and flow distribution, discrepancies in transient behavior highlight differences in modeling approaches, nodalizations, and heat transfer models. The HTTF lower plenum CFD benchmark employs nekRS to simulate flow mixing phenomena, successfully capturing relevant flow physics and demonstrating mesh independence in complex geometries. Preliminary results suggest a relatively uniform temperature field but significant unsteadiness in the flow, requiring time-averaging analyses. The GPBR200 system analysis code benchmark uses SAM's core channel and porous media models, incorporating an RCCS loop for decay heat removal. During steady-state and transient conditions, including protected de-pressurized and pressurized loss of forced cooling (DLOFC and PLOFC), both models show good agreement in predicting temperature profiles and key parameters. Notably, while the core channel model underpredicts convective heat transfer effects, both models maintain temperatures well below the TRISO fuel safety limit. These benchmarking efforts collectively enhance the predictive capabilities of the tools used in HTGR design and safety analysis, guiding developments to improve their accuracy and applicability.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A Comparison of Rotor Disk Modeling and Blade-Resolved CFD Simulations for NASA's Tiltwing Air Taxi

A multi-fidelity computational fluid dynamics analysis is carried out for NASA’s tiltwing air taxi concept operating in airplane and helicopter mode. High-fidelity simulations are computationally expensive due to individual rotor blade modeling in a time-dependent computational domain with rotating grids. The mid-fidelity rotor disk option, in its source term implementation, is explored as a more affordable alternative. Computations are performed with NASA’s OVERFLOW flow solver loosely-coupled with the comprehensive code CAMRAD II for appropriate rotor trim. Detailed comparisons are shown for the trim solution, airloads, wake geometry, and rotor performance. While the rotor disk model is able to capture the flow field with satisfactory agreement in airplane mode, it faces difficulties in helicopter mode due to the three-dimensional effects of the wake. Although this study is limited to a specific vehicle geometry, it is expected that the results are somewhat generalizable to the analysis of multi-rotor configurations.

ARMD↗

Development of an Improved RELAP5-3D Model for the High Temperature Test Facility

High-temperature gas-cooled reactors (HTGRs) are rapidly approaching deployment. Confidence in transient analysis of these systems for design, optimization, and licensing calculations requires modeling and simulation tools that have been validated against data relevant to HTGR conditions. The High Temperature Test Facility (HTTF) is an integral effects thermal hydraulics test facility for prismatic HTGRs. In spring and summer of 2019, HTTF was used for a series of experiments that now serve as the basis for the OECD/NEA Thermal Hydraulic Code Validation Benchmark for High Temperature Gas-Cooled Reactors using HTTF Data (HTGR T/H Benchmark). This benchmark contains problems for systems code, computational fluid dynamics (CFD), and coupled systems code/CFD modeling representing lower plenum mixing and both the depressurized and pressurized conduction cooldown (DCC and PCC respectively) transients. Benchmark problems include exercises for code-to-code and code-to-data comparisons as well as an exercise for error scaling between HTTF and the Modular High Temperature Gas-Cooled Reactor, which serves as the basis for the HTTF design. Previous analysis as part of the HTGR T/H benchmark used a RELAP5-3D model developed at Idaho National Laboratory (INL) and demonstrated an ability to reproduce trends in the measured data but difficulties reproducing experimental values within their uncertainty. These difficulties were largely attributed to assumptions made during the development of the initial RELAP5-3D model, which predated the HTTF experiments. A significant cause of difficulty reproducing the measured temperatures may be the radial nodalization of the previous RELAP5-3D model. The new model provides a finer nodalization to assess the impact of radial nodalization and allows for asymmetric heating within the core, which was a feature of multiple HTTF experiments. In this paper, we present the new RELAP5-3D model of HTTF. In addition to describing the new model, this paper compares the new and old models and provides results for a full-power steady state, a DCC, and a PCC in HTTF. These analyses are based on the code-to-code comparison exercises for the DCC and PCC problems of the HTGR T/H benchmark. We present the results of these exercises from the new model and compare them to the results of the old model.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Aerodynamic analysis for aircraft with nacelles, pylons, and winglets at transonic speeds

A computational method has been developed to provide an analysis for complex realistic aircraft configurations at transonic speeds. Wing-fuselage configurations with various combinations of pods, pylons, nacelles, and winglets can be analyzed along with simpler shapes such as airfoils, isolated wings, and isolated bodies. The flexibility required for the treatment of such diverse geometries is obtained by using a multiple nested grid approach in the finite-difference relaxation scheme. Aircraft components (and their grid systems) can be added or removed as required. As a result, the computational method can be used in the same manner as a wind tunnel to study high-speed aerodynamic interference effects. The multiple grid approach also provides high boundary point density/cost ratio. High resolution pressure distributions can be obtained. Computed results are correlated with wind tunnel and flight data using four different transport configurations. Experimental/computational component interference effects are included for cases where data are available. The computer code used for these comparisons is described in the appendices.

Boppe, Charles W.↗

A detailed description of the uncertainty analysis for High Area Ratio Rocket Nozzle tests at the NASA Lewis Research Center

A preliminary uncertainty analysis has been performed for the High Area Ratio Rocket Nozzle test program which took place at the altitude test capsule of the Rocket Engine Test Facility at the NASA Lewis Research Center. Results from the study establish the uncertainty of measured and calculated parameters required for the calculation of rocket engine specific impulse. A generalized description of the uncertainty methodology used is provided. Specific equations and a detailed description of the analysis are presented. Verification of the uncertainty analysis model was performed by comparison with results from the experimental program's data reduction code. Final results include an uncertainty for specific impulse of 1.30 percent. The largest contributors to this uncertainty were calibration errors from the test capsule pressure and thrust measurement devices.

Kenneth J. Davidian↗

A detailed description of the uncertainty analysis for high area ratio rocket nozzle tests at the NASA Lewis Research Center

A preliminary uncertainty analysis was performed for the High Area Ratio Rocket Nozzle test program which took place at the altitude test capsule of the Rocket Engine Test Facility at the NASA Lewis Research Center. Results from the study establish the uncertainty of measured and calculated parameters required for the calculation of rocket engine specific impulse. A generalized description of the uncertainty methodology used is provided. Specific equations and a detailed description of the analysis is presented. Verification of the uncertainty analysis model was performed by comparison with results from the experimental program's data reduction code. Final results include an uncertainty for specific impulse of 1.30 percent. The largest contributors to this uncertainty were calibration errors from the test capsule pressure and thrust measurement devices.

Davidian, Kenneth J.↗

The Effect of Delamination on Damage Path and Failure Load Prediction for Notched Composite Laminates

The influence of delamination on the progressing damage path and initial failure load in composite laminates is investigated. Results are presented from a numerical and an experimental study of center-notched tensile-loaded coupons. The numerical study includes two approaches. The first approach considers only intralaminar (fiber breakage and matrix cracking) damage modes in calculating the progression of the damage path. In the second approach, the model is extended to consider the effect of interlaminar (delamination) damage modes in addition to the intralaminar damage modes. The intralaminar damage is modeled using progressive damage analysis (PDA) methodology implemented with the VUMAT subroutine in the ABAQUS finite element code. The interlaminar damage mode has been simulated using cohesive elements in ABAQUS. In the experimental study, 2-3 specimens each of two different stacking sequences of center-notched laminates are tensile loaded. The numerical results from the two different modeling approaches are compared with each other and the experimentally observed results for both laminate types. The comparisons reveal that the second modeling approach, where the delamination damage mode is included together with the intralaminar damage modes, better simulates the experimentally observed damage modes and damage paths, which were characterized by splitting failures perpendicular to the notch tips in one or more layers. Additionally, the inclusion of the delamination mode resulted in a better prediction of the loads at which the failure took place, which were higher than those predicted by the first modeling approach which did not include delaminations.

Satyanarayana, Arunkumar↗

Depth Dose Distribution Study within a Phantom Torso after Irradiation with a Simulated Solar Particle Event at NSRL

The adequate knowledge of the radiation environment and the doses incurred during a space mission is essential for estimating an astronaut's health risk. The space radiation environment is complex and variable, and exposures inside the spacecraft and the astronaut's body are compounded by the interactions of the primary particles with the atoms of the structural materials and with the body itself Astronauts' radiation exposures are measured by means of personal dosimetry, but there remains substantial uncertainty associated with the computational extrapolation of skin dose to organ dose, which can lead to over- or underestimation of the health risk. Comparisons of models to data showed that the astronaut's Effective dose (E) can be predicted to within about a +10% accuracy using space radiation transport models for galactic cosmic rays (GCR) and trapped radiation behind shielding. However for solar particle event (SPE) with steep energy spectra and for extra-vehicular activities on the surface of the moon where only tissue shielding is present, transport models predict that there are large differences in model assumptions in projecting organ doses. Therefore experimental verification of SPE induced organ doses may be crucial for the design of lunar missions. In the research experiment "Depth dose distribution study within a phantom torso" at the NASA Space Radiation Laboratory (NSRL) at BNL, Brookhaven, USA the large 1972 SPE spectrum was simulated using seven different proton energies from 50 up to 450 MeV. A phantom torso constructed of natural bones and realistic distributions of human tissue equivalent materials, which is comparable to the torso of the MATROSHKA phantom currently on the ISS, was equipped with a comprehensive set of thermoluminescence detectors and human cells. The detectors are applied to assess the depth dose distribution and radiation transport codes (e.g. GEANT4) are used to assess the radiation field and interactions of the radiation field with the phantom torso. Lymphocyte cells are strategically embedded at selected locations at the skin and internal organs and are processed after irradiation to assess the effects of shielding on the yield of chromosome damage. The initial focus of the present experiment is to correlate biological results with physical dosimetry measurements in the phantom torso. Further on, the results of the passive dosimetry within the anthropomorphic phantoms represent the best tool to generate reliable data to benchmark computational radiation transport models in a radiation field of interest. The presentation will give first results of the physical dose distribution, the comparison with GEANT4 computer simulations based on a Voxel model of the phantom, and a comparison with the data from the chromosome aberration study.

Berger, Thomas↗

Evaluation of simulated HPGe detector efficiencies in OpenMC compared to MCNP

In late 2023, the open-source radiation transport code OpenMC introduced a pulse-height tally (PHT) feature, enabling users to track the total energy deposited by individual photons in OpenMC cells. This function represents an important improvement of the OpenMC code because it allows users to simulate the response function of a gamma detector without relying on closed-source alternatives. Despite this, limited work has been published evaluating OpenMC as a radiation transport code for simulating gamma spectroscopy experiments. This study attempts to demonstrate the usefulness of OpenMC in this space by directly comparing its PHT output to MCNP (a trusted industry-standard Monte Carlo code). In the first half of this study, a radiation transport experiment is described in which a detector with a complex internal geometry is exposed to various gamma-emitting isotopes ( 133 Ba, 137 Cs, and 60 Co) over a range of distances. The setup of this experiment was modeled in MCNP with sufficient detail to capture the efficiency characteristics of a high purity germanium (HPGe) detector for the primary gammas of each isotope. After good agreement between the MCNP model and experiment was achieved, an identical model was produced in OpenMC to allow direct comparison between the full energy peak (FEP) values produced in OpenMC and MCNP. The results show strong agreement between OpenMC and MCNP across the full range of tested energies, with each model’s FEP values typically within 2% of each other and most FEP areas within 4% of experimental data. Worse agreement was seen between the Monte Carlo codes and experiment below 300 keV (an expected result). For the 662 keV line of 137 Cs, both codes were found in poor agreement with experiment and each other over the full range of distances tested (possibly indicating an error with the 137 Cs experimental data). Ignoring the anomalous results of the 662 keV line, all other data show good qualitative and quantitative agreement between MCNP and OpenMC. This result demonstrates the accuracy of OpenMC’s PHT feature for spectroscopic applications in which detector efficiency is a primary concern.

07 - ISOTOPES AND RADIATION SOURCES↗

Comparison of laser-produced plasma spatio-temporal electron density evolution measured using interferometry with simulation results

Due to the difficulties associated with experimental measurements of laser-produced plasma (LPP) properties during the earliest stages of plasma evolution, radiation hydrodynamic codes are often used. However, although these codes have been extensively validated in the higher intensity regimes, validation at low to moderate intensities has been limited. In this study, the spatio-temporal electron density evolution in an LPP generated at moderate laser intensities and at various laser wavelengths was validated against the FLASH code for times up to 20 ns. The LPP was generated by focusing the fundamental and various harmonics radiation (1064, 532, and 266 nm) from a 6 ns full width half maximum Nd:YAG laser, at a laser intensity of 10 GW cm -2 , onto a copper target. The spatio-temporal density evolution of the expanding plasma was analyzed using Nomarski interferometry. Experimental measurements were found to be consistent with FLASH simulations, and the dependence of electron densities on wavelength was found to be in agreement with analytical models, varying as n e ∝ λ -0.7 . However, slight differences were noted in the widths and shape functions of the experimental and simulated electron density profiles.

2D fast imaging↗

An experimental comparison of nonswirling and swirling flow in a circular-to-rectangular transition duct

Circular-to-rectangular transition ducts are used as exhaust system components of aircraft with rectangular exhaust nozzles. Often, the incoming flow of these transition ducts includes a swirling velocity component remaining from the gas turbine engine. Previous transition duct studies have either not included inlet swirl or when inlet swirl was considered, only overall performance parameters were evaluated. Circular-to-rectangular transition duct flows with and without inlet swirl were explored in order to understand the effect of inlet swirl on the transition duct flow field and to provide detailed duct flow data for comparison with numerical code predictions. A method was devised to create a swirling, solid body rotational flow with minimal associated disturbances. Coefficients based on velocities and total and static pressures measured incross stream planes at four axial locations within the transition duct, along with surface static pressure measurements and surface oil film visualization, are presented for both nonswirling and swirling incoming flow. In both cases the inlet centerline Mach number was 0.35. The Reynolds number based on the inlet centerline velocity and duct inlet diameter was 1,547,000 for nonswirling and 1,366,000 for swirling flow. The maximum swirl angle was 15.6 deg. Two pair of counter-rotating side wall vortices appeared in the duct flow without inlet swirl. These vortices were absent in the swirling incoming flow cases.

Reichert, B. A.↗

Unstructured CFD and Noise Prediction Methods for Propulsion Airframe Aeroacoustics

Using unstructured mesh CFD methods for Propulsion Airframe Aeroacoustics (PAA) analysis has the distinct advantage of precise and fast computational mesh generation for complex propulsion and airframe integration arrangements that include engine inlet, exhaust nozzles, pylon, wing, flaps, and flap deployment mechanical parts. However, accurate solution values of shear layer velocity, temperature and turbulence are extremely important for evaluating the usually small noise differentials of potential applications to commercial transport aircraft propulsion integration. This paper describes a set of calibration computations for an isolated separate flow bypass ratio five engine nozzle model and the same nozzle system with a pylon. These configurations have measured data along with prior CFD solutions and noise predictions using a proven structured mesh method, which can be used for comparison to the unstructured mesh solutions obtained in this investigation. This numerical investigation utilized the TetrUSS system that includes a Navier-Stokes solver, the associated unstructured mesh generation tools, post-processing utilities, plus some recently added enhancements to the system. New features necessary for this study include the addition of two equation turbulence models to the USM3D code, an h-refinement utility to enhance mesh density in the shear mixing region, and a flow adaptive mesh redistribution method. In addition, a computational procedure was developed to optimize both solution accuracy and mesh economy. Noise predictions were completed using an unstructured mesh version of the JeT3D code.

Pao, S. Paul↗

Preliminary Structural Sizing and Alternative Material Trade Study of CEV Crew Module

This paper presents the results of a preliminary structural sizing and alternate material trade study for NASA s Crew Exploration Vehicle (CEV) Crew Module (CM). This critical CEV component will house the astronauts during ascent, docking with the International Space Station, reentry, and landing. The alternate material design study considers three materials beyond the standard metallic (aluminum alloy) design that resulted from an earlier NASA Smart Buyer Team analysis. These materials are graphite/epoxy composite laminates, discontinuously reinforced SiC/Al (DRA) composites, and a novel integrated panel material/concept known as WebCore. Using the HyperSizer (Collier Research and Development Corporation) structural sizing software and NASTRAN finite element analysis code, a comparison is made among these materials for the three composite CM concepts considered by the 2006 NASA Engineering and Safety Center Composite Crew Module project.

Bednarcyk, Brett A.↗