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,027 records · Page 57

A Reduced Dimension Static, Linearized Kalman Filter and Smoother

An approximate Kalman filter and smoother, based on approximations of the state estimation error covariance matrix, is described. Approximations include a reduction of the effective state dimension, use of a static asymptotic error limit, and a time-invariant linearization of the dynamic model for error integration. The approximations lead to dramatic computational savings in applying estimation theory to large complex systems. Examples of use come from TOPEX/POSEIDON.

Kalman filter smoother modeling ocean simulation↗

Analysis of the Feasibility of UAS-Based EMI Sensing for Underground Utilities Detection and Mapping

This paper investigates the feasibility of using a linear current sensing (LCS) technique integrated on an unmanned aerial system (UAS) for detecting and mapping underground infrastructure rapidly and cost-effectively. The LCS technique is based on data from a wide band of electromagnetic induction frequencies (50 kHz to 2 MHz) using a vector magnetic field gradiometer. This technique takes advantage of a slowly decaying secondary magnetic field in order to achieve greater standoff detection distance ($\frac{1}{R^2}$ vs. $\frac{1}{R^6}$ for compact metallic targets during EMI sensing, where R is the distance from a target to the sensor). These secondary magnetic fields are produced by the excite current on long conductors, allowing detection at a distance of 10 meters or more. The system operates between tens of kHz to a few MHz and uses either an active EMI source or existing EM fields to excite this linear current on a long metallic subsurface target. Once excited, these linear currents produce a secondary magnetic field that is detected with an above ground triaxial magnetic field gradiometer. By moving and tracking its geolocation, the system outputs rich datasets sufficient to support high-fidelity forward and inverse EMI models for estimating the depth and orientation of deep underground long linear metallic infrastructure. The system’s hardware and its integration to a UAS system are outlined, along with the formulation of LCS theory, and numerical and experimental data are presented. The results illustrate that the LCS technique offers large standoff detection, is adaptable to UAS, and could be used effectively for detecting deep underground infrastructure such as wires and pipes.

54 ENVIRONMENTAL SCIENCES↗

First-Order System Least Squares for the Stokes Equations, with Application to Linear Elasticity

Following our earlier work on general second-order scalar equations, here we develop a least-squares functional for the two- and three-dimensional Stokes equations, generalized slightly by allowing a pressure term in the continuity equation. By introducing a velocity flux variable and associated curl and trace equations, we are able to establish ellipticity in an H(exp 1) product norm appropriately weighted by the Reynolds number. This immediately yields optimal discretization error estimates for finite element spaces in this norm and optimal algebraic convergence estimates for multiplicative and additive multigrid methods applied to the resulting discrete systems. Both estimates are uniform in the Reynolds number. Moreover, our pressure-perturbed form of the generalized Stokes equations allows us to develop an analogous result for the Dirichlet problem for linear elasticity with estimates that are uniform in the Lame constants.

Cai, Z.↗

Estimation and Analysis of Nonlinear Stochastic Systems

The algebraic and geometric structures of certain classes of nonlinear stochastic systems were exploited in order to obtain useful stability and estimation results. The class of bilinear stochastic systems (or linear systems with multiplicative noise) was discussed. The stochastic stability of bilinear systems driven by colored noise was considered. Approximate methods for obtaining sufficient conditions for the stochastic stability of bilinear systems evolving on general Lie groups were discussed. Two classes of estimation problems involving bilinear systems were considered. It was proved that, for systems described by certain types of Volterra series expansions or by certain bilinear equations evolving on nilpotent or solvable Lie groups, the optimal conditional mean estimator consists of a finite dimensional nonlinear set of equations. The theory of harmonic analysis was used to derive suboptimal estimators for bilinear systems driven by white noise which evolve on compact Lie groups or homogeneous spaces.

Marcus, S. I.↗

Measuring Epoxy-Curing Kinetics

Key reaction parameters estimated from single run. Single DSC Curve used to construct linearized plot from which kinetic reaction parameters are obtained: vertical axis intercept proportional to activation energy, while slope equals reaction order. DSC methods show promise for rapid screening and estimation of cure kinetic parameters of epoxy resins.

Cizmecioglu, M.↗

On the Relationship Between Global Land-Ocean Temperature and Various Descriptors of Solar-Geomagnetic Activity and Climate

Examined are sunspot cycle- (SC-) length averages of the annual January-December values of the Global Land-Ocean Temperature Index ( ) in relation to SC-length averages of annual values of various descriptors of solar-geomagnetic activity and climate, incorporating lags of 0-5 yr. For the overall interval SC12-SC23, the is inferred to correlate best against the parameter incorporating lag = 5 yr, where the parameter refers to the resultant aa value having removed that portion of the annual aa average value due to the yearly variation of sunspot number (SSN). The inferred correlation between the and is statistically important at confidence level cl > 99.9%, having a coefficient of linear correlation r = 0.865 and standard error of estimate se = 0.149 degC. Excluding the most recent cycles SC22 and SC23, the inferred correlation is stronger, having r = 0.969 and se = 0.048 degC. With respect to the overall trend in the , which has been upwards towards warmer temperatures since SC12 (1878-1888), solar-geomagnetic activity parameters are now trending downwards (since SC19). For SC20-SC23, in contrast, comparison of the against SC-length averages of the annual value of the Mauna Loa carbon dioxide ( ) index is found to be highly statistically important (cl >> 99.9%), having r = 0.9994 and se = 0.012 degC for lag = 2 yr. On the basis of the inferred preferential linear correlation between the and , the current ongoing SC24 is inferred to have warmer than was seen in SC23 (i.e., >0.526 degC), probably in excess of 0.68 degC (relative to the 1951-1980 base period).

Wilson, Robert M.↗

A modal parameter extraction procedure applicable to linear time-invariant dynamic systems

Modal analysis has emerged as a valuable tool in many phases of the engineering design process. Complex vibration and acoustic problems in new designs can often be remedied through use of the method. Moreover, the technique has been used to enhance the conceptual understanding of structures by serving to verify analytical models. A new modal parameter estimation procedure is presented. The technique is applicable to linear, time-invariant systems and accommodates multiple input excitations. In order to provide a background for the derivation of the method, some modal parameter extraction procedures currently in use are described. Key features implemented in the new technique are elaborated upon.

Kurdila, A. J.↗

Identification of motion parameters of a rigid body from its orthogonal and perspective projections

An estimate is made of the motion parameters, namely, linear and angular velocities, of a rigid body rotating and translating in three-space. The authors assume that the velocities are constant and that the motion is not completely observable. They consider two separate cases of partial observations corresponding to the orthogonal and the perspective projections, respectively. If (x, y, z) is the Cartesian coordinate of the three-space, the authors assume in the first case that the projection of the motion on the x-y plane is observed. If (r, theta, phi) is the polar coordinates of the three-space, they assume in the second case that the parameter vector (theta, phi) is observed. The use of both of these cases to estimate the motion parameters is discussed.

Ganguly, S.↗

Systems and methods for direct estimation of battery parameters using only charge/discharge curves

Electrochemical models for the lithium-ion battery are useful in predicting and controlling its performance. The values of the parameters in these models are vital to their accuracy. However, not all parameters can be measured precisely, especially when destructive methods are prohibited. In some embodiments of the present disclosure, a parameter estimation approach is used to estimate the open circuit potential of the positive electrode (Up) using piecewise linear approximation together with all the other parameters of a single particle model. Up and 10 more parameters may be estimated from a single discharge curve without knowledge of the electrode chemistry using a technique such as a genetic algorithm. Different case studies were presented for estimating Up with different types of parameters of the battery model. The estimated parameters were then validated by comparing simulations at different discharge rates with experimental data.

Qi, Yanbo↗

The Cauchy Problem for Boltzmann Bi-linear Systems: The Mixing of Monatomic and Polyatomic Gases

Abstract From a unified vision of vector valued solutions in weighted Banach spaces, this paper establishes the existence and uniqueness for space homogeneous Boltzmann bi-linear systems with conservative collisional forms arising in complex gas dynamical structures. This broader vision is directly applied to dilute multi-component gas mixtures composed of both monatomic and polyatomic gases. Such models can be viewed as extensions of scalar Boltzmann binary elastic flows, as much as monatomic gas mixtures with disparate masses and single polyatomic gases, providing a unified approach for vector valued solutions in weighted Banach spaces. Novel aspects of this work include developing the extension of a general ODE theory in vector valued weighted Banach spaces, precise lower bounds for the collision frequency in terms of the weighted Banach norm, energy identities, angular or compact manifold averaging lemmas which provide coerciveness resulting into global in time stability, a new combinatorics estimate for p -binomial forms producing sharper estimates for the k -moments of bi-linear collisional forms. These techniques enable the Cauchy problem improvement that resolves the model with initial data corresponding to strictly positive and bounded initial vector valued mass and total energy, in addition to only a $$2^+$$ 2 + moment determined by the hard potential rates discrepancy, a result comparable in generality to the classical Cauchy theory of the scalar homogeneous Boltzmann equation.

Physics↗

SLR Station Recovery, Center of Frame Motion, and Time Varying Gravity

Weekly station position estimates, beginning with 1993, are derived from the ITRF2008-based SLR processing of up to four satellites: Lageos 1, Lageos2, Starlette, and Stella. Helmert parameters obtained from c omparison of weekly SLR station positions and the a-priori SLRF2008 station complement are evaluated for geocenter motion and scale. Two me thods for modeling time varying gravity are employed in the SLR satel lite POD processing, with GGM03S serving as the static gravity field. Both methods forward model atmosphere gravity derived from 6-hour ECM WF pressure data. The standard approach applies an annual 20x20 field estimated from 4 years of GRACE data, and the IERS2003 recommended linear rates for C20, C30, C40, C21, and S21. The alternate approach us es a new set of low-order/degree 4x4 coefficients estimated weekly fr om SLR & DORIS processing to 10 satellites from 1993-2012. This exper imental tvg4x4 model has been shown to improve the TOPEX, Jason-1, and Jason-2 altimeter satellite orbits,. In this paper we apply the more detailed time-variable gravity modeling to the SLR satellite POD pro cessing and subsequent reference frame analyses. For this study we will evaluate the orbit differences (periodic and secular) for the satel lites concerned, characterize the impact on the station coordinate solutions, and the impact on reference frame parameters (geocenter and s cale).

Zelensky, Nikita P.↗

Stability analysis of WEST L-mode discharges with improved confinement from boron powder injection

WEST L-mode plasmas with dominant electron heating and no core torque source have observed improvements in confinement during boron (B) powder injection. These results are reminiscent of previous powder injection experiments on other devices and gaseous impurity seeding experiments on WEST. During powder injection, the stored energy increased up to 25% due to enhanced ion and electron heat and particle confinement. The improvements in confinement were not indicative of an L-H transition. To identify the dominant mechanisms and the causality chain behind these improvements in confinement, we employ interpretative modeling using METIS, predictive integrated modeling using a high-fidelity plasma simulator (HFPS), and stand-alone gyrokinetic simulations using quasi-linear gyrokinetic code. Interpretative modeling with METIS allowed for the estimation of missing data while maintaining good overall consistency with experiment. These results provided the initial conditions for fully predictive flux driven simulations using the HFPS. From these simulations, quasi-linear gyrokinetic analysis was performed at ρ = 0.5 and ρ = 0.65. Collisionality was found to be a strong candidate for the turbulence suppression mechanism at ρ = 0.5, while a combination of collisionality and the T e /T i ratio was found to be the likely mechanism at ρ = 0.65. The results additionally suggested that increased Z eff (through main ion dilution) could play a role in the improved confinement, but this could not be confirmed due to a lack of experimental measurements. The modeling framework established here can now be used to evaluate and exploit a variety of future powder injection experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A multiple linear regression analysis of hot corrosion attack on a series of nickel base turbine alloys

Multiple linear regression analysis was used to determine an equation for estimating hot corrosion attack for a series of Ni base cast turbine alloys. The U transform (i.e., 1/sin (% A/100) to the 1/2) was shown to give the best estimate of the dependent variable, y. A complete second degree equation is described for the centered" weight chemistries for the elements Cr, Al, Ti, Mo, W, Cb, Ta, and Co. In addition linear terms for the minor elements C, B, and Zr were added for a basic 47 term equation. The best reduced equation was determined by the stepwise selection method with essentially 13 terms. The Cr term was found to be the most important accounting for 60 percent of the explained variability hot corrosion attack.

Barrett, C. A.↗

Impact damage of composites

A model is described for estimating the impact damage of fiber reinforced composite plates. The displacements and stresses are calculated by a three dimensional transient, finite element method of solution of the governing equations applicable to a linearly elastic body. The region in which damage occurs is estimated using the Tsai-Wu failure criterion. A computer code was developed which can be used to calculate the impact force, displacements and velocities of the plate and the impact body, stresses and strains in the plate, and the damage area. Sample numerical results are presented illustrating the type of information provided by the code. Comparisons between measured and calculated damage areas are also given.

Wu, Hsi-Young T.↗

Inverse rendering of fusion plasmas: Inferring plasma composition from imaging systems

In this work, we develop a differentiable rendering pipeline for visualising plasma emission within tokamaks, and estimating the gradients of the emission and estimating other physical quantities. Unlike prior work, we are able to leverage arbitrary representations of plasma quantities and easily incorporate them into a non-linear optimisation framework. The efficiency of our method enables not only estimation of a physically plausible image of plasma, but also recovery of the neutral Deuterium distribution from imaging and midplane measurements alone. We demonstrate our method with three different levels of complexity showing first that a poloidal neutrals density distribution can be recovered from imaging alone, second that the distributions of neutral Deuterium, electron density and electron temperature can be recovered jointly, and finally, that this can be done in the presence of realistic imaging systems that incorporate sensor cropping and quantisation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Spatial and temporal prediction of radiation dose rates near Fukushima Daiichi Nuclear Power Plant

In this paper, we have developed a methodology to estimate the spatiotemporal distribution of radiation air dose rates around the Fukushima Daiichi Nuclear Power Plant (FDNPP). In our exploratory data analysis, we found that (1) the temporal evolution of dose rates is composed of a log-linear decay trend and fluctuations of air dose rates that are spatially correlated among adjacent monitoring posts; and (2) the slope of the log-linear environmental decay trend can be represented as a function of the apparent initial dose rates, coordinate position, land-use type, and soil type. From these observations, we first estimated the log-linear decay trend at each location based on these predictors, using the random forest method. We then developed a modified Kalman filter coupled with a Gaussian process model to estimate the dose-rate time series at a given location and time. We applied this method to the Fukushima evacuation zone (as of March 2017), which included 17 monitoring post locations (with monitoring datasets collected between 2014 and 2018) and generated a time series of dose-rate maps. Our results show that this approach allows us to produce accurate spatial and temporal predictions of radiation dose-rate maps using limited spatiotemporal measurements.

61 RADIATION PROTECTION AND DOSIMETRY↗

Increasing the power of weak lensing data with multipole-based intrinsic alignment estimators

ABSTRACT It has long been known that galaxy shapes align coherently with the large-scale density field. Characterizing this effect is essential to interpreting measurements of weak gravitational lensing, the deflection of light from distant galaxies by matter overdensities along the line of sight, as it also produces coherent galaxy alignments that we wish to interpret in terms of a cosmological model. Existing direct measurements of intrinsic alignments using galaxy samples with high-quality shape and redshift measurements typically use well-understood but sub-optimal projected estimators, which do not make good use of the information in the data when comparing those estimators to theoretical models. We demonstrate a more optimal estimator, based on a multipole expansion of the correlation functions or power spectra, for direct measurements of galaxy intrinsic alignments. We show that even using the lowest order multipole alone increases the significance of inferred model parameters using simulated and real data, without any additional modelling complexity. We apply this estimator to measurements of parameters of the non-linear alignment model using data from the Sloan Digital Sky survey, demonstrating consistent results with a factor of ∼2 greater precision in parameter fits to intrinsic alignments models. This result is functionally equivalent to quadrupling the survey area, but without the attendant costs – thereby demonstrating the value in using this new estimator in current and future intrinsic alignments measurements using spectroscopic galaxy samples.

Astronomy & Astrophysics↗

Coupled Groups of g-Modes in a Sun with Mixed Core

Groups of linear g-modes can sum to create long-lived nonlinear oscillations in small "hot volumes" very deep in the Sun that help drive the modes. In these volumes (dimensions -10 Mm), the time average rate of He-3 burning doubles as temperature fluctuations exceed 10% and rises by an order of magnitude for fluctuations of 25%. To be consistent with locally large motions, we impose a mixed shell on an otherwise standard solar model before computing g-mode solutions. Mixing in the assumed shell r = (0.10+/-0.03) R(sub sun) is rapid (<<10(exp 6) yr) with slower mixing somewhat beyond. If l is the principal spherical harmonic index, a set of g-modes for any single l less than or equal to 15 with five consecutive radial harmonics can be excited with nearly linear thermal amplitudes, A(sub T) less than or equal to 0.053, throughout the star and a fractional temperature fluctuation in its hot volume of (Delta)T/T less than or equal to 0.18. These thresholds for excitation will become smaller when sets for several values of l are computed simultaneously. There is some evidence for the rotation of g-mode sets in the long solar activity record and g-mode upward wave flux has been suggested to explain the 1.3 yr reversing flows tentatively detected below the Sun's convective envelope (CE). The large local amplitudes needed for excitation implies that g-modes may transport a non-negligible fraction of the solar luminosity, yet their near linear amplitudes outside the hot volume suggests amplitudes over most of the solar surface that would be barely detectable for l > 3. A formalism is presented for summing the g-modes and estimating growth rates under the approximation that modes are strictly linear except in a hot volume which holds only a few percent of mode kinetic energy. Finally over the range 2 less than or equal to l less than or equal to 30, we summed all zonal harmonics, m, for a given l and computed the relative angular orientations that would release the most nuclear energy. This should be close to the physically preferred angular state of such a family and a few examples were displayed.

Wolff, Charles L.↗