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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Characterization of sub-picosecond laser-produced fast electrons by modeling angularly resolved bremsstrahlung measurements with 3D hybrid particle-in-cell code

Hard x-rays produced by intense laser-produced fast electrons interacting with solids are a vital source for producing radiographs of high-density objects and implosion cores for inertial confinement fusion. Accurate calculation of hard x-ray sources requires a three-dimensional (3D) simulation geometry that fully models the electron transport dynamics, including electron recirculation and the generation of absolute photon yields. To date, 3D simulations of laser-produced bremsstrahlung photons over tens of picoseconds and code benchmarking have not been performed definitively. In this study, we characterize sub-picosecond laser-produced fast electrons by modeling angularly resolved bremsstrahlung measurements for refluxing and non-refluxing targets using the 3D hybrid particle-in-cell (PIC), Large Scale Plasma code. Bremsstrahlung radiation and escaped electron data were obtained by focusing a 50-TW Leopard laser (15 J, 0.35 ps, 2 × 10 19 W/cm 2 ) on a 100-μm-thick Cu foil and a Cu with a large plastic backing (Cu–CH target). Data for both the Cu and Cu–CH targets were reproduced for simulations with a given set of electron parameters. Comparison of the simulations revealed that the hard x-ray emission from the Cu target was significantly longer in duration than that from the Cu–CH target. Here, the benchmarked hybrid PIC code could prove to be a powerful tool in the design and optimization of time- and angular-dependent bremsstrahlung sources for flash x-ray and gamma-ray radiography.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2001

Hypersonic Flows About a 25 degree Sharp Cone

This paper presents the results of a numerical study that examines the surface heating discrepancies observed between computed and measured values along a sharp cone. With Mach numbers of an order of 10 and the freestream length Reynolds number of an order of 10 000, the present computations have been made with the direct simulation Monte Carlo (DSMC) method by using the G2 code of Bird. The flow conditions are those specified for two experiments conducted in the Veridian 48-inch Hypersonic Shock Tunnel. Axisymmetric simulations are made since the test model was assumed to be at zero incidence. Details of the current calculations are presented, along with comparisons between the experimental data, for surface heating and pressure distributions. Results of the comparisons show major differences in measured and calculated results for heating distributions, with differences in excess of 25 percent for the two cases examined.

Moss, James N.↗
DOE OSTI2025

Comparison of interlaminar damage modeling strategies for hybrid composite/aluminum laminates subjected to low-velocity impact

Low-velocity impact of hybrid metal-composite structures was investigated experimentally and computationally. Composite laminates consisting of 2D woven glass fiber reinforced polymer (GFRP) and carbon fiber reinforced polymer (CFRP) were joined with a 6061-T6 aluminum plate using an epoxy adhesive. Two variations of the structure were studied; one consisting of all plies oriented at 0° and one consisting of all plies oriented at 45°. A drop tower was used to impact structures at a range of energies, including energies above and below the threshold at which the aluminum layer was perforated. Numerical simulations were implemented using Sierra/SM, an in-house transient dynamics finite element code developed at Sandia National Laboratories. A Hosford plasticity model was used to describe the response of the aluminum layer. A newly implemented orthotropic continuum damage mechanics (CDM) constitutive model was used to represent the composite laminate. This 3D-CDM model was compared to a cohesive zone model (2D-CDM/CZM) to investigate efficacy of aluminum perforation energy prediction, delamination prediction, and computational cost. Accuracy of each model was evaluated using the experimental results. Each showed good agreement with the tests for both the force and velocity histories, as well as the observed damage mechanisms. The 2D-CDM/CZM model was marginally more accurate in capturing both the composite and aluminum behavior — this model averaged error percentages of -11.2% and 10.8% for residual velocity and peak force, respectively. Meanwhile, the 3D-CDM model predictions yielded average error percentages of -35.5% (velocity) and 22.6% (force). However, the 3D-CDM model generally resulted in a decreased computational cost; the average run time was 14% shorter than the 2D-CDM/CZM model and 3x as many timesteps per hour were computed using the same computational resources. In conclusion, new experimental data on the impact and perforation resistance of metal-composite laminates is presented in addition to numerical predictions of the impact behavior.

Carbon fiber↗
DOE OSTI2021

GOTHIC Aerosol Source Depletion Studies

Pacific Northwest National Laboratory has continued work to develop an aerosol-laden flow modeling capability with the Generation of Thermal Hydraulic Information in Containment (GOTHIC TM ) computer code to perform simulations for thermal hydraulic conditions and aerosol transport and deposition in spent fuel casks. This report describes our recent work to expand our model that was originally developed in 2019 to allow for thermal characterization, carrier gas flow determination, and tracking of particulate behavior throughout the entire canister volume. This was achieved through conversion of decay heat source in all fuel tubes within the model from heaters to thermal conductors in GOTHIC as well as remeshing the internal volume of the canister. These model improvements allowed for simulations of base case scenarios for comparison of similar efforts by our collaborators with alternative tools at Sandia and Oak Ridge. Several characteristics of the GOTHIC code were elucidated by this effort and were documented for consideration in how the code is used for this effort moving forward. For instance, treatment of particle size distributions, initial spatial distributions, mesh gradients, and minimum volume concentrations were identified as important issues for consideration. Ultimately, the current version of the code is capable of tracking temperatures, flow rates, and particle behavior throughout the canister internal volume. This report presents predicted depletion times and preferred deposition patterns within the canister. Additionally, this report documents code performance and insights for further model development and use for comparisons to other codes and for predicting system behavior in experiments being performed and planned by our collaborators.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1999

Evolution of Plasmaspheric Refilling: A Comparison of Measurements with an Interhemispheric Plasmasphere Model

The evolution of plasmasphere is closely coupled to processes governing flux tube refilling and transport. In this study we follow this evolution through a series of sequential ion observations made by RIMS on the DE 1 satellite and with the intervening time intervals filled by simulations made with the FLIP model. The FLIP model solves the continuity and momentum equations for the major ion species as well as the energy equations for ions and electrons along entire flux tubes rom 100 k altitude in both hemispheres. Convection has recently been included, driven a dawn-dusk electric field model. Since the observations are necessarily restricted to two local time passes through the plasmasphere per orbit, the FLIP code will follow the observed flux tubes as they corotate and convect through the other local times. The period of DE 1 is such that nearly the same flux tubes are observed every second day or every sixth orbit. The composite of observations and simulations allows assessment of the degree to which the degree to which the FLIP simulations accurately describe the processes which govern flux tube refilling over the two-day period between observations of the same flux tubes.

Comfort, Richard H.↗
DOE OSTI2025

Validation of 3D MHD simulations of Ne/D 2 mixed shattered pellet injection in JET

Nonlinear 3D magnetohydrodynamic (MHD) modeling of shattered pellet injection (SPI) in JET tokamak plasmas is performed with the JOREK code. The study focuses on the validation of simulation result with respect to experimental observations, addressing in particular figures of merit for the efficiency of the SPI, as a technique to mitigate thermal loads by injecting radiative impurities like neon during the thermal quench (TQ) phase of disruptions. A set of JET pulses with neon/deuterium atomic mixture ratio in the shattered pellet varying from 10% to 100% provides the experimental data to be compared with. Simulations using different models for the ablation of solid fragments and the radiation of neon impurities are considered. Synthetic diagnostics are employed for a direct quantitative comparison with key JET diagnostic systems providing in particular radiation, electron temperature and density, and magnetic measurements. 3D radiation structures from MHD simulations are analyzed to estimate the toroidal asymmetry of radiation, which is difficult to measure in JET due to the presence of only two toroidally displaced multi-channel bolometry systems. JOREK simulations show that, before the TQ, the bulk of the radiation is concentrated in a toroidal region enclosing the radiating fragments that is missed by both multi-channel bolometers, suggesting a possible underestimation of the total radiation by such diagnostic systems in the pre-TQ phase. On the other hand, JOREK predicts that during the TQ the radiation becomes more toroidally symmetric, with positive implications for the heat load to plasma facing components.

JOREK↗
NASA NTRS2000

Evolution of Plasmaspheric Refilling from Comparisons of Satellite Observations with Simulations by an Interhemispheric Plasmasphere Model

The evolution of the plasmasphere is closely coupled to processes governing flux tube refilling and transport. In this study we follow this evolution through a series of sequential ion observations made by the Retarding Ton Mass Spectrometer (RIMS) on the Dynamics Explorer 1 (DE1) satellite and with the intervening time intervals filled in by simulations made with the Field Line Interhemispheric Plasma (FLIP) model. The period of DE1 (approx. 7 hours) is such that nearly the same flux tubes are observed every second day, or every sixth orbit. The intervening observations provide additional information regarding longitudinal variations. Since the observations are necessarily restricted to two local time passes through the plasmasphere per orbit, the FLIP code will follow the observed flux tubes as they corotate and convect through the other local times. The simulations use the initial plasmasphere observations to define initial conditions of where the plasmapause is located and the extent of filling in the plasma trough. The composite of observations and simulations allows assessment of the degree to which the FLIP simulations accurately describe the processes which govern flux tube refilling over the two-day period between observations on the same flux tubes. To the extent that the convection electric field model in the FLIP code is accurate, this composite also provides a global picture of what is happening in and near the plasmasphere during the period. Of particular interest are changes in plasma composition as refilling proceeds.

Comfort, Richard H.↗
NASA NTRS1985

Studies on effects of boundary conditions in confined turbulent flow predictions

The differences in k epsilon model predictions of plane and axisymmetric expansion flows is investigated. The prediction of the coaxial jet for different velocity ratios of the annular to central jet is presented. The effects of inlet kinetic energy and the energy dissipation rate profiles are investigated for swirling and nonswirling flows. The effects of expansion ration and Reynolds number on the reattachment length are also presented. The results show that the inlet k and epsilon profiles have the most significant effect on the reattachment length and flow redevelopment for the case of coaxial jet of high velocity ratio. A comparison of k epsilon model predictions for the pipe expansion flow by the PHOENICS and TEACH codes reveals some discrepancies in the predicted results. TEACH prediction seems to produce unrealistic kinetic energy profiles in some regions of the flow. PHOENICS code produces a long tail in the recirculation region under certain conditions.

Nallasamy, M.↗
DOE OSTI2020

Benchmark probabilistic solar forecasts: Characteristics and recommendations

We illustrate and compare commonly used benchmark, or reference, methods for probabilistic solar forecasting that researchers use to measure the performance of their proposed techniques. A thorough review of the literature indicates wide variation in the benchmarks implemented in probabilistic solar forecast studies. To promote consistent and sensible methodological comparisons, we implement and compare ten variants from six common benchmark classes at two temporal scales: intra-hourly forecasts and hourly resolution forecasts. Using open-source Surface Radiation Budget Network (SURFRAD) data from 2018, these benchmark methods are compared using proper probabilistic metrics and common diagnostic tools. Practical implementation issues, such as the impact of missing data and applicability for operational forecasting, are also discussed. Furthermore, we make recommendations for practitioners on the appropriate selection of benchmark methods to properly showcase state-of-the-art improvements in forecast reliability and sharpness. All code and open-source data are available on Github for reproducibility and for other researchers to apply the same benchmark methods to their own data.

14 SOLAR ENERGY↗
NASA NTRS1993

Calibration of radiation codes in climate models: Comparison of calculations with observations from the SPECtral Radiation Experiment (SPECTRE)

The primary goal of SPECTRE is to: close the loopholes by which longwave radiation models have eluded incisive comparisons with measurements. Likewise, the experimental approach was quite simple in concept, namely: accurately measure the zenith infrared radiance at high spectral resolution while simultaneously profiling the radiatively important atmospheric properties with conventional and remote sensing devices. The field phase of SPECTRE was carried out as part of FIRE Cirrus II, and detailed spectra of the down welling radiance were obtained by several interferometers simultaneous to the measurement of the optical properties of the atmosphere. We are now well along in the process of analyzing the data and calibrating radiation codes so that they may be used more effectively in climate related studies. The calibration is being done with models ranging from the most detailed (line-by-line) to the broad-band parameterizations used in climate models. This paper summarizes our progress in the calibration for clear-sky conditions. When this stage is completed, we will move on to the calibration for cirrus conditions.

Ellingson, R. G.↗
NASA NTRS2002

Buckling Design Studies of Inverted, Oblate Bulkheads for a Propellant Tank

An investigation of the deformation and buckling characteristics of a composite, oblate bulkhead that has an inverted geometry and is subjected to pressure-only loading is presented for three bulkhead geometries and thicknesses. The effects of a stiffening support ring at the bulkhead to cylinder interface are also evaluated. Buckling analyses conducted using the axisymmetric shell code BOSOR4 are discussed for several bulkhead configurations. These results are analytically verified using results from the Structural Analysis of General Shells (STAGS) code for a selected bulkhead configuration. The buckling characterization of an inverted, oblate bulkhead requires careful attention as small changes in bulkhead parameters can have a significant effect on the critical buckling load. Comparison of BOSOR4 and STAGS results provided a very good correlation between the two analysis methods. In addition, the analysis code BOSOR4 was found to be an efficient sizing tool that is useful during the preliminary design stage of a practical shell structure. Together, these two aspects should give the design engineer confidence in sizing these stability critical structures. Additional characterization is warranted, especially for a composite tank structure, since only one bulkhead configuration was examined closely.

Smeltzer, Stanley S., III↗
NASA NTRS1999

Parametric Weight Comparison of Current and Proposed Thermal Protection System (TPS) Concepts

A parametric weight assessment of advanced metallic panel, ceramic blanket, and ceramic tile thermal protection systems (TPS) was conducted using an implicit, one-dimensional (1 -D) thermal finite element sizing code. This sizing code contained models to ac- count for coatings, fasteners, adhesives, and strain isolation pads. Atmospheric entry heating profiles for two vehicles, the Access to Space (ATS) rocket-powered single-stage-to-orbit (SSTO) vehicle and a proposed Reusable Launch Vehicle (RLV), were used to ensure that the trends were not unique to a particular trajectory. Eight TPS concepts were compared for a range of applied heat loads and substructural heat capacities to identify general trends. This study found the blanket TPS concepts have the lightest weights over the majority of their applicable ranges, and current technology ceramic tiles and metallic TPS concepts have similar weights. A proposed, state-of-the-art metallic system which uses a higher temperature alloy and efficient multilayer insulation was predicted to be significantly lighter than the ceramic tile systems and approaches blanket TPS weights for higher integrated heat loads.

Myers, David E.↗
NASA NTRS2004

A Computational Study of an Oscillating VR-12 Airfoil with a Gurney Flap

Computations of the flow over an oscillating airfoil with a Gurney-flap are performed using a Reynolds Averaged Navier-Stokes code and compared with recent experimental data. The experimental results have been generated for different sizes of the Gurney flaps. The computations are focused mainly on a configuration. The baseline airfoil without a Gurney flap is computed and compared with the experiments in both steady and unsteady cases for the purpose of initial testing of the code performance. The are carried out with different turbulence models. Effects of the grid refinement are also examined and unsteady cases, in addition to the assessment of solver effects. The results of the comparisons of steady lift and drag computations indicate that the code is reasonably accurate in the attached flow of the steady condition but largely overpredicts the lift and underpredicts the drag in the higher angle steady flow.

Rhee, Myung↗
NASA NTRS2014

Analysis of a Shock-Associated Noise Prediction Model Using Measured Jet Far-Field Noise Data

A code for predicting supersonic jet broadband shock-associated noise was assessed using a database containing noise measurements of a jet issuing from a convergent nozzle. The jet was operated at 24 conditions covering six fully expanded Mach numbers with four total temperature ratios. To enable comparisons of the predicted shock-associated noise component spectra with data, the measured total jet noise spectra were separated into mixing noise and shock-associated noise component spectra. Comparisons between predicted and measured shock-associated noise component spectra were used to identify deficiencies in the prediction model. Proposed revisions to the model, based on a study of the overall sound pressure levels for the shock-associated noise component of the measured data, a sensitivity analysis of the model parameters with emphasis on the definition of the convection velocity parameter, and a least-squares fit of the predicted to the measured shock-associated noise component spectra, resulted in a new definition for the source strength spectrum in the model. An error analysis showed that the average error in the predicted spectra was reduced by as much as 3.5 dB for the revised model relative to the average error for the original model.

Dahl, Milo D.↗
DOE OSTI2024

High- and Mid-Fidelity Modeling Comparison for a Floating Marine Turbine System

There is a lack of suitable numerical tools, particularly open-source tools, that can be used for designing and optimizing marine turbine systems. The National Renewable Energy Laboratory has added features to their widely used mid-fidelity wind turbine modeling code, OpenFAST, to enable modeling of axial-flow marine turbines. This necessitated the addition of several physical effects relevant to marine turbines that are neglected for wind turbines. These include buoyancy, added mass and inertial loads, wave-current superposition, and changes to the coordinate systems. This updated version of OpenFAST allows for the modeling of both fixed and floating marine turbines at a speed comparable to real time. While efficient for large sets of load cases and design studies, mid-fidelity codes make simplifying assumptions that may impact their accuracy. High-fidelity computational fluid dynamics (CFD) simulations can capture more flow effects with fewer assumptions and provide detailed body pressure mapping and flow-field information. It is important to compare predictions between mid-fidelity and high-fidelity codes, both to verify the models and to understand the limitations. A floating marine turbine system was modeled both with OpenFAST and with the commercial CFD code STAR-CCM+. The CFD model used a three-dimensional unsteady Reynolds-averaged Navier-Stokes solver for a volume-of-fluid numerical wave and current tank. The blade-resolved simulations used the sliding-interface technique for the spinning rotor and an overset grid to accommodate the rigid-body motion of the floating system. The mooring system was modeled with a custom coupling of the CFD solver with the open-source code MoorDyn. This improves upon the existing quasi-static catenary solver in STAR-CCM+, which lacks seabed contact or line-to-line connections. Simulation results for a floating marine turbine are compared between OpenFAST and CFD, highlighting the capabilities of the mid-fidelity code and identifying the areas where a high-fidelity approach is needed.

CFD↗
DOE OSTI2025

Geometric Interpretation of the Cluster Location Problem Part II: Application to the Pahala, Hawaii, Earthquake Sequence

In the companion “Theory” article, we presented a new framing of the seismic location problem in terms of differential geometry (Harris et al., 2025). From that viewpoint, we developed a “project and correct” approach for estimating the relative locations of earthquakes. Here, in this study, we use project and correct to estimate high-precision relative locations of events from an earthquake sequence beneath the town of Pahala, Hawaii, using high-precision correlation-derived picks. The sequence was active from 2020 through 2022 and produced many highly correlated signals at Hawaii Volcano Observatory (HVO) stations on the island of Hawaii. The data we inverted consisted of 2882 events with observations at 5 HVO stations. For comparison with the travel-time image, we also produced conventional hypocenter solutions using both the Bayesloc program (Myers et al., 2007, 2009) and a purpose-built double-difference code. There were obvious structural elements in the resulting image, the resolution of which we used to test the performance of the project and the correct algorithm. For the projection step, we first produced a 3D local basis using an singular value decomposition (SVD) of the 2882 groups of times. Projection of the travel-time vectors into this basis resulted in an image with structures similar to those produced by our conventional locators, but with distortion as predicted by theory. Removing the distortion requires an inverse operator generated from the metric tensor at the geometric centroid of the events. We compared two approaches to obtaining such an inverse operator. The first uses an estimate of the geographic centroid of the event cloud from the centroid of the travel-time data. The second approach uses the centroid of the conventionally produced locations. The first approach produces a corrected image very similar to the conventional results, but with a rotation. The corrected image produced using the conventionally derived centroid is a near-exact match to the conventional locations.

Dodge, Douglas A. [Lawrence Livermore National Lab↗
NASA NTRS1994

Computational Fluid Dynamic Modeling of Rocket Based Combined Cycle Engine Flowfields

Computational Fluid Dynamic techniques are used to study the flowfield of a fixed geometry Rocket Based Combined Cycle engine operating in rocket ejector mode. Heat addition resulting from the combustion of injected fuel causes the subsonic engine flow to choke and go supersonic in the slightly divergent combustor-mixer section. Reacting flow computations are undertaken to predict the characteristics of solutions where the heat addition is determined by the flowfield. Here, adaptive gridding is used to improve resolution in the shear layers. Results show that the sonic speed is reached in the unheated portions of the flow first, while the heated portions become supersonic later. Comparison with results from another code show reasonable agreement. The coupled solutions show that the character of the combustion-based thermal choking phenomenon can be controlled reasonably well such that there is opportunity to optimize the length and expansion ratio of the combustor-mixer.

Daines, Russell L.↗
NASA NTRS1995

Stellar models with microscopic diffusion and rotational mixing. 1: Application to the Sun

The Yale stellar evolution code has been modified to include the combined effects of diffusion and rotational mixing on H-1, H-4, and the trace elements He, Li-6, Li-7, and Be-9. The interaction between rotational mixing and diffusion is studied by calculating a number of calibrated solar models. The rotational mixing inhibits the diffusion in the outer parts of the models, leading to a decrease in the envelope diffusion by 25%-50%. Conversely, diffusion leads to gradients in mean molecular weight which can inhibit the rotational mixing. The degree to which gradients in mean molecular weight inhibit the rotational mixing is somewhat uncertain. A comparison with the observed solar oblateness suggests that gradients in the mean molecular weight play a smaller role in inhibiting the rotational mixing than previously believed. This is reinforced by the fact that the model with the standard value for the inhibiting effect of mean molecular weight on the rotational mixing depletes no Li on the main sequence. This is clear in contrast to the observations. A reduction in the inhibiting effect of mean molecular weight gradients by a factor of 10 loads to noticeable main-sequence Li depletion.

Chaboyer, Brian↗