Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Prior distribution”

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 253 records · Page 14

Data-driven prediction of scaling and ignition of inertial confinement fusion experiments

Recent advances in inertial confinement fusion (ICF) at the National Ignition Facility (NIF), including ignition and energy gain, are enabled by a close coupling between experiments and high-fidelity simulations. Neither simulations nor experiments can fully constrain the behavior of ICF implosions on their own, meaning pre- and postshot simulation studies must incorporate experimental data to be reliable. Linking past data with simulations to make predictions for upcoming designs and quantifying the uncertainty in those predictions has been an ongoing challenge in ICF research. We have developed a data-driven approach to prediction and uncertainty quantification that combines large ensembles of simulations with Bayesian inference and deep learning. The approach builds a predictive model for the statistical distribution of key performance parameters, which is jointly informed by past experiments and physics simulations. The prediction distribution captures the impact of experimental uncertainty, expert priors, design changes, and shot-to-shot variations. We have used this new capability to predict a 10× increase in ignition probability between Hybrid-E shots driven with 2.05 MJ compared to 1.9 MJ, and validated our predictions against subsequent experiments. We describe our new Bayesian postshot and prediction capabilities, discuss their application to NIF ignition and validate the results, and finally investigate the impact of data sparsity on our prediction results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A mixed, unified forward/inverse framework for earthquake problems: fault implementation and coseismic slip estimate

SUMMARY We introduce a new finite-element (FE) based computational framework to solve forward and inverse elastic deformation problems for earthquake faulting via the adjoint method. Based on two advanced computational libraries, FEniCS and hIPPYlib for the forward and inverse problems, respectively, this framework is flexible, transparent and easily extensible. We represent a fault discontinuity through a mixed FE elasticity formulation, which approximates the stress with higher order accuracy and exposes the prescribed slip explicitly in the variational form without using conventional split node and decomposition discrete approaches. This also allows the first order optimality condition, that is the vanishing of the gradient, to be expressed in continuous form, which leads to consistent discretizations of all field variables, including the slip. We show comparisons with the standard, pure displacement formulation and a model containing an in-plane mode II crack, whose slip is prescribed via the split node technique. We demonstrate the potential of this new computational framework by performing a linear coseismic slip inversion through adjoint-based optimization methods, without requiring computation of elastic Green’s functions. Specifically, we consider a penalized least squares formulation, which in a Bayesian setting—under the assumption of Gaussian noise and prior—reflects the negative log of the posterior distribution. The comparison of the inversion results with a standard, linear inverse theory approach based on Okada’s solutions shows analogous results. Preliminary uncertainties are estimated via eigenvalue analysis of the Hessian of the penalized least squares objective function. Our implementation is fully open-source and Jupyter notebooks to reproduce our results are provided. The extension to a fully Bayesian framework for detailed uncertainty quantification and non-linear inversions, including for heterogeneous media earthquake problems, will be analysed in a forthcoming paper.

58 GEOSCIENCES↗

Robust error calibration for serial crystallography

Serial crystallography is an important technique with unique abilities to resolve enzymatic transition states, minimize radiation damage to sensitive metalloenzymes and perform de novo structure determination from micrometre-sized crystals. This technique requires the merging of data from thousands of crystals, making manual identification of errant crystals unfeasible. cctbx.xfel.merge uses filtering to remove problematic data. However, this process is imperfect, and data reduction must be robust to outliers. We add robustness to cctbx.xfel.merge at the step of uncertainty determination for reflection intensities. This step is a critical point for robustness because it is the first step where the data sets are considered as a whole, as opposed to individual lattices. Robustness is conferred by reformulating the error-calibration procedure to have fewer and less stringent statistical assumptions and incorporating the ability to down-weight low-quality lattices. We then apply this method to five macromolecular XFEL data sets and observe the improvements to each. The appropriateness of the intensity uncertainties is demonstrated through internal consistency. This is performed through theoretical CC 1/2 and I /σ relationships and by weighted second moments, which use Wilson's prior to connect intensity uncertainties with their expected distribution. This work presents new mathematical tools to analyze intensity statistics and demonstrates their effectiveness through the often underappreciated process of uncertainty analysis.

Mittan-Moreau, David W.↗

A derecho climatology (2004–2021) in the United States based on machine learning identification of bow echoes

Due to their persistent widespread severe winds, derechos pose significant threats to human safety and property, with impacts comparable to many tornadoes and hurricanes. Yet, automated detection of derechos remains challenging due to the absence of spatiotemporally continuous observations and the complex criteria employed to define the phenomenon. This study presents an objective derecho detection approach capable of automatically identifying derechos through both observations and model results. The approach is grounded in a physically based definition of derechos and integrates three algorithms: (1) the Python Flexible Object Tracker (PyFLEXTRKR) algorithm to track mesoscale convective systems (MCSs), (2) a semantic segmentation convolutional neural network to identify bow echoes, and (3) a comprehensive classification algorithm to detect derechos within MCS life cycles and distinguish derecho-producing from non-derecho-producing MCSs. Using this approach, we developed a novel high-resolution (4 km and hourly) observational dataset of derechos and accompanying derecho-producing MCSs over the United States east of the Rocky Mountains from 2004 to 2021. The dataset consists of two subsets based on different gust speed data sources and is analyzed to document the climatology of derechos in the United States. On average, 12–15 derechos are identified per year, aligning with previous estimations (∼6–21 events annually). The spatial distribution and seasonal variation patterns are consistent with prior studies, showing peak occurrences in the Great Plains and the Midwest during the warm season. Additionally, during the study period, derechos account for approximately 3.1 % of measured damaging gusts (≥25.93 m s−1) over the eastern United States. The dataset is publicly available at https://doi.org/10.5281/zenodo.14835362 (Li et al., 2025).

54 ENVIRONMENTAL SCIENCES↗

Role of energetic particles in charging/discharging of spacecraft dielectrics

The role that energetic particles in the substorm plasma have on the charging and discharging of typical dielectric layers used on spacecraft was investigated using spectra and pitch angle distributions measured in situ on the SCATHA spacecraft prior to and during a few kilovolt differential charging event in eclipse conditions. The particle spectra was input to deposition codes that determine the dose rate as a function of depth in Kapton and Teflon layers used in an experiment on SCATHA. The calculated ambient dose rates of a few rads/sec throughout the bulk of the samples are sufficiently high that radiation damage levels are reached on the time scale of 1 year. Surface dose is a factor of 100 higher. Bulk conductivity profiles were obtained from the dose rates using empirical relationships. The radiation induced bulk conductivities calculated at the peak charging time are found to be smaller than the intrinsic dark conductivity range of solar conditioned Kapton but higher than the corresponding value for Teflon.

Reagan, J. B.↗

A Search for Vector Magnetic Field Variations Associated with the M-Class Flares of 1991 June 10 in AR 6659

A careful analysis of a 6-hour time sequence of vector magnetograms of AR 6659, observed on 1991 June 10 with the MSFC vector magnetograph, has revealed only minor changes in the vector magnetic field azimuths in the vicinity of two M-class flares, and the association of these changes with the flares is not unambiguous. In this paper we present our analysis of the data which includes comparison of vector magnetograms prior to and during the flares, calculation of distributions of the rms variation of the azimuth at each pixel in the field of view of the active region, and examination of the variation with time of the azimuths at every pixel covered by the main flare emissions as observed with the H-alpha telescope coaligned with the vector magnetograph. We also present results of an analysis of evolutionary changes in the azimuth over the field of view of the active region.

Hagyard, Mona J.↗

A Search for Vector Magnetic Field Variations Associated with the M-Class Flares of 1991 June 10 in AR 6659

A careful analysis of a 6-hour time sequence of vector magnetograms of AR 6659, observed on 1991 June 10 with the MSFC vector magnetograph, has revealed only minor changes in the vector magnetic field azimuths in the vicinity of two M-class flares, and the association of these changes with the flares is not unambiguous. In this paper we present our analysis of the data which includes comparison of vector magnetograms prior to and during the flares, calculation of distributions of the rms variation of the azimuth at each pixel in the field of view of the active region, and examination of the variation with time of the azimuths at every pixel covered by the main flare emissions as observed with the H-alpha telescope coaligned with the vector magnetograph. We also present results of an analysis of evolutionary changes in the azimuth over the field of view of the active region.

Hagyard, Mona J.↗

A Physical Model to Estimate Snowfall over Land using AMSU-B Observations

In this study, we present an improved physical model to retrieve snowfall rate over land using brightness temperature observations from the National Oceanic and Atmospheric Administration's (NOAA) Advanced Microwave Sounder Unit-B (AMSU-B) at 89 GHz, 150 GHz, 183.3 +/- 1 GHz, 183.3 +/- 3 GHz, and 183.3 +/- 7 GHz. The retrieval model is applied to the New England blizzard of March 5, 2001 which deposited about 75 cm of snow over much of Vermont, New Hampshire, and northern New York. In this improved physical model, prior retrieval assumptions about snowflake shape, particle size distributions, environmental conditions, and optimization methodology have been updated. Here, single scattering parameters for snow particles are calculated with the Discrete-Dipole Approximation (DDA) method instead of assuming spherical shapes. Five different snow particle models (hexagonal columns, hexagonal plates, and three different kinds of aggregates) are considered. Snow particle size distributions are assumed to vary with air temperature and to follow aircraft measurements described by previous studies. Brightness temperatures at AMSU-B frequencies for the New England blizzard are calculated using these DDA calculated single scattering parameters and particle size distributions. The vertical profiles of pressure, temperature, relative humidity and hydrometeors are provided by MM5 model simulations. These profiles are treated as the a priori data base in the Bayesian retrieval algorithm. In algorithm applications to the blizzard data, calculated brightness temperatures associated with selected database profiles agree with AMSU-B observations to within about +/- 5 K at all five frequencies. Retrieved snowfall rates compare favorably with the near-concurrent National Weather Service (NWS) radar reflectivity measurements. The relationships between the NWS radar measured reflectivities Z(sub e) and retrieved snowfall rate R for a given snow particle model are derived by a histogram matching technique. All of these Z(sub e)-R relationships fall in the range of previously established Z(sub e)-R relationships for snowfall. This suggests that the current physical model developed in this study can reliably estimate the snowfall rate over land using the AMSU-B measured brightness temperatures.

Kim, Min-Jeong↗

Likelihood-Based Climate Model Evaluation

Climate models are deterministic, mathematical descriptions of the physics of climate. Confidence in predictions of future climate is increased if the physics are verifiably correct. A necessary, (but not sufficient) condition is that past and present climate be simulated well. Quantify the likelihood that a (summary statistic computed from a) set of observations arises from a physical system with the characteristics captured by a model generated time series. Given a prior on models, we can go further: posterior distribution of model given observations.

climate models↗

System and process for purification of astatine-211 from target materials

A new column-based purification system and approach are described for rapid separation and purification of the alpha-emitting therapeutic radioisotope 211At from dissolved cyclotron targets that provide highly reproducible product results with excellent 211At species distributions and high antibody labeling yields compared with prior art manual extraction results of the prior art that can be expected to enable enhanced production of purified 211At isotope products suitable for therapeutic medical applications such as treatment of cancer in human patients.

O'Hara, Matthew J.↗

System and process for purification of Astatine-211 from target materials

A new column-based purification system and approach are described for rapid separation and purification of the alpha-emitting therapeutic radioisotope 211 At from dissolved cyclotron targets that provide highly reproducible product results with excellent 211 At species distributions and high antibody labeling yields compared with prior art manual extraction results of the prior art that can be expected to enable enhanced production of purified 211 At isotope products suitable for therapeutic medical applications such as treatment of cancer in human patients.

O'Hara, Matthew J.↗

A Large Fraction of Hydrogen-rich Supernova Progenitors Experience Elevated Mass Loss Shortly Prior to Explosion

Spectroscopic detection of narrow emission lines traces the presence of circumstellar mass distributions around massive stars exploding as core-collapse supernovae. Transient emission lines disappearing shortly after the supernova explosion suggest that the material spatial extent is compact and implies an increased mass loss shortly prior to explosion. Here, we present a systematic survey for such transient emission lines (Flash Spectroscopy) among Type II supernovae detected in the first year of the Zwicky Transient Facility survey. We find that at least six out of ten events for which a spectrum was obtained within two days of the estimated explosion time show evidence for such transient flash lines. Our measured flash event fraction (>30% at 95% confidence level) indicates that elevated mass loss is a common process occurring in massive stars that are about to explode as supernovae.

Rachel J. Bruch↗

Effect of sample size on the maximum value distribution of fatigue driving forces in metals and alloys

An analytical framework is presented to predict the effects of sample size on the maximum value distribution (MVD) of the driving forces for fatigue crack formation in metals and alloys. The distribution of the maximum driving force for fatigue crack formation over the domain follows the generalized extreme value theory in the limit as the domain size increases to infinity. Here, a simulation-based analysis of microstructure influences on fatigue resistance for polycrystalline metals and alloys is very costly, and reaching those limits is intractable. This work models the MVD of Fatigue Indicator Parameters (FIPs), which serve as surrogate measures for the driving force for fatigue crack formation, at finite sample sizes prior to their convergence to a limiting extreme value distribution. Large-scale crystal plasticity finite element (CPFE) simulations of FCC Al 7075-T6 with microstructure realizations of various sizes are incorporated to calibrate and evaluate the developed framework, and a total of ∼6.5 million grains of Al 7075-T6 are examined. The calibrated analytical solution agrees well with the brute force Monte Carlo simulation results extracted from the CPFE simulations. Furthermore, the developed formulation can predict the MVD of FIPs for different sample sizes using a size-dependent parameter, and it is capable of accurately extrapolating the MVD of FIPs for much larger microstructure sample sizes than the size used for its calibration.

Crystal plasticity↗

Experimental decoy-state Bennett-Brassard 1984 quantum key distribution through a turbulent channel

In free-space quantum key distribution (QKD) in turbulent conditions, scattering and beam wandering cause intensity fluctuations which decrease the detected signal-to-noise ratio. This effect can be mitigated by rejecting received bits when the channel's transmittance is below a threshold. Thus, the overall error rate is reduced and the secure key rate increases despite the deletion of bits. In this work, we implement recently proposed selection methods focusing on the prefixed-threshold real-time selection (P-RTS) where a cutoff can be chosen prior to data collection and independently of the transmittance distribution. We perform finite-size decoy-state Bennett-Brassard 1984 QKD in a laboratory setting where we simulate the atmospheric turbulence using an acousto-optical modulator. We show that P-RTS can yield considerably higher secure key rates for a wide range of the atmospheric channel parameters. In addition, we evaluate the performance of the P-RTS method for a realistically finite sample size. We demonstrate that a near-optimal selection threshold can be predetermined even with imperfect knowledge of the channel transmittance distribution parameters.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Overwintering Distribution and Postspawn Survival of Steelhead in the Upper Columbia River Basin

Abstract Adult summer-run steelhead Oncorhynchus mykiss overwinter in freshwater for several months prior to spawning. In systems supporting mixed populations of fisheries and conservation importance, understanding the distribution and survival of pre- and postspawn fish is necessary for informed management. The upper Columbia River supports hatchery-origin components valued by anglers, natural-origin components of conservation concern, and temporary strays from downstream populations. We used radiotelemetry and PIT telemetry to monitor the behavior, distribution, and survival of adult steelhead during fall, overwintering, spawning, and postspawn periods, with a focus on use of the main stem versus four major tributaries. Adult steelhead (N = 807) were tagged at Priest Rapids Dam in 2015 and 2016. One-fifth of steelhead fell back below Priest Rapids Dam and did not reascend. A slight majority of tagged steelhead that overwintered upstream of Priest Rapids Dam did so in main-stem reservoirs (54%; N = 548). Overwintering in the main-stem Columbia River was more likely for later-arriving steelhead and was concentrated in the upstream-most reservoir. Winter tributary use was highest in the Wenatchee (26%; 2016) and Methow (18%; 2015) rivers, whereas no steelhead overwintered in the Entiat River. Harvest of hatchery-origin steelhead was 18% in 2015 and was near zero in 2016, when the fishery was suspended due to low adult returns. After accounting for reported harvest, annual overwinter survival did not differ between main-stem and tributary habitats, and relatively low adjusted survival of hatchery-origin steelhead in 2015 suggested unreported harvest. In contrast to low iteroparity rates (<3%), the majority of postspawn steelhead (56.5%) exited tributaries as kelts; kelt survival to Bonneville Dam was 65% in 2016 and 23% in 2017. Collectively, the results highlight the importance of understanding patterns of habitat use and mortality in steelhead populations when managers are faced with balancing harvest and conservation goals.

Fuchs, Nathaniel T.↗

Regularized Differentiation for Bioburden Density Estimation in Planetary Protection

In this paper, we propose and investigate the performance of two novel shrinkage estimators for bioburden density estimation in planetary protection. The estimators are based on the regularized differentiation of a cumulative count of colony forming units collected throughout the data collecting session or the life cycle of the entire mission. The regularized differentiation recasts the problem of bioburden density estimation as a linear least squares problem. The least squares problem is then solved through regularization techniques, such as truncated singular value decomposition and penalized least squares. The regularization is necessary to avoid noise amplification during the differentiation of noisy data. The two regularization estimators are compared with four other commonly used estimators to simultaneously evaluate the means of multivariable independent Poisson distributions: the maximum likelihood, noninformative Bayes estimator with Jeffreys prior, Empirical Bayes using conjugate gamma-Poisson model with gamma parameters selected by method of moments, and the Clevenson-Zidek estimator. It is shown through computer-simulated data that the regularized differentiation based on ridge regression has the smallest mean-squared error among all estimators. The analysis of shrinkage mechanism implemented by regularized differentiation is performed, and it is shown that the regularized differentiation amounts to performing a weighted averaging of all the samples. The weights are determined by the regularization parameter automatically selected by the L-curve technique. Since the method of least squares makes no distributional assumptions about the data, it presents an attractive technique for bioburden density estimation when there are concerns about the misspecification of the distributional model. The paper concludes with the analysis of the bioburden data collected during InSight mission and directions for future work.

97 - MATHEMATICS AND COMPUTING↗

Economic Evaluation of Modernization Expenditures for Electric Utility Distribution Systems: A Guide for Utility Regulators

Cost-effectiveness evaluation of potential grid modernization solutions is integral to integrated distribution system planning (IDSP). This report synthesizes and updates prior cost-effectiveness approaches identified by the U.S. Department of Energy for grid modernization investments, incorporating multi-objective decision-making. The structure of the report follows the IDSP process, from identifying objectives and priorities to prioritizing grid solutions. The report details two methods for conducting initial cost-effectiveness screening of individual and interdependent grid components: (1) Lowest Reasonable Cost and (2) Benefit-Cost Analysis. Drawing on solutions that passed the initial cost-effectiveness screen, the utility prioritizes grid solutions based on policy objectives and priorities, regulatory compliance, operational efficiencies, and other factors to develop a portfolio of distribution solutions. The utility uses multi-objective decision-making to assess each proposed expenditure against each objective and respective metric. The utility applies a weighting factor, reflecting the priority ranking of the objective, to the total numerical “score” based on the contribution of the proposed solution to addressing each objective. Ultimately, the utility ranks each expenditure from highest to lowest priority using the final score for each solution as well as its cost. The report includes examples of emerging best practices for states and utilities and a cost-effectiveness process checklist.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Cloud Modeling Using Field Project Data for the Study of Precipitation Processes

The use of cloud-resolving models (CRMs) in the study of precipitation process and their relation to the large-scale environment can be generally categorized into two approaches. The first approach is so called "cloud ensemble modeling". In this approach, many clouds of different size in various stages of their lifecycles can be present at any model simulation time. Large-scale effects are derived from observations and imposed into the model as the main forcing. The advantage of this approach is that the modeled convection will be forced to have the same intensity, thermodynamic budget and organization as the obserations.This approach will also allow CRMs to perform multi-day or multi-week time integrations. The second approach usually requires initial temperature and water vapor profiles that have a medium to large CAPE, and open lateral boundary conditions are used. The modeled clouds could be termed "self-forced convection". Model improvements, such as in the microphysics, are achieved using the second approach. In cloud ensemble modeling, accurate large-scale advective tendencies for temperature and water vapor are the main forcing for the CRMs. We found that the large-scale advective terms for temperature and water vapor are not always consistent, For example, large-scale forcing could indicate strong drying which would produce cooling in the model through evaporation but not contain large-scale advective heating to compensate. This discrepancy in forcing would cause differences between the observed and modeled latent heating profiles. Good measurements of other quantities (i.e., surface fluxes and radiation) are also required to perform variational objective analysis that computes and minimizes a "cost function" that constrains the difference between the large-scale advective forcing in temperature and water vapor. With self-forced convection, accurate vertical distributions of temperature, moisture (water vapor), and horizontal winds are required. The timing of the measurements relative to cloud development is crucial (i.e., prior to cloud triggering). Microphysical measurements (i.e., the cloud number concentration and size distribution) can also be used in this second approach but are of secondary importance with cloud ensemble modeling. In this paper, data collected during TRMM field campaigns (FCs; i.e., SCSMEX, LBA and KWAJEX) which were aimed at validating TRMM products (i.e., rainfall and the vertical distribution of latent heating) will be used to examine the impact of errors in the initial conditions (e.g.., soundings an large-scale forcing) on simulated rainfall distributions and brightness. Rainfall and precipitations simulated from a CRM will also be compared with those estimated by a schocastic model.

Tao, W.-K.↗