Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “algorithmic differentiation”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

Spectral element methods - Algorithms and architectures

Spectral element methods are high-order weighted residual techniques for partial differential equations that combine the geometric flexibility of finite element methods with the rapid convergence of spectral techniques. Spectral element methods are described for the simulation of incompressible fluid flows, with special emphasis on implementation of spectral element techniques on medium-grained parallel processors. Two parallel architectures are considered; the first, a commercially available message-passing hypercube system; the second, a developmental reconfigurable architecture based on Geometry-Defining Processors. High parallel efficiency is obtained in hypercube spectral element computations, indicating that load balancing and communication issues can be successfully addressed by a high-order technique/medium-grained processor algorithm-architecture coupling.

Fischer, Paul↗
NASA NTRS1978

An interaction solution algorithm for viscous aerodynamic performance prediction

A weak-interaction solution algorithm is established for aerodynamic flow field prediction about an isolated airfoil. It requires numerical solution of differential equations governing potential flow, viscous and turbulent boundary layer flow, and the turbulent wake flow downstream of a trailing edge. The algorithm accounts for computed viscous displacement effects on the potential flow. These in turn alter the viscous flow through imposed pressure gradients. Closure for turbulence is accomplished using a second order model. Numerical evaluations assess factors affecting solution accuracy, convergence and stability for the combined potential, boundary layer, and parabolic Navier-Stokes equation systems as solved using a finite element algorithm.

Baker, A. J.↗
NASA NTRS2003

Total Ozone from the Ozone Monitoring System (OMI) using TOMS and DOAS Methods

The Ozone Monitoring Instrument (OMI) is the Dutch-Finnish contribution to NASA's EOS-Aura satellite scheduled for launch in January 2004. OMI is an imaging spectrometer that will measure the back-scattered Solar radiance in the wavelength range of 270 to 500 nm. The instrument provides near global coverage in one day with a spatial resolution of 13x24 square kilometers. OMI is a new instrument, with a heritage from TOMS, SBW, GOME, GOMOS and SCIAMACHY. OMI'S unique capabilities for measuring important trace gases and aerosols with a small footprint and daily global coverage, in conjunction with the other Aura instruments, will make a major contribution to our understanding of stratospheric and tropospheric chemistry and climate change. OMI will provide data continuity with the 23-year ozone record of TOMS. There are three ozone products planned for OMI: total column ozone, ozone profile and tropospheric column ozone. We are developing two different algorithms for total column ozone: one similar to the algorithm currently being used to process the TOMS data, and the other an improved version of the differential optical absorption spectroscopy (DOAS) method, which has been applied to GOME and SCIAMACHY data. The main reasons for starting with two algorithms for total ozone have to do with heritage and past experience; our long-term goal is to combine the two to develop a more accurate and reliable total ozone product for OMI. We will compare the performance of these two algorithms by applying both of them to the GOME data. We will examine where and how the results differ, and use the extensive TOMS-Dobson comparison studies to assess the performance of the DOAS algorithm.

Veefkind, J. P.↗
NASA NTRS1983

Optimal trajectories for aeroassisted orbital transfer

Consideration is given to classical and minimax problems involved in aeroassisted transfer from high earth orbit (HEO) to low earth orbit (LEO). The transfer is restricted to coplanar operation, with trajectory control effected by means of lift modulation. The performance of the maneuver is indexed to the energy expenditure or, alternatively, the time integral of the heating rate. Firist-order optimality conditions are defined for the classical approach, as are a sequential gradient-restoration algorithm and a combined gradient-restoration algorithm. Minimization techniques are presented for the aeroassisted transfer energy consumption and time-delay integral of the heating rate, as well as minimization of the pressure. It is shown that the eigenvalues of the Jacobian matrix of the differential system is both stiff and unstable, implying that the sequential gradient restoration algorithm in its present version is unsuitable. A new method, involving a multipoint approach to the two-poing boundary value problem, is recommended.

Miele, A.↗
NASA NTRS1989

Dynamic grid refinement for partial differential equations on parallel computers

The fast adaptive composite grid method (FAC) is an algorithm that uses various levels of uniform grids to provide adaptive resolution and fast solution of PDEs. An asynchronous version of FAC, called AFAC, that completely eliminates the bottleneck to parallelism is presented. This paper describes the advantage that this algorithm has in adaptive refinement for moving singularities on multiprocessor computers. This work is applicable to the parallel solution of two- and three-dimensional shock tracking problems.

Mccormick, S.↗
NASA NTRS2001

Application of multivariate statistics to vestibular testing: discriminating between Meniere's disease and migraine associated dizziness

Meniere's disease (MD) and migraine associated dizziness (MAD) are two disorders that can have similar symptomatologies, but differ vastly in treatment. Vestibular testing is sometimes used to help differentiate between these disorders, but the inefficiency of a human interpreter analyzing a multitude of variables independently decreases its utility. Our hypothesis was that we could objectively discriminate between patients with MD and those with MAD using select variables from the vestibular test battery. Sinusoidal harmonic acceleration test variables were reduced to three vestibulo-ocular reflex physiologic parameters: gain, time constant, and asymmetry. A combination of these parameters plus a measurement of reduced vestibular response from caloric testing allowed us to achieve a joint classification rate of 91%, independent quadratic classification algorithm. Data from posturography were not useful for this type of differentiation. Overall, our classification function can be used as an unbiased assistant to discriminate between MD and MAD and gave us insight into the pathophysiologic differences between the two disorders.

NASA Discipline Neuroscience↗
NASA NTRS1990

An adaptive grid algorithm for one-dimensional nonlinear equations

Richards' equation, which models the flow of liquid through unsaturated porous media, is highly nonlinear and difficult to solve. Step gradients in the field variables require the use of fine grids and small time step sizes. The numerical instabilities caused by the nonlinearities often require the use of iterative methods such as Picard or Newton interation. These difficulties result in large CPU requirements in solving Richards equation. With this in mind, adaptive and multigrid methods are investigated for use with nonlinear equations such as Richards' equation. Attention is focused on one-dimensional transient problems. To investigate the use of multigrid and adaptive grid methods, a series of problems are studied. First, a multigrid program is developed and used to solve an ordinary differential equation, demonstrating the efficiency with which low and high frequency errors are smoothed out. The multigrid algorithm and an adaptive grid algorithm is used to solve one-dimensional transient partial differential equations, such as the diffusive and convective-diffusion equations. The performance of these programs are compared to that of the Gauss-Seidel and tridiagonal methods. The adaptive and multigrid schemes outperformed the Gauss-Seidel algorithm, but were not as fast as the tridiagonal method. The adaptive grid scheme solved the problems slightly faster than the multigrid method. To solve nonlinear problems, Picard iterations are introduced into the adaptive grid and tridiagonal methods. Burgers' equation is used as a test problem for the two algorithms. Both methods obtain solutions of comparable accuracy for similar time increments. For the Burgers' equation, the adaptive grid method finds the solution approximately three times faster than the tridiagonal method. Finally, both schemes are used to solve the water content formulation of the Richards' equation. For this problem, the adaptive grid method obtains a more accurate solution in fewer work units and less computation time than required by the tridiagonal method. The performance of the adaptive grid method tends to degrade as the solution process proceeds in time, but still remains faster than the tridiagonal scheme.

Gutierrez, William E.↗
NASA NTRS2019

Real-Time Modeling of Thermal Interactions in Cryogenic Ball Bearings

Based on the integration of classical differential equations of motion of bearing elements, a real-time dynamic simulation of ball bearing performance is coupled with simultaneous modeling of thermal interactions between the bearing elements. The transient heat generations are time averaged over a significantly larger time step before they are used to compute the changes in bearing temperature fields. Such an averaging algorithm eliminates the numerical difficulties associated with simultaneous integration of mechanical thermal differential equations with vastly different time scales. Under stable operation, the step change in temperatures at each thermal time step converges to a steady-state solution. Following is a list of accomplishments in the present investigation: (1) A real-time dynamic performance of both 440C and hybrid ball bearings in a LOX environment is modeled in both transient and steady-state time domains, as the inner race accelerates to a prescribed operating speed. (2) The frictional interactions at the ball to race contacts are modeled by experimentally measured traction coefficients as a function of slide-to-roll ratios. (3) Bearing heat generation, as defined by the heat transferred to circulating LOX, is validated against experimental values determined by measured LOX temperatures under prescribed pressures and flow rates. The predicted heat generation is in good agreement with the experimental values. (4) The steady-state solutions are shown to be independent of initial conditions. This establishes convergence of the numerical integration of a differential equation of motion of the bearing elements. (5) Churning and drag losses in the circulating LOX constitutes the majority of bearing heat generation under the experimental test conditions. (6) The heat dissipated in frictional interactions in hybrid bearings is significantly lower in hybrid bearings in comparison to that simulated in all-steel bearings, while the churning and drag losses are relatively unchanged.

Gupta, P. K.↗
NASA NTRS1974

Equations of condition for high order Runge-Kutta-Nystrom formulae

Derivation of the equations of condition of order eight for a general system of second-order differential equations approximated by the basic Runge-Kutta-Nystrom algorithm. For this general case, the number of equations of condition is considerably larger than for the special case where the first derivative is not present. Specifically, it is shown that, for orders two through eight, the number of equations for each order is 1, 1, 1, 2, 3, 5, and 9 for the special case and is 1, 1, 2, 5, 13, 34, and 95 for the general case.

Bettis, D. G.↗
NASA NTRS1980

An approximate factorization solution of the Navier-Stokes equations for transonic flow using body-fitted coordinates with application to NACA 64A010 airfoils

The implementation of the approximate factorization algorithm and its ability to efficiently and accurately describe transonic flow about an NACA 64A010 airfoil section is examined. The approximate factorization algorithm is developed from the nondimensional, conservative, vectorized Navier-Stokes equations expressed in curvilinear coordinates. Equations of state and transport coefficient relations appropriate to atmospheric air are appended to close the system of partial differential equations. An algebraic turbulence model is also incorporated into the equation set. This algorithm was verified by investigating the flow about an NACA 64A010 airfoil at 0, 2, and 3.5 deg angle of attack for free-stream conditions of 2,000,000 Reynolds number and 0.8 Mach number. Overall results were in good qualitative agreement with wind tunnel data sets. However, while nondimensional times of six were attained, numerical difficulties prevented any case from reaching a true steady state.

Copper, G. K.↗
NASA NTRS1982

Survey of new vector computers: The CRAY 1S from CRAY research; the CYBER 205 from CDC and the parallel computer from ICL - architecture and programming

Problems which can arise with vector and parallel computers are discussed in a user oriented context. Emphasis is placed on the algorithms used and the programming techniques adopted. Three recently developed supercomputers are examined and typical application examples are given in CRAY FORTRAN, CYBER 205 FORTRAN and DAP (distributed array processor) FORTRAN. The systems performance is compared. The addition of parts of two N x N arrays is considered. The influence of the architecture on the algorithms and programming language is demonstrated. Numerical analysis of magnetohydrodynamic differential equations by an explicit difference method is illustrated, showing very good results for all three systems. The prognosis for supercomputer development is assessed.

Gentzsch, W.↗
NASA NTRS1983

Implicit, nonswitching, vector-oriented algorithm for steady transonic flow

A rapid computation of a sequence of transonic flow solutions has to be performed in many areas of aerodynamic technology. The employment of low-cost vector array processors makes the conduction of such calculations economically feasible. However, for a full utilization of the new hardware, the developed algorithms must take advantage of the special characteristics of the vector array processor. The present investigation has the objective to develop an efficient algorithm for solving transonic flow problems governed by mixed partial differential equations on an array processor.

Lottati, I.↗
NASA NTRS1987

Parallel processing architecture for computing inverse differential kinematic equations of the PUMA arm

In advanced robot control problems, on-line computation of inverse Jacobian solution is frequently required. Parallel processing architecture is an effective way to reduce computation time. A parallel processing architecture is developed for the inverse Jacobian (inverse differential kinematic equation) of the PUMA arm. The proposed pipeline/parallel algorithm can be inplemented on an IC chip using systolic linear arrays. This implementation requires 27 processing cells and 25 time units. Computation time is thus significantly reduced.

Hsia, T. C.↗
NASA NTRS1992

Scheme For Finite-Difference Computations Of Waves

Compact algorithms generating and solving finite-difference approximations of partial differential equations for propagation of waves obtained by new method. Based on concept of discrete dispersion relation. Used in wave propagation to relate frequency to wavelength and is key measure of wave fidelity.

Davis, Sanford↗
NASA NTRS2003

Analysis Method for Inelastic, Adhesively Bonded Joints with Anisotropic Adherends

A one-dimensional analysis method for evaluating adhesively bonded joints composed of anisotropic adherends and adhesives with nonlinear material behavior is presented in the proposed paper. The strain and resulting stress field in a general, bonded joint overlap are determined by using a variable-step, finite-difference solution algorithm to iteratively solve a system of first-order differential equations. Applied loading is given by a system of combined extensional, bending, and shear forces that are applied to the edge of the joint overlap. Adherends are assumed to behave as linear, cylindrically bent plates using classical laminated plate theory that includes the effects of first-order transverse shear deformation. Using the deformation theory of plasticity and a modified von-Mises yield criterion, inelastic material behavior is modeled in the adhesive layer. Results for the proposed method are verified against previous results from the literature and shown to be in excellent agreement. An additional case that highlights the effects of transverse shear deformation between similar adherends is also presented.

Smeltzer, Stanley S., III↗
NASA NTRS2004

Multiobjective Aerodynamic Shape Optimization Using Pareto Differential Evolution and Generalized Response Surface Metamodels

Differential Evolution (DE) is a simple, fast, and robust evolutionary algorithm that has proven effective in determining the global optimum for several difficult single-objective optimization problems. The DE algorithm has been recently extended to multiobjective optimization problem by using a Pareto-based approach. In this paper, a Pareto DE algorithm is applied to multiobjective aerodynamic shape optimization problems that are characterized by computationally expensive objective function evaluations. To improve computational expensive the algorithm is coupled with generalized response surface meta-models based on artificial neural networks. Results are presented for some test optimization problems from the literature to demonstrate the capabilities of the method.

Madavan, Nateri K.↗
NASA NTRS2007

White-light Interferometry using a Channeled Spectrum: II. Calibration Methods, Numerical and Experimental Results

In the companion paper, [Appl. Opt. 46, 5853 (2007)] a highly accurate white light interference model was developed from just a few key parameters characterized in terms of various moments of the source and instrument transmission function. We develop and implement the end-to-end process of calibrating these moment parameters together with the differential dispersion of the instrument and applying them to the algorithms developed in the companion paper. The calibration procedure developed herein is based on first obtaining the standard monochromatic parameters at the pixel level: wavenumber, phase, intensity, and visibility parameters via a nonlinear least-squares procedure that exploits the structure of the model. The pixel level parameters are then combined to obtain the required 'global' moment and dispersion parameters. The process is applied to both simulated scenarios of astrometric observations and to data from the microarcsecond metrology testbed (MAM), an interferometer testbed that has played a prominent role in the development of this technology.

phase correction↗