SEARCH · Search NASA
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.
Extension of problem of singular perturbation for linear scalar constant coefficient differential- difference equation with single retardation to several retardations, noting degenerate equation solution
Singular perturbation techniques are used in this paper to develop a two time scale procedure for designing static gain output feedback controllers. It is shown that, if certain control spillover conditions are satisfied, control designs based on reduced-order models will stabilize the corrgsponding full system. Optimal output feedback control theory is used to derive the necessary conditions for the stabilizing gain matrix. The problem of stabilizing a model for a large space structure is used to illustrate the practicality of the approach.
A theoretical model is developed for the sound generated when a convected vortical or entropic gust encounters an airfoil at nonzero angle of attack. The theory is based on a linearization of the Euler equations about the steady, subsonic flow past the airfoil. High-frequency gusts, whose wavelengths are short compared to the airfoil chord, but long compared to the displacement of the mean-flow stagnation point from the leading edge, are considered. The analysis utilizes singular perturbation techniques and involves four asymptotic regions. Local regions, which scale on the gust wavelength, are present at the airfoil leading and trailing edges. Behind the airfoil a "transition" region, which is similar to the transition zone between illuminated and shadow zones in optical problems, is present. In the outer region, far away from the airfoil edges and wake, the solution has a geometric-acoustics form. The primary sound generation is found to be concentrated in the local leading-edge region. The trailing edge plays a secondary role as a scatterer of the sound generated in the leading-edge region. Parametric calculations are presented which illustrate that moderate levels of airfoil steady loading can significantly affect the sound field produced by airfoil-gust interactions.
Note is made of an earlier paper which defined finite difference operators for the Hilbert space L2(m), and gave the eigenvalues for these operators. The present work examines eigenvalues for higher order singular differential operators by using finite difference methods. The two self-adjoint operators investigated are defined by a particular value in the same Hilbert space, L2(m), and are strictly positive with compact inverses. A class of finite difference operators is considered, with the idea of application to the theory of Toeplitz matrices. The approximating operators consist of a good approximation plus a perturbing operator.
The theoretical development and the comparison of results with data of a thick wing and pylon-fuselage-fanpod-nacelle analysis are presented. The analysis utilizes potential flow theory to compute the surface velocities and pressures, section lift and center of pressure, and the total configuration lift, moment, and vortex drag. The skin friction drag is also estimated in the analysis. The perturbation velocities induced by the wing and pylon, fuselage and fanpod, and nacelle are represented by source and vortex lattices, quadrilateral vortices, and source frustums, respectively. The strengths of these singularities are solved for simultaneously including all interference effects. The wing and pylon planforms, twists, cambers, and thickness distributions, and the fuselage and fanpod geometries can be arbitrary in shape, provided the surface gradients are smooth. The flow through nacelle is assumed to be axisymmetric. An axisymmetric center engine hub can also be included. The pylon and nacelle can be attached to the wing, fuselage, or fanpod.
Bridgman crystal growth can be conducted in the so-called "detached" solidification regime, where the growing crystal is detached from the crucible wall. A small gap between the growing crystal and the crucible wall, of the order of 100 micrometers or less, can be maintained during the process. A meniscus is formed at the bottom of the melt between the crystal and crucible wall. Under proper conditions, growth can proceed without collapsing the meniscus. The meniscus shape plays a key role in stabilizing the process. Thermal and other process parameters can also affect the geometrical steady-state stability conditions of solidification. The dynamic stability theory of the shaped crystal growth process has been developed by Tatarchenko. It consists of finding a simplified autonomous set of differential equations for the radius, height, and possibly other process parameters. The problem then reduces to analyzing a system of first order linear differential equations for stability. Here we apply a modified version of this theory for a particular case of detached solidification. Approximate analytical formulas as well as accurate numerical values for the capillary stability coefficients are presented. They display an unexpected singularity as a function of pressure differential. A novel approach to study the thermal field effects on the crystal shape stability has been proposed. In essence, it rectifies the unphysical assumption of the model that utilizes a perturbation of the crystal radius along the axis as being instantaneous. It consists of introducing time delay effects into the mathematical description and leads, in general, to stability over a broader parameter range. We believe that this novel treatment can be advantageously implemented in stability analyses of other crystal growth techniques such as Czochralski and float zone methods.
A previous small perturbation analysis of circumferential waves in circumferential compression systems, assuming inviscid flow, is shown to be consistent with observations that narrow diffusers are more stable than wide ones, when boundary layer displacement effect is included. The Moore-Greitzer analysis for finite strength transients containing both surge and rotating stall in axial machines is adapted for a centrifugal compression system. Under certain assumptions, and except for a new second order swirl, the diffuser velocity field, including resonant singularities, can be carried over from the previous inviscid linear analysis. Nonlinear transient equations are derived and applied in a simple example to show that throttling through a resonant value of flow coefficient must occur in a sudden surge-like drop, accompanied by a transient rotating wave. This inner solution is superseded by an outer surge response on a longer time scale. Surge may occur purely as result of circumferential wave resonance. Numerical results are shown for various parametric choices relating to throttle schedule and the characteristic slope. A number of circumferential modes considered simultaneously is briefly discussed.
This paper investigates the effect of airfoil steady loading on the sound generated by the interaction of an airfoil with a convected disturbance. A previous theory, which included only the incidence angle contribution to the mean loading, is extended to include camber. The theory is based on a linearization of the Euler equations about a nonuniform, 0(1) Mach number subsonic mean flow. The discussion concentrates on the case of a slightly cambered airfoil at small incidence angle, interacting with a gust whose wavelength is short compared to the airfoil chord. The small parameter representing the amount of camber and incidence, and the large parameter representing the ratio of airfoil chord to disturbance wavelength, are utilized in a singular perturbation solution to the governing equations. Acoustic power calculations reveal that the amount of sound generated increases significantly with increased loading. More importantly, it is shown that the radiated acoustic power correlates very well with the strength of the mean flow around the leading edge.
In fluid theory, the ordering of low-frequency modes in a homogeneous plasma is based on the phase velocity, since modes do not intersect each other in dispersion diagrams as a function of wavenumber or other parameters. In linear kinetic theory, modes cross each other. Thus a consistent and useful classification should be based on the physical properties of the modes instead. This paper attempts such a classification by documeting the dispersion and general mode properties of the low-frequency waves (omega much less than (OMEGA(sub ci) OMEGA(sub ce) (exp 1/2)), where OMEGA(sub ci), OMEGA(sub ce) are the cyclotron frequencies of the ions and electrons, respectively) in kinetic theory, and by comparing them to the results of two-fluid theory. Kinetic theory gives a seperate Alfven/ion-cyclotron (A/IC) wave with phase speed Omega/k approximately = v(sub A) cos theta for omega much less than OMEGA(sub ci), where v(sub A) is the Alfven velocity and theta the angle of propagation between wave vector k and background magnetic field B(sub o). For a given wavenumber, the magnetosonic mode is a double-valued solution with a singular point in theta, beta parameter space, where beta is the ratio of thermal pressure to magnetic pressure. It is shown that a branch cut starting at the singular point theta approximately 30 deg, beta approximately 3 and leading to larger beta gives a practical and consitent seperation of this double-valued magnetosonic solution. Selection of this branch cut results in a moderately damped fast/magnetos onic and a heavily damped slow/sound wave. A comprehensive review of the polarization, compressibility and other mode properties is given and shown to be consistent with the selected branch cut. At small wavenumbers, the kinetic mode properties typically start to deviate significantly from their fluid counterparts at beta approximately 0.5. At larger beta, there is no longer a consistent correspondence between the fluid and kinetic modes. Kinetic theory also dictates the use of different mode properties to distinguish between them in observational data. For example, the phase between the density and magnetic field perturbation may become useless at high beta, whereas the direction of the magnetic field perturbations with respect to k and B(sub o) remains a useful characteristic. Two quantities based on this characteristic are suggested and are shown to be useful also to distinguish between the mirror mode and A/IC waves in a plasma with temperature anisotropy.
The general theory of potential aerodynamic flow around a lifting body having arbitrary shape and motion is presented. By using the Green function method, an integral representation for the potential is obtained for both supersonic and subsonic flow. Under small perturbation assumption, the potential at any point, P, in the field depends only upon the values of the potential and its normal derivative on the surface, sigma, of the body. Hence, if the point P approaches the surface of the body, the representation reduces to an integro-differential equation relating the potential and its normal derivative (which is known from the boundary conditions) on the surface sigma. For the important practical case of small harmonic oscillation around a rest position, the equation reduces to a two-dimensional Fredholm integral equation of second-type. It is shown that this equation reduces properly to the lifting surface theories as well as other classical mathematical formulas. The question of uniqueness is examined and it is shown that, for thin wings, the operator becomes singular as the thickness approaches zero. This fact may yield numerical problems for very thin wings.
The general method for analyzing steady subsonic potential aerodynamic flow around a lifting body having arbitrary shape is presented. By using the Green function method, an integral representation for the potential is obtained. Under small perturbation assumption, the potential at any point, P, in the field depends only upon the values of the potential and its normal derivative on the surface of the body. Hence if the point P approaches the surface of the body, the representation reduces to an integral equation relating the potential and its normal derivative (which is known from the boundary conditions) on the surface. The question of uniqueness is examined and it is shown that, for thin wings, the operator becomes singular as the thickness approaches zero. This fact may yield numerical problems for very thin wings. However, numerical results obtained for a rectangular wing in subsonic flow show that these problems do not appear even for thickness ratio tau = .001. Comparison with existing results shows that the proposed method is at least as fast and accurate as the lifting surface theories.
Similarity solutions were found which give the adiabatic flow of an ideal gas about two-dimensional and axisymmetric power-law bodies at infinite Mach number to second order in the body slenderness parameter. The flow variables were expressed as a sum of zero-order and perturbation similarity functions for which the axial variations in the flow equations separated out. The resulting similarity equations were integrated numerically. The solutions, which are universal functions, are presented in graphic and tabular form. To avoid a singularity in the calculations, the results are limited to body power-law exponents greater than about 0.85 for the two-dimensional case and 0.75 for the axisymmetric case. Because of the entropy layer induced by the nose bluntness (for power-law bodies other than cones and wedges), only the pressure function is valid at the body surface. The similarity results give excellent agreement with the exact solutions for inviscid flow over wedges and cones having half-angles up to about 20 deg. They give good agreement with experimental shock-wave shapes and surface-pressure distributions for 3/4-power axisymmetric bodies, considering that Mach number and boundary-layer displacement effects are not included in the theory.
During this research period we have made significant progress in the four proposed areas: (1) design of robust controllers via H infinity optimization; (2) design of robust controllers via mixed H2/H infinity optimization; (3) M-delta structure and robust stability analysis for structured uncertainties; and (4) a study on controllability and observability of perturbed plant. It is well known now that the two-Riccati-equation solution to the H infinity control problem can be used to characterize all possible stabilizing optimal or suboptimal H infinity controllers if the optimal H infinity norm or gamma, an upper bound of a suboptimal H infinity norm, is given. In this research, we discovered some useful properties of these H infinity Riccati solutions. Among them, the most prominent one is that the spectral radius of the product of these two Riccati solutions is a continuous, nonincreasing, convex function of gamma in the domain of interest. Based on these properties, quadratically convergent algorithms are developed to compute the optimal H infinity norm. We also set up a detailed procedure for applying the H infinity theory to robust control systems design. The desire to design controllers with H infinity robustness but H(exp 2) performance has recently resulted in mixed H(exp 2) and H infinity control problem formulation. The mixed H(exp 2)/H infinity problem have drawn the attention of many investigators. However, solution is only available for special cases of this problem. We formulated a relatively realistic control problem with H(exp 2) performance index and H infinity robustness constraint into a more general mixed H(exp 2)/H infinity problem. No optimal solution yet is available for this more general mixed H(exp 2)/H infinity problem. Although the optimal solution for this mixed H(exp 2)/H infinity control has not yet been found, we proposed a design approach which can be used through proper choice of the available design parameters to influence both robustness and performance. For a large class of linear time-invariant systems with real parametric perturbations, the coefficient vector of the characteristic polynomial is a multilinear function of the real parameter vector. Based on this multilinear mapping relationship together with the recent developments for polytopic polynomials and parameter domain partition technique, we proposed an iterative algorithm for coupling the real structured singular value.