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,351 records · Page 75

Application of a full potential method for analysis of complex aircraft geometries

A supersonic potential flow solver was developed to analyze the flow over complex realistic aircraft geometries. Enhancements to the method were made to accommodate regions of subsonic flow, the effect of trailing wakes on other aircraft components, and the modeling/gridding of complete configurations. Validation of the method was demonstrated by comparisons with experimental aerodynamic force and surface pressure measurements. The predicted results are in very good agreement with the experimental data. The bibliography contains additional information on the use of the potential flow code to predict the aerodynamics of high-speed wing/body configurations, waverider concepts, TAV, and the Space Shuttle orbiter package.

Jones, Kenneth M.↗

LU-SGS implicit algorithm for three-dimensional incompressible Navier-Stokes equations with source term

A numerical method is developed for solving the incompressible Navier-Stokes equations using the concept of pseudocompressibility. A lower-upper symmetric-Gauss-Seidel implicit scheme is developed for three-dimensional incompressible viscous flow computations. The present algorithm offers additional advantages when solving the flow equations with source terms. Complete vectorizability of the algorithm on oblique planes of sweep in three-dimensions is accomplished in a new flow solver, INS3D-LU code. Spatial differencing is a second-order accurate semi-discrete finite-volume method augmented by a third-order accurate numerical dissipation model which is based on spectral-radii. Comparison of numerical solutions for a curved duct with experimental data shows good agreement. The method is applied to calculate the inducer flow of the Space Shuttle Main Engine turbopump.

Yoon, Seokkwan↗

Modeling of jet-induced geyser formation in a reduced gravity environment

Flow patterns predicted by a computational model of jet-induced geyser formation in a reduced gravity environment are presented and comparison is made to patterns predicted by experimentally based correlations. The configuration studied is an idealization of a forthcoming flight experiment to examine cryogenic propellant management issues. A transitional version of the ECLIPSE code used as a computational tool for the analyses is described. It is shown that computationally predicted flow patterns are in qualitative agreement with the correlation-based predictions, and some details of the predicted flow fields are given.

Wendl, M. C.↗

Burnett solutions along the stagnation line of a cooled cylinder in low-density hypersonic flows

The Burnett equations in the cylindrical coordinate system have been applied to the hypersonic flow over a cylinder in low density environments. Through a local similarity transformation, the Burnett equations are degenerated into nonlinear third-order ordinary differential equations which are only valid along the stagnation line. These equations are solved by the multi-staged Runge-Kutta integration technique. Steady state solutions are asymptotically reached by the time-marching procedure. A newly-developed computer code is used to solve both the Navier-Stokes and the Burnett equations within the shock layer. The computed heat transfer coefficients are consistent with wind-tunnel and STS flight data. Comparisons with the DSMC solutions and the Boltzmann solutions show that the Burnett equations are more appropriate as the governing equations than the Navier-Stokes equations in the transitional flow regime.

Liaw, G. S.↗

On the optimality of a universal noiseless coder

Rice developed a universal noiseless coding structure that provides efficient performance over an extremely broad range of source entropy. This is accomplished by adaptively selecting the best of several easily implemented variable length coding algorithms. Variations of such noiseless coders have been used in many NASA applications. Custom VLSI coder and decoder modules capable of processing over 50 million samples per second have been fabricated and tested. In this study, the first of the code options used in this module development is shown to be equivalent to a class of Huffman code under the Humblet condition, for source symbol sets having a Laplacian distribution. Except for the default option, other options are shown to be equivalent to the Huffman codes of a modified Laplacian symbol set, at specified symbol entropy values. Simulation results are obtained on actual aerial imagery over a wide entropy range, and they confirm the optimality of the scheme. Comparison with other known techniques are performed on several widely used images and the results further validate the coder's optimality.

Yeh, Pen-Shu↗

Computer Simulation Results and Analysis for a Root-Raised Cosine Filter Design using Canonical Signed Digits

Filters implemented using Canonical Signed-Digit (CSD) number representation for the coefficients have been studied for many years as an efficient way to increase the speed and reduce the hardware complexity. A simple and very effective method for designing a CSD filter with a specific set of filter parameters has been investigated and its simulation results are presented here. In order to optimize filter coefficients into the corresponding CSD numbers, the Minimum Mean Square Error (MMSE) criterion is used. Furthermore, an attempt is made to improve frequency response of the CSD filter by allocating an extra non-zero digit for normalized coefficients exceeding one-half. Due to limited filter aperture width in examples presented here, the frequency response of CSD filter as the number of non-zero digits in a CSD code increases is not affected much. A root-raised cosine filter model is employed as a base line for this design approach. Finally, the CSD filter simulation results and a hardware complexity comparison with a conventional filter are also shown along with the eye diagrams an Bit-Error-Rate (BER) performance curves.

Kim, Heechul↗

Progressive Failure And Life Prediction of Ceramic and Textile Composites

An engineering approach to predict the fatigue life and progressive failure of multilayered composite and textile laminates is presented. Analytical models which account for matrix cracking, statistical fiber failures and nonlinear stress-strain behavior have been developed for both composites and textiles. The analysis method is based on a combined micromechanics, fracture mechanics and failure statistics analysis. Experimentally derived empirical coefficients are used to account for the interface of fiber and matrix, fiber strength, and fiber-matrix stiffness reductions. Similar approaches were applied to textiles using Repeating Unit Cells. In composite fatigue analysis, Walker's equation is applied for matrix fatigue cracking and Heywood's formulation is used for fiber strength fatigue degradation. The analysis has been compared with experiment with good agreement. Comparisons were made with Graphite-Epoxy, C/SiC and Nicalon/CAS composite materials. For textile materials, comparisons were made with triaxial braided and plain weave materials under biaxial or uniaxial tension. Fatigue predictions were compared with test data obtained from plain weave C/SiC materials tested at AS&M. Computer codes were developed to perform the analysis. Composite Progressive Failure Analysis for Laminates is contained in the code CPFail. Micromechanics Analysis for Textile Composites is contained in the code MicroTex. Both codes were adapted to run as subroutines for the finite element code ABAQUS and CPFail-ABAQUS and MicroTex-ABAQUS. Graphic user interface (GUI) was developed to connect CPFail and MicroTex with ABAQUS.

Xue, David Y.↗

Integration of Libration Point Orbit Dynamics into a Universal 3-D Autonomous Formation Flying Algorithm

The autonomous formation flying control algorithm developed by the Goddard Space Flight Center (GSFC) for the New Millennium Program (NMP) Earth Observing-1 (EO-1) mission is investigated for applicability to libration point orbit formations. In the EO-1 formation-flying algorithm, control is accomplished via linearization about a reference transfer orbit with a state transition matrix (STM) computed from state inputs. The effect of libration point orbit dynamics on this algorithm architecture is explored via computation of STMs using the flight proven code, a monodromy matrix developed from a N-body model of a libration orbit, and a standard STM developed from the gravitational and coriolis effects as measured at the libration point. A comparison of formation flying Delta-Vs calculated from these methods is made to a standard linear quadratic regulator (LQR) method. The universal 3-D approach is optimal in the sense that it can be accommodated as an open-loop or closed-loop control using only state information.

Folta, David↗

Numerical Simulation of Protoplanetary Vortices

The fluid dynamics within a protoplanetary disk has been attracting the attention of many researchers for a few decades. Previous works include, to list only a few among many others, the well-known prescription of Shakura & Sunyaev, the convective and instability study of Stone & Balbus and Hawley et al., the Rossby wave approach of Lovelace et al., as well as a recent work by Klahr & Bodenheimer, which attempted to identify turbulent flow within the disk. The disk is commonly understood to be a thin gas disk rotating around a central star with differential rotation (the Keplerian velocity), and the central quest remains as how the flow behavior deviates (albeit by a small amount) from a strong balance established between gravitational and centrifugal forces, transfers mass and momentum inward, and eventually forms planetesimals and planets. In earlier works we have briefly described the possible physical processes involved in the disk; we have proposed the existence of long-lasting, coherent vortices as an efficient agent for mass and momentum transport. In particular, Barranco et al. provided a general mathematical framework that is suitable for the asymptotic regime of the disk; Barranco & Marcus (2000) addressed a proposed vortex-dust interaction mechanism which might lead to planetesimal formation; and Lin et al. (2002), as inspired by general geophysical vortex dynamics, proposed basic mechanisms by which vortices can transport mass and angular momentum. The current work follows up on our previous effort. We shall focus on the detailed numerical implementation of our problem. We have developed a parallel, pseudo-spectral code to simulate the full three-dimensional vortex dynamics in a stably-stratified, differentially rotating frame, which represents the environment of the disk. Our simulation is validated with full diagnostics and comparisons, and we present our results on a family of three-dimensional, coherent equilibrium vortices.

Lin, H.↗

Validity of the Aluminum Equivalent Approximation in Space Radiation Shielding

The origin of the aluminum equivalent shield approximation in space radiation analysis can be traced back to its roots in the early years of the NASA space programs (Mercury, Gemini and Apollo) wherein the primary radiobiological concern was the intense sources of ionizing radiation causing short term effects which was thought to jeopardize the safety of the crew and hence the mission. Herein, it is shown that the aluminum equivalent shield approximation, although reasonably well suited for that time period and to the application for which it was developed, is of questionable usefulness to the radiobiological concerns of routine space operations of the 21 st century which will include long stays onboard the International Space Station (ISS) and perhaps the moon. This is especially true for a risk based protection system, as appears imminent for deep space exploration where the long-term effects of Galactic Cosmic Ray (GCR) exposure is of primary concern. The present analysis demonstrates that sufficiently large errors in the interior particle environment of a spacecraft result from the use of the aluminum equivalent approximation, and such approximations should be avoided in future astronaut risk estimates. In this study, the aluminum equivalent approximation is evaluated as a means for estimating the particle environment within a spacecraft structure induced by the GCR radiation field. For comparison, the two extremes of the GCR environment, the 1977 solar minimum and the 2001 solar maximum, are considered. These environments are coupled to the Langley Research Center (LaRC) deterministic ionized particle transport code High charge (Z) and Energy TRaNsport (HZETRN), which propagates the GCR spectra for elements with charges (Z) in the range I <= Z <= 28 (H -- Ni) and secondary neutrons through selected target materials. The coupling of the GCR extremes to HZETRN allows for the examination of the induced environment within the interior' of an idealized spacecraft as approximated by a spherical shell shield, and the effects of the aluminum equivalent approximation for a good polymeric shield material such as genetic polyethylene (PE). The shield thickness is represented by a 25 g/cm spherical shell. Although one could imagine the progression to greater thickness, the current range will be sufficient to evaluate the qualitative usefulness of the aluminum equivalent approximation. Upon establishing the inaccuracies of the aluminum equivalent approximation through numerical simulations of the GCR radiation field attenuation for PE and aluminum equivalent PE spherical shells, we Anther present results for a limited set of commercially available, hydrogen rich, multifunctional polymeric constituents to assess the effect of the aluminum equivalent approximation on their radiation attenuation response as compared to the generic PE.

Badavi, Francis F.↗

A New Analysis Tool Assessment for Rotordynamic Modeling of Gas Foil Bearings

Gas foil bearings offer several advantages over traditional bearing types that make them attractive for use in high-speed turbomachinery. They can operate at very high temperatures, require no lubrication supply (oil pumps, seals, etc.), exhibit very long life with no maintenance, and once operating airborne, have very low power loss. The use of gas foil bearings in high-speed turbomachinery has been accelerating in recent years, although the pace has been slow. One of the contributing factors to the slow growth has been a lack of analysis tools, benchmarked to measurements, to predict gas foil bearing behavior in rotating machinery. To address this shortcoming, NASA Glenn Research Center (GRC) has supported the development of analytical tools to predict gas foil bearing performance. One of the codes has the capability to predict rotordynamic coefficients, power loss, film thickness, structural deformation, and more. The current paper presents an assessment of the predictive capability of the code, named XLGFBTH (Texas A&M University). A test rig at GRC is used as a simulated case study to compare rotordynamic analysis using output from the code to actual rotor response as measured in the test rig. The test rig rotor is supported on two gas foil journal bearings manufactured at GRC, with all pertinent geometry disclosed. The resulting comparison shows that the rotordynamic coefficients calculated using XLGFBTH represent the dynamics of the system reasonably well, especially as they pertain to predicting critical speeds.

Howard, Samuel A.↗

Numerical Modeling of an Integrated Vehicle Fluids System Loop for Pressurizing a Cryogenic Tank

This paper presents a numerical model of the pressurization loop of the Integrated Vehicle Fluids (IVF) system using the Generalized Fluid System Simulation Program (GFSSP). The IVF propulsion system, being developed by United Launch Alliance to reduce system weight and enhance reliability, uses boiloff propellants to drive thrusters for the reaction control system as well as to run internal combustion engines to develop power and drive compressors to pressurize propellant tanks. NASA Marshall Space Flight Center (MSFC) conducted tests to verify the functioning of the IVF system using a flight-like tank. GFSSP, a finite volume based flow network analysis software developed at MSFC, has been used to support the test program. This paper presents the simulation of three different test series, comparison of numerical prediction and test data and a novel method of presenting data in a dimensionless form. The paper also presents a methodology of implementing a compressor map in a system level code.

LeClair, A. C.↗

A comparison of computed and measured aerodynamic characteristics of a proposed aeroassist flight experiment configuration

The use of experimental or computed data to evaluate the performance of aeroassist orbital transfer vehicles (AOTVs) is discussed. Aerodynamic and surface pressure data are derived from computed flowfield solutions (the HALIS inviscid flowfield code), the Newtonian theory, and ground-based, wind tunnel studies (hypersonic He tunnel, 31-inch Mach 10 tunnel, and a hypersonic CF4 tunnel). The wind tunnel and HALIS models were tested at angles-of-attack that ranged from -10 to 10 deg. The effects of the ellipticity of the nose, the radius of the circular arc, and the angle through which the arc passes on the aerodynamic characteristics of the vehicles are examined. It is observed that the HALIS and CF4 tunnel generated surface pressure on the AOTVs produced the most useful aerodynamic data. Good correlation is obtained for the experimental and computational data; however, the Newtonian results do not correspond to the tunnel/HALIS data.

Weilmuenster, K. J.↗

Comparisons of Cloud In-Situ Microphysical Properties of Deep Convective Clouds to Appendix D/P using Data from the HAIC-HIWC and HIWC-RADAR I Flight Campaigns

In-situ cloud data from three international flight campaigns are compared to the Federal Aviation Administration Title 14 Code of Federal Regulations Part 33 Appendix D mixed-phase/glaciated environmental envelope, and the corresponding identical European Aviation Safety Agency CS-25 Appendix P envelope. The appendices consist of a temperature-altitude envelope, a 99th percentile total water content envelope at the 17.4 Nm distance scale, a distance factor for estimation at other distance scales, ice crystal median mass diameter, and recommended liquid water content levels in mixed-phase icing conditions. The data were collected during 54 flights out of one subtropical and two tropical locations, with 472 runs from about 17,000’ to 39,000’ in approximately 115 clouds. The campaigns provide about 29,600 Nm of in situ data in deep convection over four targeted temperature intervals: -10, -30, -40, and -50, all ± 5 C. The dataset is a modern and unique documentation of the ice crystal icing environment, and results described in this article will contribute to regulatory and industry assessment of Appendices D and P.

Ice Crystal Icing↗

Comparison of Measured Dark Current Distributions with Calculated Damage Energy Distributions in HgCdTe

This paper presents a combined Monte Carlo and analytic approach to the calculation of the pixel-to-pixel distribution of proton-induced damage in a HgCdTe sensor array and compares the results to measured dark current distributions after damage by 63 MeV protons. The moments of the Coulombic, nuclear elastic and nuclear inelastic damage distributions were extracted from Monte Carlo simulations and combined to form a damage distribution using the analytic techniques first described in [1]. The calculations show that the high energy recoils from the nuclear inelastic reactions (calculated using the Monte Carlo code MCNPX [2]) produce a pronounced skewing of the damage energy distribution. While the nuclear elastic component (also calculated using the MCNPX) contributes only a small fraction of the total nonionizing damage energy, its inclusion in the shape of the damage across the array is significant. The Coulombic contribution was calculated using MRED [3-5], a Geant4 [4,6] application. The comparison with the dark current distribution strongly suggests that mechanisms which are not linearly correlated with nonionizing damage produced according to collision kinematics are responsible for the observed dark current increases. This has important implications for the process of predicting the on-orbit dark current response of the HgCdTe sensor array.

Marshall, C. J.↗

Second-order-accurate spatial differencing for the transonic small-disturbance equation

Current methods for calculating transonic flows with the small-disturbance potential equation are typically only first-order accurate in the supersonic regions of the flow. However, calculations using the full-potential equation show significant improvements in accuracy when second-order methods are used instead of first-order methods. In this paper, algorithms for the small-disturbance equations are presented, for both steady and unsteady flows, with spatial differencing that is second-order accurate in both the subsonic and supersonic regions of the flow. These algorithms are stable, simple extensions of implicit monotone algorithms; the former are only first-order accurate spatially in the supersonic regions. Several calculations in one and two dimensions have been made, and the first- and second-order calculations are compared. In one dimension, the calculations include the four types of shock-wave motion. In two dimensions, the comparisons include steady flow over a Korn airfoil and unsteady flow over a pitching airfoil. All the comparisons in two dimensions show that the new methods are just as stable numerically as the old methods, but improve the accuracy of the solutions. The algorithm improvements can be implemented in present computer codes by making minor coding modifications.

Goorjian, P. M.↗

Characterizing Computers And Predicting Computing Times

Improved method for evaluation and comparison of computers running same or different FORTRAN programs devised. Enables one to predict time necessary to run given "benchmark" or other standard program on given computer, in scalar mode and without optimization of codes generated by compiler. Such "benchmark" running times are principal measures used to characterize performances of computers; of interest to designers, manufacturers, programmers, and users.

Saavedra-Barrera, Rafael H.↗

An assumed pdf approach for the calculation of supersonic mixing layers

In an effort to predict the effect that turbulent mixing has on the extent of combustion, a one-equation turbulence model is added to an existing Navier-Stokes solver with finite-rate chemistry. To average the chemical-source terms appearing in the species-continuity equations, an assumed pdf approach is also used. This code was used to analyze the mixing and combustion caused by the mixing layer formed by supersonic coaxial H2-air streams. The chemistry model employed allows for the formation of H2O2 and HO2. Comparisons are made with recent measurements using laser Raman diagnostics. Comparisons include temperature and its rms, and concentrations of H2, O2, N2, H2O, and OH. In general, good agreement with experiment was noted.

Baurle, R. A.↗