Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “approximation algorithm”

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 1,009 records · Page 56

Rapid near-optimal aerospace plane trajectory generation and guidance

Problems associated with onboard trajectory optimization, propulsion system cycle selection, and the synthesis of guidance laws are addressed for ascent to low earth orbit of an airbreathing, single-stage-to-orbit vehicle. A multicycle propulsion system is assumed that incorporates turbojet, ramjet, scramjet, and rocket engines. An energy state approximation is applied to a singularly perturbed, four-state dynamic model for flight of a point mass over a spherical nonrotating earth. An algorithm is then derived for generating both the fuel-optimal climb profile and the guidance commands required to follow that profile. In particular, analytic switching conditions are derived that, under appropriate assumptions, efficiently govern optimal transition from one propulsion cycle to another. The algorithm proves to be computationally efficient and suitable for real-time implementation. The paper concludes with the presentation of representative numerical results that illustrate the nature of the fuel-optimal climb paths and the tracking performance of the guidance algorithm.

Corban, J. E.↗

Miniature Optical Wide-Angle-Lens Startracker (Mini-OWLS)

This paper provides a brief overview of the design considerations and the current status of the Miniature Optical Wide-Angle Lens Startracker Program. Mini-OWLS offers a revolutionary alternative to the conventional startracker. It is a small, lightweight, low cost, high performance startracker that can be used in a variety of applications including calibration and alignment of Inertial Measurement Units (IMU's) Mini-OWLS makes use of a strap down design incorporating Holographic Optical Elements (HOES) in place of conventional optics. HOES can be multiplexed so that the same aperture can be used for multiple separate optical paths looking in several directions simultaneously without startracker rotation. Additionally, separate Schmidt corrector plates are not required to compensate for spherical aberration. The optical assembly, or what would normally be considered as the telescope, is less than 20 cc in volume, weighs less than 55 grams, and contains the equivalent of three individual telescopes. Each one has a 4 deg Field of View (FOV) with a field of regard of 48 square degrees. Mini-OWLS has a bandwidth of approximately 300 nm in or near the visible wavelength. The projected resolution of the startracker is 5 to 10 arcseconds, depending on the centroiding algorithm used. The Mini-OWLS program was initiated last year and represents a miniaturized version of a similar design for aeronautical applications. The contract is managed by Wright Laboratory, Air Force Systems Command, Wright-Patterson AFB, Ohio, with funding from the Strategic Defense Initiative Organization through Eglin AFB. The initial phase of the program is to build and test a development unit. The second phase is to integrate the startracker with the Charles Stark Draper Laboratory Micromechanical Inertial Guidance System (MIGS) and the Signal Processing Packaging Design (SPPD) being developed by Texas Instruments. The preliminary design review was conducted in November 1991. Three-axes prototype telescope assemblies have been built and design evaluation tests initiated.

Miller, Rick↗

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.↗

Aeroheating Predictions for X-34 Using an Inviscid-Boundary Layer Method

Radiative equilibrium surface temperatures and surface heating rates from a combined inviscid-boundary layer method are presented for the X-34 Reusable Launch Vehicle for several points along the hypersonic descent portion of its trajectory. Inviscid, perfect-gas solutions are generated with the Langley Aerothermodynamic Upwind Relaxation Algorithm (LAURA) and the Data-Parallel Lower-Upper Relaxation (DPLUR) code. Surface temperatures and heating rates are then computed using the Langley Approximate Three-Dimensional Convective Heating (LATCH) engineering code employing both laminar and turbulent flow models. The combined inviscid-boundary layer method provides accurate predictions of surface temperatures over most of the vehicle and requires much less computational effort than a Navier-Stokes code. This enables the generation of a more thorough aerothermal database which is necessary to design the thermal protection system and specify the vehicle's flight limits.

Riley, Christropher J.↗

Aeroheating Predictions for X-34 Using an Inviscid-Boundary Layer Method

Radiative equilibrium surface temperatures and surface heating rates from a combined inviscid-boundary layer method are presented for the X-34 Reusable Launch Vehicle for several points along the hypersonic descent portion of its trajectory. Inviscid, perfect-gas solutions are generated with the Langley Aerothermodynamic Upwind Relaxation Algorithm (LAURA) and the Data-Parallel Lower-Upper Relaxation (DPLUR) code. Surface temperatures and heating rates are then computed using the Langley Approximate Three-Dimensional Convective Heating (LATCH) engineering code employing both laminar and turbulent flow models. The combined inviscid-boundary layer method provides accurate predictions of surface temperatures over most of the vehicle and requires much less computational effort than a Navier-Stokes code. This enables the generation of a more thorough aerothermal database which is necessary to design the thermal protection system and specify the vehicle's flight limits.

Riley, Christopher J.↗

Event Pileup in AXAF's ACIS CCD Camera

AXAF's high resolution mirrors will focus a point source near the optical axis to a spot that is contained within a radius of about two pixels on the ACIS Charge Coupled Devices (CCD) camera. Because of the small spot size, the accuracy to which fluxes and spectral energy distributions of bright point sources can be measured will be degrad3ed by event pileup. Event pileup occurs when two or more X-ray photons arrive simultaneously in a single detection cell on a CCD readout frame. When pileup occurs, ACIS's event detection algorithm registers the photons as a single X-ray event. The pulse height channel of the event will correspond to an energy E approximately E-1 + E-2...E-n, where n is the number of photons registered per detection cell per readout frame. As a result, pileup artificially hardens the observed spectral energy distribution. I will discuss the effort at the AXAF Science Center Lo calibrate pileup in ACIS using focused, nearly monochromatic X-ray source. I will discuss techniques for modeling and correcting pileup effects in polychromatic spectra.

McNamara, Brian R.↗

Fast Fuzzy Arithmetic Operations

In engineering applications of fuzzy logic, the main goal is not to simulate the way the experts really think, but to come up with a good engineering solution that would (ideally) be better than the expert's control, In such applications, it makes perfect sense to restrict ourselves to simplified approximate expressions for membership functions. If we need to perform arithmetic operations with the resulting fuzzy numbers, then we can use simple and fast algorithms that are known for operations with simple membership functions. In other applications, especially the ones that are related to humanities, simulating experts is one of the main goals. In such applications, we must use membership functions that capture every nuance of the expert's opinion; these functions are therefore complicated, and fuzzy arithmetic operations with the corresponding fuzzy numbers become a computational problem. In this paper, we design a new algorithm for performing such operations. This algorithm is applicable in the case when negative logarithms - log(u(x)) of membership functions u(x) are convex, and reduces computation time from O(n(exp 2))to O(n log(n)) (where n is the number of points x at which we know the membership functions u(x)).

Hampton, Michael↗

A Climate-Data Record of the "Clear-Sky" Surface Temperature of the Greenland Ice Sheet

We are developing a climate-data record (CDR of daily "clear-sky" ice-surface temperature (IST) of the Greenland Ice Sheet, from 1982 to the present using Advanced Very High Resolution Radiometer (AVHRR) (1982 - present) and Moderate-Resolution Imaging Spectroradiometer (MODIS) data (2000 - present) at a resolution of approximately 5 km. The CDR will be continued in the National Polar-orbiting Operational Environmental Satellite System Visible/Infrared Imager Radiometer Suite era. Two algorithms remain under consideration. One algorithm under consideration is based on the split-window technique used in the Polar Pathfinder dataset (Fowler et al., 2000 & 21007). Another algorithm under consideration, developed by Comiso (2006), uses a single channel of AVHRR data (channel 4) in conjunction with meteorological-station data to account for atmospheric effects and drift between AVHRR instruments. Known issues being addressed in the production of the CDR are: tune-series bias caused by cloud cover (surface temperatures can be different under clouds vs. clear areas) and cross-calibration in the overlap period between AVHRR instruments, and between AVHRR and MODIS instruments. Because of uncertainties, mainly due to clouds (Stroeve & Steffen, 1998; Wang and Key, 2005; Hall et al., 2008 and Koenig and Hall, submitted), time-series of satellite 1S'1" do not necessarily correspond to actual surface temperatures. The CDR will be validated by comparing results with automatic-,",eather station (AWS) data and with satellite-derived surface-temperature products. Regional "clear-sky" surface temperature increases in the Arctic, measured from AVHRR infrared data, range from 0.57+/-0.02 deg C (Wang and Key, 2005) to 0.72+/-0.10 deg C (Comiso, 2006) per decade since the early 1980s. Arctic warming has important implications for ice-sheet mass balance because much of the periphery of the Greenland Ice Sheet is already near 0 deg C during the melt season, and is thus vulnerable to rapid melting if temperatures continue to increase. References

Hall, D. K.↗

Computing the effective elasticity of anisotropic porous media from X-ray computed micro-tomography images

Development and optimization of composite materials designed for thermal protection of NASA’s spacecraft requires the understanding of their physical response to high-enthalpy environments. To predict their macro-scale properties and behavior, high-fidelity 3D simulations are performed at the micro-scale on realistic representations of these composites. The digital micro-structures are generated either synthetically or through X-ray micro-tomography reconstructions. One of the main challenges in the prediction of heatshield material structural response is the computation of the effective elasticity of the fibrous composite, as well as the understanding of the deformation and stresses generated at the micro-scale. These are driven by the fiber layout within the micro-structure and the distribution of the infused matrix. In this effort, the micro-mechanical linear elastic behavior of fibrous ablators is modeled through the use of a numerical method based on the Multi-Point Stress Approximation (MPSA) finite volume scheme. The MPSA, a generalization of the more commonly used Multi-Point Flux Approximation (MPFA), was discussed in a previous presentation at the 15th USCCN and in reference. To predict the behavior of fibrous and woven architectures, algorithms that compute the local fiber orientation are used. The implementation of the MPSA was verified using analytical solutions, engineering test cases and compared against legacy Finite Element Analysis (FEA) software. The stress analysis models were then applied to real geometries used by NASA in thermal protection systems such as fibrous preforms and woven materials and the results were compared to experimental data.

Federico Semeraro↗

Path Planning: Differential Dynamic Programming and Model Predictive Path Integral Control on VTOL Aircraft

This paper explores two optimal control approaches, widely used in robotics, to establish their viability as real-time trajectory planners for vehicle configurations envisioned for the emerging aviation sector of Urban Air Mobility (UAM). Differential Dynamic Programming (DDP) enables planning over highly nonlinear dynamics using second-order approximations along a nominal trajectory, and displays quadratic convergence to a local solution. Model Predictive Path Integral (MPPI) is a stochastic sampling-based algorithm that can optimize for general cost criteria, including potentially highly nonlinear formulations, and supports parallel computation through the use of modern GPU hardware. In this work, DDP and MPPI were implemented using model predictive control (MPC), and the results indicate they are able to successfully transition the aircraft over different flight envelopes and generate trajectories unique to UAM vehicles.

Differential Dynamic Programming↗

Implementation of Pin Point Landing Vision Components in an FPGA System

Pin-point landing is required to enable missions to land close, typically within 10 meters, to scientifically important targets in generally hazardous terrain. In Pin Point Landing both high accuracy and high speed estimation of position and orientation is needed to provide input to the control system to safely choose and navigate to a safe landing site. A proposed algorithm called VISion aided Inertial NAVigation (VISINAV) has shown that the accuracy requirements can be met. [2][3] VISINAV was shown in software only, and was expected to use FPGA enhancements in the future to improve the computational speed needed for pin point landing during Entry Descent and Landing (EDL). Homography, feature detection and spatial correlation are computationally intensive parts of VISINAV. Homography aligns the map image with the descent image so that small correlation windows can be used, and feature detection provides regions that spatial correlation can track from frame to frame in order to estimate vehicle motion. On MER the image Homography, Feature Detection and Correlation would take approximately 650ms tracking 75 features between frames. We implemented Homography, Feature detection and Correlation on a Virtex 4 LX160 FPGA to run in under 25ms while tracking 500 features to improve algorithm reliability and throughput.

aircraft landing guidance↗

Physically-Derived Dynamical Cores in Atmospheric General Circulation Models

The algorithm chosen to represent the advection in atmospheric models is often used as the primary attribute to classify the model. Meteorological models are generally classified as spectral or grid point, with the term grid point implying discretization using finite differences. These traditional approaches have a number of shortcomings that render them non-physical. That is, they provide approximate solutions to the conservation equations that do not obey the fundamental laws of physics. The most commonly discussed shortcomings are overshoots and undershoots which manifest themselves most overtly in the constituent continuity equation. For this reason many climate models have special algorithms to model water vapor advection. This talk focuses on the development of an atmospheric general circulation model which uses a consistent physically-based advection algorithm in all aspects of the model formulation. The shallow-water model is generalized to three dimensions and combined with the physics parameterizations of NCAR's Community Climate Model. The scientific motivation for the development is to increase the integrity of the underlying fluid dynamics so that the physics terms can be more effectively isolated, examined, and improved. The expected benefits of the new model are discussed and results from the initial integrations will be presented.

Rood, Richard B.↗

Physically-Derived Dynamical Cores in Atmospheric General Circulation Models

The algorithm chosen to represent the advection in atmospheric models is often used as the primary attribute to classify the model. Meteorological models are generally classified as spectral or grid point, with the term grid point implying discretization using finite differences. These traditional approaches have a number of shortcomings that render them non-physical. That is, they provide approximate solutions to the conservation equations that do not obey the fundamental laws of physics. The most commonly discussed shortcomings are overshoots and undershoots which manifest themselves most overtly in the constituent continuity equation. For this reason many climate models have special algorithms to model water vapor advection. This talk focuses on the development of an atmospheric general circulation model which uses a consistent physically-based advection algorithm in all aspects of the model formulation. The shallow-water model of Lin and Rood (QJRMS, 1997) is generalized to three dimensions and combined with the physics parameterizations of NCAR's Community Climate Model. The scientific motivation for the development is to increase the integrity of the underlying fluid dynamics so that the physics terms can be more effectively isolated, examined, and improved. The expected benefits of the new model are discussed and results from the initial integrations will be presented.

Rood, Richard B.↗

Analysis of long term heart rate variability: methods, 1/f scaling and implications

The use of spectral techniques to quantify short term heart rate fluctuations on the order of seconds to minutes has helped define the autonomic contributions to beat-to-beat control of heart rate. We used similar techniques to quantify the entire spectrum (0.00003-1.0 Hz) of heart rate variability during 24 hour ambulatory ECG monitoring. The ECG from standard Holter monitor recordings from normal subjects was sampled with the use of a phase locked loop, and a heart rate time series was constructed at 3 Hz. Frequency analysis of the heart rate signal was performed after a nonlinear filtering algorithm was used to eliminate artifacts. A power spectrum of the entire 24 hour record revealed power that was inversely proportional to frequency, 1/f, over 4 decades from 0.00003 to 0.1 Hz (period approximately 10 hours to 10 seconds). Displaying consecutive spectra calculated at 5 minute intervals revealed marked variability in the peaks at all frequencies throughout the 24 hours, probably accounting for the lack of distinct peaks in the spectra of the entire records.

Non-NASA Center↗

Design and Optimization of Low-thrust Orbit Transfers Using Q-law and Evolutionary Algorithms

Future space missions will depend more on low-thrust propulsion (such as ion engines) thanks to its high specific impulse. Yet, the design of low-thrust trajectories is complex and challenging. Third-body perturbations often dominate the thrust, and a significant change to the orbit requires a long duration of thrust. In order to guide the early design phases, we have developed an efficient and efficacious method to obtain approximate propellant and flight-time requirements (i.e., the Pareto front) for orbit transfers. A search for the Pareto-optimal trajectories is done in two levels: optimal thrust angles and locations are determined by Q-law, while the Q-law is optimized with two evolutionary algorithms: a genetic algorithm and a simulated-annealing-related algorithm. The examples considered are several types of orbit transfers around the Earth and the asteroid Vesta.

low-thrust orbit transfer↗

Adaptive Control of Small Outboard-Powered Boats for Survey Applications

Four autopilot controllers have been developed in this work that can both hold a desired heading and follow a straight line. These PID, adaptive PID, neuro-adaptive, and adaptive augmenting control algorithms have all been implemented into a numerical simulation of a 33-foot center console vessel with wind, waves, and current disturbances acting in the perpendicular (across-track) direction of the boat s desired trajectory. Each controller is tested for its ability to follow a desired heading in the presence of these disturbances and then to follow a straight line at two different throttle settings for the same disturbances. These controllers were tuned for an input thrust of 2000 N and all four controllers showed good performance with none of the controllers significantly outperforming the others when holding a constant heading and following a straight line at this engine thrust. Each controller was then tested for a reduced engine thrust of 1200 N per engine where each of the three adaptive controllers reduced heading error and across-track error by approximately 50% after a 300 second tuning period when compared to the fixed gain PID, showing that significant robustness to changes in throttle setting was gained by using an adaptive algorithm.

VanZwieten, T.S.↗

Sparse matrix techniques applied to modal analysis of multi-section duct liners

A simplified procedure is presented for analysis of ducts with discretely nonuniform properties. The analysis uses basis functions as the generalized coordinates. The duct eigenfunctions are approximated by finite series of these functions. The emphasis is on solution of the resulting large sparse set of linear equations. Characteristics of sparse matrix algorithms are outlined and some criteria for application are established. Analogies with structural methods are used to illustrate variations which can increase efficiency in generating values for design optimization routines. The effects of basis function selection, number of eigenfunctions and identification and ordering of equations on the sparsity and solution stability are included.

Arnold, W. R.↗

Adaptive Bayes classifiers for remotely sensed data

An algorithm is developed for a learning, adaptive, statistical pattern classifier for remotely sensed data. The estimation procedure consists of two steps: (1) an optimal stochastic approximation of the parameters of interest, and (2) a projection of the parameters in time and space. The results reported are for Gaussian data in which the mean vector of each class may vary with time or position after the classifier is trained.

Raulston, H. S.↗