SEARCH · Engineering Papers

Results for “Incremental Computing”

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 55 records · Page 3

All sources NASA NTRS DOE OSTI

NASA NTRS1983

Design, implementation and flight testing of PIF autopilots for general aviation aircraft

The designs of Proportional-Integrated-Filter (PIF) auto-pilots for a General Aviation (NAVION) aircraft are presented. The PIF autopilot uses the sampled-data regulator and command generator tracking to determine roll select, pitch select, heading select, altitude select and localizer/glideslope capture and hold autopilot modes. The PIF control law uses typical General Aviation sensors for state feedback, command error integration for command tracking, digital complementary filtering and analog prefiltering for sensor noise suppression, a control filter for computation delay accommodation and the incremental form to eliminate trim values in implementation. Theoretical developments described in detail, were needed to combine the sampled-data regulator with command generator tracking for use as a digital flight control system. The digital PIF autopilots are evaluated using closed-loop eigenvalues and linear simulations. The implementation of the PIF autopilots in a digital flight computer using a high order language (FORTRAN) is briefly described. The successful flight test results for each PIF autopilot mode is presented.

Broussard, J. R.↗

NASA NTRS1982

Design and flight testing of a digital optimal control general aviation autopilot

This paper presents the designs of Proportional-Integral-Filter (PIF) autopilots for a General Aviation (NAVION) aircraft. The PIF autopilots use modern control theory to determine heading select and altitude select and hold autopilot modes. The PIF control law uses typical General Aviation sensors for state feedback; command error integration for command tracking; digital complimentary filtering and analog prefiltering for sensor noise suppression; a control filter for computation delay accommodation; and the incremental form to eliminate trim values in implementation. Theoretical developments for the control law are described which combine the sampled-data regulator with command generator tracking for use as a digital flight control system. The digital PIF autopilots are evaluated using closed-loop eigenvalues and simulations. Successful flight test results for the PIF autopilots are presented for different turbulence conditions and quadratic weights.

Broussard, J. R.↗

NASA NTRS1999

A Computational/Experimental Study of Two Optimized Supersonic Transport Designs and the Reference H Baseline

Two supersonic transport configurations designed by use of non-linear aerodynamic optimization methods are compared with a linearly designed baseline configuration. One optimized configuration, designated Ames 7-04, was designed at NASA Ames Research Center using an Euler flow solver, and the other, designated Boeing W27, was designed at Boeing using a full-potential method. The two optimized configurations and the baseline were tested in the NASA Langley Unitary Plan Supersonic Wind Tunnel to evaluate the non-linear design optimization methodologies. In addition, the experimental results are compared with computational predictions for each of the three configurations from the Enter flow solver, AIRPLANE. The computational and experimental results both indicate moderate to substantial performance gains for the optimized configurations over the baseline configuration. The computed performance changes with and without diverters and nacelles were in excellent agreement with experiment for all three models. Comparisons of the computational and experimental cruise drag increments for the optimized configurations relative to the baseline show excellent agreement for the model designed by the Euler method, but poorer comparisons were found for the configuration designed by the full-potential code.

Cliff, Susan E.↗

NASA NTRS2001

Temperature Humidity and Sea Level Pressure Increments Induced by 1DVAR Analysis of GPS Refractivity

The Global Positioning System (GPS) transmitted signals are affected by the atmosphere. Using the radio occultation technique, where a receiver is placed on a low-Earth-orbiting platform. it is possible to perform soundings. by measuring the bending angles of the rays. The information can be converted into atmospheric refractivity. We have developed a one dimensional variational (1DVAR) analysis that uses GPS/MET 1995 refractivity and 6-hour FVDAS (Finite Volume Data Assimilation System) forecasts as background information to constrain the retrievals. The analysis increments are defined as 1DVAR minus background temperature, humidity and sea level pressure. Before assimilating the 1DVAR profiles into the FVDAS. the increments need to be understood. First, some bias could be induced in the retrievals when confronted with actual biased data: second. bias in the back-round could create undesired bias in the retrievals. Anv bias in the analyses will ultimately change the climatology of the model the retrievals will be assimilated into. We relate the increments to the reduction of the difference between observed minus computed refractivity profiles. We also point out the difference in the mean increments using backgrounds which have assimilated either NESDIS TIROS Operational Vertical Sounder (TOVS) operational retrievals or Data Assimilation Office (DAO) TOVS interactive retrievals. The climatology of the model in terms of refractivity is significantly different and this impacts the GPS 1DVAR increments. This testifies that changing the basic load of assimilated data has an influence on the impact the GPS data may have in a DAS.

Poli, Paul↗

NASA NTRS1974

Analog computer implementation of four instantaneous distortion indices

Dynamic distortion data obtained from inlet engine compatibility test on a J85-GE-13 engine are used to establish the feasibility of using on-line generated indices as control signals. These data are analyzed over time increments of 1.3 seconds on an analog computer. The analog program of four instantaneous distortion indices and their results are presented. A modified approach in determining the extent of distortion is also presented.

Costakis, W. G.↗

NASA NTRS1972

VISCEL: A general-purpose computer program for analysis of linear viscoelastic structures (user's manual), volume 1

This program, an extension of the linear equilibrium problem solver ELAS, is an updated and extended version of its earlier form (written in FORTRAN 2 for the IBM 7094 computer). A synchronized material property concept utilizing incremental time steps and the finite element matrix displacement approach has been adopted for the current analysis. A special option enables employment of constant time steps in the logarithmic scale, thereby reducing computational efforts resulting from accumulative material memory effects. A wide variety of structures with elastic or viscoelastic material properties can be analyzed by VISCEL. The program is written in FORTRAN 5 language for the Univac 1108 computer operating under the EXEC 8 system. Dynamic storage allocation is automatically effected by the program, and the user may request up to 195K core memory in a 260K Univac 1108/EXEC 8 machine. The physical program VISCEL, consisting of about 7200 instructions, has four distinct links (segments), and the compiled program occupies a maximum of about 11700 words decimal of core storage.

Gupta, K. K.↗

NASA NTRS1974

User's guide for analysis of finite elastoplastic deformation: The FIPDEF and FIPAX programs for the CDC 6600

Computer programs are presented which provide incremental finite-element analysis capability for problems of quasi-static, finite, elastoplastic deformation in two spatial dimensions (plane strain, plane stress, axisymmetric). Monotonic or cyclic loading of isotropic hardening materials is considered. The only restriction on the form of the stress-strain curve is that the rate of work hardening exceed some small positive value. The user's guide assumes familiarity with both finite-element analysis and FORTRAN IV programming for the CDC 6600. Sufficient information is provided to support problem solving ultization of the programs.

Osias, J. R.↗

NASA NTRS1975

The influence of crack closure and elasto-plastic flow on the bending of a cracked plate

The influence of crack closure and elastoplastic flow on the bending behavior of thin cracked plates is investigated using an incremental elastoplastic plate bending finite element computer program. The finite element program was developed using assumptions consistent with Kirchhoff fourth-order plate theory while the material property treatment permits general isotropic work hardening with local elastic unloading. This technique is applied to the problem of a large centrally through cracked plate subject to remote circular bending. Comparison is drawn between two cases of the bending problem. The first neglects the possibility of crack face interference with bending, and the second includes a kinematic prohibition against the crack face from passing through the symmetry plane. Results are reported which isolate the effects of elastoplastic flow and crack closure.

Jones, D. P.↗

NASA NTRS1975

Circular earth orbits attainable with fixed two-impulse expendable tug

The performance of solid rocket motors in expendable stages as an alternative concept for Shuttle interim upper stage applications is assessed from the standpoint of mission flexibility. The system comprises two stages and a payload, released from a circular Shuttle orbit at 160 naut miles altitude. The stages are sized to deliver a given payload mass to synchronous equatorial orbit. Only nonreturn missions are considered. A parametric sweep was made with small increments of pitch and yaw to compute the possible transfer orbits with an apogee greater than 160 naut miles. At several points on each transfer orbit the possibility of attaining a circular orbit with the second burn impulse of 6000 fps was assessed. Constraints on the vertical component of the second impulse and on the total velocity after firing determine the pitch angle for second burn. Results indicate a fixed two-impulse system designed to place a payload into synchronous equatorial orbit from the Shuttle orbit can deliver the payload mass into any lower circular orbit at any inclination within about 30 deg of the Shuttle orbit. The only apparent disadvantage of the system is an inability to achieve circular orbits above geosynchronous altitudes with the design payload.

Porter, R. F.↗

NASA NTRS1982

An experimental investigation of hypersonic flow over biconics at incidence and comparison to prediction

Pressure distributions, aerodynamic coefficients, shock shapes, and oil-flow patterns were measured on spherically blunted, 12.84 deg/7 deg on-axis and bent biconics at Mach 6 in air. The angle of attack, referenced to the axis of the aft-cone, was varied from -10 to 40 deg in nominal 5 deg increments. Predictions from an inviscid flowfield computer code referred to as STEIN and codes which solve the three-dimensional parabolized Navier Stokes (PNS) equations were compared with measurement. Three PNS codes were found which provided accurate predictions of windward and leeward pressure distributions for angles of attack to 25 deg.

Miller, C. G.↗

NASA NTRS1987

The matrix exponential in transient structural analysis

The primary usefulness of the presented theory is in the ability to represent the effects of high frequency linear response with accuracy, without requiring very small time steps in the analysis of dynamic response. The matrix exponential contains a series approximation to the dynamic model. However, unlike the usual analysis procedure which truncates the high frequency response, the approximation in the exponential matrix solution is in the time domain. By truncating the series solution to the matrix exponential short, the solution is made inaccurate after a certain time. Yet, up to that time the solution is extremely accurate, including all high frequency effects. By taking finite time increments, the exponential matrix solution can compute the response very accurately. Use of the exponential matrix in structural dynamics is demonstrated by simulating the free vibration response of multi degree of freedom models of cantilever beams.

Minnetyan, Levon↗

NASA NTRS1980

Three dimensional inelastic finite element analysis of laminated composites

Formulations of the inelastic response of laminated composites to thermal and mechanical loading are used as the basis for development of the computer NALCOM (Nonlinear Analysis of Laminated Composites) computer program which uses a fully three dimensional isoparametric finite element with 24 nodes and 72 degrees of freedom. An incremental solution is performed with nonlinearities introduced as pseudoloads computed for initial strains. Equilibrium iteration may be performed at every step. Elastic and elastic-plastic response of boron/epoxy and graphite/epoxy graphite/epoxy and problems of curing 0/90 sub s Gr/Ep laminates with and without circular holes are analyzed. Mechanical loading of + or - 45sub s Gr/Ep laminates is modeled and symmetry conditions which exist in angle-ply laminates are discussed. Results are compared to experiments and other analytical models when possible. All models are seen to agree reasonably well with experimetnal results for off-axis tensile coupons. The laminate analyses show the three dimensional effects which are present near holes and free corners.

Griffin, O. H., Jr.↗

NASA NTRS1998

Statistical Properties of Line Centroid Velocity Increments in the rho Ophiuchi Cloud

We present a comparison of histograms of CO (2-1) line centroid velocity increments in the rho Ophiuchi molecular cloud with those computed for spectra synthesized from a three-dimensional, compressible, but non-starforming and non-gravitating hydrodynamic simulation. Histograms of centroid velocity increments in the rho Ophiuchi cloud show clearly non-Gaussian wings, similar to those found in histograms of velocity increments and derivatives in experimental studies of laboratory and atmospheric flows, as well as numerical simulations of turbulence. The magnitude of these wings increases monotonically with decreasing separation, down to the angular resolution of the data. This behavior is consistent with that found in the phase of the simulation which has most of the properties of incompressible turbulence. The time evolution of the magnitude of the non-Gaussian wings in the histograms of centroid velocity increments in the simulation is consistent with the evolution of the vorticity in the flow. However, we cannot exclude the possibility that the wings are associated with the shock interaction regions. Moreover, in an active starforming region like the rho Ophiuchi cloud, the effects of shocks may be more important than in the simulation. However, being able to identify shock interaction regions in the interstellar medium is also important, since numerical simulations show that vorticity is generated in shock interactions.

Lis, D. C.↗

NASA NTRS1996

The incremental unknowns-a multilevel scheme for the simulation of turbulent channel flows

In numerical simulation of complex flows, it is important to identify different length scales of the flow and treat them differently. In this report, we introduce a new multilevel scheme for simulating turbulent channel flows. Two different versions of the scheme, namely the spectral and finite difference versions, are presented. The spectral version of the scheme is based on a spectral-Galerkin formulation which provides a natural decomposition of the flow into small and large wavelength parts, and which leads to linear systems that can be solved with quasi-optimal computational complexity. In the finite difference version, the Incremental Unknown (IU) is used to separate the length scales. Preliminary numerical results indicate that the scheme is well suited for turbulence computations and provides results which are comparable to that by Direct Numerical Simulation (DNS) but with significantly less CPU time.

Chen, M.↗

NASA NTRS2019

Investigation of Reduced-Order Modeling for Aircraft Stability and Control Prediction

High-fidelity computational fluid dynamics tools offer the potential to approximate increments for ground-to-flight scaling effects, as well as to augment the dynamic damping derivative data for motion-based flight simulators. Unfortunately, the computational expense is currently prohibitive for populating a complete simulator database. This work investigates an existing surrogate-based, indicial response reduced-order model methodology as a means to efficiently augment a flight simulator database with high-fidelity nonlinear aerodynamic damping derivatives. Creation of the reduced-order model is based on the superposition integrals of the step response with the derivative of its corresponding input signal. Step responses are calculated using a computational grid motion approach that separates the effects of angle of attack and sideslip angle from angular rates, and rates from angle of attack and sideslip. It is demonstrated that the transients produced during the start of a forced-oscillation motion are captured by the reduced-order model to the level of fidelity of a comparable computational solution. Aerodynamic coefficients computed within minutes by the reduced-order model for an aircraft undergoing an 18-second half Lazy-8 maneuver and a 25-second Immelmann turn maneuver are compared with those from full computational flight solutions that required days to complete. Finally, a cost-benefit assessment is included that demonstrates a compelling advantage for this approach. d for maneuvering, flexible vehicles.

Frink, Neal T.↗

NASA NTRS1993

Decision theory for computing variable and value ordering decisions for scheduling problems

Heuristics that guide search are critical when solving large planning and scheduling problems, but most variable and value ordering heuristics are sensitive to only one feature of the search state. One wants to combine evidence from all features of the search state into a subjective probability that a value choice is best, but there has been no solid semantics for merging evidence when it is conceived in these terms. Instead, variable and value ordering decisions should be viewed as problems in decision theory. This led to two key insights: (1) The fundamental concept that allows heuristic evidence to be merged is the net incremental utility that will be achieved by assigning a value to a variable. Probability distributions about net incremental utility can merge evidence from the utility function, binary constraints, resource constraints, and other problem features. The subjective probability that a value is the best choice is then derived from probability distributions about net incremental utility. (2) The methods used for rumor control in Bayesian Networks are the primary way to prevent cycling in the computation of probable net incremental utility. These insights lead to semantically justifiable ways to compute heuristic variable and value ordering decisions that merge evidence from all available features of the search state.

Linden, Theodore A.↗

NASA NTRS1996

A General Interface Method for Aeroelastic Analysis of Aircraft

The aeroelastic analysis of an aircraft requires an accurate and efficient procedure to couple aerodynamics and structures. The procedure needs an interface method to bridge the gap between the aerodynamic and structural models in order to transform loads and displacements. Such an interface method is described in this report. This interface method transforms loads computed by any aerodynamic code to a structural finite element (FE) model and converts the displacements from the FE model to the aerodynamic model. The approach is based on FE technology in which virtual work is employed to transform the aerodynamic pressures into FE nodal forces. The displacements at the FE nodes are then converted back to aerodynamic grid points on the aircraft surface through the reciprocal theorem in structural engineering. The method allows both high and crude fidelities of both models and does not require an intermediate modeling. In addition, the method performs the conversion of loads and displacements directly between individual aerodynamic grid point and its corresponding structural finite element and, hence, is very efficient for large aircraft models. This report also describes the application of this aero-structure interface method to a simple wing and an MD-90 wing. The results show that the aeroelastic effect is very important. For the simple wing, both linear and nonlinear approaches are used. In the linear approach, the deformation of the structural model is considered small, and the loads from the deformed aerodynamic model are applied to the original geometry of the structure. In the nonlinear approach, the geometry of the structure and its stiffness matrix are updated in every iteration and the increments of loads from the previous iteration are applied to the new structural geometry in order to compute the displacement increments. Additional studies to apply the aero-structure interaction procedure to more complicated geometry will be conducted in the second phase of the present contract.

Tzong, T.↗

NASA NTRS2013

Investigations of Fluid-Structure-Coupling and Turbulence Model Effects on the DLR Results of the Fifth AIAA CFD Drag Prediction Workshop

The accurate calculation of aerodynamic forces and moments is of significant importance during the design phase of an aircraft. Reynolds-averaged Navier-Stokes (RANS) based Computational Fluid Dynamics (CFD) has been strongly developed over the last two decades regarding robustness, efficiency, and capabilities for aerodynamically complex configurations. Incremental aerodynamic coefficients of different designs can be calculated with an acceptable reliability at the cruise design point of transonic aircraft for non-separated flows. But regarding absolute values as well as increments at off-design significant challenges still exist to compute aerodynamic data and the underlying flow physics with the accuracy required. In addition to drag, pitching moments are difficult to predict because small deviations of the pressure distributions, e.g. due to neglecting wing bending and twisting caused by the aerodynamic loads can result in large discrepancies compared to experimental data. Flow separations that start to develop at off-design conditions, e.g. in corner-flows, at trailing edges, or shock induced, can have a strong impact on the predictions of aerodynamic coefficients too. Based on these challenges faced by the CFD community a working group of the AIAA Applied Aerodynamics Technical Committee initiated in 2001 the CFD Drag Prediction Workshop (DPW) series resulting in five international workshops. The results of the participants and the committee are summarized in more than 120 papers. The latest, fifth workshop took place in June 2012 in conjunction with the 30th AIAA Applied Aerodynamics Conference. The results in this paper will evaluate the influence of static aeroelastic wing deformations onto pressure distributions and overall aerodynamic coefficients based on the NASA finite element structural model and the common grids.

Keye, Stefan↗