Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Distribution parameters”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Heating and acceleration of ions with Kappa distribution functions by low‐frequency Alfvén wave

Abstract Heating and acceleration of ions with Kappa distribution functions (with parameter ) in low‐beta plasmas, by a low‐frequency Alfvén wave, is investigated using test‐particle simulations, yielding interesting new results. As long as the Alfvén wave amplitude is sufficiently large, the computed net heating energy of ions becomes independent of the wave frequency and amplitude, always approaching the same value of . The eventual energy of ions is dictated only by the initial ion energy and the ratio of the magnetic field energy density to the plasma density. The heating effect of the Kappa ions increases with . During the heating process, the ions are picked up by the Alfvén wave and pitch angle scattered, forming a quasi‐isotropic spherical shell velocity distribution. The Kappa ions are accelerated in the parallel direction, reaching a bulk flow speed roughly equal to the local Alfvén speed. Higher ‐value in the initial Kappa distribution leads to faster saturation. The above results may explain certain features of the ion heating and acceleration in the solar wind and corona.

Li, Kehua↗
DOE OSTI2020

A Bayesian Approach for Statistical–Physical Bulk Parameterization of Rain Microphysics. Part I: Scheme Description

A novel framework is proposed for the bulk parameterization of rain microphysics: the Bayesian Observationally Constrained Statistical–Physical Scheme (BOSS). It is designed to facilitate direct constraint by observations using Bayesian inference. BOSS combines existing process-level microphysical knowledge with flexible process rate formulations and parameters constrained by observations within a Bayesian framework. Furthermore, using a raindrop size distribution (DSD) normalization method that relates DSD moments to one another via generalized power series, generalized multivariate power expressions are derived for the microphysical process rates as functions of a set of prognostic DSD moments. The scheme is flexible and can utilize any number and combination of prognostic moments and any number of terms in the process rate formulations. This means that both uncertainty in parameter values and structural uncertainty associated with the process rate formulations can be investigated systematically, which is not possible using traditional schemes. In this paper, BOSS is compared to two- and three-moment versions of a traditional bulk rain microphysics scheme (denoted as MORR). It is demonstrated that some process formulations in MORR are analytically equivalent to the generalized power expressions in BOSS using one or two terms, while others are not. BOSS is able to replicate the behavior of MORR in idealized one-dimensional rainshaft tests, but with a much more flexible and systematic design. Part II of this study describes the application of BOSS to derive rain microphysical process rates and posterior parameter distributions in Bayesian experiments using Markov chain Monte Carlo sampling constrained by synthetic observations.

58 GEOSCIENCES↗
DOE OSTI2024

Impact of Pressure Distribution and Magnitude on the Performance of Lithium Metal Anodes

Li metal anodes are a critical battery technology due to their ability to substantially increase the energy density of Li-based batteries. It is well known that pressure greatly impacts the performance of a Li-metal anode. However, precisely how the pressure value and distribution of pressure affect performance is unclear. Furthermore, the solid-electrolyte interphase composition that forms under varying pressure distributions remains a key parameter for practical lithium metal anodes. In this work, different pressure distributions were employed by using differently shaped and oriented mechanical springs in the coin cells, resulting in varying contact points. Pressure-sensitive films were used to spatially map the pressure and correlate it to the performance. It was found that higher average pressure does not necessarily have a positive effect on performance. When high pressure is paired with poor pressure uniformity, the performance is in fact worse likely due to the current focusing effect, rendering unsatisfied cycling stability. This work points to the importance of controlling the relationship between average pressure and pressure uniformity.

25 ENERGY STORAGE↗
DOE OSTI2021

Characterizing patterns and variability of building electric load profiles in time and frequency domains

The rapid development of advanced metering infrastructure provides a new data source—building electrical load profiles with high temporal resolution. Electric load profile characterization can generate useful information to enhance building energy modeling and provide metrics to represent patterns and variability of load profiles. Such characterizations can be used to identify changes to building electricity demand due to operations or faulty equipment and controls. In this study, we proposed a two-path approach to analyze high temporal resolution building electrical load profiles: (1) time-domain analysis and (2) frequency-domain analysis. Furthermore, the commonly adopted time-domain analysis can extract and quantify the distribution of key parameters characterizing load shape such as peak-base load ratio and morning rise time, while a frequency-domain analysis can identify major periodic fluctuations and quantify load variability. We implemented and evaluated both paths using whole-year 15-minute interval smart meter data of 188 commercial office building in Northern California. The results from these two paths are consistent with each other and complementary to represent full dynamics of load profiles. The time- and frequency-domain analyses can be used to enhance building energy modeling by: (1) providing more realistic assumptions about building operation schedules, and (2) validating the simulated electric load profiles using the developed variability metrics against the real building load data.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Personalized and uncertainty-aware coronary hemodynamics simulations: From Bayesian estimation to improved multi-fidelity uncertainty quantification

Non-invasive simulations of coronary hemodynamics have improved clinical risk stratification and treatment outcomes for coronary artery disease, compared to relying on anatomical imaging alone. However, simulations typically use empirical approaches to distribute total coronary flow amongst the arteries in the coronary tree, which ignores patient variability, the presence of disease, and other clinical factors. Further, uncertainty in the clinical data often remains unaccounted for in the modeling pipeline. We present an end-to-end uncertainty-aware pipeline to (1) personalize coronary flow simulations by incorporating vessel-specific coronary flows as well as cardiac function; and (2) predict clinical and biomechanical quantities of interest with improved precision, while accounting for uncertainty in the clinical data. We assimilate patient-specific measurements of myocardial blood flow from clinical CT myocardial perfusion imaging to estimate branch-specific coronary artery flows. Simulated noise in the clinical data is used to estimate the joint posterior distributions of the model parameters using adaptive Markov Chain Monte Carlo sampling. Additionally, the posterior predictive distribution for the relevant quantities of interest is determined using a new approach combining multi-fidelity Monte Carlo estimation with non-linear, data-driven dimensionality reduction. This leads to improved correlations between high- and low-fidelity model outputs. Our framework accurately recapitulates clinically measured cardiac function as well as branch-specific coronary flows under measurement noise uncertainty. We observe substantial reductions in confidence intervals for estimated quantities of interest compared to single-fidelity Monte Carlo estimation and state-of-the-art multi-fidelity Monte Carlo methods. This holds especially true for quantities of interest that showed limited correlation between the low- and high-fidelity model predictions. In addition, the proposed multi-fidelity Monte Carlo estimators are significantly cheaper to compute than traditional estimators, under a specified confidence level or variance. The proposed pipeline for personalized and uncertainty-aware predictions of coronary hemodynamics is based on routine clinical measurements and recently developed techniques for CT myocardial perfusion imaging. The proposed pipeline offers significant improvements in precision and reduction in computational cost.

Bayesian parameter estimation↗
DOE OSTI2023

Plasma‐catalytic synthesis of acrylonitrile from methane and nitrogen

Here, in this work, we demonstrate plasma-catalytic synthesis of hydrogen and acrylonitrile (AN) from CH 4 and N 2 . The process involves two steps: (1) plasma synthesis of C 2 H 2 and HCN in a nominally 1:1 stoichiometric ratio with high yield up to 90% and (2) downstream thermocatalytic reaction of these intermediates to make AN. The effect of process parameters on product distributions and specific energy requirements are reported. If the catalytic conversion of C 2 H 2 and HCN in the downstream thermocatalytic step to AN were perfect, which will require further improvements in the thermocatalytic reactor, then at the maximum output of our 1 kW radiofrequency 13.56 MHz transformer, a specific energy requirement of 73 kWh kg AN −1 was determined. The expectation is that scaling up the process to higher throughputs would result in decreases in specific energy requirement into the predicted economically viable range less than 10 kWh kg AN −1 .

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

Coupled Thermo-Hydro-Mechanical Processes in Fractured Rocks: Some Past Scientific Highlights and Future Research Directions

Abstract Coupled thermo-hydro-mechanical (THM) processes in fractured rocks have been a topic of intense scientific research for more than 30 years. The present paper takes a look into the past and highlights some scientific advances which are of an unusual “out-of-the-box” nature, and then looks forward and discusses possible directions of future research in this interesting field of study. Concerning future research directions, we see a trend from a focus on coupled THM processes in single fractures or a few interacting fractures, to the study of coupled THM behavior in complex fracture network systems where the fractures act collectively giving rise to local stress concentration points and points of large pressure gradients. Three examples of future research directions are presented. First is an effort towards identifying characterizing parameters of a fracture network that play a direct controlling role in major coupled THM phenomena (such as induced seismicity and flow channeling), rather than parameters of stochastic distributions of fractures in the network. The second example of research direction is accounting for the heterogeneity and hierarchy of fractures in a fault or fracture zone which has been associated with major THM events in a number of geo-energy projects. The third example is at the opposite end of the first; here it is recognized that in some cases, the coupled THM processes in fractured rocks may be controlled dominantly by only a few key bridges. Identification, characterization, and evaluation of these key bridges should be one of the important research directions in the coming days.

58 GEOSCIENCES↗
DOE OSTI2023

Impact of artificial topological changes on flow and transport through fractured media due to mesh resolution

Abstract We performed a set of numerical simulations to characterize the interplay of fracture network topology, upscaling, and mesh refinement on flow and transport properties in fractured porous media. We generated a set of generic three-dimensional discrete fracture networks at various densities, where the radii of the fractures were sampled from a truncated power-law distribution, and whose parameters were loosely based on field site characterizations. We also considered five network densities, which were defined using a dimensionless version of density based on percolation theory. Once the networks were generated, we upscaled them into a single continuum model using the upscaled discrete fracture matrix model presented by Sweeney et al. (2019). We considered steady, isothermal pressure-driven flow through each domain and then simulated conservative, decaying, and adsorbing tracers using a pulse injection into the domain. For each simulation, we calculated the effective permeability and solute breakthrough curves as quantities of interest to compare between network realizations. We found that selecting a mesh resolution such that the global topology of the upscaled mesh matches the fracture network is essential. If the upscaled mesh has a connected pathway of fracture (higher permeability) cells but the fracture network does not, then the estimates for effective permeability and solute breakthrough will be incorrect. False connections cannot be eliminated entirely, but they can be managed by choosing appropriate mesh resolution and refinement for a given network. Adopting octree meshing to obtain sufficient levels of refinement leads to fewer computational cells (up to a 90% reduction in overall cell count) when compared to using a uniform resolution grid and can result in a more accurate continuum representation of the true fracture network.

58 GEOSCIENCES↗
DOE OSTI2020

Modeling strength and failure variability due to porosity in additively manufactured metals

To model and quantify the variability in plasticity and failure of additively manufactured metals due to imperfections in their microstructure, we have developed uncertainty quantification methodology based on pseudo marginal likelihood and embedded variability techniques. We account for both the porosity resolvable in computed tomography scans of the initial material and the sub-threshold distribution of voids through a physically motivated model. Calibration of the model indicates that the sub-threshold population of defects dominates the yield and failure response. Finally, the technique also allows us to quantify the distribution of material parameters connected to microstructural variability created by the manufacturing process, and, thereby, make assessments of material quality and process control.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Generating realistic building electrical load profiles through the Generative Adversarial Network (GAN)

Building electrical load profiles can improve understanding of building energy efficiency, demand flexibility, and building-grid interactions. Current approaches to generating load profiles are time-consuming and not capable of reflecting the dynamic and stochastic behaviors of real buildings; some approaches also trigger data privacy concerns. In this study, we proposed a novel approach for generating realistic electrical load profiles of buildings through the Generative Adversarial Network (GAN), a machine learning technique that is capable of revealing an unknown probability distribution purely from data. The proposed approach has three main steps: (1) normalizing the daily 24-hour load profiles, (2) clustering the daily load profiles with the k-means algorithm, and (3) using GAN to generate daily load profiles for each cluster. The approach was tested with an open-source database – the Building Data Genome Project. We validated the proposed method by comparing the mean, standard deviation, and distribution of key parameters of the generated load profiles with those of the real ones. The KL divergence of the generated and real load profiles are within 0.3 for majority of parameters and clusters. Additionally, results showed the load profiles generated by GAN can capture not only the general trend but also the random variations of the actual electrical loads in buildings. We report the proposed GAN approach can be used to generate building electrical load profiles, verify other load profile generation models, detect changes to load profiles, and more importantly, anonymize smart meter data for sharing, to support research and applications of grid-interactive efficient buildings.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Detailed chemistry modeling of rotating detonations with dilute n -heptane sprays and preheated air

Utilization of liquid fuels is crucial to enabling commercialization of rotating detonation engines (RDE) in the near future. In this study, Eulerian-Lagrangian simulations are conducted for rotating detonative combustion with dilute n-heptane sprays and preheated air. Two-dimensional flattened configuration is used and a skeletal chemical mechanism with 44 species and 112 elementary reactions for n-heptane combustion is adopted. The flow structure, droplet distribution, and thermochemical parameters in the refill zone are first analyzed. It is shown that the mixture in the refill zone is heterogeneous, including evaporating droplets, vapor, and air. When the total temperature is below 950 K, the average equivalence ratio increases with the total temperature. When it is higher than 950 K, the average equivalence ratio is almost constant. Subsequently, the chemical explosive mode analysis is applied to identify the controlling reactions and dominant combustion modes in the fuel refill zone and reaction fronts. Results demonstrate that the initiation reaction (R104: n-C 7 H 16 + O 2 2-C 7 H 15 + HO 2 ) and low-temperature reaction (R107: RO 2 R'O 2 H) are dominant in the upstream and downstream of the refill zone, respectively. The intermediate species from low-temperature chemistry, R'O 2 H, is found to be important for the chemical explosive mode in the undetonated mixture. The influence of species diffusion and dispersed droplets is further analyzed. Results show that vapor autoignition facilitated by droplet evaporation occurs in the refill zone. Finally, the effects of the air total temperature on the detonation propagation speed and RDE propulsion performance are investigated. Further, it is found that the detonation propagation speed and specific impulse increase with air total temperature. The total pressure ratio first increases and then decreases as the air total temperature increases. Moreover, when the total temperature of the preheated air is above 1,300 K, the effects of low-temperature chemistry are negligible.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2025

Accelerating Embedding Potential Optimization by Reconstructing the Pseudo-Valence Electron Density

Density functional embedding theory (DFET) enables use of electronic structure methods with higher accuracy than density functional theory in a local region, with applications thus far ranging from (photo/electro)catalysis to reactions in solution. DFET partitions a large collection of atoms into smaller groups that interact via a shared embedding (interaction) potential V emb , determined via functional optimization. The optimized effective potential (OEP) process used to optimize V emb is time-consuming and becomes a computational bottleneck due to sharp, oscillating features of V emb near nuclei. Here, similar to pseudopotential theory, by reconstructing electron densities used in the OEP process from smoother pseudo-valence-only (PVO) electron densities as proxies for total densities of the full system and subsystems, we can retain accuracy in the embedded electronic structure calculations while potentially reducing the overhead of V emb construction, within the projector augmented-wave (PAW) formalism. We explore three different chemical reactions as exemplars to test PVO–DFET, namely, H 2 dissociative adsorption on a Cu(111) surface, H 2 O adsorption on a Pt(111) surface, and aqueous [Ca 2+ –SO 4 2– ] ion-pair formation. The PVO approximation works well for all three systems with minimal loss of accuracy (∼10–70 meV error relative to the original exact-derivative (ED) approach) while accelerating V emb generation for the Cu and Pt systems respectively by 20× and 5×. Given proper numerical convergence parameters, the spatial distributions of differences between PVO- and ED-based V emb outside the core regions are small, explaining the exceptional agreement between the two approaches. Finally, we anticipate that this more efficient PVO–DFET approximation will be useful whenever computation of V emb is much more expensive than subsequent embedded high-level electron correlation calculations.

approximation↗
DOE OSTI2024

Anion and Cation Migration at 2D/3D Halide Perovskite Interfaces

Here, this study explores the ionic dynamics in 2D/3Dperovskite solar cells, which are known for their improved efficiency and stability. The focus is on the impact of halide choice in 3D perovskites treated with phenethylammonium halide salts (PEAX, X = Br and I). Our findings reveal that light and heat drive ionic migration in these structures, with PEA + species diffusing into the 3D film in PEABr-treated samples. Mixed-halide 3D perovskites show halide interdiffusion, with bromine migrating to the surface and iodine diffusing into the film. Cathodoluminescence microscopy reveals localized 2Dphases on the 3D perovskite, which become more evenly distributed after thermal treatment. Both PEAX salts enhance the performance of photovoltaic devices. This improvement is attributed to the passivation capabilities of the salts themselves and their respective Ruddlesden–Popper (RP) phases. Annealed PEAI-treated devices show a better balance between efficiency and statistical distribution of photovoltaic parameters.

25 ENERGY STORAGE↗
DOE OSTI2025

A Probabilistic Model for Global EMIC Wave Activity Using Van Allen Probes Observations

Electromagnetic ion cyclotron (EMIC) waves play a key role in radiation belt dynamics through resonant interactions. However, their low occurrence probability, high variability, and spatial intermittency pose challenges for accurate modeling. In this study, we present a machine learning (ML)-based global EMIC wave model built on the entire data set from the Van Allen Probes mission. To capture the distinct statistical characteristics of wave occurrence and amplitude, the model is separated into two modules: an occurrence model trained using ML techniques, and a wave amplitude model sampled from observed probability distributions. The input parameters are limited to real-time or predictable variables to ensure practical applicability. Our model shows strong performance across the entire test set and demonstrates improved predictive capability over a baseline random occurrence model, particularly during quiet geomagnetic conditions. Evaluation during both quiet and active periods confirms the model's ability to represent the clustered and intermittent nature of EMIC wave activity. Furthermore, the model provides global estimates of wave power, enabling integration with radiation belt electron data and showing signatures consistent with wave-induced scattering. We found a good correlation between the global wave activity from the model and relativistic electron observation by Van Allen Probes, regardless of the availability of in situ wave observations. The modular structure of the model also allows for straightforward expansion for additional wave properties, such as wave frequency, which can be modeled independently. This flexible, event-sensitive approach offers a promising framework for data-driven radiation belt simulations and space weather applications.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Ultrafast terahertz emission from emerging symmetry-broken materials

Abstract Nonlinear optical spectroscopies are powerful tools for investigating both static material properties and light-induced dynamics. Terahertz (THz) emission spectroscopy has emerged in the past several decades as a versatile method for directly tracking the ultrafast evolution of physical properties, quasiparticle distributions, and order parameters within bulk materials and nanoscale interfaces. Ultrafast optically-induced THz radiation is often analyzed mechanistically in terms of relative contributions from nonlinear polarization, magnetization, and various transient free charge currents. While this offers material-specific insights, more fundamental symmetry considerations enable the generalization of measured nonlinear tensors to much broader classes of systems. We thus frame the present discussion in terms of underlying broken symmetries, which enable THz emission by defining a system directionality in space and/or time, as well as more detailed point group symmetries that determine the nonlinear response tensors. Within this framework, we survey a selection of recent studies that utilize THz emission spectroscopy to uncover basic properties and complex behaviors of emerging materials, including strongly correlated, magnetic, multiferroic, and topological systems. We then turn to low-dimensional systems to explore the role of designer nanoscale structuring and corresponding symmetries that enable or enhance THz emission. This serves as a promising route for probing nanoscale physics and ultrafast light-matter interactions, as well as facilitating advances in integrated THz systems. Furthermore, the interplay between intrinsic and extrinsic material symmetries, in addition to hybrid structuring, may stimulate the discovery of exotic properties and phenomena beyond existing material paradigms.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Uncertainty quantification in multivariable regression for material property prediction with Bayesian neural networks

With the increased use of data-driven approaches and machine learning-based methods in material science, the importance of reliable uncertainty quantification (UQ) of the predicted variables for informed decision-making cannot be overstated. UQ in material property prediction poses unique challenges, including multi-scale and multi-physics nature of materials, intricate interactions between numerous factors, limited availability of large curated datasets, etc. In this work, we introduce a physics-informed Bayesian Neural Networks (BNNs) approach for UQ, which integrates knowledge from governing laws in materials to guide the models toward physically consistent predictions. To evaluate the approach, we present case studies for predicting the creep rupture life of steel alloys. Experimental validation with three datasets of creep tests demonstrates that this method produces point predictions and uncertainty estimations that are competitive or exceed the performance of conventional UQ methods such as Gaussian Process Regression. Additionally, we evaluate the suitability of employing UQ in an active learning scenario and report competitive performance. The most promising framework for creep life prediction is BNNs based on Markov Chain Monte Carlo approximation of the posterior distribution of network parameters, as it provided more reliable results in comparison to BNNs based on variational inference approximation or related NNs with probabilistic outputs.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Skyrmion lattice formation and destruction mechanisms probed with TR-SANS

Magnetic skyrmions are topologically protected, nanoscale whirls of the spin configuration that tend to form hexagonally ordered arrays. As a topologically non-trivial structure, the nucleation and annihilation of the skyrmion, as well as the interaction between skyrmions, varies from conventional magnetic systems. Recent works have suggested that the ordering kinetics in these materials occur over millisecond or longer timescales, which is unusually slow for magnetic dynamics. The current work investigates the skyrmion ordering kinetics, particularly during lattice formation and destruction, using time-resolved small angle neutron scattering (TR-SANS). Evaluating the time-resolved structure and intensity of the neutron diffraction pattern reveals the evolving real-space structure of the skyrmion lattice and the timeframe of the formation. Measurements were performed on three prototypical skyrmion materials: MnSi, (Fe,Co)Si, and Cu 2 OSeO 3 . To probe lattice formation and destruction kinetics, the systems were prepared in the stable skyrmion state, and then a square-wave magnetic field modulation was applied. The measurements show that the skyrmions quickly form ordered domains, with a significant distribution in lattice parameters, which then converge to the final structure; the results confirm the slow kinetics, with formation times between 10 ms and 99 ms. Comparisons are made between the measured formation times and the fundamental material properties, suggesting the ordering temperature, saturation magnetization and magnetocrystalline anisotropy may be driving the timeframes. Micromagnetic simulations were also performed and support a scaling of the kinetics with sample volume, a behavior which is caused by the reconciling of misaligned domains.

Liyanage, W. L. N. C.↗
DOE OSTI2023

In-situ study of rules of nanostructure evolution, severe plastic deformations, and friction under high pressure

Severe plastic deformations under high pressure are used to produce nanostructured materials but were studied ex-situ. Rough diamond anvils are introduced to reach maximum friction equal to yield strength in shear and the first in-situ study of the evolution of the pressure-dependent yield strength and radial distribution of nanostructural parameters are performed for severely pre-deformed Zr. ω-Zr behaves like perfectly plastic, isotropic, and strain-path-independent and reaches steady values of the crystallite size and dislocation density, which are pressure-, strain- and strain-path-independent. However, steady states for a-Zr obtained with smooth and rough anvils are different, causing major challenge in plasticity theory. In-situ study of severe plastic deformation of α-Zr with rough diamond anvils revealed that pressure-dependent yield strength, crystallite size, and dislocation density are getting steady and plastic strain- and strain-path-independent.

36 MATERIALS SCIENCE↗