SEARCH · Engineering Papers

Results for “code comparisons”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 757 records · Page 42

All sources NASA NTRS DOE OSTI

DOE OSTI2021

A deterministic verification strategy for electrostatic particle-in-cell algorithms in arbitrary spatial dimensions using the method of manufactured solutions

As simulations of kinetic plasmas continue to increase in scope and complexity, a rigorous and straightforward method for verifying particle-in-cell (PIC) implementations is necessary to ensure their correctness. Here, in this paper, we present a deterministic method for the rigorous verification of multidimensional, multispecies, electrostatic particle-in-cell codes based on the method of manufactured solutions. Specifically, we prove that rigorous verification is possible through the exclusive examination of errors of grid quantities (i.e., moments and/or fields), allowing for a very light-weight and non-intrusive implementation in existing PIC codes. This is a marked improvement over earlier PIC verification approaches (only demonstrated with one species in 1D-1V), which rely on the comparison of cumulative distribution functions, and are invasive on the code base, introduce additional stochastic noise, are significantly more computationally expensive, and lack rigorous convergence properties. Interestingly, we show that different grid quantities feature different rates of convergence with the number of particles and mesh size, impacting the verification process. These theoretical results are confirmed numerically with a multi-species 2D-2V particle-in-cell code, including the ability of the method to detect order reduction due to an incorrect implementation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1985

New integration techniques for chemical kinetic rate equations. 2: Accuracy comparison

A comparison of the accuracy of several techniques recently developed for solving stiff differential equations is presented. The techniques examined include two general purpose codes EEPISODE and LSODE developed for an arbitrary system of ordinary differential equations, and three specialized codes CHEMEQ, CREKID, and GCKP84 developed specifically to solve chemical kinetic rate equations. The accuracy comparisons are made by applying these solution procedures to two practical combustion kinetics problems. Both problems describe adiabatic, homogeneous, gas phase chemical reactions at constant pressure, and include all three combustion regimes: induction, heat release, and equilibration. The comparisons show that LSODE is the most efficient code - in the sense that it requires the least computational work to attain a specified accuracy level. An important finding is that an iterative solution of the algebraic enthalpy conservation equation for the temperature can be more accurate and efficient than computing the temperature by integrating its time derivative.

Radhakrishnan, K.↗

NASA NTRS1986

New integration techniques for chemical kinetic rate equations. II - Accuracy comparison

A comparison of the accuracy of several techniques recently developed for solving stiff differential equations is presented. The techniques examined include two general purpose codes EEPISODE and LSODE developed for an arbitrary system of ordinary differential equations, and three specialized codes CHEMEQ, CREKID, and GCKP84 developed specifically to solve chemical kinetic rate equations. The accuracy comparisons are made by applying these solution procedures to two practical combustion kinetics problems. Both problems describe adiabatic, homogeneous, gas phase chemical reactions at constant pressure, and include all three combustion regimes: induction heat release, and equilibration. The comparisons show that LSODE is the most efficient code - in the sense that it requires the least computational work to attain a specified accuracy level. An important finding is that an iterative solution of the algebraic enthalpy conservation equation for the temperature can be more accurate and efficient than computing the temperature by integrating its time derivative.

Radhakrishnan, K.↗

NASA NTRS2017

Comparing HZETRN, SHIELD, FLUKA and GEANT Transport Codes

For the first time, the American (NASA) and Russian (ROSCOSMOS) space radiation transport codes, HZETRN and SHIELD respectively, are directly compared to each other. Calculations are presented for Galactic Cosmic Ray (GCR) minimum Hydrogen, Oxygen and Iron projectiles incident on a uniform Aluminum cylinder of varying thickness. Comparisons are made for the flux spectra of neutrons, light ions, heavy ions and pions emitted from the back of the Aluminum cylinder. In order to provide more benchmark comparisons, some calculations with the GEANT and FLUKA transport codes are also shown.

John W Norbury↗

NASA NTRS1989

Calculation of steady and unsteady pressures at supersonic speeds with CAP-TSD

A finite difference technique is used to solve the transonic small disturbance flow equation making use of shock capturing to treat wave discontinuities. Thus the nonlinear effects of thickness and angle of attack are considered. Such an approach is made feasible by the development of a new code called CAP-TSD (Computational Aeroelasticity Program - Transonic Small Disturbance), and is based on a fully implicit approximate factorization (AF) finite difference method to solve the time dependent transonic small disturbance equation. The application of the CAP-TSD code to the calculation of low to moderate supersonic steady and unsteady flows is presented. In particular, comparisons with exact linear theory solutions are made for steady and unsteady cases to evaluate shock capturing and other features of the current method. In addition, steady solutions obtained from an Euler code are used to evaluate the small disturbance aspects of the code. Steady and unsteady pressure comparisons are made with measurements for an F-15 wing model and for the RAE tailplane model.

Bennett, Robert M.↗

NASA NTRS2008

Implementation of a Blowing Boundary Condition in the LAURA Code

Preliminary steps toward modeling a coupled ablation problem using a finite-volume Navier-Stokes code (LAURA) are presented in this paper. Implementation of a surface boundary condition with mass transfer (blowing) is described followed by verification and validation through comparisons with analytic results and experimental data. Application of the code to a carbon-nosetip ablation problem is demonstrated and the results are compared with previously published data. It is concluded that the code and coupled procedure are suitable to support further ablation analyses and studies.

Thompson, Richard a.↗

DOE OSTI2025

Demonstration of Monte Carlo Reference Equilibrium Core Solutions for Pebble Bed Reactors Using Shift in kugelpy

This report describes how to produce reference solutions for the equilibrium core search of pebble bed reactors using the pyrates Python library. The pyrates library uses the kugelpy methodology to perform running-in calculations to reach an equilibrium core. This approach relies on full-core Monte Carlo neutron transport calculations and the shifting of fuel pebble inventory through the reactor to simulate pebble movement in the core from one step to the next of a running-in scenario. In addition to the use of Serpent for the Monte Carlo neutron transport calculations, as part of this work, updates for the use of the Shift Monte Carlo code were integrated into the official pyrates GitHub repository. Comparisons of running-in simulations for a generic pebble bed reactor produced by using pyrates with Serpent and Shift are shown. Although consistent trends are shown between the use of the two Monte Carlo codes, the limiting factor in either calculation is the computation time due to the need to perform several hundred full-core neutron transport calculations before an equilibrium core is reached.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2025

Demonstration of Monte Carlo Reference Equilibrium Core Solutions for Pebble Bed Reactors Using Shift in kugelpy

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRSDate not supplied

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↗

DOE OSTI2025

Simulations of edge and SOL turbulence in diverted negative and positive triangularity plasmas

Optimizing the performance of magnetic confinement fusion devices is critical to achieving an attractive fusion reactor design. Negative triangularity (NT) scenarios have been shown to achieve excellent levels of energy confinement, while avoiding edge localized modes. Modeling turbulent transport in the edge and SOL is key in understanding the impact of NT on turbulence and extrapolating the results to future devices and regimes. Previous gyrokinetic turbulence studies have reported beneficial effects of NT across a broad range of parameters. However, most simulations have focused on the inner plasma region, neglecting the impact of NT on the outermost edge. In this work, we investigate the effect of NT in edge and scrape-off layer simulations, including the magnetic X-point and separatrix. For the first time, we employ a multi-fidelity approach, combining global, non-linear gyrokinetic simulations with drift-reduced fluid simulations, to gain a deeper understanding of the underlying physics at play. First-principles simulations using the GENE-X code demonstrate that in comparable NT and PT geometries, similar profiles are achieved, while the turbulent heat flux is reduced by more than 50% in NT. Comparisons with results from the drift-reduced fluid turbulence code GRILLIX suggest that the turbulence is driven by trapped electron modes. The parallel heat flux width on the divertor targets is reduced in NT, primarily due to a lower spreading factor S.

GENE-X↗

NASA NTRS1994

Design oriented structural analysis

Desirable characteristics and benefits of design oriented analysis methods are described and illustrated by presenting a synoptic description of the development and uses of the Equivalent Laminated Plate Solution (ELAPS) computer code. ELAPS is a design oriented structural analysis method which is intended for use in the early design of aircraft wing structures. Model preparation is minimized by using a few large plate segments to model the wing box structure. Computational efficiency is achieved by using a limited number of global displacement functions that encompass all segments over the wing planform. Coupling with other codes is facilitated since the output quantities such as deflections and stresses are calculated as continuous functions over the plate segments. Various aspects of the ELAPS development are discussed including the analytical formulation, verification of results by comparison with finite element analysis results, coupling with other codes, and calculation of sensitivity derivatives. The effectiveness of ELAPS for multidisciplinary design application is illustrated by describing its use in design studies of high speed civil transport wing structures.

Giles, Gary L.↗

NASA NTRS1989

Investigation of the transient fuel preburner manifold and combustor

A computational fluid dynamics (CFD) model with finite rate reactions, FDNS, was developed to study the start transient of the Space Shuttle Main Engine (SSME) fuel preburner (FPB). FDNS is a time accurate, pressure based CFD code. An upwind scheme was employed for spatial discretization. The upwind scheme was based on second and fourth order central differencing with adaptive artificial dissipation. A state of the art two-equation k-epsilon (T) turbulence model was employed for the turbulence calculation. A Pade' Rational Solution (PARASOL) chemistry algorithm was coupled with the point implicit procedure. FDNS was benchmarked with three well documented experiments: a confined swirling coaxial jet, a non-reactive ramjet dump combustor, and a reactive ramjet dump combustor. Excellent comparisons were obtained for the benchmark cases. The code was then used to study the start transient of an axisymmetric SSME fuel preburner. Predicted transient operation of the preburner agrees well with experiment. Furthermore, it was also found that an appreciable amount of unburned oxygen entered the turbine stages.

Wang, Ten-See↗

NASA NTRS1992

Stagnation-point heat-transfer rate predictions at aeroassist flight conditions

The results are presented for the stagnation-point heat-transfer rates used in the design process of the Aeroassist Flight Experiment (AFE) vehicle over its entire aeropass trajectory. The prediction methods used in this investigation demonstrate the application of computational fluid dynamics (CFD) techniques to a wide range of flight conditions and their usefulness in a design process. The heating rates were computed by a viscous-shock-layer (VSL) code at the lower altitudes and by a Navier-Stokes (N-S) code for the higher altitude cases. For both methods, finite-rate chemically reacting gas was considered, and a temperature-dependent wall-catalysis model was used. The wall temperature for each case was assumed to be radiative equilibrium temperature, based on total heating. The radiative heating was estimated by using a correlation equation. Wall slip was included in the N-S calculation method, and this method implicitly accounts for shock slip. The N-S/VSL combination of projection methods was established by comparison with the published benchmark flow-field code LAURA results at lower altitudes, and the direct simulation Monte Carlo results at higher altitude cases. To obtain the design heating rate over the entire forward face of the vehicle, a boundary-layer method (BLIMP code) that employs reacting chemistry and surface catalysis was used. The ratio of the VSL or N-S method prediction to that obtained from the boundary-layer method code at the stagnation point is used to define an adjustment factor, which accounts for the errors involved in using the boundary-layer method.

Gupta, Roop N.↗

NASA NTRS2003

Large-Eddy Simulation of Conductive Flows at Low Magnetic Reynolds Number

In this paper we study the LES method with dynamic procedure in the context of conductive flows subject to an applied external magnetic field at low magnetic Reynolds number R(sub m). These kind of flows are encountered in many industrial applications. For example, in the steel industry, applied magnetic fields can be used to damp turbulence in the casting process. In nuclear fusion devices (Tokamaks), liquid-lithium flows are used as coolant blankets and interact with the surrounding magnetic field that drives and confines the fusion plasma. Also, in experimental facilities investigating the dynamo effect, the flow consists of liquid-sodium for which the Prandtl number and, as a consequence, the magnetic Reynolds number is low. Our attention is focused here on the case of homogeneous (initially isotropic) decaying turbulence. The numerical simulations performed mimic the thought experiment described in Moffatt in which an initially homogeneous isotropic conductive flow is suddenly subjected to an applied magnetic field and freely decays without any forcing. Note that this flow was first studied numerically by Schumann. It is well known that in that case, extra damping of turbulence occurs due to the Joule effect and that the flow tends to become progressively independent of the coordinate along the direction of the magnetic field. Our comparison of filtered direct numerical simulation (DNS) predictions and LES predictions show that the dynamic Smagorinsky model enables one to capture successfully the flow with LES, and that it automatically incorporates the effect of the magnetic field on the turbulence. Our paper is organized as follows. In the next section we summarize the LES approach in the case of MHD turbulence at low R(sub m) and recall the definition of the dynamic Smagorinsky model. In Sec. 3 we describe the parameters of the numerical experiments performed and the code used. Section 4 is devoted to the comparison of filtered DNS results and LES results. Conclusions are presented in Sec. 5.

Knaepen, B.↗

NASA NTRS1983

Numerical computation of transonic flow about wing-fuselage configurations on a vector computer

The transonic wing analysis code TWING, which uses the AF2 relaxation algorithm, has been vectorized to run on the Cray-1S computer. Vectorization of this code improved computational efficiency over that of the CDC 7600 computer by factors of 11 to 13. The improvement compares favorably with the prediction of a theoretical performance model. A convenient generalization now permits the treatment of rudimentary wing-fuselage combinations. Flow predictions for a transport configuration in both isolated-wing and wing-fuselage modes show the expected trends in shock strength and position when compared with wind-tunnel results. An isolated fighter wing is examined in terms of execution time on three different computers and in comparison with experimental data. The computational fluid dynamics code produced during this study is a careful union of an efficient three-dimensional, transonic, numerical algorithm and the vector features presently available on modern computers.

Thomas, S. D.↗

NASA NTRS1995

Ocean observations with EOS/MODIS: Algorithm development and post launch studies

An investigation of the influence of stratospheric aerosol on the performance of the atmospheric correction algorithm was carried out. The results indicate how the performance of the algorithm is degraded if the stratospheric aerosol is ignored. Use of the MODIS 1380 nm band to effect a correction for stratospheric aerosols was also studied. The development of a multi-layer Monte Carlo radiative transfer code that includes polarization by molecular and aerosol scattering and wind-induced sea surface roughness has been completed. Comparison tests with an existing two-layer successive order of scattering code suggests that both codes are capable of producing top-of-atmosphere radiances with errors usually less than 0.1 percent. An initial set of simulations to study the effects of ignoring the polarization of the the ocean-atmosphere light field, in both the development of the atmospheric correction algorithm and the generation of the lookup tables used for operation of the algorithm, have been completed. An algorithm was developed that can be used to invert the radiance exiting the top and bottom of the atmosphere to yield the columnar optical properties of the atmospheric aerosol under clear sky conditions over the ocean, for aerosol optical thicknesses as large as 2. The algorithm is capable of retrievals with such large optical thicknesses because all significant orders of multiple scattering are included.

Gordon, Howard R.↗

NASA NTRS1996

Computation of Thermally Perfect Compressible Flow Properties

A set of compressible flow relations for a thermally perfect, calorically imperfect gas are derived for a value of c(sub p) (specific heat at constant pressure) expressed as a polynomial function of temperature and developed into a computer program, referred to as the Thermally Perfect Gas (TPG) code. The code is available free from the NASA Langley Software Server at URL http://www.larc.nasa.gov/LSS. The code produces tables of compressible flow properties similar to those found in NACA Report 1135. Unlike the NACA Report 1135 tables which are valid only in the calorically perfect temperature regime the TPG code results are also valid in the thermally perfect, calorically imperfect temperature regime, giving the TPG code a considerably larger range of temperature application. Accuracy of the TPG code in the calorically perfect and in the thermally perfect, calorically imperfect temperature regimes are verified by comparisons with the methods of NACA Report 1135. The advantages of the TPG code compared to the thermally perfect, calorically imperfect method of NACA Report 1135 are its applicability to any type of gas (monatomic, diatomic, triatomic, or polyatomic) or any specified mixture of gases, ease-of-use, and tabulated results.

Witte, David W.↗

NASA NTRS1983

Numerical solution to the glancing sidewall oblique shock wave/turbulent boundary layer interaction in three dimension

A supersonic three-dimensional viscous forward-marching computer design code called PEPSIS is used to obtain a numerical solution of the three-dimensional problem of the interaction of a glancing sidewall oblique shock wave and a turbulent boundary layer. Very good results are obtained for a test case that was run to investigate the use of the wall-function boundary-condition approximation for a highly complex three-dimensional shock-boundary layer interaction. Two additional test cases (coarse mesh and medium mesh) are run to examine the question of near-wall resolution when no-slip boundary conditions are applied. A comparison with experimental data shows that the PEPSIS code gives excellent results in general and is practical for three-dimensional supersonic inlet calculations.

Anderson, B. H.↗