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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Cross-Validation of Computational and Experimental Distributed Surface Pressures on the Space Launch System

This paper presents a new workflow for comparing experimental pressure-sensitive paint (PSP) data to computational fluid dynamic (CFD) simulations by way of mapping data from corresponding grids utilizing interpolation methods. In addition to generating quantitative and qualitative point-to-point comparisons between PSP and CFD data, this workflow extracts sectional loading data from both grids and generates lineload comparison charts for corresponding PSP and CFD runs. Experimental PSP data presented in this paper were taken from a 2016 NASA Ames Research Center Unitary Plan Wind Tunnel 11- by 11-Foot Transonic WindTunnel Facility test of the NASA Space Launch System. CFD simulation data for comparison purposes were generated using the FUN3D code. Overall, interpolation onto PSP grids versus CFD grids yields comparable surface pressure fields. However, lineload comparisons are easier to make on the CFD grid-mapped data due to the grid topology and the current capabilities of the lineload analysis tools at NASA Langley Research Center. This workflow is written using contemporary software (Python, Tecplot, PyTecplot), is compatible with existing tools at NASA Langley, and is developed to be adaptable depending on the situation.

SLS↗
DOE OSTI2024

Code Benchmark of Depressurized Conduction Cooldown Transient in the High Temperature Test Facility

This paper presents results from modeling of a depressurized conduction cooldown (DCC) transient at the High Temperature Test Facility (HTTF) as part of the OECD-NEA Thermal Hydraulics Code Validation Benchmark for High-Temperature Gas-Cooled Reactors using HTTF Data . This paper briefly describes the benchmark and the models being used. It then presents a comparison of steady state and transient results based on the Problem 2 Exercise 1A and 1B definitions. We compare block and helium temperature distributions, mass flow distribution, and energy balance in steady state. All models show comparable mass flow distributions and energy balances. The temperatures within the core and outer regions are comparable in all models too, but inner reflector temperatures can vary significantly. Despite that, we find that the models are in good agreement for the full-power steady state. In the DCC, we look at block temperature at the core midplane and RCCS water exit temperature. The INL and ANL models are found to be in excellent agreement with one another on block temperature over time, while the agreement when the KAERI and NRG models are added into consideration is good. Differences in the transient heat removal from the RCCS cause the differences in block temperature over time in these models. The CNL models show similar trends to the INL, ANL, KAERI, and NRG models, but the temperatures are high because the volumes used in calculating the average temperature include the heater rods in the CNL models only. The HUN-REN model shows results that suggest significantly lower heat removal in the RCCS which merit further investigation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1995

Off-design computer code for calculating the aerodynamic performance of axial-flow fans and compressors

An off-design axial-flow compressor code is presented and is available from COSMIC for predicting the aerodynamic performance maps of fans and compressors. Steady axisymmetric flow is assumed and the aerodynamic solution reduces to solving the two-dimensional flow field in the meridional plane. A streamline curvature method is used for calculating this flow-field outside the blade rows. This code allows for bleed flows and the first five stators can be reset for each rotational speed, capabilities which are necessary for large multistage compressors. The accuracy of the off-design performance predictions depend upon the validity of the flow loss and deviation correlation models. These empirical correlations for the flow loss and deviation are used to model the real flow effects and the off-design code will compute through small reverse flow regions. The input to this off-design code is fully described and a user's example case for a two-stage fan is included with complete input and output data sets. Also, a comparison of the off-design code predictions with experimental data is included which generally shows good agreement.

Schmidt, James F.↗
NASA NTRS1987

High speed simulation of flexible multibody dynamics

A multiflexible body dynamics code intended for fast turnaround control design trades is described. Nonlinear rigid body dynamics and linearized flexible dynamics combine to provide efficient solution of the equations of motion. Comparison with results from the DISCOS code provide verification of accuracy.

Jacot, A. D.↗
DOE OSTI2022

Comparison of measurement and simulation of ATLAS cavern radiation background

Sixteen Medipix2 pixel detector based (MPX) devices were operated at various positions within the ATLAS detector and cavern continuously from early 2008 up to 2013. In addition to photons, each MPX detector is capable to detect charged particles, and neutrons as it is covered with a mask of converter materials dividing its area into regions sensitive to thermal or fast neutrons. The MPX detector network was effectively used for real-time measurements of the spectral characteristics and composition of complex radiation fields in ATLAS. This article reports comparison of the results of measurements performed with MPX detectors during the LHC operation period in 2010 and 2011 with Monte Carlo simulations results from the FLUGG and GCALOR codes. For the purpose of this comparison, the MPX detectors were operated in tracking mode with low threshold (8–10 keV) allowing one to distinguish among particle categories based on the recognition of track patterns left by the particles in the MPX sensitive layer. The comparison of measurements with simulations shows that the agreement between measured and simulated data is satisfactory in most cases within a factor of two.

47 OTHER INSTRUMENTATION↗
NASA NTRS1993

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 airfoil configurations. Comparisons are shown for a landing configuration with an advanced-technology flap. Grid convergence studies are 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.

Anderson, W. Kyle↗
DOE OSTI2026

Evaluation of Sandia NCS Benchmark Suite Updates

The Sandia Nuclear Criticality Safety (NCS) program’s benchmark suite was recently updated. This suite is used to ensure that NCS calculations using computer-based neutron transportation codes have an established baseline comparison of calculated versus known experimental results. The Evaluated Nuclear Data File (ENDF) version used for the MCNP models in the suite was changed from ENDF/B-VII.1 to ENDF/B-VIII.0. Additionally, relevant thermal scattering law data libraries (TSLs) were updated. The sensitivity of the calculational bias of each benchmark model to these changes is discussed. Implementation of the ENDF/B-VIII.0 library and updated TSLs results in improvements to bias distribution in the intermediate enriched uranium, plutonium, and mixed uranium–plutonium (IEU, PU, and MIX) fissionable material benchmark categories, but a small bias increase in low- and high-enriched uranium categories (LEU and HEU, respectively). The results also highlight the sensitivity of the benchmarks, with average lethargy of neutrons causing fission energies (EALF) in the intermediate energy range to ENDF/B library changes. The most numerous bias changes were observed in the thermal energy region when transitioning from ENDF/B-VII.1 to ENDF/B-VIII.0. In conclusion, most of the unique bias changes observed in MCNP 6.3.0 between the two nuclear data libraries were in the LEU-COMP-THERM evaluation subset.

ICSBEP↗
DOE OSTI2021

Feedback density and causal complexity of simulation model structure

Measures of simulation model complexity generally focus on outputs; we propose measuring the complexity of a model’s causal structure to gain insight into its fundamental character. This article introduces tools for measuring causal complexity. First, we introduce a method for developing a model’s causal structure diagram, which characterises the causal interactions present in the code. Causal structure diagrams facilitate comparison of simulation models, including those from different paradigms. Next, we develop metrics for evaluating a model’s causal complexity using its causal structure diagram. We discuss cyclomatic complexity as a measure of the intricacy of causal structure and introduce two new metrics that incorporate the concept of feedback, a fundamental component of causal structure. The first new metric introduced here is feedback density, a measure of the cycle-based interconnectedness of causal structure. The second metric combines cyclomatic complexity and feedback density into a comprehensive causal complexity measure. Finally, we demonstrate these complexity metrics on simulation models from multiple paradigms and discuss potential uses and interpretations. These tools enable direct comparison of models across paradigms and provide a mechanism for measuring and discussing complexity based on a model’s fundamental assumptions and design.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1981

Relations for the thermodynamic and transport properties in the testing environment of the Langley hypersonic CF4 tunnel

Thermodynamic and transport properties of gaseous CF4 that can be used in flow field computer codes and theoretical analyses for comparison of results with experimental data from the Langley Hypersonic CF4 Tunnel are presented. The thermodynamic relations which are based on thermally perfect but calorifically imperfect gas are adequate for the testing environment of the CF4 tunnel and are simpler to use than the imperfect gas relations required to define the complete expansion from the tunnel reservoir. Relations for the transport properties are based on the kinetic theory of gases in which published experimental data are used in the derivation of the relations. Extensive experimental data were located for viscosity and the derived relation should provide values for viscosity with errors of less than 1 percent. The experimental data for thermal conductivity were limited with significant disagreement between the various sources. The derived relation will probably provide values for the thermal conductivity with errors of no more than 5 percent which is within the accuracy of the experimental data.

Sutton, K.↗
NASA NTRS1984

Modeling material failure with a vectorized routine

The computational aspects of modelling material failure in structural wood members are presented with particular reference to vector processing aspects. Wood members are considered to be highly orthotropic, inhomogeneous, and discontinuous due to the complex microstructure of wood material and the presence of natural growth characteristics such as knots, cracks and cross grain in wood members. The simulation of strength behavior of wood members is accomplished through the use of a special purpose finite element/fracture mechanics routine, program STARW (Strength Analysis Routine for Wood). Program STARW employs quadratic finite elements combined with singular crack tip elements in a finite element mesh. Vector processing techniques are employed in mesh generation, stiffness matrix formation, simultaneous equation solution, and material failure calculations. The paper addresses these techniques along with the time and effort requirements needed to convert existing finite element code to a vectorized version. Comparisons in execution time between vectorized and nonvectorized routines are provided.

Cramer, S. M.↗
NASA NTRS1984

Investigation of mixing in a turbofan exhaust duct. II Computer code application and verification

A three-dimensional analysis of turbofan forced mixer nozzle aerodynamics demonstrates that the complex flow structure is dominated by geometrically induced secondary flow rather than by turbulence. The test apparatus consisted of a fixed upstream model section and a rotating shroud. The Mach number of the fan and core streams at the mixing plane (lobe exit) was 0.45, the bypass ratio was about 4, and the Reynolds number based on the shroud radius was 1,100,000. The three velocity components near the exit plane of the lobes were measured using flow angularity probes to provide information about the mixer inflow conditions for turbulent computations. The validity of a previous computer code was demonstrated in a comparison of the nozzle exit temperature data with the computed temperature distributions. The mechanism most responsible for the generation of secondary flow within the lobes is due to the turning of the fan and core streams in opposite radial directions.

Povinelli, L. A.↗
NASA NTRS1986

Effect of subcooling on the on-orbit pressurization rate of cryogenic propellant tankage

The SOLA-ECLIPSE code is being developed to enable prediction of the behavior of cryogenic propellants in spacecraft tankage. A brief description of the formulations used for modeling heat transfer and for determining thermodynamic state is presented. Code performance is verified through comparison to experimental data for the self-pressurization of scale model liquid hydrogen tanks. SOLA-ECLIPSE is used to examine the effect of initial subcooling of the liquid phase on the self-pressurization rate of an on-orbit full scale liquid hydrogen tank typical for a chemical propulsion Orbital Transfer Vehicle. The computational predictions show that even small amounts of subcooling will significantly decrease the self-pressurization rate. Further, if the cooling is provided by a Thermodynamic Vent System, it is concluded that small levels of subcooling will maximize propellant conservation.

Hochstein, J. I.↗
NASA NTRS1987

Temperature fields due to jet induced mixing in a typical OTV tank

The Eclipse Code is being developed as a general tool for analysis of cryogenic propellant behavior in spacecraft tankage. The focus of the work being reported is on prediction of temperature fields due to introduction of a cold jet along the centerline of a typical Orbit Transfer Vehicle tank. A brief description of the formulations used for modeling heat transfer and turbulent flow is presented. Code performance is verified through comparison to experimental data for mixing in small scale tanks. An unexpected difficulty in computing long duration flows is reviewed. Preliminary results for a partially filled full scale tank are obtained by approximating the free surface by a spherical solid boundary.

Hochstein, J. I.↗
NASA NTRS1987

Ice measurements by Geosat radar altimetry

Radar altimetry for ice-covered ocean and land is more complex and variable than open ocean radar altimetry; attention is presently given to Geosat ice-sheet topography for the Greenland and Antarctic ice sheets between 72 deg N and 72 deg S which owes its excellent accuracy to the well separated spacing of the orbital tracks and an 18-month geodetic mission duration. A surface elevation map of southern Greenland, produced from 110 days of retracked Geosat data, is presented in color-coded three-dimensional perspective. Comparisons are made between Seasat and Geosat data for ice mass elevations in Greenland.

Zwally, H. Jay↗
NASA NTRS1988

Modeling of impulsive propellant reorientation

The impulsive propellant reorientation process is modeled using the (Energy Calculations for Liquid Propellants in a Space Environment (ECLIPSE) code. A brief description of the process and the computational model is presented. Code validation is documented via comparison to experimentally derived data for small scale tanks. Predictions of reorientation performance are presented for two tanks designed for use in flight experiments and for a proposed full scale OTV tank. A new dimensionless parameter is developed to correlate reorientation performance in geometrically similar tanks. Its success is demonstrated.

Hochstein, John I.↗
NASA NTRS1989

Modeling of impulsive propellant reorientation

The impulsive propellant reorientation process is modeled using the Energy Calculations for Liquid Propellants in a Space Environment (ECLIPSE) code. A brief description of the process and the computational model is presented. Code validation is documented via comparison to experimentally derived data for small scale tanks. Predictions of reorientation performance are presented for two tanks designed for use in flight experiments and for a proposed full scale OTV tank. A new dimensionless parameter is developed to correlate reorientation performance in geometrically similar tanks. Its success is demonstrated.

Hochstein, John I.↗
NASA NTRS1989

Application of a full-potential solver to bending-torsion flutter in cascades

A two-dimensional, unsteady, full-potential cascade aerodynamic code has been used for flutter calculations. The code, which uses a time-marching implicit scheme with internal Newton-iterations at each time step, allows oscillatory blade motions with phase lag between adjacent blades to be modeled using a single inter-blade passage. This aerodynamic code has been used with a two-degrees-of-freedom typical section structural model to determine the aeroelastic stability of a bladed disk. The unsteady pressure results calculated from the aerodynamic code have been verified by comparison with results from the linear and linearized potential theories. Flutter results have been presented for two cases. A comparison with linear theory shows identical trends; however, the present calculations show a lower flutter frequency.

Bakhle, Milind A.↗
NASA NTRS1990

Numerical investigation of laminar-turbulent transition in a flat plate wake

Lamina-turbulent transition of high-deficit flat plate wakes is investigated by direct numerical simulations using the complete Navier Stokes equations. The simulations are based on a spatial model so that both the base flow and the disturbance flow can develop in the downstream direction. The Navier Stokes equations are used in a vorticity-velocity form and are solved using a combination of finite difference and spectral approximations. Fourier series are used in the spanwise direction. Second-order finite-differences are used to approximate the spatial derivatives in the streamwise and transverse directions. For the temporal discretion, a combination of ADI, Crank-Nicolson, and Adams-Bashforth methods is employed. The discretized velocity equations are solved using fast Helmholtz solvers. Code validation is accomplished by comparison of the numerical results to both linear stability and to experiments. Calculations of two- and/or three-dimensional sinuous and mode disturbances in the wake of flat plate are undertaken. For calculations of two-dimensional disturbances, the wake is forced at an amplitude level so that nonlinear disturbance development may be observed. In addition, the forcing amplitude is varied in order to determine its effect on the disturbance behavior. To investigate the onset of three-dimensionality, the wake is forced with a small-amplitude three-dimensional disturbance and a larger amplitude two-dimensional disturbance. The two-dimensional forcing amplitude is varied in order to determine its influence on the three-dimensional flow field.

Fasel, Hermann F.↗