Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “STRENGTH CALCULATIONS”

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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
NASA NTRS1980

Comparison between optical measurements and a numerical solution of the flow field within a transonic axial-flow compressor rotor

A comparison between numerical and experimental results is presented for the flowfield within a transonic axial-flow compressor rotor. The rotor was tested at design speed and a wide open throttle discharge condition. The relative tip Mach number was 1.4. A laser anemometer system was used to measure velocity and flow angle upstream, within, and downstream of the rotor. A holographic interferometer was used to visualize the rotor shock system near the tip. The computational procedure solves the full three-dimensional Euler equations using a time-marching technique. Shock location and shape determined from the two optical systems are compared. Calculated relative Mach number and flow angle contours, shock locations, and shock strength are compared to values measured with the laser anemometer.

Strazisar, A. J.↗
NASA NTRS1982

Thermal radiation view factor: Methods, accuracy and computer-aided procedures

The computer aided thermal analysis programs which predicts the result of predetermined acceptable temperature range prior to stationing of these orbiting equipment in various attitudes with respect to the Sun and the Earth was examined. Complexity of the surface geometries suggests the use of numerical schemes for the determination of these viewfactors. Basic definitions and standard methods which form the basis for various digital computer methods and various numerical methods are presented. The physical model and the mathematical methods on which a number of available programs are built are summarized. The strength and the weaknesses of the methods employed, the accuracy of the calculations and the time required for computations are evaluated. The situations where accuracies are important for energy calculations are identified and methods to save computational times are proposed. Guide to best use of the available programs at several centers and the future choices for efficient use of digital computers are included in the recommendations.

Kadaba, P. V.↗
NASA NTRS1983

Crater ejecta scaling laws - Fundamental forms based on dimensional analysis

Self-consistent scaling laws are developed for meteoroid impact crater ejecta. Attention is given to the ejection velocity of material as a function of the impact point, the volume of ejecta with a threshold velocity, and the thickness of ejecta deposit in terms of the distance from the impact. Use is made of recently developed equations for energy and momentum coupling in cratering events. Consideration is given to scaling of laboratory trials up to real-world events and formulations are developed for calculating the ejection velocities and ejecta blanket profiles in the gravity and strength regimes of crater formation. It is concluded that, in the gravity regime, the thickness of an ejecta blanket is the same in all directions if the thickness and range are expressed in terms of the crater radius. In the strength regime, however, the ejecta velocities are independent of crater size, thereby allowing for asymmetric ejecta blankets. Controlled experiments are recommended for the gravity/strength transition.

Housen, K. R.↗
NASA NTRS1984

Pair production rates in mildly relativistic, magnetized plasmas

Electron-positron pairs may be produced by either one or two photons in the presence of a strong magnetic field. In magnetized plasmas with temperatures kT approximately sq mc, both of these processes may be important and could be competitive. The rates of one-photon and two-photon pair production by photons with Maxwellian, thermal bremsstrahlung, thermal synchrotron and power law spectra are calculated as a function of temperature or power law index and field strength. This allows a comparison of the two rates and a determination of the conditions under which each process may be a significant source of pairs in astrophysical plasmas. It is found that for photon densities n(gamma) or = 10 to the 25th power/cu cm and magnetic field strengths B or = 10 to the 12th power G, one-photon pair production dominates at kT approximately sq mc for a Maxwellian, at kT approximately 2 sq mc for a thermal bremsstrahlung spectrum, at all temperatures for a thermal synchrotron spectrum, and for power law spectra with indices s approximately 4.

Burns, M. L.↗
NASA NTRS1984

Pair production rates in mildly relativistic magnetized plasmas

Electron-positron pairs may be produced by either one or two photons in the presence of a strong magnetic field. In magnetized plasmas with temperatures kT approximately sq mc, both of these processes may be important and could be competitive. The rates of one-photon and two-photon pair production by photons with Maxwellian, thermal bremsstrahlung, thermal synchrotron and power law spectra are calculated as a function of temperature or power law index and field strength. This allows a comparison of the two rates and a determination of the conditions under which each process may be a significant source of pairs in astrophysical plasmas. It is found that for photon densities n(gamma) or = 10 to the 25th power/cu cm and magnetic field strengths B or = 10 to the 12th power G, one-photon pair production dominates at kT approximately sq mc for a Maxwellian, at kT approximately 2 sq mc for a thermal bremsstrahlung spectrum, at all temperatures for a thermal synchrotron spectrum, and for power law spectra with indices s approximately 4.

Burns, M. L.↗
NASA NTRS1985

Spectral analysis of groove spacing on Ganymede

A quantitative analysis of groove spacing on Ganymede is described. Fourier transforms of a large number of photometric profiles across groove sets are calculated and the resulting power spectra are examined for the position and strength of peaks representing topographic periodicities. The geographic and global statistical distribution of groove wavelengths are examined, and these data are related to models of groove tectonism. It is found that groove spacing on Ganymede shows an approximately long-normal distribution with a minimum of about 3.5 km, a maximum of about 17 km, and a mean of 8.4 km. Groove spacing tends to be quite regular within a single groove set but can vary substantially from one groove set to another within a single geographic region.

Grimm, R. E.↗
NASA NTRS1986

Calculations of cosmogenic nuclide production rates in the Earth's atmosphere and their inventories

The production rates of cosmogenic isotopes in the Earth's atmosphere and their resulting terrestrial abundances have been calculated, taking into account both geomagnetic and solar-modulatory effects. The local interstellar flux was assumed to be that of Garcia-Munoz, et al. Solar modulation was accounted for using the heliocentric potential model and expressed in terms of the Deep River neutron monitor count rates. The geomagnetic field was presented by vertical cutoffs calculated by Shea and Smart and the non-vertical cutoffs calculated using ANGRI. The local interstellar particle flux was first modulated using the heliocentric potential field. The modulated cosmic-ray fluxes reaching the earth's orbit then interacted with the geomagnetic field as though it were a high-pass filter. The interaction of the cosmic radiation with the Earth's atmosphere was calculated utilizing the Bolztmann transport equation. Spallation cross sections for isotope production were calculated using the formalism of Silberberg and Tsao and other cross sections were taken from standard sources. Inventories were calculated by accounting from the variation in solar modulation and geomagnetic field strength with time. Results for many isotope, including C-14, Be-7 and Be-10 are in generally good agreement with existing data. The C-14 inventory, for instance, amounts to 1.75/sq cm(e)/s, in excellent agreement with direct estimates.

Obrien, K.↗
NASA NTRS1989

Modeling the classical nova outburst. I - Exploring the physics of a new mechanism

Model calculations were performed to describe a mechanism that produces classical nova outbursts on white dwarfs of 1 solar mass or less and for accretion rates of 4 x 10 to the -10th solar mass/yr or greater, i.e., the parameters corresponding to observed data of nova systems. Calculations point to four factors that can induce nuclear runaways of sufficient strength to eject about 0.0001 solar mass at speeds of several hundred to a few thousand km per second, as is observed in classical novae. These are (1) the effects of storage of angular momentum in the star's envelope during the accretion phase; (2) the reduction of centrifugal forces in the star's outer layers during the early nuclear runaway phase, through the inward transport of angular momentum; (3) the inward movement of the zone of peak nuclear burning through the convectively induced shear instability during the runaway phase; and (4) the mixing of original CO stellar matter and H-rich matter, also through the convectively induced shear instability.

Kutter, G. Siegfried↗
NASA NTRS1990

The free-wake computation of rotor-body flows

In this paper a method is described for predicting the compressible, free-wake, flow about a lifting rotor-body configuration. The method is an extension of a unique vorticity embedded full-potential method used to calculate free-wake rotor hover performance. An unusual feature of this method is that it obviates the requirement for multiple grids to treat the rotor-body problem. The approach used to treat the body is similar to that used to include the rotor wake in the full potential calculation. The body is modeled as a structured circulation sheet and the strength of this sheet is determined in an iterative manner. Initially the method is tested to compute the flow past simple isolated bodies like cylinders and spheres. After a comparison of these simple computations with exact solutions this procedure is included into the HELIX-I, free-wake rotor code, to compute the flow around a rotor mounted on a large whirl tower. The effects of the tower on rotor wake geometry and load distribution are presented.

Ramachandran, K.↗
NASA NTRS1992

Runge-Kutta methods combined with compact difference schemes for the unsteady Euler equations

Recent development using compact difference schemes to solve the Navier-Stokes equations show spectral-like accuracy. A study was made of the numerical characteristics of various combinations of the Runge-Kutta (RK) methods and compact difference schemes to calculate the unsteady Euler equations. The accuracy of finite difference schemes is assessed based on the evaluations of dissipative error. The objectives are reducing the numerical damping and, at the same time, preserving numerical stability. While this approach has tremendous success solving steady flows, numerical characteristics of unsteady calculations remain largely unclear. For unsteady flows, in addition to the dissipative errors, phase velocity and harmonic content of the numerical results are of concern. As a result of the discretization procedure, the simulated unsteady flow motions actually propagate in a dispersive numerical medium. Consequently, the dispersion characteristics of the numerical schemes which relate the phase velocity and wave number may greatly impact the numerical accuracy. The aim is to assess the numerical accuracy of the simulated results. To this end, the Fourier analysis is to provide the dispersive correlations of various numerical schemes. First, a detailed investigation of the existing RK methods is carried out. A generalized form of an N-step RK method is derived. With this generalized form, the criteria are derived for the three and four-step RK methods to be third and fourth-order time accurate for the non-linear equations, e.g., flow equations. These criteria are then applied to commonly used RK methods such as Jameson's 3-step and 4-step schemes and Wray's algorithm to identify the accuracy of the methods. For the spatial discretization, compact difference schemes are presented. The schemes are formulated in the operator-type to render themselves suitable for the Fourier analyses. The performance of the numerical methods is shown by numerical examples. These examples are detailed. described. The third case is a two-dimensional simulation of a Lamb vortex in an uniform flow. This calculation provides a realistic assessment of various finite difference schemes in terms of the conservation of the vortex strength and the harmonic content after travelling a substantial distance. The numerical implementation of Giles' non-refelctive equations coupled with the characteristic equations as the boundary condition is discussed in detail. Finally, the single vortex calculation is extended to simulate vortex pairing. For the distance between two vortices less than a threshold value, numerical results show crisp resolution of the vortex merging.

Yu, Sheng-Tao↗
NASA NTRS1994

Methodology for sensitivity analysis, approximate analysis, and design optimization in CFD for multidisciplinary applications

The straightforward automatic-differentiation and the hand-differentiated incremental iterative methods are interwoven to produce a hybrid scheme that captures some of the strengths of each strategy. With this compromise, discrete aerodynamic sensitivity derivatives are calculated with the efficient incremental iterative solution algorithm of the original flow code. Moreover, the principal advantage of automatic differentiation is retained (i.e., all complicated source code for the derivative calculations is constructed quickly with accuracy). The basic equations for second-order sensitivity derivatives are presented; four methods are compared. Each scheme requires that large systems are solved first for the first-order derivatives and, in all but one method, for the first-order adjoint variables. Of these latter three schemes, two require no solutions of large systems thereafter. For the other two for which additional systems are solved, the equations and solution procedures are analogous to those for the first order derivatives. From a practical viewpoint, implementation of the second-order methods is feasible only with software tools such as automatic differentiation, because of the extreme complexity and large number of terms. First- and second-order sensitivities are calculated accurately for two airfoil problems, including a turbulent flow example; both geometric-shape and flow-condition design variables are considered. Several methods are tested; results are compared on the basis of accuracy, computational time, and computer memory. For first-order derivatives, the hybrid incremental iterative scheme obtained with automatic differentiation is competitive with the best hand-differentiated method; for six independent variables, it is at least two to four times faster than central finite differences and requires only 60 percent more memory than the original code; the performance is expected to improve further in the future.

Taylor, Arthur C., III↗
NASA NTRS1993

Three-dimensional Navier-Stokes analysis and redesign of an imbedded bellmouth nozzle in a turbine cascade inlet section

Verification of proposed turbopump blading performance will involve evaluation of candidate blades in cascade test facilities. It is necessary to be able to predict the flow fields within these cascades for the results to be applicable to actual engine environments. This work presents the results of a study to predict the flow field for the NASA Lewis Transonic Turbine Blade Cascade Facility, which is similar to those used to evaluate rocket propulsion turbines. A pitchwise nonuniform total pressure distribution was observed at the blade row leading edge plane. A CFD analysis was used to show that the cause of the flow nonuniformity was a pair of vortices that originated in an embedded bellmouth inlet. Further CFD analysis was used to verify that a redesigned inlet section resulted in a flow with acceptable uniformity. A computational analysis was chosen because physical accessibility to the inlet section was limited, and because a computational approach also allows one to examine design changes cheaper and more quickly than an experimental approach would. The PARC code, a general purpose, three-dimensional, Navier-Stokes code with multiblock solution capability, was chosen for the present study. Results are presented detailing the computational requirements needed to accurately predict flows of this nature. Calculations of the original geometry showed total pressure loss regions consistent in strength and in location to experimental measurements. An examination of the results shows that the distortions are caused by a pair of vortices that originate as a result of the interaction of the flow with the imbedded bellmouth. Computations were performed for an inlet geometry which eliminated the imbedded bellmouth by bridging the region between it and the upstream wall. This analysis indicated that eliminating the imbedded bellmouth eliminates the troublesome pair of vortices, resulting in a flow with much greater pitchwise uniformity.

Giel, P. W.↗
NASA NTRS1998

Locations of Halide Ions in Tetragonal Lysozyme Crystals

Anions play an important role in the crystallization of lysozyme, and are known to bind to the crystalline protein. Previous studies employing X-ray crystallography had found one chloride ion binding site in the tetragonal crystal form of the protein and four nitrate ion binding sites in the monoclinic form. Studies using other approaches have reported more chloride ion binding sites, but their locations were not known. Knowing the precise location of these anions is also useful in determining the correct electrostatic fields surrounding the protein. In the first part of this study the anion positions in the tetragonal form were determined from the difference Fourier map obtained from the lysozyme crystals grown in bromide and chloride solutions under identical conditions. The anion locations were then obtained from standard crystallographic methods and five possible anion binding sites were found in this manner. The sole chloride ion binding site found in previous studies was confirmed. The remaining four sites were new ones for tetragonal lysozyme crystals. However, three of these new sites and the previously found one corresponded to the four unique binding sites found for nitrate ions in monoclinic crystals. This suggests that most of the anion binding sites in lysozyme remain unchanged, even when different anions and different crystal forms of lysozyme are employed. It is unlikely that there are many more anions in the tetragonal lysozyme crystal structure. Assuming osmotic equilibrium it can be shown that there are at most three more anions in the crystal channels. Some of the new anion binding sites found in this study were, as expected, in pockets containing basic residues. However, some of them were near neutral, but polar, residues. Thus, the study also showed the importance of uncharged, but polar groups, on the protein surface in determining its electrostatic field. This was important for the second part of this study where the electrostatic field surrounding the protein was accurately determined. This was achieved by solving the linearized version of the Poisson-Boltzmann equation for the protein in solution. The solution was computed employing the commercial code Delphi which uses a finite difference technique. This has recently become available as a module in the general protein visualization code Insight II. Partial charges were assigned to the polar groups of lysozyme for the calculations done here. The calculations showed the complexity of the electrostatic field surrounding the protein. Although most of the region near the protein surface had a positive field strength, the active site cleft was negatively charged and this was projected a considerable distance. This might explain the occurrence of "head-to-side" interactions in the formation of lysozyme aggregates in solution. Pockets of high positive field strength were also found in the vicinity of the anion locations obtained from the crystallographic part of this study, confirming the validity of these calculations. This study clearly shows not only the importance of determining the counterion locations in protein crystals and the electrostatic fields surrounding the protein, but also the advantage of performing them together.

Lim, Kap↗
NASA NTRS2002

Time Dependent Density Functional Theory Calculations of Large Compact PAH Cations: Implications for the Diffuse Interstellar Bands

We investigate the electronic absorption spectra of several maximally pericondensed polycyclic aromatic hydrocarbon radical cations with time dependent density functional theory calculations. We find interesting trends in the vertical excitation energies and oscillator strengths for this series containing pyrene through circumcoronene, the largest species containing more than 50 carbon atoms. We discuss the implications of these new results for the size and structure distribution of the diffuse interstellar band carriers.

Weisman, Jennifer L.↗
NASA NTRS2007

Analytical Approach for Estimating Preliminary Mass of ARES I Crew Launch Vehicle Upper Stage Structural Components

In January 2004, President Bush gave the National Aeronautics and Space Administration (NASA) a vision for Space Exploration by setting our sight on a bold new path to go back to the Moon, then to Mars and beyond. In response to this vision, NASA started the Constellation Program, which is a new exploration launch vehicle program. The primary mission for the Constellation Program is to carry out a series of human expeditions ranging from Low Earth Orbit to the surface of Mars and beyond for the purposes of conducting human exploration of space, as specified by the Vision for Space Exploration (VSE). The intent is that the information and technology developed by this program will provide the foundation for broader exploration activities as our operational experience grows. The ARES I Crew Launch Vehicle (CLV) has been designated as the launch vehicle that will be developed as a "first step" to facilitate the aforementioned human expeditions. The CLV Project is broken into four major elements: First Stage, Upper Stage Engine, Upper Stage (US), and the Crew Exploration Vehicle (CEV). NASA's Marshall Space Flight Center (MSFC) is responsible for the design of the CLV and has the prime responsibility to design the upper stage of the vehicle. The US is the second propulsive stage of the CLV and provides CEV insertion into low Earth orbit (LEO) after separation from the First Stage of the Crew Launch Vehicle. The fully integrated Upper Stage is a mix of modified existing heritage hardware (J-2X Engine) and new development (primary structure, subsystems, and avionics). The Upper Stage assembly is a structurally stabilized cylindrical structure, which is powered by a single J-2X engine which is developed as a separate Element of the CLV. The primary structure includes the load bearing liquid hydrogen (LH2) and liquid oxygen (LOX) propellant tanks, a Forward Skirt, the Intertank structure, the Aft Skirt and the Thrust Structure. A Systems Tunnel, which carries fluid and electrical power functions to other Elements of the CLV, is included as secondary structure. The MSFC has an overall responsibility for the integrated US element as well as structural design an thermal control of the fuel tanks, intertank, interstage, avionics, main propulsion system, Reaction Control System (RCS) for both the Upper Stage and the First Stage. MSFC's Spacecraft and Vehicle Department, Structural and Analysis Design Division is developing a set of predicted mass of these elements. This paper details the methodology, criterion and tools used for the preliminary mass predictions of the upper stage structural assembly components. In general, weight of the cylindrical barrel sections are estimated using the commercial code Hypersizer, whereas, weight of the domes are developed using classical solutions. HyperSizer is software that performs automated structural analysis and sizing optimization based on aerospace methods for strength, stability, and stiffness. Analysis methods range from closed form, traditional hand calculations repeated every day in industry to more advanced panel buckling algorithms. Margin-of-safety reporting for every potential failure provides the engineer with a powerful insight into the structural problem. Optimization capabilities include finding minimum weight panel or beam concepts, material selections, cross sectional dimensions, thicknesses, and lay-ups from a library of 40 different stiffened and sandwich designs and a database of composite, metallic, honeycomb, and foam materials. Multiple different concepts (orthogrid, isogrid, and skin stiffener) were run for multiple loading combinations of ascent design load with and with out tank pressure as well as proof pressure condition. Subsequently, selected optimized concept obtained from Hypersizer runs was translated into a computer aid design (CAD) model to account for the wall thickness tolerance, weld land etc for developing the most probable weight of the components. The flow diram summarizes the analysis steps used in developing these predicted mass.

Aggarwal, Pravin↗
NASA NTRS2010

Dst Index in the 2008 GEM Modeling Challenge - Model Performance for Moderate and Strong Magnetic Storms

The GEM 2008 modeling challenge efforts are expanding beyond comparing in-situ measurements in the magnetosphere and ionosphere to include the computation of indices to be compared. The Dst index measures the largest deviations of the horizontal magnetic field at 4 equatorial magnetometers from the quiet-time background field and is commonly used to track the strength of the magnetic disturbance of the magnetosphere during storms. Models can calculate a proxy Dst index in various ways, including using the Dessler-Parker Sckopke relation and the energy of the ring current and Biot-Savart integration of electric currents in the magnetosphere. The GEM modeling challenge investigates 4 space weather events and we compare models available at CCMC against each other and the observed values of Ost. Models used include SWMF/BATSRUS, OpenGGCM, LFM, GUMICS (3D magnetosphere MHD models), Fok-RC, CRCM, RAM-SCB (kinetic drift models of the ring current), WINDMI (magnetosphere-ionosphere electric circuit model), and predictions based on an impulse response function (IRF) model and analytic coupling functions with inputs of solar wind data. In addition to the analysis of model-observation comparisons we look at the way Dst is computed in global magnetosphere models. The default value of Dst computed by the SWMF model is for Bz the Earth's center. In addition to this, we present results obtained at different locations on the Earth's surface. We choose equatorial locations at local noon, dusk (18:00 hours), midnight and dawn (6:00 hours). The different virtual observatory locations reveal the variation around the earth-centered Dst value resulting from the distribution of electric currents in the magnetosphere during different phases of a storm.

Rastaetter, Lutz↗
NASA NTRS2014

An Interactive Method of Characteristics Java Applet to Design and Analyze Supersonic Aircraft Nozzles

The Method of Characteristics (MOC) is a classic technique for designing supersonic nozzles. An interactive computer program using MOC has been developed to allow engineers to design and analyze supersonic nozzle flow fields. The program calculates the internal flow for many classic designs, such as a supersonic wind tunnel nozzle, an ideal 2D or axisymmetric nozzle, or a variety of plug nozzles. The program also calculates the plume flow produced by the nozzle and the external flow leading to the nozzle exit. The program can be used to assess the interactions between the internal, external and plume flows. By proper design and operation of the nozzle, it may be possible to lessen the strength of the sonic boom produced at the rear of supersonic aircraft. The program can also calculate non-ideal nozzles, such as simple cone flows, to determine flow divergence and nonuniformities at the exit, and its effect on the plume shape. The computer program is written in Java and is provided as free-ware from the NASA Glenn central software server.

Analysis↗