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,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

In-flight polarization angle calibration for LiteBIRD: blind challenge and cosmological implications

We present a demonstration of the in-flight polarization angle calibration for the JAXA/ISAS second strategic large class mission, LiteBIRD, and estimate its impact on the measurement of the tensor-to-scalar ratio parameter, r, using simulated data. We generate a set of simulated sky maps with CMB and polarized foreground emission, and inject instrumental noise and polarization angle offsets to the 22 (partially overlapping) LiteBIRD frequency channels. Our in-flight angle calibration relies on nulling the EB cross correlation of the polarized signal in each channel. This calibration step has been carried out by two independent groups with a blind analysis, allowing an accuracy of the order of a few arc-minutes to be reached on the estimate of the angle offsets. Both the corrected and uncorrected multi-frequency maps are propagated through the foreground cleaning step, with the goal of computing clean CMB maps. We employ two component separation algorithms, the Bayesian-Separation of Components and Residuals Estimate Tool (B-SeCRET), and the Needlet Internal Linear Combination (NILC). We find that the recovered CMB maps obtained with algorithms that do not make any assumptions about the foreground properties, such as NILC, are only mildly affected by the angle miscalibration. However, polarization angle offsets strongly bias results obtained with the parametric fitting method. Once the miscalibration angles are corrected by EB nulling prior to the component separation, both component separation algorithms result in an unbiased estimation of the r parameter. While this work is motivated by the conceptual design study for LiteBIRD, its framework can be broadly applied to any CMB polarization experiment. In particular, the combination of simulation plus blind analysis provides a robust forecast by taking into account not only detector sensitivity but also systematic effects.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1984

The temporal behavior of upper stratospheric ozone at low latitudes - Evidence from Nimbus 4 BUV data for short-term responses to solar ultraviolet variability

Ozone mixing ratios at pressure levels near 2 mbar are analyzed for the purpose of estimating the average response of the upper stratospheric ozone to solar UV variability on the time scale of the solar rotation period. The data were obtained from observations made with the Nimbus 4 backscattering UV radiometer in 13 latitude zones between 65 degrees N and 65 degrees S. The anaysis showed that temporal variations are negatively correlated with changes in zonally averaged equivalent temperature measured simultaneously by Nimbus 4 radiometer. A linear regression analysis is performed in order to obtain estimates of the average percent change of ozone at low latitudes, on the considered time scale for given changes in 10.7-cm flux, and in the UV flux model developed by Lean et al. (1982). Reproductions of the ozone profiles are provided.

Hood, L. L.↗
DOE OSTI2020

On the inference of viscoelastic constants from stress relaxation experiments

Several constitutive theories have been proposed in the literature to model the viscoelastic response of materials, including widely used rheological constitutive models. These models are characterized by certain parameters ("time constants") that define the time scales over which the material relaxes. These parameters are primarily obtained from stress relaxation experiments using curve-fitting techniques. However, the question of how best to estimate these time constants remains open. As a step towards answering this question, we propose an optimal experimental design approach based on ideas from information geometry, namely Fisher information and Kullback-Leibler divergence. The material is modeled as a spring element in parallel with multiple Maxwell elements and described using a one- or two-term Prony series. Treating the time constants as unknowns, we develop expressions for the Fisher information and Kullback-Leibler divergence that allow us to maximize information gain from experimental data. Based on the results of this study, we propose that the largest time constant estimated from a stress relaxation experiment for a linear viscoelastic material should be at most one-fifth of the total time of the experiment in order to maximize information gain. Our results also provide confirmation that the equilibrium modulus of the material cannot be reliably determined from curve-fitting to data from a stress relaxation experiment.

Fisher information↗
DOE OSTI2017

Optimal design and dispatch of a system of diesel generators, photovoltaics and batteries for remote locations

Renewable energy technologies, specifically, solar photovoltaic cells, combined with battery storage and diesel generators, form a hybrid system capable of independently powering remote locations, i.e., those isolated from larger grids. If sized correctly, hybrid systems reduce fuel consumption compared to diesel generator-only alternatives. We present an optimization model for establishing a hybrid power design and dispatch strategy for remote locations, such as a military forward operating base, that models the acquisition of different power technologies as integer variables and their operation using nonlinear expressions. Our cost-minimizing, nonconvex, mixed-integer, nonlinear program contains a detailed battery model. Due to its complexities, we present linearizations, which include exact and convex under-estimation techniques, and a heuristic, which determines an initial feasible solution to serve as a “warm start” for the solver. We determine, in a few hours at most, solutions within 5% of optimality for a candidate set of technologies; these solutions closely resemble those from the nonlinear model. Lastly, our instances contain real data spanning a yearly horizon at hour fidelity and demonstrate that a hybrid system could reduce fuel consumption by as much as 50% compared to a generator-only solution.

14 SOLAR ENERGY↗
DOE OSTI2018

Numerical integration method for triple integrals using the second kind Chebyshev wavelets and Gauss–Legendre quadrature

In this paper, a computational method combining the second kind Chebyshev wavelets and Gauss–Legendre quadrature is proposed for numerical integrations of arbitrary functions over regions like cuboid, tetrahedron, cylinder, cone, paraboloid and ellipsoid. Gauss–Legendre quadrature is used to convert a triple integral into a double integral and integral regions are transformed to the standard integration region by linear and nonlinear transformation. Moreover, convergence and accuracy estimation of the second kind Chebyshev wavelets expansion of two dimensions is given. Illustrative examples have been demonstrated to show the applicability and accuracy of the present method.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Trunk muscle forces and spinal loads while walking in persons with lower limb amputation: Influences of chronic low back pain

Persons with lower limb amputation (LLA) are at high risk for developing chronic low back pain (LBP), often with biomechanical factors considered as likely contributors. Here, trunk and pelvis kinematics, muscle forces, and resultant spinal loads were characterized in persons with LLA, with and without chronic LBP. Thirty-five persons with unilateral LLA – 19 with chronic LBP (“LLA-cLBP”), 16 without LBP (“LLA-nLBP”) – and 15 (uninjured) persons without LBP (“CTR-nLBP”) walked overground (1.3 m/s) while thorax and pelvis kinematics were tracked (and ranges of motion [ROM] computed), and used as inputs for a non-linear finite element model of the spine to estimate global and local muscle forces, and resultant spinal loads. In the frontal and transverse planes, thorax ROM were up to 66.6% smaller in LLA-nLBP versus LLA-cLBP (P < 0.001) and CTR-nLBP (P < 0.001). In the sagittal plane, pelvis ROM was 50.4% smaller in LLA-nLBP versus LLA-cLBP (P = 0.014). LLA-cLBP exhibited 45.5% and 34.2% greater peak local and global muscle forces, respectively, versus CTR-nLBP (P < 0.011). Up to 48.1% greater spinal loads were observed in LLA-cLBP versus CTR-nLBP (P < 0.013); peak compression and local muscle forces were respectively 20.2% and 41.0% larger in LLA-nLBP versus CTR-nLBP (P < 0.005). Despite differences in trunk and pelvis kinematics between LLA-cLBP and LLA-nLBP, trunk muscle forces and spinal loads were similar (P > 0.101) between these groups. Similar loading parameters regardless of LBP presence, while highly dependent on trunk muscle activation strategies, may mitigate further accumulation of mechanical fatigue. Finally, it remains important to understand the temporality of loading with respect to LBP onset following LLA.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2017

Numerical Analysis of Experiments on Thermally Induced Dissociation of Methane Hydrates in Porous Media

Numerical simulation is essential for the prediction and evaluation of hydrocarbon reservoir performance. Numerical simulators developed for the description of the behavior of hydrates under production and the corresponding flow of fluids and heat accounting for all known processes are powerful, but they need validation through comparison to field or experimental data in order to instill confidence in their predictions. In this study, we analyze by means of numerical simulation the results of an experiment of methane hydrate dissociation by thermal stimulation in unconsolidated porous media heated through the vessel walls. The physics captured by the model include multicomponent heat and mass transfer, multiphase flow through porous media, and the phase behavior of the CH 4 + H 2 O system involved in methane hydrate formation and dissociation. The set of governing equations consists of the mass and energy conservation equations coupled with constitutive relationships, i.e., the dissolution of gas in H 2 O, relative permeability and capillary pressure models, composite thermal conductivity models, and methane hydrate phase equilibria. The model geometry describes accurately the hydrate reactor used in a recent experimental study investigating methane hydrate dissociation behavior [Chong et. al. Appl. Energy 2016, 177, 409-421]. The cumulative gas production is estimated and validated against three tests of experimental data involving different boundary temperatures, showing a good agreement between observations and numerical predictions. The predicted evolution of the spatial distributions of different phases over time shows that hydrate dissociation progresses inward from the reactor boundary to the center, methane gas accumulates to the top of the reactor because of buoyancy, and water migrates down to the bottom of the reactor because of gravity. A sharp hydrate dissociation front is predicted, and the estimated location of hydrate dissociation front suggests a linear relationship with the square root of time. A sensitivity analysis on the thermal conductivity of sand under fully saturated conditions is conducted to elucidate its effect on the gas production behavior. In addition, the energy efficiency ratio computed from the simulation of this boundary-wall heating technique varies from 14.0 to 16.2. Deviations between observations and predictions of the evolution of the temperature profile are attributed to initial heterogeneous distribution of the hydrate phase in the hydrate reactor.

42 ENGINEERING↗
DOE OSTI2016

Impact of post-Born lensing on the CMB

Lensing of the CMB is affected by post-Born lensing, producing corrections to the convergence power spectrum and introducing field rotation. We show numerically that the lensing convergence power spectrum is affected at the ∼< 0.2% level on accessible scales, and that this correction and the field rotation are negligible for observations with arcminute beam and noise levels ∼> 1 μK arcmin. The field rotation generates ∼ 2.5% of the total lensing B-mode polarization amplitude (0.2% in power on small scales), but has a blue spectrum on large scales, making it highly subdominant to the convergence B modes on scales where they are a source of confusion for the signal from primordial gravitational waves. Since the post-Born signal is non-linear, it also generates a bispectrum with the convergence. We show that the post-Born contributions to the bispectrum substantially change the shape predicted from large-scale structure non-linearities alone, and hence must be included to estimate the expected total signal and impact of bispectrum biases on CMB lensing reconstruction quadratic estimators and other observables. The field-rotation power spectrum only becomes potentially detectable for noise levels || 1 μK arcmin, but its bispectrum with the convergence may be observable at ∼ 3σ with Stage IV observations. Rotation-induced and convergence-induced B modes are slightly correlated by the bispectrum, and the bispectrum also produces additional contributions to the lensed BB power spectrum.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Study of ion-gyroscale fluctuations in low-density L-mode plasmas heated by NBI on KSTAR

Broadband density fluctuations with peak frequency ranging from 150 to 400 kHz were measured using a multichannel microwave imaging reflectometer in core region of the low-density L-mode plasmas heated by neutral beam injection on KSTAR. These fluctuations have been studied by comparing the dominant mode scales estimated from the measurement with those predicted from linear gyrokinetic simulation. The measured poloidal wavenumbers are qualitatively comparable to those of the 'fastest growing modes' from simulations, whereas they are larger than those of the 'transport-dominant modes' by about a factor of three. The agreement on wavenumbers between the measurement and linear simulation (for the fastest growing modes) is probably due to sufficiently weak $E \times B$ flow shear compared to the maximum linear growth rate. Meanwhile, the transport-dominant modes seem to be related to the fluctuations in lower frequencies (~80–150 kHz) observed in some of the measurement.

Physics↗
DOE OSTI2021

Influence of inhomogeneous stochasticity on the falsifiability of mean-field theories and examples from accretion disc modelling

ABSTRACT Despite spatial and temporal fluctuations in turbulent astrophysical systems, mean-field theories can be used to describe their secular evolution. However, observations taken over time scales much shorter than dynamical time scales capture a system in a single state of its turbulence ensemble. Comparing with mean-field theory can falsify the latter only if the theory is additionally supplied with a quantified precision. The central limit theorem provides appropriate estimates to the precision only when fluctuations contribute linearly to an observable and with constant coherent scales. Here, we introduce an error propagation formula that relaxes both limitations, allowing for non-linear functional forms of observables and inhomogeneous coherent scales and amplitudes of fluctuations. The method is exemplified in the context of accretion disc theories, where inhomogeneous fluctuations in the surface temperature are propagated to the disc emission spectrum – the latter being a non-linear and non-local function of the former. The derived precision depends non-monotonically on emission frequency. Using the same method, we investigate how binned spectral fluctuations in telescope data change with the spectral resolving power. We discuss the broader implications for falsifiability of a mean-field theory.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Nonlinear Photovoltaic Degradation Rates: Modeling and Comparison Against Conventional Methods

Although common practice for estimating photovoltaic (PV) degradation rate (RD) assumes a linear behavior, field data have shown that degradation rates are frequently nonlinear. This article presents a new methodology to detect and calculate nonlinear RD based on PV performance time-series from nine different systems over an eight-year period. Prior to performing the analysis and in order to adjust model parameters to reflect actual PV operation, synthetic datasets were utilized for calibration purposes. A change-point analysis is then applied to detect changes in the slopes of PV trends, which are extracted from constructed performance ratio (PR) time-series. Once the number and location of change points is found, the ordinary least squares method is applied to the different segments to compute the corresponding rates. The obtained results verified that the extracted trends from the PR time-series may not always be linear and therefore, “nonconventional” models need to be applied. All thin-film technologies demonstrated nonlinear behavior whereas nonlinearity detected in the crystalline silicon systems is thought to be due to a maintenance event. A comparative analysis between the new methodology and other conventional methods demonstrated levelized cost of energy differences of up to 6.14%, highlighting the importance of considering nonlinear degradation behavior

14 SOLAR ENERGY↗
DOE OSTI2024

CONSTRAINT-INDEPENDENT CONSTANT CTOA DETERMINATION FOR DUCTILE STABLE CRACK GROWTH

Crack tip opening angle (CTOA) has been used as a reliable fracture toughness parameter for decades to characterize stable ductile crack growth for thin-walled aerospace structures in the low-constraint conditions. Recently, the CTOA parameter was also applied to the pipeline industry, and a CTOA test standard ASTM E3039 was thus developed for testing a critical constant CTOA. Research showed that the constant CTOA can reasonably describe fracture toughness required to arrest a dynamic crack propagation for a modern gas pipeline. However, the CTOA fracture criterion requires constraint-independent CTOA toughness against stable ductile crack growth. ASTM E3039 recommends a drop weight tearing test (DWTT) specimen for CTOA testing. Since a shallow crack is used, DWTT measured CTOA may depend on constraint level at the crack tip. To understand if it is the case, this paper evaluates the critical CTOA for a set of fracture toughness tests on single edge notched bend (SENB) specimens with shallow and deep cracks based on four CTOA estimation models. In which, the Ln(P)-LLD linear fit model is similar to that used by ASTM E3039 in the CTOA calculation. Fracture test data for X80 pipeline steel and HY80 structural steel are considered in the CTOA evaluation. The results show that the four CTOA models can determine a crack size-independent constant CTOA over stable ductile crack growth for the SENB specimens. As a result, CTOA determined by ASTM E3039 is constraint&#x2;independent and transferable to use for an actual crack propagating in a gas pipeline.

Zhu, Xian-Kui↗
DOE OSTI2022

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↗
DOE OSTI2024

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↗
DOE OSTI2018

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↗
DOE OSTI2025

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↗
DOE OSTI2022

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↗