SEARCH · Engineering Papers

Results for “conditional random fields”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 199 records · Page 11

All sources NASA NTRS DOE OSTI

DOE OSTI2021

PastoralScape: An Environment-Driven Model of Vaccination Decision Making Within Pastoralist Groups in East Africa

Economic and cultural resilience among pastoralists in East Africa is threatened by the interconnected forces of climate change, contagious diseases spread and evolving national and international trade. A key factor in the resilience of livestock that communities depend on is human decision making regarding vaccination against prevalent diseases such as Rift Valley fever and Contagious Bovine Pleuropneumonia. This paper describes an agent-based model that couples models of disease propagation, animal health, human decision making, and external GIS data sources capturing measures of foraging condition. We describe the design of the sub-models, their coupling, and demonstrate the sensitivity of the model to parameters that relate to controllable factors such as government and NGO information sources that can influence human decision making patterns. This model is intended to form the basis upon which richer economic and human factor models can be built.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Randomized probe imaging through deep k-learning

Randomized probe imaging (RPI) is a single-frame diffractive imaging method that uses highly randomized light to reconstruct the spatial features of a scattering object. The reconstruction process, known as phase retrieval, aims to recover a unique solution for the object without measuring the far-field phase information. Typically, reconstruction is done via time-consuming iterative algorithms. In this work, we propose a fast and efficient deep learning based method to reconstruct phase objects from RPI data. The method, which we call deep k-learning, applies the physical propagation operator to generate an approximation of the object as an input to the neural network. This way, the network no longer needs to parametrize the far-field diffraction physics, dramatically improving the results. Deep k-learning is shown to be computationally efficient and robust to Poisson noise. The advantages provided by our method may enable the analysis of far larger datasets in photon starved conditions, with important applications to the study of dynamic phenomena in physical science and biological engineering.

Guo, Zhen (ORCID:0000000213473451)↗

DOE OSTI2020

Effective entropy of quantum fields coupled with gravity

Entanglement entropy, or von Neumann entropy, quantifies the amount of uncertainty of a quantum state. For quantum fields in curved space, entanglement entropy of the quantum field theory degrees of freedom is well-defined for a fixed background geometry. In this paper, we propose a generalization of the quantum field theory entanglement entropy by including dynamical gravity. The generalized quantity named effective entropy, and its Renyi entropy generalizations, are defined by analytic continuation of a replica calculation. The replicated theory is defined as a gravitational path integral with multiple copies of the original boundary conditions, with a co-dimension-2 brane at the boundary of region we are studying. We discuss different approaches to define the region in a gauge invariant way, and show that the effective entropy satisfies the quantum extremal surface formula. When the quantum fields carry a significant amount of entanglement, the quantum extremal surface can have a topology transition, after which an entanglement island region appears. Our result generalizes the Hubeny-Rangamani-Takayanagi formula of holographic entropy (with quantum corrections) to general geometries without asymptotic AdS boundary, and provides a more solid framework for addressing problems such as the Page curve of evaporating black holes in asymptotic flat spacetime. We apply the formula to two example systems, a closed two-dimensional universe and a four-dimensional maximally extended Schwarzchild black hole. We discuss the analog of the effective entropy in random tensor network models, which provides more concrete understanding of quantum information properties in general dynamical geometries. We show that, in absence of a large boundary like in AdS space case, it is essential to introduce ancilla that couples to the original system, in order for correctly characterizing quantum states and correlation functions in the random tensor network. Using the superdensity operator formalism, we study the system with ancilla and show how quantum information in the entanglement island can be reconstructed in a state-dependent and observer-dependent map. We study the closed universe (without spatial boundary) case and discuss how it is related to open universe.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2015

Enhanced focus steering abilities of multi-element therapeutic arrays operating in nonlinear regimes

Steering abilities of a typical HIFU therapeutic array operated in linear and nonlinear regimes were compared using numerical simulation with the 3D Westervelt equation. The array included 256 elements of 1.2 MHz frequency and 6.6 mm diameter distributed in a quasi-random pattern over a spherical shell with a 130 mm aperture and a focal length of 120 mm. In the case of linear focusing, thermal effects are proportional to the intensity level and the criterion for safe array operation is that the intensity in the grating lobes should be less than 10% of the intensity in the main focus. In the case of nonlinear focusing, the heating effect is no longer proportional to intensity; therefore the heat deposition rate was chosen as the relevant metric, using the same 10% threshold for the secondary lobe in comparison with the focal maximum. When steering the focus, the same linearly predicted intensity level at the main focus was maintained by increasing the array power. Numerical simulations of the acoustic field were performed for nonlinear propagation both in water and in tissue. It was shown that for shock-forming conditions in the main focus, the steering range of safe electronic focusing is larger than that for linear propagation conditions. Nonlinear sonication regimes therefore can be used to enlarge tissue volumes that can be sonicated using electronic steering of the focus of HIFU arrays.

42 ENGINEERING↗

DOE OSTI2016

SU-F-T-68: Characterizes of Microdetectors in Electron Beam Dosimetry

Purpose: Electron beam dosimetry requires high resolution data due to finite range that can be accomplished with small volume detectors. The small-field used in advance technologies in photon beam has created a market for microdetectors, however characteristics are significantly variable in photon beams and relatively unknown in electron beam that is investigated in this study. Methods: Among nearly 2 dozen microdetectors that have been investigated in small fields of photon beam, two popular detectors (microDiamond 60019 (PTW)) and W1 plastic scintillator detector (Standard Imaging)) that are tissue equivalent and have very small sensitive volume are selected. Electron beams from Varian linear accelerators were used to investigate dose linearity dose rate dependence, energy dependence, depth dose and profiles in a reference condition in a water phantom. For W1 that has its own Supermax electrometer point by point measurements were performed. For microDiamond, a PTW-scanning tank was used for both scanning and point dose measurements. Results: W1 detector showed excellent dose linearity (r{sup 2} =1.0) from 5–500 MU either with variation of dose rate or beam energy. Similar findings were also observed for microdiamond with r{sup 2}=1.0. Percent variations in dose/MU for W1 and microDiamond were 0.2–1.1% and 0.4–1.2%, respectively among dose rate and beam energy. This variation was random for microDiamond, whereas it decreased with beam energy and dose rate for W1. The depth dose and profiles were within ±1 mm for both detectors. Both detectors did not show any energy dependence in electron beams. Conclusion: Both microDiamond and W1 detectors provided superior characteristics of beam parameters in electron beam including dose, dose rate linearity and energy independence. Both can be used in electron beam except W1 require point by point measurements and microdiamond requires 1500 MU for initial quenching.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2016

Simulating the large-scale structure of HI intensity maps

Intensity mapping of neutral hydrogen (HI) is a promising observational probe of cosmology and large-scale structure. We present wide field simulations of HI intensity maps based on N-body simulations of a 2.6 Gpc / h box with 2048{sup 3} particles (particle mass 1.6 × 10{sup 11} M{sub ⊙} / h). Using a conditional mass function to populate the simulated dark matter density field with halos below the mass resolution of the simulation (10{sup 8} M{sub ⊙} / h < M{sub halo} < 10{sup 13} M{sub ⊙} / h), we assign HI to those halos according to a phenomenological halo to HI mass relation. The simulations span a redshift range of 0.35 ∼< z ∼< 0.9 in redshift bins of width Δ z ≈ 0.05 and cover a quarter of the sky at an angular resolution of about 7'. We use the simulated intensity maps to study the impact of non-linear effects and redshift space distortions on the angular clustering of HI. Focusing on the autocorrelations of the maps, we apply and compare several estimators for the angular power spectrum and its covariance. We verify that these estimators agree with analytic predictions on large scales and study the validity of approximations based on Gaussian random fields, particularly in the context of the covariance. We discuss how our results and the simulated maps can be useful for planning and interpreting future HI intensity mapping surveys.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Karhunen–Loève deep learning method for surrogate modeling and approximate Bayesian parameter estimation

We evaluate the performance of the Karhunen-Loève Deep Neural Network (KL-DNN) framework for surrogate modeling and approximate Bayesian parameter estimation in partial differential equation models. In the surrogate model, the Karhunen-Loève (KL) expansions are used for the dimensionality reduction of the number of unknown parameters and variables, and a deep neural network is employed to relate the reduced space of parameters to that of the state variables. The KL-DNN surrogate model is used to formulate a maximum-a-posteriori-like least-squares problem, which is randomized to draw samples of the posterior distribution of the parameters. We test the proposed framework for a hypothetical unconfined aquifer via comparison with the forward MODFLOW and inverse PEST++ iterative ensemble smoother (IES) solutions as well as the state-of-the-art Fourier neural operator (FNO) and deep operator networks (DeepONets) operator learning surrogate models. Our results show that the KL-DNN surrogate model outperforms FNO and DeepONet for forward predictions. For solving inverse problems, the randomized algorithm provides the same or more accurate Bayesian predictions of the parameters than IES as evidenced by the higher log-predictive probability of both the estimated parameter field and the forecast hydraulic head. The posterior mean obtained from the randomized algorithm is closer to the reference parameter field than that obtained with FNO as the maximum a posteriori estimate.

Approximate Bayesian inference↗

DOE OSTI2015

Compression induced phase transition of nematic brush: A mean-field theory study

Responsive behavior of polymer brush to the external compression is one of the most important characters for its application. For the flexible polymer brush, in the case of low grafting density, which is widely studied by the Gaussian chain model based theory, the compression leads to a uniform deformation of the chain. However, in the case of high grafting density, the brush becomes anisotropic and the nematic phase will be formed. The normal compression tends to destroy the nematic order, which leads to a complex responsive behaviors. Under weak compression, chains in the nematic brush are buckled, and the bending energy and Onsager interaction give rise to the elasticity. Under deep compression, the responsive behaviors of the nematic polymer brush depend on the chain rigidity. For the compressed rigid polymer brush, the chains incline to re-orientate randomly to maximize the orientational entropy and its nematic order is destroyed. For the compressed flexible polymer brush, the chains incline to fold back to keep the nematic order. A buckling-folding transition takes place during the compressing process. For the compressed semiflexible brush, the chains are collectively tilted to a certain direction, which leads to the breaking of the rotational symmetry in the lateral plane. These responsive behaviors of nematic brush relate to the properties of highly frustrated worm-like chain, which is hard to be studied by the traditional self-consistent field theory due to the difficulty to solve the modified diffusion equation. To overcome this difficulty, a single chain in mean-field theory incorporating Monte Carlo simulation and mean-field theory for the worm-like chain model is developed in present work. This method shows high performance for entire region of chain rigidity in the confined condition.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

A survey of unsupervised learning methods for high-dimensional uncertainty quantification in black-box-type problems

Constructing surrogate models for uncertainty quantification (UQ) on complex partial differential equations (PDEs) having inherently high-dimensional O(10 n ), n ≥ 2, stochastic inputs (e.g., forcing terms, boundary conditions, initial conditions) poses tremendous challenges. The “curse of dimensionality” can be addressed with suitable unsupervised learning techniques used as a pre-processing tool to encode inputs onto lower-dimensional subspaces while retaining its structural information and meaningful properties. In this work, we review and investigate thirteen dimension reduction methods including linear and nonlinear, spectral, blind source separation, convex and non-convex methods and utilize the resulting embeddings to construct a mapping to quantities of interest via polynomial chaos expansions (PCE). Here, we refer to the general proposed approach as manifold PCE (m-PCE), where manifold corresponds to the latent space resulting from any of the studied dimension reduction methods. To investigate the capabilities and limitations of these methods we conduct numerical tests for three physics-based systems (treated as black-boxes) having high-dimensional stochastic inputs of varying complexity modeled as both Gaussian and non-Gaussian random fields to investigate the effect of the intrinsic dimensionality of input data. We demonstrate both the advantages and limitations of the unsupervised learning methods and we conclude that a suitable m-PCE model provides a cost-effective approach compared to alternative algorithms proposed in the literature, including recently proposed expensive deep neural network-based surrogates and can be readily applied for high-dimensional UQ in stochastic PDEs.

42 ENGINEERING↗

DOE OSTI2016

Clinical Experience and Evaluation of Patient Treatment Verification With a Transit Dosimeter

Purpose: To prospectively evaluate a protocol for transit dosimetry on a patient population undergoing intensity modulated radiation therapy (IMRT) and to assess the issues in clinical implementation of electronic portal imaging devices (EPIDs) for treatment verification. Methods and Materials: Fifty-eight patients were enrolled in the study. Amorphous silicon EPIDs were calibrated for dose and used to acquire images of delivered fields. Measured EPID dose maps were back-projected using the planning computed tomographic (CT) images to calculate dose at prespecified points within the patient and compared with treatment planning system dose offline using point dose difference and point γ analysis. The deviation of the results was used to inform future action levels. Results: Two hundred twenty-five transit images were analyzed, composed of breast, prostate, and head and neck IMRT fields. Patient measurements demonstrated the potential of the dose verification protocol to model dose well under complex conditions: 83.8% of all delivered beams achieved the initial set tolerance level of Δ{sub D} of 0 ± 5 cGy or %Δ{sub D} of 0% ± 5%. Importantly, the protocol was also sensitive to anatomic changes and spotted that 3 patients from 20 measured prostate patients had undergone anatomic change in comparison with the planning CT. Patient data suggested an EPID-reconstructed versus treatment planning system dose difference action level of 0% ± 7% for breast fields. Asymmetric action levels were more appropriate for inversed IMRT fields, using absolute dose difference (−2 ± 5 cGy) or summed field percentage dose difference (−6% ± 7%). Conclusions: The in vivo dose verification method was easy to use and simple to implement, and it could detect patient anatomic changes that impacted dose delivery. The system required no extra dose to the patient or treatment time delay and so could be used throughout the course of treatment to identify and limit systematic and random errors in dose delivery for patient groups.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2026

Measurement of Local Impedance Characteristics of In-Plane PEM Electrolyzer Component Features via Segmented Cell

Segmented cells are a class of advanced diagnostic devices that enable the measurement of current distribution within the active area of an electrochemical cell, allowing for the evaluation of localized impacts from operating conditions, flow fields and component inhomogeneities. Segmented cell devices have been applied to support R&D efforts for the advancement of polymer electrolyte membrane fuel cell technologies and more recently proton exchange membrane water electrolysis (PEMWE). Beyond current density distributions, segmented cells have been successfully deployed to collect the distribution of other relevant parameters such as high frequency resistance and electrode potentials. Very recently, a promising impedance method has been reported for water electrolysis where a one-dimensional segmented cell is coupled to a multichannel potentiostat, to enable the measurement of local cell impedances, further expanding the diagnostic capabilities of the segmented cell. In this study, we will present recent efforts to accurately measure local impedance characteristics for a PEMWE cell using a two-dimensional segmented cell device coupled with a multichannel potentiostat. The distribution of local resistances, i.e. high frequency, charge transfer, mass transport and catalyst layer, and iR-free overpotentials were successfully measured after eliminating sources of random and systematic error. We will discuss the application of this diagnostic on a cell containing a 1 cm feature simulating an ionomer skin. The results from these tests give insights if such an irregularity needs to be classified as a defect, and highlight the usefulness of coupling impedance spectroscopy with segmented cells for holistic spatial diagnostics of electrochemical devices.

08 HYDROGEN↗

DOE OSTI2019

1D PIC-DSMC analysis of a high-pressure nanosecond pulse discharge breakdown in helium

Kinetic simulations of plasma phenomena during and after formation of the conductive plasma channel of a nanosecond pulse discharge are analyzed and compared to existing experimental measurements. Particle-in-cell with direct simulation Monte Carlo collisions (PIC-DSMC) modeling is used to analyze a discharge in helium at 200 Torr and 300 K over a 1 cm gap. The analysis focuses on physics that would not be reproduced by fluid models commonly used at this high number density and collisionality, specifically non-local and stochastic phenomena. Similar analysis could be used to improve the predictive capability of lower fidelity or reduced order models. First, the modeling results compare favorably with experimental measurements of electron number density, temperature, and 1D electron energy distribution function at the same conditions. Second, it is shown that the ionization wave propagates in a stochastic, stepwise manner, dependent on rare, random ionization events ahead of the ionization wave when the ionization fraction in front of the ionization wave is very low, analagous to the stochastic branching of streamers in 3D. Third, analysis shows high-energy runaway electrons accelerated in the cathode layer produce electron densities in the negative glow region over an order of magnitude above those in the positive column. Finally, future work to develop reduced order models of these two phenomena would improve the accuracy of fluid plasma models without the cost of PIC-DSMC simulations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

An improved method to access initial states in relativistic heavy-ion collisions

Abstract Observables in heavy-ion collisions are generally categorized into centralities, which reflect an average over events within a range of impact parameters including a wide variety of initial-state configurations. A multiple binning method using spectator neutrons within each centrality has been shown to provide access to events with rare initial-state conditions. This work suggests an improvement in quantifying the difference between standard centrality and spectator neutron binning towards accessing the initial-state properties. A selection of events with higher initial-state density at a fixed participating nucleon number was observed to result in larger final-state particle production and smaller elliptic flow. The relative difference between observables in centrality and spectator binning shows reduced sensitivity for the observables dominated by impact parameter fluctuations in the initial state, such as triangular flow. This property renders the spectator binning method a good candidate for separating geometric contributions from random fluctuations in the initial state towards final-state observables.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2017

Binary Black Hole Mergers from Field Triples: Properties, Rates, and the Impact of Stellar Evolution

We consider the formation of binary black hole (BH) mergers through the evolution of field massive triple stars. In this scenario, favorable conditions for the inspiral of a BH binary are initiated by its gravitational interaction with a distant companion, rather than by a common-envelope phase invoked in standard binary evolution models. We use a code that follows self-consistently the evolution of massive triple stars, combining the secular triple dynamics (Lidov–Kozai cycles) with stellar evolution. After a BH triple is formed, its dynamical evolution is computed using either the orbit-averaged equations of motion, or a high-precision direct integrator for triples with weaker hierarchies for which the secular perturbation theory breaks down. Most BH mergers in our models are produced in the latter non-secular dynamical regime. We derive the properties of the merging binaries and compute a BH merger rate in the range (0.3–1.3) Gpc{sup −3} yr{sup −1}, or up to ≈2.5 Gpc{sup −3} yr{sup −1} if the BH orbital planes have initially random orientation. Finally, we show that BH mergers from the triple channel have significantly higher eccentricities than those formed through the evolution of massive binaries or in dense star clusters. Measured eccentricities could therefore be used to uniquely identify binary mergers formed through the evolution of triple stars. While our results suggest up to ≈10 detections per year with Advanced-LIGO, the high eccentricities could render the merging binaries harder to detect with planned space based interferometers such as LISA.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2017

Adiabatic phase-conserving processes for executing quantum operations with ultracold atoms

We have studied the regimes of deterministic single-atom Rydberg excitation in the conditions of Rydberg blockade and the methods of compensation for the dynamic phase of the wave function during the adiabatic passage. Using these methods, we have proposed schemes of single-qubit and two-qubit quantum states with mesoscopic atomic ensembles containing a random number of atoms, considred as quibits. The double adiabatic passage of the Förster resonance for two interacting atoms with a deterministic phase shift can be used for the implementation of two-qubit gates with reduced sensitivity of the gate fidelity to the fluctuations of the interatomic distance. (paper)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2018

A Galerkin method with two-dimensional Haar basis functions for the computation of the Karhunen–Loève expansion

We study the numerical approximation of a homogeneous Fredholm integral equation of second kind associated with the Karhunen–Loève expansion of Gaussian random fields. We develop, validate, and discuss an algorithm based on the Galerkin method with two-dimensional Haar wavelets as basis functions. The shape functions are constructed from the orthogonal decomposition of tensor-product spaces of one-dimensional Haar functions, and a recursive algorithm is employed to compute the matrix of the discrete eigenvalue system without the explicit calculation of integrals, allowing the implementation of a fast and efficient algorithm that provides considerable reduction in CPU time, when compared with classical Galerkin methods. Numerical experiments confirm the convergence rate of the method and assess the approximation error and the sparsity of the eigenvalue system when the wavelet expansion is truncated. We illustrate the numerical solution of a diffusion problem with random input data with the present method. In this problem, accuracy was retained after dropping the coefficients below a threshold value that was numerically determined. A similar method with scaling functions rather than wavelet functions does not need a discrete wavelet transform and leads to eigenvalue systems with better conditioning but lower sparsity.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2017

High Pressure Hydrogen Pressure Relief Devices: Accelerated Life Testing and Application Best Practices

Pressure relief devices (PRDs ) are used to protect high pressure systems from burst failure caused by overpressurization. Codes and standards require the use of PRDs for the safe design of many pressurized systems. These systems require high reliability due to the risks associated with a burst failure. Hydrogen service can increase the risk of PRD failure due to material property degradation caused by hydrogen attack. The National Renewable Energy Laboratory (NREL) has conducted an accelerated life test on a conventional spring loaded PRD. Based on previous failures in the field, the nozzles specific to these PRDs are of particular interest. A nozzle in a PRD is a small part that directs the flow of fluid toward the sealing surface to maintain the open state of the valve once the spring force is overcome. The nozzle in this specific PRD is subjected to the full tensile force of the fluid pressure. These nozzles are made from 440C material, which is a type of hardened steel that is commonly chosen for high pressure applications because of its high strength properties. In a hydrogen environment, however, 440C is considered a worst case material since hydrogen attack results in a loss of almost all ductility and thus 440C is prone to fatigue and material failure. Accordingly, 440C is not recommended for hydrogen service. Conducting an accelerated life test on a PRD with 440C material provides information on necessary and sufficient conditions required to produce crack initiation and failure. The accelerated life test also provides information on other PRD failure modes that are somewhat statistically random in nature.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2017

Stellar clusterings around 'isolated' Massive YSOs in the LMC

Observations suggest that there is a significant fraction of O stars in the field of the Milky Way that appear to have formed in isolation or in low-mass clusters (<100 M{sub ⊙}). The existence of these high-mass stars that apparently formed in the field challenges the generally accepted paradigm, which requires star formation to occur in clustered environments. In order to understand the physical conditions for the formation of these stars, it is necessary to observe isolated high-mass stars while they are still forming. With the Hubble Space Telescope, we observe the seven most isolated massive (>8 M{sub ⊙}) young stellar objects (MYSOs) in the Large Magellanic Cloud. The observations show that while these MYSOs are remote from other MYSOs, OB associations, and even known giant molecular clouds, they are actually not isolated at all. Imaging reveals ∼100 to several hundred pre-main-sequence (PMS) stars in the vicinity of each MYSO. These previously undetected PMS stars form prominent compact clusters around the MYSOs, and in most cases they are also distributed sparsely across the observed regions. Contrary to what previous high-mass field star studies show, these observations suggest that high-mass stars may not be able to form in clusters with masses less than 100 M{sub ⊙}. If these MYSOs are indeed the best candidates for isolated high-mass star formation, then the lack of isolation is at odds with random sampling of the initial mass function. Moreover, while isolated MYSOs may not exist, we find evidence that isolated clusters containing O stars can exist, which in itself is rare.

79 ASTRONOMY AND ASTROPHYSICS↗