Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “variational algorithms”

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 127 records · Page 7

Optimum structural design with plate bending elements - A survey

A survey is presented of recently published papers in the field of optimum structural design of plates, largely with respect to the minimum-weight design of plates subject to such constraints as fundamental frequency maximization. It is shown that, due to the availability of powerful computers, the trend in optimum plate design is away from methods tailored to specific geometry and loads and toward methods that can be easily programmed for any kind of plate, such as finite element methods. A corresponding shift is seen in optimization from variational techniques to numerical optimization algorithms. Among the topics covered are fully stressed design and optimality criteria, mathematical programming, smooth and ribbed designs, design against plastic collapse, buckling constraints, and vibration constraints.

Haftka, R. T.↗

Precipitating electron energy flux and auroral zone conductances - An empirical model

Data from the low energy electron (LEE) experiments on the Atmosphere Explorer C and D satellites have been used to determine the average global distribution of the energy flux of precipitating auroral electrons and their average energy for different levels of geomagnetic activity. Measurements from the Atmosphere Explorer unified abstract file (15-s resolution) have been binned according to invariant latitude (in the range 50-90 deg), magnetic local time, and geomagnetic activity as measured by the Kp and auroral electrojet (AE) indices, separately. Bin-averaged values of precipitating energy flux and average energy have been calculated, and a smoothing algorithm used to reduce stochastic variations in the raw data. The results indicate that, for the parameters studied, the AE inces does a superior job of ordering the data with regard to geomagnetic activity. The global distribution of the auroral enhancement porition of the Pedersen and Hall conductances were inferred from the data by means of an empirical fit to detailed energy deposition calculations.

Spiro, R. W.↗

A Strassen-Newton algorithm for high-speed parallelizable matrix inversion

Techniques are described for computing matrix inverses by algorithms that are highly suited to massively parallel computation. The techniques are based on an algorithm suggested by Strassen (1969). Variations of this scheme use matrix Newton iterations and other methods to improve the numerical stability while at the same time preserving a very high level of parallelism. One-processor Cray-2 implementations of these schemes range from one that is up to 55 percent faster than a conventional library routine to one that is slower than a library routine but achieves excellent numerical stability. The problem of computing the solution to a single set of linear equations is discussed, and it is shown that this problem can also be solved efficiently using these techniques.

Bailey, David H.↗

Measurement of a zonal wind profile on Titan by Doppler tracking of the Cassini entry probe

A program, called the Cassini mission, intended to study the Saturn system by utilizing a Saturn orbiter and a probe descending to the surface of Titan, is discussed. Winds are expected to cause perturbations to the probe local horizontal velocity, resulting in an anomalous drift in the probe location and a shift in the frequency of the probe telemetry, due to the Doppler effect. By using an iterative algorithm, in which the time variation of the probe telemetry frequency is monitored throughout the descent, and the probe trajectory is updated to reflect the effect of wind on the probe location, a highly accurate relative wind profile can be recovered. By adding a single wind velocity, measured by independent means, an absolute wind profile can be obtained. However, the accuracy of the zonal winds recovery is limited by errors in trajectory, and frequency.

Atkinson, D. H.↗

Measurement of subresolution terrain displacements using SPOT panchromatic imagery

It is shown that satellite-derived imagery can be used to measure subresolution horizontal terrain displacements associated with present-day earthquakes, sand dune migration, glacial motion, coastal sediment transport, and pre-eruptive volcanic processes. This use of these data allows the detection of change and the determination of rates of many environmental processes worldwide. The premise of the method is that areally extensive, subresolution spatial differences in ground patterns between images acquired at different times can be accurately measured to high precision and can be distinguished from systematic image differences, such as those due to sensing-system attitude variations. Thus far, a preliminary algorithm has been tested that shows that the approach is reasonable.

Crippen, Robert E.↗

New Global Electron Density Observations from GPS-RO in the D- and E-Region Ionosphere

A novel retrieval technique is developed for electron density (N(sub e)) in the D- and E-region (80-120 km) using the high-quality 50-Hz GPS radio occultation (GPS-RO) phase measurements. The new algorithm assumes a slow, linear variation in the F-region background when the GPS-RO passes through the D- and E-region, and extracts the N(sub e) profiles at 80-130 km from the phase advance signal caused by N(sub e). Unlike the conventional Abel function, the new approach produces a sharp N(sub e) weighting function in the lower ionosphere, and the N(sub e) retrievals are in good agreement with the IRI (International Reference Ionosphere) model in terms of monthly maps, zonal means and diurnal variations. The daytime GPS-RO N(sub e) profiles can be well characterized by the alpha-Chapman function of three parameters (N(sub mE), h(sub mE) and H), showing that the bottom of E-region is deepening and sharpening towards the summer pole. At high latitudes the monthly GPS-RO N(sub e) maps at 80-120 km reveal clear enhancement in the auroral zones, more prominent at night, as a result of energetic electron precipitation (EEP) from the outer radiation belt. The D-/E-region auroral N(sub e) is strongly correlated with K(sub p) on a daily basis. The new N(sub e) data allow further comprehensive analyses of the sporadic E (E(sub s)) phenomena in connection with the background N(sub e) in the E-region. The layered (2-10 km) and fluctuated (less than 2 km) E(sub s) components, namely N(sub e_Layer) than N(sub e_Pert), are extracted with respect to the background N( sub e_Region) on a profile-by-profile basis. The N(sub e_Layer) component has a strong but highly-refined peak at approximately 105 km, with an amplitude smaller than N(sub e_Region) approximately by an order of magnitude. The N(sub e_Pert) component, which was studied extensively in the past, is approximately 2 orders of magnitude weaker than N(sub e_Layer). Both N(sub e_Layer) and N(sub e_Pert) are subject to significant diurnal and semidiurnal variations, showing downward progression with local time in amplitude. The 11-year solar cycle dominates the N(sub e) interannual variations, showing larger N(sub e_Region) and N(sub e_Layer) but smaller N(sub e_Pert) amplitudes in the solar maximum years. Enhanced Ne profiles are often observed in the polar winter, showing good correlation with solar proton events (SPEs) and geomagnetic activity. The new methodology offers great potential for retrieving low N(sub e) in the D-region, where radio propagation and communication blackouts can occur due to enhanced ionization. For space weather applications it is recommended for GPSRO operations to raise the top of high-rate data acquisition to approximately 140 km in the future.

GSP-R↗

Computer-aided analysis of Landsat data for surveying Texas coastal zone environments

The feasibility of using machine-aided processing of Landsat data to inventory environmental units was studied by analyzing geometrically corrected and spatially registered Landsat data collected over the Matagorda Bay area of the Texas coastal estuarine system. A clustering algorithm (nonsupervised processor) was used to divide the data into groups of sample points of similar spectral characteristics, and correlation of spectral classes with reference data on a point-to-point basis showed the coastal features exhibit unique spectral variations. Use of a maximum likelihood algorithm permitted discrimination of 13 terrestrial and aquatic environments.

Kristof, S. J.↗

An advancing front Delaunay triangulation algorithm designed for robustness

A new algorithm is described for generating an unstructured mesh about an arbitrary two-dimensional configuration. Mesh points are generated automatically by the algorithm in a manner which ensures a smooth variation of elements, and the resulting triangulation constitutes the Delaunay triangulation of these points. The algorithm combines the mathematical elegance and efficiency of Delaunay triangulation algorithms with the desirable point placement features, boundary integrity, and robustness traditionally associated with advancing-front-type mesh generation strategies. The method offers increased robustness over previous algorithms in that it cannot fail regardless of the initial boundary point distribution and the prescribed cell size distribution throughout the flow-field.

Mavriplis, D. J.↗

An advancing front Delaunay triangulation algorithm designed for robustness

A new algorithm is described for generating an unstructured mesh about an arbitrary two-dimensional configuration. Mesh points are generated automatically by the algorithm in a manner which ensures a smooth variation of elements, and the resulting triangulation constitutes the Delaunay triangulation of these points. The algorithm combines the mathematical elegance and efficiency of Delaunay triangulation algorithms with the desirable point placement features, boundary integrity, and robustness traditionally associated with advancing-front-type mesh generation strategies. The method offers increased robustness over previous algorithms in that it cannot fail regardless of the initial boundary point distribution and the prescribed cell size distribution throughout the flow-field.

Mavriplis, D. J.↗

Fine guidance sensor design optimization for Space Infrared Telescope Facility (SIRTF)

The present paper is concerned with the general issues of design optimization of a staring mosaic sensor for guide-star tracking, taking into account the choice of the proper interpolation algorithm. It was found that the final choice of algorithm depends on the severity of the positional accuracy requirement and the expected variation in the blur spot shape. The considered track algorithms estimate the position either of the peak or of the centroid of the image of a point source. Attention is given to aspects of problem definition, the total effect of error sources on algorithm accuracy, a parametric analysis of performance behavior, Monte Carlo results, and the design of the SIRTF Fine Guidance Sensor (FGS). The algorithm error, the signal-to-noise error, and the nonuniformity error are examined.

Grossman, S. B.↗

Variable Sampling Mapping

The performance of an optical system (for example, a telescope) is limited by the misalignments and manufacturing imperfections of the optical elements in the system. The impact of these misalignments and imperfections can be quantified by the phase variations imparted on light traveling through the system. Phase retrieval is a methodology for determining these variations. Phase retrieval uses images taken with the optical system and using a light source of known shape and characteristics. Unlike interferometric methods, which require an optical reference for comparison, and unlike Shack-Hartmann wavefront sensors that require special optical hardware at the optical system's exit pupil, phase retrieval is an in situ, image-based method for determining the phase variations of light at the system s exit pupil. Phase retrieval can be used both as an optical metrology tool (during fabrication of optical surfaces and assembly of optical systems) and as a sensor used in active, closed-loop control of an optical system, to optimize performance. One class of phase-retrieval algorithms is the iterative transform algorithm (ITA). ITAs estimate the phase variations by iteratively enforcing known constraints in the exit pupil and at the detector, determined from modeled or measured data. The Variable Sampling Mapping (VSM) technique is a new method for enforcing these constraints in ITAs. VSM is an open framework for addressing a wide range of issues that have previously been considered detrimental to high-accuracy phase retrieval, including undersampled images, broadband illumination, images taken at or near best focus, chromatic aberrations, jitter or vibration of the optical system or detector, and dead or noisy detector pixels. The VSM is a model-to-data mapping procedure. In VSM, fully sampled electric fields at multiple wavelengths are modeled inside the phase-retrieval algorithm, and then these fields are mapped to intensities on the light detector, using the properties of the detector and optical system, for comparison with measured data. Ultimately, this model-to-data mapping procedure enables a more robust and accurate way of incorporating the exit-pupil and image detector constraints, which are fundamental to the general class of ITA phase retrieval algorithms.

Smith, Jeffrey, S.↗

Improvements to Lunar BRDF-Corrected Nighttime Satellite Imagery: Uses and Applications

Observations made by the VIIRS day/night band (DNB) provide daily, nighttime measurements to monitor Earth surface processes.However, these observations are impacted by variations in reflected solar radiation on the moon's surface. As the moon transitions from new to full phase, increasing radiance is reflected to the Earth's surface and contributes additional reflected moonlight from clouds and land surface, in addition to emissions from other light sources observed by the DNB. The introduction of a bi-directional reflectance distribution function (BRDF) algorithm serves to remove these lunar variations and normalize observed radiances. Provided by the Terrestrial Information Systems Laboratory at Goddard Space Flight Center, a 1 km gridded lunar BRDF-corrected DNB product and VIIRS cloud mask can be used for a multitude of nighttime applications without influence from the moon. Such applications include the detection of power outages following severe weather events using pre-and post-event DNB imagery, as well as the identification of boat features to curtail illegal fishing practices. This presentation will provide context on the importance of the lunar BRDF correction algorithm and explore the aforementioned uses of this improved DNB product for applied science applications.

NWS↗

Angle-of-Attack-Modulated Terminal Point Control for Neptune Aerocapture

An aerocapture guidance algorithm based on a calculus of variations approach is developed, using angle of attack as the primary control variable. Bank angle is used as a secondary control to alleviate angle of attack extremes and to control inclination. The guidance equations are derived in detail. The controller has very small onboard computational requirements and is robust to atmospheric and aerodynamic dispersions. The algorithm is applied to aerocapture at Neptune. Three versions of the controller are considered with varying angle of attack authority. The three versions of the controller are evaluated using Monte Carlo simulations with expected dispersions.

Queen, Eric M.↗

Atmospheric Ascent Guidance for Rocket-Powered Launch Vehicles

An advanced ascent guidance algorithm for rocket- powered launch vehicles is developed. This algorithm cyclically solves the calculus-of-variations two-point boundary-value problem starting at vertical rise completion through main engine cutoff. This is different from traditional ascent guidance algorithms which operate in a simple open-loop mode until high dynamic pressure (including the critical max-Q) portion of the trajectory is over, at which time guidance operates under the assumption of negligible aerodynamic acceleration (i.e., vacuum dynamics). The initial costate guess is corrected based on errors in the terminal state constraints and the transversality conditions. Judicious approximations are made to reduce the order and complexity of the state/costate system. Results comparing guided launch vehicle trajectories with POST open-loop trajectories are given verifying the basic formulation of the algorithm. Multiple shooting is shown to be a very effective numerical technique for this application. In particular, just one intermediate shooting point, in addition to the initial shooting point, is sufficient to significantly reduce sensitivity to the guessed initial costates. Simulation results from a high-fidelity trajectory simulation are given for the case of launch to sub-orbital cutoff conditions as well as launch to orbit conditions. An abort to downrange landing site formulation of the algorithm is presented.

Dukeman, Greg A.↗

High-Lift Optimization Design Using Neural Networks on a Multi-Element Airfoil

The high-lift performance of a multi-element airfoil was optimized by using neural-net predictions that were trained using a computational data set. The numerical data was generated using a two-dimensional, incompressible, Navier-Stokes algorithm with the Spalart-Allmaras turbulence model. Because it is difficult to predict maximum lift for high-lift systems, an empirically-based maximum lift criteria was used in this study to determine both the maximum lift and the angle at which it occurs. Multiple input, single output networks were trained using the NASA Ames variation of the Levenberg-Marquardt algorithm for each of the aerodynamic coefficients (lift, drag, and moment). The artificial neural networks were integrated with a gradient-based optimizer. Using independent numerical simulations and experimental data for this high-lift configuration, it was shown that this design process successfully optimized flap deflection, gap, overlap, and angle of attack to maximize lift. Once the neural networks were trained and integrated with the optimizer, minimal additional computer resources were required to perform optimization runs with different initial conditions and parameters. Applying the neural networks within the high-lift rigging optimization process reduced the amount of computational time and resources by 83% compared with traditional gradient-based optimization procedures for multiple optimization runs.

Greenman, Roxana M.↗

An historical survey of computational methods in optimal control.

Review of some of the salient theoretical developments in the specific area of optimal control algorithms. The first algorithms for optimal control were aimed at unconstrained problems and were derived by using first- and second-variation methods of the calculus of variations. These methods have subsequently been recognized as gradient, Newton-Raphson, or Gauss-Newton methods in function space. A much more recent addition to the arsenal of unconstrained optimal control algorithms are several variations of conjugate-gradient methods. At first, constrained optimal control problems could only be solved by exterior penalty function methods. Later algorithms specifically designed for constrained problems have appeared. Among these are methods for solving the unconstrained linear quadratic regulator problem, as well as certain constrained minimum-time and minimum-energy problems. Differential-dynamic programming was developed from dynamic programming considerations. The conditional-gradient method, the gradient-projection method, and a couple of feasible directions methods were obtained as extensions or adaptations of related algorithms for finite-dimensional problems. Finally, the so-called epsilon-methods combine the Ritz method with penalty function techniques.

Polak, E.↗

An automatic editing algorithm for GPS data

An algorithm has been developed to edit automatically Global Positioning System data such that outlier deletion, cycle slip identification, and correction are independent of clock instability, selective availability, receiver-satellite kinematics, and tropospheric conditions. This algorithm, called TurboEdit, operates on undifferenced, dual frequency carrier phase data, and requires the use of P code pseudorange data and a smoothly varying ionospheric electron content. TurboEdit was tested on the large data set from the CASA Uno experiment, which contained over 2500 cycle slips.Analyst intervention was required on 1 percent of the station-satellite passes, almost all of these problems being due to difficulties in extrapolating variations in the ionospheric delay. The algorithm is presently being adapted for real time data editing in the Rogue receiver for continuous monitoring applications.

Blewitt, Geoffrey↗