Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “yield 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 595 records · Page 33

Welch Method and Bootstrapping Applied to Subcritical Gamma Noise

We measured the prompt neutron decay constant 𝛼 of the CROCUS zero-power reactor at the Swiss Federal Institute of Technology Lausanne using cross-power spectral density (CPSD) analysis of gamma-gamma correlations from two trans-stilbene organic scintillators positioned near the reactor core. We measured critical and subcritical states, with water levels ranging from 960 mm (critical) to 800 mm (𝜌=−1.4 $ subcritical). Our analysis used the Welch method, dividing signal segments for fast Fourier transform (FFT) frequency analysis and applying bootstrapping uncertainty quantification that uses Welch-defined segments. Results demonstrated a clear increase in the measured 𝛼 as reactor reactivity decreased, distinguishing critical from subcritical conditions. At the 960-mm critical level, 𝛼 was estimated at 155.9 ± 0.7 s −1 , and for the 800-mm subcritical level, 𝛼 increased significantly to 367.3 ± 6.9 s –1 . A linear regression of subcritical states yielded a critical estimate of 154.0 ± 3.1 s –1 , aligning with the static 𝛼 estimate at critical. The bootstrapping method produced normally distributed 𝛼 estimates, confirming data consistency. The gamma CPSD 𝛼 estimates clearly distinguish reactor states and improve monitoring of zero-power reactors. The future deployment of modular and microreactors as potential candidates for noise analysis is demonstrated in CROCUS, particularly zero-power mock-ups of new designs. The improvement of noise analysis in the subcritical domain from this work will support experimental data for reactor deployment and procedure.

CROCUS↗

Computer Model Calibration with Time Series Data Using Deep Learning and Quantile Regression

Computer models play a key role in many scientific and engineering problems. One major source of uncertainty in computer model experiments is input parameter uncertainty. Computer model calibration is a formal statistical procedure to infer input parameters by combining information from model runs and observational data. The existing standard calibration framework suffers from inferential issues when the model output and observational data are high-dimensional dependent data, such as large time series, due to the difficulty in building an emulator and the nonidentifiability between effects from input parameters and data-model discrepancy. To overcome these challenges, we, in this work, propose a new calibration framework based on a deep neural network (DNN) with long short-term memory layers that directly emulates the inverse relationship between the model output and input parameters. Adopting the “learning with noise” idea, we train our DNN model to filter out the effects from data-model discrepancy on input parameter inference. We also formulate a new way to construct interval predictions for DNN using quantile regression to quantify the uncertainty in input parameter estimates. Through a simulation study and real data application with the Weather Research and Forecasting Model Hydrological modeling system (WRF-Hydro), we show our approach can yield accurate point estimates and well-calibrated interval estimates for input parameters.

97 MATHEMATICS AND COMPUTING↗

Notes on Regression Analysis for Radar Parameter Estimation

A fundamental task of radar, beyond merely detecting a target, is to estimate some parameters associated with it. For example, this might include range, direction, velocity, etc. In any case, multiple measurements, often noisy, need to be processed to yield a ‘best estimate’ of the parameter. A common mathematical method for doing so is called “Regression” analysis. The goal is to minimize the expected squared error in the estimate. Even when alternate algorithms are considered, the least squared-error regression analysis is the benchmark against which alternatives are compared.

47 OTHER INSTRUMENTATION↗

High deuteron and neutron yields from the interaction of a petawatt laser with a cryogenic deuterium jet

A compact high-flux, short-pulse neutron source would have applications from nuclear astrophysics to cancer therapy. Laser-driven neutron sources can achieve fluxes much higher than spallation and reactor neutron sources by reducing the volume and time in which the neutron-producing reactions occur by orders of magnitude. We report progress towards an efficient laser-driven neutron source in experiments with a cryogenic deuterium jet on the Texas Petawatt laser. Neutrons were produced both by laser-accelerated multi-MeV deuterons colliding with Be and mixed metallic catchers and by d (d,n) 3 He fusion reactions within the jet. We observed deuteron yields of 10 13 /shot in quasi-Maxwellian distributions carrying ~8-10% of the input laser energy. We obtained neutron yields greater than 10 10 /shot and found indications of a deuteron-deuteron fusion neutron source with high peak flux (>10 22 cm -2 s -1 ). The estimated fusion neutron yield in our experiment is one order of magnitude higher than any previous laser-induced dd fusion reaction. Though many technical challenges will have to be overcome to convert this proof-of-principle experiment into a consistent ultra-high flux neutron source, the neutron fluxes achieved here suggest laser-driven neutron sources can support laboratory study of the rapid neutron-capture process, which is otherwise thought to occur only in astrophysical sites such as core-collapse supernova, and binary neutron star mergers.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Blowing Snow at McMurdo Station, Antarctica During the AWARE Field Campaign: Multi-Instrument Observations of Blowing Snow

In polar regions, blowing snow (BLSN) can play a substantial role in regional thermodynamics and radiation properties. The extent of the impact depends upon the depth of the BLSN layer, a property that is difficult to measure and not well understood. The West Antarctic Radiation Experiment (AWARE) Field Campaign saw the deployment of a large suite of observational instruments to McMurdo Station, Antarctica, allowing for in-depth investigation of BLSN. During the year of the campaign, BLSN occurred ~7.4% of the time. In this study, additional remote sensing observations are used to supplement an existing ceilometer-based BLSN characterization algorithm creating a multi-instrument (MI) depth estimation algorithm to yield layer depths with more certainty. Backscatter coefficient and linear depolarization ratios from the micropulse and high spectral resolution lidar (HSRL) are incorporated along with color ratio derived from HSRL backscatter and Ka-band radar reflectivity observations. This provided estimates of BLSN depth that were less impacted by noise/artifacts in any one set of observations. When applied to the data from AWARE, the MI algorithm yielded an average depth of 168.2 m with ~60% of plumes less than 200 m in depth. Events occurring with precipitation saw the most value in the additional instruments, with the suite of observations allowing for more separation of BLSN particles from fallstreaks.

54 ENVIRONMENTAL SCIENCES↗

Using spatial coherence to retrieve cloud properties - ERBE experience and CERES applications

Cloud properties for individual ERBE scanner fields of view are obtained from simultaneous, collocated imagery data. The cloud properties are used to assess the performance of the ERBE maximum likelihood estimate of cloud conditions which ultimately influences estimates of the radiative fluxes at the top of the atmosphere. During EOS the analysis of MODIS observations within the fields of view of the CERES instruments promises to yield, in addition, estimates of surface radiative fluxes as well as the vertical distribution of radiative heating within the atmosphere.

Coakley, James A., Jr.↗

Determination of the Thermal Offset of the Eppley Precision Spectral Pyranometer

Eppley's precision spectral pyranometer (PSP) is used in networks around the world to measure downwelling diffuse and global solar irradiance at the surface of the Earth. In recent years several studies have shown significant discrepancy between irradiances measured by pyranometers and those computed by atmospheric radiative transfer models. Pyranometer measurements have been questioned because observed diffuse irradiances sometimes are below theoretical minimum values for a pure molecular atmosphere, and at night the instruments often produce nonzero signals ranging between + 5 and - 10 W/sq m. We install thermistor sondes in the body of a PSP as well as on its inner dome to monitor the temperature gradients within the instrument, and we operate a pyrgeometer (PIR) instrument side by side with the PSP. We derive a relationship between the PSP output and thermal radiative exchange by the dome and the detector and a relationship between the PSP output and the PIR thermopile output (net-IR). We determine the true PSP offset by quickly capping the instrument at set time intervals. For a ventilated and shaded PSP, the thermal offset can reach - 15 W/sq m under clear skies, whereas it remains close to zero for low overcast clouds. We estimate the PSP thermal offset by two methods: (1) using the PSP temperatures and (2) using the PIR net-IR signal. The offset computed from the PSP temperatures yields a reliable estimate of the true offset (+/- 1 W/sq m). The offset computed from net-IR is consistent with the true offset at night and under overcast skies but predicts only part of the true range under clear skies.

Haeffelin, Martial↗

The helium abundance of Jupiter from Voyager

Full disk measurements recorded 31 days before the Voyager 1 encounter with Jupiter by the radiometer of the infrared instrument, IRIS, indicate a geometric albedo of 0.274 + or - 0.013. Combining this measurement with the Pioneer derived phase integral of 1.25 and our error estimate of 0.1 yields a Jovian Bond albedo of 0.343 + or - 0.032. Infrared spectra recorded at the same time by the Michelson interferometer, along with a model extrapolation to low wave numbers not covered by the instrument, yield a thermal emission of (1.359 + or - 0.014) .001 W cm to the (-2) power. As in the case of the albedo measurement, the quoted errors in the emission measurement reflect estimates of systematic effects and are uncertain while the random component is negligible. From these measurements the internal heat flux of Jupiter is estimated to be (5.444 + or - 0.425) .0001 W cm to the (-2) power, and the energy balance defined as the ratio of emitted thermal to absorbed solar energy is 1.668 + or - 0.085.

Gautier, D.↗

Estimation of new production in the ocean by compound remote sensing

Compound remote sensing, a technique for estimating oceanic new production using remotely sensed data on ocean color and temperature, is discussed. This approach, which depends on parameterizations developed from ship observations as well as on satellite data, yields more representative estimates of the large-scale average new production than those calculated from ship data alone. The approach is illustated with data for the Georges Bank.

Sathyendranath, Shubha↗

Improvements to GNSS Positioning in Challenging Environments by 3DMA Lidar Informed Selective Satellites Usage

The use of global navigation satellite systems (GNSS) for position estimation tends to yield poor results when operating inside of an urban canyon due to large obstructions (e.g., buildings) that disrupt signals as they travel from a satellite to a receiver resulting in a position estimate that may significantly fluctuate in magnitude and direction. Identifying and removing signals that are non-line-of-sight (NLOS) to the receiver and only using signals that are line-of-sight (LOS) can improve the estimated position. However, quickly and accurately determining the LOS status of each measurement can be challenging without additional information about the operating environment. Use of publicly available lidar data can be used to incorporate techniques, such as 3D-mapping-aided (3DMA), to estimate the LOS status of satellites and augment the position solution accordingly. To complicate the issue, the error on the GNSS position estimate in an urban canyon is often so large that is it not sufficient to use as an approximate location for LOS prediction. That is, at times the calculated GNSS solution is not representative of the true location and cannot be used to accurately predict which satellites are within LOS due to the difference in the physical geometry associated with the two locations. This paper explores the use of a GNSS/inertial fused position solution as the initial position estimate for predicting which satellites are within LOS in an urban environment and the impact that removal of predicted NLOS satellites has on the GNSS position solution.

GPS↗

Developing “inverted-corona” fusion targets as high-fluence neutron sources

We present experimental studies of inverted-corona targets as neutron sources at the OMEGA Laser Facility and the National Ignition Facility (NIF). Laser beams are directed onto the inner walls of a capsule via laser-entrance holes (LEHs), heating the target interior to fusion conditions. The fusion fuel is provided either as a wall liner, e.g., deuterated plastic (CD), or as a gas fill, e.g., D 2 gas. Such targets are robust to low-mode drive asymmetries, allowing for single-sided laser drive. On OMEGA, 1.8-mm-diameter targets with either a 10-μm CD liner or up to 2 atm of D 2 -gas fill were driven with up to 18 kJ of laser energy in a 1-ns square pulse. Neutron yields of up to 1.5 × 1010 generally followed expected trends with fill pressure or laser energy, although the data imply some mix of the CH wall into the fusion fuel for either design. Comparable performance was observed with single-sided (1x LEH) or double-sided (2x LEH) drive. NIF experiments tested the platform at scaled up dimensions and energies, combining a 15-μm CD liner and a 3-atm D2-gas fill in a 4.5-mm diameter target, laser-driven with up to 330 kJ. Neutron yields up to 2.6 × 10 12 were measured, exceeding the scaled yield expectation from the OMEGA data. The observed energy scaling on the NIF implies that the neutron production is gas dominated, suggesting a performance boost from using deuterium–tritium (DT) gas. In this work, we estimate that neutron yields exceeding 10 14 should be readily achievable using a modest laser drive of ~300 kJ with a DT fill.

47 OTHER INSTRUMENTATION↗

Iterative-Transform Phase Retrieval Using Adaptive Diversity

A phase-diverse iterative-transform phase-retrieval algorithm enables high spatial-frequency, high-dynamic-range, image-based wavefront sensing. [The terms phase-diverse, phase retrieval, image-based, and wavefront sensing are defined in the first of the two immediately preceding articles, Broadband Phase Retrieval for Image-Based Wavefront Sensing (GSC-14899-1).] As described below, no prior phase-retrieval algorithm has offered both high dynamic range and the capability to recover high spatial-frequency components. Each of the previously developed image-based phase-retrieval techniques can be classified into one of two categories: iterative transform or parametric. Among the modifications of the original iterative-transform approach has been the introduction of a defocus diversity function (also defined in the cited companion article). Modifications of the original parametric approach have included minimizing alternative objective functions as well as implementing a variety of nonlinear optimization methods. The iterative-transform approach offers the advantage of ability to recover low, middle, and high spatial frequencies, but has disadvantage of having a limited dynamic range to one wavelength or less. In contrast, parametric phase retrieval offers the advantage of high dynamic range, but is poorly suited for recovering higher spatial frequency aberrations. The present phase-diverse iterative transform phase-retrieval algorithm offers both the high-spatial-frequency capability of the iterative-transform approach and the high dynamic range of parametric phase-recovery techniques. In implementation, this is a focus-diverse iterative-transform phaseretrieval algorithm that incorporates an adaptive diversity function, which makes it possible to avoid phase unwrapping while preserving high-spatial-frequency recovery. The algorithm includes an inner and an outer loop (see figure). An initial estimate of phase is used to start the algorithm on the inner loop, wherein multiple intensity images are processed, each using a different defocus value. The processing is done by an iterative-transform method, yielding individual phase estimates corresponding to each image of the defocus-diversity data set. These individual phase estimates are combined in a weighted average to form a new phase estimate, which serves as the initial phase estimate for either the next iteration of the iterative-transform method or, if the maximum number of iterations has been reached, for the next several steps, which constitute the outerloop portion of the algorithm. The details of the next several steps must be omitted here for the sake of brevity. The overall effect of these steps is to adaptively update the diversity defocus values according to recovery of global defocus in the phase estimate. Aberration recovery varies with differing amounts as the amount of diversity defocus is updated in each image; thus, feedback is incorporated into the recovery process. This process is iterated until the global defocus error is driven to zero during the recovery process. The amplitude of aberration may far exceed one wavelength after completion of the inner-loop portion of the algorithm, and the classical iterative transform method does not, by itself, enable recovery of multi-wavelength aberrations. Hence, in the absence of a means of off-loading the multi-wavelength portion of the aberration, the algorithm would produce a wrapped phase map. However, a special aberration-fitting procedure can be applied to the wrapped phase data to transfer at least some portion of the multi-wavelength aberration to the diversity function, wherein the data are treated as known phase values. In this way, a multiwavelength aberration can be recovered incrementally by successively applying the aberration-fitting procedure to intermediate wrapped phase maps. During recovery, as more of the aberration is transferred to the diversity function following successive iterations around the ter loop, the estimated phase ceases to wrap in places where the aberration values become incorporated as part of the diversity function. As a result, as the aberration content is transferred to the diversity function, the phase estimate resembles that of a reference flat.

Dean, Bruce H.↗

Constraining the masses of high-redshift clusters with weak lensing: Revised shape calibration testing for the impact of stronger shears and increased blending

Weak lensing measurements suffer from well-known shear estimation biases, which can be partially corrected for with the use of image simulations. Here we present an analysis of simulated images that mimic Hubble Space Telescope/Advance Camera for Surveys observations of high-redshift galaxy clusters, including cluster specific issues such as non-weak shear and increased blending. Our synthetic galaxies have been generated to have similar observed properties as the background-selected source samples studied in the real images. First, we used simulations with galaxies placed on a grid to determine a revised signal-to-noise-dependent (S/N KSB ) correction for multiplicative shear measurement bias, and to quantify the sensitivity of our KSB+ bias calibration to mismatches of galaxy or PSF properties between the real data and the simulations. Next, we studied the impact of increased blending and light contamination from cluster and foreground galaxies, finding it to be negligible for high-redshift (z > 0.7) clusters, whereas shear measurements can be affected at the ~1% level for lower redshift clusters given their brighter member galaxies. Finally, we studied the impact of fainter neighbours and selection bias using a set of simulated images that mimic the positions and magnitudes of galaxies in Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS) data, thereby including realistic clustering. While the initial SExtractor object detection causes a multiplicative shear selection bias of –0.028 ± 0.002, this is reduced to –0.016 ± 0.002 by further cuts applied in our pipeline. Given the limited depth of the CANDELS data, we compared our CANDELS-based estimate for the impact of faint neighbours on the multiplicative shear measurement bias to a grid-based analysis, to which we added clustered galaxies to even fainter magnitudes based on Hubble Ultra Deep Field data, yielding a refined estimate of ~ –0.013. Our sensitivity analysis suggests that our pipeline is calibrated to an accuracy of ~0.015 once all corrections are applied, which is fully sufficient for current and near-future weak lensing studies of high-redshift clusters. As an application, we used it for a refined analysis of three highly relaxed clusters from the South Pole Telescope Sunyaev-Zeldovich survey, where we now included measurements down to the cluster core (r > 200 kpc) as enabled by our work. Compared to previously employed scales (r > 500 kpc), this tightens the cluster mass constraints by a factor 1.38 on average.

79 ASTRONOMY AND ASTROPHYSICS↗

Computational investigation of incident ion angles and material erosion at rough graphite and silicon carbide divertor surfaces

We present a computational investigation of the dependence of material erosion on the incident ion angle at rough graphite and silicon carbide divertor surfaces. Ion angle distributions (IADs) for D plasmas at NSTX-U and DIII-D divertors were calculated by an equation-of-motion model that traces the ion trajectories in the sheath. Then, the effective sputtering yields and ion shadowed area fractions were calculated by a Monte Carlo micro-patterning and roughness code that applied the calculated IADs to surface topographic data that were obtained from optical confocal microscopy of rough graphite and SiC divertor surfaces from NSTX-U and DIII-D experiments. The calculations found that the effective sputtering yields, the sputtering pattern, and the shadowed area are determined by the detailed surface topology rather than the root mean square roughness R RMS , which represents deviations from a flat surface. The suppression of the effective sputtering yields for rough surfaces compared to the yield for a smooth surface was accounted for by the change of the mean local incident ion angle (LIIA) (θ'). The mean surface inclination angle distribution (SIAD) (δ) was found to be a useful parameter to estimate the LIIA from the calculated IADs. We report global empirical formulas for the mean LIIA and fraction of the area shadowed from the main ions for D plasmas for rough surfaces with B-field incident angles α = 85°–89° as a function of the mean SIAD (δ). We propose the use of the mean LIIA (θ') to estimate the sputtering yield for rough surfaces from the angular dependence of the sputtering yield.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Fluid volcanism on Miranda and Ariel - Flow morphology and composition

Several types of volcanic units have been recognized on the icy Uranian satellites Miranda and Ariel. On Miranda, ridges characterized by crest grooves are up to 10 km wide and 500 m high. A continuous flat-topped flow band also 10 km wide and 500 m high forms the outer southern margin of Elsinore Corona, which appears to comprise coalesced flow bands and ridges. On Ariel, in addition to at least one ridge unit similar to those on Miranda, flood plains material has covered the floors of deep chasmata (grabens) and an irregular depression. Flows on both satellites are characterized by linear vent geometries and distinct topographic margins, which indicate extrusion of a relatively viscous material. The topography of the flows can be used to estimate flow viscosity or yield strength using a Bingham plastic model. Extrusion viscosity estimates, incorporating plausible volcanologically based emplacement time scales and a rigid crust correction, range from 10 MP to 1 GP (10 TP in the unlikely absence of a chilled crust). Viscosity estimates are dependent on the assumed emplacement time scale, however, and could be as high as 10 PP, if a solid-state-based time scale is assumed.

Schenk, Paul M.↗

Spectrum Modal Analysis for the Detection of Low-Altitude Windshear with Airborne Doppler Radar

A major obstacle in the estimation of windspeed patterns associated with low-altitude windshear with an airborne pulsed Doppler radar system is the presence of strong levels of ground clutter which can strongly bias a windspeed estimate. Typical solutions attempt to remove the clutter energy from the return through clutter rejection filtering. Proposed is a method whereby both the weather and clutter modes present in a return spectrum can be identified to yield an unbiased estimate of the weather mode without the need for clutter rejection filtering. An attempt will be made to show that modeling through a second order extended Prony approach is sufficient for the identification of the weather mode. A pattern recognition approach to windspeed estimation from the identified modes is derived and applied to both simulated and actual flight data. Comparisons between windspeed estimates derived from modal analysis and the pulse-pair estimator are included as well as associated hazard factors. Also included is a computationally attractive method for estimating windspeeds directly from the coefficients of a second-order autoregressive model. Extensions and recommendations for further study are included.

Kunkel, Matthew W.↗

Evaluation of the Event Driven Phenology Model Coupled with the VegET Evapotranspiration Model Through Comparisons with Reference Datasets in a Spatially Explicit Manner

A new model coupling scheme with remote sensing data assimilation was developed for estimation of daily actual evapotranspiration (ET). The scheme represents a mix of the VegET, a physically based model to estimate ET from a water balance, and an event driven phenology model (EDPM), where the EDPM is an empirically derived crop specific model capable of producing seasonal trajectories of canopy attributes. In this experiment, the scheme was deployed in a spatially explicit manner within the croplands of the Northern Great Plains. The evaluation was carried out using 2007-2009 land surface forcing data from the North American Land Data Assimilation System (NLDAS) and crop maps derived from remotely sensed data of NASA's Moderate Resolution Imaging Spectroradiometer (MODIS). We compared the canopy parameters produced by the phenology model with normalized difference vegetation index (NDVI) data derived from the MODIS nadir bi-directional reflectance distribution function (BRDF) adjusted reflectance (NBAR) product. The expectations of the EDPM performance in prognostic mode were met, producing determination coefficient (r2) of 0.8 +/-.0.15. Model estimates of NDVI yielded root mean square error (RMSE) of 0.1 +/-.0.035 for the entire study area. Retrospective correction of canopy dynamics with MODIS NDVI brought the errors down to just below 10% of observed data range. The ET estimates produced by the coupled scheme were compared with ones from the MODIS land product suite. The expected r2=0.7 +/-.15 and RMSE = 11.2 +/-.4 mm per 8 days were met and even exceeded by the coupling scheme0 functioning in both prognostic and retrospective modes. Minor setbacks of the EDPM and VegET performance (r2 about 0.5 and additional 30 % of RMSR) were found on the peripheries of the study area and attributed to the insufficient EDPM training and to spatially varying accuracy of crop maps. Overall the experiment provided sufficient evidence of soundness and robustness of the EDPM and VegET coupling scheme, assuring its potential for spatially explicit applications.

Kovalskyy, V.↗

Likelihood Analysis for Mega Pixel Maps

The derivation of cosmological parameters from astrophysical data sets routinely involves operations counts which scale as O(N(exp 3) where N is the number of data points. Currently planned missions, including MAP and Planck, will generate sky maps with N(sub d) = 10(exp 6) or more pixels. Simple "brute force" analysis, applied to such mega-pixel data, would require years of computing even on the fastest computers. We describe an algorithm which allows estimation of the likelihood function in the direct pixel basis. The algorithm uses a conjugate gradient approach to evaluate X2 and a geometric approximation to evaluate the determinant. Monte Carlo simulations provide a correction to the determinant, yielding an unbiased estimate of the likelihood surface in an arbitrary region surrounding the likelihood peak. The algorithm requires O(N(sub d)(exp 3/2) operations and O(Nd) storage for each likelihood evaluation, and allows for significant parallel computation.

Kogut, Alan J.↗