Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical integration”

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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
NASA NTRS1998

Polar Hydra Data Analysis

The science activities are: 1) Hydra is still operating successfully on orbit. 2) A large amount of analysis and discovery has occurred with the Hydra ground data processing this past year. 3) Full interdetector calibration has been implemented and documented. This intercalibration was necessitated by the incorrect installation of bias resistors in the pre-acceleration stage to the electron channeltrons. This had the effect of making the counting efficiency for electrons energy dependent as well as channeltron specific. The nature of the error had no impact on the ion detection efficiency since they have a different bias arrangement. This intercalibration is so effective, that the electron and ion moment densities are routinely produced with a level of agreement better than 20%. 4) The data processing routinely removes glint in the sensors and produces public energy time spectrograms on the web overnight. 6) Routine, but more intensive computer processing codes are operational that determine for electrons and ions, the density, the flow vector, the pressure tensor and the heat flux by numerical integration. These codes use the magnetic field to sustain the quality of their output. To gain access to this high quality magnetic field within our data stream we have monitored Russell's web page for zero levels and timing files (since his data acquisition is not telemetry synchronous) and have a local reconstruction of B for our use. We have also detected a routine anomaly in the magnetometer data stream that we have documented to Chris Russell and developed an editing algorithm to intercept these "hits" and remove them from the geophysical analysis.

Scudder, J. D.↗
NASA NTRS1960

Temperature Histories in Ceramic-Insulated Heat-Sink Nozzle

Temperature histories were calculated for a composite nozzle wall by a simplified numerical integration calculation procedure. These calculations indicated that there is a unique ratio of insulation and metal heat-sink thickness that will minimize total wall thickness for a given operating condition and required running time. The optimum insulation and metal thickness will vary throughout the nozzle as a result of the variation in heat-transfer rate. The use of low chamber pressure results in a significant increase in the maximum running time of a given weight nozzle. Experimentally measured wall temperatures were lower than those calculated. This was due in part to the assumption of one-dimensional or slab heat flow in the calculation procedure.

Ciepluch, Carl C.↗
NASA NTRS1959

Structure of Weak Shock Waves in a Monatomic Gas

The profiles and thicknesses of normal shock waves in argon at Mach numbers of 1.335, 1.454, 1.576, and 1-713 were determined experimentally by means of a free-molecule probe whose equilibrium temperature is related by kinetic theory to the local flow properties and their gradients. Comparisons were made between the experimental shock profiles and the theoretical profiles calculated from the Navier-Stokes equations, the Grad 13-moment equations, and the Burnett equations. New, very accurate numerical integrations of the Burnett equations were obtained for this purpose with results quite different from those found by Zoller, to whom the solution of this problem is frequently attributed. The experimental shock profiles were predicted with approximately equal success by the Navier-Stokes and Burnett theories, while the 13-moment method was definitely less satisfactory. A surprising feature of the theoretical results is the relatively small difference in predictions between the Navier-Stokes and Burnett theories in the present range of shock strengths and the contrastingly large difference between predictions of Burnett and the 13-moment theories. It is concluded that the Navier-Stokes equations are correct for weak shocks and that within the present shock strength range the Burnett equations make no improvement which merits the trouble of solving them. For shocks of noticeably greater strength, say with a shock Mach number of more than 2.5, it remains fundamentally doubtful that any of these theories can be correct.

Sherman, F. S.↗
NASA NTRS1961

On the Long Period Luni-Solar Effect in the Motion of an Artificial Satellite

Two systems of formulas are presented for the determination of the long period perturbations caused by the Sun and the Moon in the motion of an artificial satellite. The first system can be used to determine the lunar effect for all satellites. The second method is more convenient for finding the lunar effect for close satellites and the solar effect for all satellites. Knowledge of these effects is essential for determining the stability of the satellite orbit. The basic equations of both systems are arranged in a form which permits the use of numerical integration. The two theories are more accurate and more adaptable to the use of electronic machines than the analytical developments obtained previously.

Musen, Peter↗
NASA NTRS1959

Approximate Analytical Solutions for Hypersonic Flow Over Slender Power Law Bodies

Approximate analytical solutions are presented for two-dimensional and axisymmetric hypersonic flow over slender power law bodies. Both zero order (M approaches infinity) and first order (small but nonvanishing values of 1/(M(Delta)(sup 2) solutions are presented, where M is free-stream Mach number and Delta is a characteristic slope. These solutions are compared with exact numerical integration of the equations of motion and appear to be accurate particularly when the shock is relatively close to the body.

Mirels, Harold↗
NASA NTRS1976

Nitric Oxide Formation by Meteoroids in the Upper Atmosphere

The process of nitric oxide formation during atmospheric entry of meteoroids is analyzed theoretically. An ablating meteoroid is assumed to be a point source in a uniform flow with a continuum regime evolving in its wake. The amount of nitric oxide produced by high-temperature reactions of air in the continuum regime is calculated by numerical integration of chemical-rate equations. This is accomplished by assuming that flow properties are constant across the reacting region, the radius of the region being determined from considerations of shock-wave formation and molecular diffusion. The results, when summed over the observed mass, velocity, and entry-angle distributions of meteoroids, provide annual global production rates of nitric oxide as a function of altitude. The peak production of nitric oxide is found to occur at altitudes between 9 x 10(exp 4) and 10(exp 5) m, the total annual rate being about 4 x 10(exp 7) kg. The present results suggest that the large concentration of nitric oxide observed below 9.5 x 10(exp 4) m could be attributed to meteoroids instead of photodissociation of nitrogen into metastable, 2D-state atoms, as has been previously hypothesized.

Menees, Gene P.↗
NASA NTRS1994

Analytic Sensitivities for Shape Optimization in Equivalent Plate Structural Wing Models

Equivalent plate modeling techniques based on Ritz analysis with simple polynomials prove to be efficient tools for structural modeling of wings in the preliminary design stage. Accuracy problems are encountered, however, when these models are used to obtain finite difference behavior sensitivities with respect to planform shape. The accuracy problems are associated with the poor numerical conditioning of static and eigenvalue equations. As higher-order polynomials are being used to Improve the analysis itself, the more sensitive is the finite difference derivative to the step size used. This article describes a formulation of wing equivalent plate modeling in which it is simple to obtain analytic, explicit expressions for stiffness and mass matrix elements without the need to perform numerical integration. This formulation leads naturally to analytic expressions for the derivatives of displacements, stresses, and natural frequencies with respect to shape design variables. This article examines the accuracy of finite difference derivatives compared with the analytic derivatives, and shows that In some cases it is impossible to obtain any information of value by finite differences. Analytic sensitivities, in this case, are still sufficiently accurate for design optimization.

Livne, Eli↗
NASA NTRS1999

A Parallel Compact Multi-Dimensional Numerical Algorithm with Aeroacoustics Applications

In this study we propose a novel method to parallelize high-order compact numerical algorithms for the solution of three-dimensional PDEs (Partial Differential Equations) in a space-time domain. For this numerical integration most of the computer time is spent in computation of spatial derivatives at each stage of the Runge-Kutta temporal update. The most efficient direct method to compute spatial derivatives on a serial computer is a version of Gaussian elimination for narrow linear banded systems known as the Thomas algorithm. In a straightforward pipelined implementation of the Thomas algorithm processors are idle due to the forward and backward recurrences of the Thomas algorithm. To utilize processors during this time, we propose to use them for either non-local data independent computations, solving lines in the next spatial direction, or local data-dependent computations by the Runge-Kutta method. To achieve this goal, control of processor communication and computations by a static schedule is adopted. Thus, our parallel code is driven by a communication and computation schedule instead of the usual "creative, programming" approach. The obtained parallelization speed-up of the novel algorithm is about twice as much as that for the standard pipelined algorithm and close to that for the explicit DRP algorithm.

Povitsky, Alex↗
NASA NTRS2000

A Real-Time Method for Estimating Viscous Forebody Drag Coefficients

This paper develops a real-time method based on the law of the wake for estimating forebody skin-friction coefficients. The incompressible law-of-the-wake equations are numerically integrated across the boundary layer depth to develop an engineering model that relates longitudinally averaged skin-friction coefficients to local boundary layer thickness. Solutions applicable to smooth surfaces with pressure gradients and rough surfaces with negligible pressure gradients are presented. Model accuracy is evaluated by comparing model predictions with previously measured flight data. This integral law procedure is beneficial in that skin-friction coefficients can be indirectly evaluated in real-time using a single boundary layer height measurement. In this concept a reference pitot probe is inserted into the flow, well above the anticipated maximum thickness of the local boundary layer. Another probe is servomechanism-driven and floats within the boundary layer. A controller regulates the position of the floating probe. The measured servomechanism position of this second probe provides an indirect measurement of both local and longitudinally averaged skin friction. Simulation results showing the performance of the control law for a noisy boundary layer are then presented.

Whitmore, Stephen A.↗
NASA NTRS1999

Taxiing, Take-Off, and Landing Simulation of the High Speed Civil Transport Aircraft

The aircraft industry jointly with NASA is studying enabling technologies for higher speed, longer range aircraft configurations. Higher speeds, higher temperatures, and aerodynamics are driving these newer aircraft configurations towards long, slender, flexible fuselages. Aircraft response during ground operations, although often overlooked, is a concern due to the increased fuselage flexibility. This paper discusses modeling and simulation of the High Speed Civil Transport aircraft during taxiing, take-off, and landing. Finite element models of the airframe for various configurations are used and combined with nonlinear landing gear models to provide a simulation tool to study responses to different ground input conditions. A commercial computer simulation program is used to numerically integrate the equations of motion and to compute estimates of the responses using an existing runway profile. Results show aircraft responses exceeding safe acceptable human response levels.

Reaves, Mercedes C.↗
NASA NTRS1999

A Sectored-One-Dimensional Model for Simulating Combustion Instabilities in Premix Combustors

A one-dimensional, CFD based combustor simulation has been developed that exhibits self-excited, thermoacoustic oscillations in premixed combustor geometries that typically have large, abrupt changes in cross sectional area. The combustor geometry is approximated by dividing it into a finite number of one-dimensional sectors. Within each sector, the equations of motion are integrated numerically, along with a species transport and a reaction equation. Across the sectors, mass and energy are conserved, and momentum loss is prescribed using appropriately compatible boundary conditions that account for the area change. The resulting simulation and associated boundary conditions essentially represent a one-dimensional, multi-block technique. Details of the simulation code are presented herein. Results are then shown comparing experimentally observed and simulated operation of a particular combustor rig that exhibited different instabilities at different operating points. It will be shown that the simulation closely matched the rig data in oscillation amplitudes, frequencies, and operating points at which the instabilities occurred. Finally, advantages and limitations of the simulation technique are discussed.

Paxson, Daniel E.↗
NASA NTRS2000

LSENS, The NASA Lewis Kinetics and Sensitivity Analysis Code

A general chemical kinetics and sensitivity analysis code for complex, homogeneous, gas-phase reactions is described. The main features of the code, LSENS (the NASA Lewis kinetics and sensitivity analysis code), are its flexibility, efficiency and convenience in treating many different chemical reaction models. The models include: static system; steady, one-dimensional, inviscid flow; incident-shock initiated reaction in a shock tube; and a perfectly stirred reactor. In addition, equilibrium computations can be performed for several assigned states. An implicit numerical integration method (LSODE, the Livermore Solver for Ordinary Differential Equations), which works efficiently for the extremes of very fast and very slow reactions, is used to solve the "stiff" ordinary differential equation systems that arise in chemical kinetics. For static reactions, the code uses the decoupled direct method to calculate sensitivity coefficients of the dependent variables and their temporal derivatives with respect to the initial values of dependent variables and/or the rate coefficient parameters. Solution methods for the equilibrium and post-shock conditions and for perfectly stirred reactor problems are either adapted from or based on the procedures built into the NASA code CEA (Chemical Equilibrium and Applications).

Radhakrishnan, K.↗
NASA NTRS2000

Radiation Heat Transfer Between Diffuse-Gray Surfaces Using Higher Order Finite Elements

This paper presents recent work on developing methods for analyzing radiation heat transfer between diffuse-gray surfaces using p-version finite elements. The work was motivated by a thermal analysis of a High Speed Civil Transport (HSCT) wing structure which showed the importance of radiation heat transfer throughout the structure. The analysis also showed that refining the finite element mesh to accurately capture the temperature distribution on the internal structure led to very large meshes with unacceptably long execution times. Traditional methods for calculating surface-to-surface radiation are based on assumptions that are not appropriate for p-version finite elements. Two methods for determining internal radiation heat transfer are developed for one and two-dimensional p-version finite elements. In the first method, higher-order elements are divided into a number of sub-elements. Traditional methods are used to determine radiation heat flux along each sub-element and then mapped back to the parent element. In the second method, the radiation heat transfer equations are numerically integrated over the higher-order element. Comparisons with analytical solutions show that the integration scheme is generally more accurate than the sub-element method. Comparison to results from traditional finite elements shows that significant reduction in the number of elements in the mesh is possible using higher-order (p-version) finite elements.

Gould, Dana C.↗
NASA NTRS1999

Recent Progress on the VASIMR

Experimental and theoretical studies on the Variable Specific Impulse Magnetoplasma Rocket (VASIMR) have continued through a NASA led collaborative program involving several research groups. In the experimental area, performance characterization of the VASIMR helicon plasma source has been obtained over a portion of the parameter space, with helium and hydrogen propellant. Density (10(exp 18) - 10(exp 19)/ cubic meter) and temperature (5 eV) were measured at moderate degree of ionization in two separate experimental devices. Helicon design improvement and optimization will be discussed. Experiments with the ion cyclotron resonance heating (ICRH) subsection have begun and preliminary results will be discussed. Theoretical picture and integrated numerical simulation continue to be refined to account for the main physics elements of the VASIMR, including RF absorption and particle acceleration with subsequent detachment in the magnetic nozzle.

ChangDiaz, Franklin R.↗
NASA NTRS1997

Issues in Defining Software Architectures in a GIS Environment

The primary mission of the Pan-American Center for Earth and Environmental Studies (PACES) is to advance the research areas that are relevant to NASA's Mission to Planet Earth program. One of the activities at PACES is the establishment of a repository for geographical, geological and environmental information that covers various regions of Mexico and the southwest region of the U.S. and that is acquired from NASA and other sources through remote sensing, ground studies or paper-based maps. The center will be providing access of this information to other government entities in the U.S. and Mexico, and research groups from universities, national laboratories and industry. Geographical Information Systems(GIS) provide the means to manage, manipulate, analyze and display geographically referenced information that will be managed by PACES. Excellent off-the-shelf software exists for a complete GIS as well as software for storing and managing spatial databases, processing images, networking and viewing maps with layered information. This allows the user flexibility in combining systems to create a GIS or to mix these software packages with custom-built application programs. Software architectural languages provide the ability to specify the computational components and interactions among these components, an important topic in the domain of GIS because of the need to integrate numerous software packages. This paper discusses the characteristics that architectural languages address with respect to the issues relating to the data that must be communicated between software systems and components when systems interact. The paper presents a background on GIS in section 2. Section 3 gives an overview of software architecture and architectural languages. Section 4 suggests issues that may be of concern when defining the software architecture of a GIS. The last section discusses the future research effort and finishes with a summary.

Acosta, Jesus↗
NASA NTRS2001

Numerical Simulation of Delamination Growth in Composite Materials

The use of decohesion elements for the simulation of delamination in composite materials is reviewed. The test methods available to measure the interfacial fracture toughness used in the formulation of decohesion elements are described initially. After a brief presentation of the virtual crack closure technique, the technique most widely used to simulate delamination growth, the formulation of interfacial decohesion elements is described. Problems related with decohesion element constitutive equations, mixed-mode crack growth, element numerical integration and solution procedures are discussed. Based on these investigations, it is concluded that the use of interfacial decohesion elements is a promising technique that avoids the need for a pre-existing crack and pre-defined crack paths, and that these elements can be used to simulate both delamination onset and growth.

Camanho, P. P.↗
NASA NTRS2002

The Fate of Ejecta from Hyperion

Ejecta from Saturn's moon Hyperion is subject to powerful perturbations from nearby Titan, which control its ultimate fate. We have performed numerical integrations to simulate a simplified system consisting of Saturn (including oblateness), Tethys, Dione, Titan, Hyperion, Iapetus, and the Sun (treated simply as a massive satellite). In addition, 1050 massless particles were ejected from Hyperion at five different points in its orbit. These particles started more or less evenly distributed over latitude and longitude, 1 km above Hyperion's mean radius, and were ejected radially outward at speeds 10\% faster than its escape speed. Only about 4\% of the particles survived for the 100,000-year course of the integration, while $\sim$8/% escaped from the Saturnian system. Titan accreted $\sim$77\% of all the particles, while Hyperion reaccreted only $\sim$5\%. This may help to account for Hyperion's rugged shape. Three particles hit Rhea and 2 hit Dione, but $\sim$5\% of the particles were removed when they penetrated within 150,000 km of Saturn. Most removals occurred within the first few thousand years. In general, ejecta from Hyperion are much more widely scattered than previously thought, and cross the orbits of all of the other classical satellites.

Lissauer, Jack J.↗
NASA NTRS2002

Formation and Migration of Trans-Neptunian Objects

We consider that trans-Neptunian objects (TNOs) with diameter d greater than 100 kilometers moving now in not very eccentric orbits could be formed directly by the compression of large rarefied dust condensations (with a greater than 30 AU), but not by the accretion of smaller solid planetesimals. Probably, some planetesimals with d approximately 100- 1000 kilometers in the feeding zone of the giant planets and even large main-belt asteroids also could be formed directly by such compression. Some smaller objects (TNOs, planetesimals, asteroids) could be debris of larger objects, and other such objects could be formed directly by compression of condensations. A small portion of planetesimals from the feeding zone of the giant planets that entered into the trans-Neptunian region could left in eccentrical orbits beyond Neptune and became so called "scattered objects". Our computer runs in which gravitational interactions of bodies were taken into account with the use of the spheres method, showed that the embryos of Uranus and Neptune could increase their semimajor axes from less than 10 AU to their present values, moving permanently in orbits with small eccentricities, due to gravitational interactions with the migrating planetesimals. Later on, Thommes et al. considered a similar model using direct numerical integration. The total amount of water delivered to the Earth during the formation of the giant, planets was about the mass of water in the Earth oceans. The end of such bombardment could be caused mainly by the planetesimals that became scattered objects. For the present mass of the trans-Neptunian belt, the collisional lifetime of 1-kilometer TNO is about the age of the Solar system, but only a small portion (less than 1%) of 100-km TNOs could be destroyed during this age. The probability of destruction of a typical TNO (with 30 less than a,less than 50 AU) by scattered objects can be of the same order of magnitude as that by typical TNOs. TNOs could be even more often destroyed during planet formation than during last 4 Gyr. The analysis of the results of the orbital evolution of Jupiter-crossing objects (JCOs) and TNOs showed that, in principle, the trans-Neptunian belt can provide up to 100% of Earth-crossing objects, but, of course, some of them came from the main asteroid belt. Most of the collisions of former JCOs with the Earth were from orbits with aphelia inside Jupiter's orbit and belonged to a small portion of objects which moved in such orbits for a long time. About 1 of 300 JCOs collided with the Sun. The ratio of the total mass of icy planetesimals that migrated from the feeding zone of the giant planets and collided with the planet to the mass of this planet was greater (by a factor of 3 in our runs) for Mars than that for Earth and Venus.

Ipatov, S. I.↗