Search NASA⌕ Search

SEARCH · Search NASA

Results for “linear estimation”

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 757 records · Page 42

Solution Irregularity Remediation for Spatial Discretization Error Estimation for S N Transport Solutions

The discrete ordinates linear Boltzmann transport equation is typically solved in its spatially discretized form, incurring spatial discretization error. Quantification of this error for purposes such as adaptive mesh refinement or error analysis requires an a posteriori estimator, which utilizes the numerical solution to the spatially discretized equation to compute an estimate. Because the quality of the numerical solution informs the error estimate, irregularities, present in the true solution for any realistic problem configuration, tend to cause the largest deviation in the error estimate vis-a-vis the true error. In this paper, an analytical partial singular characteristic tracking (pSCT) procedure for reducing the estimator’s error is implemented within our novel residual source estimator for a zeroth-order discontinuous Galerkin scheme, at the additional cost of a single inner iteration. Here, a metric-based evaluation of the pSCT scheme versus the standard residual source estimator is performed over the parameter range of a Method of Manufactured Solutions test suite. The pSCT scheme generates near-ideal accuracy in the estimate in problems where the dominant source of the estimator’s error is the solution irregularity, namely, problems where the true solution is discontinuous and problems where the true solution’s first derivative is discontinuous and the scattering ratio is low. In problems where the scattering ratio is high and the true solution is discontinuous in the first derivative, the error in the scattering source, which is not converged by the pSCT scheme, is greater than the error incurred due to the irregularity. Ultimately, a pSCT scheme is judged to be useful for error estimation in problems where the computational cost of the scheme is justified. In the presence of many irregularities, such a scheme may be intractable for general use, but in benchmarks, as an analytical tool, or in problems that have nondissipative discontinuities, the scheme may prove invaluable.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Comparison of stability and control parameters for a light, single-engine, high-winged aircraft using different flight test and parameter estimation techniques

The longitudinal and lateral stability and control parameters for a high wing, general aviation, airplane are examined. Estimations using flight data obtained at various flight conditions within the normal range of the aircraft are presented. The estimations techniques, an output error technique (maximum likelihood) and an equation error technique (linear regression), are presented. The longitudinal static parameters are estimated from climbing, descending, and quasi steady state flight data. The lateral excitations involve a combination of rudder and ailerons. The sensitivity of the aircraft modes of motion to variations in the parameter estimates are discussed.

Suit, W. T.↗

Data-Driven Distributed Algorithms for Estimating Eigenvalues and Eigenvectors of Interconnected Dynamical Systems

Here, the paper presents data-driven algorithms to estimate in a distributed manner the eigenvalues, right and left eigenvectors of an unknown linear (or linearized) interconnected dynamic system. In particular, the proposed algorithms do not require the identification of the system model in advance before performing the estimation. As a first step, we consider interconnected dynamical system with distinct eigenvalues. The proposed strategy first estimates the eigenvalues using the well-known Prony method. The right and left eigenvectors are then estimated by solving distributively a set of linear equations. One important feature of the proposed algorithms is that the topology of communication network used to perform the distributed estimation can be chosen arbitrarily, given that it is connected, and is also independent of the structure or sparsity of the system (state) matrix. The proposed distributed algorithms are demonstrated via a numerical example.

97 MATHEMATICS AND COMPUTING↗

Integrated system identification and modal state estimation for control of flexible space structures

A novel approach of integrated system identification and modal state estimation is proposed for control of linear dynamical systems including flexible space structures. There are four steps involved in this approach. First, the relation between a stochastic state space model of a dynamical system and the coefficients of its autoregressive model with exogenous input is derived. Second, an adaptive least-squares transversal predictor is used to estimate the coefficients of the model. Third, a state space model and a steady state Kalman filter gain of the dynamical system are then identified from the coefficients of the model by using the eigensystem realization algorithm. Fourth, a modal state estimator is constructed using the modal parameters of the identified model. On-line implementation of this algorithm can continually improve the modal parameters and the filter gain. It can also gradually update the system model when the system characteristics are slowly changing. A numerical example is used to illustrate the feasibility of the new approach.

Chen, Chung-Wen↗

Potential application of multipolarization SAR for pine-plantation biomass estimation

This paper presents the technique and the potential utility of multipolarization Synthetic Aperture Radar (SAR) data for pine-plantation biomass estimation. Three channels of SAR data, one from the Shuttle Imaging Radar SIR-A and the other two from the aircraft SAR, were acquired over the Baldwin County, Alabama, study area. The SIR-A data were acquired with HH polarization and the aircraft SAR data with VV and VH polarizations. Linear regression techniques are used to estimate the pine-plantation biomass, tree height, and age using 21 test plots. The results indicate that the multipolarization data are highly related to the plantation biomass. The results suggest a potential application of multipolarization SAR for pine-plantation biomass estimation.

Wu, Shih-Tseng↗

Anomalous cross-field transport in a Hall thruster inferred from direct measurement of instability growth rates

The contribution of the electron drift instability to anomalous electron transport is experimentally assessed in a Hall effect discharge. The transport is represented by an anomalous collision frequency, which is related through quasilinear theory to the energy and growth rate of the instability. The wave energy is measured directly with ion saturation probes, while estimates of the growth rate are employed based on both linearized theory and direct measurement. The latter measurement is performed with a bispectral analysis method. The wave-driven collision frequency is compared to measurements of the actual collision frequency inferred from a method based on laser- induced fluorescence. It is found that estimates for transport using linearized theory for the growth differ by over an order of magnitude from the actual anomalous collision frequency in the plasma. The wave-driven anomalous collision frequency with measured growth, however, is shown to agree with the electron collision frequency in magnitude and capture aspects of the trends in spatial variation. This result demonstrates experimentally that wave-driven effects ultimately can explain the observed cross-field transport in these devices. As a result, the implications of this finding are discussed in the context of the key lengthscales that drive the transport as well as the implications identifying reduced fidelity models that could be used to predict anomalous collision frequency.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Validation of a computer code for analysis of subsonic aerodynamic performance of wings with flaps in combination with a canard or horizontal tail and an application to optimization

Extensive correlations of computer code results with experimental data are employed to illustrate the use of a linearized theory, attached flow method for the estimation and optimization of the longitudinal aerodynamic performance of wing-canard and wing-horizontal tail configurations which may employ simple hinged flap systems. Use of an attached flow method is based on the premise that high levels of aerodynamic efficiency require a flow that is as nearly attached as circumstances permit. The results indicate that linearized theory, attached flow, computer code methods (modified to include estimated attainable leading-edge thrust and an approximate representation of vortex forces) provide a rational basis for the estimation and optimization of aerodynamic performance at subsonic speeds below the drag rise Mach number. Generally, good prediction of aerodynamic performance, as measured by the suction parameter, can be expected for near optimum combinations of canard or horizontal tail incidence and leading- and trailing-edge flap deflections at a given lift coefficient (conditions which tend to produce a predominantly attached flow).

Carlson, Harry W.↗

Technical Evaluation of the NASA Model for Cancer Risk to Astronauts Due to Space Radiation

At the request of NASA, the National Research Council's (NRC's) Committee for Evaluation of Space Radiation Cancer Risk Model1 reviewed a number of changes that NASA proposes to make to its model for estimating the risk of radiation-induced cancer in astronauts. The NASA model in current use was last updated in 2005, and the proposed model would incorporate recent research directed at improving the quantification and understanding of the health risks posed by the space radiation environment. NASA's proposed model is defined by the 2011 NASA report Space Radiation Cancer Risk Projections and Uncertainties--2010 . The committee's evaluation is based primarily on this source, which is referred to hereafter as the 2011 NASA report, with mention of specific sections or tables. The overall process for estimating cancer risks due to low linear energy transfer (LET) radiation exposure has been fully described in reports by a number of organizations. The approaches described in the reports from all of these expert groups are quite similar. NASA's proposed space radiation cancer risk assessment model calculates, as its main output, age- and gender-specific risk of exposure-induced death (REID) for use in the estimation of mission and astronaut-specific cancer risk. The model also calculates the associated uncertainties in REID. The general approach for estimating risk and uncertainty in the proposed model is broadly similar to that used for the current (2005) NASA model and is based on recommendations by the National Council on Radiation Protection and Measurements. However, NASA's proposed model has significant changes with respect to the following: the integration of new findings and methods into its components by taking into account newer epidemiological data and analyses, new radiobiological data indicating that quality factors differ for leukemia and solid cancers, an improved method for specifying quality factors in terms of radiation track structure concepts as opposed to the previous approach based on linear energy transfer, the development of a new solar particle event (SPE) model, and the updates to galactic cosmic ray (GCR) and shielding transport models. The newer epidemiological information includes updates to the cancer incidence rates from the life span study (LSS) of the Japanese atomic bomb survivors, transferred to the U.S. population and converted to cancer mortality rates from U.S. population statistics. In addition, the proposed model provides an alternative analysis applicable to lifetime never-smokers (NSs). Details of the uncertainty analysis in the model have also been updated and revised. NASA's proposed model and associated uncertainties are complex in their formulation and as such require a very clear and precise set of descriptions. The committee found the 2011 NASA report challenging to review largely because of the lack of clarity in the model descriptions and derivation of the various parameters used. The committee requested some clarifications from NASA throughout its review and was able to resolve many, but not all, of the ambiguities in the written description.

Source record↗

Assessment of existing analytic methods for prediction of high angle-of-attack loads on delta wings at supersonic speeds

An assessment of the applicability of four loading prediction methods to high angle-of-attack conditions for simplified wing-body configurations is provided. The methods are: The tangent wedge approximation, the linear theory methods of Middleton and Woodward, and a shock-fitting finite-difference technique. Estimates obtained by these methods were compared with experimental pressure data on delta wings to examine the effects of Mach number, camber, sweep angle, and angle of attack. Results indicate that all of the methods provided reasonable estimates at moderate angles of attack. At these moderate angles of attack, the methods of Middleton and Woodward provided good estimates at Mach numbers higher than those usually associated with linear theory. Only the finite-difference method provided reasonable load estimates at high angles of attack.

Landrum, E. J.↗

Water impact analysis of space shuttle solid rocket motor by the finite element method

Preliminary analysis showed that the doubly curved triangular shell elements were too stiff for these shell structures. The doubly curved quadrilateral shell elements were found to give much improved results. A total of six load cases were analyzed in this study. The load cases were either those resulting from a static test using reaction straps to simulate the drop conditions or under assumed hydrodynamic conditions resulting from a drop test. The latter hydrodynamic conditions were obtained through an emperical fit of available data. Results obtained from a linear analysis were found to be consistent with results obtained elsewhere with NASTRAN and BOSOR. The nonlinear analysis showed that the originally assumed loads would result in failure of the shell structures. The nonlinear analysis also showed that it was useful to apply internal pressure as a stabilizing influence on collapse. A final analysis with an updated estimate of load conditions resulted in linear behavior up to full load.

Buyukozturk, O.↗

Gradient methods for identification of distributed parameter systems

Parameter-dependence properties are developed in the context of a linear abstract Cauchy problem governed by a parameter-dependent operator. Of particular interest are problems in which the parameter induces an unbounded perturbation of the evolution operator. Conditions are stated under which the derivative of the state with respect to the parameter possesses certain smoothness properties. These properties lead to local convergence results for a gradient-based estimation algorithm based on quasi-linearization. Numerical results concerning a delay-differential equation indicate the effect of a forcing function on parameter sensitivity.

Brewer, Dennis W.↗

LS47: A DSN station location set compatible with JPL development ephemeris DE108

An updated Deep Space Network station location set, LS47, is presented which is compatible with JPL Development Ephemeris DE108. Analytic procedures for linearly correcting station spin axis and longitude estimates for an ephemeris update based on Brouwer-Clemence Set III parameters are briefly discussed. The validity of this technique is demonstrated by a comparison of a linearly corrected solution with one explicitly determined by reprocessing the data. A mission data base, including Viking 1 and 2 encounter data, is first used to obtain an updated DE96 compatible station location solution, LS46, which in turn is adjusted to form the DE108 solution, LS47. Improved station Z-heights are estimated by using available very long baseline interferometry data. Spin axis differences between LS46 and LS47 are relatively insignificant; however, the ephemeris change introduces a -0.8 x 10 to the -5 power degree rotation in the DE96 longitude ephemeris.

Ellis, J.↗

Parameterized Linear Longitudinal Airship Model

A parameterized linear mathematical model of the longitudinal dynamics of an airship is undergoing development. This model is intended to be used in designing control systems for future airships that would operate in the atmospheres of Earth and remote planets. Heretofore, the development of linearized models of the longitudinal dynamics of airships has been costly in that it has been necessary to perform extensive flight testing and to use system-identification techniques to construct models that fit the flight-test data. The present model is a generic one that can be relatively easily specialized to approximate the dynamics of specific airships at specific operating points, without need for further system identification, and with significantly less flight testing. The approach taken in the present development is to merge the linearized dynamical equations of an airship with techniques for estimation of aircraft stability derivatives, and to thereby make it possible to construct a linearized dynamical model of the longitudinal dynamics of a specific airship from geometric and aerodynamic data pertaining to that airship. (It is also planned to develop a model of the lateral dynamics by use of the same methods.) All of the aerodynamic data needed to construct the model of a specific airship can be obtained from wind-tunnel testing and computational fluid dynamics

Kulczycki, Eric↗

Theoretical Characterization of C{sub 6}N, C{sub 6}N{sup −}, and C{sub 6}N{sup +}

This work emphasizes the stability of C{sub 6}N linear carbon chains and carbon clusters containing three-body rings. C{sub 6}N possesses at least 44 neutral isomers and 38 and 35 isomers with a negative or a positive charge. The lowest-energy structures, which can be candidates for laboratory and astrophysical detection, were studied with RCCSD(T)-F12 and MRCI/CASSCF, specifying properties for various electronic states. Neutral C{sub 6}N displays two prominent equilibrium structures, a nitrogen-terminated linear form (X{sup 2}Π) and a C{sub 2v} form (X{sup 2}B{sub 2}) containing a three-carbon ring. They are separated by 0.21 eV. For the linear one, Renner–Teller and spin–orbit effects are expected. Its equilibrium spin–orbit constant ∣ A{sub so,e} ∣was predicted to be 29.09 cm{sup −1}, the dipole moment was computed to be 0.2059 D, and B{sub e} was computed to be 868.20 MHz, very close to the experimental value (B{sub 0}(exp)-B{sub e}(calc) ∼ 4.9 MHz). The anion presents two stable forms, a linear triplet (X{sup 3}Σ{sup −}) and a C{sub 2v}(X{sup 1}A{sub 1}) isomer, whereas for the cation most stable structures are both linear. For the linear anion and the two linear cations B{sub 0} was estimated to be 870.05 MHz, 875.72 MHz, and 917.36 MHz and μ was estimated to be 4.4200 D, 4.2013D, and 0.7828D. For the neutral C{sub 2V} form, Ae = 45048.74 MHz, Be = 1062.78 MHz, Ce = 1038.30 MHz, and μ = 0.8840 D, whereas for the negative C{sub 2V} form, Ae = 40087.48 MHz, Be = 1067.35 MHz, Ce = 1039.67 MHz, and μ = 6.2854 D. Below the electron affinity of the most stable isomer (EA = −3.42 eV), the linear anion displays three probably existing electronic states. Detectability is discussed in terms of the symmetry and spin multiplicity of the ground electronic states.

79 ASTRONOMY AND ASTROPHYSICS↗

Compact Polarimetry in a Low Frequency Spaceborne Context

Compact polarimetry has been shown to be an interesting alternative mode to full polarimetry when global coverage and revisit time are key issues. It consists on transmitting a single polarization, while receiving on two. Several critical points have been identified, one being the Faraday rotation (FR) correction and the other the calibration. When a low frequency electromagnetic wave travels through the ionosphere, it undergoes a rotation of the polarization plane about the radar line of sight for a linearly polarized wave, and a simple phase shift for a circularly polarized wave. In a low frequency radar, the only possible choice of the transmit polarization is the circular one, in order to guaranty that the scattering element on the ground is illuminated with a constant polarization independently of the ionosphere state. This will allow meaningful time series analysis, interferometry as long as the Faraday rotation effect is corrected for the return path. In full-polarimetric (FP) mode, two techniques allow to estimate the FR: Freeman method using linearly polarized data, and Bickel and Bates theory based on the transformation of the measured scattering matrix to a circular basis. In CP mode, an alternate procedure is presented which relies on the bare surface scattering properties. These bare surfaces are selected by the conformity coefficient, invariant with FR. This coefficient is compared to other published classifications to show its potential in distinguishing three different scattering types: surface, doublebounce and volume. The performances of the bare surfaces selection and FR estimation are evaluated on PALSAR and airborne data. Once the bare surfaces are selected and Faraday angle estimated over them, the correction can be applied over the whole scene. The algorithm is compared with both FP techniques. In the last part of the paper, the calibration of a CP system from the point of view of classical matrix transformation methods in polarimetry is proposed.

radar polarization diversity↗

Estimation of abatement potentials and costs of air pollution emissions in China

Understanding the air pollution emission abatement potential and associated control cost is a prerequisite to design cost ef?cient control policies. In this study, a linear programming algorithm model, International Control Cost Estimate Tool, was updated with cost data for applications of 56 types of end-of-pipe technologies and ?ve types of renewable energy in 10 major sectors namely power generation, industry combustion, cement pro-duction, iron and steel production, other industry processes, domestic combustion, transportation, solvent use, livestock rearing, and fertilizer use. The updated model was implemented to estimate the abatement potential and marginal cost of multiple pollutants in China. The total maximum abatement potentials of sulfur dioxide (SO2), nitrogen oxides (NOx), primary particulate matter (PM2.5), non-volatile organic compounds (NMVOCs), and ammonia (NH3) in China were estimated to be 19.2, 20.8, 9.1, 17.2 and 8.6 Mt, respectively, which accounted for 89.7%, 89.9%, 94.6%, 74.0%, and 80.2% of their total emissions in 2014, respectively. The associated control cost of such reductions was estimated as 92.5, 469.7, 75.7, 449.0, and 361.8 billion CNY in SO2, NOx, primary PM2.5, NMVOCs and NH3, respectively. Shandong, Jiangsu, Henan, Zhejiang, and Guangdong provinces exhibited large abatement potentials for all pollutants. Provincial disparity analysis shows that high GDP regions tend to have higher reduction potential and total abatement costs. End-of-pipe technologies tended be a cost-ef?cient way to control pollution in industries processes (i.e., cement plants, iron and steel plants, lime production, building ceramic production, glass and brick production), whereas such technologies were less cost- effective in fossil fuel-related sectors (i.e., power plants, industry combustion, domestic combustion, and transportation) compared with renewable energy. The abatement potentials and marginal abatement cost curves developed in this study can further be used as a crucial component in an integrated model to design optimized cost-ef?cient control policies.

Zhang, Fenfen↗

Model-Based Engine Control Architecture with an Extended Kalman Filter

This paper discusses the design and implementation of an extended Kalman filter (EKF) for model-based engine control (MBEC). Previously proposed MBEC architectures feature an optimal tuner Kalman Filter (OTKF) to produce estimates of both unmeasured engine parameters and estimates for the health of the engine. The success of this approach relies on the accuracy of the linear model and the ability of the optimal tuner to update its tuner estimates based on only a few sensors. Advances in computer processing are making it possible to replace the piece-wise linear model, developed off-line, with an on-board nonlinear model running in real-time. This will reduce the estimation errors associated with the linearization process, and is typically referred to as an extended Kalman filter. The non-linear extended Kalman filter approach is applied to the Commercial Modular Aero-Propulsion System Simulation 40,000 (C-MAPSS40k) and compared to the previously proposed MBEC architecture. The results show that the EKF reduces the estimation error, especially during transient operation.

engine control↗

Input-output characterization of an ultrasonic testing system by digital signal analysis

Ultrasonic test system input-output characteristics were investigated by directly coupling the transmitting and receiving transducers face to face without a test specimen. Some of the fundamentals of digital signal processing were summarized. Input and output signals were digitized by using a digital oscilloscope, and the digitized data were processed in a microcomputer by using digital signal-processing techniques. The continuous-time test system was modeled as a discrete-time, linear, shift-invariant system. In estimating the unit-sample response and frequency response of the discrete-time system, it was necessary to use digital filtering to remove low-amplitude noise, which interfered with deconvolution calculations. A digital bandpass filter constructed with the assistance of a Blackman window and a rectangular time window were used. Approximations of the impulse response and the frequency response of the continuous-time test system were obtained by linearly interpolating the defining points of the unit-sample response and the frequency response of the discrete-time system. The test system behaved as a linear-phase bandpass filter in the frequency range 0.6 to 2.3 MHz. These frequencies were selected in accordance with the criterion that they were 6 dB below the maximum peak of the amplitude of the frequency response. The output of the system to various inputs was predicted and the results were compared with the corresponding measurements on the system.

Williams, J. H., Jr.↗