Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fixed point method”

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 73 records · Page 4

A system-approach to the elastohydrodynamic lubrication point-contact problem

The classical EHL (elastohydrodynamic lubrication) point contact problem is solved using a new system-approach, similar to that introduced by Houpert and Hamrock for the line-contact problem. Introducing a body-fitted coordinate system, the troublesome free-boundary is transformed to a fixed domain. The Newton-Raphson method can then be used to determine the pressure distribution and the cavitation boundary subject to the Reynolds boundary condition. This method provides an efficient and rigorous way of solving the EHL point contact problem with the aid of a supercomputer and a promising method to deal with the transient EHL point contact problem. A typical pressure distribution and film thickness profile are presented and the minimum film thicknesses are compared with the solution of Hamrock and Dowson. The details of the cavitation boundaries for various operating parameters are discussed.

Lim, Sang Gyu↗

A system-approach to the elastohydrodynamic lubrication point-contact problem

The classical EHL (elastohydrodynamic lubrication) point contact problem is solved using a new system-approach, similar to that introduced by Houpert and Hamrock for the line-contact problem. Introducing a body-fitted coordinate system, the troublesome free-boundary is transformed to a fixed domain. The Newton-Raphson method can then be used to determine the pressure distribution and the cavitation boundary subject to the Reynolds boundary condition. This method provides an efficient and rigorous way of solving the EHL point contact problem with the aid of a supercomputer and a promising method to deal with the transient EHL point contact problem. A typical pressure distribution and film thickness profile are presented and the minimum film thicknesses are compared with the solution of Hamrock and Dowson. The details of the cavitation boundaries for various operating parameters are discussed.

Lim, Sang G.↗

Parameter estimation for slit-type scanning sensors

The Infrared Astronomical Satellite, scheduled for launch into a 900 km near-polar orbit in August 1982, will perform an infrared point source survey by scanning the sky with slit-type sensors. The description of position information is shown to require the use of a non-Gaussian random variable. Methods are described for deciding whether separate detections stem from a single common source, and a formulism is developed for the scan-to-scan problems of identifying multiple sightings of inertially fixed point sources for combining their individual measurements into a refined estimate. Several cases are given where the general theory yields results which are quite different from the corresponding Gaussian applications, showing that argument by Gaussian analogy would lead to error.

Fowler, J. W.↗

An adaptive identification and control scheme for large space structures

A unified identification and control scheme capable of achieving space at form performance objectives under nominal or failure conditions is described. Preliminary results are also presented, showing that the methodology offers much promise for effective robust control of large space structures. The control method is a multivariable, adaptive, output predictive controller called Model Predictive Control (MPC). MPC uses a state space model and input reference trajectories of set or tracking points to adaptively generate optimum commands. For a fixed model, MPC processes commands with great efficiency, and is also highly robust. A key feature of MPC is its ability to control either nonminimum phase or open loop unstable systems. As an output controller, MPC does not explicitly require full state feedback, as do most multivariable (e.g., Linear Quadratic) methods. Its features are very useful in LSS operations, as they allow non-collocated actuators and sensors. The identification scheme is based on canonical variate analysis (CVA) of input and output data. The CVA technique is particularly suited for the measurement and identification of structural dynamic processes - that is, unsteady transient or dynamically interacting processes such as between aerodynamics and structural deformation - from short, noisy data. CVA is structured so that the identification can be done in real or near real time, using computationally stable algorithms. Modeling LSS dynamics in 1-g laboratories has always been a major impediment not only to understanding their behavior in orbit, but also to controlling it. In cases where the theoretical model is not confirmed, current methods provide few clues concerning additional dynamical relationships that are not included in the theoretical models. CVA needs no a priori model data, or structure; all statistically significant dynamical states are determined using natural, entropy-based methods. Heretofore, a major limitation in applying adaptive controllers to LSS problems has been the requirement for accurate models. We feel that MPC coupled with the CVA identification scheme will resolve this limitation. The results are based on research performed under an SBIR grant, administered by the Air Force Office of Scientific Research.

Carroll, J. V.↗

Modified Interior Distance Functions (Theory and Methods)

In this paper we introduced and developed the theory of Modified Interior Distance Functions (MIDF's). The MIDF is a Classical Lagrangian (CL) for a constrained optimization problem which is equivalent to the initial one and can be obtained from the latter by monotone transformation both the objective function and constraints. In contrast to the Interior Distance Functions (IDF's), which played a fundamental role in Interior Point Methods (IPM's), the MIDF's are defined on an extended feasible set and along with center, have two extra tools, which control the computational process: the barrier parameter and the vector of Lagrange multipliers. The extra tools allow to attach to the MEDF's very important properties of Augmented Lagrangeans. One can consider the MIDFs as Interior Augmented Lagrangeans. It makes MIDF's similar in spirit to Modified Barrier Functions (MBF's), although there is a fundamental difference between them both in theory and methods. Based on MIDF's theory, Modified Center Methods (MCM's) have been developed and analyzed. The MCM's find an unconstrained minimizer in primal space and update the Lagrange multipliers, while both the center and the barrier parameter can be fixed or updated at each step. The MCM's convergence was investigated, and their rate of convergence was estimated. The extension of the feasible set and the special role of the Lagrange multipliers allow to develop MCM's, which produce, in case of nondegenerate constrained optimization, a primal and dual sequences that converge to the primal-dual solutions with linear rate, even when both the center and the barrier parameter are fixed. Moreover, every Lagrange multipliers update shrinks the distance to the primal dual solution by a factor 0 less than gamma less than 1 which can be made as small as one wants by choosing a fixed interior point as a 'center' and a fixed but large enough barrier parameter. The numericai realization of MCM leads to the Newton MCM (NMCM). The approximation for the primal minimizer one finds by Newton Method followed by the Lagrange multipliers update. Due to the MCM convergence, when both the center and the barrier parameter are fixed, the condition of the MDF Hessism and the neighborhood of the primal ninimizer where Newton method is 'well' defined remains stable. It contributes to both the complexity and the numerical stability of the NMCM.

Polyak, Roman A.↗

An Efficient Ray-Tracing Method for Determining Terrain Intercepts in EDL Simulations

The calculation of a ray's intercept from an arbitrary point in space to a prescribed surface is a common task in computer simulations. The arbitrary point often represents an object that is moving according to the simulation, while the prescribed surface is fixed in a defined frame. For detailed simulations, this surface becomes complex, taking the form of real-world objects such as mountains, craters or valleys which require more advanced methods to accurately calculate a ray's intercept location. Incorporation of these complex surfaces has commonly been implemented in graphics systems that utilize highly optimized graphics processing units to analyze such features. This paper proposes a simplified method that does not require computationally intensive graphics solutions, but rather an optimized ray-tracing method for an assumed terrain dataset. This approach was developed for the Mars Science Laboratory mission which landed on the complex terrain of Gale Crater. First, this paper begins with a discussion of the simulation used to implement the model and the applicability of finding surface intercepts with respect to atmosphere modeling, altitude determination, radar modeling, and contact forces influencing vehicle dynamics. Next, the derivation and assumptions of the intercept finding method are presented. Key assumptions are noted making the routines specific to only certain types of surface data sets that are equidistantly spaced in longitude and latitude. The derivation of the method relies on ray-tracing, requiring discussion on the formulation of the ray with respect to the terrain datasets. Further discussion includes techniques for ray initialization in order to optimize the intercept search. Then, the model implementation for various new applications in the simulation are demonstrated. Finally, a validation of the accuracy is presented along with the corresponding data sets used in the validation. A performance summary of the method will be shown using the analysis from the Mars Science Laboratory's terminal descent sensing model. Alternate uses will also be shown for determining horizon maps and orbiter set times.

Shidner, Jeremy D.↗

Error modeling of precision orientation sensors in a fixed base simulation

Models of noise and dynamic characteristics of gyro and autocollimator for very small signal levels are presented. Measurements were evaluated using spectral techniques for identifying noise from base motion. The experiment was constructed to measure the precession, due to relativistic effects, of an extremely precise earth-orbiting gyroscope. The design goal for nonrelativistic gyro drift is 0.001 arcsec per year. An analogous fixed base simulator was used in developing methods of instrument error modeling and performance evaluation applicable to the relativity experiment sensors and other precision pointing instruments. Analysis of autocollimator spectra uncovered the presence of a platform gimbal resonance. The source of resonance was isolated to gimbal bearing elastic restraint properties most apparent at very small levels of motion. A model of these properties which include both elastic and coulomb friction characteristics is discussed, and a describing function developed.

Klinger, D. L.↗

Desingularization of periodic vortex sheet roll-up

An analytical approach is used in an attempt to model the evolution of a vortex sheet past the critical time by means of a desingularization method. Evolution of the sheet, which is embedded in a two-dimensional flow, is described by approximating the total circulation between a fixed material point and an imaginary point on the curve of the flow. A linear dispersion relation is defined which shows that the short wavelength modes of the sheet are not unstable and therefore do adversely effect the computations as artifacts of round-off errors. The desingularization approach is demonstrated to converge beyond the critical time for the vortex sheet. Application of the technique for the study of the vortex sheet shed from an elliptically loaded wing is indicated.

Krasny, R.↗

Predictive Lateral Logic for Numerical Entry Guidance Algorithms

Recent entry guidance algorithm development123 has tended to focus on numerical integration of trajectories onboard in order to evaluate candidate bank profiles. Such methods enjoy benefits such as flexibility to varying mission profiles and improved robustness to large dispersions. A common element across many of these modern entry guidance algorithms is a reliance upon the concept of Apollo heritage lateral error (or azimuth error) deadbands in which the number of bank reversals to be performed is non-deterministic. This paper presents a closed-loop bank reversal method that operates with a fixed number of bank reversals defined prior to flight. However, this number of bank reversals can be modified at any point, including in flight, based on contingencies such as fuel leaks where propellant usage must be minimized.

Smith, Kelly M.↗

NASTRAN hydroelastic modal studies. Volume 3: NASTRAN 3-D hydroelastic analysis and MESHGEN user's manuals

The method used to formulate the fluid/structure equations is described in the updates to the NASTRAN Theoretical Manual. The basis for defining the fluid are three dimensional finite elements connected to fluid grid points defining the Eulerian pressure at a point fixed in space. The use of a pressure single degree of freedom at each point rather than three displacements allows a finer mesh of elements with a reasonable matrix order. In the formulation of the fluid/structure system the interior fluid degrees of freedom are transformed and removed from the solution matrices. The eigenvalues of the combination are extracted from small, fully dense, symmetric mass and stiffness matrices, efficiently processed with the Givens method. The solution matrices are defined by the free surface displacements and the reduced structure coordinates.

Source record↗

Optimization of simple structures with higher mode frequency constraints.

Results of a study of least weight optimization of simple structures with a single natural frequency constraint. Variational techniques were used to derive the necessary equations, and numerical methods were used, where required, to find solutions to the resulting nonlinear, two-point boundary value problems. The numerical results suggest a similarity between the fixed fundamental frequency solution and solutions in which a single frequency other than fundamental is held fixed. Such a similarity is shown to occur under special circumstances, and the knowledge of the fundamental solution in these cases makes it possible to calculate the solution to problems with constraints on frequencies other than the fundamental.

Weisshaar, T. A.↗

Electron temperature measurements using a 12-channel array probe

The most common technique for determining the mean kinetic energy of electrons in low-temperature plasmas utilizes the so-called Langmuir probe. The present investigation is concerned with the study of the electron temperature in a pulsed high-beta plasma as a function of time by means of a miniature array of 12 planar Langmuir probes. By using a probe with many individually collecting surfaces, each biased at a different fixed voltage, an approximation to the true probe characteristic is obtained when monitoring the collected currents for each independent subprobe. The employed method, by using many points, provides enough information to reconstruct the entire I-V curve. Attention is given to the principle of operation of the new probe, calibration, accuracy, time resolution, and applications of the new method. By employing two probes it is possible to perform correlation measurements to study heat flow and temperature fluctuations.

Wild, N.↗

PEGASUS 5: An Automated Pre-Processor for Overset-Grid CFD

An all new, automated version of the PEGASUS software has been developed and tested. PEGASUS provides the hole-cutting and connectivity information between overlapping grids, and is used as the final part of the grid generation process for overset-grid computational fluid dynamics approaches. The new PEGASUS code (Version 5) has many new features: automated hole cutting; a projection scheme for fixing gaps in overset surfaces; more efficient interpolation search methods using an alternating digital tree; hole-size optimization based on adding additional layers of fringe points; and an automatic restart capability. The new code has also been parallelized using the Message Passing Interface standard. The parallelization performance provides efficient speed-up of the execution time by an order of magnitude, and up to a factor of 30 for very large problems. The results of three example cases are presented: a three-element high-lift airfoil, a generic business jet configuration, and a complete Boeing 777-200 aircraft in a high-lift landing configuration. Comparisons of the computed flow fields for the airfoil and 777 test cases between the old and new versions of the PEGASUS codes show excellent agreement with each other and with experimental results.

Suhs, Norman E.↗

Results of continuous synchronous orbit testing of sealed nickel-cadmium cells

Test results from continuous synchronous orbit testing of sealed nickel cadmium cells are presented. The synchronous orbit regime simulates a space satellite maintaining a position over a fixed point on earth as the earth rotates on its axis and revolves about the sun. Characteristics of each lot of cells, test conditions, and charge control methods are described.

Harkness, J. D.↗

A microcomputer-based position updating system for general aviation utilizing Loran-C

Modern digital electronic technology is used to produce a device to convert LORAN C to useful pilot information using a simple software algebra and low cost microprocessor devices. Results indicate that the processor based LORAN C navigator has an accuracy of 1.0 nm or less over an area typically covered by a triad of Loran C stations and can execute a position update in less than 0.2 seconds. The system was tested in 30 hours of flight and proved that it can give reliable and accurate navigation information. Methods of converting time differences to position, design considerations for the microcomputer system, and the system for coordinate conversion are discussed. Testing with predetermined points and possible fixes for errors are also considered.

Fischer, J. P.↗

Search-based optimization

The problem of determining the minimum cost hypothetical ancestral sequences for a given cladogram is known to be NP-complete (Wang and Jiang, 1994). Traditionally, point estimations of hypothetical ancestral sequences have been used to gain heuristic, upper bounds on cladogram cost. These include procedures with such diverse approaches as non-additive optimization of multiple sequence alignment, direct optimization (Wheeler, 1996), and fixed-state character optimization (Wheeler, 1999). A method is proposed here which, by extending fixed-state character optimization, replaces the estimation process with a search. This form of optimization examines a diversity of potential state solutions for cost-efficient hypothetical ancestral sequences and can result in greatly more parsimonious cladograms. Additionally, such an approach can be applied to other NP-complete phylogenetic optimization problems such as genomic break-point analysis. c2003 The Willi Hennig Society. Published by Elsevier Science (USA). All rights reserved.

Non-NASA Center↗

Remote sensing of biomass and annual net aerial primary productivity of a salt marsh

Net aerial primary productivity is the rate of storage of organic matter in above-ground plant issues exceeding the respiratory use by the plants during the period of measurement. It is pointed out that this plant tissue represents the fixed carbon available for transfer to and consumption by the heterotrophic organisms in a salt marsh or the estuary. One method of estimating annual net aerial primary productivity (NAPP) required multiple harvesting of the marsh vegetation. A rapid nondestructive remote sensing technique for estimating biomass and NAPP would, therefore, be a significant asset. The present investigation was designed to employ simple regression models, equating spectral radiance indices with Spartina alterniflora biomass to nondestructively estimate salt marsh biomass. The results of the study showed that the considered approach can be successfully used to estimate salt marsh biomass.

Hardisky, M. A.↗

Analysis of a method for precisely relating a seafloor point to a distant point on land

A study of the environmental constraints and engineering aspects of the acoustic portion of a system for making geodetic ties between undersea reference points and others on land is described. Important areas in which to make such observations initially would be from the California mainland out to oceanic points seaward of the San Andreas fault, and across the Aleutian Trench. The overall approach would be to operate a GPS receiver in a relative positioning (interferometric) mode to provide the long range element of the baseline determination (10 to 1,000 km) and an array of precision sea floor acoustic transponders to link the locally moving sea surface GPS antenna location to a fixed sea floor point. Analyses of various environmental constrants (tides, waves, currents, sound velocity variations) lead to the conclusion that, if one uses a properly designed transponder having a remotely controllable precise retransmission time delay, and is careful with regard to methods for installing these on the sea floor, one should, in many ocean locations, be able to achieve sub-decimeter overall system accuracy. Achievements of cm accuracy or better will require additional understanding of time and space scales of variation of sound velocity structure in the ocean at relevant locations.

Spiess, F. N.↗