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 37 records · Page 2

An investigation of using an RQP based method to calculate parameter sensitivity derivatives

Estimation of the sensitivity of problem functions with respect to problem variables forms the basis for many of our modern day algorithms for engineering optimization. The most common application of problem sensitivities has been in the calculation of objective function and constraint partial derivatives for determining search directions and optimality conditions. A second form of sensitivity analysis, parameter sensitivity, has also become an important topic in recent years. By parameter sensitivity, researchers refer to the estimation of changes in the modeling functions and current design point due to small changes in the fixed parameters of the formulation. Methods for calculating these derivatives have been proposed by several authors (Armacost and Fiacco 1974, Sobieski et al 1981, Schmit and Chang 1984, and Vanderplaats and Yoshida 1985). Two drawbacks to estimating parameter sensitivities by current methods have been: (1) the need for second order information about the Lagrangian at the current point, and (2) the estimates assume no change in the active set of constraints. The first of these two problems is addressed here and a new algorithm is proposed that does not require explicit calculation of second order information.

Beltracchi, Todd J.↗

Preliminary investigation of pilot scanning techniques of dial pointing instruments

Two pilots' methods of looking at instruments with needle pointers in a fixed base helicopter simulation were observed. A total of 45 runs were analyzed for each pilot. The data indicated that two apparently different techniques were being used; one looking at the needle point, the other looking at a fixed spot on the instrument and reading the needle direction parafoveally. The latter technique is found to be somewhat faster with both pilots accomplishing the flying task.

Harris, R. L., Sr.↗

Calculation of free-fall trajectories using numerical optimization methods.

An important problem in space flight is the calculation of trajectories for nonthrusting vehicles between fixed points in a given time. A new procedure based on Hamilton's principle for solving such two-point boundary-value problems is presented. It employs numerical optimization methods to perform the extremization required by Hamilton's principle. This procedure is applied to the calculation of an Earth-Moon trajectory. The results show that the initial guesses required to obtain an iteration procedure which converges are not critical and that convergence can be obtained to any predetermined degree of accuracy.

Hull, D. G.↗

An exploratory investigation of weight estimation techniques for hypersonic flight vehicles

The three basic methods of weight prediction (fixed-fraction, statistical correlation, and point stress analysis) and some of the computer programs that have been developed to implement them are discussed. A modified version of the WAATS (Weights Analysis of Advanced Transportation Systems) program is presented, along with input data forms and an example problem.

Cook, E. L.↗

Optimal active control for Burgers equations

A method for active fluid flow control based on control theory is discussed. Dynamic programming and fixed point successive approximations are used to accommodate the nonlinear control problem. The long-term goal of this project is to establish an effective method applicable to complex flows such as turbulence and jets. However, in this report, the method is applied to stochastic Burgers equation as an intermediate step towards this goal. Numerical results are compared with those obtained by gradient search methods.

Ikeda, Yutaka↗

Some estimation formulae for continuous time-invariant linear systems

In this brief paper we examine a Riccati equation decomposition due to Reid and Lainiotis and apply the result to the continuous time-invariant linear filtering problem. Exploitation of the time-invariant structure leads to integration-free covariance recursions which are of use in covariance analyses and in filter implementations. A super-linearly convergent iterative solution to the algebraic Riccati equation (ARE) is developed. The resulting algorithm, arranged in a square-root form, is thought to be numerically stable and competitive with other ARE solution methods. Certain covariance relations that are relevant to the fixed-point and fixed-lag smoothing problems are also discussed.

Bierman, G. J.↗

Crustal deformation measurements in central Japan determined by a Global Positioning System fixed-point network

Results are presented from temporally dense measurements of crustal deformation associated with the convergence of the Eurasian (EUR), Pacific, North American, and Philippine Sea (PHS) plates, carried out in April 1988 by a 10-station GPS fixed-point network established in central Japan. Using regional orbit relaxation methods, the analysis of the first 17-month data revealed significant horizontal deformation across the Suruga trough. Namely, it was found that a site in the northern tip of PHS plate moved nearly westward with a velocity of 28 +/-5 mm per year, and a site at the southeastern tip of EUR plate moved south-southwestward with a velocity of 18 +/-5 mm per year. A significant vertical uplift with a velocity of 20 mm/yr was detected at a site inland of the Tokai district located in the Akaishi uplift zone and at a site on the Hatsushima Island in Sagami Bay.

Shimada, Seiichi↗

Fast Multigrid Techniques in Total Variation-Based Image Reconstruction

Existing multigrid techniques are used to effect an efficient method for reconstructing an image from noisy, blurred data. Total Variation minimization yields a nonlinear integro-differential equation which, when discretized using cell-centered finite differences, yields a full matrix equation. A fixed point iteration is applied with the intermediate matrix equations solved via a preconditioned conjugate gradient method which utilizes multi-level quadrature (due to Brandt and Lubrecht) to apply the integral operator and a multigrid scheme (due to Ewing and Shen) to invert the differential operator. With effective preconditioning, the method presented seems to require Omicron(n) operations. Numerical results are given for a two-dimensional example.

Oman, Mary Ellen↗

Hybrid state vector methods for structural dynamic and aeroelastic boundary value problems

A computational technique is developed that is suitable for performing preliminary design aeroelastic and structural dynamic analyses of large aspect ratio lifting surfaces. The method proves to be quite general and can be adapted to solving various two point boundary value problems. The solution method, which is applicable to both fixed and rotating wing configurations, is based upon a formulation of the structural equilibrium equations in terms of a hybrid state vector containing generalized force and displacement variables. A mixed variational formulation is presented that conveniently yields a useful form for these state vector differential equations. Solutions to these equations are obtained by employing an integrating matrix method. The application of an integrating matrix provides a discretization of the differential equations that only requires solutions of standard linear matrix systems. It is demonstrated that matrix partitioning can be used to reduce the order of the required solutions. Results are presented for several example problems in structural dynamics and aeroelasticity to verify the technique and to demonstrate its use. These problems examine various types of loading and boundary conditions and include aeroelastic analyses of lifting surfaces constructed from anisotropic composite materials.

Lehman, L. L.↗

Sinusoidal transform coding

It has been shown that an analysis/synthesis system based on a sinusoidal representation of speech leads to synthetic speech that is essentially perceptually indistinguishable from the original. Strategies for coding the amplitudes, frequencies and phases of the sine waves have been developed that have led to a multirate coder operating at rates from 2400 to 9600 bps. The encoded speech is highly intelligible at all rates with a uniformly improving quality as the data rate is increased. A real-time fixed-point implementation has been developed using two ADSP2100 DSP chips. The methods used for coding and quantizing the sine-wave parameters for operation at the various frame rates are described.

Mcaulay, Robert J.↗

Automated Parameter Studies Using a Cartesian Method

Computational Fluid Dynamics (CFD) is now routinely used to analyze isolated points in a design space by performing steady-state computations at fixed flight conditions (Mach number, angle of attack, sideslip), for a fixed geometric configuration of interest. This "point analysis" provides detailed information about the flowfield, which aides an engineer in understanding, or correcting, a design. A point analysis is typically performed using high fidelity methods at a handful of critical design points, e.g. a cruise or landing configuration, or a sample of points along a flight trajectory.

Murman, Scott M.↗

A strawman targeting technique for spinner spacecraft in the presence of pointing uncertainty

A graphical method that allows the selection of the attitude orientation of a spinner spacecraft prior to the application of a single fixed magnitude velocity change to effect a desired orbit transfer is presented. The selection of the final pointing vector along or opposite to which the desired Delta-V is applied, takes into account pointing errors due to initial attitude determination uncertainty and execution errors introduced during the rhumb line precession of the spin axis. The characterization of the target space in terms of any achieved orbit parameters of interest as well as the final pointing uncertainty is sufficient to determine the right final attitude that insures satisfaction of mission requirements.

Kechichian, J. A.↗

First Terrestrial Soft X-ray Aurora Observations by Chandra

Northern polar "auroral" regions of Earth was observed by High-Resolution Camera in imaging mode (T32C-I) aboard Chandra X-Ray Observatory (CXO) during mid December 2003 - mid April 2004. Ten CXO observations, each approximately 20 min duration, were made in a non-conventional method (due to CXO technical issues), such that Chandra was aimed at a fixed point in sky and the Earth's polar cusp was allowed to drift through the HRC-I field-of-view. The observations were performed when CXO was near apogee and timed during northern winter mostly near midnight (6 hr), except two observations which occurred around 1200 UT, so that northern polar region is entirely in dark and solar fluoresced x-ray contamination can be avoided. These observations were aimed at searching the Earth's soft x-ray aurora and to do a comparative study with Jupiter's x-ray aurora, where a pulsating x-ray hot-spot near the northern magnetic pole has been observed by Chandra that implies a particle source region near Jupiter's magnetopause, and entry of heavy solar wind ions due to high-latitude reconnection as a viable explanation for the soft x-ray emissions. The first Chandra soft (0.1-2 keV) x-ray observations of Earth's aurora show that it is highly variable (intense arc, multiple arcs, diffuse, at times almost absent). In at least one of the observations an isolated blob of emission is observed where we expect cusp to be: giving indication of solar wind charge-exchange signature in x-rays. We are comparing the Chandra x-ray observations with observations at other wavelengths and particle data from Earth-orbiting satellites and solar wind measurements from near-Earth ACE and SOH0 spacecraft. Preliminary results from these unique CXO-Earth observations will be presented and discussed.

Bhardwaj, Anil↗

Bifurcation Analysis Using Rigorous Branch and Bound Methods

For the study of nonlinear dynamic systems, it is important to locate the equilibria and bifurcations occurring within a specified computational domain. This paper proposes a new approach for solving these problems and compares it to the numerical continuation method. The new approach is based upon branch and bound and utilizes rigorous enclosure techniques to yield outer bounding sets of both the equilibrium and local bifurcation manifolds. These sets, which comprise the union of hyper-rectangles, can be made to be as tight as desired. Sufficient conditions for the existence of equilibrium and bifurcation points taking the form of algebraic inequality constraints in the state-parameter space are used to calculate their enclosures directly. The enclosures for the bifurcation sets can be computed independently of the equilibrium manifold, and are guaranteed to contain all solutions within the computational domain. A further advantage of this method is the ability to compute a near-maximally sized hyper-rectangle of high dimension centered at a fixed parameter-state point whose elements are guaranteed to exclude all bifurcation points. This hyper-rectangle, which requires a global description of the bifurcation manifold within the computational domain, cannot be obtained otherwise. A test case, based on the dynamics of a UAV subject to uncertain center of gravity location, is used to illustrate the efficacy of the method by comparing it with numerical continuation and to evaluate its computational complexity.

Smith, Andrew P.↗

Recent Advances to Estimation of Fixed-Interface Modal Models Using Dynamic Substructuring

In 2010, Allen & Mayes proposed to estimate the fixed-interface modes of a structure by measuring the modes of the structure bolted to a fixture and then applying constraints to the fixture using the transmission simulator method. While the method proved useful, and has indeed been used in studies since that point, a few peculiarities were noted. First, in some cases the estimated fixed-base natural frequencies were observed to converge very slowly to the true values (in simulated experiments) as the number of constraints was increased. To formulate these constraints, prior studies used only the free-interface modes of the fixture or the measured modes of the assembly. This work extends that to consider other sets of constraints, showing improved results. Furthermore, in some prior studies it has been observed that there were errors of more than 10% in the natural frequencies even when the fixture motion was hundreds of times smaller than the motion of the structure of interest (and so it had presumably been removed). This work explores this phenomenon, seeking to use the strain energy in the fixture, to the extent that it can be estimated using a test-analysis model for the fixture, as a metric to predict frequency error. The proposed methods are explored by applying them to simulated measurements from a beam and from the NASA Space Launch System coupled to the Mobile Launcher.

Allen, Mathew S.↗

Momentum Flux Determination Using the Multi-beam Poker Flat Incoherent Scatter Radar

In this paper, we develop an estimator for the vertical flux of horizontal momentum with arbitrary beam pointing, applicable to the case of arbitrary but fixed beam pointing with systems such as the Poker Flat Incoherent Scatter Radar (PFISR). This method uses information from all available beams to resolve the variances of the wind field in addition to the vertical flux of both meridional and zonal momentum, targeted for high-frequency wave motions. The estimator utilises the full covariance of the distributed measurements, which provides a significant reduction in errors over the direct extension of previously developed techniques and allows for the calculation of an error covariance matrix of the estimated quantities. We find that for the PFISR experiment, we can construct an unbiased and robust estimator of the momentum flux if sufficient and proper beam orientations are chosen, which can in the future be optimized for the expected frequency distribution of momentum-containing scales. However, there is a potential trade-off between biases and standard errors introduced with the new approach, which must be taken into account when assessing the momentum fluxes. We apply the estimator to PFISR measurements on 23 April 2008 and 21 December 2007, from 60-85 km altitude, and show expected results as compared to mean winds and in relation to the measured vertical velocity variances.

Momentum flux↗

Simultaneous quaternion estimation (QUEST) and bias determination

Tests of a new method for the simultaneous estimation of spacecraft attitude and sensor biases, based on a quaternion estimation algorithm minimizing Wahba's loss function are presented. The new method is compared with a conventional batch least-squares differential correction algorithm. The estimates are based on data from strapdown gyros and star trackers, simulated with varying levels of Gaussian noise for both inertially-fixed and Earth-pointing reference attitudes. Both algorithms solve for the spacecraft attitude and the gyro drift rate biases. They converge to the same estimates at the same rate for inertially-fixed attitude, but the new algorithm converges more slowly than the differential correction for Earth-pointing attitude. The slower convergence of the new method for non-zero attitude rates is believed to be due to the use of an inadequate approximation for a partial derivative matrix. The new method requires about twice the computational effort of the differential correction. Improving the approximation for the partial derivative matrix in the new method is expected to improve its convergence at the cost of increased computational effort.

Markley, F. Landis↗