Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Perturbative methods”

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,351 records · Page 75

Direct Fault Tolerant RLV Altitude Control: A Singular Perturbation Approach

In this paper, we present a direct fault tolerant control (DFTC) technique, where by "direct" we mean that no explicit fault identification is used. The technique will be presented for the attitude controller (autopilot) for a reusable launch vehicle (RLV), although in principle it can be applied to many other applications. Any partial or complete failure of control actuators and effectors will be inferred from saturation of one or more commanded control signals generated by the controller. The saturation causes a reduction in the effective gain, or bandwidth of the feedback loop, which can be modeled as an increase in singular perturbation in the loop. In order to maintain stability, the bandwidth of the nominal (reduced-order) system will be reduced proportionally according to the singular perturbation theory. The presented DFTC technique automatically handles momentary saturations and integrator windup caused by excessive disturbances, guidance command or dispersions under normal vehicle conditions. For multi-input, multi-output (MIMO) systems with redundant control effectors, such as the RLV attitude control system, an algorithm is presented for determining the direction of bandwidth cutback using the method of minimum-time optimal control with constrained control in order to maintain the best performance that is possible with the reduced control authority. Other bandwidth cutback logic, such as one that preserves the commanded direction of the bandwidth or favors a preferred direction when the commanded direction cannot be achieved, is also discussed. In this extended abstract, a simplistic example is proved to demonstrate the idea. In the final paper, test results on the high fidelity 6-DOF X-33 model with severe dispersions will be presented.

Zhu, J. J.↗

A Versatile Platform for Nanotechnology Based on Circular Permutation of a Chaperonin Protein

A number of protein complexes have been developed as nanoscale templates. These templates can be functionalized using the peptide sequences that bind inorganic materials. However, it is difficult to integrate peptides into a specific position within a protein template. Integrating intact proteins with desirable binding or catalytic activities is an even greater challenge. We present a general method for modifying protein templates using circular permutation so that additional peptide sequence can be added in a wide variety of specific locations. Circular permutation is a reordering of the polypeptide chain such that the original termini are joined and new termini are created elsewhere in the protein. New sequence can be joined to the protein termini without perturbing the protein structure and with minimal limitation on the size and conformation of the added sequence. We have used this approach to modify a chaperonin protein template, placing termini at five different locations distributed across the surface of the protein complex. These permutants are competent to form the double-ring structures typical of chaperonin proteins. The permuted double-rings also form the same assemblies as the unmodified protein. We fused a fluorescent protein to two representative permutants and demonstrated that it assumes its active structure and does not interfere with assembly of chaperonin double-rings.

Paavola, Chad↗

A spectral method determination of the first critical Rayleigh number for a low-Prandtl number crystal melt in a cylindrical container

The onset of laminar axisymmetric Rayleigh-Benard convection is investigated for a low-Prandtl number liquid metal in a cylindrical container. All surfaces are considered to be solid and no-slip. Two separate cases are examined for the thermal boundary conditions at the side wall, one with conducting and the other with insulated surface. The governing Boussinesq system is first perturbed and then simplified by introducing a Stokes stream function. Subsequently, a Chebyshev Galerkin spectral model is employed to reduce the simplified system to a system of first-order nonlinear ordinary differential equations. A local stability analysis determines the two values of the first critical Rayleigh number, Ra(sub cl), for the insulated and conducting side walls. As expected, the conducting Ra(sub cl) value of 2882.5 obtained from the present approach exceeded the corresponding insulated Ra(sub cl) value of 2331.6. For the insulated case, an earlier study using a different numerical approach suggests that Ra(sub cl) = 2261.9, while an experimental study measured Ra(sub cl) = 2700.

Dietz, C. M., Jr.↗

Computational methods for frictional contact with applications to the Space Shuttle orbiter nose-gear tire: Comparisons of experimental measurements and analytical predictions

A computational procedure is presented for the solution of frictional contact problems for aircraft tires. A Space Shuttle nose-gear tire is modeled using a two-dimensional laminated anisotropic shell theory which includes the effects of variations in material and geometric parameters, transverse-shear deformation, and geometric nonlinearities. Contact conditions are incorporated into the formulation by using a perturbed Lagrangian approach with the fundamental unknowns consisting of the stress resultants, the generalized displacements, and the Lagrange multipliers associated with both contact and friction conditions. The contact-friction algorithm is based on a modified Coulomb friction law. A modified two-field, mixed-variational principle is used to obtain elemental arrays. This modification consists of augmenting the functional of that principle by two terms: the Lagrange multiplier vector associated with normal and tangential node contact-load intensities and a regularization term that is quadratic in the Lagrange multiplier vector. These capabilities and computational features are incorporated into an in-house computer code. Experimental measurements were taken to define the response of the Space Shuttle nose-gear tire to inflation-pressure loads and to inflation-pressure loads combined with normal static loads against a rigid flat plate. These experimental results describe the meridional growth of the tire cross section caused by inflation loading, the static load-deflection characteristics of the tire, the geometry of the tire footprint under static loading conditions, and the normal and tangential load-intensity distributions in the tire footprint for the various static vertical-loading conditions. Numerical results were obtained for the Space Shuttle nose-gear tire subjected to inflation pressure loads and combined inflation pressure and contact loads against a rigid flat plate. The experimental measurements and the numerical results are compared.

Tanner, John A.↗

Global Ionosphere Perturbations Monitored by the Worldwide GPS Network

For the first time, measurements from the Global Positioning System (GPS) worldwide network are employed to study the global ionospheric total electron content(TEC) changes during a magnetic storm (November 26, 1994). These measurements are obtained from more than 60 world-wide GPS stations which continuously receive dual-frequency signals. Based on the delays of the signals, we have generated high resolution global ionospheric maps (GIM) of TEC at 15 minute intervals. Using a differential method comparing storm time maps with quiet time maps, we find that significant TEC increases (the positive effect ) are the major feature in the winter hemisphere during this storm (the maximum percent change relative to quiet times is about 150 percent).

GPS global positioning system TEC total electron c↗

Simulation of the TDRS multipath environment

Design principles and implementation methods are discussed for simulating the propagation path between a tracking and data relay satellite and a mission spacecraft. The emphasis is on multipath and Doppler simulation but additive disturbances are also considered. The recommended form of the simulator is fed separately with the unmodulated carrier, the unmodulated subcarriers (or spread-spectrum components) and the data signals. The perturbations are also introduced separately; then successive modulation operations are performed. The simulator is segmented into elements that perform the various functions of direct and specular multipath, diffuse fading, Doppler shift and delay spread. Delay spreads are realized by discrete delays operating on baseband signals. Doppler simulation and ionospheric or diffuse multipath fading are applied to individual paths before or after modulation of the carrier by delayed baseband signals. Block diagrams are presented on how the different elements are combined to create a complete channel simulator.

Baghdady, E. J.↗

Performance analysis of the multichannel astrometric photometer

It has been proposed that extrasolar planetary systems may be observed if perturbations in star position due to the orbit of Jupiter-type planets could be detected. To see this motion, high accuracy measurements of 0.01 milliarcsecond are required over a relatively large field of view. Techniques using a moving Ronchi grating have been proposed for this application and have been successful in ground-based lower resolution tests. The method may have application to other precision angular measurement problems. This paper explores the theoretical description of the method, considers certain of the error sources, and presents a preliminary calculation of the performance which may be achieved.

Huang, Chunsheng↗

High resolution frequency to time domain transformations applied to the stepped carrier MRIS measurements

Two narrow-band radar systems are developed for high resolution target range estimation in inhomogeneous media. They are reformulations of two presently existing systems such that high resolution target range estimates may be achieved despite the use of narrow bandwidth radar pulses. A double sideband suppressed carrier radar technique originally derived in 1962, and later abandoned due to its inability to accurately measure target range in the presence of an interfering reflection, is rederived to incorporate the presence of an interfering reflection. The new derivation shows that the interfering reflection causes a period perturbation in the measured phase response. A high resolution spectral estimation technique is used to extract the period of this perturbation leading to accurate target range estimates independent of the signal-to-interference ratio. A non-linear optimal signal processing algorithm is derived for a frequency-stepped continuous wave radar system. The resolution enhancement offered by optimal signal processing of the data over the conventional Fourier Transform technique is clearly demonstrated using measured radar data. A method for modeling plane wave propagation in inhomogeneous media based on transmission line theory is derived and studied. Several simulation results including measurement of non-uniform electron plasma densities that develop near the heat tiles of a space re-entry vehicle are presented which verify the validity of the model.

Ardalan, Sasan H.↗

Convection in Ice I with Composite Newtonian/Non-Newtonian Rheology: Application to the Icy Galilean Satellites

Ice I exhibits a complex rheology at temperature and pressure conditions appropriate for the interiors of the outer ice I shells of Europa, Ganymede, and Callisto. We use numerical methods to determine the conditions required to trigger convection in an ice I shell with a stress-, temperature-, and grain-size-dependent rheology measured in laboratory experiments by Goldsby and Kohlstedt [2001] (henceforth GK2001). Triggering convection from an initially conductive ice shell with a non-Newtonian rheology for ice I requires that a finite-amplitude temperature perturbation be issued to the ice shell [2]. Here, we characterize the amplitude and wavelength of temperature perturbation required to initiate convection in the outer ice I shells of Europa, Ganymede, and Callisto using the GK2001 rheology for a range of ice grain sizes.

Barr, Amy C.↗

An Empirical Method for Determining the Lunar Gravity Field

A method has been devised to determine the spherical harmonic coefficients of the lunar gravity field. This method consists of a two-step data reduction and estimation process. In the first step, a weighted least-squares empirical orbit determination scheme is applied to Doppler tracking data from lunar orbits to estimate long-period Kepler elements and rates. Each of the Kepler elements is represented by an independent function of time. The long-period perturbing effects of the earth, sun, and solar radiation are explicitly modeled in this scheme. Kepler element variations estimated by this empirical processor are ascribed to the non-central lunar gravitation features. Doppler data are reduced in this manner for as many orbits as are available. In the second step, the Kepler element rates are used as input to a second least-squares processor that estimates lunar gravity coefficients using the long-period Lagrange perturbation equations.

Ferrari, A. J.↗

Lunar gravity derived from long-period satellite motion - A proposed method.

A new method has been devised to determine the spherical harmonic coefficients of the lunar gravity field. This method consists of a two-step data reduction and estimation process. In the first step, a weighted least-squares empirical orbit determination scheme is applied to Doppler tracking data from lunar orbits to estimate long-period Kepler elements and rates. Each of the Kepler elements is represented by an independent function of time. The long-period perturbing effects of the earth, sun, and solar radiation are explicitly modeled in this scheme. Kepler element variations estimated by this empirical processor are then ascribed to the non-central lunar gravitation features. Doppler data are reduced in this manner for as many orbits as are available. In the second step, the Kepler element rates are used as input to a second least-squares processor that estimates lunar gravity coefficients using the long-period Lagrange perturbation equations.

Ferrari, A. J.↗

Performance analysis of flexible aircraft with active control

Small perturbation equations of motion of a flexible aircraft with an active control technology (ACT) system were developed to evaluate the stability and performance of the controlled aircraft. The total aircraft system was formulated in state vector format and the system of equations was completed with fully unsteady and low frequency aerodynamics for arbitrary, complex configurations based on a potential aerodynamic method. The ACT system equations were incorporated in the digital computer program FCAP (Flight Control Analysis Program) which can be used for the analysis of complete aircraft configurations, including control system, with either low frequency or fully unsteady aerodynamics. The application of classical performance analyses including frequency response, poles and zeros, mean square response, and time response in FCAP in state vector format was discussed.

Noll, R. B.↗

Aerodynamic parameter estimation via Fourier modulating function techniques

Parameter estimation algorithms are developed in the frequency domain for systems modeled by input/output ordinary differential equations. The approach is based on Shinbrot's method of moment functionals utilizing Fourier based modulating functions. Assuming white measurement noises for linear multivariable system models, an adaptive weighted least squares algorithm is developed which approximates a maximum likelihood estimate and cannot be biased by unknown initial or boundary conditions in the data owing to a special property attending Shinbrot-type modulating functions. Application is made to perturbation equation modeling of the longitudinal and lateral dynamics of a high performance aircraft using flight-test data. Comparative studies are included which demonstrate potential advantages of the algorithm relative to some well established techniques for parameter identification. Deterministic least squares extensions of the approach are made to the frequency transfer function identification problem for linear systems and to the parameter identification problem for a class of nonlinear-time-varying differential system models.

Pearson, A. E.↗

Towards Robust Designs Via Multiple-Objective Optimization Methods

Fabricating and operating complex systems involves dealing with uncertainty in the relevant variables. In the case of aircraft, flow conditions are subject to change during operation. Efficiency and engine noise may be different from the expected values because of manufacturing tolerances and normal wear and tear. Engine components may have a shorter life than expected because of manufacturing tolerances. In spite of the important effect of operating- and manufacturing-uncertainty on the performance and expected life of the component or system, traditional aerodynamic shape optimization has focused on obtaining the best design given a set of deterministic flow conditions. Clearly it is important to both maintain near-optimal performance levels at off-design operating conditions, and, ensure that performance does not degrade appreciably when the component shape differs from the optimal shape due to manufacturing tolerances and normal wear and tear. These requirements naturally lead to the idea of robust optimal design wherein the concept of robustness to various perturbations is built into the design optimization procedure. The basic ideas involved in robust optimal design will be included in this lecture. The imposition of the additional requirement of robustness results in a multiple-objective optimization problem requiring appropriate solution procedures. Typically the costs associated with multiple-objective optimization are substantial. Therefore efficient multiple-objective optimization procedures are crucial to the rapid deployment of the principles of robust design in industry. Hence the companion set of lecture notes (Single- and Multiple-Objective Optimization with Differential Evolution and Neural Networks ) deals with methodology for solving multiple-objective Optimization problems efficiently, reliably and with little user intervention. Applications of the methodologies presented in the companion lecture to robust design will be included here. The evolutionary method (DE) is first used to solve a relatively difficult problem in extended surface heat transfer wherein optimal fin geometries are obtained for different safe operating base temperatures. The objective of maximizing the safe operating base temperature range is in direct conflict with the objective of maximizing fin heat transfer. This problem is a good example of achieving robustness in the context of changing operating conditions. The evolutionary method is then used to design a turbine airfoil; the two objectives being reduced sensitivity of the pressure distribution to small changes in the airfoil shape and the maximization of the trailing edge wedge angle with the consequent increase in airfoil thickness and strength. This is a relevant example of achieving robustness to manufacturing tolerances and wear and tear in the presence of other objectives.

Man Mohan, Rai↗

Response analysis of flexible aircraft with active control

The small-perturbation equations of motion of a flexible aircraft with an active control technology (ACT) system are developed to evaluate the stability and performance of the controlled aircraft. The total aircraft system is formulated in state vector format and the system of equations is completed with fully unsteady and low-frequency aerodynamics for arbitrary, complex configurations based on a potential aerodynamic method. The ACT system equations are incorporated in the digital computer program FCAP (Flight Control Analysis Program) which can be used for the analysis of complete aircraft configurations, including control system, with either low-frequency or fully unsteady aerodynamics. The application of classical performance analyses including frequency response, poles and zeros, mean-square response, and time response in FCAP in state vector format is discussed.

Noll, R. B.↗

Aerodynamic sound generation due to vortex-aerofoil interaction. Part 2: Analysis of the acoustic field

The Lighthill method was the basic procedure used to analyze the sound field associated with a vortex of modified strength interacting with an airfoil. A free vortex interacting with an airfoil in uniform motion was modeled in order to determine the sound field due to all the acoustic sources, not only on the airfoil surfaces (dipoles), but also the ones distributed on the perturbed flow field (quadrupoles) due to the vortex-airfoil interaction. Because inviscid flow is assumed in the study of the interaction, the quadrupoles considered in the perturbed flow field are entirely due to an unsteady flow field. The effects of airfoil thickness on the second radiation are examined by using a symmetric Joukowski airfoil for the vortex-airfoil interaction. Sound radiation in a plane, far field simplification, and computation of the sound field are discussed.

Parasarathy, R.↗

Calculation of sensitivity derivatives in thermal problems by finite differences

The optimum design of a structure subject to temperature constraints is considered. When mathematical optimization techniques are used, derivatives of the temperature constraints with respect to the design variables are usually required. In the case of large aerospace structures, such as the Space Shuttle, the computation of these derivatives can become prohibitively expensive. Analytical methods and a finite difference approach have been considered in studies conducted to improve the efficiency of the calculation of the derivatives. The present investigation explores two possibilities for enhancing the effectiveness of the finite difference approach. One procedure involves the simultaneous solution of temperatures and derivatives. The second procedure makes use of the optimum selection of the magnitude of the perturbations of the design variables to achieve maximum accuracy.

Haftka, R. T.↗

Do stars that lose mass expand or contract? A semianalytical approach

Semianalytic methods are used to determine whether a star of polytropic index n expands or contracts following an insentropic loss of mass. It is found that most stars shrink, including radiative stars which can be approximated by polytropes of index roughly three. However, stars that have a ratio of specific heats gamma only slightly larger than 1 + 1/n do expand if the perturbation of the core density is sufficiently large. For stars that shrink, the fractional change in the radius is not a linear function of the fractional change in mass. Using Lane's theorem, it is shown that a star with gamma = 4/3 cannot lose mass isentropically, although the radius increases if the core density is reduced.

Heisler, J.↗