Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ensemble methods”

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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
NASA NTRS2006

Comet Mineralogy as Inferred from Infrared Spectra of Comets

For most comets, infrared (IR) spectroscopy (remote sensing) is the method through which we diagnose the mineralogy and size distribution of dust in their comae. The shape and contrast of the IR spectral features depend on the particle size: optically active minerals (absorbing of visible and near-IR solar photons) and submicron solid grains or highly porous (> 90% vacuum) grains primarily contribute to the shapes of the observed resonances. Comet mineralogies typically are determined by fitting thermal emission models of ensembles of discrete mineral grains to observed IR spectral energy distributions. The absorptivities (Q-abs) and scattering efficiencies (Q-scat) of the discrete mineral grains are computed using Mie scattering, Maxwell-Garnet mixing, Discrete Dipole Approximation, and Multi-Layered Sphere codes. These techniques when applied to crystalline minerals, specifically olivine (Mg_x, Fe_1-x)2 Si04, x>0.9, require the use of ellipsoidal shaped particles with elongated axial ratios or hollow spheres to produce the shapes of the resonances observed both from comet comae and laboratory samples. The wavelength positions of the distinct resonances from submicron-radii crystalline silicates, as well as their thermal equilibrium temperatures, constrain the crystalline olivine to have a relatively high Mg-content (x>0.9, or Fo>90). Only resonances computed for submicron Mg-rich crystalline olivine and crystalline orthopyroxene match the observed IR spectral features. However, this has led to the interpretation that micron-radii and larger crystals are absent from comet comae. Furthermore, the mass fraction of silicate crystals is dependent upon whether just the submicron portion of the size distribution is being compared or the submicron crystals compare to the aggregates of porous amorphous silicates that are computationally tractable as porous spheres. We will discuss the Deep Impact results as examples of these challenges to interpreting mid-IR spectra of short-period comets.

Wooden, Diane H.↗
DOE OSTI2025

Evaluating Sea Breezes and Associated Convective Cloud Evolution in the Model Gray Zone

We characterize convective clouds associated with sea‐breeze circulations (SBC) using multi‐agency observations and multi‐case ensemble model simulations. The focus is on assessing convective cloud lifecycle properties and their merging behavior, as well as the environmental conditions they are embedded in, particularly SBC features. In total, 46 SBC days over the Houston‐Galveston region are selected and simulated using the Weather Research and Forecasting (WRF) model at a gray zone scale with a forecast‐like parameterization setup. Advanced techniques, including change‐point detection, a Lagrangian cloud tracking method, and a newly developed cell merging and splitting detection algorithm, are applied and/or developed for this study. Our findings indicate that the WRF model at 1 km grid spacing well represents the thermodynamic conditions over the region, as well as SBC timing and intensity. However, for the associated convective cells, WRF overestimates the 30‐dBZ echo top height, cell area, and maximum radar reflectivity compared to radar observations. This overestimation is potentially due to under‐resolved entrainment processes, an overestimated merging frequency, and the overestimation of updraft intensity. Furthermore, the model exhibits a deficiency in simulating congestus clouds, showing a more rapid transition from shallow to deep convection compared to observed behavior. Moreover, observations indicate stronger, deeper, and wider clouds when merging happens. Conversely, in simulations, the merging process does not necessarily lead to higher or longer‐lived cells, as many cases experience rapid and frequent merging and splitting which may result in more variance in convective updraft velocity during the convection lifetime.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Path integral factorization and the gravitational effective action

Abstract We discuss the factorization and continuity properties of fields in the Euclidean gravitational path integral with higher dimension operators constructed from powers of the Riemann tensor. We construct the boundary terms corresponding to the microcanonical ensemble and show that the saddle point approximation to the path integral with a quasilocal energy constraint generally yields a saddle point with discontinuous temperature. This extends a previous result for the Euclidean Schwarzschild–de Sitter geometry in Einstein gravity and shows that it is robust against at least some types of quantum corrections from heavy fields. As an application, we compute the entropy of SdS in D = 4 using the BTZ method. Our result matches the entropy calculated using Wald’s formula.

Astronomy & Astrophysics↗
DOE OSTI2021

Comparing different generations of idealized solar geoengineering simulations in the Geoengineering Model Intercomparison Project (GeoMIP)

Abstract. Solar geoengineering has been receiving increased attention in recent years as a potential temporary solution to offset global warming. One method of approximating global-scale solar geoengineering in climate models is via solar reduction experiments. Two generations of models in the Geoengineering Model Intercomparison Project (GeoMIP) have now simulated offsetting a quadrupling of the CO2 concentration with solar reduction. This simulation is idealized and designed to elicit large responses in the models. Here, we show that energetics, temperature, and hydrological cycle changes in this experiment are statistically indistinguishable between the two ensembles. Of the variables analyzed here, the only major differences involve highly parameterized and uncertain processes, such as cloud forcing or terrestrial net primary productivity. We conclude that despite numerous structural differences and uncertainties in models over the past two generations of models, including an increase in climate sensitivity in the latest generation of models, the models are consistent in their aggregate climate response to global solar dimming.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1989

Determination of longitudinal aerodynamic derivatives using flight data from an icing research aircraft

A flight test was performed with the NASA Lewis Research Center's DH-6 icing research aircraft. The purpose was to employ a flight test procedure and data analysis method, to determine the accuracy with which the effects of ice on aircraft stability and control could be measured. For simplicity, flight testing was restricted to the short period longitudinal mode. Two flights were flown in a clean (baseline) configuration, and two flights were flown with simulated horizontal tail ice. Forty-five repeat doublet maneuvers were performed in each of four test configurations, at a given trim speed, to determine the ensemble variation of the estimated stability and control derivatives. Additional maneuvers were also performed in each configuration, to determine the variation in the longitudinal derivative estimates over a wide range of trim speeds. Stability and control derivatives were estimated by a Modified Stepwise Regression (MSR) technique. A measure of the confidence in the derivative estimates was obtained by comparing the standard error for the ensemble of repeat maneuvers, to the average of the estimated standard errors predicted by the MSR program. A multiplicative relationship was determined between the ensemble standard error, and the averaged program standard errors. In addition, a 95 percent confidence interval analysis was performed for the elevator effectiveness estimates, C sub m sub delta e. This analysis identified the speed range where changes in C sub m sub delta e could be attributed to icing effects. The magnitude of icing effects on the derivative estimates were strongly dependent on flight speed and aircraft wing flap configuration. With wing flaps up, the estimated derivatives were degraded most at lower speeds corresponding to that configuration. With wing flaps extended to 10 degrees, the estimated derivatives were degraded most at the higher corresponding speeds. The effects of icing on the changes in longitudinal stability and control derivat

Ranaudo, R. J.↗
NASA NTRS2002

Protein Kinase Classification with 2866 Hidden Markov Models and One Support Vector Machine

The main application considered in this paper is predicting true kinases from randomly permuted kinases that share the same length and amino acid distributions as the true kinases. Numerous methods already exist for this classification task, such as HMMs, motif-matchers, and sequence comparison algorithms. We build on some of these efforts by creating a vector from the output of thousands of structurally based HMMs, created offline with Pfam-A seed alignments using SAM-T99, which then must be combined into an overall classification for the protein. Then we use a Support Vector Machine for classifying this large ensemble Pfam-Vector, with a polynomial and chisquared kernel. In particular, the chi-squared kernel SVM performs better than the HMMs and better than the BLAST pairwise comparisons, when predicting true from false kinases in some respects, but no one algorithm is best for all purposes or in all instances so we consider the particular strengths and weaknesses of each.

Weber, Ryan↗
NASA NTRS2011

Lessons Learned from Assimilating Altimeter Data into a Coupled General Circulation Model with the GMAO Augmented Ensemble Kalman Filter

Satellite altimetry measurements have provided global, evenly distributed observations of the ocean surface since 1993. However, the difficulties introduced by the presence of model biases and the requirement that data assimilation systems extrapolate the sea surface height (SSH) information to the subsurface in order to estimate the temperature, salinity and currents make it difficult to optimally exploit these measurements. This talk investigates the potential of the altimetry data assimilation once the biases are accounted for with an ad hoc bias estimation scheme. Either steady-state or state-dependent multivariate background-error covariances from an ensemble of model integrations are used to address the problem of extrapolating the information to the sub-surface. The GMAO ocean data assimilation system applied to an ensemble of coupled model instances using the GEOS-5 AGCM coupled to MOM4 is used in the investigation. To model the background error covariances, the system relies on a hybrid ensemble approach in which a small number of dynamically evolved model trajectories is augmented on the one hand with past instances of the state vector along each trajectory and, on the other, with a steady state ensemble of error estimates from a time series of short-term model forecasts. A state-dependent adaptive error-covariance localization and inflation algorithm controls how the SSH information is extrapolated to the sub-surface. A two-step predictor corrector approach is used to assimilate future information. Independent (not-assimilated) temperature and salinity observations from Argo floats are used to validate the assimilation. A two-step projection method in which the system first calculates a SSH increment and then projects this increment vertically onto the temperature, salt and current fields is found to be most effective in reconstructing the sub-surface information. The performance of the system in reconstructing the sub-surface fields is particularly impressive for temperature, but not as satisfactory for salt.

Keppenne, Christian↗
NASA NTRS2016

Prediction of Weather Impacts on Airport Arrival Meter Fix Capacity

This paper introduces a data driven model for predicting airport arrival capacity with a look-ahead time 2-8 hour forecast. The model is suitable for air traffic flow management by explicitly investigating the impact of convective weather on airport arrival meter fix throughput. Estimation of the arrival airport capacity under arrival meter fix flow constraints due to severe weather is an important part of Air Traffic Management (ATM). Airport arrival capacity can be reduced if one or more airport arrival meter fixes are partially or completely blocked by convective weather. When the predicted airport arrival demands exceed the predicted available airport's arrival capacity for a sustained period, Ground Delay Program (GDP) operations will be triggered by ATM system. Serious imbalances between demand and capacity occur most frequently when the airport capacity is severely degraded due to either bad airport terminal surface weather or inclement convective weather around airport arrival fixes. A model that predicts the weather-impacted airport arrival meter fix throughput may help ATM personnel to plan GDP operations more efficiently. This paper identifies the characteristics of air traffic flow across arrival meter fixes at Newark Liberty International Airport (EWR). The proposed approach, based on machine-learning methods, is developed to predict the weather impacted EWR arrival Meter Fix (MF) throughput. Sector forecast coverage is used to envision the weather impact on airport arrival MF flow, and the validation is accomplished by using Convective Weather Avoidance Model (CWAM) 0.5 to 2-hour and Collaborative Convective Forecast Product (CCFP) 4 to 8-hour look-ahead forecast data for the period of April-September in 2014. Furthermore, the regression tree ensemble learning of random forests approach for translating a sector forecast coverage model to an EWR arrival meter fix throughput model is examined. The results suggest that ATM decision makers in charge of MF flow control and GDP planning may benefit from adopting the airport arrival meter capacity prediction models to estimate the inclement weather impacts.

Wang, Yao X.↗
DOE OSTI2022

Deep Learning Provides Substantial Improvements to County‐Level Fire Weather Forecasting Over the Western United States

Abstract The recent wildfires in the western United States during 2018 and 2020 caused record‐breaking fire damage and casualties. Despite remarkable advances in fire modeling and weather forecasting, it remains challenging to anticipate catastrophic wildfire events and associated damage. One key missing component is a fire weather prediction system with sufficiently long lead time capable of providing useful regional details. Here, we develop a hybrid prediction model of wildfire danger called CFS with super resolution (CFS‐SR) as a proof of concept to fill that void. The CFS‐SR model is constructed by integrating the Climate Forecast System version 2 with a deep learning (DL) technique from Single Image Super Resolution, a method widely used in enhancing image resolution. We show that for the 2018–2019 fire season, the CFS‐SR model significantly improves accuracy in forecasting fire weather at lead times of up to 7 days with an enhanced spatial resolution up to 4 km. This level of high resolution provides county‐level fire weather forecast, making it more practical for allocating resources to mitigate wildfire danger. Our study demonstrates that a proper combination of ensemble climate predictions with DL techniques can boost predictability at finer spatial scales, increasing the utility of fire weather forecasts for practical applications.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Peculiarities of the formation of photonic nanojets by a matrix of dielectric microtoroids

We report the results of theoretical modelling of photonic nanojets (PNJs) formed by laser radiation scattering on a single-layer ordered assembly of dielectric microtoroids placed on the surface of a transparent matrix (silicone film). Using the method of computational electrodynamics (FDTD), the main PNJ parameters (length, width, peak intensity) are analysed under the conditions of mutual influence of the light fields of neighbouring microparticles. It is shown that the main factor affecting the PNJ characteristics under study is the spatial configuration of the radiation-scattering particle, namely, the internal diameter of its cross section. It is found that for certain configurations of toroid placement in a cluster, a PNJ ensemble is implemented with parameters significantly better than those for a single toroid. (paper)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Influence of inhomogeneous stochasticity on the falsifiability of mean-field theories and examples from accretion disc modelling

ABSTRACT Despite spatial and temporal fluctuations in turbulent astrophysical systems, mean-field theories can be used to describe their secular evolution. However, observations taken over time scales much shorter than dynamical time scales capture a system in a single state of its turbulence ensemble. Comparing with mean-field theory can falsify the latter only if the theory is additionally supplied with a quantified precision. The central limit theorem provides appropriate estimates to the precision only when fluctuations contribute linearly to an observable and with constant coherent scales. Here, we introduce an error propagation formula that relaxes both limitations, allowing for non-linear functional forms of observables and inhomogeneous coherent scales and amplitudes of fluctuations. The method is exemplified in the context of accretion disc theories, where inhomogeneous fluctuations in the surface temperature are propagated to the disc emission spectrum – the latter being a non-linear and non-local function of the former. The derived precision depends non-monotonically on emission frequency. Using the same method, we investigate how binned spectral fluctuations in telescope data change with the spectral resolving power. We discuss the broader implications for falsifiability of a mean-field theory.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Unraveling the Dynamics of Nucleosome Arrays

The organization of genomic DNA into chromatin is a fundamental determinant of genome stability, regulation, and cellular function. Nucleosomes, the basic repeating units of chromatin, assemble into higher-order structures whose organization and heterogeneity remain difficult to characterize using conventional ensemble-averaged techniques. A key need in the field is the development of experimental approaches capable of directly visualizing nucleosome assemblies and their structural variability at the single-molecule level. This LDRD Lab-Wide project focused on establishing and evaluating atomic force microscopy (AFM)–based approaches for the characterization of nucleosome assemblies. The work emphasized experimental workflows for preparing, imaging, and assessing multi-nucleosome systems, rather than isolated single nucleosomes. Through method development and exploratory measurements, the project demonstrated the feasibility of applying scanning probe microscopy to investigate chromatin-relevant assemblies and provided preliminary insight into the strengths and limitations of this approach for future quantitative studies. Results and lessons learned from this effort were disseminated to the broader scientific community through multiple national conference presentations, helping to position LLNL for continued work in chromatin and genome organization research.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2020

Subseasonal-to-Seasonal Hindcast Skill Assessment of Ridging Events Related to Drought Over the Western United States

Persistent atmospheric ridging events centered near the western United States are associated with widespread precipitation deficits and meteorological drought. Due to the relatively low skill of- dynamical models in forecasting precipitation on subseasonal‐to‐seasonal (S2S) time scales across this region, forecasts of ridging are explored in this study as a potential bridge for early warning drought prediction. To assess skill, we evaluate deterministic and probabilistic S2S hindcasts (out to 6 week lead time) of ridging events and geopotential height anomalies over the western United States in five ensemble hindcast systems. Prediction skill for ridging events is shown to be highly variable across models. For some models, longer‐time‐averaged patterns of geopotential height anomalies across the first 6 weeks can be skillfully simulated, when evaluated probabilistically against climatology. The most skillful models show modest skill in forecasting above normal ridging occurrences at lead times of Weeks 3–4 and Weeks 5–6,with some sensitivity to the specific ridge location and method for determining skill. Using the European Centre for Medium‐Range Weather Forecasts (ECMWF) model as a case study, longer lead time forecast busts are shown to often occur under extended periods of highly active La Nina‐like tropical convection. Model errors at simulating these tropical convection features may have a disproportionately large impact and degrade downstream forecasts over the western United States. Despite the documented model shortcomings, our results highlight an opportunity to target skillful ridging forecasts from dynamical models for improving early warning drought forecasting at S2S lead times.

S2S↗
DOE OSTI2022

Hybrid Kohn-Sham+Thomas-Fermi scheme for high-temperature density functional theory

We report density functional theory has revolutionised our understanding of dense, strongly coupled systems in the condensed-matter limit. At elevated temperatures, such as those found in warm dense matter and dense plasma, the computational cost of calculating enough states to accurately describe a thermal ensemble of electrons is a major obstacle to practical Kohn-Sham calculations. This problem has recently been tackled from a number of approaches, including a plane-wave approximation for electrons at high energy. In this work, we demonstrate how the Thomas-Fermi density of states, equivalent to a local plane-wave approximation, can yield improved results for thermodynamic quantities, electron density profiles, and interatomic forces while requiring fewer states than existing methods.

36 MATERIALS SCIENCE↗
NASA NTRS2000

Low-dimensional Representation of Error Covariance

Ensemble and reduced-rank approaches to prediction and assimilation rely on low-dimensional approximations of the estimation error covariances. Here stability properties of the forecast/analysis cycle for linear, time-independent systems are used to identify factors that cause the steady-state analysis error covariance to admit a low-dimensional representation. A useful measure of forecast/analysis cycle stability is the bound matrix, a function of the dynamics, observation operator and assimilation method. Upper and lower estimates for the steady-state analysis error covariance matrix eigenvalues are derived from the bound matrix. The estimates generalize to time-dependent systems. If much of the steady-state analysis error variance is due to a few dominant modes, the leading eigenvectors of the bound matrix approximate those of the steady-state analysis error covariance matrix. The analytical results are illustrated in two numerical examples where the Kalman filter is carried to steady state. The first example uses the dynamics of a generalized advection equation exhibiting nonmodal transient growth. Failure to observe growing modes leads to increased steady-state analysis error variances. Leading eigenvectors of the steady-state analysis error covariance matrix are well approximated by leading eigenvectors of the bound matrix. The second example uses the dynamics of a damped baroclinic wave model. The leading eigenvectors of a lowest-order approximation of the bound matrix are shown to approximate well the leading eigenvectors of the steady-state analysis error covariance matrix.

Tippett, Michael K.↗
DOE OSTI2020

PbS/CdS Quantum Dot Room-Temperature Single-Emitter Spectroscopy Reaches the Telecom O and S Bands via an Engineered Stability

We synthesized PbS/CdS core/shell quantum dots (QDs) to have functional single-emitter properties for room-temperature, solid-state operation in the telecom O and S bands. Two shell-growth methods–cation exchange and successive ionic layer adsorption and reaction (SILAR)–were employed to prepare QD heterostructures with shells of 2–16 monolayers. PbS/CdS QDs were sufficiently bright and stable to resolve photoluminescence (PL) spectra representing both bands from single nanocrystals using standard detection methods, and for a QD emitting in the O-band a second-order correlation function showed strong photon antibunching, important steps toward demonstrating the utility of lead chalcogenide QDs as single-photon emitters (SPEs). Irrespective of type, few telecom-SPEs exist that are capable of such room-temperature operation. Access to single-QD spectra enabled a direct assessment of spectral line width, which was ~70–90 meV compared to much broader ensemble spectra (~300 meV). We show inhomogeneous broadening results from dispersity in PbS core sizes that increases dramatically with extended cation exchange. Quantum yields (QYs) are negatively impacted at thick shells (>6 monolayers) and, especially, by SILAR-growth conditions. Time-resolved PL measurements revealed that, with SILAR, initially single-exponential PL-decays transition to biexponential, with opening of nonradiative carrier-recombination channels. Radiative decay times are, overall, longer for core/shell QDs compared to PbS cores, which we demonstrate can be partially attributed to some core/shell sizes occupying a quasi-type II electron–hole localization regime. Finally, we demonstrate that shell engineering and the use of lower laser-excitation powers can afford significantly suppressed blinking and photobleaching. However, dependence on shell thickness comes at a cost of less-than-optimal brightness, with implications for both materials and experimental design

36 MATERIALS SCIENCE↗
DOE OSTI2025

Efficient Measurement of Length Distribution of 1D Nanoparticles in Solution via Optical Polarimetry

The efficient measurement of the length distribution of nanotubes, nanowires, and other one-dimensional (1D) nanoparticles in solution is important to enable their incorporation into materials and devices and to optimize their processing for properties of interest, such as thermal/electrical conductivity or mechanical strength, in suspensions and composites. We report an electric-field (E-field)-assisted optical-polarimetry technique to measure the length distribution of ensembles of high-aspect-ratio particles in dilute suspension. The degree of alignment of polarizable 1D particles suspended in a fluid under Brownian motion explicitly depends on the E-field strength and the particle length. We show that it is possible to extract the length distribution of 1D nanoparticles suspended in an insulating fluid by applying a range of E-fields and using optical polarimetry to measure the corresponding alignment order parameter. Notably, the method is relatively insensitive to the diameter of the 1D particles, which can be poorly known or vary within a sample. The technique is validated with silver nanowires and carbon nanotubes of known lengths, as well as polymer-depletion-length-separated single-wall carbon nanotube samples with length distributions independently measured with analytical ultracentrifugation. Finally, we demonstrate the ability of the optical-polarimetry technique to quantify changes in the length distribution of ultranarrow, sub-nanometer-diameter single-wall carbon nanotubes under different types and durations of ultrasonication. Within its range of applicability (polarizable 1D nanoparticles in the 0.5 to 15 μm length range, constrained by the voltage stability of the media and the suspended particles), the E-field-assisted optical-polarimetry method is a particularly efficient and accurate method to measure the length distribution of nanowires and nanotubes in suspension.

1D nanoparticles↗
DOE OSTI2021

AFM-IR and s-SNOM-IR measurements of chemically addressable monolayers on Au nanoparticles

The performance of catalysts depends on their nanoscale properties, and local variations in structure and composition can have a dramatic impact on the catalytic reactivity. Therefore, probing the localized reactivity of catalytic surfaces using high spatial resolution vibrational spectroscopy, such as infrared (IR) nanospectroscopy and tip-enhanced Raman spectroscopy, is essential for mapping their reactivity pattern. Two fundamentally different scanning probe IR nanospectroscopy techniques, namely, scattering-type scanning near-field optical microscopy (sSNOM) and atomic force microscopy-infrared spectroscopy (AFM-IR), provide the capabilities for mapping the reactivity pattern of catalytic surfaces with a spatial resolution of ~20 nm. Herein, we compare these two techniques with regard to their applicability for probing the vibrational signature of reactive molecules on catalytic nanoparticles. For this purpose, we use chemically addressable self-assembled molecules on Au nanoparticles as model systems. We identified significant spectral differences depending on the measurement technique, which originate from the fundamentally different working principles of the applied methods. While AFM-IR spectra provided information from all the molecules that were positioned underneath the tip, the s-SNOM spectra were more orientation-sensitive. Due to its field-enhancement factor, the s-SNOM spectra showed higher vibrational signals for dipoles that were perpendicularly oriented to the surface. The s-SNOM sensitivity to the molecular orientation influenced the amplitude, position, and signal-to-noise ratio of the collected spectra. Ensemble-based IR measurements verified that differences in the localized IR spectra stem from the enhanced sensitivity of s-SNOM measurements to the adsorption geometry of the probed molecules.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗