SEARCH · Engineering Papers

Results for “code comparisons”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 685 records · Page 38

All sources NASA NTRS DOE OSTI

NASA NTRS1993

CFD Modeling of Turbulent Duct Flows for Coolant Channel Analysis

The design of modern liquid rocket engines requires the analysis of chamber coolant channels to maximize the heat transfer while minimizing the coolant flow. Coolant channels often do not remain at a constant cross section or at uniform curvature. New designs require higher aspect ratio coolant channels than previously used. To broaden the analysis capability and to complement standard analysis tools an investigation on the accuracy of CFD predictions for coolant channel flow has been initiated. Validation of CFD capabilities for coolant channel analysis will enhance the capabilities for optimizing design parameters without resorting to extensive experimental testing. The eventual goal is to use CFD to determine the flow fields of unique coolant channel designs and therefore determine critical heat transfer coefficients. In this presentation the accuracy of a particular CFD code is evaluated for turbulent flows. The first part of the presentation is a comparison of numerical results to existing cold flow data for square curved ducts (NASA CR-3367, 'Measurements of Laminar and Turbulent Flow in a Curved Duct with Thin Inlet Boundary Layers'). The results of this comparison show good agreement with the relatively coarse experimental data. The second part of the presentation compares two cases of higher aspect ratio channels (AR=2.5,10) to show changes in axial and secondary flow strength. These cases match experimental work presently in progress and will be used for future validation. The comparison shows increased secondary flow strength of the higher aspect ratio case due to the change in radius of curvature. The presentation includes a test case with a heated wall to demonstrate the program's capability. The presentation concludes with an outline of the procedure used to validate the CFD code for future design analysis.

Ungewitter, Ronald J.↗

NASA NTRS2000

Experimental Results from the Thermal Energy Storage-2 (TES-2) Flight Experiment

Thermal Energy Storage-2 (TES-2) is a flight experiment that flew on the Space Shuttle Endeavour (STS-72), in January 1996. TES-2 originally flew with TES-1 as part of the OAST-2 Hitchhiker payload on the Space Shuttle Columbia (STS-62) in early 1994. The two experiments, TES-1 and TES-2 were identical except for the fluoride salts to be characterized. TES-1 provided data on lithium fluoride (LiF), TES-2 provided data on a fluoride eutectic (LiF/CaF2). Each experiment was a complex autonomous payload in a Get-Away-Special payload canister. TES-1 operated flawlessly for 22 hr. Results were reported in a paper entitled, Effect of Microgravity on Materials Undergoing Melting and Freezing-The TES Experiment, by David Namkoong et al. A software failure in TES-2 caused its shutdown after 4 sec of operation. TES-1 and 2 were the first experiments in a four experiment suite designed to provide data for understanding the long duration microgravity behavior of thermal energy storage salts that undergo repeated melting and freezing. Such data have never been obtained before and have direct application for the development of space-based solar dynamic (SD) power systems. These power systems will store energy in a thermal energy salt such as lithium fluoride or a eutectic of lithium fluoride/calcium difluoride. The stored energy is extracted during the shade portion of the orbit. This enables the solar dynamic power system to provide constant electrical power over the entire orbit. Analytical computer codes were developed for predicting performance of a space-based solar dynamic power system. Experimental verification of the analytical predictions were needed prior to using the analytical results for future space power design applications. The four TES flight experiments were to be used to obtain the needed experimental data. This paper will address the flight results from the first and second experiments, TES-1 and 2, in comparison to the predicted results from the Thermal Energy Storage Simulation (TESSIM) analytical computer code. An analysis of the TES-2 data was conducted by Cleveland State University Professor, Mounir Ibrahim. TESSIM validation was based on two types of results; temperature history of various points on the containment vessel and TES material distribution within the vessel upon return from flight. The TESSIM prediction showed close comparison with the flight data. Distribution of the TES material within the vessel was obtained by a tomography imaging process. The frozen TES material was concentrated toward the colder end of the canister. The TESSIM prediction indicated a similar pattern. With agreement between TESSIM and the flight data, a computerized representation was produced to show the movement and behavior of the void during the entire melting and freezing cycles.

Tolbert, Carol↗

NASA NTRS1975

NASTRAN nonlinear vibration analysis of beam and frame structures

A capability for the nonlinear vibration analysis of beam and frame structures suitable for use with NASTRAN level 15.5 is described. The nonlinearity considered is due to the presence of axial loads induced by longitudinal end restraints and lateral displacements that are large compared to the beam height. A brief discussion is included of the mathematical analysis and the geometrical stiffness matrix for a prismatic beam (BAR) element. Also included are a brief discussion of the equivalent linearization iterative process used to determine the nonlinear frequency, the required modifications to subroutines DBAR and XMPLBD of the NASTRAN code, and the appropriate vibration capability, four example problems are presented. Comparisons with existing experimental and analytical results show that excellent accuracy is achieved with NASTRAN in all cases.

Mei, C.↗

NASA NTRS1975

A modified Henyey method for computing radiative transfer hydrodynamics

The implicit hydrodynamic code of Kutter and Sparks (1972), which is limited to optically thick regions and employs the diffusion approximation for radiative transfer, is modified to include radiative transfer effects in the optically thin regions of a model star. A modified Henyey method is used to include the solution of the radiative transfer equation in this implicit code, and the convergence properties of this method are proven. A comparison is made between two hydrodynamic models of a classical Cepheid with a 12-day period, one of which was computed with the diffusion approximation and the other with the modified Henyey method. It is found that the two models produce nearly identical light and velocity curves, but differ in the fact that the former never has temperature inversions in the atmosphere while the latter does when sufficiently strong shocks are present.

Karp, A. H.↗

NASA NTRS1978

Nonlinear behavior of circular plates with work hardening

Tests were performed on two simply supported plates of aluminum alloy 2024-0, under a central concentrated load, with peak deflections up to 2.6 times the thickness. The load was provided by a small-diameter hard-steel rod. The plates had diameter-to-thickness ratios (D/h) of 20 and 41. Measurements were made of load, deflections and strains; membrane and bending strains were calculated from the test data. The test data are presented in comparison with theoretical predictions generated by the finite-element-computer code PLANS, which includes material and geometric nonlinearities. The theoretical prediction was excellent for deflections, and generally good for strains, when the central force was represented by a line load around the loading rod's contact circle.

Winter, R.↗

NASA NTRS1983

Thermomechanical loading of multilayered cylindrical geometries in thermal cycling from 300 to 1300 K

The principal of multiple material layering is well known as an effective method of reducing heat transfer; however, thermal gradients can impose significant loads and lead to delamination and subsequent component failure. An analysis is developed and experimental data are discussed for the thermomechanical effects of multilayered materials on a heat sink substrate of cylindrical geometry subject to thermal cycling. The geometry is heated in cross-flow by a high-velocity flame and cooled in crossflow by ambient-temperature air from a critical flow orifice. Each layer of material possesses a threshold beyond which small changes in temperature or mechanical loading greatly influence the life in thermal cycling of the layered materials. Comparisons are made between the thermomechanical loads predicted by various numerical codes for the linear case and by a simplified analytic model.

Hendricks, R. C.↗

NASA NTRS1984

Potential flow analysis of glaze ice accretions on an airfoil

The results of an analytical/experimental study of the flow fields about an airfoil with leading edge glaze ice accretion shapes are presented. Tests were conducted in the Icing Research Tunnel to measure surface pressure distributions and boundary layer separation reattachment characteristics on a general aviation wing section to which was affixed wooden ice shapes which approximated typical glaze ice accretions. Comparisons were made with predicted pressure distributions using current airfoil analysis codes as well as the Bristow mixed analysis/design airfoil panel code. The Bristow code was also used to predict the separation reattachment dividing streamline by inputting the appropriate experimental surface pressure distribution.

Zaguli, R. J.↗

NASA NTRS1984

Implicit treatment of the unsteady full potential equation in conservation form

An implicit, conservative treatment for the unsteady full potential equation in two-dimensions is presented. The method employs a local time linearization for density, and introduces flux biasing concepts based on sonic conditions for the generation of artificial viscosity to capture shocks without any overshoots. The boundary condition is treated implicitly using a splitting procedure consistent with the approximate factorization scheme. This allows for extremely large Courant numbers, even for nonorthogonal grid at the body. The method has application not only to unsteady problems, but also to generate the starting blunt body solution for a supersonic full potential marching code. Results are presented for flows over cylinders, spheres and airfoils. Comparisons are made with available Euler and full potential results, and are in excellent agreement.

Shankar, V.↗

NASA NTRS1984

The prediction of transonic flows on advancing rotors

The object of this investigation was to find and demonstrate a means of performing efficient finite-difference computations of rotor loading for a trimmed rotor in high-speed, advancing flight. The essence of the scheme that was developed is a loose-coupled iteration procedure between a finite difference and a comprehensive integral rotor code. The coupling involves a transfer of appropriate load and inflow data between the two codes such that consistency is maintained. Sample computations, including a limited comparison with model rotor data, are presented. The scheme converges rapidly. The finite-difference computation of one rotor revolution with this scheme may provide sufficient accuracy for many purposes.

Tung, C.↗

NASA NTRS1989

Manipulation strategies for massive space payloads

Control for the bracing strategy is being examined. It was concluded earlier that trajectory planning must be improved to best achieve the bracing motion. Very interesting results were achieved which enable the inverse dynamics of flexible arms to be calculated for linearized motion in a more efficient manner than previously published. The desired motion of the end point beginning at t=0 and ending at t=t sub f is used to calculate the required torque at the joint. The solution is separated into a causal function that is zero for t is less than 0 and an accusal function which is zero for t is greater than t sub f. A number of alternative end point trajectories were explored in terms of the peak torque required, the amount of anticipatory action, and other issues. The single link case is the immediate subject and an experimental verification of that case is being performed. Modeling with experimental verification of closed chain dynamics continues. Modeling effort has pointed out inaccuracies that result from the choice of numerical techniques used to incorporate the closed chain constraints when modeling our experimental prototype RALF (Robotic Arm Large and Flexible). Results were compared to TREETOPS, a multi body code. The experimental verification work is suggesting new ways to make comparisons with systems having structural linearity and joint and geometric nonlinearity. The generation of inertial forces was studied with a small arm that will damp the large arm's vibration.

Book, Wayne J.↗

NASA NTRS1992

Direct numerical simulation of a combusting droplet with convection

The evaporation and combustion of a single droplet under forced and natural convection was studied numerically from first principles using a numerical scheme that solves the time-dependent multiphase and multispecies Navier-Stokes equations and tracks the sharp gas-liquid interface cutting across an arbitrary Eulerian grid. The flow fields both inside and outside of the droplet are resolved in a unified fashion. Additional governing equations model the interphase mass, energy, and momentum exchange. Test cases involving iso-octane, n-hexane, and n-propanol droplets show reasonable comparison rate, and flame stand-off distance. The partially validated code is, thus, readied to be applied to more demanding droplet combustion situations where substantial drop deformation render classical models inadequate.

Liang, Pak-Yan↗

NASA NTRS1993

HZEFRG1: An energy-dependent semiempirical nuclear fragmentation model

Methods for calculating cross sections for the breakup of high-energy heavy ions by the combined nuclear and coulomb fields of the interacting nuclei are presented. The nuclear breakup contributions are estimated with an abrasion-ablation model of heavy ion fragmentation that includes an energy-dependent, mean free path. The electromagnetic dissociation contributions arising from the interacting coulomb fields are estimated by using Weizsacker-Williams theory extended to include electric dipole and electric quadrupole contributions. The complete computer code that implements the model is included as an appendix. Extensive comparisons of cross section predictions with available experimental data are made.

Townsend, Lawrence W.↗

NASA NTRS1992

A NASA study of the impact of technology on future sea based attack aircraft

A conceptual aircraft design study was recently completed evaluating carrier-based, subsonic attack aircraft using contemporary and future technology assumptions. The study examined a configuration matrix that was made up of light and medium bomb loads, one and two man crews, internal and external weapons carriage, as well as conventional and flying wing planforms. Use of common technology assumptions, engine cycle simulation code, design mission, and consistent application of methods allow for direct comparison of the aircraft. This paper describes the design study ground rules and the aircraft designed. The aircraft descriptions include weights, dimensions, layout, design mission, design constraints, maneuver performance, and fallout mission performance. The strengths, and weaknesses of each aircraft are highlighted.

Hahn, Andrew S.↗

NASA NTRS1995

A NASTRAN-based computer program for structural dynamic analysis of Horizontal Axis Wind Turbines

This paper describes a computer program developed for structural dynamic analysis of horizontal axis wind turbines (HAWT's). It is based on the finite element method through its reliance on NASTRAN for the development of mass, stiffness, and damping matrices of the tower end rotor, which are treated in NASTRAN as separate structures. The tower is modeled in a stationary frame and the rotor in one rotating at a constant angular velocity. The two structures are subsequently joined together (external to NASTRAN) using a time-dependent transformation consistent with the hub configuration. Aerodynamic loads are computed with an established flow model based on strip theory. Aeroelastic effects are included by incorporating the local velocity and twisting deformation of the blade in the load computation. The turbulent nature of the wind, both in space and time, is modeled by adding in stochastic wind increments. The resulting equations of motion are solved in the time domain using the implicit Newmark-Beta integrator. Preliminary comparisons with data from the Boeing/NASA MOD2 HAWT indicate that the code is capable of accurately and efficiently predicting the response of HAWT's driven by turbulent winds.

Lobitz, Don W.↗

NASA NTRS1996

Spectroscopic Diagnostics and an Arc Jet Heated Air Plasma

Spectral radiation measurements were made in the range of 200 to 900 nm across a section of the plenum of an arc jet wind tunnel using a series of optical fibers. The spectra contained line radiation from Oxygen and Nitrogen atoms and molecular radiation from N2(+), N2, and NO. Abel inversion technique is used to obtain radial distribution of the spectra. The analysis yielded radial profiles of the electronic excitation, vibrational and rotational temperatures of the flow field. Spectral fitting yielded branching ratios for different vibrational and rotational bands. Relatively mild flow conditions, i.e. enthalpy and mass flow rate, were used for prolonged measurements of up to and over two hours to establish the best experimental methods of temperature determinations. Signal to noise was improved by at least an order of magnitude enabling the molecular vibrational band heads of N2(+) (first negative system), N2 (second positive system), and NO (beta, gamma, delta, and epsilon systems) to be resolved in the lower ultraviolet wavelength regions. The increased signal to noise ratio also enabled partial resolution of the rotational lines of N2(+) and N2 in certain regions of minimal overlap. Comparison of the spectra with theoretical models such as the NEQAIR2 code are presented and show potential for fitting the spectra when reliable calibration is performed for the complete wavelength range.

Mack, Larry Howard, Jr.↗

NASA NTRS2000

An Overview of the NCC Spray/Monte-Carlo-PDF Computations

This paper advances the state-of-the-art in spray computations with some of our recent contributions involving scalar Monte Carlo PDF (Probability Density Function), unstructured grids and parallel computing. It provides a complete overview of the scalar Monte Carlo PDF and Lagrangian spray computer codes developed for application with unstructured grids and parallel computing. Detailed comparisons for the case of a reacting non-swirling spray clearly highlight the important role that chemistry/turbulence interactions play in the modeling of reacting sprays. The results from the PDF and non-PDF methods were found to be markedly different and the PDF solution is closer to the reported experimental data. The PDF computations predict that some of the combustion occurs in a predominantly premixed-flame environment and the rest in a predominantly diffusion-flame environment. However, the non-PDF solution predicts wrongly for the combustion to occur in a vaporization-controlled regime. Near the premixed flame, the Monte Carlo particle temperature distribution shows two distinct peaks: one centered around the flame temperature and the other around the surrounding-gas temperature. Near the diffusion flame, the Monte Carlo particle temperature distribution shows a single peak. In both cases, the computed PDF's shape and strength are found to vary substantially depending upon the proximity to the flame surface. The results bring to the fore some of the deficiencies associated with the use of assumed-shape PDF methods in spray computations. Finally, we end the paper by demonstrating the computational viability of the present solution procedure for its use in 3D combustor calculations by summarizing the results of a 3D test case with periodic boundary conditions. For the 3D case, the parallel performance of all the three solvers (CFD, PDF, and spray) has been found to be good when the computations were performed on a 24-processor SGI Origin work-station.

Raju, M. S.↗

NASA NTRS2002

Albedo in the ATIC Experiment

ATIC(Advanced Thin Ionization Calorimeter) is a balloon borne experiment designed to measure the cosmic ray composition for elements from hydrogen to iron and their energy spectra from approx.50 GeV to near 100 TeV. It consists of a Si-matrix detector to determine the charge of a CR particle, a scintillator hodoscope for tracking, carbon interaction targets and a fully active BGO calorimeter. ATIC had its first 16-day flight from McMurdo, Antarctica from 28/12/2000 to 13/01/2000. The ATIC flight collected approximately 25 million events. To measure charge of primary particle in presence of radiation scattered back from the interaction and subsequent shower development in the calorimeter a charge detector must be a mosaic of small detector pads so that the pad containing the signal from the incident particle has no additional signal from albedo particles. Therefore the silicon matrix was built of 4480 individual silicon pads each 2 cm x 1.5 cm. The matrix consists of four planes of detectors and the active detector area, in these planes are partially overlapped to completely cover the aperture. The lateral and amplitude distributions of albedo signals in Si-matrix are analyzed for different primary nuclei and different energy deposits in BGO calorimeter. The greater part of albedo signals has Q near 1, where Q = square root of Amplitude(MIP). The albedo distribution exponentially decreases up to Q near 8. These high values are produced by slow protons and plans. There are also a small number of signals of Q > 8, mainly for heavy nucleus primaries. These signals are apparently generated by neutrons. The comparison of the experimental data and simulations with GEANT 3-21 code using QGSM generator for nucleus-nucleus interactions is presented.

Sokolskaya, N. V.↗

NASA NTRS2004

Evaluation of a CFD Method for Aerodynamic Database Development using the Hyper-X Stack Configuration

A computational fluid dynamic (CFD) study is performed on the Hyper-X (X-43A) Launch Vehicle stack configuration in support of the aerodynamic database generation in the transonic to hypersonic flow regime. The main aim of the study is the evaluation of a CFD method that can be used to support aerodynamic database development for similar future configurations. The CFD method uses the NASA Langley Research Center developed TetrUSS software, which is based on tetrahedral, unstructured grids. The Navier-Stokes computational method is first evaluated against a set of wind tunnel test data to gain confidence in the code s application to hypersonic Mach number flows. The evaluation includes comparison of the longitudinal stability derivatives on the complete stack configuration (which includes the X-43A/Hyper-X Research Vehicle, the launch vehicle and an adapter connecting the two), detailed surface pressure distributions at selected locations on the stack body and component (rudder, elevons) forces and moments. The CFD method is further used to predict the stack aerodynamic performance at flow conditions where no experimental data is available as well as for component loads for mechanical design and aero-elastic analyses. An excellent match between the computed and the test data over a range of flow conditions provides a computational tool that may be used for future similar hypersonic configurations with confidence.

Parikh, Paresh↗