Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “robust 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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Mapping Permafrost Variability and Degradation Using Seismic Surface Waves, Electrical Resistivity, and Temperature Sensing: A Case Study in Arctic Alaska

Subsurface processes significantly influence surface dynamics in permafrost regions, necessitating utilizing diverse geophysical methods to reliably constrain permafrost characteristics. This research uses multiple geophysical techniques to explore the spatial variability of permafrost in undisturbed tundra and its degradation in disturbed tundra in Utqiagvik, Alaska. Here, we integrate multiple quantitative techniques, including multichannel analysis of surface waves (MASW), electrical resistivity tomography (ERT), and ground temperature sensing, to study heterogeneity in permafrost’s geophysical characteristics. MASW results reveal active layer shear wave velocities (V s ) between 240 and 370 m/s, and permafrost V s between 450 and 1,700 m/s, typically showing a low-high-low velocity pattern. Additionally, we find an inverse relationship between in situ V s and ground temperature measurements. The V s profiles along with electrical resistivity profiles reveal cryostructures such as cryopeg and ice-rich zones in the permafrost layer. The integrated results of MASW and ERT provide valuable information for characterizing permafrost heterogeneity and cryostructure. Corroboration of these geophysical observations with permafrost core samples’ stratigraphies and salinity measurements further validates these findings. This combination of geophysical and temperature sensing methods along with permafrost core sampling confirms a robust approach for assessing permafrost’s spatial variability in coastal environments. Our results also indicate that civil infrastructure systems such as gravel roads and pile foundations affect permafrost by thickening the active layer, lowering the V s , and reducing heterogeneity. We show how the resulting V s profiles can be used to estimate key parameters for designing buildings in permafrost regions and maintaining existing infrastructure in polar regions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Pairing properties of the $t$–$t$'–$t$"–$J$ model

We study the pairing properties of the two-dimensional t–t'–t"–J model, where t' and t" are second and third neighbor hoppings, at a doping level x ≈ 0.1. Recent studies of the t–t'–J model find strong pairing for t' > 0, associated with electron doping, but an absence of pairing for t'<0, associated with hole doping. This is in contrast to the cuprates, where the highest transition temperatures appear for hole doping. Model parametrizations for the cuprates estimate a t" comparable to t', which, in principle, might fix this discrepancy. However, we find that it does not; we observe a suppression of pairing for the hole-doped system (t' < 0,t" > 0), while for the electron-doped system (t' > 0,t" < 0) d-wave pairing is robust. As a result, extended hoppings appear to be insufficient to make the one-band t–t'–J model capable of describing the pairing in the hole-doped system.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Sensitive detection of structural dynamics using a statistical framework for comparative crystallography

Chemical and conformational changes are crucial to protein function and its pharmacological control. X-ray crystallography can reveal these changes in atomic detail, but standard analysis methods, which refine separate datasets, often overlook differences that are subtle or arise in only a subset of molecules. Direct comparison of crystallographic datasets is, in principle, more powerful, but systematic errors (“scales”) often mask changes in the crystallographic observables (“structure factors”). Machine learning algorithms that jointly estimate scales and structure factors can address this limitation. Here, we augment this approach with multivariate, structured priors derived from crystallographic theory, implemented in the variational deep learning framework Careless. Doing so strongly improves the detection of protein dynamics, element-specific anomalous signals, and the binding of drug candidates, offering a robust approach to comparative crystallography and, potentially, to detection of protein dynamics by other structure determination methods.

Hekstra, Doeke R. [Harvard Univ., Cambridge, MA (U↗
DOE OSTI2021

Characterizing the magnetic fields of nearby molecular clouds using submillimeter polarization observations

Of all the factors that influence star formation, magnetic fields are perhaps the least well understood. The goal of this paper is to characterize the 3D magnetic field properties of nearby molecular clouds through various methods of statistically analysing maps of polarized dust emission. Our study focuses on nine clouds, with data taken from the Planck Sky Survey as well as data from the Balloon-borne Large Aperture Submillimeter Telescope for Polarimetry observations of Vela C. We compare the distributions of polarization fraction (p), dispersion in polarization angles ($\mathcal {S}$), and hydrogen column density (N H ) for each of our targeted clouds. To broaden the scope of our analysis, we compare the distributions of our clouds’ polarization observables with measurements from synthetic polarization maps generated from numerical simulations. We also use the distribution of polarization fraction measurements to estimate the inclination angle of each cloud’s cloud-scale magnetic field. We obtain a range of inclination angles associated with our clouds, varying from 16° to 69°. We establish inverse correlations between p and both $\mathcal {S}$ and NH in almost every cloud, but we are unable to establish a statistically robust $\mathcal {S}$ versus NH trend. By comparing the results of these different statistical analysis techniques, we are able to propose a more comprehensive view of each cloud’s 3D magnetic field properties. These detailed cloud analyses will be useful in the continued studies of cloud-scale magnetic fields and the ways in which they affect star formation within these molecular clouds.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Sparsifying priors for Bayesian uncertainty quantification in model discovery

We propose a probabilistic model discovery method for identifying ordinary differential equations governing the dynamics of observed multivariate data. Our method is based on the sparse identification of nonlinear dynamics (SINDy) framework, where models are expressed as sparse linear combinations of pre-specified candidate functions. Promoting parsimony through sparsity leads to interpretable models that generalize to unknown data. Instead of targeting point estimates of the SINDy coefficients, we estimate these coefficients via sparse Bayesian inference. The resulting method, uncertainty quantification SINDy (UQ-SINDy), quantifies not only the uncertainty in the values of the SINDy coefficients due to observation errors and limited data, but also the probability of inclusion of each candidate function in the linear combination. UQ-SINDy promotes robustness against observation noise and limited data, interpretability (in terms of model selection and inclusion probabilities) and generalization capacity for out-of-sample forecast. Sparse inference for UQ-SINDy employs Markov chain Monte Carlo, and we explore two sparsifying priors: the spike and slab prior, and the regularized horseshoe prior. UQ-SINDy is shown to discover accurate models in the presence of noise and with orders-of-magnitude less data than current model discovery methods, thus providing a transformative method for real-world applications which have limited data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Reweighting Monte Carlo predictions and automated fragmentation variations in Pythia 8

This work reports on a method for uncertainty estimation in simulated collider-event predictions. The method is based on a Monte Carlo-veto algorithm, and extends previous work on uncertainty estimates in parton showers by including uncertainty estimates for the Lund string-fragmentation model. This method is advantageous from the perspective of simulation costs: a single ensemble of generated events can be reinterpreted as though it was obtained using a different set of input parameters, where each event now is accompanied with a corresponding weight. This allows for a robust exploration of the uncertainties arising from the choice of input model parameters, without the need to rerun full simulation pipelines for each input parameter choice. Such explorations are important when determining the sensitivities of precision physics measurements. Accompanying code is available at https://gitlab.com/uchep/mlhad-weights-validation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Heliostat Sizing Methodology for Solar Heat for Industrial Processes

This study presents a method to obtain a heliostat size that minimizes the levelized cost of a heliostat-based concentrating solar thermal system for industrial process heat (IPH) applications at operating temperatures from 565 to 1550 degrees Celsius. The method extends prior work by embedding a routine for system design that obtains near-optimal subsystem sizes, increasing the fidelity of drive cost functions, and adding an optical performance model. An illustrative business case is developed for Daggett, California, targeting specified annual thermal energy outputs of 50 to 400 GWhth. Optical performance is modeled using verified estimates from the literature. A surrogate heliostat cost model, derived from commercial heliostat designs and scaled for production volume, installation, and operations and maintenance costs, is used to develop cost functions. Results show that heliostat size strongly affects the levelized cost of heat (LCOH), producing a characteristic U-shaped trend with a robust near-optimal window of 8 - 12 m2; the heliostat size producing the lowest project cost in our study grows slightly as the project size increases, and is reduced as the operating temperature increases. The findings in this study are consistent with the general trend of smaller heliostats under deployment at existing projects for high-temperature industrial process heat and reflect the significant reduction in power electronics and other per-heliostat costs. The methodology we propose is general and can be tailored to revised cost curves as the technology continues to evolve.

14 SOLAR ENERGY↗
DOE OSTI2021

Heavy Oil Polymer EOR in the Challenging Alaskan Arctic - It Works!

Under co-sponsorship of the US Department of Energy and Hilcorp Alaska LLC the first ever polymer field pilot commenced on 8/28/2018 in the Schrader Bluff heavy oil reservoir at the Milne Point Field on Alaska North Slope (ANS). The primary objective of the pilot is to prove the efficacy of polymer Enhanced Oil Recovery (EOR) to unlock the vast heavy oil resources on ANS. More than two and half years after startup, the polymer injection, supporting laboratory experiments and simulation studies steadily continue. The pilot started injecting hydrolyzed polyacrylamide (HPAM), at a concentration of 1,750 ppm to achieve a target viscosity of 45 cP, into the two horizontal injectors in the flood pattern. Production is monitored in the two horizontal producers. Based on laboratory measurements of polymer viscosity at reservoir conditions, the team decided to reduce polymer concentration to 1,200 ppm since July 2020 in an effort to control injection pressure and optimize polymer utilization. Quality control (QC) on the field ensures uniform polymer solution properties. Representative rock and fluid systems and test conditions are utilized in the corefloods on polymer retention, effect of injection water salinity, polymer loading, and their combinations on oil recovery. A history matched reservoir simulation model for forecasting oil recovery was developed on the basis of all the available field data. Field concerns related to the post-polymer breakthrough impact on flow assurance is addressed via specialized laboratory tests. Notwithstanding early operational disruptions and hydration issues, continuous polymer injection in both injectors has been achieved. To date, 950,000 lbs of polymer or 2 million barrels of polymer solution, equating to 13% of total pore volume (PV), has been placed in the flood pattern, serving as an effective indicator of adequate polymer injectivity. So far, the success of polymer EOR is evident from drastically reduced water cut in the producers, an estimated incremental 1,000 bopd over waterflood, and a favorable polymer utilization of 1.7 lbs/barrel of incremental oil. Polymer breakthrough was observed 26 months after the start of polymer injection. Main observations from corefloods are unusually high polymer retention values in some cases and a positive response to low salinity water. Although the heterogeneity in the flood pattern and exceptionally low water cut pose some challenges, persistent novel and justifiable simulation approaches have resulted in a robust history matched model. Experimental results on produced fluids treatment provide operational guidance for improved oil-water separation and mitigation of heater tube fouling. The ongoing success of the pilot is a key indicator of bringing the team ever closer to meeting the project’s primary objective. The encouraging results of the pilot is one of the drivers that has provided the impetus to apply polymer EOR throughout the Milne Point Field, which would increase oil recovery and extend the economic life of the Trans Alaska Pipeline System. The scientific knowledge, including the many lessons learned from this pilot also has referential value for other potential heavy oil EOR projects throughout the world. The collected data, operational lessons learned, and the overall success of the pilot are summarized in the paper.

None, None↗
DOE OSTI2021

Forecasting Multi-Wave Epidemics Through Bayesian Inference

We present a simple, near-real-time Bayesian method to infer and forecast a multiwave outbreak, and demonstrate it on the COVID-19 pandemic. The approach uses timely epidemiological data that has been widely available for COVID-19. It provides short-term forecasts of the outbreak’s evolution, which can then be used for medical resource planning. The method postulates one- and multiwave infection models, which are convolved with the incubation-period distribution to yield competing disease models. The disease models’ parameters are estimated via Markov chain Monte Carlo sampling and information-theoretic criteria are used to select between them for use in forecasting. The method is demonstrated on two- and three-wave COVID-19 outbreaks in California, New Mexico and Florida, as observed during Summer-Winter 2020. We find that the method is robust to noise, provides useful forecasts (along with uncertainty bounds) and that it reliably detected when the initial single-wave COVID-19 outbreaks transformed into successive surges as containment efforts in these states failed by the end of Spring 2020.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Data-Efficient Methods for Determining Flory–Huggins χ Parameters in Multicomponent Polymer Formulations

Polymer formulations are essential in diverse applications including personal care products, coatings, paints, adhesives, and plastic materials. Designing these formulations requires navigating large, complex design spaces, where phase and self-assembly behavior critically impact performance. The Flory–Huggins χ parameter, which quantifies segmental miscibility, is widely used to parametrize the excess free energy of mixing in formulation models. In this work, we introduce two data-efficient, top-down methods for estimating χ parameters using the Random Phase Approximation (RPA): (i) Boundary Nonlinear Regression (Boundary-NLR), which fits theoretical spinodal boundaries to experimental phase boundaries, and (ii) Surrogate Model Inverse Parameter Estimation (SMIPE), which uses a Gaussian Process Classifier to fit sparse phase maps via a surrogate model. Both methods allow rapid parametrization of polymer field-theoretic models without the need for additional experiments. We evaluate these approaches on data sets involving polymer–solvent–nonsolvent ternary mixtures and block copolymer–solvent systems, demonstrating their robustness to experimental noise and their relevance for real-world formulation design.

copolymers↗
DOE OSTI2023

An optimized search for dark matter in the galactic halo with HAWC

The Galactic Halo is a key target for indirect dark matter detection. The High Altitude Water Cherenkov (HAWC) observatory is a high-energy (~300 GeV to >100 TeV) gamma-ray detector located in central Mexico. HAWC operates via the water Cherenkov technique and has both a wide field of view of ~ 2 sr and a >95% duty cycle, making it ideal for analyses of highly extended sources. We made use of these properties of HAWC and a new background-estimation technique optimized for extended sources to probe a large region of the Galactic Halo for dark matter signals. With this approach, we set improved constraints on dark matter annihilation and decay between masses of 10 and 100 TeV. Due to the large spatial extent of the HAWC field of view, these constraints are robust against uncertainties in the Galactic dark matter spatial profile.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Using an Atom Interferometer to Infer Gravitational Entanglement Generation

If gravitational perturbations are quantized into gravitons in analogy with the electromagnetic field and photons, the resulting graviton interactions should lead to an entangling interaction between massive objects. We suggest a test of this prediction. To do this, we introduce the concept of interactive quantum information sensing. This novel sensing protocol is tailored to provable verification of weak dynamical entanglement generation between a pair of systems. We show that this protocol is highly robust to typical thermal noise sources. The sensitivity can moreover be increased both using an initial thermal state and/or an initial phase of entangling via a non-gravitational interaction. We outline a concrete implementation testing the ability of the gravitational field to generate entanglement between an atomic interferometer and mechanical oscillator. Preliminary numerical estimates suggest that near-term devices could feasibly be used to perform the experiment.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Universal dwell time optimization for deterministic optics fabrication

Computer-Controlled Optical Surfacing (CCOS) has been greatly developed and widely used for precision optical fabrication in the past three decades. It relies on robust dwell time solutions to determine how long the polishing tools must dwell at certain points over the surfaces to achieve the expected forms. However, as dwell time calculations are modeled as ill-posed deconvolution, it is always non-trivial to reach a reliable solution that 1) is non-negative, since CCOS systems are not capable of adding materials, 2) minimizes the residual in the clear aperture 3) minimizes the total dwell time to guarantee the stability and efficiency of CCOS processes, 4) can be flexibly adapted to different tool paths, 5) the parameter tuning of the algorithm is simple, and 6) the computational cost is reasonable. In this study, we propose a novel Universal Dwell time Optimization (UDO) model that universally satisfies these criteria. First, the matrix-based discretization of the convolutional polishing model is employed so that dwell time can be flexibly calculated for arbitrary dwell points. Second, UDO simplifies the inverse deconvolution as a forward scalar optimization for the first time, which drastically increases the solution stability and the computational efficiency. Finally, the dwell time solution is improved by a robust iterative refinement and a total dwell time reduction scheme. The superiority and general applicability of the proposed algorithm are verified on the simulations of different CCOS processes. A real application of UDO in improving a synchrotron X-ray mirror using Ion Beam Figuring (IBF) is then demonstrated. The simulation indicates that the estimated residual in the 92.3 mm × 15.7 mm CA can be reduced from 6.32 nm Root Mean Square (RMS) to 0.20 nm RMS in 3.37 min. After one IBF process, the measured residual in the CA converges to 0.19 nm RMS, which coincides with the simulation.

36 MATERIALS SCIENCE↗
DOE OSTI2024

In-vivo neuronal dysfunction by Aβ and tau overlaps with brain-wide inflammatory mechanisms in Alzheimer’s disease

The molecular mechanisms underlying neuronal dysfunction in Alzheimer’s disease (AD) remain uncharacterized. Here, we identify genes, molecular pathways and cellular components associated with whole-brain dysregulation caused by amyloid-beta (Aβ) and tau deposits in the living human brain. We obtained in-vivo resting-state functional MRI (rs-fMRI), Aβ- and tau-PET for 47 cognitively unimpaired and 16 AD participants from the Translational Biomarkers in Aging and Dementia cohort. Adverse neuronal activity impacts by Aβ and tau were quantified with personalized dynamical models by fitting pathology-mediated computational signals to the participant’s real rs-fMRIs. Then, we detected robust brain-wide associations between the spatial profiles of Aβ-tau impacts and gene expression in the neurotypical transcriptome (Allen Human Brain Atlas). Within the obtained distinctive signature of in-vivo neuronal dysfunction, several genes have prominent roles in microglial activation and in interactions with Aβ and tau. Moreover, cellular vulnerability estimations revealed strong association of microglial expression patterns with Aβ and tau’s synergistic impact on neuronal activity (q < 0.001). These results further support the central role of the immune system and neuroinflammatory pathways in AD pathogenesis. Neuronal dysregulation by AD pathologies also associated with neurotypical synaptic and developmental processes. In addition, we identified drug candidates from the vast LINCS library to halt or reduce the observed Aβ-tau effects on neuronal activity. Top-ranked pharmacological interventions target inflammatory, cancer and cardiovascular pathways, including specific medications undergoing clinical evaluation in AD. Our findings, based on the examination of molecular-pathological-functional interactions in humans, may accelerate the process of bringing effective therapies into clinical practice.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

On detrending stream velocity time series for robust tidal flow turbulence characterization

We investigated the impact of detrending techniques on turbulence quantities from tidal stream flow data, focusing on the autocorrelation function, $ρ_{uu}$, and velocity spectrum, $Φ(f)$. Standard detrending methods, including high-pass frequency-based and polynomial-based techniques, are examined, alongside a proposed alternative method, the empirical mode decomposition (EMD). Our results highlight that intervals of flow acceleration and deceleration, typical in tidal and riverine flows, significantly affect the estimation of turbulence quantities using high-pass frequency filtering and polynomial detrending of varying orders. These methods can strongly influence $ρ_{uu}$ and $Φ(f)$, thereby affecting the accurate estimation of derived quantities. We examine two variations of detrending data using EMD; the first removes only the EMD residue, and the second removes both the residue and the largest scale intrinsic mode function (IMF). By comparing the detrended spectra with the modeled von Kármán spectra, we demonstrate that the second variation (i.e., removing the residue and the largest scale IMF) successfully removed the large-scale trend of the data while retaining the energy of other scales.

Detrend↗
DOE OSTI2020

Planck 2018 results

We analyse the Planck full-mission cosmic microwave background (CMB) temperature and E -mode polarization maps to obtain constraints on primordial non-Gaussianity (NG). We compare estimates obtained from separable template-fitting, binned, and optimal modal bispectrum estimators, finding consistent values for the local, equilateral, and orthogonal bispectrum amplitudes. Our combined temperature and polarization analysis produces the following final results: f NL local = -0.9 ± 5.1; f NL equil = -26 ± 47; and f NL ortho = -38 ± 24 (68% CL, statistical). These results include low-multipole (4 ≤ ℓ < 40) polarization data that are not included in our previous analysis. The results also pass an extensive battery of tests (with additional tests regarding foreground residuals compared to 2015), and they are stable with respect to our 2015 measurements (with small fluctuations, at the level of a fraction of a standard deviation, which is consistent with changes in data processing). Polarization-only bispectra display a significant improvement in robustness; they can now be used independently to set primordial NG constraints with a sensitivity comparable to WMAP temperature-based results and they give excellent agreement. In addition to the analysis of the standard local, equilateral, and orthogonal bispectrum shapes, we consider a large number of additional cases, such as scale-dependent feature and resonance bispectra, isocurvature primordial NG, and parity-breaking models, where we also place tight constraints but do not detect any signal. The non-primordial lensing bispectrum is, however, detected with an improved significance compared to 2015, excluding the null hypothesis at 3.5 σ . Beyond estimates of individual shape amplitudes, we also present model-independent reconstructions and analyses of the Planck CMB bispectrum. Our final constraint on the local primordial trispectrum shape is g NL local = (-5.8 ± 6.5) × 10 4 (68% CL, statistical), while constraints for other trispectrum shapes are also determined. Exploiting the tight limits on various bispectrum and trispectrum shapes, we constrain the parameter space of different early-Universe scenarios that generate primordial NG, including general single-field models of inflation, multi-field models (e.g. curvaton models), models of inflation with axion fields producing parity-violation bispectra in the tensor sector, and inflationary models involving vector-like fields with directionally-dependent bispectra. Our results provide a high-precision test for structure-formation scenarios, showing complete agreement with the basic picture of the ΛCDM cosmology regarding the statistics of the initial conditions, with cosmic structures arising from adiabatic, passive, Gaussian, and primordial seed perturbations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Evaluation of Converter Performance Considering Static and Dynamic Device Part-to-Part Variability

This paper presents a methodology to incorporate and analyze the impact of semiconductor device part-to-part variation on power converter performance. By integrating extensive static and dynamic device characterization data with an automated compact model generation process that reflects manufacturing variability, device models with inherent variability features are utilized in converter simulations for a comprehensive assessment of performance impacts. The traditional converter performance evaluation process typically yields fixed efficiency values, often dismissing the inherent part-to-part variability caused by the manufacturing process of semiconductor devices. To address this limitation, a large population of devices was characterized to capture variations in static parameters-such as transfer, output, and capacitance characteristics-as well as dynamic behaviors, including switching losses. This data-driven approach enables the development of individual compact models, which were then integrated into converter simulations to evaluate efficiency ranges rather than single point estimated values. The converter simulation results show that part-to-part component variation can lead to significant efficiency deviations, exceeding several percentage points in high-power conversion applications. By offering a more accurate representation of converter behavior under real-world manufacturing conditions, this methodology enables designers to anticipate performance variability, improving the robustness of power converter designs.

device characterization↗
DOE OSTI2021

Distributional Deep Reinforcement Learning-Based Emergency Frequency Control

Emergency frequency control is one of the most critical approaches to maintain power system stability after major disturbances. With the increasing number of grid-connected renewable energy sources, existing model-based methods of frequency control are facing up with challenges of computational speed and scalability for large-scale systems. In this paper, the emergency frequency control problem is formulated as a Markov Decision Process (MDP) and solved through a novel Distributional Deep Reinforcement Learning (DDRL) method, namely the distributional soft actor critic (DSAC) method. Compared with other RL methods that only estimate the mean value, the proposed DSAC model estimates the distribution of value function over returns. This advancement can lead to more insights and knowledge for the agent, with the benefit of a much faster and more stable learning process, and the improved frequency control performance. Here, the simulation results on IEEE 39-bus and IEEE 118-bus systems demonstrate the effectiveness and robustness of proposed models, as well as the advantage compared to other state-of-the-art DRL algorithms.

24 POWER TRANSMISSION AND DISTRIBUTION↗