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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

Transient rotor dynamic rub phenomena - Theory and test

This paper develops an implicit integration scheme for transient rotor dynamic rub prediction and includes a correlation study with actual test results. A Nordsieck-like numerical integration scheme is applied directly to the second-order equations of motion. The assumption that forces and torques on the rotor are functions of the position and velocity at the point of application and its nearest neighbor is made in order to make the computational time proportional to the number of elements in the rotor dynamics model rather than the cube of the number. The test rig consists of a turbine driven, flexible shaft supported by squeeze film dampers. The blade loss event occurs due to collision of a balance bolt on one of the disks with a high speed plunger. The rotor is seen to spiral outward and contact against a stationary assemblage of seal shoes.

Kascak, A. F.↗
NASA NTRS1989

Modeling of the Earth's gravity field using the New Global Earth Model (NEWGEM)

Traditionally, the global gravity field was described by representations based on the spherical harmonics (SH) expansion of the geopotential. The SH expansion coefficients were determined by fitting the Earth's gravity data as measured by many different methods including the use of artificial satellites. As gravity data have accumulated with increasingly better accuracies, more of the higher order SH expansion coefficients were determined. The SH representation is useful for describing the gravity field exterior to the Earth but is theoretically invalid on the Earth's surface and in the Earth's interior. A new global Earth model (NEWGEM) (KIM, 1987 and 1988a) was recently proposed to provide a unified description of the Earth's gravity field inside, on, and outside the Earth's surface using the Earth's mass density profile as deduced from seismic studies, elevation and bathymetric information, and local and global gravity data. Using NEWGEM, it is possible to determine the constraints on the mass distribution of the Earth imposed by gravity, topography, and seismic data. NEWGEM is useful in investigating a variety of geophysical phenomena. It is currently being utilized to develop a geophysical interpretation of Kaula's rule. The zeroth order NEWGEM is being used to numerically integrate spherical harmonic expansion coefficients and simultaneously determine the contribution of each layer in the model to a given coefficient. The numerically determined SH expansion coefficients are also being used to test the validity of SH expansions at the surface of the Earth by comparing the resulting SH expansion gravity model with exact calculations of the gravity at the Earth's surface.

Kim, Yeong E.↗
NASA NTRS2009

Mean Element Propagations Using Numerical Averaging

The long-term evolution characteristics (and stability) of an orbit are best characterized using a mean element propagation of the perturbed two body variational equations of motion. The averaging process eliminates short period terms leaving only secular and long period effects. In this study, a non-traditional approach is taken that averages the variational equations using adaptive numerical techniques and then numerically integrating the resulting EOMs. Doing this avoids the Fourier series expansions and truncations required by the traditional analytic methods. The resultant numerical techniques can be easily adapted to propagations at most solar system bodies.

orbit propagations↗
NASA NTRS1984

A computer program for predicting nonlinear uniaxial material responses using viscoplastic models

A computer program was developed for predicting nonlinear uniaxial material responses using viscoplastic constitutive models. Four specific models, i.e., those due to Miller, Walker, Krieg-Swearengen-Rhode, and Robinson, are included. Any other unified model is easily implemented into the program in the form of subroutines. Analysis features include stress-strain cycling, creep response, stress relaxation, thermomechanical fatigue loop, or any combination of these responses. An outline is given on the theoretical background of uniaxial constitutive models, analysis procedure, and numerical integration methods for solving the nonlinear constitutive equations. In addition, a discussion on the computer program implementation is also given. Finally, seven numerical examples are included to demonstrate the versatility of the computer program developed.

Chang, T. Y.↗
NASA NTRS1981

A study of the viscous and nonadiabatic flow in radial turbines

A method for analyzing the viscous nonadiabatic flow within turbomachine rotors is presented. The field analysis is based upon the numerical integration of the incompressible Navier-Stokes equations together with the energy equation over the rotors blade-to-blade stream channels. The numerical code used to solve the governing equations employs a nonorthogonal boundary fitted coordinate system that suits the most complicated blade geometries. Effects of turbulence are modeled with two equations; one expressing the development of the turbulence kinetic energy and the other its dissipation rate. The method of analysis is applied to a radial inflow turbine. The solution obtained indicates the severity of the complex interaction mechanism that occurs between different flow regimes (i.e., boundary layers, recirculating eddies, separation zones, etc.). Comparison with nonviscous flow solutions tend to justify strongly the inadequacy of using the latter with standard boundary layer techniques to obtain viscous flow details within turbomachine rotors. Capabilities and limitations of the present method of analysis are discussed.

Khalil, I.↗
NASA NTRS1973

Quantum theory of an optical maser. VI - Transient behavior.

The transient behavior of a laser is discussed using the quantum theory as did Scully and Lamb. The formal solution of the density-matrix equation is expressed in terms of exponentially decaying eigenmodes. Some of the lower decay constants are obtained numerically. The equations for the moments of the density matrix are then derived and solved by a truncation method. The equations of motion are integrated numerically for the case where the average number of photons in a laser cavity has the realistically large value 1.3 x 100,000. An alternative Fokker-Planck-equation approach is discussed.

Wang, Y. K.↗
NASA NTRS1986

Quantitative modeling of the 'convection surge' mechanism of ion acceleration

To explore the properties of the 'convection surge' mechanism postulated by Quinn and Southwood (1982), a quantitative computer model of the mechanism is constructed. The model incorporates the dipolarization of the field line shapes as well as the results of recent transient electric field measurements by the AST-6 satellite. The computer model numerically integrates the equations of motion of ensembles of ions populating a single flux tube that goes through the convection surge or 'dipolarization' transformation. The numerical results confirm that the mechanism easily generates the bounce-phase-bunched ion distributions in question, and show that the convection surge can generate dramatically field-aligned ion distributions. It is hypothesized that the convection surge mechanism is fundamentally associated with the processes responsible for transporting tail ion populations to the geosynchronous regions of the earth's magnetosphere.

Mauk, B. H.↗
NASA NTRS1993

Multigrid method for integral equations and automatic programs

Several iterative algorithms based on multigrid methods are introduced for solving linear Fredholm integral equations of the second kind. Automatic programs based on these algorithms are introduced using Simpson's rule and the piecewise Gaussian rule for numerical integration.

Lee, Hosae↗
NASA NTRS2015

Numerical Computation of a Continuous-thrust State Transition Matrix Incorporating Accurate Hardware and Ephemeris Models

A significant body of work exists showing that providing a nonlinear programming (NLP) solver with expressions for the problem constraint gradient substantially increases the speed of program execution and can also improve the robustness of convergence, especially for local optimizers. Calculation of these derivatives is often accomplished through the computation of spacecraft's state transition matrix (STM). If the two-body gravitational model is employed as is often done in the context of preliminary design, closed form expressions for these derivatives may be provided. If a high fidelity dynamics model, that might include perturbing forces such as the gravitational effect from multiple third bodies and solar radiation pressure is used then these STM's must be computed numerically. We present a method for the power hardward model and a full ephemeris model. An adaptive-step embedded eight order Dormand-Prince numerical integrator is discussed and a method for the computation of the time of flight derivatives in this framework is presented. The use of these numerically calculated derivatieves offer a substantial improvement over finite differencing in the context of a global optimizer. Specifically the inclusion of these STM's into the low thrust missiondesign tool chain in use at NASA Goddard Spaceflight Center allows for an increased preliminary mission design cadence.

Low Thrust↗
NASA NTRS1973

The particle resonance in spiral galaxies - Nonlinear effects.

A theory is developed to account for the nonlinear effects, near the particle resonance, found by numerical integration. There are four equilibrium points in the rotating frame of reference: two of them unstable, at the minima of potential (L1, L2), and two stable, at the maxima of potential (L4, L5). Many particles are trapped in librating orbits around L4, L5. Using the lowest-order terms of a 'third integral' of motion near L4, L5, the behavior of the trapped orbits is determined, and the approximate theoretical results are compared with the orbits found numerically by computer. An estimate of the trapped mass is given. It amounts to about 30% of the total mass between 9 and 11 kpc when the particle resonance is at r sub s = 10 kpc and the maximum force due to the spiral is 4% of the axisymmetric force. Some effects due to these mass concentrations are found numerically and theoretically.

Contopoulos, G.↗
NASA NTRS1992

Radiative transfer in the electrode region of an arc heater

Radiation from hot gases in the electrode region of an arc heater are discussed. Major original contributions are made: (1) an analytical expression for the radiative heat flux, for an axisymmetric geometry and a grey reflecting wall, is derived, (2) the expression is rendered analytically differentiable/integrable by using appropriate techniques to remove the singularities, and (3), finally a two band model is devised to compute the radiative heat load at the walls in NASA Ames' jet facility. The advantage of using these techniques is that the time taken for the numerical integration is considerably reduced and computation of a 39 x 39 axisymmetric geometry becomes manageable. Fluid equations, representative of a high current arc heater, are numerically solved for an axisymmetric geometry. The computed average heat load of 400 Watts/sq cm using the two band model agrees well with the experimental estimates made at NASA Ames' facility.

Durgapal, P.↗
NASA NTRS2002

Two- and Three-Dimensional Numerical Experiments Representing Two Limiting Cases of an In-Line Pair of Finger Seal Components

The work presented here concerns the numerical development and simulation of the flow, pressure patterns and motion of a pair of fingers arranged behind each other and axially aligned in-line. The fingers represent the basic elemental component of a Finger Seal (FS) and form a tight seal around the rotor. Yet their flexibility allows compliance with rotor motion and in a passive-adaptive mode complies also with the hydrodynamic forces induced by the flowing fluid. While the paper does not treat the actual staggered configuration of a finger seal, the inline arrangement represents a first step towards that final goal. The numerical 2-D (axial-radial) and 3-D results presented herein were obtained using a commercial package (CFD-ACE+). Both models use an integrated numerical approach, which couples the hydrodynamic fluid model (Navier-Stokes based) to the solid mechanics code that models the compliance of the fingers.

Braun, M. J.↗
NASA NTRS1991

Numerical simulation of spatially-evolving instability in plane channel flow

The spatial stability of plane channel flow is analyzed using a three-dimensional, time-dependent spectral/finite difference code (Danabasoglu et al., 1990) which integrates numerically the Navier-Stokes equations. The study centers on inflow disturbance amplitudes effects on the secondary instability. The resolution requirements along the spacewise direction, which become critical before the breakdown stage, are of particular interest. A direct comparison is made with the experiments of Nishioka et al. (1980).

Danabasoglu, G.↗
NASA NTRS1987

Effect of atmospheric turbulence on the bit error probability of a space to ground near infrared laser communications link using binary pulse position modulation and an avalanche photodiode detector

The effect of atmospheric turbulence on the bit error rate of a space-to-ground near infrared laser communications link is investigated, for a link using binary pulse position modulation and an avalanche photodiode detector. Formulas are presented for the mean and variance of the bit error rate as a function of signal strength. Because these formulas require numerical integration, they are of limited practical use. Approximate formulas are derived which are easy to compute and sufficiently accurate for system feasibility studies, as shown by numerical comparison with the exact formulas. A very simple formula is derived for the bit error rate as a function of signal strength, which requires only the evaluation of an error function. It is shown by numerical calculations that, for realistic values of the system parameters, the increase in the bit error rate due to turbulence does not exceed about thirty percent for signal strengths of four hundred photons per bit or less. The increase in signal strength required to maintain an error rate of one in 10 million is about one or two tenths of a db.

Safren, H. G.↗
NASA NTRS2015

Experiments with Julia for Astrodynamics Applications

Julia's potential for solving complex astrodynamics problems is studied. Julia is a high-level, new, dynamic programming language with performance approaching C/Fortran and has features like inbuilt parallelism, variable accuracy, integrated numerical libraries and direct C and Fortran interfaces. Two astrodynamics problems are solved in Julia: 1) Lambert's problem, using the vercosine formulation and 2) trajectory integration. Implemented algorithms are compared with C and Fortran based counterparts on key performance parameters (speed, development effort, etc). Using Julia for fast and reliable astrodynamics software development is also discussed.

Arora, Nitin↗
NASA NTRS1974

Surface electrical properties experiment. Part 2: Theory of radio-frequency interferometry in geophysical subsurface probing

The radiation fields due to a horizontal electric dipole laid on the surface of a stratified medium were calculated using a geometrical optics approximation, a modal approach, and direct numerical integration. The solutions were obtained from the reflection coefficient formulation and written in integral forms. The calculated interference patterns are compared in terms of the usefulness of the methods used to obtain them. Scattering effects are also discussed and all numerical results for anisotropic and isotropic cases are presented.

Kong, J. A.↗
NASA NTRS1986

Turbulent two-phase flow in annular seals

Steady, turbulent two-phase fluid flow in a rotating annular seal with no eccentricity is analyzed. The fluid is assumed to be a homogeneous mixture of liquid and vapor in thermodynamic equilibrum. Further, the flow is assumed to be adiabatic, but the effects due to heat generation by viscous dissipation are accounted for fully. Solution of the model governing differential equations is accomplished by use of a fourth-order Runge-Kutta numerical integration scheme. The calculation of mass leakage rates under choked and unchoked conditions are discussed and the phenomenon of all-liquid choked flow is explained. Several numerical examples are presented supposing cyrogenic oxygen as the sealed fluid.

Beatty, P. A.↗