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.
Report shows how spatial-operator algebra developed by authors used to derive linearized mathematical models of dynamics of robotic manipulators. Aspects of spatial-operator algebra also described in "Robot Control Based On Spatial-Operator Algebra" (NPO-17918), "Unified Formulation Of Dynamics Of Robot Arms" (NPO-18040), and "Spatial-Operator Algebra For Robotic Manipulators" (NPO-17770).
A standard problem in gamma-ray astronomy data analysis is the decomposition of a set of observed counts, described by Poisson statistics, according to a given multicomponent linear model, with underlying physical count rates or fluxes which are to be estimated from the data. Despite its conceptual simplicity, the linear least-squares (LLSQ) method for solving this problem has generally been limited to situations in which the number n(sub i) of counts in each bin i is not too small, conventionally more than 5-30. It seems to be widely believed that the failure of the LLSQ method for small counts is due to the failure of the Poisson distribution to be even approximately normal for small numbers. The cause is more accurately the strong anticorrelation between the data and the wieghts w(sub i) in the weighted LLSQ method when square root of n(sub i) instead of square root of bar-n(sub i) is used to approximate the uncertainties, sigma(sub i), in the data, where bar-n(sub i) = E(n(sub i)), the expected value of N(sub i). We show in an appendix that, avoiding this approximation, the correct equations for the Poisson LLSQ (PLLSQ) problems are actually identical to those for the maximum likelihood estimate using the exact Poisson distribution. We apply the method to solve a problem in high-resolution gamma-ray spectroscopy for the JPL High-Resolution Gamma-Ray Spectrometer flown on HEAO 3. Systematic error in subtracting the strong, highly variable background encountered in the low-energy gamma-ray region can be significantly reduced by closely pairing source and background data in short segments. Significant results can be built up by weighted averaging of the net fluxes obtained from the subtraction of many individual source/background pairs. Extension of the approach to complex situations, with multiple cosmic sources and realistic background parameterizations, requires a means of efficiently fitting to data from single scans in the narrow (approximately = 1.2 keV, HEAO 3) energy channels of a Ge spectrometer, where the expected number of counts obtained per scan may be very low. Such an analysis system is discussed and compared to the method previously used.
Since numerical weather prediction models are unable to accurately forecast the severity and the location of the storm cells several hours into the future when compared with observation data, there has been a growing interest in probabilistic description of convective weather. The classical approach for generating uncertainty bounds consists of integrating the state equations and covariance propagation equations forward in time. This step is readily recognized as the process update step of the Kalman Filter algorithm. The second well known method, known as the Monte Carlo method, consists of generating output samples by driving the forecast algorithm with input samples selected from distributions. The statistical properties of the distributions of the output samples are then used for defining the uncertainty bounds of the output variables. This method is computationally expensive for a complex model compared to the covariance propagation method. The main advantage of the Monte Carlo method is that a complex non-linear model can be easily handled. Recently, a few different methods for probabilistic forecasting have appeared in the literature. A method for computing probability of convection in a region using forecast data is described in Ref. 5. Probability at a grid location is computed as the fraction of grid points, within a box of specified dimensions around the grid location, with forecast convection precipitation exceeding a specified threshold. The main limitation of this method is that the results are dependent on the chosen dimensions of the box. The examples presented Ref. 5 show that this process is equivalent to low-pass filtering of the forecast data with a finite support spatial filter. References 6 and 7 describe the technique for computing percentage coverage within a 92 x 92 square-kilometer box and assigning the value to the center 4 x 4 square-kilometer box. This technique is same as that described in Ref. 5. Characterizing the forecast, following the process described in Refs. 5 through 7, in terms of percentage coverage or confidence level is notionally sound compared to characterizing in terms of probabilities because the probability of the forecast being correct can only be determined using actual observations. References 5 through 7 only use the forecast data and not the observations. The method for computing the probability of detection, false alarm ratio and several forecast quality metrics (Skill Scores) using both the forecast and observation data are given in Ref. 2. This paper extends the statistical verification method in Ref. 2 to determine co-occurrence probabilities. The method consists of computing the probability that a severe weather cell (grid location) is detected in the observation data in the neighborhood of the severe weather cell in the forecast data. Probabilities of occurrence at the grid location and in its neighborhood with higher severity, and with lower severity in the observation data compared to that in the forecast data are examined. The method proposed in Refs. 5 through 7 is used for computing the probability that a certain number of cells in the neighborhood of severe weather cells in the forecast data are seen as severe weather cells in the observation data. Finally, the probability of existence of gaps in the observation data in the neighborhood of severe weather cells in forecast data is computed. Gaps are defined as openings between severe weather cells through which an aircraft can safely fly to its intended destination. The rest of the paper is organized as follows. Section II summarizes the statistical verification method described in Ref. 2. The extension of this method for computing the co-occurrence probabilities in discussed in Section HI. Numerical examples using NCWF forecast data and NCWD observation data are presented in Section III to elucidate the characteristics of the co-occurrence probabilities. This section also discusses the procedure for computing throbabilities that the severity of convection in the observation data will be higher or lower in the neighborhood of grid locations compared to that indicated at the grid locations in the forecast data. The probability of coverage of neighborhood grid cells is also described via examples in this section. Section IV discusses the gap detection algorithm and presents a numerical example to illustrate the method. The locations of the detected gaps in the observation data are used along with the locations of convective weather cells in the forecast data to determine the probability of existence of gaps in the neighborhood of these cells. Finally, the paper is concluded in Section V.
Linear multivariable regression models for predicting day and night Eddy Dissipation Rate (EDR) from available meteorological data sources are defined and validated. Model definition is based on a combination of 1997-2000 Dallas/Fort Worth (DFW) data sources, EDR from Aircraft Vortex Spacing System (AVOSS) deployment data, and regression variables primarily from corresponding Automated Surface Observation System (ASOS) data. Model validation is accomplished through EDR predictions on a similar combination of 1994-1995 Memphis (MEM) AVOSS and ASOS data. Model forms include an intercept plus a single term of fixed optimal power for each of these regression variables; 30-minute forward averaged mean and variance of near-surface wind speed and temperature, variance of wind direction, and a discrete cloud cover metric. Distinct day and night models, regressing on EDR and the natural log of EDR respectively, yield best performance and avoid model discontinuity over day/night data boundaries.
A review is presented of various rotordynamic problems which have been encountered and eliminated in developing the current flight engines, and continuing subsynchronous problems which are being encountered in developing a 109% power level engine. The basic model for the High Pressure Oxygen Turbopump (HPOTP) of the SSME including the structural dynamic model for the rotor and housing and component models for the liquid and gas seals, turbine-clearance excitation forces, and impeller-diffuser forces are discussed. Results from a linear model are used to examine the synchronous response and stability characteristics of the HPOTP, examining bearing load and stability problems associated with the second critical speed. Various seal modifications are examined and shown to have favorable consequences with respect to bearing reactions and stability.
It is argued that if the nonlinearities in a system are mild, and the controller is sufficiently stabilizing, the inaccuracies of a linear model, which is often taken to be sufficient for a controller design, can be safely neglected. For systems with severe nonlinearity a linearizing technique is described, based on the change of state coordinates and nonlinear feedback; in the total context of a stable feedback design the linearization technique is considered robust. Furthermore, attention is paid to partial linearization using the same transformations. It is found that there always exist maximally linearizing transformations which are not necessarily unique.
Techniques are described for constructing analytic models of the chemical evolution of galaxies subject to infall of metal-poor material onto a maturing disk. A class of linear models is discussed which takes the star-formation rate within a defined region to be proportional to the mass of interstellar gas within that region, and the instantaneous recycling approximation is adopted. The solutions are obtained by approximately matching the infall rate to parametrized familiies of functions for which the equations are exactly soluble. The masses, the primary and secondary metallicities, and the gas concentrations of radioactive chronometers can all then be analytically expressed. Surveys of galactic abundances in location and in time can be compared to the parameter spaces of the analytic representations.
This paper analyzes the dynamics of a long tether connecting two spacecraft in earth orbit, one of the spacecraft having dominant mass. In particular, it considers the material damping of the tether. The nominal position of the tether is stabilized by the gravity gradient such that it is aligned with the local vertical. The tether is modeled as a viscoelastic flexible continuum. Modal frequencies are derived in an analytical approximation form. Damping ratios are estimated according to the linear model calibrated by ground measurements. The results show that, with properly chosen tether material and braiding structure, longitudinal vibrations of the tethered system are well damped.
This paper discusses an approach to the modeling and performance for the preliminary design phase of a large (6.2 MW) horizontal axis wind turbine generator (WTG). Two control philosophies are presented, both of which are based on linearized models of the WT mechanical and electrical systems. The control designs are compared by showing the performance through detailed non-linear time simulation. The disturbances considered are wind gusts, and electrical faults near the WT terminals.
Forebody Vortex Control (FVC) is an emerging technology that has received widespread and concentrated attention by many researchers for application on fighter aircraft to enhance aerodynamic controllability at high angles of attack. This research explores potential application of FVC to a NASP-type configuration. The configuration investigated is characterized by a slender, circular cross-section forebody and a 78 deg swept delta wing. A man-in-the-loop, six-degress-of-freedom, high-fidelity simulation was developed that demonstrates the implementation and advantages of pneumatic forebody vortex control. Static wind tunnel tests were used as the basis for the aerodynamic characteristics modeled in the simulation. Dynamic free-to-roll wind tunnel tests were analyzed and the wing rock motion investigated. A non-linear model of the dynamic effects of the bare airframe and the forebody vortex control system were developed that closely represented the observed behavior. Multiple state-of-the-art digital flight control systems were developed that included different utilizations of pneumatic vortex control. These were evaluated through manned simulation. Design parameters for a pneumatic forebody vortex control system were based on data collected regarding the use of blowing and the mass flow required during realistic flight maneuvers.
The increasing electrification of aircraft propulsion systems is leading to new control architectures being developed to address integration between electric machines and gas-based turbine engines. For hybrid-electric propulsion systems, current conceptual architectures often couple electric machines with the shafts of gas turbine engines and introduce energy storage. Leveraging the electrical power system of hybridized engines, Turbine Electrified Energy Management (TEEM) is a recent control approach that improves transient operability in an effort to enable more efficient and lighter weight turbomachinery. This study seeks to expand TEEM’s application beyond traditional proportional integral (PI) control by presenting linear model predictive control (MPC) schemes to execute the TEEM concept. Through constraint selection and cost function design, transient operability goals for TEEM are considered with no external logic or saturation. Unique to the designs are the use of a washout filter, which simplifies transient detection and motor activation logic. The proposed architectures are implemented with both centralized MPC and distributed MPC approaches, and comparisons are drawn to a benchmark PI controller simulated on a nonlinear turbofan engine model at one ground condition and one cruise condition. Performance is evaluated using compressor maps, stall margin performance, and two novel metrics: transient stack usage and transient excursion integral. Results reveal the linear MPC scheme performs comparably to the baseline controller and can be implemented in at least two distinct configurations with potential for further modifications, thus establishing the groundwork for future investigations.
The increasing electrification of aircraft propulsion systems is leading to new control architectures being developed to address integration between electric machines and gas-based turbine engines. For hybrid-electric propulsion systems, current conceptual architectures often couple electric machines with the shafts of gas turbine engines and introduce energy storage. Leveraging the electrical power system of hybridized engines, Turbine Electrified Energy Management (TEEM) is a recent control approach that improves transient operability in an effort to enable more efficient and lighter weight turbomachinery. This study seeks to expand TEEM’s application beyond traditional proportional-integral (PI) control by presenting linear model predictive control (MPC) schemes to execute the TEEM concept. Through constraint selection and cost function design, transient operability goals for TEEM are considered with no external logic or saturation. Unique to the designs are the use of a washout filter, which simplifies transient detection and motor activation logic. The proposed architectures are implemented with both centralized MPC and distributed MPC approaches, and comparisons are drawn to a benchmark PI controller simulated on a nonlinear turbofan engine model at one ground condition and one cruise condition. Performance is evaluated using compressor maps, stall margin performance, and two novel metrics: transient stack usage and transient excursion integral. Results reveal the linear MPC scheme performs comparably to the baseline controller and can be implemented in at least two distinct configurations with potential for further modifications, thus establishing the groundwork for future investigations.
The increasing electrification of aircraft propulsion systems is leading to new control architectures being developed to address integration between electric machines and gas-based turbine engines. For hybrid-electric propulsion systems, current conceptual architectures often couple electric machines with the shafts of gas turbine engines and introduce energy storage. Leveraging the electrical power system of hybridized engines, Turbine Electrified Energy Management (TEEM) is a recent control approach that improves transient operability in an effort to enable more efficient and lighter weight turbomachinery. This study seeks to expand TEEM’s application beyond traditional proportional-integral (PI) control by presenting linear model predictive control (MPC) schemes to execute the TEEM concept. Through constraint selection and cost function design, transient operability goals for TEEM are considered with no external logic or saturation. Unique to the designs are the use of a washout filter, which simplifies transient detection and motor activation logic. The proposed architectures are implemented with both centralized MPC and distributed MPC approaches, and comparisons are drawn to a benchmark PI controller simulated on a nonlinear turbofan engine model at one ground condition and one cruise condition. Performance is evaluated using compressor maps, stall margin performance, and two novel metrics: transient stack usage and transient excursion integral. Results reveal the linear MPC scheme performs comparably to the baseline controller and can be implemented in at least two distinct configurations with potential for further modifications, thus establishing the groundwork for future investigations.
Minimum-time accelerations of the F100 turbofan engine are presented. A piecewise-linear engine model, having three state variables and four control variables, is used to obtain the minimum-time solutions. The linear model which applies at a given time in the trajectory is determined by calculating a normalized distance from the current state to the equilibrium state associated with each linear model. The linear model associated with the closest equilibrium point is then used. The control histories for the minimum-time solutions are used as input to a nonlinear simulation of the F100 engine to verify the accuracy of the piecewise linear solutions.
Minimum-time accelerations of the F100 turbofan engine are presented. A piecewise-linear engine model, having three state variables and four control variables, is used to obtain the minimum-time solutions. The linear model which applies at a given time in the trajectory is determined by calculating a normalized 'distance' from the current state to the equilibrium state associated with each linear model. The linear model associated with the closest equilibrium point is then used. The control histories for the minimum-time solutions are used as input to a nonlinear simulation of the F100 engine to verify the accuracy of the piecewise-linear solutions.
This paper derives a symbolic multi-body rigid nonlinear model for a space telescope using Stoneking’s implementation of Kane’s method. This symbolic nonlinear model is linearized using Matlab symbolic functions diff and inv because the analytic linearization is intractable for manual derivation. The linearized system model is then used to design the controllers using both linear quadratic regulator (LQR) and robust pole assignment methods. The closed-loop systems for the two designs are simulated using both the rigid model as well as a second model containing flexible modes. The performances of the two designs are compared based on the simulation testing results. Our conclusion is that the robust pole assignment design offers better performance than that of the LQR system in terms of actuator usage and pointing accuracy. However, the LQR approach remains an effective first design step that can inform the selection of real eigenvalues for robust pole assignment. The proposed method may be used for the modeling and controller designs for various multi-body systems.
Statistical model, which consists of on/off binary system, assumes Poisson detection process and log-normal atmospheric scintillation. Based upon detection process and atmospheric ionization, piecewise linear model for adaptive threshold system is developed.
A model of the nonlinear dynamics of a deformable maneuvering multibody system is described, whereby elastic deformation are modeled by restoring forces and dissipative forces at point mass appendages. This model is brought into bilinear form. Estimation of deformations occasioned by rapid slewing maneuvers is carried out by a filter based on a globally equivalent linear model of the bilinear dynamics, and is shown to be an improvement over the extended Kalman filter. To further alleviate the computational burden, the estimated deformation state is propagated between observations by a low dimensional operator spline interpolator of bilinear system Volterra series, which is easily implemented.