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.

SourceResearch connectionFrom yearThrough year

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

At least 379 records · Page 21

All sources NASA NTRS DOE OSTI

NASA NTRS1991

Surface albedo from bidirectional reflectance

The validity of integrating over discrete wavelength bands is examined to estimate total shortwave bidirectional reflectance of vegetated and bare soil surfaces. Methods for estimating albedo from multiple angle, discrete wavelength band radiometer measurements are studied. These methods include a numerical integration technique and the integration of an empirically derived equation for bidirectional reflectance. It is concluded that shortwave albedos estimated through both techniques agree favorably with the independent pyranometer measurements. Absolute rms errors are found to be 0.5 percent or less for both grass sod and bare soil surfaces.

Ranson, K. J.↗

NASA NTRS1972

A study of nozzle and ejector flow problems by the method of integral relations

The application of the method of integral relations to nozzle and ejector flow problems was examined. For nozzle flow problems, the general formulation is that the approaching flow may be rotational. Particular attention was given to the phenomenon of choking under nonuniform flow conditions. Numerical integration of the governing ordinary differential equations was also investigated. This scheme of analysis was also applied to study the interacting flow field within an ejector system.

Howlett, L. D.↗

NASA NTRS1979

Aperture amplitude and phase control of offset dual reflectors

An approximate synthesis method is presented for solving the problem of offset dual-shaped reflectors. It is shown that the equations obtained by this method are partial differential equations that can be integrated numerically, except that these equations are generally not total and therefore do not have a 'smooth' solution. However, it is demonstrated that these equations form a nearly total differential that enables one to generate a smooth subreflector by integrating the differential equations and then synthesizing a main reflector which gives an exact solution for the specified aperture phase distribution. Several important illustrative examples are treated, which validate the adequacy of the proposed method.

Galindo-Israel, V.↗

NASA NTRS1990

Numerical derivation of forced nutation terms for a rigid earth

The results of a numerical integration of the Euler equations for a rigid earth model covering a time span of 250 years are compared with Kinoshita's (1977) theory for the forced nutations and with a new nutation series by Kinoshita and Souchay. Numerical corrections to some of the analytically derived nutation terms are presented.

Schastok, J.↗

NASA NTRS1992

A general numerical model for wave rotor analysis

Wave rotors represent one of the promising technologies for achieving very high core temperatures and pressures in future gas turbine engines. Their operation depends upon unsteady gas dynamics and as such, their analysis is quite difficult. This report describes a numerical model which has been developed to perform such an analysis. Following a brief introduction, a summary of the wave rotor concept is given. The governing equations are then presented, along with a summary of the assumptions used to obtain them. Next, the numerical integration technique is described. This is an explicit finite volume technique based on the method of Roe. The discussion then focuses on the implementation of appropriate boundary conditions. Following this, some results are presented which first compare the numerical approximation to the governing differential equations and then compare the overall model to an actual wave rotor experiment. Finally, some concluding remarks are presented concerning the limitations of the simplifying assumptions and areas where the model may be improved.

Paxson, Daniel W.↗

NASA NTRS1992

Chaotic evolution of the solar system

The evolution of the entire planetary system has been numerically integrated for a time span of nearly 100 million years. This calculation confirms that the evolution of the solar system as a whole is chaotic, with a time scale of exponential divergence of about 4 million years. Additional numerical experiments indicate that the Jovian planet subsystem is chaotic, although some small variations in the model can yield quasi-periodic motion. The motion of Pluto is independently and robustly chaotic.

Sussman, Gerald J.↗

NASA NTRS1991

Distributed Minimal Residual (DMR) method for acceleration of iterative algorithms

A new method for enhancing the convergence rate of iterative algorithms for the numerical integration of systems of partial differential equations was developed. It is termed the Distributed Minimal Residual (DMR) method and it is based on general Krylov subspace methods. The DMR method differs from the Krylov subspace methods by the fact that the iterative acceleration factors are different from equation to equation in the system. At the same time, the DMR method can be viewed as an incomplete Newton iteration method. The DMR method was applied to Euler equations of gas dynamics and incompressible Navier-Stokes equations. All numerical test cases were obtained using either explicit four stage Runge-Kutta or Euler implicit time integration. The formulation for the DMR method is general in nature and can be applied to explicit and implicit iterative algorithms for arbitrary systems of partial differential equations.

Lee, Seungsoo↗

NASA NTRS1991

A Numerical Study of Hypersonic Forebody/Inlet Integration Problem

A numerical study of hypersonic forebody/inlet integration problem is presented in the form of the view-graphs. The following topics are covered: physical/chemical modeling; solution procedure; flow conditions; mass flow rate at inlet face; heating and skin friction loads; 3-D forebogy/inlet integration model; and sensitivity studies.

Kumar, Ajay↗

NASA NTRS1979

Study of periodic motions of a satellite with a magnetic damper

The motion of a satellite with a magnetic damper in the plane of a circular polar orbit is studied. The asymptotics of periodic solutions are constructed for a satellite close to axisymmetric and the radius of convergence is evaluated for the power series obtained. In a broad range of values of parameters, a periodic solution is obtained by numerical integration of equations of motion of the satellite. The asymptotics of a bifurcated curve obtained (the curve on which origin of a pair of periodic solutions occurs) in the space of the parameters agrees well with the results of numerical computation with all physical values of these parameters. A breakdown is made of the space of the initial data of phase variables in the field of effect of different types of periodic motion.

Sadov, Y. A.↗

NASA NTRS1990

Accurate calculation and instability of supersonic wake flows

This study is concerned with the computation and linear stability of a class of laminar compressible wake flows. The emphasis is on correct basic flow profiles that satisfy the steady equations of motion, and to this end the unperturbed state is obtained through numerical integration of the compressible boundary-layer equations. The linear stability of the flow is examined via the Rayleigh equation that describes evolution of inviscid disturbances. Analytical results are given for short- and long-wavelength disturbances and some numerical results of the general eigenvalue problem are also reported.

Papageorgiou, Demetrius T.↗

NASA NTRS1995

A case study of view-factor rectification procedures for diffuse-gray radiation enclosure computations

The view factors which are used in diffuse-gray radiation enclosure calculations are often computed by approximate numerical integrations. These approximately calculated view factors will usually not satisfy the important physical constraints of reciprocity and closure. In this paper several view-factor rectification algorithms are reviewed and a rectification algorithm based on a least-squares numerical filtering scheme is proposed with both weighted and unweighted classes. A Monte-Carlo investigation is undertaken to study the propagation of view-factor and surface-area uncertainties into the heat transfer results of the diffuse-gray enclosure calculations. It is found that the weighted least-squares algorithm is vastly superior to the other rectification schemes for the reduction of the heat-flux sensitivities to view-factor uncertainties. In a sample problem, which has proven to be very sensitive to uncertainties in view factor, the heat transfer calculations with weighted least-squares rectified view factors are very good with an original view-factor matrix computed to only one-digit accuracy. All of the algorithms had roughly equivalent effects on the reduction in sensitivity to area uncertainty in this case study.

Taylor, Robert P.↗

DOE OSTI2025

$\mathrm{SageNet}$: Fast Neural Network Emulation of the Stiff-amplified Gravitational Waves from Inflation

Accurate modeling of the inflationary gravitational waves (GWs) requires time-consuming, iterative numerical integrations of differential equations to take into account their backreaction on the expansion history. To improve computational efficiency while preserving accuracy, we present the Stiff-amplified Gravitational-wave Emulator Network (SageNet), a deep learning framework designed to replace conventional numerical solvers (code available at https://github.com/YifangLuo/SageNet). SageNet employs a long short-term memory architecture to emulate the present-day energy density spectrum of the inflationary GWs with possible stiff amplification, Ω GW (f). Trained on a data set of 25,689 numerically generated solutions, SageNet allows accurate reconstructions of Ω GW (f) and generalizes well to a wide range of cosmological parameters; 90.9% of the test emulations with randomly distributed parameters exhibit errors of under 4%. In addition, SageNet demonstrates its ability to learn and reproduce the artificial, adaptive sampling patterns in numerical calculations, which implement denser sampling of frequencies around changes in spectral indices in Ω GW (f). The dual capability of learning both physical and artificial features of the numerical GW spectra establishes SageNet as a robust alternative to exact numerical methods. Finally, our benchmark tests show that SageNet reduces the computation time from tens of seconds to milliseconds, achieving a speedup of ∼10 4 times over standard CPU-based numerical solvers with the potential for further acceleration on GPU hardware. These capabilities make SageNet a powerful tool for accelerating Bayesian inference procedures for extended cosmological models. In a broad sense, the SageNet framework offers a fast, accurate, and generalizable solution to modeling cosmological observables whose theoretical predictions demand costly differential equation solvers.

Astronomy data modeling↗

NASA NTRS1974

Program manual for HILTOP, a heliocentric interplanetary low thrust trajectory optimization program. Part 1: User's guide

A performance-analysis computer program, that was developed explicitly to generate optimum electric propulsion trajectory data for missions of interest in the exploration of the solar system is presented. The program was primarily designed to evaluate the performance capabilities of electric propulsion systems, and in the simulation of a wide variety of interplanetary missions. A numerical integration of the two-body, three-dimensional equations of motion and the Euler-Lagrange equations was used in the program. Transversality conditions which permit the rapid generation of converged maximum-payload trajectory data, and the optimization of numerous other performance indices for which no transversality conditions exist are included. The ability to simulate constrained optimum solutions, including trajectories having specified propulsion time and constant thrust cone angle, is also in the program. The program was designed to handle multiple-target missions with various types of encounters, such as rendezvous, stopover, orbital capture, and flyby. Performance requirements for a variety of launch vehicles can be determined.

Mann, F. I.↗

NASA NTRS1992

Membrane triangles with corner drilling freedoms. III - Implementation and performance evaluation

This paper completes a three-part series on the formulation of 3-node, 9-dof membrane triangles with corner drilling freedoms based on parametrized variational principles. The first four sections cover element implementation details including determination of optimal parameters and treatment of distributed loads. Then three elements of this type, labeled ALL, FF and EFF-ANDES, are tested on standard plane stress problems. ALL represents numerically integrated versions of Allman's 1988 triangle; FF is based on the free formulation triangle presented by Bergan and Felippa in 1985; and EFF-ANDES represent two different formulations of the optimal triangle derived in Parts I and II. The numerical studies indicate that the ALL, FF and EFF-ANDES elements are comparable in accuracy for elements of unitary aspect ratios. The ALL elements are found to stiffen rapidly in inplane bending for high aspect ratios, whereas the FF and EFF elements maintain accuracy. The EFF and ANDES implementations have a moderate edge in formation speed over the FF.

Felippa, Carlos A.↗

NASA NTRS2004

Higher-Order Finite Elements for Computing Thermal Radiation

Two variants of the finite-element method have been developed for use in computational simulations of radiative transfers of heat among diffuse gray surfaces. Both variants involve the use of higher-order finite elements, across which temperatures and radiative quantities are assumed to vary according to certain approximations. In this and other applications, higher-order finite elements are used to increase (relative to classical finite elements, which are assumed to be isothermal) the accuracies of final numerical results without having to refine computational meshes excessively and thereby incur excessive computation times. One of the variants is termed the radiation sub-element (RSE) method, which, itself, is subject to a number of variations. This is the simplest and most straightforward approach to representation of spatially variable surface radiation. Any computer code that, heretofore, could model surface-to-surface radiation can incorporate the RSE method without major modifications. In the basic form of the RSE method, each finite element selected for use in computing radiative heat transfer is considered to be a parent element and is divided into sub-elements for the purpose of solving the surface-to-surface radiation-exchange problem. The sub-elements are then treated as classical finite elements; that is, they are assumed to be isothermal, and their view factors and absorbed heat fluxes are calculated accordingly. The heat fluxes absorbed by the sub-elements are then transferred back to the parent element to obtain a radiative heat flux that varies spatially across the parent element. Variants of the RSE method involve the use of polynomials to interpolate and/or extrapolate to approximate spatial variations of physical quantities. The other variant of the finite-element method is termed the integration method (IM). Unlike in the RSE methods, the parent finite elements are not subdivided into smaller elements, and neither isothermality nor other unrealistic physical conditions are assumed. Instead, the equations of radiative heat transfer are integrated numerically over the parent finite elements by use of a computationally efficient Gaussian integration scheme.

Gould, Dana C.↗

NASA NTRS1985

Steady inviscid three-dimensional flows

The present analysis combines some of the theoretical concepts suggested by Hawthorne (1955) with a numerical integration procedure suggested by Martin (1978). The resulting algorithm is for inviscid subsonic flows. Thus, it is restricted to high Reynolds number flows. Chang and Adamczyk (1983) have provided a detailed derivation of the present algorithm along with a discussion of its stability bounds. The present paper represents a summary of this work. The integration of the continuity equation is considered along with an evaluation of the entropy, total temperature, and vorticity field. Attention is given to the shear-flow algorithm construction, and an application to a shear flow in a turning channel. A description of numerical results is also provided. The discussed algorithm represents a new procedure for solving inviscid subsonic three-dimensional rotational flows.

Adamczyk, J. J.↗

NASA NTRS1993

Improved Algorithm Computes Incompressible Flow

Algorithm numerically integrates Navier-Stokes equations of time-dependent or steady flow of incompressible, viscous fluid. Simulates realistic three-dimensional flows bounded by stationary or moving surface(s) of complicated shape(s). Only two arbitrary parameters specified by user. Algorithm described in "Numerical Simulation Of Flow Through An Artificial Heart" (ARC-12478). Algorithm tested by application to number of steady and unsteady flows, including flow in artificial heart discussed in noted prior article.

Rogers, Stuart E.↗

NASA NTRS2001

Human Exploration Missions Study: Space Surveillance Telescope Transfer to and Station at a Halo Orbit at the Earth-Sun Libration Point L2

This study was undertaken to determine mission profile and delta velocity requirements to place a telescope at the Earth-Sun libration point L2. The program, Integrated Mission Program (IMP), was selected to be used in the investigation. A description of IMP and its capabilities may be found in the Addenda. The Addenda also contains the libration halo equations, constants and other parameters. Comments regarding the chaotic nature of numerical integration near the libration points are also attached in the Addenda. A basic two stage S/C with a simple mission profile was selected. This profile is shown.

Dauro, Vincent A., Sr.↗