Search NASA⌕ Search

SEARCH · Search NASA

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.

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

At least 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Towards robust surrogate models: Benchmarking machine learning approaches to expediting phase field simulations of brittle fracture

Data-driven approaches have the potential to make modeling complex, nonlinear physical phenomena significantly more computationally tractable. For example, computational modeling of fracture is a core challenge where machine learning techniques have the potential to provide a much needed speedup that would enable progress in areas such as multi-scale modeling and uncertainty quantification. Currently, phase field modeling (PFM) of fracture is one such approach that offers a convenient variational formulation to model crack nucleation, branching and propagation. To date, machine learning techniques have shown promise in approximating PFM simulations. While standard fracture benchmarks represent realistic scenarios frequently observed in practice, they typically do not provide sufficiently challenging tests for data-driven methods. Here, to address this gap, we introduce a challenging dataset based on PFM simulations designed to benchmark and advance ML methods for fracture modeling. This dataset includes three energy decomposition methods, two boundary conditions, and 1000 random initial crack configurations for a total of 6000 simulations. Each sample contains 100 time steps capturing the temporal evolution of the crack field. Alongside this dataset, we also implement and evaluate Physics Informed Neural Networks (PINN), Fourier Neural Operators (FNO), and UNet models as baselines, and explore the impact of ensembling strategies on prediction accuracy. With this combination of our dataset and baseline models drawn from the literature we aim to provide a standardized and challenging benchmark for evaluating machine learning approaches to solid mechanics. Our results highlight both the promise and limitations of popular current models, and demonstrate the utility of this dataset as a testbed for advancing machine learning in fracture mechanics research.

Benchmark dataset↗
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 OSTI2026

Public Reference Data for Megawatt-Scale Hydrogen Electrolysis - Simulated Marine Hydrokinetic Tidal Turbine

The U.S. Department of Energy and National Laboratory of the Rockies (NLR) demonstrate hydrogen electrolysis, hydrogen compression and storage, and variable hydrogen fuel cell power production using megawatt-scale equipment at NLR’s Flatirons Campus as part of the Advanced Research on Integrated Energy Systems (ARIES) initiative. This dataset is part of that effort and is intended for academic, national laboratory, industrial, and other stakeholders to plan, design, and validate models of megawatt-scale hydrogen technologies and diverse energy infrastructure. These data provide a baseline for how existing hydrogen electrolysis technologies perform when coupled with other energy technologies. This dataset contains inputs and outputs from simulations of a floating marine hydrokinetic turbine over approximately half a tidal cycle (~6.6 hours). Inflow conditions were derived from field measurements in Alaska’s Cook Inlet and represent a tidal environment in which the current speed ramps from near 0 m/s to a peak of 3 m/s and back. The original acoustic doppler current profiler dataset is publicly available on the Marine and Hydrokinetic Data Repository. In a full tidal cycle, the flow reverses and the rotor would reorient; this reversal was not modeled. In the Cook Inlet campaign , turbulence intensity was similar in both directions. Two inflow cases are included. In the first case, labeled “raw” in the files, the measured current time series was used directly in the InflowWind module of OpenFAST. Speed and direction were applied as a function of time and elevation, uniformly in the horizontal direction. With full spatial coherence, this approach captures high turbulent variability and results in pronounced power fluctuations, so it is considered a conservative, near-worst-case representation of loading. In the second case, labeled “average” in the files, a 30-minute moving average was applied to extract the slowly varying mean speed. The residual fluctuations about this mean were used to generate spatially varying, full-field turbulence inputs with TurbSim, giving a more physically realistic representation of the inflow across the rotor disk. Two random realizations were used to produce distinct inflow conditions for two OpenFAST simulations representing a two-turbine array. The same turbulence intensity is applied across the full time series, producing larger fluctuations at the start and end, where the mean speed is low. The second case is the more appropriate framework for performance and power assessment but overpredicts turbulence at lower flow speeds and underpredicts it at higher speeds. As the floating platform moves and the rotor changes its x-position, Taylor’s frozen turbulence hypothesis used by InflowWind assumes a constant rather than a time-varying mean velocity, introducing some inaccuracy in the velocity plane sampling. The turbine modeled is the 500-kW Reference Model 1, a horizontal-axis two-bladed hydrokinetic turbine on a four-column floating semisubmersible substructure . Simulations were performed using OpenFAST v4.1 with the Reference Open Source Controller (ROSCO) v2.10. All input files required to reproduce the simulations are included. The electrolyzer is a 1.25-MW proton exchange membrane type MC250 system manufactured by Nel . This unit supports up to 2.5 MW, but NLR has only a single 1.25-MW stack. The datasets report hydrogen balance-of-plant and system data, all captured at 1 Hz, including hydrogen mass production measured with an Emerson Coriolis flow meter. The system controls hydrogen production by varying direct current applied to the stack, from a maximum of 3,000 A to a minimum safe operating current of 300 A, or 10%. Because the current–voltage characteristic changes as the stack ages and efficiency degrades, the actual minimum safe operating power changes over time. The simulated tidal turbine time series data was translated from power (kilowatts) to current (amperes) using a curve fit with calibration data and sent to the electrolyzer power supply at 1-Hz. Each zip file represents a single tidal electrolysis experiment and is named: {technology}_{inflow method}_{number of 500 kW tidal turbines connected} For instance, “tidal-500kW-RM1_average_2.zip” is a 6-hour experiment using the 500-kW tidal reference model, scaled by 2x (1-MW) to better match the electrolyzer maximum of 1.25MW, fed with the 30-minute moving average current case. Each zip folder contains the following files: A .csv file of raw data. An .xlsx file explaining all the fields in the raw data. A .png plot showing the time series of hydrogen production in kilograms per hour, electrolysis power consumption, and input wave power. A .csv file combines all tidal profiles as "combined_tidal_experiments.csv." A separate experiment, “characterization_200.zip,” shows the MC250 electrolyzer steady-state response with 30-minute load steps over 5 hours and is accessible with this entry.

08 HYDROGEN↗
DOE OSTI2016

Aggregate modeling of fast-acting demand response and control under real-time pricing

This paper develops and assesses the performance of a short-term demand response (DR) model for utility load control with applications to resource planning and control design. Long term response models tend to underestimate short-term demand response when induced by prices. This has two important consequences. First, planning studies tend to undervalue DR and often overlook its benefits in utility demand management program development. Second, when DR is not overlooked, the open-loop DR control gain estimate may be too low. This can result in overuse of load resources, control instability and excessive price volatility. Our objective is therefore to develop a more accurate and better performing short-term demand response model. We construct the model from first principles about the nature of thermostatic load control and show that the resulting formulation corresponds exactly to the Random Utility Model employed in economics to study consumer choice. The model is tested against empirical data collected from field demonstration projects and is shown to perform better than alternative models commonly used to forecast demand in normal operating conditions. The results suggest that (1) existing utility tariffs appear to be inadequate to incentivize demand response, particularly in the presence of high renewables, and (2) existing load control systems run the risk of becoming unstable if utilities close the loop on real-time prices.

Aggregate Modeling↗
DOE OSTI2025

Detecting outbreaks using a spatial latent field

In this paper, we present a method for estimating the infection-rate of a disease as a spatial-temporal field. Our data comprises time-series case-counts of symptomatic patients in various areal units of a region. We extend an epidemiological model, originally designed for a single areal unit, to accommodate multiple units. The field estimation is framed within a Bayesian context, utilizing a parameterized Gaussian random field as a spatial prior. We apply an adaptive Markov chain Monte Carlo method to sample the posterior distribution of the model parameters condition on COVID-19 case-count data from three adjacent counties in New Mexico, USA. Our results suggest that the correlation between epidemiological dynamics in neighboring regions helps regularize estimations in areas with high variance (i.e., poor quality) data. Using the calibrated epidemic model, we forecast the infection-rate over each areal unit and develop a simple anomaly detector to signal new epidemic waves. Our findings show that anomaly detector based on estimated infection-rates outperforms a conventional algorithm that relies solely on case-counts.

Safta, Cosmin [Sandia National Laboratories (SNL-C↗
DOE OSTI2018

Quantitative Boltzmann–Gibbs Principles via Orthogonal Polynomial Duality

We study fluctuation fields of orthogonal polynomials in the context of particle systems with duality. We thereby obtain a systematic orthogonal decomposition of the fluctuation fields of local functions, where the order of every term can be quantified. This implies a quantitative generalization of the Boltzmann–Gibbs principle. In the context of independent random walkers, we complete this program, including also fluctuation fields in non-stationary context (local equilibrium). For other interacting particle systems with duality such as the symmetric exclusion process, similar results can be obtained, under precise conditions on the n particle dynamics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Conditional Karhunen-Loève expansion for uncertainty quantification and active learning in partial differential equation models

We use a conditional Karhunen-Lo` eve (KL) model to quantify and reduce uncertainty in a stochastic partial differential equation (SPDE) problem with partially-known space-dependent coefficient, Y (x). We assume that a small number of Y (x) measurements are available and model Y (x) with a KL expansion. We achieve reduction in uncertainty by conditioning the KL expansion coefficients on measurements. We consider two approaches for conditioning the KL expansion: In Approach 1, we condition the KL model first and then truncate it. In Approach 2, we first truncate the KL expansion and then condition it. We employ the conditional KL expansion together with Monte Carlo and sparse grid collocation methods to compute the moments of the solution of the SPDE problem. Uncertainty of the problem is further reduced by adaptively selecting additional observation locations using two active learning methods. Method 1 minimizes the variance of the PDE coefficient, while Method 2 minimizes the variance of the solution of the PDE. We demonstrate that conditioning leads to dimension reduction of the KL representation of Y (x). For a linear diffusion SPDE with uncertain log-normal coefficient, we show that Approach 1 provides a more accurate approximation of the conditional log-normal coefficient and solution of the SPDE than Approach 2 for the same number of random dimensions in a conditional KL expansion. Furthermore, Approach 2 provides a good estimate for the number of terms of the truncated KL expansion of the conditional field of Approach 1. Finally, we demonstrate that active learning based on Method 2 is more efficient for uncertainty reduction in the SPDE’s states (i.e., it leads to a larger reduction of the variance) than active learning using Method 2.

Conditioned Karhunen-Lo` eve expanion, machine lea↗
DOE OSTI2016

Spatiotemporal modeling of node temperatures in supercomputers

Los Alamos National Laboratory (LANL) is home to many large supercomputing clusters. These clusters require an enormous amount of power (~500-2000 kW each), and most of this energy is converted into heat. Thus, cooling the components of the supercomputer becomes a critical and expensive endeavor. Recently a project was initiated to investigate the effect that changes to the cooling system in a machine room had on three large machines that were housed there. Coupled with this goal was the aim to develop a general good-practice for characterizing the effect of cooling changes and monitoring machine node temperatures in this and other machine rooms. This paper focuses on the statistical approach used to quantify the effect that several cooling changes to the room had on the temperatures of the individual nodes of the computers. The largest cluster in the room has 1,600 nodes that run a variety of jobs during general use. Since extremes temperatures are important, a Normal distribution plus generalized Pareto distribution for the upper tail is used to model the marginal distribution, along with a Gaussian process copula to account for spatio-temporal dependence. A Gaussian Markov random field (GMRF) model is used to model the spatial effects on the node temperatures as the cooling changes take place. This model is then used to assess the condition of the node temperatures after each change to the room. The analysis approach was used to uncover the cause of a problematic episode of overheating nodes on one of the supercomputing clusters. Lastly, this same approach can easily be applied to monitor and investigate cooling systems at other data centers, as well.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Computational materials reliability assessment of hydrogen fueled gas turbine power generation engines

The use of blended fuel sources in land based gas turbine engines drives variations in the resulting operational profile (temperatures and pressures) which can impact engine reliability. Furthermore, variability in the manufacture of components affects the resulting microstructure which directly impacts material performance and reliability. Currently, data-driven models are typically used for maintaining and inspecting fleets of engines. Without explicitly capturing material and operational sources of variability conservatism must be used in developing component-level reliability models. Therefore, there exists an opportunity to use information from materials-scale physics models to better inform reliability modeling and reduce conservatism; the impact is more cost-efficient operation and maintenance of current and future fleets. Specifically, this work establishes a computational framework for evaluating the probabilistic high temperature creep performance of hot-section Ni-based superalloys where uncertainty comes from both microstructural and operational variability. A novel high-fidelity physics model which phenomenologically captures grain-boundary sensitive phenomena has been established. A probabilistic calibration procedure was used to calibrate the model and capture uncertainty in the parameterized model coefficients. A design of experiments methodology was established for identifying informative microstructural digital representations for suitable for forward model evaluation. Results show that training a machine-learning surrogate using this design criteria outperforms random selection of microstructural representations. Finally, two surrogate models were developed: (1) a deterministic surrogate model which predicts the local field response given microstructure, constitutive model parameters, and operating conditions (stress, temperature) and (2) a probabilistic model, where uncertainty comes from constitutive law uncertainty, built using denoising diffusion probabilistic models which samples responses given (1) microstructure and (2) operating conditions. These surrogate models enable partner Siemens Energy to rapidly perform UQ analysis specific to creep deformation across a range of microstructures and operating conditions. The impact is that these ML and physics codes can be used to establish more advanced reliability models for the inspection, servicing, and maintenance of land based gas turbine engines.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Experimental Investigation of Nanosecond and Subnanosecond Pulsed DBD in Atmospheric Air: Fast Imaging and Spectroscopy

Dielectric barrier discharge (DBD), as an easy and simple way of generation of non-thermal plasma, has found a number of applications in variuos fields. However, development of a homogeneous, or uniform, DBD that would operate at atmospheric pressure conditions in atmospheric air, would open a number of new applications in various fields from thin film coatings to plasma medicine. Unfortunately, at atmospheric pressure, a uniform DBD can be easily transformed into a filamentary dielectric DBD; therefore some serious issues arise, such as gas heating due to strong discharges in the random microdischarge channel and non-uniform energy distribution, which adversely affect applications. These issues traditionally are solved by the use of an appropriate working gas composition, an alternating current driving frequency, lowering of gas pressure, etc. The transitions between discharge modes in the same experimental conditions have been thoroughly investigated in nitrogen, rare gases and their mixtures with air and other gases. In many cases (for example, in plasma medicine), these methods, especially those related to gas composition and pressure, may not be applied in a convenient manner. Recent advances in pulsed power technology permitted application of much faster voltage rise times (including the subnanosecond range) and short (few nanoseconds) pulses, and revealed that uniform DBD can, in fact, be generated in atmospheric air. Such discharges are in the great interest for many applications, however to date there still little understanding of the mechanisms of their operation and characteristics. Currently, there is no adequate model of the uniform dielectric barrier discharge development in atmospheric air. Although extensive studies have been performed on understanding of the nature of the pulsed DBD uniformity, until now there is still little understanding of the mechanism of the DBD transition from the filamentary mode to uniform mode. One of the reasons for this is that development of streamers, and later – filaments, occurs on the sub-nanosecond and nanosecond time scales, and therefore requires imaging and other diagnostic techniques with corresponding speed of registration. Recently developed technologies allow such studies. Here we demonstrate that DBD uniformity strongly depends on applied electric field in the discharge gap. More specifically, the discharge uniformity may be achieved in the case when two conditions are satisfied: (1) stong overvoltage in the discharge gap (provided by fast rise times), when anode-directed streamers are formed, and (2) short pulse duration that prevents discharge overheating due to rising conductivity (current) which leads to formation of filaments. We show that by controlling the applied (global) electric field in ns-pulsed DBD, it is possible to control the uniformity of the discharge. In addition, this offers better control of the discharge chemistry due to local changes of electric fields and therefore electron energy distribution function.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

4D ML reconstruction as a tool for volumetric PET-based treatment verification in ion beam radiotherapy

Purpose: An innovative strategy to improve the sensitivity of positron emission tomography (PET)-based treatment verification in ion beam radiotherapy is proposed. Methods: Low counting statistics PET images acquired during or shortly after the treatment (Measured PET) and a Monte Carlo estimate of the same PET images derived from the treatment plan (Expected PET) are considered as two frames of a 4D dataset. A 4D maximum likelihood reconstruction strategy was adapted to iteratively estimate the annihilation events distribution in a reference frame and the deformation motion fields that map it in the Expected PET and Measured PET frames. The outputs generated by the proposed strategy are as follows: (1) an estimate of the Measured PET with an image quality comparable to the Expected PET and (2) an estimate of the motion field mapping Expected PET to Measured PET. The details of the algorithm are presented and the strategy is preliminarily tested on analytically simulated datasets. Results: The algorithm demonstrates (1) robustness against noise, even in the worst conditions where 1.5 × 10{sup 4} true coincidences and a random fraction of 73% are simulated; (2) a proper sensitivity to different kind and grade of mismatches ranging between 1 and 10 mm; (3) robustness against bias due to incorrect washout modeling in the Monte Carlo simulation up to 1/3 of the original signal amplitude; and (4) an ability to describe the mismatch even in presence of complex annihilation distributions such as those induced by two perpendicular superimposed ion fields. Conclusions: The promising results obtained in this work suggest the applicability of the method as a quantification tool for PET-based treatment verification in ion beam radiotherapy. An extensive assessment of the proposed strategy on real treatment verification data is planned.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2025

Thermo-Mechanical Phase-Field Modeling of Fracture in High-Burnup UO2 Fuels Under Transient Conditions

This study presents a novel multiphysics phase-field fracture model to analyze high-burnup uranium dioxide (UO2) fuel behavior under transient reactor conditions. Fracture is treated as a stochastic phase transition, which inherently accounts for the random microstructural effects that lead to variations in the value of fracture strength. Moreover, the model takes into consideration the effects of temperature and burnup on thermal conductivity. Therefore, the model is able to predict crack initiation, propagation, and complex morphologies in response to thermal gradients and stress distributions. Several simulations were conducted to investigate the effects of operational and transient conditions on fracture behavior and the resulting cracking patterns. High-burnup fuels exhibit reduced thermal conductivity, elevating temperature gradients and resulting in extensive radial and circumferential cracks. Transient heating rates and temperatures significantly affect fracture patterns, with higher heating rates generating steeper gradients and more irregular crack trajectories. This approach provides critical insights into fuel integrity during accident scenarios and supports the safety evaluation of extended burnup limits.

Chemistry↗
DOE OSTI2017

Local properties of the large-scale peaks of the CMB temperature

In the present work, we study the largest structures of the CMB temperature measured by Planck in terms of the most prominent peaks on the sky, which, in particular, are located in the southern galactic hemisphere. Besides these large-scale features, the well-known Cold Spot anomaly is included in the analysis. All these peaks would contribute significantly to some of the CMB large-scale anomalies, as the parity and hemispherical asymmetries, the dipole modulation, the alignment between the quadrupole and the octopole, or in the case of the Cold Spot, to the non-Gaussianity of the field. The analysis of the peaks is performed by using their multipolar profiles, which characterize the local shape of the peaks in terms of the discrete Fourier transform of the azimuthal angle. In order to quantify the local anisotropy of the peaks, the distribution of the phases of the multipolar profiles is studied by using the Rayleigh random walk methodology. Finally, a direct analysis of the 2-dimensional field around the peaks is performed in order to take into account the effect of the galactic mask. The results of the analysis conclude that, once the peak amplitude and its first and second order derivatives at the centre are conditioned, the rest of the field is compatible with the standard model. In particular, it is observed that the Cold Spot anomaly is caused by the large value of curvature at the centre.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Considerations for the temperature stratification in a pre-burn constant-volume combustion chamber

In recent years, the Engine Combustion Network (ECN) has developed as a worldwide reference for understanding and describing engine combustion processes, successfully bringing together experimental and numerical efforts. Since experiments and numerical simulations both target the same boundary conditions, an accurate characterization of the stratified environment that is inevitably present in experimental facilities is required. The difference between the core-, and pressure-derived bulk-temperature of pre-burn combustion vessels has been addressed in various previous publications. Additionally, thermocouple measurements have provided initial data on the boundary layer close to the injector nozzle, showing a transition to reduced ambient temperatures. The conditions at the start of fuel injection influence physicochemical properties of a fuel spray, including near nozzle mixing, heat release computations, and combustion parameters. To address the temperature stratification in more detail, thermocouple measurements at larger distances from the spray axis have been conducted. Both the temperature field prior to the pre-combustion event that preconditions the high-temperature, high-pressure ambient, as well as the stratification at the moment of fuel injection were studied. To reveal the cold boundary layer near the injector with a better spatial resolution, Rayleigh scattering experiments and thermocouple measurements at various distances close to the nozzle have been carried out. The impact of the boundary layers and temperature stratification are illustrated and quantified using numerical simulations at Spray A conditions. Next to a reference simulation with a uniform temperature field, six different stratified temperature distributions have been generated. These distributions were based on the mean experimental temperature superimposed by a randomized variance, again derived from the experiments. The results showed that an asymmetric flame structure arises in the computed results when the temperature stratification input is used. In these predictions, first-stage ignition is advanced by 24 $\mu$s , while second-stage ignition is delayed by 11 $\mu$s. At the same time a lift-off length difference between the top and the bottom of up to 1.1 mm is observed. Furthermore, the lift-off length is less stable over time. Given the shown dependency, the temperature data is made available along with the vessel geometry data as a recommended basis for future numerical simulations.

42 ENGINEERING↗
DOE OSTI2018

On the Effect of Small-Angle Scattering by Density Fluctuations on the Efficiency of Linear Transformation of Ordinary and Extraordinary Waves in a Toroidally Inhomogeneous Plasma

We have analyzed the efficiency of tunneling of quasi-optical wave beams through the evanescent region in the vicinity of plasma cutoff in a randomly inhomogeneous magnetoactive plasma. A new theoretical model proposed here makes it possible to study the effect of a random phase modulation induced by density fluctuations in the wave beam path on the efficiency of the linear transformation of waves in the 2D geometry corresponding to the experimental conditions for heating the overdense plasma in toroidal magnetic systems. We have derived a general analytic expression connecting the wave beam coefficient of transformation averaged over the ensemble of random wave field realizations with its phase correlation function. We have analyzed the dependence of the coefficient of transformation on the correlation length for the random phase distribution in the beam and on the traversed path length to the interaction region. The threshold value of the path length above which the fluctuations produce a dominating effect has been determined. The importance of taking into account the 2D effects has been demonstrated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Effect of Mach number and volume fraction in air-shock interacting with a bed of randomly distributed spherical particles

Shock-particle interaction is a fundamental problem in many engineering applications, with the dynamics being heavily influenced by the incident-shock Mach number and the particle volume fraction. In this paper we present fully resolved inviscid simulations of an incident-shock wave traveling through a bed of randomly distributed spherical particles. We vary the strength of the incident shock along with the particle volume fraction in order to study the complex wave interaction during shock-particle interaction. In this study we are interested in the early-time behavior during which the particles do not move and hence in our simulations all the particles are fixed in space. We compute the streamwise average of flow field quantities to generate the x - t contour plots to study the unsteady oscillations inside the particle bed. We observe that the transmitted shock slows down under certain conditions and it is partly due to tortuosity and partly due to weakening caused by energy dissipation. We also present the force histories of the streamwise drag and lift forces for all the particles. The random distribution of particles leads to high variability in the drag force experienced by the particles. We compute the mean peak drag force as a function of the streamwise location to study the mean behavior of the transmitted shock. Based on our findings, we propose simple modifications to improve the current point-particle models used in Euler-Lagrange simulations of shock interacting with a bed of particles.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2015

Influence of cooling rate in planar thermally assisted magnetic random access memory: Improved writeability due to spin-transfer-torque influence

This paper investigates the effect of a controlled cooling rate on magnetic field reversal assisted by spin transfer torque (STT) in thermally assisted magnetic random access memory. By using a gradual linear decrease of the voltage at the end of the write pulse, the STT decays more slowly or at least at the same rate as the temperature. This condition is necessary to make sure that the storage layer magnetization remains in the desired written direction during cooling of the cell. The influence of the write current pulse decay rate was investigated on two exchange biased synthetic ferrimagnet (SyF) electrodes. For a NiFe based electrode, a significant improvement in writing reproducibility was observed using a gradual linear voltage transition. The write error rate decreases by a factor of 10 when increasing the write pulse fall-time from ∼3 ns to 70 ns. For comparison, a second CoFe/NiFe based electrode was also reversed by magnetic field assisted by STT. In this case, no difference between sharp and linear write pulse fall shape was observed. We attribute this observation to the higher thermal stability of the CoFe/NiFe electrode during cooling. In real-time measurements of the magnetization reversal, it was found that Ruderman-Kittel-Kasuya-Yosida (RKKY) coupling in the SyF electrode vanishes for the highest pulse voltages that were used due to the high temperature reached during write. As a result, during the cooling phase, the final state is reached through a spin-flop transition of the SyF storage layer.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Efficient and robust quantum random number generation by photon number detection

We present an efficient and robust quantum random number generator based upon high-rate room temperature photon number detection. We employ an electric field-modulated silicon avalanche photodiode, a type of device particularly suited to high-rate photon number detection with excellent photon number resolution to detect, without an applied dead-time, up to 4 photons from the optical pulses emitted by a laser. By both measuring and modeling the response of the detector to the incident photons, we are able to determine the illumination conditions that achieve an optimal bit rate that we show is robust against variation in the photon flux. We extract random bits from the detected photon numbers with an efficiency of 99% corresponding to 1.97 bits per detected photon number yielding a bit rate of 143 Mbit/s, and verify that the extracted bits pass stringent statistical tests for randomness. Our scheme is highly scalable and has the potential of multi-Gbit/s bit rates.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗