Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Singularity”

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 685 records · Page 38

Hypersonic vehicle control law development using H infinity and mu-synthesis

Applicability and effectivness of robust control techniques to a single-stage-to-orbit (SSTO) airbreathing hypersonic vehicle on an ascent accelerating path and their effectivness are explored in this paper. An SSTO control system design problem, requiring high accuracy tracking of velocity and altitude commands while limiting an angle of attack oscillations, minimizing control power usage and stabilizing the vehicle all in the presence of atmospheric turbulence and uncertainty in the system, was formulated to compare results of the control designs using H infinity and mu-synthesis procedures. The math model, an integrated flight/propulsion dynamic model of a conical accelerator class vehicle, was linearized as the vehicle accelerated through Mach 8. Controller analysis was conducted using the singular value technique and the mu-analysis approach. Analysis results were obtained in both the frequency and the time domains. The results clearly demonstrate the inherent advantages of the structured singular value framework for this class of problems. Since payload performance margins are so critical for the SSTO mission, it is crucial that adequate stability margins be provided without sacrificing any payload mass.

Gregory, Irene M.↗

Non-iterative computation of the infimum in H-infinity-optimization for plants with invariant zeros on the j-omega axis

A simple and noniterative procedure for the computation of the exact value of the infimum in the singular H-infinity optimization problem is presented. The problem formulation is general, and no restrictions are placed on the finite and infinite zero structures of the system and on the direct feedthrough terms between the control input and the controlled output variables and between the disturbance input and the measurement output variables. The method is applicable to a class of singular H-infinity optimization problems for which the transfer functions from the control input to the controlled output and from the disturbance input to the measurement output satisfy certain geometric conditions. In particular, the authors extend the result of B. M. Chen et al. (1991) by allowing these two transfer functions to have invariant zeros on the j-omega axis.

Chen, Ben M.↗

Attitude control with realization of linear error dynamics

An attitude control law is derived to realize linear unforced error dynamics with the attitude error defined in terms of rotation group algebra (rather than vector algebra). Euler parameters are used in the rotational dynamics model because they are globally nonsingular, but only the minimal three Euler parameters are used in the error dynamics model because they have no nonlinear mathematical constraints to prevent the realization of linear error dynamics. The control law is singular only when the attitude error angle is exactly pi rad about any eigenaxis, and a simple intuitive modification at the singularity allows the control law to be used globally. The forced error dynamics are nonlinear but stable. Numerical simulation tests show that the control law performs robustly for both initial attitude acquisition and attitude control.

Paielli, Russell A.↗

Nonlinear vortex trail dynamics. II - Analytic solutions

Spatially periodic large amplitude solutions of the von Karman model are obtained in the neighborhood of singularities. These singularities correspond to vortex clusters in the physical plane. The quasi-periodic and unbounded solutions found analytically confirm earlier numerical work and show qualitative agreement with experimental observations of large-scale phenomena of vortex trails. Separatrices or heteroclinic orbits were explicitly found for an integrable approximate equation, which indicate that the von Karman model itself supports chaotic solutions.

Lim, Chjan C.↗

A note on the distortion of a flat-plate boundary layer by free-stream vorticity normal to the plate

The purpose of this note is to construct a local solution that eliminates a residual velocity discontinuity in the inviscid portion of a solution obtained in a recent paper by Goldstein, Leib and Cowley (1992). This result is of importance because it shows that the solution obtained in that paper is entirely non-singular outside the viscous wall boundary layer and that any singularity in the problem will have to arise in the usual way through a breakdown in the viscous boundary layer.

Goldstein, M. E.↗

Mutual coupling between parallel columns of periodic slots in a ground plane surrounded by dielectric slabs

Arrays of slots with arbitrary orientation in a conducting plane which are infinitely periodic in one dimension and finitely periodic in another dimension are considered. The perfect conducting plane extends infinitely and is bounded on each side by dielectric slabs of finite thickness and infinite extent. Single columns of slots are represented by equivalent magnetic scattering currents, which are solved for via the moment method. The mutual coupling (admittance) between slot columns in the presence of the stratified media is found by the array scanning method (ASM), which expresses the admittance as the average of the scan admittance of an artificially constructed doubly infinite array of slots over all real scan angles. The technique avoids the use of Sommerfeld integrals, but still gives rise to singularities at scan angles corresponding to the resonant excitation of surface waves. An analytical approximation is made to remove these surface wave singularities, thus making numerical implementation of the method practical. The method is quite general and may be extended to thin slot elements of arbitrary shape and orientation.

Skinner, J. P.↗

Three-dimensional kinematic reconnection of plasmoids with nulls

The global nonlinear dynamics of magnetic field lines in plasmoids with a pair of nulls, where B = 0, is studied. The aim of this analysis is to describe the separatrix surfaces on which singularities can occur in ideal magnetohydrodynamics because of topological changes in the field. These separatrix surfaces should locate the boundary layers associated with 3D reconnection in the presence of resistivity or inertia. It is found that the field lines exhibit chaotic scattering with several properties in common with plasmoid models without nulls (in which one component of the magnetic field never changes sign). In particular, the singular surfaces can be fractal, implying complex current density structures down to the dissipation scale. These generic features are expected to exist in typical coronal magnetic geometries exhibiting three-dimensional reconnection and the formation of current sheets.

Lau, Yun-Tung↗

Confining potential in momentum space

A method is presented for the solution in momentum space of the bound state problem with a linear potential in r space. The potential is unbounded at large r leading to a singularity at small q. The singularity is integrable, when regulated by exponentially screening the r-space potential, and is removed by a subtraction technique. The limit of zero screening is taken analytically, and the numerical solution of the subtracted integral equation gives eigenvalues and wave functions in good agreement with position space calculations.

Norbury, John W.↗

Asteroidal collision probabilities

Several past calculations of collision probabilities between pairs of bodies on independent orbits have yielded inconsistent results. We review the methodologies and identify their various problems. Greenberg's (1982) collision probability formalism (now with a corrected symmetry assumption) is equivalent to Wetherill's (1967) approach, except that it includes a way to avoid singularities near apsides. That method shows that the procedure by Namiki and Binzel (1991) was accurate for those cases where singularities did not arise.

Bottke, William F., Jr.↗

Matched asymptotic expansion of the Hamilton-Jacobi-Bellman equation for aeroassisted plane-change maneuvers

In this paper we develop a general procedure for constructing a matched asymptotic expansion of the Hamilton-Jacobi-Bellman equation based on the method of characteristics. The development is for a class of perturbation problems whose solution exhibits two-time-scale behavior. A regular expansion for problems of this type is inappropriate since it is not uniformly valid over a narrow range of the independent variable. Of particular interest here is the manner in which matching and boundary conditions are enforced when the expansion is carried out to first order. Two cases are distinguished - one where the left boundary condition coincides with or lies to the right of the singular region and one where the left boundary condition lies to the left of the singular region. A simple example is used to illustrate the procedure, and its potential application to aeroassisted plane change is described.

Calise, Anthony J.↗

A comparative study of multivariable robustness analysis methods as applied to integrated flight and propulsion control

Three multivariable robustness analysis methods are compared and contrasted. The focus of the analysis is on system stability and performance robustness to uncertainty in the coupling dynamics between two interacting subsystems. Of particular interest is interacting airframe and engine subsystems, and an example airframe/engine vehicle configuration is utilized in the demonstration of these approaches. The singular value (SV) and structured singular value (SSV) analysis methods are compared to a method especially well suited for analysis of robustness to uncertainties in subsystem interactions. This approach is referred to here as the interacting subsystem (IS) analysis method. This method has been used previously to analyze airframe/engine systems, emphasizing the study of stability robustness. However, performance robustness is also investigated here, and a new measure of allowable uncertainty for acceptable performance robustness is introduced. The IS methodology does not require plant uncertainty models to measure the robustness of the system, and is shown to yield valuable information regarding the effects of subsystem interactions. In contrast, the SV and SSV methods allow for the evaluation of the robustness of the system to particular models of uncertainty, and do not directly indicate how the airframe (engine) subsystem interacts with the engine (airframe) subsystem.

Schierman, John D.↗

Statistics of arcs in clusters of galaxies

Samples of gravitational lens events in clusters show many large arcs compared to arclets, relative to what can be obtained by idealized singular lens models. We describe the probability of image magnification for point sources and for simple but more realistic gravitational lensing models that include a finite core size and an ellipticity. In addition, we explore the changes in the probability distribution of image magnifications, distortions, and angular extents for sources of different sizes as the parameters of the lenses are varied. A finite core in spherically symmetric lens models introduces a discontinuity in the probability distribution at which the relative number of highly magnified images is increased. In elliptical lenses, this discontinuity and its effect are replaced by a continuous increase in the probability of obtaining high-magnification images relative to singular spherically symmetric models. We also find that the finite size of the source causes a further increase in the expected number of images just below the maximum possible magnification.

Bergmann, Anton G.↗

Symmetric Galerkin boundary formulations employing curved elements

Accounts of the symmetric Galerkin approach to boundary element analysis (BEA) have recently been published. This paper attempts to add to the understanding of this method by addressing a series of fundamental issues associated with its potential computational efficiency. A new symmetric Galerkin theoretical formulation for both the (harmonic) heat conduction and the (biharmonic) elasticity problem that employs regularized singular and hypersingular boundary integral equations (BIEs) is presented. The novel use of regularized BIEs in the Galerkin context is shown to allow straightforward incorporation of curved, isoparametric elements. A symmetric reusable intrinsic sample point (RISP) numerical integration algorithm is shown to produce a Galerkin (i.e., double) integration strategy that is competitive with its counterpart (i.e., singular) integration procedure in the collocation BEA approach when the time saved in the symmetric equation solution phase is also taken into account. This new formulation is shown to be capable of employing hypersingular BIEs while obviating the requirement of C 1 continuity, a fact that allows the employment of the popular continuous element technology. The behavior of the symmetric Galerkin BEA method with regard to both direct and iterative equation solution operations is also addressed. A series of example problems are presented to quantify the performance of this symmetric approach, relative to the more conventional unsymmetric BEA, in terms of both accuracy and efficiency. It is concluded that appropriate implementations of the symmetric Galerkin approach to BEA indeed have the potential to be competitive with, if not superior to, collocation-based BEA, for large-scale problems.

Kane, J. H.↗

Dynamic structural correlation via nonlinear programming techniques

A solution to the correlation between structural dynamic test results and finite element analyses of the same components is presented in this paper. Basically, the method can be categorized as a Levenberg-Marquardt type Gauss-Newton method which requires only the differences between FE modal analyses and test results and their first derivatives with respect to preassigned design variables. With proper variable normalization and equation scaling, the method has been made numerically better-conditioned and the inclusion of the Levenberg-Marquardt technique overcomes any remaining difficulty encountered in inverting singular or near-singular matrices. An important feature is that each iteration requires only one function evaluation along with the associated design sensitivity analysis and so the procedure is computationally efficient.

Ting, T.↗

Integrals and integral equations in linearized wing theory

The formulas of subsonic and supersonic wing theory for source, doublet, and vortex distributions are reviewed and a systematic presentation is provided which relates these distributions to the pressure and to the vertical induced velocity in the plane of the wing. It is shown that care must be used in treating the singularities involved in the analysis and that the order of integration is not always reversible. Concepts suggested by the irreversibility of order of integration are shown to be useful in the inversion of singular integral equations when operational techniques are used. A number of examples are given to illustrate the methods presented, attention being directed to supersonic flight speed.

Lomax, Harvard↗

A special method for finding body distortions that reduce the wave drag of wing and body combinations at supersonic speeds

For a given wing and supersonic Mach number, the problem of shaping an adjoining fuselage so that the combination will have a low wave drag is considered. Only fuselages that can be simulated by singularities (multipoles) distributed along the body axis are studied. However, the optimum variations of such singularities are completely specified in terms of the given wing geometry. An application is made to an elliptic wing having a biconvex section, a thickness-chord ratio equal to 0.05 at the root, and an aspect ratio equal to 3. A comparison of the theoretical results with a wind-tunnel experiment is also presented.

Lomax, Harvard↗

A wall interference assessment/correction system

A Wall Signature method originally developed by Hackett was selected to be adapted for the Ames 12-ft Wind Tunnel Wall Interference Assessment/Correction (WIAC) System in the project. This method uses limited measurements of the static pressure at the wall, in conjunction with the solid wall boundary condition, to determine the strength and distribution of singularities representing the test article. The singularities are used in turn for estimating wall interference at the model location. The lifting interference will be treated separately by representing in a horseshoe vortex system for the model's lifting effects. The development and implementation of a working prototype will be completed, delivered, and documented with a software manual. The WIAC code will be validated by conducting numerically simulated experiments rather than actual wind tunnel experiments. The simulations will be used to generate both free-air and confined wind-tunnel flow fields for each of the test articles over a range of test configurations. Specifically, the pressure signature at the test section wall will be computed for the tunnel case to provide the simulated 'measured' data. These data will serve as the input for the WIAC method - Wall Signature method. The performance of the WIAC method then may be evaluated by comparing the corrected parameters with those for the free-air simulation. The following two additional tasks are included in the supplement No. 1 to the basic Grant. On-line wall interference calculation: The developed wall signature method (modified Hackett's method) for Ames 12-ft Tunnel will be the pre-computed coefficients which facilitate the on-line calculation of wall interference; and support system effects estimation: The effects on the wall pressure measurements due to the presence of the model support systems will be evaluated.

Lo, Ching F.↗

Guidance law development for aeroassisted transfer vehicles using matched asymptotic expansions

This report addresses and clarifies a number of issues related to the Matched Asymptotic Expansion (MAE) analysis of skip trajectories, or any class of problems that give rise to inner layers that are not associated directly with satisfying boundary conditions. The procedure for matching inner and outer solutions, and using the composite solution to satisfy boundary conditions is developed and rigorously followed to obtain a set of algebraic equations for the problem of inclination change with minimum energy loss. A detailed evaluation of the zeroth order guidance algorithm for aeroassisted orbit transfer is performed. It is shown that by exploiting the structure of the MAE solution procedure, the original problem, which requires the solution of a set of 20 implicit algebraic equations, can be reduced to a problem of 6 implicit equations in 6 unknowns. A solution that is near optimal, requires a minimum of computation, and thus can be implemented in real time and on-board the vehicle, has been obtained. Guidance law implementation entails treating the current state as a new initial state and repetitively solving the zeroth order MAE problem to obtain the feedback controls. Finally, a general procedure is developed for constructing a MAE solution up to first order, of the Hamilton-Jacobi-Bellman equation based on the method of characteristics. The development is valid for a class of perturbation problems whose solution exhibits two-time-scale behavior. A regular expansion for problems of this type is shown to be inappropriate since it is not valid over a narrow range of the independent variable. That is, it is not uniformly valid. Of particular interest here is the manner in which matching and boundary conditions are enforced when the expansion is carried out to first order. Two cases are distinguished-one where the left boundary condition coincides with, or lies to the right of, the singular region, and another one where the left boundary condition lies to the left of the singular region. A simple example is used to illustrate the procedure where the obtained solution is uniformly valid to O(Epsilon(exp 2)). The potential application of this procedure to aeroassisted plane change is also described and partially evaluated.

Calise, Anthony J.↗