SEARCH · Engineering Papers
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.
The analysis and the FORTRAN program are presented for the determination of attitude of a spin-stabilized spacecraft. The use of telemetry data that provide information about two reference vectors and their relation to the spin is outlined. A technique for the determination of the spin-axis orientation that employs only simple calculations is described.
Permanent magnet excited linear alternators rated tens of kW and coupled to free-piston Stirling engines are presently viewed as promising candidates for long term generation of electric power in both space and terrestrial applications. Series capacitive cancellation of the internal inductive reactance of such alternators was considered a viable way to both increase power extraction and to suppress unstable modes of the thermodynamic oscillation. Idealized toroidal and cylindrical alternator geometries are used for a comparative study of the issues of specific mass and capacitive tuning, subject to stability criteria. The analysis shows that the stator mass of an alternator designed to be capacitively tuned is always greater than the minimum achievable stator mass of an alternator designed with no capacitors, assuming equal utilization of materials ratings and the same frequency and power to a resistive load. This conclusion is not substantially altered when the usually lesser masses of the magnets and of any capacitors are added. Within the reported stability requirements and under circumstances of normal materials ratings, this study finds no clear advantage to capacitive tuning. Comparative plots of the various constituent masses are presented versus the internal power factor taken as a design degree of freedom. The explicit formulas developed for stator core, coil, capacitor, and magnet masses and for the degree of magnet utilization provide useful estimates of scaling effects.
Coupled flap-lag-torsion equations of motion of an isolated horizontal axis wind turbine (HAWT) blade have been formulated. The analysis neglects blade-tower coupling. The final nonlinear equations have periodic coefficients. A new and convenient method of generating an appropriate time-dependent equilibrium position, required for the stability analysis, has been implemented and found to be computationally efficient. Steady-state response and stability boundaries for an existing (typical) HAWT blade are presented. Such stability boundaries have never been published in the literature. The results show that the isolated blade under study is basically stable. The tower shadow (wake) has a considerable effect on the out-of-plane response but leaves blade stability unchanged. Nonlinear terms can significantly affect linearized stability boundaries; however, they have a negligible effect on response, thus implying that a time-dependent equilibrium position (or steady-state response), based completely on the linear system, is appropriate for the type of HAWT blades under study.
Experimental aerodynamic investigations were made on a .006 scale model 040-A delta wing space shuttle orbiter configuration. These tests were conducted to determine six-degree-of-freedom force and moment data for preliminary stability and control analysis. Data were obtained over a Mach number range from 0.6 to 4.96 at angles of attack from -10 deg to 50 deg at zero degrees sideslip and at angles of sideslip from -10 deg to 10 deg at constants angles of attack of 0 deg, 15 deg, 30 deg, and 45 deg. Various aileron, elevator, (elevon) rudder and rudder flare deflection angles were tested to establish the control effectiveness and vehicle stability. Model component buildup data were also obtained to provide a data base for future configuration modifications. Plotted data results are presented in both the body and stability axis system.
FeO x , SrFeO 3-y (SF), and Ba 0.5 Sr 0.5 FeO 3–y (BSF) were studied for application as fuel oxidation catalysts in solid oxide fuel cells (SOFC) anodes. Electrodes were prepared by impregnation into porous yttria-stabilized zirconia (YSZ), with La 0.3 Sr 0.7 TiO 3–y (LST) added for electronic conductivity. The electrode impedances decreased dramatically upon addition of SF and BSF and much less when only SrO or FeO x were added. Temperature Programmed Desorption (TPD) of O 2 from oxidized BSF showed O 2 desorbing between 200 °C and 700 °C, while no O 2 desorbed from Fe 2 O 3 below 900 °C. The results, together with thermodynamic analysis, suggest that stabilization of Fe +4 in the perovskite lattice plays an important role in enhancing the catalytic activity of SF and BSF by providing access to a Fe 3+ /Fe 4+ redox couple that can accept oxygen anions from the yttria-stabilized zirconia (YSZ) electrolyte. Oxygen is weakly bound to these sites facilitating its reaction with adsorbed H 2 .
Analytical method for predicting stability and control characteristics of large elastic airplanes at subsonic and supersonic speeds
A method is proposed to determine directly the amount of artificial viscosity needed for stability using an eigenvalue analysis for a finite difference representation of the Navier-Stokes equations. The stability and growth of small perturbations about a steady flow over the airfoils are analyzed for various amounts of artificial viscosity. The eigenvalues were determined for a small perturbation about a steady inviscid flow over a NACA 0012 airfoil at a Mach number of 0.8 and angle of attack of 0 degrees. The movement of the eigenvalue constellation with respect to the amount of artificial viscosity is studied. The stability boundries as a function of the amount of artificial viscosity from both the eigenvalue analysis and the time marching scheme are also presented. This procedure not only allows for determining the effect of varying amounts of artificial viscosity, but also for the effects of different forms of terms for artificial viscosity.
A method is proposed to determine directly the amount of artificial viscosity needed for stability using an eigenvalue analysis for a finite difference representation of the Navier-Stokes equations. The stability and growth of small perturbations about a steady flow over airfoils are analyzed for various amounts of artificial viscosity. The eigenvalues were determined for a small time-dependent perturbation about a steady inviscid flow over an NACA 0012 airfoil at a Mach number of 0.8 and angle of attack of 0 deg. The method has been applied to inviscid flows here, but as discussed is also applicable to viscous flows. The movement of the eigenvalue constellation with respect to the amount of artificial viscosity is studied. The stability boundaries as a function of the amount of artificial viscosity from both the eigenvalue analysis and the time-marching scheme are also presented. The eigenvalue procedure not only allows for determining the effect of varying amounts of artificial viscosity, but also for the effects of different forms of artificial viscosity.
Various tools for the aeroelastic stability and response analysis of rotor blades in hover and forward flight were developed and incorporated in a comprehensive package capable of performing aeroelastic tailoring of rotor blades in forward flight. The results indicate that substantial vibration reductions, of order 15-40%, in the vibratory hub shears can be achieved by relatively small modifications of the initial design. Furthermore the optimized blade can be up to 20% lighter than the original design. Accomplishments are reported for the following tasks: (1) finite element modeling of rotary-wing aeroelastic problems in hover and forward flight; (2) development of numerical methods for calculating the aeroelastic response and stability of rotor blades in forward fight; (3) formulation of the helicopter air resonance problem in hover with active controls; and (4) optimum design of rotor blades for vibration reduction in forward flight.
Spatial operators have been used to analyze the dynamics of robotic multibody systems and to develop novel computational dynamics algorithms. Mass matrix factorization, inversion, diagonalization, and linearization are among several new insights obtained using such operators. While initially developed for serial rigid body manipulators, the spatial operators and the related mathematical analysis have been shown to extend very broadly including to tree and closed topology systems, to systems with flexible joints, links, etc. This work uses concepts from graph theory to explore the mathematical foundations of spatial operators. The goal is to study and characterize the properties of the spatial operators at an abstract level so that they can be applied to a broader range of dynamics problems. The rich mathematical properties of the kinematics and dynamics of robotic multibody systems has been an area of strong research interest for several decades. These properties are important to understand the inherent physical behavior of systems, for stability and control analysis, for the development of computational algorithms, and for model development of faithful models. Recurring patterns in spatial operators leads one to ask the more abstract question about the properties and characteristics of spatial operators that make them so broadly applicable. The idea is to step back from the specific application systems, and understand more deeply the generic requirements and properties of spatial operators, so that the insights and techniques are readily available across different kinematics and dynamics problems. In this work, techniques from graph theory were used to explore the abstract basis for the spatial operators. The close relationship between the mathematical properties of adjacency matrices for graphs and those of spatial operators and their kernels were established. The connections hold across very basic requirements on the system topology, the nature of the component bodies, the indexing schemes, etc. The relationship of the underlying structure is intimately connected with efficient, recursive computational algorithms. The results provide the foundational groundwork for a much broader look at the key problems in kinematics and dynamics. The properties of general graphs and trees of nodes and edge were examined, as well as the properties of adjacency matrices that are used to describe graph connectivity. The nilpotency property of such matrices for directed trees was reviewed, and the adjacency matrices were generalized to the notion of block weighted adjacency matrices that support block matrix elements. This leads us to the development of the notion of Spatial Kernel Operator SKO kernels. These kernels provide the basis for the development of SKO resolvent operators.
Consideration of the role of stability of the mean jet flow of an edgetone system in the edgetone operation. A brief description of the main features of the edgetone operation and a brief review of the pertinent background are presented. This is followed by a discussion of the inadequacy of parallel flow stability considerations in explaining edgetone features. An analysis is made of the stability characteristics of a slightly nonparallel jet flow, and the role of jet stability based on these considerations in the edgetone operation is discussed. It is shown that the characteristics obtained by considering the nonparallel nature of the jet flowfield lead to a satisfactory description of the phenomenon, including a frequency selection criterion based on the concept of the most amplified disturbance frequency at each edge position.
An analysis and calculations on the stability of the free laminar boundary layer between parallel streams were made for an incompressible fluid using the Tollmien-Schlichting theory of small disturbances. Because the boundary conditions are at infinity, two solutions of the Orr-Sommerfeld stability equations need not be considered, and the remaining two solutions are exponential in character at the infinite boundaries. The calculations show that the flow is unstable except for very low Reynolds numbers.
This computational study investigates the effects of discrete roughness elements on a blunt cone at zero degrees angle of attack in a Mach 6 flow. Motivation was provided by experiments conducted in the Air Force Research Laboratory Mach 6 High Reynolds Number facility on a 7-degree half-angle cone with a roughness array located at 45 degrees from the apex on two nose tips of different blutness but equivalent roughness Reynolds number. Transition was only affected on the blunter cone, indicating that the transition onset is associated with the combined effects of bluntness and roughness. The present study investigates the 15.24 mm nose radius, 420 azimuthal wavenumber case via Navier-Stokes computations of the laminar base flow and instability analysis. Plane-marching parabolized stability equations (PSE) and inflow-resolvent analysis based on the three-dimensional, harmonic linearized Navier-Stokes equations (HLNSE) are used to calculate the amplification of disturbances along the roughness wake as well as over the roughness nearfield. Results show that the roughness shape can have a great impact on the characteristics of the most amplified wake instabilities. For the experimental configuration with cubic roughness elements of 15 μ m height, the flow is marginally unstable. For prismatic elements of 20 μ m height, the PSE predicts a logarithmic disturbance amplification ratio of N = 5.6 along its wake, but this ratio increases to N = 9.6 when the amplification over the roughness and separation regions is included in the inflow-resolvent analysis.
Areas relating to the development and improvement of the single-fuselage, yawed-wing transonic transport concept were investigated. These included: (1) developing an alternate configuration with a simplified engine installation;(2) determining a structural design speed placard that would allow the engine-airframe match for optimum airplane performance; and (3) conducting an aeroelastic stability and control analysis of the yawed-wing configuration with a flexible wing. A two-engine, single-fuselage, yawed-wing configuration was developed that achieved the Mach 1.2 design mission at 5560 km (3000 nmi) and payload of 18,140 kg (40,000 lb) with a gross weight of 217,700 kg (480,000 lb). This airplane was slightly heavier than the aft-integrated four-engine configuration that had been developed in a previous study. A modified structural design speed placard, which was determined, resulted in a 6% to 8% reduction in the gross weight of the yawed-wing configurations. The dynamic stability characteristics of the single-fuselage yawed-wing configuration were found to be very dependent on the magnitude of the pitch/roll coupling, the static longitudinal stability, and the dihedral effect.
The classical theory of flight dynamics for airplane longitudinal stability and control analysis was extended to the case of a hypervelocity reentry vehicle. This includes the elements inherent in supersonic and hypersonic flight such as the influence of the Mach number on aerodynamic characteristics, and the effect of the reaction control system and aerodynamic controls on the trim condition through a wide range of speed. Phugoid motion and angle of attack oscillation for typical cases of cruising flight, ballistic entry, and glide entry are investigated. In each case, closed form solutions for the variations in altitude, flight path angle, speed and angle of attack are obtained. The solutions explicitly display the influence of different regions design parameters and trajectory variables on the stability of the motion.
The preliminary conceptual design of the Korean fusion demonstration reactor (K-DEMO) with a major radius of 6.8 m and the fusion power of 2200 MW has been studied since 2012. The overall configuration of the K-DEMO divetor system based on the ITER-like water-cooled tungsten technology is a double-null type symmetric divertor subdivided into 32 toroidal modules for the vertical maintenance. A detached divertor scenario with impurity seeding was considered as the primary approach for the power exhaust to reduce the peak heat flux lower than the engineering limit of 10 MW/m 2 . The power exhaust performance at the scrap off layer was estimated by using UEDGE-2D code, a two-dimensional fluid transport code for collisional edge plasma and neutral species like N, Ne, and Ar. Particle and heat flux on inboard and outboard divertor targets were calculated for the detached cased depending on parameters such as the impurity seeding rate, pumping rate, and the pedestal density. On the other hand, a magnetic solution like X-divertor, snowflake divertor, and super X-divertor to expand the plasma wet area was considered for K-DEMO since the detached divertor increasing a radiation fraction by impurity seeding might be able to be unstable. However, the extremely high current of poloidal coils was required more than the engineering limit, 20 MA, to form magnetic field lines for the alternative divertors. Based on the physical calculation of the edge plasma, engineering analyses were carried out to find out the thermal and structural reliability. The thermo-hydraulic analysis confirmed thermal stability, whether all comprising materials are operating within their allowable temperature windows when the case of the peak heat flux is set to 10 MW/m 2 on the outboard divertor target. The response surface optimization method derived two optimal design candidates employing two kinds of heat sink materials, respectively: the reduced activation ferritic martensitic (RAFM) steel and CuCrZr alloy. Additionally, the drawbacks and merits of the two materials were definite. The optimal design with applying RAFM steel was vulnerable to withstand thermal and mechanical loads since low thermal conductivity caused too thin thickness of the heat sink. On the other hand, the CuCrZr alloy has critical drawbacks in terms of activation and radioactive waste despite its high thermal conductivity. Meanwhile, preliminary electromagnetic (EM) analysis was carried out to estimate the EM loads caused by the abnormal behaviors of plasma since EM loads are one of the most critical external loads for designing a DEMO divertor.
This paper presents an aeroservoelastic flight dynamic modeling approach for NASA Lift+Cruise eVTOL aircraft in forward flight. A conceptual aerodynamic analysis for forward flight is performed in CHARM, a rotorcraft aeromechanics unsteady vortex-lattice flow solver. A trim analysis is performed to obtain the steady state aerodynamic information at a trim point to be used for the aeroelastic analysis. A NASTRAN finite element model for the Lift+Cruise vehicle is available. The structural dynamic mode shapes are extracted from the NASTRAN for the aeroelastic analysis. An aeroservoelastic model is developed using the Theodorsen’s theory. A stability and control derivative analysis is performed to formulate a flight dynamic model of the Lift+Cruise vehicle. An aeroservoelastic flight dynamic model is formulated by coupling these two models together. Flight dynamic and flutter analyses are performed to determine vehicle stability and flutter-free operation.
Dynamic modeling of a low cost single axis LST system and the pointing stability of that system are investigated. The effects of nonlinear friction of the bearings of the reaction wheels, quantization, and sensor noise on the pointing error are covered. In addition, data are given on self sustained oscillations of the system induced by quantization and methods of evaluating attitude error of the digital LST.