Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Randomized 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 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Stochastic Framework for Optimal Control of Planetary Reentry Trajectories Under Multilevel Uncertainties

We present a novel stochastic optimal control framework that accounts for various types of uncertainties, with application to reentry trajectory planning. The formulation of the optimal trajectory control problem is presented in the context of an indirect method where a functional objective associated with the terminal vehicle speed is to be minimized. Uncertain input parameters in the optimal trajectory control model, including aerodynamic parameters and initial and terminal conditions, are modeled as aleatory random variables, while the statistical parameters of these aleatory distributions are themselves random variables. The parametric and model uncertainties are simultaneously propagated through an extended polynomial chaos expansion (EPCE) formalism. Several metrics are described to evaluate response statistics and presented as insightful tools for robust decision making. Specifically, the response probability density function (PDF) reflecting influence of both epistemic and aleatory uncertainties is obtained. By sampling over the random variables representing model error, an ensemble of response PDFs is generated and the associated failure probability is estimated as a random variable with its own polynomial chaos expansion. Besides, the sensitivity index functions of response PDF with respect to the statistical parameters are evaluated. Coupling parametric and model uncertainties within the EPCE framework leads to a robust and efficient paradigm for multilevel uncertainty propagation and PDF characterization in general optimal control problems.

Engineering↗
DOE OSTI2021

A Robust Approach to Sequential Information Theoretic Planning

In many sequential planning applications a natural approach to generating high quality plans is to maximize an information reward such as mutual information (MI). Unfortunately, MI lacks a closed form in all but trivial models, and so must be estimated. In applications where the cost of plan execution is expensive, one desires planning estimates which admit theoretical guarantees. Through the use of robust M-estimators we obtain bounds on absolute deviation of estimated MI. Moreover, we propose a sequential algorithm which integrates inference and planning by maximally reusing particles in each stage. We validate the utility of using robust estimators in the sequential approach on a Gaussian Markov Random Field wherein information measures have a closed form. Lastly, we demonstrate the benefits of our integrated approach in the context of sequential experiment design for inferring causal regulatory networks from gene expression levels. Here, our method shows improvements over a recent method which selects intervention experiments based on the same MI objective.

Zheng, Sue↗
DOE OSTI2021

A new open‐path eddy covariance method for nitrous oxide and other trace gases that minimizes temperature corrections

Abstract Low‐power, open‐path gas sensors enable eddy covariance (EC) flux measurements in remote areas without line power. However, open‐path flux measurements are sensitive to fluctuations in air temperature, pressure, and humidity. Laser‐based, open‐path sensors with the needed sensitivity for trace gases like methane (CH 4 ) and nitrous oxide (N 2 O) are impacted by additional spectroscopic effects. Corrections for these effects, especially those related to temperature fluctuations, often exceed the flux of gases, leading to large uncertainties in the associated fluxes. For example, the density and spectroscopic corrections arising from temperature fluctuations can be one or two orders of magnitude greater than background N 2 O fluxes. Consequently, measuring background fluxes with laser‐based, open‐path sensors is extremely challenging, particularly for N 2 O and gases with similar high‐precision requirements. We demonstrate a new laser‐based, open‐path N 2 O sensor and a general approach applicable to other gases that minimizes temperature‐related corrections for EC flux measurements. The method identifies absorption lines with spectroscopic effects in the opposite direction of density effects from temperature and, thus, density and spectroscopic effects nearly cancel one another. The new open‐path N 2 O sensor was tested at a corn ( Zea mays L.) field in Southwestern Michigan, United States. The sensor had an optimal precision of 0.1 ppbv at 10 Hz and power consumption of 50 W. Field trials showed that temperature‐related corrections were 6% of density corrections, reducing EC random errors by 20‐fold compared to previously examined lines. Measured open‐path N 2 O EC fluxes showed excellent agreement with those made with static chambers ( m = 1.0 ± 0.3; r 2 = .96). More generally, we identified absorption lines for CO 2 and CH 4 flux measurements that can reduce the temperature‐related corrections by 10–100 times compared to existing open‐path sensors. The proposed method provides a new direction for future open‐path sensors, facilitating the expansion of accurate EC flux measurements.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Scattering of radio frequency waves by randomly modulated density interfaces in the edge of fusion plasmas

In the scrape-off layer and the edge region of a tokamak, the plasma is strongly turbulent and scatters the radio frequency (RF) electromagnetic waves that propagate through this region. It is important to know, whether used for diagnostics or for heating and current drive, the spectral properties of these scattered RF waves. The spectral changes influences the interpretation of the diagnostic-data obtained and the current and heating profiles. A full-wave, 3D electromagnetic code ScaRF (see Papadopoulos et al. 2019) has been developed for studying the RF wave propagation through turbulent plasma. ScaRF is a finite-difference frequency-domain (FDFD) method for solving Maxwell's equations. The magnetized plasma is defined through the cold plasma, anisotropic permittivity tensor. As a result, ScaRF can be used to study the scattering of any cold plasma RF wave. It can be for the study of scattering of electron cyclotron waves in ITER-type and medium-sized tokamaks such as TCV, ASDEX-U, DIII-D. For the case of medium-sized tokamaks, there's experimental evidence that drift waves and rippling modes are present in the edge region (see Ritz et al. 1984). Hence, we study the scattering of RF waves by periodic density interfaces (plasma gratings) in the form of a superposition of spatial modes with varying periodicity and random amplitudes (see Papadopoulos et al. 2019). The power reflection coefficient (a random variable) is calculated for different realizations of the density interface. In this work, the uncertainty of the power reflection coefficient is rigorously quantified by use of the Polynomial Chaos Expansion (see Xiu & Karniadakis 2002) method in conjunction with the Smolyak sparse grid integration (see Papadopoulos et al. 2018) (PCE-SG). The PCE-SG method is proven accurate and much more efficient (roughly 2-orders of magnitude shorter execution time) compared to alternative methods such as the Monte Carlo (MC) approach.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Grassmannian Diffusion Maps--Based Dimension Reduction and Classification for High-Dimensional Data

This work introduces the Grassmannian diffusion maps (GDMaps), a novel nonlinear dimensionality reduction technique that defines the affinity between points through their representation as low-dimensional subspaces corresponding to points on the Grassmann manifold. Here, the method is designed for applications, such as image recognition and data-based classification of constrained high-dimensional data where each data point itself is a high-dimensional object (i.e., a large matrix) that can be compactly represented in a lower-dimensional subspace. The GDMaps is composed of two stages. The first is a pointwise linear dimensionality reduction wherein each high-dimensional object is mapped onto the Grassmann manifold representing the low-dimensional subspace on which it resides. The second stage is a multipoint nonlinear kernel-based dimension reduction using diffusion maps to identify the subspace structure of the points on the Grassmann manifold. To this end, an appropriate Grassmannian kernel is used to construct the transition matrix of a random walk on a graph connecting points on the Grassmann manifold. Spectral analysis of the transition matrix yields low-dimensional Grassmannian diffusion coordinates embedding the data into a low-dimensional reproducing kernel Hilbert space. Further, a novel data classification/recognition technique is developed based on the construction of an overcomplete dictionary of reduced dimension whose atoms are given by the Grassmannian diffusion coordinates. Three examples are considered. First, a "toy" example shows that the GDMaps can identify an appropriate parametrization of structured points on the unit sphere. The second example demonstrates the ability of the GDMaps to revealing the intrinsic subspace structure of high-dimensional random field data. In the last ex- ample, a face recognition problem is solved considering face images subject to varying illumination conditions, changes in face expressions, and occurrence of occlusions. The technique presented high recognition rates (i.e., 95% in the best case) using a fraction of the data required by conventional methods.

42 ENGINEERING↗
DOE OSTI2021

Global Sensitivity Analysis of a Reactive Transport Model for Mineral Scale Formation During Hydraulic Fracturing

Injection of water-based hydraulic fracturing fluid (HFF) into tight shale gas/oil formations can increase formation permeability and enhance production rates, but this process frequently causes mineral scale formation that can occlude pore space and hinder flow. To identify the most important factors that control the formation of mineral scales, we applied a novel global sensitivity analysis method—distance-based generalized sensitivity analysis (DGSA)—to a reactive transport model (RTM) that was previously built and calibrated to simulate precipitation of barite [BaSO4] and iron (hydr)oxide [Fe(OH) 3 ] in shale matrices and on fracture surfaces. Reactive transport simulations were run with model parameters randomly sampled based on assigned uncertainties. Modeling results for barite and Fe(OH)3 formation were clustered using machine-learning algorithms. A list of ranked critical input parameters was obtained after statistical quantification of cumulative distribution functions of input parameters. We found that barite formation is most sensitive to the rate of sulfate ion generation, which is determined by the pyrite dissolution rate coefficient and oxidant availability. In addition, barite formation is sensitive to the initial amounts of barite in HFF and shale, followed by barite thermodynamics/kinetics. For Fe(OH) 3 formation, the ranked factors are Fe(OH)3 precipitation rate coefficients, initial HFF pH, initial Fe(OH) 3 amount in HFF, and oxidant availability. Overall, our results provide insights into managing mineral scale formation during hydraulic fracturing to enhance production. Meanwhile, this study serves as an example of global sensitivity analysis of RTMs using the efficient, straightforward, and open-source DGSA method.

58 GEOSCIENCES↗
DOE OSTI2021

Thermal conductivity modeling of monodispersed microspheres using discrete element method

Particle beds are widely used in various systems and processes, such as particle heat exchangers, granular flow reactors, and additive manufacturing. Accurate modeling of the thermal conductivity of particle beds and understanding of their heat transfer mechanisms are important. However, previous models were based on simple cubic packing of particles, which could not accurately represent the actual heat transfer processes under certain conditions. Here, we examine the effect of the packing structure on the thermal conductivity of particle beds. We use monodispersed silica microspheres with average particle sizes ranging from 23 to 330 μm as a model material. We employ a transient hot-wire technique to measure the thermal conductivity of the particle beds with packing density of 43%–57% within a temperature range of room temperature to 500 °C and under N2 gaseous pressures of 20–760 Torr. We then use a discrete element method (DEM) to obtain the realistic packing structure of the particles, which is then fed into a finite-element model (FEM) to calculate the thermal conductivity, with the consideration of solid conduction, gas conduction, and radiation heat transfer. Our results show that the thermal conductivity model based on the more realistic random packing structure derived from the DEM shows better agreement with the experimental data compared to that based on the simple cubic-packing structure. The combined DEM and FEM methodology can serve as a useful tool to predict the effective thermal conductivity of particle beds and to quantify different heat transfer mechanisms under various conditions.

14 SOLAR ENERGY↗
DOE OSTI2022

How to estimate soil organic carbon stocks of agricultural fields? perspectives using ex-ante evaluation

Estimating soil organic carbon (SOC) stocks of agricultural fields has a range of important applications from development of sustainable management practices to monitoring carbon stocks. There are many estimation strategies with the potential for more reliable estimates of SOC stock and more efficient use of soil sampling and analysis resources, especially by leveraging readily available auxiliary information such as remote sensing. However, concrete guidance for strategy selection is lacking. This study narrows this gap with a comparison of strategies for estimating deep SOC stock (0–60 cm) in a prototypical field. Using high density SOC stock measurements and simulation, we built on past studies by 1) ex-ante evaluating a large number of strategy options, 2) using a Bayesian approach to quantify the uncertainty of the comparison, and 3) considering multiple Bayesian models to assess sensitivity to this modeling choice. We found that, using readily available auxiliary information, both balanced and stratified sampling offer substantial improvements over simple random sampling. The auxiliary information most important for this improvement is a Sentinel-2 SOC index = blue / (green × red), followed by the topographic wetness index. We found that these results are robust to the choice of mapping method, but that there is uncertainty in the magnitude of improvement. Here, we recommend future studies implement this Bayesian approach for simulated ex-ante evaluation of SOC stock estimation strategies across more fields to investigate the generalizability of these findings.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Self-Assembled Vesicles from Mixed Brush Nanoparticles in Solution

The self-assembly of binary polymer-grafted nanoparticles (NPs) in a selective solvent is investigated using coarse-grained simulations. Simulations are performed using theoretically informed Langevin dynamics (TILD), a particle-based method that employs a particle-to-mesh scheme to efficiently calculate the nonbonded interactions. Here, the particles are densely grafted with two immiscible polymers, A and B, that are permanently bound to the NP either at random grafting sites (random-grafted) or with all the A chains on one hemisphere of the NP and all the B chains on the other hemisphere (Janus-grafted). For NPs with random grafting, the polymers phase-separate on the surface of the NP to form Janus-type structures in dilute solution, even though some of the chains have to stretch around the particle to form the Janus structure. When the solvent quality is sufficiently poor for the solvophobic chains, the binary grafted NPs assemble into various structures, including double-walled vesicles. In particular, vesicles are formed when the solvophilic volume fraction is between 0.2 and 0.3, in a similar range to that required for vesicle formation in diblock copolymers in a selective solvent. For mixed-grafted NPs, there is considerable variation in the structure of each individual NP, but nevertheless, these NPs form ordered vesicles, similar to those formed by Janus-grafted NPs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

amPEPpy 1.0: a portable and accurate antimicrobial peptide prediction tool

Antimicrobial peptides (AMPs) are promising alternative antimicrobial agents. Currently, however, portable, user-friendly and efficient methods for predicting AMP sequences from genome-scale data are not readily available. In this work, we present amPEPpy, an open-source, multi-threaded command-line application for predicting AMP sequences using a random forest classifier.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Grid-based diffusion Monte Carlo for fermions without the fixed-node approximation

A diffusion Monte Carlo algorithm is introduced that can determine the correct nodal structure of the wave function of a few-fermion system and its ground-state energy without an uncontrolled bias. This is achieved by confining signed random walkers to the points of a uniform infinite spatial grid, allowing them to meet and annihilate one another to establish the nodal structure without the fixed-node approximation. An imaginary-time propagator is derived rigorously from a discretized Hamiltonian, governing a non-Gaussian, sign-flipping, branching, and mutually annihilating random walk of particles. The accuracy of the resulting stochastic representations of a fermion wave function is limited only by the grid and imaginary-time resolutions and can be improved in a controlled manner. Here, the method is tested for a series of model problems including fermions in a harmonic trap as well as the He atom in its singlet or triplet ground state. For the latter case, the energies approach from above with increasing grid resolution and converge within 0.015 E h of the exact basis-set-limit value for the grid spacing of 0.08 a.u. with a statistical uncertainty of 10 –5 E h without an importance sampling or Jastrow factor.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Deep-learning-based image registration for nano-resolution tomographic reconstruction

Nano-resolution full-field transmission X-ray microscopy has been successfully applied to a wide range of research fields thanks to its capability of non-destructively reconstructing the 3D structure with high resolution. Due to constraints in the practical implementations, the nano-tomography data is often associated with a random image jitter, resulting from imperfections in the hardware setup. Without a proper image registration process prior to the reconstruction, the quality of the result will be compromised. Here a deep-learning-based image jitter correction method is presented, which registers the projective images with high efficiency and accuracy, facilitating a high-quality tomographic reconstruction. This development is demonstrated and validated using synthetic and experimental datasets. We report the method is effective and readily applicable to a broad range of applications. Together with this paper, the source code is published and adoptions and improvements from our colleagues in this field are welcomed.

deep learning↗
DOE OSTI2020

Distinctiveness of genes contributing to growth of Pseudomonas syringae in diverse host plant species

A variety of traits are necessary for bacterial colonization of the interior of plant hosts, including well-studied virulence effectors as well as other phenotypes contributing to bacterial growth and survival within the apoplast. High-throughput methods such as transposon sequencing (TnSeq) are powerful tools to identify such genes in bacterial pathogens. However, there is little information as to the distinctiveness of traits required for bacterial colonization of different hosts. Here, we utilize randomly barcoded TnSeq (RB-TnSeq) to identify the genes that contribute to the ability of Pseudomonas syringae strain B728a to grow within common bean (Phaseolus vulgaris), lima bean (Phaseolus lunatus), and pepper (Capsicum annuum); species representing two different plant families. The magnitude of contribution of most genes to apoplastic fitness in each of the plant hosts was similar. However, 50 genes significantly differed in their fitness contributions to growth within these species. These genes encoded proteins in various functional categories including polysaccharide synthesis and transport, amino acid metabolism and transport, cofactor metabolism, and phytotoxin synthesis and transport. Six genes that encoded unannotated, hypothetical proteins also contributed differentially to growth in these hosts. The genetic repertoire of a relatively promiscuous pathogen such as P. syringae may thus be shaped, at least in part, by the conditional contribution of some fitness determinants.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Model Exploration of an Information-Based Healthcare Intervention Using Parallelization and Active Learning

This paper describes the application of a large-scale active learning method to characterize the parameter space of a computational agent-based model developed to investigate the impact of CommunityRx, a clinical information-based health intervention that provides patients with personalized information about local community resources to meet basic and self-care needs. Additionally, the diffusion of information about community resources and their use is modeled via networked interactions and their subsequent effect on agents' use of community resources across an urban population. A random forest model is iteratively fitted to model evaluations to characterize the model parameter space with respect to observed empirical data. We demonstrate the feasibility of using high-performance computing and active learning model exploration techniques to characterize large parameter spaces; by partitioning the parameter space into potentially viable and non-viable regions, we rule out regions of space where simulation output is implausible to observed empirical data. We argue that such methods are necessary to enable model exploration in complex computational models that incorporate increasingly available micro-level behavior data. We provide public access to the model and high-performance computing experimentation code.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

S AP F LOWER : an automated tool for sap flow data preprocessing, gap-filling, and analysis using deep learning

Sap flow, a critical process in plant water use and ecosystem water cycles, is often measured using thermal dissipation probes (TDP) due to their ease of installation and continuous data collection. However, sap flow data frequently include noise, outliers, and gaps, creating challenges for analysis and requiring substantial manual processing. We developed S AP F LOWER , a tool that automates data preprocessing, model training, gap-filling, sapwood area scaling and modeling, and water use analysis. It integrates autocleaning, machine learning and deep learning models (e.g. random forest, Gaussian process regression, long short-term memory (LSTM), bidirectional LSTM (BiLSTM)), and efficient workflows to process sap flow data. S AP F LOWER can remove over 90% of noisy data while preserving legitimate variations and achieve high accuracy in gap-filling based on user-determined parameters. Random forest, LSTM, and BiLSTM models reduced root mean square error to 10% or less for long-term gaps. Model training and prediction can be performed efficiently within seconds. S AP F LOWER significantly enhances the efficiency and accessibility of TDP data analysis by automating complex tasks, enabling researchers without programming expertise to employ advanced techniques. Future improvements will focus on species-specific corrections for TDP and support for additional measurement methods. S AP F LOWER is openly available on GitHub (https://github.com/JiaxinWang123/SapFlower) and Zenodo (doi: 10.5281/zenodo.13665919).

ecosystem water balance↗
DOE OSTI2024

Nowcasting Earthquakes With Stochastic Simulations: Information Entropy of Earthquake Catalogs

Earthquake nowcasting has been proposed as a means of tracking the change in large earthquake potential in a seismically active area. The method was developed using observable seismic data, in which probabilities of future large earthquakes can be computed using Receiver Operating Characteristic methods. Furthermore, analysis of the Shannon information content of the earthquake catalogs has been used to show that there is information contained in the catalogs, and that it can vary in time. So an important question remains, where does the information originate? In this paper, we examine this question using stochastic simulations of earthquake catalogs. Our catalog simulations are computed using an Earthquake Rescaled Aftershock Seismicity (“ERAS”) stochastic model. This model is similar in many ways to other stochastic seismicity simulations, but has the advantage that the model has only 2 free parameters to be set, one for the aftershock (Omori-Utsu) time decay, and one for the aftershock spatial migration away from the epicenter. Generating a simulation catalog and fitting the two parameters to the observed catalog such as California takes only a few minutes of wall clock time. While clustering can arise from random, Poisson statistics, we show that significant information in the simulation catalogs arises from the “non-Poisson” power-law aftershock clustering, implying that the practice of de-clustering observed catalogs may remove information that would otherwise be useful in forecasting and nowcasting. We also show that the nowcasting method provides similar results with the ERAS model as it does with observed seismicity.

58 GEOSCIENCES↗
DOE OSTI2020

Robustness of Optical Response for Self-Assembled Plasmonic Metamaterials with Morphological Disorder and Surface Roughness

Bottom-up fabrication of metallized biotemplated nanostructures to form specific plasmonic nanoresonators holds promise as a means of achieving large-scale optical metamaterials. However, in contrast to top-down methods, the stochastic growth of self-assembled nanoresonators is prone to significant disorder and surface roughness, which naturally raise an important question about the robustness of their resonant properties in terms of structural imperfections. An aggregated-random-sphere model is developed to mimic the nucleated growth of metallized DNA origami assembly, leading to meta-atoms with realistic, experimentally observed morphological disorder and surface roughness. Using the well-known split-ring-resonator (SRR) motif as an example, the resonant properties of meta-atoms under different levels of roughness are investigated and a strong tolerance of optical response against morphological disorder is revealed. It is found that in SRRs, even with dramatic roughness introduced, the expected resonances are still observed, despite broadening line shapes compared to ideal smooth structure. Only for extreme disorder, which causes drastic segmentation of SRRs, does the resonant response disappear. Finally, the demonstrations are very encouraging for the prospects of bottom-up fabrication toward versatile functional metamaterials and metadevices.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Interlaboratory Study of Digital Volume Correlation Error Due to X-Ray Computed Tomography Equipment and Scan Parameters: an Update from the DVC Challenge

Background: The quality of Digital Volume Correlation (DVC) full-field displacement measurements depends directly on the characteristics of the X-ray Computed Tomography (XCT) equipment, and scan procedures used to acquire the tomographic images. Objective: In this work, we seek to experimentally study the effects of XCT equipment and tomographic scan procedures on the quality of these images for DVC analysis, and to survey the level of DVC error that may be achieved using standard XCT operating procedures. Methods: Six participants in an interlaboratory study acquired high-quality XCT scans of a syntactic foam before and after rigid body motion. The resulting images were correlated using commercial DVC software to quantify error sources due to random image noise, reconstruction artifacts, as well as systematic spatial or temporal distortion. Results: In the absence of rigid body motion, the standard deviation of the displacement measurements ranged from 0.012 to 0.043 voxels using a moderate subvolume size, indicating that subvoxel measurement resolution could readily be achieved with a variety of XCT equipment and scan recipes. Comparison of consecutive scans without rigid body motion showed transient dilatational displacement gradients due to self-heating of the X-ray source and/or thermal expansion of the foam. Evaluation of the scans after rigid body motion showed significant, machine-specific spatial distortion in the displacement fields of up to 0.5 voxels; new approaches to remove this error need to be developed. Conclusions: Analysis of the scan protocols used in the interlaboratory study, as well as a complementary parametric sensitivity study, showed that the DVC error was strongly influenced by the XCT equipment, but could be mitigated by adjusting the total scan duration.

47 OTHER INSTRUMENTATION↗