Engineering Papers⌕ Search

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.

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

At least 541 records · Page 30

Verification and Validation of the STAT7 Code

This document presents the verification and validation of the STAT7 Version 1.1 code, which has been developed to perform the steady-state thermal hydraulics (TH) safety calculations required to support conversion of the Massachusetts Institute of Technology Research Reactor from highly enriched uranium fuel to low- enriched uranium fuel. The key capabilities of the code, including statistical processing, water property generation, and TH solution, are verified independently. While the first two capabilities are checked using statistical tests and NIST data tables, each component of the TH solution is verified using simple geometry test problems. The comparison between results obtained by the code and by hand calculations on these problems shows that each of these components works correctly. Additionally, a comprehensive comparison between STAT7 Version 1.1 and PLTEMP/ANL Version 4.3 Revision 162 calculations is performed for the validation. Good agreement is achieved between the results obtained with the two codes for all TH parameters belonging to the fuel core region.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Performance Characteristics of a Boundary Layer Ingesting Inlet - Distortion Tolerant Fan Through CFD Code Coupling

A method of constructing a high-fidelity simulation using two NASA CFD codes, Launch Ascent and Vehicle Aerodynamics (LAVA), which is the general-purpose solver, and a turbomachinery specific code NASA-TURBO (TURBO), will be presented. Validation of the code coupling is performed through comparison to experimental data obtained in the 2017 wind tunnel test of NASA’s Boundary Layer Ingesting Inlet - Distortion Tolerant Fan (BLI2DTF). For the simulation, the external flow field is computed by LAVA utilizing a steady-state Reynolds-averaged Navier-Stokes (RANS) methodology. Data from LAVA is exchanged with TURBO, which simulated the propulsor utilizing an unsteady-RANS methodology. A sensitivity analysis is performed using averaged and instantaneous data transfer between the two codes, as well as varying the amount of simulation time TURBO computes between updates. The results of the coupled simulations appear to replicate the BLI2DTF experiment correctly. Static pressure taps leading into and around the nacelle, total pressure ratio and efficiency of the propulsor, and profiles of total temperature, pressure, and Mach number at the rear measurement plane, aft of the propulsor, are in good agreement with experimental data. These results provide confidence that the coupled simulation is behaving in a physically realizable manner.

CFD↗

Computational analysis for an advanced transport configuration with engine nacelle

A small-disturbance transonic analysis code is used to calculate the flow-field effects of adding an engine nacelle to a wing/body configuration. Analyses are performed on an advanced transport configuration with and without engine nacelles. Two nacelle shapes are analyzed and the effects of the nacelle installation on pressure distributions are compared with experimental results obtained by shifting the nacelle longitudinally and vertically relative to the wing. Effects of varying the nacelle installation yaw angle are also analyzed and compared with experimental data. These comparisons show that the analysis code is adequately sensitive to variations in nacelle shape, longitudinal and vertical location beneath the wing, and the nacelle installation yaw angle. Results indicate that the code can be used as an effective guide during the design process.

Waggoner, E. G.↗

LSENS: A General Chemical Kinetics and Sensitivity Analysis Code for homogeneous gas-phase reactions. Part 1: Theory and numerical solution procedures

LSENS, the Lewis General Chemical Kinetics and Sensitivity Analysis Code, has been developed for solving complex, homogeneous, gas-phase chemical kinetics problems and contains sensitivity analysis for a variety of problems, including nonisothermal situations. This report is part 1 of a series of three reference publications that describe LENS, provide a detailed guide to its usage, and present many example problems. Part 1 derives the governing equations and describes the numerical solution procedures for the types of problems that can be solved. The accuracy and efficiency of LSENS are examined by means of various test problems, and comparisons with other methods and codes are presented. LSENS is a flexible, convenient, accurate, and efficient solver for chemical reaction problems such as static system; steady, one-dimensional, inviscid flow; reaction behind incident shock wave, including boundary layer correction; and perfectly stirred (highly backmixed) reactor. In addition, the chemical equilibrium state can be computed for the following assigned states: temperature and pressure, enthalpy and pressure, temperature and volume, and internal energy and volume. For static problems the code computes the sensitivity coefficients of the dependent variables and their temporal derivatives with respect to the initial values of the dependent variables and/or the three rate coefficient parameters of the chemical reactions.

Radhakrishnan, Krishnan↗

Test Generator for MATLAB Simulations

MATLAB Automated Test Tool, version 3.0 (MATT 3.0) is a software package that provides automated tools that reduce the time needed for extensive testing of simulation models that have been constructed in the MATLAB programming language by use of the Simulink and Real-Time Workshop programs. MATT 3.0 runs on top of the MATLAB engine application-program interface to communicate with the Simulink engine. MATT 3.0 automatically generates source code from the models, generates custom input data for testing both the models and the source code, and generates graphs and other presentations that facilitate comparison of the outputs of the models and the source code for the same input data. Context-sensitive and fully searchable help is provided in HyperText Markup Language (HTML) format.

Henry, Joel↗

SpaceX Dragon 1 Post Flight MMOD Inspection Campaign

Introduction The SpaceX Dragon 1 spacecraft performed 20 resupply missions to the International Space Station (ISS) between 2012 and 2020. After each mission, a team from the NASA John-son Space Center Hypervelocity Impact Technology Group inspected each Dragon 1 capsule for hypervelocity impact damage features. Data from these inspections are collected into a database that includes impact feature dimensions as well as the location on the vehicle. Additional details on the type and size of particle that produced the damage site are provided when sampling data and definitive spectroscopic analysis results are available. Observation data can be compared with impact estimates from risk assessment codes as a check on the micrometeoroid and orbital debris (MMOD) environment predictions. Scope A general description of the areas of the vehicle that were inspected are provided as well as mission details such as exposure duration and launch dates. The paper documents the general inspection procedure for collection of data and the post inspection data analysis process. It also provides details of the observation data collected as well as the results of analysis of intact samples collected for spectroscopic analysis to discern the source of the impacting particle. A comparison between observed impacts and the expected number of damage features calculated by Bumper 3 with the latest MMOD environments are also presented. Findings Statistics on the >300 impact features documented in the database will provide insight into the depth to diameter ratios and other relationships. The quality of the comparison between the observations and code predictions are dependent on several factors. The paper provides de-tails of each of these variables. (1) Damage equations (2) Impact condition assumptions a. Projectile density b. Impact speed c. Impact angle (3) Analysis results a. pre-flight vs. on-orbit damage b. MMOD vs. non-MMOD Conclusions and Recommendations The ISS visiting vehicle impact database is an ongoing project. SpaceX provides crew rotations as well as resupply and cargo return with the Dragon 2 spacecraft. Sierra Space is manifested for resupply and cargo return services starting in 2023. Boeing is expected to provide crew rotations as well. All of these spacecraft will provide additional opportunities for post flight MMOD inspections of space exposed hardware. Acknowledgements Over the course of the inspection campaign various personnel at the SpaceX Texas Test Site in McGregor and at the SpaceX headquarters in Hawthorne CA provided valuable support to this activity.

James L Hyde↗

flash -light on the ring : hydrodynamic simulations of expanding supernova shells near supermassive black holes

ABSTRACT The way supermassive black holes (SMBHs) in Galactic Centres (GCs) accumulate their mass is not completely determined. At large scales, it is governed by galactic encounters, mass inflows connected to spirals arms and bars, or due to expanding shells from supernova (SN) explosions in the central parts of galaxies. The investigation of the latter process requires an extensive set of gas dynamical simulations to explore the multidimensional parameter space needed to frame the phenomenon. The aims of this paper are to extend our investigation of the importance of SNe for inducing accretion on to an SMBH and carry out a comparison between the fully hydrodynamic code flash and the much less computationally intensive code ring, which uses the thin shell approximation. We simulate 3D expanding shells in a gravitational potential similar to that of the GC with a variety of homogeneous and turbulent environments. In homogeneous media, we find convincing agreement between flash and ring in the shapes of shells and their equivalent radii throughout their whole evolution until they become subsonic. In highly inhomogeneous, turbulent media, there is also a good agreement of shapes and sizes of shells, and of the times of their first contact with the central 1-pc sphere, where we assume that they join the accretion flow. The comparison supports the proposition that an SN occurring at a galactocentric distance of 5 pc typically drives 1–3 M⊙ into the central 1 pc around the GC.

79 ASTRONOMY AND ASTROPHYSICS↗

Navier-Stokes computations and experimental comparisons for multielement airfoil configurations

A two-dimensional unstructured Navier-Stokes code is utilized for computing the flow around multielement configurations. Comparisons are shown for a landing configuration with an advanced-technology flap as well as for a takeoff configuration. For each, a grid convergence study is conducted to assess inaccuracies caused by inadequate grid resolution. Although adequate resolution is obtained for determining the pressure distributions, further refinement is needed to sufficiently resolve the velocity profiles at high angles of attack. For the advanced flap configuration, comparisons of pressure distributions and lift are made with experimental data. Here, two flap riggings and two Reynolds numbers are considered. In general, the trends caused by variations in these quantities are well predicted by the computations, although the angle of attack for maximum lift is overpredicted. For the takeoff configuration, comparisons between pressures, velocity profiles, and Reynolds stresses are made at a single angle of attack with generally good agreement.

Anderson, W. K.↗

Comparison of Temperature-Dependent Cross Section Treatments Within the Shift Monte Carlo Radiation Transport Code

Multiple calculation methods have been developed to account for Doppler broadening of cross sections in Monte Carlo radiation transport simulations. This work compares the accuracy and performance of four different broadening methods, including an on-the-fly multipole treatment, applied to three different geometries of increasing complexity, including a nuclear reactor assembly with a realistic temperature profile. The comparisons reveal how the combination of physics implementations in the Shift Monte Carlo code result in slightly different temperature-dependent behavior, and different levels of performance, for a variety of materials and situations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Terry Turbopump Expanded Operating Band Modeling and Simulation Efforts in Fiscal Year 2021 Extended Period of Performance (Final Report)

This report documents the progress made under the Terry Turbine Expanded Operating Band (TTEXOB) program's modeling and simulation (MODSIM) initiative at Sandia National Laboratories (SNL ). It describes the US Federal Fiscal Year 2021 (FY21) extended period-of-performance MODSIM work completed since the closure of FY20 with due reference to the Texas A&M University (TAMU) hybrid milestone 5/6 experimental program. This work, which falls under Milestone 7 of the program, provides a counterpart to the various experiments. The overall TTEXOB program and its milestone-based approach are described in the program's Summary Plan. Details of the individual milestone test plans can be found in the corresponding detailed test plan, e.g. the Milestone 3 and 4 Detailed Test Plan. SNL MODISM is conducted alongside experiments performed at TAMU, and SNL technical staff regularly consults with TAMU on the experimental program. In FY21, MELCOR code models and capabilities were exercised in two different contexts: experimental comparisons to the TAMU ZS-1 and GS-2, and stand-alone analyses of a station black-out (SBO) scenario in a generic boiling water reactor (BWR). Code to experiment comparisons met with fair success when turbine losses were well characterized as for the ZS-1 turbine. Both deterministic and Bayesian calibration processes were used to find a recommended turbine torque multiplier for ZS-1 type turbines. This process could be repeated for GS-2 type turbines if GS-2 losses were better understood. Stand-alone generic BWR SBO calculations revealed that three different modes of self-regulating turbopump behavior may be observed depending on certain modeling parameters and choices having to do with turbine nozzles. Aspects of this predicted behavior may have been observed in TAMU GS-2 experiments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Comparison of rate one-half, equivalent constraint length 24, binary convolutional codes for use with sequential decoding on the deep-space channel

Virtually all previously-suggested rate 1/2 binary convolutional codes with KE = 24 are compared. Their distance properties are given; and their performance, both in computation and in error probability, with sequential decoding on the deep-space channel is determined by simulation. Recommendations are made both for the choice of a specific KE = 24 code as well as for codes to be included in future coding standards for the deep-space channel. A new result given in this report is a method for determining the statistical significance of error probability data when the error probability is so small that it is not feasible to perform enough decoding simulations to obtain more than a very small number of decoding errors.

Massey, J. L.↗

Development of a CFD code for internal flows in liquid fueled engines

To support the design efforts of engines in which liquid propellants are used, one is often required to analyze incompressible, two-dimensional or axisymmetric flows within ducts and cavities with rotating walls and complicated geometries. The steady-state solution is of interest in most cases. This code is intended to provide a tool for efficient computational fluid dynamics (CFD) analysis of such flow problems, taking advantage of the artificial compressibility concept and the Beam-Warming numerical scheme modified for second order accurate, implicit boundary conditions. These concepts ensure a stable, robust, and accurate algorithm due to the reduced speed of sound and the accuracy of the boundary conditions. The code is dedicated only to two-dimensional or axisymmetric flows with or without swirl. Three-dimensional computation is excluded to increase efficiency and speed. This paper briefly presents the theory of the code, as well as several benchmark applications with comparison to well known analytical solutions. In all these test cases, the code produced accurate results.

Dakhoul, Youssef↗

Comparison of SPHC Hydrocode Results with Penetration Equations and Results of Other Codes

The SPHC hydrodynamic code was used to simulate impacts of spherical aluminum projectiles on a single-wall aluminum plate and on a generic Whipple shield. Simulations were carried out in two and three dimensions. Projectile speeds ranged from 2 kilometers per second to 10 kilometers per second for the single-wall runs, and from 3 kilometers per second to 40 kilometers per second for the Whipple shield runs. Spallation limit results of the single-wall simulations are compared with predictions from five standard penetration equations, and are shown to fall comfortably within the envelope of these analytical relations. Ballistic limit results of the Whipple shield simulations are compared with results from the AUTODYN-2D and PAM-SHOCK-3D codes presented in a paper at the Hypervelocity Impact Symposium 2000 and the Christiansen formulation of 2003.

Evans, Steven W.↗

Modelling the nonlinear plasma response to externally applied three-dimensional fields with the Stepped Pressure Equilibrium Code

Small-amplitude, symmetry-breaking magnetic field perturbations, including resonant magnetic perturbations (RMPs) and error fields, can profoundly impact plasma properties in both tokamaks and stellarators. In this work, we perform the first comparison between the Stepped Pressure Equilibrium Code (SPEC) (a comparatively fast and efficient equilibrium code based on energy-minimisation principles) and M3D-C 1 (a high-fidelity albeit computationally expensive initial-value extended-magnetohydro-dynamic (MHD) code) to assess the conditions under which SPEC can be used to model the nonlinear, non-ideal plasma response to an externally applied (m = 2,n = 1) RMP field in an experimentally relevant geometry. We find that SPEC is able to capture the plasma response in the weakly nonlinear regime – meaning perturbation amplitudes below the threshold for break up of the separatrix and onset of secondary magnetic island formation – when around half of the total toroidal flux is enclosed in the volume containing the q = 2 resonant surface. Further, the observed dependence of SPEC solutions on input parameters, including toroidal flux and the number of volumes into which the plasma is partitioned, indicates that additional exploration of the underlying Multi-Region Relaxed MHD physics model is needed to constrain the choice of parameters. Nonetheless, this work suggests promising applications of SPEC to optimisation and fusion plasma design.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

An Experimental Comparison of Cyclic Variations in Diesel–Natural Gas and POMDME–Natural Gas Dual Fuel Combustion

Abstract Cyclic variations in internal combustion engines are caused by various factors, including combustion mixture stratification, in-cylinder flows, local fluctuations in air-fuel ratio, etc. Cyclic variations have a profound impact on engine performance and emissions. In this study, cyclic variations in dual fuel combustion are analyzed, comparing diesel–natural gas (NG) and polyoxymethylene dimethyl ether (POMDME)-NG dual fuel combustion. Cyclic variability was initially quantified using the coefficient of variation of gross indicated mean effective pressure (IMEPg) computed from experimental cylinder pressure data. The cases analyzed in this study had a coefficient of variation (COV) of IMEPg greater than or around 5%, which was the lower limit of onset of instability for this engine. Experiments were performed at two fixed start of injection (SOI) of high-cetane fuel: 310 CAD and 350 CAD. For all experiments, a constant load of 5 bar IMEPg was maintained, and the intake boost pressure and rail pressure were fixed at 1.5 bar and 500 bar, respectively. For each case, 1000 cycles of cylinder pressure data were recorded, filtered, and processed using an in-house heat release analysis code for each cycle. A comparison between individual cycles and the “ensemble averaged cycle” was made for both diesel–NG and POMDME–NG combustion. For the early SOI of 310 CAD, the peak cylinder pressure fluctuations of individual cycle were found to be ± 15 bar for both fuel combinations, compared to the ensemble averaged cycle, and < 1/10th of the cycles had an IMEPg lower than 0.05 bar of the ensemble averaged cycle. However, the peak pressure fluctuations were found to be lower for POMDME–NG (±3 bar) than diesel–NG dual fuel combustion at 350 CAD SOI, indicating lower cyclic variations. The higher reactivity of POMDME helped reduce fluctuations in combustion phasing at the retarded SOI. The presence of cycles of deterioration and cycles of recovery were also observed with diesel–NG combustion for 310 CAD SOI, and the scatter in the IMEPg return map was similar for both fuel combinations. The IMEPg return map for POMDME–NG combustion was less scattered at the 350 CAD SOI.

Engineering↗

Verification and Validation of the STAT7 Code

This document presents the verification and validation of the STAT7 Version 1.0 code which has been developed to perform the steady-state thermal hydraulic safety calculation of the MIT Research Reactor required upon conversion of the reactor from high enrichment uranium fuel to low enrichment uranium fuel. The key capabilites of the code including the statistical processing, water property generation, and thermal hydraulic solution are verified independently. While the first two capabilites are checked using statistical tests and NIST data tables, each component of the thermal hydraulic solution is verified using simple geometry test problems. The comparison between results obtained by the code and hand calculations on these problems shows that each of these components work correctly. Additionally, a comprehensive comparison between STAT7 Version 1.0 and PLTEMP/ANL Version 4.2 calculations for a five-channel test problem is performed for the validation. The good agreement between the obtained results for all thermal hydraulic parameters belonging to the fuel core region is achieved. An issue in bypass flowrate calculation was detected when both fin-friction and bypass flow claculations are involved. This problem can be corrected in a later version of the code, however, does not affect the application of STAT7 code as it has been used to date where the flow distribution in the core region and bypass channel is specified as an input parameter.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Verification and Validation of the STAT7 Code

This document presents the verification and validation of the STAT7 Version 1.0 code which has been developed to perform the steady-state thermal hydraulic safety calculation of the MIT Research Reactor required upon conversion of the reactor from highly enriched uranium fuel to low- enriched uranium fuel. The key capabilities of the code including the statistical processing, water property generation, and thermal hydraulic solution are verified independently. While the first two capabilities are checked using statistical tests and NIST data tables, each component of the thermal hydraulic solution is verified using simple geometry test problems. The comparison between results obtained by the code and hand calculations on these problems shows that each of these components work correctly. Additionally, a comprehensive comparison between STAT7 Version 1.0 and PLTEMP/ANL Version 4.2 calculations for a five-channel test problem is performed for the validation. The good agreement between the obtained results for all thermal hydraulic parameters belonging to the fuel core region is achieved. An issue in bypass flowrate calculation was detected when both fin-friction and bypass flow calculations are involved. This problem can be corrected in a later version of the code. However, it does not affect the application of STAT7 as it has been used to date where the flow distribution in the core region and bypass channel is specified as an input parameter.

97 MATHEMATICS AND COMPUTING↗

Implementation and validation of a wake model for vortex-surface interactions in low speed forward flight

The wake model was implemented using a VAX 750 and a Microvax II workstation. Online graphics capability using a DISSPLA graphics package. The rotor model used by Beddoes was significantly extended to include azimuthal variations due to forward flight and a simplified scheme for locating critical points where vortex elements are placed. A test case was obtained for validation of the predictions of induced velocity. Comparison of the results indicates that the code requires some more features before satisfactory predictions can be made over the whole rotor disk. Specifically, shed vorticity due to the azimuthal variation of blade loading must be incorporated into the model. Interactions between vortices shed from the four blades of the model rotor must be included. The Scully code for calculating the velocity field is being modified in parallel with these efforts to enable comparison with experimental data. To date, some comparisons with flow visualization data obtained at Georgia Tech were performed and show good agreement for the isolated rotor case. Comparison of time-resolved velocity data obtained at Georgia Tech also shows good agreement. Modifications are being implemented to enable generation of time-averaged results for comparison with NASA data.

Komerath, Narayanan M.↗