Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “variance 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,243 records · Page 69

Aerodynamic parameters of the Navion airplane extracted from flight

An iterative method, which is characterized as a maximum-likelihood minimum-variance technique, was used to extract the aerodynamic parameters of a Navion airplane from flight data. The purposes were to compare the results with parameters obtained from wind-tunnel tests and with results obtained by analog matching the same data, and to develop techniques for application of the parameter extraction program. Results from the study showed that the parameter-extraction program can produce aerodynamic parameters which will permit close estimation of the aircraft time histories used in the extraction process. The program determined an estimate of the standard deviations of the states and parameters. These estimates were used to indicate how well the calculated states fit the flight data and the confidence in the values of the estimated parameters. The study also showed that the values of the parameters were affected by the data and mathematical model used during the extraction process. Because of the lack of confidence in the parameters extracted by use of some of the sets of data, several parameters were estimated by other methods. By using a combination of methods, a set of parameters which gave a fit to the data was obtained.

Suit, W. T.↗

Spectral Dimensionality and Scale of Urban Radiance

Characterization of urban radiance and reflectance is important for understanding the effects of solar energy flux on the urban environment as well as for satellite mapping of urban settlement patterns. Spectral mixture analyses of Landsat and Ikonos imagery suggest that the urban radiance field can very often be described with combinations of three or four spectral endmembers. Dimensionality estimates of Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) radiance measurements of urban areas reveal the existence of 30 to 60 spectral dimensions. The extent to which broadband imagery collected by operational satellites can represent the higher dimensional mixing space is a function of both the spatial and spectral resolution of the sensor. AVIRIS imagery offers the spatial and spectral resolution necessary to investigate the scale dependence of the spectral dimensionality. Dimensionality estimates derived from Minimum Noise Fraction (MNF) eigenvalue distributions show a distinct scale dependence for AVIRIS radiance measurements of Milpitas, California. Apparent dimensionality diminishes from almost 40 to less than 10 spectral dimensions between scales of 8000 m and 300 m. The 10 to 30 m scale of most features in urban mosaics results in substantial spectral mixing at the 20 m scale of high altitude AVIRIS pixels. Much of the variance at pixel scales is therefore likely to result from actual differences in surface reflectance at pixel scales. Spatial smoothing and spectral subsampling of AVIRIS spectra both result in substantial loss of information and reduction of apparent dimensionality, but the primary spectral endmembers in all cases are analogous to those found in global analyses of Landsat and Ikonos imagery of other urban areas.

Small, Christopher↗

Mitigating the effects of undersampling in weak lensing shear estimation with metacalibration

ABSTRACT metacalibration is a state-of-the-art technique for measuring weak gravitational lensing shear from well-sampled galaxy images. We investigate the accuracy of shear measured with metacalibration from fitting elliptical Gaussians to undersampled galaxy images. In this case, metacalibration introduces aliasing effects leading to an ensemble multiplicative shear bias about 0.01 for Euclid and even larger for the Roman Space Telescope, well exceeding the missions’ requirements. We find that this aliasing bias can be mitigated by computing shapes from weighted moments with wider Gaussians as weight functions, thereby trading bias for a slight increase in variance of the measurements. We show that this approach is robust to the point-spread function in consideration and meets the stringent requirements of Euclid for galaxies with moderate to high signal-to-noise ratios. We therefore advocate metacalibration as a viable shear measurement option for weak lensing from upcoming space missions.

Kannawadi, Arun↗

Deep Learning enabled spectral energy conversion for in situ exposure measurements

A detector-specific deep learning (DL) approach is presented for spectra-to-exposure conversion using large-format sodium iodide (NaI(Tl)) detectors deployed for in situ environmental radiation measurements in emergency response scenarios. Accurate determination of exposure from NaI spectra is challenging due to poor energy resolution, partial energy absorption, and the strong sensitivity of traditionally deployed analytical conversion methods to calibrated source geometry and pre-deployment assumptions. Here, to address these limitations, a multi-layer perceptron model was trained on a hybrid in situ /Monte Carlo dataset constructed to span a broad range of photon energies, spatial extents, and realistic deployment variability, representative of general in situ emergency response conditions. The DL model was evaluated against commonly fielded analytical approaches under matched simulation conditions, including a single-factor method, a G-function method, and a modeled pressurized ion chamber (PIC) baseline. This study was intentionally computational in scope to enable controlled, like-for-like comparisons between conversion techniques while minimizing confounding real-world variability. Comparison to the modeled PIC provides contextual benchmarking and is not intended as a field inter-comparison with deployed instruments. Across the evaluated 20 keV to 3 MeV energy range, the DL approach consistently exhibited higher accuracy and reduced variance relative to the analytical methods against a deterministically calculated exposure. This may indicate improved robustness to spectral complexity without reliance on source-, geometric-, or spectral region-specific optimization. While results do not represent real-world validation, the presented work demonstrates that deep learning may effectively learn the nonlinear detector response-to-exposure relationship for asymmetric NaI(Tl) detectors and offers a promising pathway for improving in situ exposure estimation using spectroscopic systems already integrated into initial real-time emergency response operations.

61 RADIATION PROTECTION AND DOSIMETRY↗

Searching for hidden mirror symmetries in CMB fluctuations from WMAP 7 year maps

We search for hidden mirror symmetries at large angular scales in the WMAP 7 year Internal Linear Combination map of CMB temperature anisotropies using global pixel based estimators introduced for this aim. Two different axes are found for which the CMB intensity pattern is anomalously symmetric (or anti-symmetric) under reflection with respect to orthogonal planes at the 99.84(99.96)% CL (confidence level), if compared to a result for an arbitrary axis in simulations without the symmetry. We have verified that our results are robust to the introduction of the galactic mask. The direction of such axes is close to the CMB kinematic dipole and nearly orthogonal to the ecliptic plane, respectively. If instead the real data are compared to those in simulations taken with respect to planes for which the maximal mirror symmetry is generated by chance, the confidence level decreases to 92.39(76.65)%. But when the effect in question translates into the anomalous alignment between normals to planes of maximal mirror (anti)-symmetry and these natural axes mentioned. We also introduce the representation of the above estimators in the harmonic domain, confirming the results obtained in the pixel one. The symmetry anomaly is shown to be almost entirely due to low multipoles, so it may have a cosmological and even primordial origin. Contrary, the anti-symmetry one is mainly due to intermediate multipoles that probably suggests its non-fundamental nature. We have demonstrated that these anomalies are not connected to the known issue of the low variance in WMAP observations and we have checked that axially symmetric parts of these anomalies are small, so that the axes are not the symmetry ones.

79 ASTRONOMY AND ASTROPHYSICS↗

Reducing uncertainties in satellite estimates of aerosol–cloud interactions over the subtropical ocean by integrating vertically resolved aerosol observations

Satellite quantification of aerosol effects on clouds relies on aerosol optical depth (AOD) as a proxy for aerosol concentration or cloud condensation nuclei (CCN). However, the lack of error characterization of satellite-based results hampers their use for the evaluation and improvement of global climate models. We show that the use of AOD for assessing aerosol–cloud interactions (ACIs) is inadequate over vast oceanic areas in the subtropics. Instead, we postulate that a more physical approach that consists of matching vertically resolved aerosol data from the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) satellite at the cloud-layer height with Moderate Resolution Imaging Spectroradiometer (MODIS) Aqua cloud retrievals reduces uncertainties in satellite-based ACI estimates. Combined aerosol extinction coefficients (σ) below cloud top (σBC) from the Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) and cloud droplet number concentrations (N d ) from MODIS Aqua yield high correlations across a broad range of σBC values, with σBC quartile correlations ≥0.78. In contrast, CALIOP-based AOD yields correlations with MODIS Nd of 0.54–0.62 for the two lower AOD quartiles. Moreover, σBC explains 41% of the spatial variance in MODIS N d , whereas AOD only explains 17%, primarily caused by the lack of spatial covariability in the eastern Pacific. Compared with σBC, near-surface σ weakly correlates in space with MODIS N d , accounting for a 16% variance. It is concluded that the linear regression calculated from ln(N d )–ln(σBC) (the standard method for quantifying ACIs) is more physically meaningful than that derived from the Nd–AOD pair.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Modeling the Spin Equilibrium of Neutron Stars in LMXBs Without Gravitational Radiation

In this paper we discuss the spin-equilibrium of accreting neutron stars in LMXBs. We demonstrate that, when combined with a naive spin-up torque, the observed data leads to inferred magnetic fields which are at variance with those of galactic millisecond radiopulsars. This indicates the need for either additional spin-down torques (eg. gravitational radiation) or an improved accretion model. We show that a simple consistent accretion model can be arrived at by accounting for radiation pressure in rapidly accreting systems (above a few percent of the Eddington accretion rate). In our model the inner disk region is thick and significantly sub-Keplerian, and the estimated equilibrium periods are such that the LMXB neutron stars have properties that accord well with the galactic millisecond radiopulsar sample. The implications for future gravitational-wave observations are also discussed briefly.

Andersson, N.↗

A Conditional Autoencoder for Galaxy Photometric Parameter Estimation

Astronomical photometric surveys routinely image billions of galaxies, and traditionally infer the parameters of a parametric model for each galaxy. This approach has served us well, but the computational expense of deriving a full posterior probability distribution function is a challenge for increasingly ambitious surveys. In this paper, we use deep learning methods to characterize galaxy images, training a conditional autoencoder on mock data. The autoencoder can reconstruct and denoise galaxy images via a latent space engineered to include semantically meaningful parameters, such as brightness, location, size, and shape. Our model recovers galaxy fluxes and shapes on mock data with a lower variance than the Hyper Suprime-Cam photometry pipeline, and returns reasonable answers even for inputs outside the range of its training data. When applied to data in the training range, the regression errors on all extracted parameters are nearly unbiased with a variance near the Cramr-Rao bound.

79 ASTRONOMY AND ASTROPHYSICS↗

Potential and limitations of inferring ecosystem photosynthetic capacity from leaf functional traits

Abstract The aim of this study was to systematically analyze the potential and limitations of using plant functional trait observations from global databases versus in situ data to improve our understanding of vegetation impacts on ecosystem functional properties ( EFP s). Using ecosystem photosynthetic capacity as an example, we first provide an objective approach to derive robust EFP estimates from gross primary productivity ( GPP ) obtained from eddy covariance flux measurements. Second, we investigate the impact of synchronizing EFP s and plant functional traits in time and space to evaluate their relationships, and the extent to which we can benefit from global plant trait databases to explain the variability of ecosystem photosynthetic capacity. Finally, we identify a set of plant functional traits controlling ecosystem photosynthetic capacity at selected sites. Suitable estimates of the ecosystem photosynthetic capacity can be derived from light response curve of GPP responding to radiation (photosynthetically active radiation or absorbed photosynthetically active radiation). Although the effect of climate is minimized in these calculations, the estimates indicate substantial interannual variation of the photosynthetic capacity, even after removing site‐years with confounding factors like disturbance such as fire events. The relationships between foliar nitrogen concentration and ecosystem photosynthetic capacity are tighter when both of the measurements are synchronized in space and time. When using multiple plant traits simultaneously as predictors for ecosystem photosynthetic capacity variation, the combination of leaf carbon to nitrogen ratio with leaf phosphorus content explains the variance of ecosystem photosynthetic capacity best (adjusted R 2 = 0.55). Overall, this study provides an objective approach to identify links between leaf level traits and canopy level processes and highlights the relevance of the dynamic nature of ecosystems. Synchronizing measurements of eddy covariance fluxes and plant traits in time and space is shown to be highly relevant to better understand the importance of intra‐ and interspecific trait variation on ecosystem functioning.

54 ENVIRONMENTAL SCIENCES↗

The NANOGrav nine-year data set: Noise budget for pulsar arrival times on intraday timescales

The use of pulsars as astrophysical clocks for gravitational wave (GW) experiments demands the highest possible timing precision. Pulse times of arrival (TOAs) are limited by stochastic processes that occur in the pulsar itself, along the line of sight through the interstellar medium, and in the measurement process. On timescales of seconds to hours, the TOA variance exceeds that from template-fitting errors due to additive noise. We assess contributions to the total variance from two additional effects: amplitude and phase jitter intrinsic to single pulses and changes in the interstellar impulse response from scattering. The three effects have different dependencies on time, frequency, and pulse signal-to-noise ratio. We use data on 37 pulsars from the North American Nanohertz Observatory for GWs to assess the individual contributions to the overall intraday noise budget for each pulsar. We detect jitter in 22 pulsars and estimate the average value of rms jitter in our pulsars to be ∼1% of pulse phase. We examine how jitter evolves as a function of frequency and find evidence for evolution. Finally, we compare our measurements with previous noise parameter estimates and discuss methods to improve GW detection pipelines.

79 ASTRONOMY AND ASTROPHYSICS↗

Measuring kinetic energy changes in the mesoscale with low acquisition rates

We report on the measurement of the average kinetic energy changes in isothermal and non-isothermal quasistatic processes in the mesoscale, realized with a Brownian particle trapped with optical tweezers. Our estimation of the kinetic energy change allows to access to the full energetic description of the Brownian particle. Kinetic energy estimates are obtained from measurements of the mean square velocity of the trapped bead sampled at frequencies several orders of magnitude smaller than the momentum relaxation frequency. The velocity is tuned applying a noisy electric field that modulates the amplitude of the fluctuations of the position and velocity of the Brownian particle, whose motion is equivalent to that of a particle in a higher temperature reservoir. Additionally, we show that the dependence of the variance of the time-averaged velocity on the sampling frequency can be used to quantify properties of the electrophoretic mobility of a charged colloid. Our method could be applied to detect temperature gradients in inhomogeneous media and to characterize the complete thermodynamics of biological motors and of artificial micro and nanoscopic heat engines.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Data-driven predictive probability density function control of fiber length stochastic distribution shaping in refining process

Pulp is the most important raw material in paper in paper industries, whose Fiber Length Stochastic Distribution (FLSD) shaping directly determines the energy consumption and paper quality of the subsequent papermaking processes. However, the mean and variance are insufficient to describe the output FLSD shaping, which display non-Gaussian distribution properties. Therefore, the traditional control, method based on the mean and variance of fiber length is difficult to control the FLSD shaping effectively. In this paper, a novel data-driven predictive probability density function (PDF) control method is proposed for output FLSD shaping in refining process. Primarily, in order to improve the approximation accuracy of the PDF employing RBF neural network, the parameters (center value and width) of each RBF basis functions are tuned via utilizing iterative learning control (ILC) law, and the corresponding estimations of weights law can be obtained. Secondly, considering that the conventional linear model of weights vector has drawbacks of low accuracy and weak generalization ability, random vector functional link networks (RVFLNs) based data driven nonlinear modelling method is employed to characterize the prediction model between the input variables and weights vector. Finally, in order to reduce the randomness of the output FLSD, a minimum entropy control method under mean constraint is employed to design the predictive PDF controller for the output FLSD shaping with the help of the established nonlinear stochastic distribution model. Both simulations and experiments are given to illustrate the effectiveness and practicability of the proposed method.

Probability debsity function, Fiber length distrib↗

Fitting dynamic models to the Geosat sea level observations in the tropical Pacific Ocean. I - A free wave model

Free, equatorially trapped sinusoidal wave solutions to a linear model on an equatorial beta plane are used to fit the Geosat altimetric sea level observations in the tropical Pacific Ocean. The Kalman filter technique is used to estimate the wave amplitude and phase from the data. The estimation is performed at each time step by combining the model forecast with the observation in an optimal fashion utilizing the respective error covariances. The model error covariance is determined such that the performance of the model forecast is optimized. It is found that the dominant observed features can be described qualitatively by basin-scale Kelvin waves and the first meridional-mode Rossby waves. Quantitatively, however, only 23 percent of the signal variance can be accounted for by this simple model.

Fu, Lee-Lueng↗

Measurement System Characterization in the Presence of Measurement Errors

In the calibration of a measurement system, data are collected in order to estimate a mathematical model between one or more factors of interest and a response. Ordinary least squares is a method employed to estimate the regression coefficients in the model. The method assumes that the factors are known without error; yet, it is implicitly known that the factors contain some uncertainty. In the literature, this uncertainty is known as measurement error. The measurement error affects both the estimates of the model coefficients and the prediction, or residual, errors. There are some methods, such as orthogonal least squares, that are employed in situations where measurement errors exist, but these methods do not directly incorporate the magnitude of the measurement errors. This research proposes a new method, known as modified least squares, that combines the principles of least squares with knowledge about the measurement errors. This knowledge is expressed in terms of the variance ratio - the ratio of response error variance to measurement error variance.

Commo, Sean A.↗

Estimation of uncertainties for subcritical benchmark measurements - 14387

Subcritical measurements are used for many applications including nonproliferation, treaty verification, in-situ monitoring for criticality, and reactor subcriticality monitoring. Most analysis methods involve correlated neutron detection, based upon the property that fission events can create multiple neutrons that are born at nearly the same time (on a scale of 10{sup -13} or 10{sup -14} seconds). Many different time-correlated methods have been used since the 1950's and are still widely utilized today. This particular work describes an uncertainty analysis of measured data using the Hage-Cifarelli formalism of the Feynman Variance-to-Mean method. The Feynman Variance-to-Mean method was introduced by Feynman et al. in 1956 and is based upon the differences between the detector count data and the data that would be expected from measurements taken from a Poissonian source. The Hage-Cifarelli formalism relates moments of the Feynman histograms to several parameters of a system: leakage multiplication (M{sub L}), spontaneous fission rate (F{sub S}), (α,n) neutron emission rate (S{sub α}), detector efficiency (ε), and the moments of the number of neutrons emitted per fission (ν). The first part of the uncertainty method described in this work will show how to determine uncertainties of the first (R{sub 1}) and second (R{sub 2}) reduced factorial moments of the Feynman histograms. The second part of the method will describe how to determine the uncertainties in the other parameters given uncertainties in the reduced factorial moments; in particular, this work will focus on determining the uncertainties in M{sub L} and F{sub S}. This work will utilize measured data from a recent benchmark evaluation with a 4.5 kg sphere of α-phase plutonium to validate the uncertainty analysis. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Image-derived and arterial blood sampled input functions for quantitative PET imaging of the angiotensin II subtype 1 receptor in the kidney

Purpose: The radioligand {sup 11}C-KR31173 has been introduced for positron emission tomography (PET) imaging of the angiotensin II subtype 1 receptor in the kidney in vivo. To study the biokinetics of {sup 11}C-KR31173 with a compartmental model, the input function is needed. Collection and analysis of arterial blood samples are the established approach to obtain the input function but they are not feasible in patients with renal diseases. The goal of this study was to develop a quantitative technique that can provide an accurate image-derived input function (ID-IF) to replace the conventional invasive arterial sampling and test the method in pigs with the goal of translation into human studies. Methods: The experimental animals were injected with [{sup 11}C]KR31173 and scanned up to 90 min with dynamic PET. Arterial blood samples were collected for the artery derived input function (AD-IF) and used as a gold standard for ID-IF. Before PET, magnetic resonance angiography of the kidneys was obtained to provide the anatomical information required for derivation of the recovery coefficients in the abdominal aorta, a requirement for partial volume correction of the ID-IF. Different image reconstruction methods, filtered back projection (FBP) and ordered subset expectation maximization (OS-EM), were investigated for the best trade-off between bias and variance of the ID-IF. The effects of kidney uptakes on the quantitative accuracy of ID-IF were also studied. Biological variables such as red blood cell binding and radioligand metabolism were also taken into consideration. A single blood sample was used for calibration in the later phase of the input function. Results: In the first 2 min after injection, the OS-EM based ID-IF was found to be biased, and the bias was found to be induced by the kidney uptake. No such bias was found with the FBP based image reconstruction method. However, the OS-EM based image reconstruction was found to reduce variance in the subsequent phase of the ID-IF. The combined use of FBP and OS-EM resulted in reduced bias and noise. After performing all the necessary corrections, the areas under the curves (AUCs) of the AD-IF were close to that of the AD-IF (average AUC ratio =1 ± 0.08) during the early phase. When applied in a two-tissue-compartmental kinetic model, the average difference between the estimated model parameters from ID-IF and AD-IF was 10% which was within the error of the estimation method. Conclusions: The bias of radioligand concentration in the aorta from the OS-EM image reconstruction is significantly affected by radioligand uptake in the adjacent kidney and cannot be neglected for quantitative evaluation. With careful calibrations and corrections, the ID-IF derived from quantitative dynamic PET images can be used as the input function of the compartmental model to quantify the renal kinetics of {sup 11}C-KR31173 in experimental animals and the authors intend to evaluate this method in future human studies.

60 APPLIED LIFE SCIENCES↗

Exploration with Scalable Gaussian Process Reinforcement Learning

Exploration is a challenging problem in reinforcement learning (RL), especially in environments with sparse rewards. Quantifying and utilizing the parametric uncertainty has been shown to be paramount for successful exploration [Osband et al., 2018]. Bayesian, or approximately Bayesian, methods present a principled means of estimating the parametric uncertainty in RL problems. Gaussian processes, nonparametric Bayesian models, are often impractical due to poor scalability and computational bottlenecks. We introduce a scalable Gaussian process RL (GPRL) method which directly induces sparsity in the covariance matrix to facilitate faster computation. This is a departure from previous GPRL methods which instead rely on data reduction and subsampling. We compare various covariance-based exploration techniques (Thompson sampling, upper confidence bound, and probabilistic maximum variance) which leverage our scalable GP framework in sparse reward environments. Finally, we show favorable comparison against the bootstrapped deep Q-Network.

97 MATHEMATICS AND COMPUTING↗

Spectral analysis of GEOS-3 altimeter data and frequency domain collocation

The mathematical background in spectral analysis as applied to geodetic applications is summarized. The resolution (cut-off frequency) of the GEOS 3 altimeter data is examined by determining the shortest wavelength (corresponding to the cut-off frequency) recoverable. The data from some 18 profiles are used. The total power (variance) in the sea surface topography with respect to the reference ellipsoid as well as with respect to the GEM-9 surface is computed. A fast inversion algorithm for matrices of simple and block Toeplitz matrices and its application to least squares collocation is explained. This algorithm yields a considerable gain in computer time and storage in comparison with conventional least squares collocation. Frequency domain least squares collocation techniques are also introduced and applied to estimating gravity anomalies from GEOS 3 altimeter data. These techniques substantially reduce the computer time and requirements in storage associated with the conventional least squares collocation. Numerical examples given demonstrate the efficiency and speed of these techniques.

Eren, K.↗