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 1,189 records · Page 66

WE-D-204-06: An Open Source ImageJ CatPhan Analysis Tool

Purpose: The CatPhan is a popular QA device for assessing CT image quality. There are a number of software options which perform analysis of the CatPhan. However, there is often little ability for the user to adjust the analysis if it isn’t running properly, and these are all expensive options. An open source tool is an effective solution. Methods: To use the software, the user imports the CT as an image sequence in ImageJ. The user then scrolls to the slice with the lateral dots. The user then runs the plugin. If tolerance constraints are not already created, the user is prompted to enter them or to use generic tolerances. Upon completion of the analysis, the plugin calls pdfLaTex to compile the pdf report. There is a csv version of the report as well. A log of the results from all CatPhan scans is kept as a csv file. The user can use this to baseline the machine. Results: The tool is capable of detecting the orientation of the phantom. If the CatPhan was scanned backwards, one can simply flip the stack of images horizontally and proceed with the analysis. The analysis includes Sensitometry (estimating the effective beam energy), HU values and linearity, Low Contrast Visibility (using LDPE & Polystyrene), Contrast Scale, Geometric Accuracy, Slice Thickness Accuracy, Spatial resolution (giving the MTF using the line pairs as well as the point spread function), CNR, Low Contrast Detectability (including the raw data), Uniformity (including the Cupping Effect). Conclusion: This is a robust tool that analyzes more components of the CatPhan than other software options (with the exception of ImageOwl). It produces an elegant pdf and keeps a log of analyses for long-term tracking of the system. Because it is open source, users are able to customize any component of it.

60 APPLIED LIFE SCIENCES↗

Scaled Vecchia Approximation for Fast Computer-Model Emulation

Many scientific phenomena are studied using computer experiments consisting of multiple runs of a computer model while varying the input settings. Gaussian processes (GPs) are a popular tool for the analysis of computer experiments, enabling interpolation between input settings, but direct GP inference is computationally infeasible for large datasets. We adapt and extend a powerful class of GP methods from spatial statistics to enable the scalable analysis and emulation of large computer experiments. Specifically, we apply Vecchia’s ordered conditional approximation in a transformed input space, with each input scaled according to how strongly it relates to the computer-model response. The scaling is learned from the data by estimating parameters in the GP covariance function using Fisher scoring. Our methods are highly scalable, enabling estimation, joint prediction, and simulation in near-linear time in the number of model runs. In several numerical examples, our approach substantially outperformed existing methods.

97 MATHEMATICS AND COMPUTING↗

Neural Networks to Find the Optimal Forcing for Offsetting the Anthropogenic Climate Change Effects

Abstract Of great relevance to climate engineering is the systematic relationship between the radiative forcing to the climate system and the response of the system, a relationship often represented by the linear response function (LRF) of the system. However, estimating the LRF often becomes an ill-posed inverse problem due to high-dimensionality and nonunique relationships between the forcing and response. Recent advances in machine learning make it possible to address the ill-posed inverse problem through regularization and sparse system fitting. Here, we develop a convolutional neural network (CNN) for regularized inversion. The CNN is trained using the surface temperature responses from a set of Green’s function perturbation experiments as imagery input data together with data sample densification. The resulting CNN model can infer the forcing pattern responsible for the temperature response from out-of-sample forcing scenarios. This promising proof of concept suggests a possible strategy for estimating the optimal forcing to negate certain undesirable effects of climate change. The limited success of this effort underscores the challenges of solving an inverse problem for a climate system with inherent nonlinearity. Significance Statement Predicting the climate response for a given climate forcing is a direct problem, while inferring the forcing for a given desired climate response is often an inverse, ill-posed, problem, posing a new challenge to the climate community. This study makes the first attempt to infer the radiative forcing for a given target pattern of global surface temperature response using a deep learning approach. The resulting deeply trained convolutional neural network inversion model shows promise in capturing the forcing pattern corresponding to a given surface temperature response, with a significant implication on the design of an optimal solar radiation management strategy for curbing global warming. This study also highlights the technical challenges that future research should prioritize in seeking feasible solutions to the inverse climate problem.

Ren, Huiying↗

Non-Linear Langevin and Fractional Fokker–Planck Equations for Anomalous Diffusion by Lévy Stable Processes

The numerical solutions to a non-linear Fractional Fokker–Planck (FFP) equation are studied estimating the generalized diffusion coefficients. The aim is to model anomalous diffusion using an FFP description with fractional velocity derivatives and Langevin dynamics where Lévy fluctuations are introduced to model the effect of non-local transport due to fractional diffusion in velocity space. Distribution functions are found using numerical means for varying degrees of fractionality of the stable Lévy distribution as solutions to the FFP equation. The statistical properties of the distribution functions are assessed by a generalized normalized expectation measure and entropy and modified transport coefficient. The transport coefficient significantly increases with decreasing fractality which is corroborated by analysis of experimental data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Determination of Constraint-Independent Crack Tip Opening Angle for Stable Crack Growth in High-Strength Ductile Steels

The crack tip opening angle (CTOA) is one of fracture toughness parameters that has been used for decades in describing large stable crack growth in thin-walled aerospace structures under the low-constraint conditions. Recently, the pipeline industry has developed a growing interest in using the CTOA parameter to serve as the minimum required fracture toughness to arrest dynamic crack propagation in modern gas transmission pipelines made of high-strength ductile steel. To meet this industrial need, the CTOA test standard ASTM E3039 was therefore developed for measuring the constant critical CTOA. ASTM E3039 recommends a drop weight tearing test (DWTT) specimen with a shallow crack for standard CTOA testing, but its CTOA may depend on the low constraint condition of the DWTT specimen at the crack tip. Verifying the constraint independence of the DWTT-measured CTOA thus becomes indispensable for applying CTOA toughness to the running fracture control in the pipeline design. For this purpose, the present paper evaluates critical CTOA values in a set of fracture toughness tests on single-edge notched bend (SENB) specimens with shallow and deep cracks, based on four CTOA estimation models. Among these, the Ln(P)-LLD linear fit model is similar to that recommended by ASTM E3039 for CTOA calculation. Fracture test data for X80 pipeline steel and HY80 structural steel were considered in the CTOA evaluation. The results showed that the four CTOA models were able to determine a constraint independent CTOA value for stable crack growth in the SENB specimens. As a result, a single, reliable, constant CTOA value could be determined regardless of the specimen geometry or the crack-tip constraint conditions. Therefore, the CTOA measured using ASTM E3039 is constraint-independent and transferable to use in cases of actual cracks propagating in gas transmission pipelines.

36 MATERIALS SCIENCE↗

Conditional H I Mass Functions and the H I -to-halo Mass Relation in the Local Universe

We present a new H I mass estimator that relates log 10 (M HI /M * ) to a linear combination of four galaxy properties: stellar surface mass density, color index u - r, stellar mass, and concentration index, with the scatter of individual galaxies around the mean H I mass modeled with a Gaussian distribution function. We calibrate the estimator using the xGASS sample, including both H I detection and nondetection, and constrain the model parameters through Bayesian inferences. Tests with mock catalogs demonstrate that our estimator provides unbiased H I masses for optical samples like SDSS. We apply our estimator to the SDSS spectroscopic sample to estimate the H I mass function (HIMF) of local galaxies, as well as the conditional H I mass function in galaxy groups and the H I –halo mass relation. Our HIMF agrees with the ALFALFA measurements at M HI ≳ 5 × 10 9 M ⊙ , but with higher amplitude and a steeper slope at lower masses. We show that this discrepancy is caused primarily by the cosmic variance, which is corrected for the SDSS sample but not for ALFALFA. The total CHIMFs for all halo masses can be described by a single Schechter function, while those of central galaxies show a double-Gaussian profile. The total H I mass in a group increases monotonically with halo mass, but for central galaxies, the H I mass shows weak dependence on halo mass when M h ≳ 10 12 M ⊙ . The observed H I –halo mass relation is not reproduced by current hydrodynamic simulations and semianalytic models of galaxy formation.

79 ASTRONOMY AND ASTROPHYSICS↗

Estimate of the Total Energy of the Beh(sub 2) Molecule

The existence of gaseous BeH 2 , predicted from simple valence theory, has still not been established experimentally. We have estimated the total energy of this molecule, assumed linear symmetric, as the sum of four negative terms; the calculated SCF energy, the SCF error to the Hartree-Fock limit, the correlation energy to the exact nonrelativistic energy, and the relativistic energy.

Molecular energy↗

Polarization tracking of a partially coherent signal using a double loop

We consider the problem of simultaneously estimating the phase and polarization angles of a linearly polarized electromagnetic wave. The solution of this problem has applications to space communications and ionospheric probing. The nonlinear stochastic characteristics of a cross-coupled double-loop implementation for the tracker are evaluated using the multidimensional Fokker-Planck equation. The exact stationary probability density is found for a special case of equal signal-to-noise ratio (SNR) in the two loops. The stationary probability density is obtained for a general class of double loops when the two loop bandwidths differ greatly.

Ohlson, J. E.↗

Computer Mapping of Water Quality in Saginaw Bay with LANDSAT Digital Data

The author has identified the following significant results. LANDSAT digital data and ground truth measurements for Saginaw Bay (Lake Huron), Michigan, for 31 July 1975 were correlated by stepwise linear regression and the resulting equations used to estimate invisible water quality parameters in nonsampled areas. Chloride, conductivity, total Kjeldahl nitrogen, total phosphorus, and chlorophyll a were best correlated with the ratio of LANDSAT Band 4 to Band 5. Temperature and Secchi depth correlate best with Band 5.

Rogers, R. H.↗

Full-scale study of the cooling system aerodynamics of an operating piston engine installed in a light aircraft wing panel

Measurements of the drag and of the nacelle internal pressures on a wing and nacelle that housed a horizontally opposed piston engine were made in the 40- by 80-Foot Wind Tunnel at Ames Research Center. These tests are follow-ons to earlier tests made with the same wing and nacelle but in which the engine was replaced with an electric motor and an adjustable orifice plate. In the initial tests the orifice plate was used to control the rate of cooling-air flow through the nacelle and thereby to simulate a range of gasoline engine types. Good agreement was found between the results of those tests and of the test reported here. Also, the upper and lower plenum pressure and cooling-air flow rate were found to be related by conventional equations used to represent the flow through orifices. Tests were run with three cooling air inlet sizes over a free-stream velocity range from 50 to 150 knots, an angle of attack range from 0 deg to 10 deg, and a cowl-flap deflection range from 0 deg to 30 deg. The data were analyzed by computing a flow coefficient similar to that used in the analysis of orifices. It was found that all of the flow coefficient values fell within a band that varied linearly with inlet area. The linear mean line through this band provides an estimate of the relationship between cooling-air flow rate and upper plenum pressure over a wide range of test conditions.

Corsiglia, V. R.↗

The role of modern control theory in the design of controls for aircraft turbine engines

The development, applications, and current research in modern control theory (MCT) are reviewed, noting the importance for fuel-efficient operation of turbines with variable inlet guide vanes, compressor stators, and exhaust nozzle area. The evolution of multivariable propulsion control design is examined, noting a basis in a matrix formulation of the differential equations defining the process, leading to state space formulations. Reports and papers which appeared from 1970-1982 which dealt with problems in MCT applications to turbine engine control design are outlined, including works on linear quadratic regulator methods, frequency domain methods, identification, estimation, and model reduction, detection, isolation, and accommodation, and state space control, adaptive control, and optimization approaches. Finally, NASA programs in frequency domain design, sensor failure detection, computer-aided control design, and plant modeling are explored

Zeller, J.↗

Options for a new Arizona observatory - Configurations, costs and site

The process of creating a major observatory, with emphasis on the configurations, costs, and site involved, is examined. The premise that cost optimization starts with a modestly priced building block, such as a 7.5-m telescope, is considered. The telescope unit and the configurations in which it is used, ranging from a single telescope or a set of separate telescopes, to an MMT in which the light of all of the telescopes can be phased when desired, are described. It is shown that options exist for using elements in a nonredundant linear array or Michelson interferometer. A means of estimating costs of different facilities, all on the same basis, is provided. The importance of placing the observatory on a site, such as Mt. Graham in Arizona, which maximizes the output per unit cost of the facility is discussed.

Woolf, N. J.↗

Configuration trade and code validation study on a conical hypersonic vehicle

A study has been conducted on a generic wing-cone transatmospheric vehicle at Mach numbers form 2.5 to 4.5. The objectives of the study were to experimentally define the aerodynamic characteristics of the vehicle and evaluate several computational aerodynamic prediction methods through comparison with the experimental results. The baseline wing-cone configuration fuselage consisted of a 5 deg half-angle cone forebody, cylindrical midbody, and 9 deg truncated cone afterbody. The 4-percent-thick diamond airfoil wing had an aspect ratio of 1. Several configuration variables were investigated to provide trade information on canard, wing-position and incidence, vertical tail, and nose bluntness effects. Results of the study show that wing-position and wing-incidence effects on the longitudinal aerodynamic characteristics can be significantly influenced by wing-body interference. The use of positive wing incidence to provide favorable forebody orientation for possible inlet performance improvement is accompanied by trim drag and lift-drag ratio penalties. The lateral-directional stability characteristics were strongly influenced by the location of the vertical tails. The higher-order full-potential method provided better estimates of the aerodynamic characteristics than either the linearized supersonic potential method or the tangent-cone/tangent-wedge/shock-expansion on method.

Covell, Peter F.↗

Approximation theory for LQG (Linear-Quadratic-Gaussian) optimal control of flexible structures

An approximation theory is presented for the LQG (Linear-Quadratic-Gaussian) optimal control problem for flexible structures whose distributed models have bounded input and output operators. The main purpose of the theory is to guide the design of finite dimensional compensators that approximate closely the optimal compensator. The optimal LQG problem separates into an optimal linear-quadratic regulator problem and an optimal state estimation problem. The solution of the former problem lies in the solution to an infinite dimensional Riccati operator equation. The approximation scheme approximates the infinite dimensional LQG problem with a sequence of finite dimensional LQG problems defined for a sequence of finite dimensional, usually finite element or modal, approximations of the distributed model of the structure. Two Riccati matrix equations determine the solution to each approximating problem. The finite dimensional equations for numerical approximation are developed, including formulas for converting matrix control and estimator gains to their functional representation to allow comparison of gains based on different orders of approximation. Convergence of the approximating control and estimator gains and of the corresponding finite dimensional compensators is studied. Also, convergence and stability of the closed-loop systems produced with the finite dimensional compensators are discussed. The convergence theory is based on the convergence of the solutions of the finite dimensional Riccati equations to the solutions of the infinite dimensional Riccati equations. A numerical example with a flexible beam, a rotating rigid body, and a lumped mass is given.

Gibson, J. S.↗

Contamination effects of GPS Navstar solar array performance

The solar arrays on Navstars 1 to 6 have been limiting the electrical power capabilities in the extended life (beyond the 5 year design life). The departure from predicted performance consists of an extra 2.5 percent per year degradation beyond the radiation model estimates. This degradation is unusual in showing a linear rather than exponential decay with time. The performance of the arrays on these satellites has been examined in order to predict future behavior and to make refined projections on the Navstar 7 to 11 solar arrays. Evidence obtained from flight experiments on Navstar 5 and 6, and from laboratory experiments, suggests that contamination of the solar arrays while on orbit may be responsible. The evidence for photo-induced contamination of spacecraft surfaces is presented, and the effect on solar array output in the case of the GPS satellites is shown to be consistent with the observed anomalies.

Marvin, Dean C.↗

Estimation of motion parameters for a rigid body from its orthogonal projection

An estimate is presented of the motion parameters, namely, linear and angular velocities of a rigid body rotating and translating in three-dimensional-space. It is assumed that the velocities are constant and that only the orthogonal projection of the motion is observable. In particular, if (x, y, z) is the Cartesian coordinate, it is assumed that the projection of the motion on the x-y plane is observed and the information along the z coordinate is lost.

Ganguly, S.↗

Fusion techniques using distributed Kalman filtering for detecting changes in systems

A comparison is made of the performances of two detection strategies that are based on different data fusion techniques. The strategies detect changes in a linear system. One detection strategy involves combining the estimates and error covariance matrices of distributed Kalman filters, generating a residual from the used estimates, comparing this residual to a threshold, and making a decision. The other detection strategy involves a distributed decision process in which estimates from distributed Kalman filters are used to generate distributed residuals which are compared locally to a threshold. Local decisions are made and these decisions are then fused into a global decision. The performances of each of these detection schemes are compared, and it is concluded that better performance is achieved when local decisions are made and then fused into a global decision.

Belcastro, Celeste M.↗