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.
Sufficient conditions for global stability of a buck regulator using a discrete control law are found. The method of paired systems and Liapunov functions are used to establish global stability and to study the convergence of the regulator. A heuristic argument is given that the optimal switching curves associated with the paired continuous systems approximate the optimal switching curves of the discrete systems.
In this paper we review the existing and develop new continuous Galerkin methods for solving time dependent partial differential equations with higher order derivatives in one and multiple space dimensions. We review local discontinuous Galerkin methods for convection diffusion equations involving second derivatives and for KdV type equations involving third derivatives. We then develop new local discontinuous Galerkin methods for the time dependent bi-harmonic type equations involving fourth derivatives, and partial differential equations involving fifth derivatives. For these new methods we present correct interface numerical fluxes and prove L(exp 2) stability for general nonlinear problems. Preliminary numerical examples are shown to illustrate these methods. Finally, we present new results on a post-processing technique, originally designed for methods with good negative-order error estimates, on the local discontinuous Galerkin methods applied to equations with higher derivatives. Numerical experiments show that this technique works as well for the new higher derivative cases, in effectively doubling the rate of convergence with negligible additional computational cost, for linear as well as some nonlinear problems, with a local uniform mesh.
Many space-science experiments need an active isolation system to provide a sufficiently quiescent microgravity environment. Modern control methods provide the potential for both high-performance and robust stability in the presence of parametric uncertainties that are characteristic of microgravity vibration isolation systems. While H2 and H(infinity) methods are well established, neither provides the levels of attenuation performance and robust stability in a compensator with low order. Mixed H2/H(infinity), controllers provide a means for maximizing robust stability for a given level of mean-square nominal performance while directly optimizing for controller order constraints. This paper demonstrates the benefit of mixed norm design from the perspective of robustness to parametric uncertainties and controller order for microgravity vibration isolation. A nominal performance metric analogous to the mu measure, for robust stability assessment is also introduced in order to define an acceptable trade space from which different control methodologies can be compared.
Many microgravity science experiments need an active isolation system to provide a sufficiently quiescent acceleration environment. The g-LIMIT vibration isolation system will provide isolation for Microgravity Science Glovebox experiments in the International Space Station. While standard control system technologies have been demonstrated for these applications, modern control methods have the potential for meeting performance requirements while providing robust stability in the presence of parametric uncertainties that are characteristic of microgravity vibration isolation systems. While H2 and H infinity methods are well established, neither provides the levels of attenuation performance and robust stability in a compensator with low order. Mixed H2/mu controllers provide a means for maximizing robust stability for a given level of mean-square nominal performance while directly optimizing for controller order constraints. This paper demonstrates the benefit of mixed norm design from the perspective of robustness to parametric uncertainties and controller order for microgravity vibration isolation. A nominal performance metric analogous to the mu measure for robust stability assessment is also introduced in order to define an acceptable trade space from which different control methodologies can be compared.
An experimental and analytical investigation has been made of the effectiveness of flares and of fins for providing aerodynamic stability of low-fineness-ratio bodies near zero angle of attack. Wind-tunnel tests were performed at Mach numbers from 0.6 to 3.8 to measure the normal force, center of pressure, and drag of bodies consisting of conical noses in combination with cylindrical midsections and finned or flared aftersections. This study included an investigation of the effects of nose bluntness, midbody length, flare angle (from 0° to 20°), and fin leading-edge bluntness. The results showed that for the same plan-form area, flares were more effective than fins, but for the same drag, fins were more effective than flares. Flow separation, which occurred ahead of a flare at certain Mach numbers, resulted in an increase in the stabilizing effectiveness of the flare and a decrease in the drag. However, flow separation led to large undesirable shifts in the center of pressure. Blunt fins were found to be more effective than sharp fins, with the vertical as well as the horizontal fins contributing to the effectiveness. Comparisons of the analytical with the experimental results indicated that the analytical method employed gave good estimates of the stabilizing effectiveness of both fins and flares throughout the supersonic Mach number range.
Schlieren methods of flow visualization and hot-wire anemometry for velocity measurements were used to investigate the wakes generated by hovering model propellers and rotors. The research program was directed toward investigating (1) the stability of the tip vortex, (2) the effects produced by various tip shapes on performance and tip vortex characteristics, and (3) the shock formation and noise characteristics associated with various tip shapes. A free-wake analysis was also conducted for comparison with the vortex stability experimental results. Schlieren photographs showing wake asymmetry, interaction, and instability are presented along with a discussion of the effects produced by the number of blades, collective pitch, and tip speed. Two hot-wire anemometer techniques, used to measure the maximum circumferential velocity in the tip vortex, are discussed.
A large set of variable coefficient linear systems of ordinary differential equations which possess two different time scales, a slow one and a fast one is considered. A small parameter epsilon characterizes the stiffness of these systems. A system of o.d.e.s. in this set is approximated by a general class of multistep discretizations which includes both one-leg and linear multistep methods. Sufficient conditions are determined under which each solution of a multistep method is uniformly bounded, with a bound which is independent of the stiffness of the system of o.d.e.s., when the step size resolves the slow time scale, but not the fast one. This property is called stability with large step sizes. The theory presented lets one compare properties of one-leg methods and linear multistep methods when they approximate variable coefficient systems of stiff o.d.e.s. In particular, it is shown that one-leg methods have better stability properties with large step sizes than their linear multistep counter parts. The theory also allows one to relate the concept of D-stability to the usual notions of stability and stability domains and to the propagation of errors for multistep methods which use large step sizes.
In most studies of stability of plates, the axial stress has been taken as uniform compression throughout flat rectangular plates. Buckling of isotropic plates under a compressive stress that varies linearly from one loaded edge to the other has been studied by Libove et al. Cases of practical interest exist, however, in which the axial stress is not uniform but varies from tension at both loaded edges to compression in the middle. An example is the stability of the crown of the hat stiffened panel, a candidate configuration of the upper and lower skin of the Blended Wing Body (BWB) Aircraft. The BWB Aircraft is an advanced long-range ultra-high-capacity airliner with the principal feature being the pressurized wide double-deck body which is blended into the wing. In the present research, analytical methods are used to investigate the local stability of the crown in order to minimize its weight while optimizing its buckling strength. The crown is modeled as a rectangular laminated composite plate subjected to a second degree parabolic variation of axial stresses in the longitudinal direction. A varying tension-compression- tension axial stresses are induced in the crown of the stiffeners due to bending. The change in axial stresses is equilibrated by nonuniform shear stresses along the plate edges and transverse normal stresses.
Developing general analytical methods for predicting planar and three dimensional attitude motion of gravity gradient spacecraft in elliptical orbit
A method is introduced that weakens the Lyapunov or in track instability of satellite equations of motion. The method utilizes a linearized energy integral of satellite motion as a constraint on solutions obtained by numerical integration. The procedure prevents local numerical error from altering the frequency associated with the fast angular variable and thereby reduces the Lyapunov instability and the global numerical error. Applications of the method to satellite motion show accuracy improvements of two to three orders of magnitude in position and velocity after 50 revolutions. A modification of the method is presented that allows the use of slowly varying integrals of motion.
Numerical methods were developed that will form the computational part of the turbulence closure scheme. A wave model was developed for the two-dimensional shear layer. This configuration is being used as a test case for the closure schemes. Various numerical schemes were examined to give efficient solutions of the Rayleigh equation for this geometry. These include both spectral and finite difference methods. Secondly, numerical methods are under development to solve the non-separable Rayleigh equation. This solution is required for the closure scheme in more complex geometries. A model problem was used to assist in the algorithm development. Two-dimensional spectral methods and a hybrid spectral/finite difference technique were developed. An analytic solution of the Rayleigh equation for a basic elliptic flow was obtained. This will be used to verify the stability codes developed for arbitrary geometries. Other numerical methods for solving the Rayleigh equation based on the boundary element technique were also examined. These solutions are forming the basis of a model for the shock structure in jets of arbitrary geometry.
Among the topics discussed are developments in structural engineering hardware and software, computation for fracture mechanics, trends in numerical analysis and parallel algorithms, mechanics of materials, advances in finite element methods, composite materials and structures, determinations of random motion and dynamic response, optimization theory, automotive tire modeling methods and contact problems, the damping and control of aircraft structures, and advanced structural applications. Specific topics covered include structural design expert systems, the evaluation of finite element system architectures, systolic arrays for finite element analyses, nonlinear finite element computations, hierarchical boundary elements, adaptive substructuring techniques in elastoplastic finite element analyses, automatic tracking of crack propagation, a theory of rate-dependent plasticity, the torsional stability of nonlinear eccentric structures, a computation method for fluid-structure interaction, the seismic analysis of three-dimensional soil-structure interaction, a stress analysis for a composite sandwich panel, toughness criterion identification for unidirectional composite laminates, the modeling of submerged cable dynamics, and damping synthesis for flexible spacecraft structures.
Launch vehicles frequently experience a reduced stability margin through the transonic Mach number range. This reduced stability margin is caused by an undamping of the aerodynamics in one of the lower frequency flexible or rigid body modes. Analysis of the behavior of a flexible vehicle is routinely performed with quasi-steady aerodynamic lineloads derived from steady rigid computational fluid dynamics (CFD). However, a quasi-steady aeroelastic stability analysis can be unconservative at the critical Mach numbers where experiment or unsteady computational aeroelastic (CAE) analysis show a reduced or even negative aerodynamic damping. This paper will present a method of enhancing the quasi-steady aeroelastic stability analysis of a launch vehicle with unsteady aerodynamics. The enhanced formulation uses unsteady CFD to compute the response of selected lower frequency modes. The response is contained in a time history of the vehicle lineloads. A proper orthogonal decomposition of the unsteady aerodynamic lineload response is used to reduce the scale of data volume and system identification is used to derive the aerodynamic stiffness, damping and mass matrices. The results of the enhanced quasi-static aeroelastic stability analysis are compared with the damping and frequency computed from unsteady CAE analysis and from a quasi-steady analysis. The results show that incorporating unsteady aerodynamics in this way brings the enhanced quasi-steady aeroelastic stability analysis into close agreement with the unsteady CAE analysis.
A natural laminar flow (NLF) wind tunnel model has been designed and analyzed for a wind tunnel test in the National Transonic Facility (NTF) at the NASA Langley Research Center. The NLF design method is built into the CDISC design module and uses a Navier-Stokes flow solver, a boundary layer profile solver, and stability analysis and transition prediction software. The NLF design method alters the pressure distribution to support laminar flow on the upper surface of wings with high sweep and flight Reynolds numbers. The method addresses transition due to attachment line contamination/transition, Gortler vortices, and crossflow and Tollmien-Schlichting modal instabilities. The design method is applied to the wing of the Common Research Model (CRM) at transonic flight conditions. Computational analysis predicts significant extents of laminar flow on the wing upper surface, which results in drag savings. A 5.2 percent scale semispan model of the CRM NLF wing will be built and tested in the NTF. This test will aim to validate the NLF design method, as well as characterize the laminar flow testing capabilities in the wind tunnel facility.
Based on the method of the energy principle, the effect of the coriolis force on the stability of rotating hydromagnetic systems was examined and the condition for instability was derived. It is shown that, in rotating systems, the coriolis force inhibits the onset of convective motion and that, once the condition for instability is fulfilled, overstable states are produced with respect to this motion. Such states seem to affect the onset of turbulent convective motion.
We investigate the structure and stability of rotationally symmetric nonhomogeneous time-independent solutions derived from a simple analytic model of self-propagating star formation. For this purpose we employ two methodologies: We use bifurcation theoretical methods to prove the existence of nonhomogeneous axisymmetric stationary solutions of an appropriate nonlinear evolution equation for the stellar density. We show that the nonhomogeneous solution branch bifurcates from the homogeneous one at a critical parameter value of the star formation rate. Further, the analytical theory allows us to show that the new solution set is stable in the weakly nonlinear regime near the bifurcation point. To follow the solution branch further, we use numerical methods. The numerical calculation shows the structure and stability of these solutions. We conclude that no periodic time-dependent solutions of this special model exist, and no further bifurcations can be found. The same results have been found in simulations of stochastic self-propagating star formation based on similar models. Therefore, our findings provide a natural explanation, why long-lived large-scale structure have not been found in those simulations.
Explore the source record for details and available documents.
This report is a compilation of PID (Proportional Integral Derivative) results for both longitudinal and lateral directional analysis that was completed during Fall 1994. It had earlier established that the maneuvers available for PID containing independent control surface inputs from OBES were not well suited for extracting the cross-coupling static (i.e., C(sub N beta)) or dynamic (i.e., C(sub Npf)) derivatives. This was due to the fact that these maneuvers were designed with the goal of minimizing any lateral directional motion during longitudinal maneuvers and vice-versa. This allows for greater simplification in the aerodynamic model as far as coupling between longitudinal and lateral directions is concerned. As a result, efforts were made to reanalyze this data and extract static and dynamic derivatives for the F/A-18 HARV (High Angle of Attack Research Vehicle) without the inclusion of the cross-coupling terms such that more accurate estimates of classical model terms could be acquired. Four longitudinal flights containing static PID maneuvers were examined. The classical state equations already available in pEst for alphadot, qdot and thetadot were used. Three lateral directional flights of PID static maneuvers were also examined. The classical state equations already available in pEst for betadot, p dot, rdot and phi dot were used. Enclosed with this document are the full set of longitudinal and lateral directional parameter estimate plots showing coefficient estimates along with Cramer-Rao bounds. In addition, a representative time history match for each type of meneuver tested at each angle of attack is also enclosed.