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 253 records · Page 14

Effect of Impact Location on the Response of Shuttle Wing Leading Edge Panel 9

The objective of this paper is to compare the results of several simulations performed to determine the worst-case location for a foam impact on the Space Shuttle wing leading edge. The simulations were performed using the commercial non-linear transient dynamic finite element code, LS-DYNA. These simulations represent the first in a series of parametric studies performed to support the selection of the worst-case impact scenario. Panel 9 was selected for this study to enable comparisons with previous simulations performed during the Columbia Accident Investigation. The projectile for this study is a 5.5-in cube of typical external tank foam weighing 0.23 lb. Seven locations spanning the panel surface were impacted with the foam cube. For each of these cases, the foam was traveling at 1000 ft/s directly aft, along the orbiter X-axis. Results compared from the parametric studies included strains, contact forces, and material energies for various simulations. The results show that the worst case impact location was on the top surface, near the apex.

Lyle, Karen H.↗

Simulation of Guided Wave Interaction with In-Plane Fiber Waviness

Reducing the timeline for certification of composite materials and enabling the expanded use of advanced composite materials for aerospace applications are two primary goals of NASA's Advanced Composites Project (ACP). A key a technical challenge area for accomplishing these goals is the development of rapid composite inspection methods with improved defect characterization capabilities. Ongoing work at NASA Langley is focused on expanding ultrasonic simulation capabilities for composite materials. Simulation tools can be used to guide the development of optimal inspection methods. Custom code based on elastodynamic finite integration technique is currently being developed and implemented to study ultrasonic wave interaction with manufacturing defects, such as in-plane fiber waviness (marcelling). This paper describes details of validation comparisons performed to enable simulation of guided wave propagation in composites containing fiber waviness. Simulation results for guided wave interaction with in-plane fiber waviness are also discussed. The results show that the wavefield is affected by the presence of waviness on both the surface containing fiber waviness, as well as the opposite surface to the location of waviness.

Leckey, Cara A. C.↗

Bison Verification and Validation Activities for TRISO

Numerical modeling and simulation (M&S) tools play a key role in the research, development, and overall safety assessments of next-generation nuclear energy systems. One such tool, Bison, is a nuclear fuel performance code that is applicable to many fuel forms (e.g., light-water reactor fuel, oxide and metallic fuel for fast reactors, tri-structural isotropic (TRISO) fuel, and plate fuel), and it uses the finite element method to model the thermo- mechanical response of nuclear fuels. One fuel form widely utilized in Generation-IV high-temperature gas-cooled and fluoride- salt-cooled nuclear reactor concepts is TRISO fuel. Recently, Bison’s capabilities were significantly expanded to enable it to model the performance of TRISO particles and compacts. It is important that Bison’s computational results be reliable and predictive, since this code is used to inform high-consequence decisions. The various processes developed to address this issue generally entail two fundamental steps: verification and validation (V&V). Verification ensures that the code functions correctly and is reliable. Code/solution verification, code benchmark, and software quality assurance exercises are examples of verification activities. On the other hand, validation is the process of assessing a code’s capability to accurately model physical problems. Comparisons between code results and experiments quantify the validation level. Application of V&V procedures is crucial to the development of computational tools that are free of coding mistakes and can accurately represent reality. The current study presents an overview of Bison V&V activities relevant to the TRISO fuel concept, which include code/solution verification exercises, CRP-6 Benchmark—a Coordinated Research Program through the International Atomic Energy Agency (IAEA)—exercises, and validation exercises with the Advanced Gas Reactor (AGR)- 1/2/3/4 experiment series.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Coupled 2-dimensional cascade theory for noise and unsteady aerodynamics of blade row interaction in turbofans. Volume 1: Theory development and parametric studies

Typical analytical models for interaction between rotor and stator in a turbofan analyze the effect of wakes from the rotor impinging on the stator, producing unsteady loading, and thereby generating noise. Reflection/transmission characteristics of the rotor are sometimes added in a separate calculation. In those models, there is a one-to-one relationship between wake harmonics and noise harmonics; that is, the BPF (blade passing frequency) wake harmonic causes only the BPF noise harmonic, etc. This report presents a more complete model in which flow tangency boundary conditions are satisfied on two cascades in relative motion for several harmonics simultaneously. By an extension of S.N. Smith's code for two dimensional flat plate cascades, the noise generation/frequency scattering/blade row reflection problem is solved in a single matrix inversion. It is found that the BPF harmonic excitation of the stator scatters considerable energy in the higher BPF harmonics due to relative motion between the blade rows. Furthermore, when swirl between the rotor and stator is modeled, a 'mode trapping' effect occurs which explains observations on fans operating at rotational speeds below BFP cuton: the BPF mode amplifies between blade rows by multiple reflections but cannot escape to the inlet and exit ducts. However, energy scattered into higher harmonics does propagate and dominates the spectrum at two and three times BPF. This report presents the complete derivation of the theory, comparison with a previous (more limited) coupled rotor/stator interaction theory due to Kaji and Okazaki, exploration of the mode trapping phenomenon, and parametric studies showing the effects of vane/blade ratio and rotor/stator interaction. For generality, the analysis applies to stages where the rotor is either upstream or downstream of the stator and to counter rotation stages. The theory has been coded in a FORTRAN program called CUP2D, documented in Volume 2 of this report. It is concluded that the new features of this analysis - unsteady coupling, frequency scattering, and flow turning between rotor and stator - have a profound effect on noise generation caused by rotor/stator interaction. Treating rotors and stators as isolated cascades is not adequate for noise analysis and prediction.

Hanson, Donald B.↗

Quantifying deformation during Zry-4 burst testing: a comparison of BISON and a combined in-situ digital image correlation and infrared thermography method

The development and verification of nuclear fuel-cladding performance codes require extensive testing to establish empirical correlations, necessitating more accelerated testing methods that can replace large validation studies. Here, in the present work, a technique was developed and implemented to experimentally quantify the relationship between temperature, cladding internal pressure, and strain in-situ during a burst test, with the specific aim of generating data for code validation and model refinement. Digital image correlation was used to measure hoop, axial, and radial strain, and infrared thermography to quantify axial temperature gradients. Several key experimental modifications to traditional LOCA burst testing were necessary, but are shown to effectively have little impact on burst temperatures. The time/temperature dependent pressure and time/spatially dependent thermal gradient data were used with BISON to simulate cladding burst and strain rates. Although the measured DIC strains reach appreciable amounts at slightly lower temperatures than BISON simulated strains, there is good agreement between the measured and simulated strains in terms of magnitude and rates leading up to burst. However, in the few seconds prior to and during burst, the BISON simulated strains are significantly lower than the measured strains, indicating the potential for model improvement. The combined experiment and simulation technique may be applied to any new cladding concept to accelerate fuel qualification.

36 MATERIALS SCIENCE↗

Genome-wide discovery of G-quadruplexes in barley

G-quadruplexes (G4s) are four-stranded nucleic acid structures with closely spaced guanine bases forming square planar G-quartets. Aberrant formation of G4 structures has been associated with genomic instability. However, most plant species are lacking comprehensive studies of G4 motifs. In this study, genome-wide identification of G4 motifs in barley was performed, followed by a comparison of genomic distribution and molecular functions to other monocot species, such as wheat, maize, and rice. Similar to the reports on human and some plants like wheat, G4 motifs peaked around the 5' untranslated region (5' UTR), the first coding domain sequence, and the first intron start sites on antisense strands. Our comparative analyses in human, Arabidopsis, maize, rice, and sorghum demonstrated that the peak points could be erroneously merged into a single peak when large window sizes are used. We also showed that the G4 distributions around genic regions are relatively similar in the species studied, except in the case of Arabidopsis. G4 containing genes in monocots showed conserved molecular functions for transcription initiation and hydrolase activity. Additionally, we provided examples of imperfect G4 motifs.

54 ENVIRONMENTAL SCIENCES↗

HTTF Benchmarking Activities in FY 2023

This report describes the work conducted in FY23 within the High Temperature Gas-Cooled Reactor (HTGR) Applications Drivers activity under the Multiphysics Applications technical area in the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program within DOE-NE. The focus of this past year’s activity was on the OECD/NEA thermal hydraulic benchmark for high temperature gas-cooled reactors using HTTF data (HTGR T/H). The activities focus on code-to-code comparisons for two benchmark problems: (1) Depressurized Conduction Cooldown using SAM, and (2) Lower Plenum Mixing using Nek5000/RS.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Evaluation of CFD Turbulent Heating Prediction Techniques and Comparison With Hypersonic Experimental Data

Results from a study to assess the accuracy of turbulent heating and skin friction prediction techniques for hypersonic applications are presented. The study uses the original and a modified Baldwin-Lomax turbulence model with a space marching code. Grid converged turbulent predictions using the wall damping formulation (original model) and local damping formulation (modified model) are compared with experimental data for several flat plates. The wall damping and local damping results are similar for hot wall conditions, but differ significantly for cold walls, i.e., T(sub w) / T(sub t) < 0.3, with the wall damping heating and skin friction 10-30% above the local damping results. Furthermore, the local damping predictions have reasonable or good agreement with the experimental heating data for all cases. The impact of the two formulations on the van Driest damping function and the turbulent eddy viscosity distribution for a cold wall case indicate the importance of including temperature gradient effects. Grid requirements for accurate turbulent heating predictions are also studied. These results indicate that a cell Reynolds number of 1 is required for grid converged heating predictions, but coarser grids with a y(sup +) less than 2 are adequate for design of hypersonic vehicles. Based on the results of this study, it is recommended that the local damping formulation be used with the Baldwin-Lomax and Cebeci-Smith turbulence models in design and analysis of Hyper-X and future hypersonic vehicles.

Dilley, Arthur D.↗

The AGORA High-resolution Galaxy Simulations Comparison Project. V. Satellite Galaxy Populations in a Cosmological Zoom-in Simulation of a Milky Way–Mass Halo

We analyze and compare the satellite halo populations at z ~ 2 in the high-resolution cosmological zoom-in simulations of a 10 12 M ⊙ target halo (z = 0 mass) carried out on eight widely used astrophysical simulation codes (Art-I, Enzo, Ramses, Changa, Gadget-3, Gear, Arepo-t, and Gizmo) for the AGORA High-resolution Galaxy Simulations Comparison Project. We use slightly different redshift epochs near z = 2 for each code (hereafter "z ~ 2") at which the eight simulations are in the same stage in the target halo's merger history. After identifying the matched pairs of halos between the CosmoRun simulations and the DMO simulations, we discover that each CosmoRun halo tends to be less massive than its DMO counterpart. When we consider only the halos containing stellar particles at z ~ 2, the number of satellite galaxies is significantly fewer than that of dark matter halos in all participating AGORA simulations and is comparable to the number of present-day satellites near the Milky Way or M31. The so-called "missing satellite problem" is fully resolved across all participating codes simply by implementing the common baryonic physics adopted in AGORA and the stellar feedback prescription commonly used in each code, with sufficient numerical resolution (≲100 proper pc at z = 2). We also compare other properties such as the stellar mass–halo mass relation and the mass–metallicity relation. Our work highlights the value of comparison studies such as AGORA, where outstanding problems in galaxy formation theory are studied simultaneously on multiple numerical platforms.

79 ASTRONOMY AND ASTROPHYSICS↗

Comparison Analysis Between Nevada and TSS Using a Jet Engine Exhaust Nozzle Test Rig

Thermal analysis in both simple and complex models can require calculating the propagation of radiant energy to and from multiple surfaces. This can be accomplished through simple estimation techniques or complex computationally intense computer modeling simulations. Currently there are a variety of computer analysis techniques used to simulate the propagation of radiant energy, each having advantages and disadvantages. The major objective of this effort was to compare two ray tracing radiation propagation analysis programs (NEVADA and TSS) Net Energy Verification and Determination Analyzer and Thermal Synthesizer System with experimental data. Results from a non-flowing, electrically heated test rig was used to verify the calculated radiant energy propagation from a nozzle geometry that represents an aircraft propulsion nozzle system. In general the programs produced comparable overall results, and results slightly higher then the experimental data. Upon inspection of individual radiation interchange factors, differences were evident and would have been magnified if a more radical model temperature profile was analyzed. Bidirectional reflectivity data (BRDF) was not used do to modeling limitations in TSS. For code comparison purposes, this nozzle geometry represents only one case for one set of analysis conditions. Since each computer code has advantages and disadvantages bases on scope, requirements, and desired accuracy, the usefulness of this single case study may be limiting.

Baumeister, Joseph F.↗

HST Studies of the Chromospheres, Wind, and Mass-Loss Rates of Cool Giant and Supergiant Stars

UV spectra of K-M giant and supergiant stars and of carbon stars have been acquired with the Goddard High Resolution Spectrograph (GHRS) on the Hubble Space Telescope (HST). These spectra have been used to measure chromospheric flow and turbulent velocities, study the acceleration of their stellar winds, acquire constraints on their outer atmospheric structure, and enable estimates of their mass-loss rates. Results from our observations of the giant stars Gamma Dra (K5 III hybrid), Alpha Tau (K5 III), Gamma Cru (M3.4 III), Mu Gem (M3 IIIab), and 30 Her (MG III), the supergiants Alpha Ori (M2 Iab) and Lambda Vel (K5 Ib), and the carbon stars TX Psc (NO; C6,2) and TW Hor (NO; C7,2) will be summarized and compared. The high resolution and wavelength accuracy of these data have allowed the direct measurement of the acceleration of the stellar winds in the chromospheres of several of these stars (from initial velocities of 3-9 km/s to upper velocities of 15-25 km/s) and of the chromospheric macroturbulence (-25-35 km/s). The high signal-to-noise and large dynamic range of these spectra have allowed the detection and identification of numerous new emission features, including weak C IV emission indicative of hot transition-region plasma in the non-coronal giant Alpha Tau, many new fluorescent lines of Fe II, and the first detection of fluorescent molecular hydrogen emission and of Ca II recombination lines in the UV spectrum of a giant star. The UV spectrum of two carbon stars have been studied with unprecedented resolution and reveal extraordinarily complicated Mg II lines nearly smothered by circumstellar absorptions. Finally, comparison of synthetic UV emission line profiles computed with the Lamers et al. (1987) Sobolev with Exact Integration (SEI) code with observations of chromospheric emission lines overlain with wind absorption features provides estimates of the mass-loss rates for four of these stars.

Carpenter, Kenneth G.↗

Timing Analysis with INTEGRAL: Comparing Different Reconstruction Algorithms

INTEGRAL is one of the few instruments capable of detecting X-rays above 20keV. It is therefore in principle well suited for studying X-ray variability in this regime. Because INTEGRAL uses coded mask instruments for imaging, the reconstruction of light curves of X-ray sources is highly non-trivial. We present results from the comparison of two commonly employed algorithms, which primarily measure flux from mask deconvolution (ii-lc-extract) and from calculating the pixel illuminated fraction (ii-light). Both methods agree well for timescales above about 10 s, the highest time resolution for which image reconstruction is possible. For higher time resolution, ii-light produces meaningful results, although the overall variance of the lightcurves is not preserved.

Grinberg, V.↗

Parametric study of potential early commercial power plants Task 3-A MHD cost analysis

The development of costs for an MHD Power Plant and the comparison of these costs to a conventional coal fired power plant are reported. The program is divided into three activities: (1) code of accounts review; (2) MHD pulverized coal power plant cost comparison; (3) operating and maintenance cost estimates. The scope of each NASA code of account item was defined to assure that the recently completed Task 3 capital cost estimates are consistent with the code of account scope. Improvement confidence in MHD plant capital cost estimates by identifying comparability with conventional pulverized coal fired (PCF) power plant systems is undertaken. The basis for estimating the MHD plant operating and maintenance costs of electricity is verified.

Source record↗

Numerical Investigation of Radiative Heat Transfer in Laser Induced Air Plasmas

Radiative heat transfer is one of the most important phenomena in the laser induced plasmas. This study is intended to develop accurate and efficient methods for predicting laser radiation absorption and plasma radiative heat transfer, and investigate the plasma radiation effects in laser propelled vehicles. To model laser radiation absorption, a ray tracing method along with the Beer's law is adopted. To solve the radiative transfer equation in the air plasmas, the discrete transfer method (DTM) is selected and explained. The air plasma radiative properties are predicted by the LORAN code. To validate the present nonequilibrium radiation model, several benchmark problems are examined and the present results are found to match the available solutions. To investigate the effects of plasma radiation in laser propelled vehicles, the present radiation code is coupled into a plasma aerodynamics code and a selected problem is considered. Comparisons of results at different cases show that plasma radiation plays a role of cooling plasma and it lowers the plasma temperature by about 10%. This change in temperature also results in a reduction of the coupling coefficient by about 10-20%. The present study indicates that plasma radiation modeling is very important for accurate modeling of aerodynamics in a laser propelled vehicle.

Liu, J.↗

Serpent and MCNP Calculations of the Energy Deposition in the Transformational Challenge Reactor

This paper focuses on the calculation of the energy deposition in the Transformational Challenge Reactor by two major Monte Carlo codes: Serpent and MCNP. The first software computation relies on Kinetic Energy Released per unit Mass (KERMA) factors while the second one relies on Q-values. The results from these two independent computation methodologies are in very good agreement; however, Serpent runs much faster than MCNP (for the same computational model) and allows for a detailed energy deposition distribution from a 1-mm-side square mesh with a relative statistical error between 0.5% and 1%. This detailed energy deposition is suitable for multiphysics analyses aimed at design optimizations. In order to calculate the energy deposition, Serpent needs enhanced ACE files (distributed by the software developers). Unlike other Monte Carlo software that uses inputs based on Python or Java languages, the Serpent input syntax is very similar to that of MCNP; a Python script can convert a MCNP input to a Serpent input in seconds. For simulations not requiring the calculation of the energy deposition, Serpent can also read nuclear data from MCNP ACE files, which eventually improves the comparison of the results of the two codes.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A Survey of Multimaterial Treatments for Thermal Radiative Transfer

Arbitrary Lagrangian-Eulerian methods are a popular choice for hydrodynamic modeling in radiation (rad-hydro) simulations. Because these methods involve a relaxation step that moves the mesh relative to material boundaries, multimaterial spatial zones are generally present. Accurate treatments of these zones are needed to resolve various physical phenomena of interest for inertial confinement fusion applications. However, these codes are often paired with single-material, deterministic thermal radiative transfer (TRT) codes that are oblivious to the material compositions of each zone. These single-material TRT codes can only accept homogenized material properties (opacities, specific heats, etc.) from the hydrodynamic code and output homogenized solutions. After each TRT time step, the multimaterial hydrodynamic code must dehomogenize the quantities computed by the TRT package in order to update subzonal material temperatures. The process by which hydrodynamic codes perform this dehomogenization has not been well documented in previous literature, and the methods can vary significantly from code to code. The purpose of this paper is to document, study, and compare existing techniques used for rad-hydro simulations as well as present a new method with potentially promising results. We summarize several methods and give comparisons on infinite-medium problems as well a finite-medium problem for two of the methods.

42 ENGINEERING↗

Involute Working Group – Validation of CFD Turbulence Models for Steady-State Safety Analysis

There are three research reactors in the world having fuel plates curved as circle-involute (a spiral generated around a circle): The Oak Ridge National Laboratory (ORNL) High Flux Isotope Reactor (HFIR) located in Tennessee, U.S.A. [ORNL, 2018]; The Laue-Langevin Institute (ILL) High Flux Reactor (RHF) located in Grenoble, France [ILL, 2018]; and, The Technical University of Munich (TUM) Research Neutron Source Heinz Maier-Leibnitz (FRM II) located in Garching, Germany [TUM, 2018]. All three reactors are currently using Highly Enriched Uranium (HEU, 235U/U ≥ 20 wt. %) as fuel and all three are actively engaged in activities to convert to Low Enriched Uranium (LEU, 235U/U < 20 wt. %) fuel. More information on these reactors can be found in Section two. For various reasons, these reactors have expressed interest in using Computational Fluid Dynamics (CFD) tools to perform their Steady-State Thermal-Hydraulic (SSTH) safety calculations. Using CFD tools for SSTH analysis represents generally a significant departure from traditional methods and requires thorough verification and validation (V&V) to be accepted by regulators. This is why Argonne National Laboratory (ANL) and the involute reactors formed an informal group to help each other in this endeavor. This so-called Involute Working Group (IWG) aims at supporting the qualification of CFD tools for SSTH safety calculations. Activities include benchmarking, code-to-code comparison, Verification and Validation (V&V) as well as technical support. More information on the IWG and the reasons why involute-plate reactors are interested in CFD can be found in Section three. The present study aims at testing three turbulence models commonly implemented in commercial software (namely K-Epsilon, K-Omega and Reynolds-Stress) and assessing their relevance for involute-plate reactor application. An experiment performed by Gambill and Bundy in 1961 was performed to measure heat transfer coefficients in conditions prototypic of the HFIR reactor [Gambill, 1961]. More details on this experiment are provided in Section four. Based on the information available, CFD models of this experiment have been created and the experimental results have been compared to the ones obtained with these CFD models. CFD models and the results of this comparison are provided in Section five.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Examination of the collision force method for analyzing the responses of simple containment/deflection structures to impact by one engine rotor blade fragment

An approximate collision analysis, termed the collision-force method, was developed for studying impact-interaction of an engine rotor blade fragment with an initially circular containment ring. This collision analysis utilizes basic mass, material property, geometry, and pre-impact velocity information for the fragment, together with any one of three postulated patterns of blade deformation behavior: (1) the elastic straight blade model, (2) the elastic-plastic straight shortening blade model, and (3) the elastic-plastic curling blade model. The collision-induced forces are used to predict the resulting motions of both the blade fragment and the containment ring. Containment ring transient responses are predicted by a finite element computer code which accommodates the large deformation, elastic-plastic planar deformation behavior of simple structures such as beams and/or rings. The effects of varying the values of certain parameters in each blade-behavior model were studied. Comparisons of predictions with experimental data indicate that of the three postulated blade-behavior models, the elastic-plastic curling blade model appears to be the most plausible and satisfactory for predicting the impact-induced motions of a ductile engine rotor blade and a containment ring against which the blade impacts.

Zirin, R. M.↗