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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
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.↗
DOE OSTI2022

Numerical methods and hypoexponential approximations for gamma distributed delay differential equations

Abstract Gamma distributed delay differential equations (DDEs) arise naturally in many modelling applications. However, appropriate numerical methods for generic gamma distributed DDEs have not previously been implemented. Modellers have therefore resorted to approximating the gamma distribution with an Erlang distribution and using the linear chain technique to derive an equivalent system of ordinary differential equations (ODEs). In this work, we address the lack of appropriate numerical tools for gamma distributed DDEs in two ways. First, we develop a functional continuous Runge–Kutta (FCRK) method to numerically integrate the gamma distributed DDE without resorting to Erlang approximation. We prove the fourth-order convergence of the FCRK method and perform numerical tests to demonstrate the accuracy of the new numerical method. Nevertheless, FCRK methods for infinite delay DDEs are not widely available in existing scientific software packages. As an alternative approach to solving gamma distributed DDEs, we also derive a hypoexponential approximation of the gamma distributed DDE. This hypoexponential approach is a more accurate approximation of the true gamma distributed DDE than the common Erlang approximation but, like the Erlang approximation, can be formulated as a system of ODEs and solved numerically using standard ODE software. Using our FCRK method to provide reference solutions, we show that the common Erlang approximation may produce solutions that are qualitatively different from the underlying gamma distributed DDE. However, the proposed hypoexponential approximations do not have this limitation. Finally, we apply our hypoexponential approximations to perform statistical inference on synthetic epidemiological data to illustrate the utility of the hypoexponential approximation.

97 MATHEMATICS AND COMPUTING↗
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.↗
NASA NTRS1991

Benchmark calculations of thermal reaction rates. I - Quantal scattering theory

The thermal rate coefficient for the prototype reaction H + H2 yields H2 + H with zero total angular momentum is calculated by summing, averaging, and numerically integrating state-to-state reaction probabilities calculated by time-independent quantum-mechanical scattering theory. The results are very carefully converged with respect to all numerical parameters in order to provide high-precision benchmark results for confirming the accuracy of new methods and testing their efficiency.

Chatfield, David C.↗
NASA NTRS1994

Hot-Gas-Slide and Coolant-Side Heat Transfer in Liquid Rocket Engine Combustors

The objectives of this article are to develop a multidisciplinary, computational methodology to predict the hot-gas-side and coolant-side heat transfer in film cooling assisted, regeneratively cooled liquid rocket engine combustors, and to use it in parametric studies to recommend optimized design of the coolant channels for a developmental combustor. An integrated numerical model which incorporates computational fluid dynamics (CFD) for the hot-gas thermal environment, and thermal analysis for the liner and coolant channels, was developed. This integrated CFD/thermal model was validated by comparing predicted heat fluxes with those of hot-firing test and industrial design methods for a 40-k calorimeter thrust chamber and the Space Shuttle Main Engine main combustion chamber. Parametric studies were performed for the advanced main combustion chamber to find a strategy for a proposed coolant channel design.

Ten-see Wang↗
NASA NTRS1985

Accumulation of the Terrestrial Planets

The validity of the two body approximation in calculations of planetary growth is examined. Use of this approximation is essential to practical 3 dimensional simulations of planetary growth. With regard to gravitational perturbations, the changes in semimajor axis, eccentricity, and inclination resulting from close planetesimal encounters (near 1 A.U.) out to 10 Tisserand sphere of influence radii were calculated by two and three dimensional numerical integration. These were compared with the results of treating the encounter as a two body problem, as is customary in Monte Carlo calculations of orbital evolution and in numerical and analytical studies of planetary accumulation. It is found that for values of (V/V sub e) approx. 0.35 (V = relative velocity, V = escape velocity of largest body), the two body approximation fails to describe the outcome of individual encounters. In this low velocity region, the two body gravitational focusing cross section is no longer valid.

Wetherill, G. W.↗