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 1,063 records · Page 59

Real-time minimal-bit-error probability decoding of convolutional codes

A recursive procedure is derived for decoding of rate R = 1/n binary convolutional codes which minimizes the probability of the individual decoding decisions for each information bit, subject to the constraint that the decoding delay be limited to Delta branches. This new decoding algorithm is similar to, but somewhat more complex than, the Viterbi decoding algorithm. A real-time, i.e., fixed decoding delay, version of the Viterbi algorithm is also developed and used for comparison to the new algorithm on simulated channels. It is shown that the new algorithm offers advantages over Viterbi decoding in soft-decision applications, such as in the inner coding system for concatenated coding.

Lee, L.-N.↗

A comparison of computational methods for three-dimensional, turbulent turbomachinery flow fields

A space-marching method and a time-marching method have been used to compute the three-dimensional turbulent flow in an end wall cascade of airfoils. Using an identical grid and turbulence model, the two codes were used to predict a variety of flow quantities. Predictions by the two methods are compared to each other and to experimental data. In general both methods predict measured quantities well, with a small edge in prediction accuracy going to the space-marching method. Secondary flow comparisons show the time-marching solution more accurately predicting the underturning of the flow in the outer portion of the end wall boundary layer while the space-marching method more accurately predicted the overturning of the flow very near the end wall. The prediction comparisons are discussed along with computational details and other attributes of the two methods.

Kirtley, K. R.↗

Comparison between computational and experimental data for a hypersonic laser propelled vehicle

A 3-D axisymmetric hypersonic engine inlet was investigated using PARC2D, an ideal gas Computational Fluid Dynamics code. The code was used to predict the results of tests conducted in the Rensselaer Polytechnic Institute Hypersonic Shock Tunnel which measured surface and pitot pressures, and shock positions (through Schlieren photography) at freestream Mach numbers of 10, 13, and 15. A strong viscous/shock interaction was observed in both the experiment and the CFD results, due to the model's parabolic compression ramp. Good agreement was found between the experimental results and the CFD solution both for surface pressures and shock positions. Agreement between pitot pressures was less reliable.

Messitt, D. G.↗

Comparison of level-1 and level-2 buckling and postbuckling solutions

The accuracy of the solutions obtained with the Level-1 computational modules of DISDECO (Delft Interactive Shell Design Code) are investigated by comparing them with solutions obtained with ANILISA, a Level-2 computational module. Whereas the Level-1 solutions are based on the use of trigonometric functions to represent the out-of-plane displacement W and involve Galerkin type approximations to reduce the governing partial differential equations to simple algebraic expressions, the Level-2 solutions employ a truncated Fourier expansion in the circumferential direction and solve the resulting ordinary differential equations forming a nonlinear eigenvalue problem by a Stodola like technique. This involves the numerical integration of the set of ordinary differential equations, which in turn makes it possible to satisfy the specified boundary conditions exactly. It is shown that in many applications the Level-1 predictions concerning the critical buckling load of the perfect structure and its imperfection sensitivity are quite accurate and reliable. For those cases where this is not true, it is shown how the pre and postprocessing facilities of DISDECO can be used efficiently to arrive at a reliable prediction. This involves always the use of the Level-2 computational module ANILISA with its rigorous prebuckling solution and its feature which makes it possible to satisfy the specified boundary conditions exactly.

Arbocz, J.↗

A comprehensive comparison of turbulence models in the far wake

In the present study, the far wake was examined numerically using an implicit, upwind, finite-volume, compressible Navier-Stokes code. The numerical grid started at 500 equivalent circular cylinder diameters in the wave, and extended to 4000 equivalent diameters. By concentrating only on the far wake, the numerical difficulties and fine mesh requirements near the wake-generating body were eliminated. At the time of this writing, results for the K-epsilon and K-omega turbulence models at low Mach number have been completed and show excellent agreement with previous incompressible results and far-wake similarity solutions. The code is presently being used to compare the performance of various other turbulence models, including Reynolds stress models and the new anisotropic two-equation turbulence models being developed at NASA Langley. By increasing our physical understanding of the deficiencies and limits of these models, it is hoped that improvements to the universality of the models can be made. Future plans include examination of two-dimensional momentumless wakes as well.

Cimbala, John M.↗

Numerical Simulation of Flow in a Whirling Annular Seal and Comparison with Experiments

The turbulent flow field in a simulated annular seal with a large clearance/radius ratio (0.015) and a whirling rotor was simulated using an advanced 3D CFD code SCISEAL. A circular whirl orbit with synchronous whirl was imposed on the rotor center. The flow field was rendered quasi-steady by making a transformation to a totaling frame. Standard k-epsilon model with wall functions was used to treat the turbulence. Experimentally measured values of flow parameters were used to specify the seal inlet and exit boundary conditions. The computed flow-field in terms of the velocity and pressure is compared with the experimental measurements inside the seal. The agreement between the numerical results and experimental data with correction is fair to good. The capability of current advanced CFD methodology to analyze this complex flow field is demonstrated. The methodology can also be extended to other whirl frequencies. Half- (or sub-) synchronous (fluid film unstable motion) and synchronous (rotor centrifugal force unbalance) whirls are the most unstable whirl modes in turbomachinery seals, and the flow code capability of simulating the flows in steady as well as whirling seals will prove to be extremely useful in the design, analyses, and performance predictions of annular as well as other types of seals.

Athavale, M. M.↗

Comparison of Organ Dosimetry for Astronaut Phantoms: Earth-Based vs. Microgravity-Based Anthropometry and Body Positioning

The purpose of this study is to use NASA radiation transport codes to compare astronaut organ dose equivalents resulting from solar particle events (SPE), geomagnetically trapped protons, and free-space galactic cosmic rays (GCR) using phantom models representing Earth-based and microgravity-based anthropometry and positioning. Methods: The Univer sity of Florida hybrid adult phantoms were scaled to represent male and female astronauts with 5th, 50th, and 95th percentile heights and weights as measured on Earth. Another set of scaled phantoms, incorporating microgravity-induced changes, such as spinal lengthening, leg volume loss, and the assumption of the neutral body position, was also created. A ray-tracer was created and used to generate body self-shielding distributions for dose points within a voxelized phantom under isotropic irradiation conditions, which closely approximates the free-space radiation environment. Simplified external shielding consisting of an aluminum spherical shell was used to consider the influence of a spacesuit or shielding of a hull. These distributions were combined with depth dose distributions generated from the NASA radiation transport codes BRYNTRN (SPE and trapped protons) and HZETRN (GCR) to yield dose equivalent. Many points were sampled per organ. Results: The organ dos e equivalent rates were on the order of 1.5-2.5 mSv per day for GCR (1977 solar minimum) and 0.4-0.8 mSv per day for trapped proton irradiation with shielding of 2 g cm-2 aluminum equivalent. The organ dose equivalents for SPE irradiation varied considerably, with the skin and eye lens having the highest organ dose equivalents and deep-seated organs, such as the bladder, liver, and stomach having the lowest. Conclus ions: The greatest differences between the Earth-based and microgravity-based phantoms are observed for smaller ray thicknesses, since the most drastic changes involved limb repositioning and not overall phantom size. Improved self-shielding models reduce the overall uncertainty in organ dosimetry for mission-risk projections and assessments for astronauts

VanBaalen, Mary↗

Computation of axisymmetric and ionized flows using particle and continuum methods

Comparisons between particle and continuum simulations of hypersonic near-continuum flows are presented. The particle approach employs the direct simulation Monte Carlo method (DSMC), and the continuum approach solves the Euler equations. Both simulations have thermochemistry models for air implemented including ionization. A new axisymmetric DSMC code which is efficiently vectorized is developed for this study. In this DSMC code, particular attention is paid to matching the relaxation rates employed in the continuum approach. This investigation represents a continuation of a previous study which considered thermochemical relaxation in one-dimensional shock waves of nitrogen. Comparison of the particle and continuum methods is first made for an axisymmetric blunt-body flow of air at 7 km/s. Very good agreement is obtained for the two solutions. The two techniques also compare well for a one-dimensional shock wave in air at 10 km/s. In both applications, the results are found to be sensitive to various aspects of the chemistry models employed.

Boyd, Iain D.↗

Application of FUN3D and CFL3D to the Third Workshop on CFD Uncertainty Analysis

Two Reynolds-averaged Navier-Stokes computer codes - one unstructured and one structured - are applied to two workshop cases (for the 3rd Workshop on CFD Uncertainty Analysis, held at Instituto Superior Tecnico, Lisbon, in October 2008) for the purpose of uncertainty analysis. The Spalart-Allmaras turbulence model is employed. The first case uses the method of manufactured solution and is intended as a verification case. In other words, the CFD solution is expected to approach the exact solution as the grid is refined. The second case is a validation case (comparison against experiment), for which modeling errors inherent in the turbulence model and errors/uncertainty in the experiment may prevent close agreement. The results from the two computer codes are also compared. This exercise verifies that the codes are consistent both with the exact manufactured solution and with each other. In terms of order property, both codes behave as expected for the manufactured solution. For the backward facing step, CFD uncertainty on the finest grid is computed and is generally very low for both codes (whose results are nearly identical). Agreement with experiment is good at some locations for particular variables, but there are also many areas where the CFD and experimental uncertainties do not overlap.

Rumsey, C. L.↗

GOME Calibration and Validation Using Backscatter UV Techniques

GOME radiance, irradiance, and ozone products were validated by NASA Goddard Space Flight Center through three tasks which included, pre-launch calibration comparisons with SBUV and TOMS radiometric standards, validation of GOME Level-1 irradiance and radiance and Level 2 total ozone data products using SBUV/2 and TOMS algorithms and data, and studies of GOME data using the Goddard radiative transfer code. The prelaunch calibration using the NASA large aperture integrating sphere was checked against that provided by TPD. Agreement in the calibration constants, derived in air, between the Goddard and TPD system were better than 3%. Validation of Level-1 irradiance data included comparison of GOME and SSBUV and the UARS solar irradiances measurements. Large wavelength dependent differences, as high as 10%, were noted between GOME and the US instruments. This discrepancy has now been attributed to radiometric sensitivity changes experienced by GOME when operating in a vacuum. GOME Earth radiance data were then compared to the NOAA-14 SBUV/2 radiances. These results show that between 340 and 400 nm the differences in GOME and SBUV/2 data are less than 5% with some wavelength dependence. At wavelengths shorter than 300 nm, differences are of the order of 10% or more where the GOME radiances are larger. To test GOME DOAS retrieved total ozone values, these values were compared with ozone amounts retrieved using GOME radiances in the TOMS version-7 algorithm. The differences showed a solar zenith angle dependence ranging from 0 to 10% where the TOMS algorithm values were higher. GOME radiances below 300 nm were further validated by selecting radiances at wavelengths normally used by SBUV and processing them through the SBUV ozone profile algorithm and then compared to climatological values. The GOME ozone profiles ranged from 10-30% lower over altitude compared to climatological values. This is consistent with the offsets detected in the SBUV/2 radiance comparisons at wavelengths shorter than 300 nm.

Hilsenrath, E.↗

A Comparison of Telemetry Signals in the Presence and Absence of a Subcarrier

The performance of deep-space telemetry signals that employ a residual carrier modulation technique is compared in the presence and absence of a subcarrier. When the subcarrier is present, the performance for the resulting pulse-coded modulation/phase-shift keyed/phase-modulated (PCM/PSK/PM) scheme is evaluated for both sine-wave and square-wave subcarriers and Non-Return-to Zero (NRZ) data. When the subcarrier is absent, the performance for the resulting PCM/PM technique is evaluated for both the NRZ and the bi-phase data format. The comparison is based on telemetry performance as well as bandwidth efficiency. The first criterion is characterized in terms of the Symbol Error Rate (SER) as a function of symbol SNR, loop bandwidth-to-data rate ratio, and modulation index. The bandwidth efficiency is characterized by the occupancy factor. The results of both the analysis and measurements show that when the Interference-to-Carrier Ratio (ICR) is less than -20 dB, the performance degradation in the absence of a subcarrier is negligible. Various combinations of loop bandwidth-to-data rate ratios and modulation indices that achieve this performance are derived and listed. Bandwidth occupancy comparison indicates that PCM/PM/NRZ is the most efficient in this regard. Therefore, by eliminating the subcarrier and using the PCM/PM/NRZ scheme, many advantages can be realized without any sacrifice in performance.

Shihabi, Mazen M.↗

Comparison of Models for Ball Bearing Dynamic Capacity and Life

Generalized formulations for dynamic capacity and life of ball bearings, based on the models introduced by Lundberg and Palmgren and Zaretsky, have been developed and implemented in the bearing dynamics computer code, ADORE. Unlike the original Lundberg-Palmgren dynamic capacity equation, where the elastic properties are part of the life constant, the generalized formulations permit variation of elastic properties of the interacting materials. The newly updated Lundberg-Palmgren model allows prediction of life as a function of elastic properties. For elastic properties similar to those of AISI 52100 bearing steel, both the original and updated Lundberg-Palmgren models provide identical results. A comparison between the Lundberg-Palmgren and the Zaretsky models shows that at relatively light loads the Zaretsky model predicts a much higher life than the Lundberg-Palmgren model. As the load increases, the Zaretsky model provides a much faster drop off in life. This is because the Zaretsky model is much more sensitive to load than the Lundberg-Palmgren model. The generalized implementation where all model parameters can be varied provides an effective tool for future model validation and enhancement in bearing life prediction capabilities.

Rolling-element Fatigue↗

Comparison of Fluka-2006 Monte Carlo Simulation and Flight Data for the ATIC Detector

We have performed a detailed Monte Carlo (MC) simulation for the Advanced Thin Ionization Calorimeter (ATIC) detector using the MC code FLUKA-2006 which is capable of simulating particles up to 10 PeV. The ATIC detector has completed two successful balloon flights from McMurdo, Antarctica lasting a total of more than 35 days. ATIC is designed as a multiple, long duration balloon flight, investigation of the cosmic ray spectra from below 50 GeV to near 100 TeV total energy; using a fully active Bismuth Germanate(BGO) calorimeter. It is equipped with a large mosaic of.silicon detector pixels capable of charge identification, and, for particle tracking, three projective layers of x-y scintillator hodoscopes, located above, in the middle and below a 0.75 nuclear interaction length graphite target. Our simulations are part of an analysis package of both nuclear (A) and energy dependences for different nuclei interacting in the ATIC detector. The MC simulates the response of different components of the detector such as the Si-matrix, the scintillator hodoscopes and the BGO calorimeter to various nuclei. We present comparisons of the FLUKA-2006 MC calculations with GEANT calculations and with the ATIC CERN data and ATIC flight data.

Gunasingha, R.M.↗

Computation of axisymmetric and ionized hypersonic flows using particle and continuum methods

Comparisons between particle and continuum simulations of hypersonic near-continuum flows are presented. The particle approach employs the direct simulation Monte Carlo (DSMC) method, and the continuum approach solves the appropriate equations of fluid flow. Both simulations have thermochemistry models for air implemented including ionization. A new axisymmetric DSMC code that is efficiently vectorized is developed for this study. In this DSMC code, particular attention is paid to matching the relaxation rates employed in the continuum approach. This investigation represents a continuum of a previous study that considered thermochemical relaxation in one-dimensional shock waves of nitrogen. Comparison of the particle and continuum methods is first made for an axisymmetric blunt-body flow of air at 7 km/s. Very good agreement is obtained for the two solutions. The two techniques also compare well for a one-dimensional shock wave in air at 10 km/s. In both applications, the results are found to be sensitive to various aspects of the chemistry model employed.

Boyd, Iain D.↗

Modeling Film-Coolant Flow Characteristics at the Exit of Shower-Head Holes

The coolant flow characteristics at the hole exits of a film-cooled blade are derived from an earlier analysis where the hole pipes and coolant plenum were also discretized. The blade chosen is the VKI rotor with three staggered rows of shower-head holes. The present analysis applies these flow characteristics at the shower-head hole exits. A multi-block three-dimensional Navier-Stokes code with Wilcox's k-omega model is used to compute the heat transfer coefficient on the film-cooled turbine blade. A reasonably good comparison with the experimental data as well as with the more complete earlier analysis where the hole pipes and coolant plenum were also gridded is obtained. If the 1/7th power law is assumed for the coolant flow characteristics at the hole exits, considerable differences in the heat transfer coefficient on the blade surface, specially in the leading-edge region, are observed even though the span-averaged values of h (heat transfer coefficient based on T(sub o)-T(sub w)) match well with the experimental data. This calls for span-resolved experimental data near film-cooling holes on a blade for better validation of the code.

Garg, Vijay K.↗

Outer satellite atmospheres: Their nature and planetary interactions

Significant progress is reported in early modeling analysis of observed sodium cloud images with our new model which includes the oscillating Io plasma torus ionization sink. Both the general w-D morphology of the region B cloud as well as the large spatial gradient seen between the region A and B clouds are found to be consistent with an isotropic flux of sodium atoms from Io. Model analysis of the spatially extended high velocity directional features provided substantial evidence for a magnetospheric wind driven gas escape mechanism from Io. In our efforts to define the source(s) of hydrogen atoms in the Saturn system, major steps were taken in order to understand the role of Titan. We have completed the comparison of the Voyager UVS data with previous Titan model results, as well as the update of the old model computer code to handle the spatially varying ionization sink for H atoms.

Smyth, W. H.↗

Computation of supersonic viscous flows around pointed bodies at large incidence

A recently reported parabolized Navier-Stokes method has been extended to compute turbulent supersonic flows around cones and an ogive-cylinder body at large incidence. The algebraic eddy-viscosity turbulence model contained in the code was modified to properly account for the large regions of cross-flow separation that occur in these flows. Extensive comparisons between computed results and experimentally measured flow fields are presented. The results show good agreement for viscous-layer profiles and details of the external leeward-side vortex structure at angles of attack up to three times the cone half angles. Details of the modified turbulence model are presented and discussed.

Degani, D.↗

Inner-shell transitions in Fe XIX-XXII in the X-ray spectra of solar flares and Tokamaks

Calculated spectra of the ions Fe XIX-XXII for various densities and temperatures are presented, thereby extending the work begun by Doschek, Feldman, and Cowan (1981). The calculations are based on a code (the Cowan code) that computes both the level structure of an ion and intensity factors for the 1s-2p satellite lines. A comparison is made between the calculated spectra and those observed in solar flares by the P78-1 and SMM instruments. The observed intensities of Fe XX lines, which are the most sensitive to density, are found to agree well with those calculated in the low-density limit. The agreement for lines arising from other ions is also very good. It is also seen that the predicted density variations in Fe XX are confirmed by the higher density Princeton Large Torus plasmas. Thus a possible useful density diagnostic is indicated for tokamak and high-density astrophysical plasmas, perhaps including some solar flares.

Phillips, K. J. H.↗