Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Function 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,279 records · Page 71

On the sea-state bias of the Geosat altimeter

The sea-state bias in a satellite altimeter's range measurement is caused by the influence of ocean waves on the radar return pulse; it results in an estimate of sea level that is too low according to some function of the wave height. This bias is here estimated for Geosat by correlating collinear differences of altimetric sea-surface heights with collinear differences of significant wave heights (H1/3). Corrections for satellite orbit error are estimated simultaneously with the sea-state bias. Based on twenty 17-day repeat cycles of the Geosat Exact Repeat Mission, the solution for the sea-state bias is 2.6 + or - 0.2 percent of H1/3. The least-squares residuals, however, show a correlation with wind speed U, so the traditional model of the bias has been supplemented with a second term: H1/3 + alpha-2H1/3U. This second term produces a small, but statistically significant, reduction in variance of the residuals. Both systematic and random errors in H1/3 and U tend to bias the estimates of alpha-1 and alpha-2, which complicates comparisons of the results with ground-based measurements of the sea-state bias.

Ray, Richard D.↗

An error covariance model for sea surface topography and velocity derived from TOPEX/POSEIDON altimetry

In order to facilitate the use of satellite-derived sea surface topography and velocity oceanographic models, methodology is presented for deriving the total error covariance and its geographic distribution from TOPEX/POSEIDON measurements. The model is formulated using a parametric model fit to the altimeter range observations. The topography and velocity modeled with spherical harmonic expansions whose coefficients are found through optimal adjustment to the altimeter range residuals using Bayesian statistics. All other parameters, including the orbit, geoid, surface models, and range corrections are provided as unadjusted parameters. The maximum likelihood estimates and errors are derived from the probability density function of the altimeter range residuals conditioned with a priori information. Estimates of model errors for the unadjusted parameters are obtained from the TOPEX/POSEIDON postlaunch verification results and the error covariances for the orbit and the geoid, except for the ocean tides. The error in the ocean tides is modeled, first, as the difference between two global tide models and, second, as the correction to the present tide model, the correction derived from the TOPEX/POSEIDON data. A formal error covariance propagation scheme is used to derive the total error. Our global total error estimate for the TOPEX/POSEIDON topography relative to the geoid for one 10-day period is found tio be 11 cm RMS. When the error in the geoid is removed, thereby providing an estimate of the time dependent error, the uncertainty in the topography is 3.5 cm root mean square (RMS). This level of accuracy is consistent with direct comparisons of TOPEX/POSEIDON altimeter heights with tide gauge measurements at 28 stations. In addition, the error correlation length scales are derived globally in both east-west and north-south directions, which should prove useful for data assimilation. The largest error correlation length scales are found in the tropics. Errors in the velocity field are smallest in midlatitude regions. For both variables the largest errors caused by uncertainty in the geoid. More accurate representations of the geoid await a dedicated geopotential satellite mission. Substantial improvements in the accuracy of ocean tide models are expected in the very near future from research with TOPEX/POSEIDON data.

Tsaoussi, Lucia S.↗

Spacecraft drag-free technology development: On-board estimation and control synthesis

Estimation and control methods for a Drag-Free spacecraft are discussed. The functional and analytical synthesis of on-board estimators and controllers for an integrated attitude and translation control system is represented. The framework for detail definition and design of the baseline drag-free system is created. The techniques for solution of self-gravity and electrostatic charging problems are applicable generally, as is the control system development.

Key, R. W.↗

Direct evaluation of dynamical large-deviation rate functions using a variational ansatz

We describe a simple form of importance sampling designed to bound and compute large-deviation rate functions for time-extensive dynamical observables in continuous-time Markov chains. We start with a model, defined by a set of rates, and a time-extensive dynamical observable. We construct a reference model, a variational ansatz for the behavior of the original model conditioned on atypical values of the observable. Direct simulation of the reference model provides an upper bound on the large-deviation rate function associated with the original model, an estimate of the tightness of the bound, and, if the ansatz is chosen well, the exact rate function. The exact rare behavior of the original model does not need to be known in advance. We use this method to calculate rate functions for currents and counting observables in a set of network- and lattice models taken from the literature. Straightforward ansätze yield bounds that are tighter than bounds obtained from Level 2.5 of large deviations via approximations that involve uniform scalings of rates. Here, we show how to correct these bounds in order to recover the rate functions exactly. Our approach is complementary to more specialized methods and offers a physically transparent framework for approximating and calculating the likelihood of dynamical large deviations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

SU-F-T-116: Predicting IQ and the Risk of Hearing Loss Following Proton Versus Photon Radiotherapy for Pediatric Brain Tumor Patients

Purpose: The increased sparing of normal tissues in intensity modulated proton therapy (IMPT) compared to photon intensity modulated radiotherapy (IMRT) in brain tumor treatments should translate into improved neurocognitive outcomes. Models were used to estimate the intelligence quotient (IQ) and the risk of hearing loss 5 years post radiotherapy and to compare outcomes of proton against photon in pediatric brain tumors. Methods: Patients who had received radical IMRT were randomly selected from our retrospective database: 10 cases each of craniopharyngioma, ependymoma and medulloblastoma, and 20 cases of glioma. The existing planning CT and contours were used to generate IMPT plans. The RBE-corrected dose to brain structures and cochleas were calculated for both IMPT and IMRT. A model was applied to estimate IQ using a Markov chain Monte Carlo technique. The reported incidence of hearing loss as a function of cochlear dose was used to estimate the rate of occurrence. Results: The average brain dose was less in all IMPT plans compared to IMRT: ranging from a 6.7% reduction (P=0.003) in the case of medulloblastoma to 38% (P=0.007) for craniopharyngioma. This dose reduction translated into a gain in IQ of 1.9 points on average for protons vs photons for the whole cohort at 5 years post-treatment (P=0.011). In terms of specific diseases, the gains in IQ ranged from 0.8 points for medulloblastoma, to 2.7 points for craniopharyngioma. Hearing loss probability was evaluated on a per-ear-basis and was found to be systematically less for proton versus photon: overall 2.9% versus 7.2% (P < 0.001). Conclusion: A novel method was developed to predict neurocognitive outcomes in pediatric brain tumor patients on a case-by-case basis. A modest gain in IQ was systematically observed for proton in all patients. Given the uncertainties within the model used and our reinterpretation, these gains may be underestimated.

60 APPLIED LIFE SCIENCES↗

Global Autocorrelation Scales of the Partial Pressure of Oceanic CO2

A global database of approximately 1.7 million observations of the partial pressure of carbon dioxide in surface ocean waters (pCO2) collected between 1970 and 2003 is used to estimate its spatial autocorrelation structure. The patterns of the lag distance where the autocorrelation exceeds 0.8 is similar to patterns in the spatial distribution of the first baroclinic Rossby radius of deformation indicating that ocean circulation processes play a significant role in determining the spatial variability of pCO2. For example, the global maximum of the distance at which autocorrelations exceed 0.8 averages about 140 km in the equatorial Pacific. Also, the lag distance at which the autocorrelation exceed 0.8 is greater in the vicinity of the Gulf Stream than it is near the Kuroshio, approximately 50 km near the Gulf Stream as opposed to 20 km near the Kuroshio. Separate calculations for times when the sun is north and south of the equator revealed no obvious seasonal dependence of the spatial autocorrelation scales. The pCO2 measurements at Ocean Weather Station (OWS) 'P', in the eastern subarctic Pacific (50 N, 145 W) is the only fixed location where an uninterrupted time series of sufficient length exists to calculate a meaningful temporal autocorrelation function for lags greater than a few days. The estimated temporal autocorrelation function at OWS 'P', is highly variable. A spectral analysis of the longest four pCO2 time series indicates a high level of variability occurring over periods from the atmospheric synoptic to the maximum length of the time series, in this case 42 days. It is likely that a relative peak in variability with a period of 3-6 days is related to atmospheric synoptic period variability and ocean mixing events due to wind stirring. However, the short length of available time series makes identifying temporal relationships between pCO2 and atmospheric or ocean processes problematic.

Li, Zhen↗

Effects of error covariance structure on estimation of model averaging weights and predictive performance

When conducting model averaging for assessing groundwater conceptual model uncertainty, the averaging weights are often evaluated using model selection criteria such as AIC, AICc, BIC, and KIC (Akaike Information Criterion, Corrected Akaike Information Criterion, Bayesian Information Criterion, and Kashyap Information Criterion, respectively). However, this method often leads to an unrealistic situation in which the best model receives overwhelmingly large averaging weight (close to 100%), which cannot be justified by available data and knowledge. It was found in this study that this problem was caused by using the covariance matrix, CE, of measurement errors for estimating the negative log likelihood function common to all the model selection criteria. This problem can be resolved by using the covariance matrix, Cek, of total errors (including model errors and measurement errors) to account for the correlation between the total errors. An iterative two-stage method was developed in the context of maximum likelihood inverse modeling to iteratively infer the unknown Cek from the residuals during model calibration. The inferred Cek was then used in the evaluation of model selection criteria and model averaging weights. While this method was limited to serial data using time series techniques in this study, it can be extended to spatial data using geostatistical techniques. The method was first evaluated in a synthetic study and then applied to an experimental study, in which alternative surface complexation models were developed to simulate column experiments of uranium reactive transport. It was found that the total errors of the alternative models were temporally correlated due to the model errors. The iterative two-stage method using Cek resolved the problem that the best model receives 100% model averaging weight, and the resulting model averaging weights were supported by the calibration results and physical understanding of the alternative models. Using Cek obtained from the iterative two-stage method also improved predictive performance of the individual models and model averaging in both synthetic and experimental studies.

error correlation↗

Centrality dependence of J /ψ and ψ (2S) production and nuclear modification in p-Pb collisions at \( \sqrt{s_{\mathrm{NN}}} \) = 8 . 16 TeV

The inclusive production of the J/ ψ and ψ (2S) charmonium states is studied as a function of centrality in p-Pb collisions at a centre-of-mass energy per nucleon pair \( \sqrt{s_{\mathrm{NN}}} \) = 8 . 16 TeV at the LHC. The measurement is performed in the dimuon decay channel with the ALICE apparatus in the centre-of-mass rapidity intervals – 4 . 46 < y cms < – 2 . 96 (Pb-going direction) and 2 . 03 < y cms < 3 . 53 (p-going direction), down to zero transverse momentum ( p T ). The J/ ψ and ψ (2S) production cross sections are evaluated as a function of the collision centrality, estimated through the energy deposited in the zero degree calorimeter located in the Pb-going direction. The p T -differential J/ ψ production cross section is measured at backward and forward rapidity for several centrality classes, together with the corresponding average \( \left\langle {p}_{\mathrm{T}}\right\rangle \) and \( \left\langle {p}_{\mathrm{T}}^2\right\rangle \) values. The nuclear effects affecting the production of both charmonium states are studied using the nuclear modification factor. In the p-going direction, a suppression of the production of both charmonium states is observed, which seems to increase from peripheral to central collisions. In the Pb-going direction, however, the centrality dependence is different for the two states: the nuclear modification factor of the J/ ψ increases from below unity in peripheral collisions to above unity in central collisions, while for the ψ (2S) it stays below or consistent with unity for all centralities with no significant centrality dependence. The results are compared with measurements in p-Pb collisions at \( \sqrt{s_{\mathrm{NN}}} \) = 5 . 02 TeV and no significant dependence on the energy of the collision is observed. Finally, the results are compared with theoretical models implementing various nuclear matter effects.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

SU(2) Flavor Asymmetry of the Proton Sea in Chiral Effective Theory

We refine the computation of the $\bar{d}$ - $\bar{u}$ flavor asymmetry in the proton sea with a complementary effort to reveal the dynamics of pion exchange in high-energy processes. In particular, we discuss the efficacy of pion exchange models to simultaneously describe leading neutron electroproduction at HERA along with the $\bar{d}$ - $\bar{u}$ flavor asymmetry in the proton. A detailed χ 2 analysis of the ZEUS and H1 data, when combined with constraints on the pion flux from Drell-Yan data, allows regions of applicability of one-pion exchange to be delineated. Based on the fit results, we also address a possible estimate for leading proton structure functions in upcoming tagged deep-inelastic scattering experiments at Jefferson Lab on the deuteron with forward protons.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Microstructural Characterization and High Strain Rate Plastic Flow Behavior of SMAW Armox500T Steel Joints from Spherical Indentation Experiments

Static indentation and uniaxial compression tests have been conducted on Armox500T steel and its weldments to predict the constraint factor (CF), i.e., ratio of Meyer’s hardness to uniaxial flow stress. Series of dynamic indentation experiments were carried out at impact velocities ranging from 5 to 300 m/s to estimate dynamic hardness as a function of average strain. Subsequently, dynamic indentation (DI) test data and CF determined under static indentation conditions have been used to study the high strain rate plastic flow behavior of Armox500T weldments in comparison with the base metal. It was observed that flow stress for Armox500T and its weldments under dynamic loading conditions (10{sup 4} s{sup −1}) are significantly higher than flow stress measured under static loading conditions (10{sup −3} s{sup −1}). The plastic flow behavior computed from DI is in good agreement with the data evaluated through conventional split Hopkinson pressure bar technique. Further, the study of the microstructure of base metal and weldments by optical microscopy and SEM revealed a considerable variation in the microstructure.

36 MATERIALS SCIENCE↗

Land cover change-induced decline in terrestrial gross primary production over the conterminous United States from 2001 to 2016

As one of the most dynamic aspects of global environmental change, land cover change (LCC) has a profound impact on terrestrial carbon sequestration. However, LCC-induced carbon fluxes are still the most uncertain terms in global and regional carbon budgets. Ecosystem gross primary production (GPP) is the total carbon uptake by vegetation through photosynthesis, serving as a major control on ecosystem function and land carbon balance during and after the modification of the land surface. However, accurately capturing LCC-induced GPP changes requires both high-quality land cover data and controlling for variation driven by other environmental factors such as climate. In this study, we comprehensively examined the effects of LCC on annual GPP trends over the conterminous United States (CONUS) from 2001 to 2016 using the USGS National Land Cover Database, a remote sensing-driven ecosystem model, and the Google Earth Engine cloud computing platform. We designed a series of model experiments to identify LCC effects on GPP by controlling climate effects. During the study period, LCC exerted a strong negative effect on total GPP across the CONUS ([-2.2, -1.8] Tg C yr -2 ), while climate had smaller positive effects ([0.17, 0. 92] Tg C yr -2 ). The LCC-induced reduction of GPP was mainly caused by net forest loss ([-1.98, -1.39] Tg C yr -2 ) and urban expansion ([-2.03, -1.92] Tg C yr -2 ), but was partially offset by increases in crop area ([+0.66, +0.79] Tg C yr -2 ). Ensemble simulations from TRENDY did not capture the strong negative LCC influences on GPP, likely due to limitations of the adopted land use/cover data. Overall, our study provides a novel perspective on LCC-induced GPP changes, which could help to improve our understanding of ecosystem function changes and constrain the estimation of land carbon balance in the context of anthropogenic activity and climate change.

54 ENVIRONMENTAL SCIENCES↗

Laser-diagnostic mapping of temperature and soot statistics in a 2-m diameter turbulent pool fire

We present spatial profiles of temperature and soot-volume-fraction statistics from a sooting 2-m base diameter turbulent pool fire, burning a 10%-toluene / 90%-methanol fuel mixture. Dual-pump coherent anti-Stokes Raman scattering and laser-induced incandescence are utilized to obtain radial profiles of temperature and soot probability density functions (pdf) as well as estimates of temperature/soot joint statistics at three vertical heights above the surface of the methanol/toluene fuel pool. Results are presented both in the fuel vapor-dome region at ¼ base diameter and in the actively burning region at ½ and ¾ diameters above the fuel surface. The spatial evolution of the soot and temperature pdfs is discussed and profiles of the temperature and soot mean and rms statistics are provided. Joint temperature/soot statistics are presented as spatially resolved conditional averages across the fire plume, and in terms of a joint pdf obtained by including measurements from multiple spatial locations.

42 ENGINEERING↗

Hedonic Price Estimates of Lake Water Quality: Valued Attribute, Instrumental Variables, and Ecological-Economic Benefits

We estimate the effect of lake water quality on residential housing prices using cross-sectional data on water clarity from 113 lakes across the United States. An instrumental-variables approach is developed to address potential endogeneity bias in the water-clarity variable. Three lake-based physical variables serve as instrumental variables: total nitrogen concentration in lake water, total phosphorus concentration in lake water, and water temperature near the lake surface. The econometric results include three methodological findings: the instruments are valid and strong; the ordinary-least-squares estimate of the coefficient on the water-clarity variable is unbiased; and water clarity is the attribute of water quality valued in the housing market. The estimated water-clarity effect shows that a one-tenth of a meter change in water clarity leads to a one-percent change in housing price, or an elasticity of 0.20 at mean clarity. Coupling this effect with estimates of an ecological production function, a lake-specific benefit index is developed that shows the effect on housing values of bringing lakes into compliance with the US Environmental Protection Agency's regional recommendations for phosphorous concentration in lake waters.

54 ENVIRONMENTAL SCIENCES↗

Kilovoltage Intrafraction Monitoring for Prostate Intensity Modulated Arc Therapy: First Clinical Results

Purpose: Most linear accelerators purchased today are equipped with a gantry-mounted kilovoltage X-ray imager which is typically used for patient imaging prior to therapy. A novel application of the X-ray system is kilovoltage intrafraction monitoring (KIM), in which the 3-dimensional (3D) tumor position is determined during treatment. In this paper, we report on the first use of KIM in a prospective clinical study of prostate cancer patients undergoing intensity modulated arc therapy (IMAT). Methods and Materials: Ten prostate cancer patients with implanted fiducial markers undergoing conventionally fractionated IMAT (RapidArc) were enrolled in an ethics-approved study of KIM. KIM involves acquiring kV images as the gantry rotates around the patient during treatment. Post-treatment, markers in these images were segmented to obtain 2D positions. From the 2D positions, a maximum likelihood estimation of a probability density function was used to obtain 3D prostate trajectories. The trajectories were analyzed to determine the motion type and the percentage of time the prostate was displaced {>=}3, 5, 7, and 10 mm. Independent verification of KIM positional accuracy was performed using kV/MV triangulation. Results: KIM was performed for 268 fractions. Various prostate trajectories were observed (ie, continuous target drift, transient excursion, stable target position, persistent excursion, high-frequency excursions, and erratic behavior). For all patients, 3D displacements of {>=}3, 5, 7, and 10 mm were observed 5.6%, 2.2%, 0.7% and 0.4% of the time, respectively. The average systematic accuracy of KIM was measured at 0.46 mm. Conclusions: KIM for prostate IMAT was successfully implemented clinically for the first time. Key advantages of this method are (1) submillimeter accuracy, (2) widespread applicability, and (3) a low barrier to clinical implementation. A disadvantage is that KIM delivers additional imaging dose to the patient.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

An evaluation of tri-valent oxide (Cr 2 O 3 ) as a grain enlarging dopant for UO 2 nuclear fuels fabricated under reducing environment

A study was performed to evaluate the microstructure and crystallography of nominally 500–2000 Cr 2 O 3 -doped UO 2 fabricated in a temperature range of 1150–1750°C under reducing experimental conditions. We observed an increase in grain size of the samples with the increase in heat treating temperature as expected. For a given sintering temperature (1700–1750°C), an increase in the grain size was also observed with the increase in Cr 2 O 3 concentration up to a value of ~1000–1200 wppm. A decrease in fission gas release as a function of grain size was estimated for the Cr 2 O 3 -doped UO 2 samples assuming specified post-irradiation annealing conditions. A nearly linear decrease was obtained in the lattice parameter of the Cr 2 O 3 -doped UO 2 fcc phase with the increase in Cr 2 O 3 concentration, especially up to a nominal value of 1000 wppm. The lattice parameter decrease was also persistent with the increase in the average grain size as a result of addition of Cr 2 O 3 into the UO 2 lattice. An increase in the crystallite size and a decrease in the microstrain of the $fcc$ phase were observed with the increase in the average grain size of the samples, indicating a higher crystallinity of the Cr 2 O 3 -doped samples than that of the undoped UO 2 sample.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Uniform-in-phase-space data selection with iterative normalizing flows

Improvements in computational and experimental capabilities are rapidly increasing the amount of scientific data that are routinely generated. In applications that are constrained by memory and computational intensity, excessively large datasets may hinder scientific discovery, making data reduction a critical component of data-driven methods. Datasets are growing in two directions: the number of data points and their dimensionality. Whereas dimension reduction typically aims at describing each data sample on lower-dimensional space, the focus here is on reducing the number of data points. A strategy is proposed to select data points such that they uniformly span the phase-space of the data. The algorithm proposed relies on estimating the probability map of the data and using it to construct an acceptance probability. An iterative method is used to accurately estimate the probability of the rare data points when only a small subset of the dataset is used to construct the probability map. Instead of binning the phase-space to estimate the probability map, its functional form is approximated with a normalizing flow. Therefore, the method naturally extends to high-dimensional datasets. The proposed framework is demonstrated as a viable pathway to enable data-efficient machine learning when abundant data are available.

97 MATHEMATICS AND COMPUTING↗

Optimization of nonlinear turbulence in stellarators

We present new stellarator equilibria that have been optimized for reduced turbulent transport using nonlinear gyrokinetic simulations within the optimization loop. The optimization routine involves coupling the pseudo-spectral GPU-native gyrokinetic code GX with the stellarator equilibrium and optimization code DESC . Since using GX allows for fast nonlinear simulations, we directly optimize for reduced nonlinear heat fluxes. To handle the noisy heat flux traces returned by these simulations, we employ the simultaneous perturbation stochastic approximation (SPSA) method that only uses two objective function evaluations for a simple estimate of the gradient. We show several examples that optimize for both reduced heat fluxes and good quasi-symmetry as a proxy for low neoclassical transport. Finally, we run full transport simulations using the T3D stellarator transport code to evaluate the changes in the macroscopic profiles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗