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,657 records · Page 92
All sourcesNASA NTRSDOE OSTI
NASA NTRS1986

Spline-based distributed system identification with application to large space antennas

A parameter and state estimation technique for distributed models is demonstrated through the solution of a problem generic to large space antenna system identification. Assuming the position of the reflective surface of the maypole (hoop/column) antenna to be approximated by the static two-dimensional, stretched-membrane partial differential equation with variable-stiffness coefficient functions, a spline-based approximation procedure is described that estimates the shape and stiffness functions from data set observations. For given stiffness functions, the Galerkin projection with linear spline-based functions is applied to project the distributed problem onto a finite-dimensional subspace wherein algebraic equations exist for determining a static shape (state) prediction. The stiffness functions are then parameterized by cubic splines and the parameters estimated by an output error technique. Numerical results are presented for data descriptive of a 100-m-diameter maypole antenna.

Banks, H. T.↗
DOE OSTI2013

Turbulent regimes in the tokamak scrape-off layer

The non-linear turbulent regimes in the tokamak scrape-off layer (SOL) are identified according to the linear instability responsible for the perpendicular transport. Four regions of the SOL operational parameters are determined where turbulence is driven by the inertial or resistive branches of the ballooning mode or of drift waves. The analysis, based on the linear electrostatic drift-reduced Braginskii equations, evaluates the pressure scale length self-consistently from the balance between plasma losses at the vessel and perpendicular turbulent transport. The latter is estimated by assuming that turbulence saturation occurs due to a local flattening of the plasma gradients and associated removal of the linear instability drive; it is also shown that transport is led by the mode that maximizes the ratio of the linear growth to the poloidal wavenumber. The methodology used to identify the turbulent regimes is confirmed by the results of non-linear simulations of SOL turbulence. The identification of the turbulent regimes, the predicted pressure scale length, and the poloidal wavenumber of the leading mode are in reasonable agreement with non-linear simulation results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Monte Carlo Perturbation Theory Estimates of Sensitivities to System Dimensions

Here, Monte Carlo methods are developed using adjoint-based perturbation theory and the differential operator method to compute the sensitivities of the k-eigenvalue, linear functions of the flux (reaction rates), and bilinear functions of the forward and adjoint flux (kinetics parameters) to system dimensions for uniform expansions or contractions. The calculation of sensitivities to system dimensions requires computing scattering and fission sources at material interfaces using collisions occurring at the interface—which is a set of events with infinitesimal probability. Kernel density estimators are used to estimate the source at interfaces using collisions occurring near the interface. The methods for computing sensitivities of linear and bilinear ratios are derived using the differential operator method and adjoint-based perturbation theory and are shown to be equivalent to methods previously developed using a collision history–based approach. The methods for determining sensitivities to system dimensions are tested on a series of fast, intermediate, and thermal critical benchmarks as well as a pressurized water reactor benchmark problem with iterated fission probability used for adjoint-weighting. The estimators are shown to agree within 5% and 3σ of reference solutions obtained using direct perturbations with central differences for the majority of test problems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1982

Robust Kalman filter design for active flutter suppression systems

Additional insight is provided into the use of the Doyle-Stein (1979, 1981) technique in aeroelastic control problems by examining the application of the method to a flutter control problem. The system to be controlled consists of a full-size wind tunnel model of a wing, plus an aileron, an actuator, and an accelerometer used to sense the motion of the wing. A full-state feedback controller was designed using linear optimal control theory, and a Kalman filter was used in the feedback loop for state estimation. The filter design procedure is explained along with that to improve closed-loop properties of the system. The locus of the poles of the filter is examined as a scalar design parameter is varied. The Doyle-Stein design procedure is shown to substantially improve the stability properties of an active flutter controller designed using the linear quadratic Gaussian control theory.

Garrard, W. L.↗
DOE OSTI2022

Multimodal parameter spaces of a complex multi-channel neuron model

One of the most common types of models that helps us to understand neuron behavior is based on the Hodgkin–Huxley ion channel formulation (HH model). A major challenge with inferring parameters in HH models is non-uniqueness: many different sets of ion channel parameter values produce similar outputs for the same input stimulus. Such phenomena result in an objective function that exhibits multiple modes (i.e., multiple local minima). This non-uniqueness of local optimality poses challenges for parameter estimation with many algorithmic optimization techniques. HH models additionally have severe non-linearities resulting in further challenges for inferring parameters in an algorithmic fashion. To address these challenges with a tractable method in high-dimensional parameter spaces, we propose using a particular Markov chain Monte Carlo (MCMC) algorithm, which has the advantage of inferring parameters in a Bayesian framework. The Bayesian approach is designed to be suitable for multimodal solutions to inverse problems. We introduce and demonstrate the method using a three-channel HH model. We then focus on the inference of nine parameters in an eight-channel HH model, which we analyze in detail. We explore how the MCMC algorithm can uncover complex relationships between inferred parameters using five injected current levels. The MCMC method provides as a result a nine-dimensional posterior distribution, which we analyze visually with solution maps or landscapes of the possible parameter sets. The visualized solution maps show new complex structures of the multimodal posteriors, and they allow for selection of locally and globally optimal value sets, and they visually expose parameter sensitivities and regions of higher model robustness. We envision these solution maps as enabling experimentalists to improve the design of future experiments, increase scientific productivity and improve on model structure and ideation when the MCMC algorithm is applied to experimental data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2013

Linearity of Climate Response to Increases in Black Carbon Aerosols

The impact of absorbing aerosols on global climate are not completely understood. Here, we present results of idealized experiments conducted with the Community Atmosphere Model (CAM4) coupled to a slab ocean model (CAM4-SOM) to simulate the climate response to increases in tropospheric black carbon aerosols (BC) by direct and semi-direct effects. CAM4-SOM was forced with 0, 1x, 2x, 5x and 10x an estimate of the present day concentration of BC while maintaining their estimated present day global spatial and vertical distribution. The top of the atmosphere (TOA) radiative forcing of BC in these experiments is positive (warming) and increases linearly as the BC burden increases. The total semi-direct effect for the 1x experiment is positive but becomes increasingly negative for higher BC concentrations. The global average surface temperature response is found to be a linear function of the TOA radiative forcing. The climate sensitivity to BC from these experiments is estimated to be 0.42 K $ W^{-1} m^{2}$ when the semi-direct effects are accounted for and 0.22 K $ W^{-1} m^{2}$ with only the direct effects considered. Global average precipitation decreases linearly as BC increases, with a precipitation sensitivity to atmospheric absorption of 0.4 $\%$ $W^{-1}m^{2}$ . The hemispheric asymmetry of BC also causes an increase in southward cross-equatorial heat transport and a resulting northward shift of the inter-tropical convergence zone in the simulations at a rate of 4$^{\circ}$N $ PW^{-1}$. Global average mid- and high-level clouds decrease, whereas the low-level clouds increase linearly with BC. The increase in marine stratocumulus cloud fraction over the south tropical Atlantic is caused by increased BC-induced diabatic heating of the free troposphere.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Detection of the pairwise kinematic Sunyaev-Zel’dovich effect and pairwise velocity with DESI DR1 galaxies and ACT DR6 and Planck CMB data

We present a 9.3⁢𝜎 detection of the pairwise kinematic Sunyaev-Zel’dovich (kSZ) effect by combining a sample of 913,286 Luminous Red Galaxies (LRGs) from the Dark Energy Spectroscopic Instrument Data Release 1 (DESI DR1) catalog and coadded Atacama Cosmology Telescope (ACT DR6) and Planck cosmic microwave background (CMB) temperature maps. This represents the highest-significance pairwise kSZ measurement to date. The analysis uses three ACT CMB temperature maps: coadded 150 GHz, total frequency maps, and a component separated Internal Linear Combination (ILC) map, all of which cover 19,000 square degrees of the sky from Advanced ACTPol observations conducted between 2017 and 2022. Comparison of the results of these three maps serves as a consistency check for potential foreground contamination that may depend on the observation frequency. An estimate of the best-fit mass-averaged optical depth is obtained by comparing the pairwise kSZ curve with the linear theory prediction of the pairwise velocity under the best-fit Planck cosmology and is compared with predictions from simulations. This estimate serves as a reference point for future comparisons with thermal SZ–derived optical depth measurements for the same DESI cluster samples, which will be presented in a companion paper. Finally, we employ a machine learning approach, trained on simulations to estimate the optical depth for 456,803 DESI LRG-identified clusters within the simulated mass range ( ≳ 10 13 ⁢𝑀 ⊙ ). These are combined with the measured kSZ signal to infer the individual cluster peculiar velocities, providing the opportunity to constrain the behavior of gravity and the dark sector over a range of cosmic scales and epochs.

Gong, Yulin [Cornell University] (ORCID:0000000346↗
NASA NTRS1981

A model describing the microwave emission from a multi-layer snowpack at 37 GHz

A multilayer emission model is described and applied to emission measurements obtained at 37 GHz and H polarization using a microwave radiometer attached to a truck-mounted boom in Steamboat Springs, Colorado in 1977. Estimated absorption and scattering coefficients and their dependence on wetness were obtained using calculated values of the dielectric constant at 37 GHz along with the model. It was found that the scattering coefficient is comparable in value to the absorption coefficient for dry snow however, the absorption coefficient increases linearly with increasing snow wetness while the scattering coefficient decreases linearly with increasing wetness. The emission from each layer of the snowpack was also calculated using the estimated coefficients. It is shown that for dry snow, the ground underneath the snowpack contributes about 45% of all measured emission while the rest is due to emission from all the layers within the snowpack. When the wetness of the top 5 cm layer of snowpack increases to about 2% by volume, this top 5 cm snowlayer contributes more than 90% of all the measured emission.

Abdelrazik, M.↗
NASA NTRS2005

Elastic Thickness Estimates for Coronae Associated with Chasmata on Venus

Coronae are large-scale circular tectonic features surrounded by annular ridges. They are generally considered unique to Venus and may offer insights into the differences in lithospheric structure or mantle convective pattern between Venus and Earth. 68% of all coronae are associated with chasmata or fracture belts. The remaining 32% are located at volcanic rises or in the plains. Chasmata are linear to arcuate troughs, with trough parallel fractures and faults which extend for 1000 s of kilometers. Estimates of the elastic thickness of the lithosphere (T(sub e)) have been calculated in a number of gravity/topography studies of Venus and for coronae specifically. None of these studies, however, have explored the dependence of T(sub e) on the tectonic history of the region, as implied from the interpretation of relative timing relationships between coronae and surrounding features. We examine the relationship between the local T(sub e) and the relative ages of coronae and chasmata with the aim of further constraining the origin and evolution of coronae and chasmata systems.

Hoogenboom, T.↗
NASA NTRS2016

Uncertainty Estimates of Psychoacoustic Thresholds Obtained from Group Tests

Adaptive psychoacoustic test methods, in which the next signal level depends on the response to the previous signal, are the most efficient for determining psychoacoustic thresholds of individual subjects. In many tests conducted in the NASA psychoacoustic labs, the goal is to determine thresholds representative of the general population. To do this economically, non-adaptive testing methods are used in which three or four subjects are tested at the same time with predetermined signal levels. This approach requires us to identify techniques for assessing the uncertainty in resulting group-average psychoacoustic thresholds. In this presentation we examine the Delta Method of frequentist statistics, the Generalized Linear Model (GLM), the Nonparametric Bootstrap, a frequentist method, and Markov Chain Monte Carlo Posterior Estimation and a Bayesian approach. Each technique is exercised on a manufactured, theoretical dataset and then on datasets from two psychoacoustics facilities at NASA. The Delta Method is the simplest to implement and accurate for the cases studied. The GLM is found to be the least robust, and the Bootstrap takes the longest to calculate. The Bayesian Posterior Estimate is the most versatile technique examined because it allows the inclusion of prior information.

Rathsam, Jonathan↗
DOE OSTI2019

Mathematical Modelling of Gas–Liquid, Two-Phase Flows in a Ladle Shroud

Differential and macroscopic models of argon–steel flows in ladle shroud have been developed. In this, argon–steel, two-phase flow phenomena have been formulated via a transient, three dimensional, turbulent flow model, based on the volume of fluid (VOF) calculation procedure. While realizable k–ε turbulence model has been applied to map turbulence, commercial, CFD software ANSYS-Fluent™ (Version 18), has been applied to carry out numerical calculations. Predictions from the model have been directly assessed against experimental measurements across the range of shroud dimensions and volumetric flow rates typically practiced in the industry. It is demonstrated that the two-phase turbulent flow model captures the general features of gas–liquid flows in ladle shroud providing estimates of free jet length and threshold gas flow rates (required to halt air ingression) which are in agreement with corresponding experimental measurements. In the absence of differential solutions, a macroscopic model has been worked out through dimensional analysis embodying multiple non-linear regression. It is shown that dimensionless free jet length in bloom and slab casting shrouds can be estimated reasonably accurately from the following correlation (in SI unit), viz., ((L{sub jet})/(D{sub CN}))=2.8×10{sup −2}(((Q{sub G})/(Q{sub L}))){sup 1.14}(((gD{sub CN}{sup 5})/(Q{sub L}{sup 2}))){sup 0.8}(((σD{sub CN}{sup 3})/(ρ{sub L}Q{sub L}{sup 2}))){sup −0.9}(((D{sub sh})/(D{sub CN}))){sup 2.0}(((ρ{sub G})/(ρ{sub L}))){sup −0.30}in which, L{sub jet} is the free liquid jet length (m), Q{sub G} is the gas flow rate (m{sup 3}/s), Q{sub L} is the liquid flow rate (m{sup 3}/s), D{sub sh} is shroud diameter (m), D{sub CN} is the collector nozzle diameter (m), σ is the interfacial tension (N/m), and ρ{sub G} as well as ρ{sub L} are respectively density of gas and liquid (kg/m{sup 3}). It is demonstrated that the proposed correlation is consistent with the laws of physical modeling and leads to estimates that are in good agreement with predictions from the differential models, for both air-water as well as argon–steel systems. Numerical simulations as well as macroscopic modeling have indicated that thermo-physical properties of the gas–liquid system are important and exert some influences on the gas–liquid, two-phase, flow in ladle shrouds, albeit not to a large extent. Despite dissimilar thermo-physical properties, full scale water modeling appears to be sufficiently predictive and provides reasonable macroscopic descriptions of the two-phase flow phenomena in industrial ladle shroud systems.

42 ENGINEERING↗
DOE OSTI2020

Projecting Antarctica's contribution to future sea level rise from basal ice shelf melt using linear response functions of 16 ice sheet models (LARMIP-2)

Abstract. The sea level contribution of the Antarctic ice sheet constitutes a large uncertainty in future sea level projections. Here we apply a linear response theory approach to 16 state-of-the-art ice sheet models to estimate the Antarctic ice sheet contribution from basal ice shelf melting within the 21st century. The purpose of this computation is to estimate the uncertainty of Antarctica's future contribution to global sea level rise that arises from large uncertainty in the oceanic forcing and the associated ice shelf melting. Ice shelf melting is considered to be a major if not the largest perturbation of the ice sheet's flow into the ocean. However, by computing only the sea level contribution in response to ice shelf melting, our study is neglecting a number of processes such as surface-mass-balance-related contributions. In assuming linear response theory, we are able to capture complex temporal responses of the ice sheets, but we neglect any self-dampening or self-amplifying processes. This is particularly relevant in situations in which an instability is dominating the ice loss. The results obtained here are thus relevant, in particular wherever the ice loss is dominated by the forcing as opposed to an internal instability, for example in strong ocean warming scenarios. In order to allow for comparison the methodology was chosen to be exactly the same as in an earlier study (Levermann et al., 2014) but with 16 instead of 5 ice sheet models. We include uncertainty in the atmospheric warming response to carbon emissions (full range of CMIP5 climate model sensitivities), uncertainty in the oceanic transport to the Southern Ocean (obtained from the time-delayed and scaled oceanic subsurface warming in CMIP5 models in relation to the global mean surface warming), and the observed range of responses of basal ice shelf melting to oceanic warming outside the ice shelf cavity. This uncertainty in basal ice shelf melting is then convoluted with the linear response functions of each of the 16 ice sheet models to obtain the ice flow response to the individual global warming path. The model median for the observational period from 1992 to 2017 of the ice loss due to basal ice shelf melting is 10.2 mm, with a likely range between 5.2 and 21.3 mm. For the same period the Antarctic ice sheet lost mass equivalent to 7.4 mm of global sea level rise, with a standard deviation of 3.7 mm (Shepherd et al., 2018) including all processes, especially surface-mass-balance changes. For the unabated warming path, Representative Concentration Pathway 8.5 (RCP8.5), we obtain a median contribution of the Antarctic ice sheet to global mean sea level rise from basal ice shelf melting within the 21st century of 17 cm, with a likely range (66th percentile around the mean) between 9 and 36 cm and a very likely range (90th percentile around the mean) between 6 and 58 cm. For the RCP2.6 warming path, which will keep the global mean temperature below 2 ∘C of global warming and is thus consistent with the Paris Climate Agreement, the procedure yields a median of 13 cm of global mean sea level contribution. The likely range for the RCP2.6 scenario is between 7 and 24 cm, and the very likely range is between 4 and 37 cm. The structural uncertainties in the method do not allow for an interpretation of any higher uncertainty percentiles. We provide projections for the five Antarctic regions and for each model and each scenario separately. The rate of sea level contribution is highest under the RCP8.5 scenario. The maximum within the 21st century of the median value is 4 cm per decade, with a likely range between 2 and 9 cm per decade and a very likely range between 1 and 14 cm per decade.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

The value of continuity: Refined isogeometric analysis and fast direct solvers

Here, we propose the use of highly continuous finite element spaces interconnected with low continuity hyperplanes to maximize the performance of direct solvers. Starting from a highly continuous Isogeometric Analysis (IGA) discretization, we introduce C0-separators to reduce the interconnection between degrees of freedom in the mesh. By doing so, both the solution time and best approximation errors are simultaneously improved. We call the resulting method “refined Isogeometric Analysis (rIGA)”. To illustrate the impact of the continuity reduction, we analyze the number of Floating Point Operations (FLOPs), computational times, and memory required to solve the linear system obtained by discretizing the Laplace problem with structured meshes and uniform polynomial orders. Theoretical estimates demonstrate that an optimal continuity reduction may decrease the total computational time by a factor between p 2 and p 3 , with pp being the polynomial order of the discretization. Numerical results indicate that our proposed refined isogeometric analysis delivers a speed-up factor proportional to p 2 . In a 2D mesh with four million elements and p=5, the linear system resulting from rIGA is solved 22 times faster than the one from highly continuous IGA. In a 3D mesh with one million elements and p=3, the linear system is solved 15 times faster for the refined than the maximum continuity isogeometric analysis.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

The increase in animal mortality risk following exposure to sparsely ionizing radiation is not linear quadratic with dose

The US government regulates allowable radiation exposures relying, in large part, on the seventh report from the committee to estimate the Biological Effect of Ionizing Radiation (BEIR VII), which estimated that most contemporary exposures- protracted or low-dose, carry 1.5 fold less risk of carcinogenesis and mortality per Gy than acute exposures of atomic bomb survivors. This correction is known as the dose and dose rate effectiveness factor for the life span study of atomic bomb survivors (DDREF LSS ). As a result, it was calculated by applying a linear-quadratic dose response model to data from Japanese atomic bomb survivors and a limited number of animal studies.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗
DOE OSTI2016

Enhanced Performance Controller Design for Stochastic Systems by Adding Extra State Estimation onto the Existing Closed Loop Control

To enhance the performance of the tracking property , this paper presents a novel control algorithm for a class of linear dynamic stochastic systems with unmeasurable states, where the performance enhancement loop is established based on Kalman filter. Without changing the existing closed loop with the PI controller, the compensative controller is designed to minimize the variances of the tracking errors using the estimated states and the propagation of state variances. Moreover, the stability of the closed-loop systems has been analyzed in the mean-square sense. A simulated example is included to show the effectiveness of the presented control algorithm, where encouraging results have been obtained.

Stochastic coupling attenuation↗
NASA NTRS1979

Maximum likelihood estimation for life distributions with competing failure modes

The general model for the competing failure modes assuming that location parameters for each mode are expressible as linear functions of the stress variables and the failure modes act independently is presented. The general form of the likelihood function and the likelihood equations are derived for the extreme value distributions, and solving these equations using nonlinear least squares techniques provides an estimate of the asymptotic covariance matrix of the estimators. Monte-Carlo results indicate that, under appropriate conditions, the location parameters are nearly unbiased, the scale parameter is slightly biased, and the asymptotic covariances are rapidly approached.

Sidik, S. M.↗
NASA NTRS1979

A modification to linearized theory for prediction of pressure loadings on lifting surfaces at high supersonic Mach numbers and large angles of attack

A new linearized-theory pressure-coefficient formulation was studied. The new formulation is intended to provide more accurate estimates of detailed pressure loadings for improved stability analysis and for analysis of critical structural design conditions. The approach is based on the use of oblique-shock and Prandtl-Meyer expansion relationships for accurate representation of the variation of pressures with surface slopes in two-dimensional flow and linearized-theory perturbation velocities for evaluation of local three-dimensional aerodynamic interference effects. The applicability and limitations of the modification to linearized theory are illustrated through comparisons with experimental pressure distributions for delta wings covering a Mach number range from 1.45 to 4.60 and angles of attack from 0 to 25 degrees.

Carlson, H. W.↗