Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Distribution function”

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 1,225 records · Page 68
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Simulating pitch angle scattering using an explicitly solvable energy-conserving algorithm

Particle distribution functions evolving under the Lorentz operator can be simulated with the Langevin equation for pitch angle scattering. This approach is frequently used in particle based Monte-Carlo simulations of plasma collisions, among others. However, most numerical treatments do not guarantee energy conservation, which may lead to unphysical artifacts such as numerical heating and spectra distortions. We present a novel structure-preserving numerical algorithm for the Langevin equation for pitch angle scattering. Similar to the well-known Boris algorithm, the proposed numerical scheme takes advantage of the structure-preserving properties of the Cayley transform when calculating the velocity-space rotations. The resulting algorithm is explicitly solvable, while preserving the norm of velocities down to machine precision. We demonstrate that the method has the same order of numerical convergence as the traditional stochastic Euler-Maruyama method.

Zhang, Xin↗
DOE OSTI2021

Observation of synergy between lower hybrid waves at two frequencies in EAST

Synergistic effects between two frequencies of lower hybrid (LH) waves—operating at 2.45 and 4.6 GHz—were observed in experiment on EAST for the first time. At low density (n_e,lin ≈ 2.0 × 10^19m^−3), simultaneous injection of a 65/35 mix of 2.45 GHz/4.6 GHz power achieved an LHCD efficiency that was 25% higher than what should be expected from the linear combination of the two sources. The experiment was interpreted with time-dependent simulations, using the equilibrium and transport solver TRANSP, coupled with the ray-tracing code GENRAY and the Fokker-Planck solver CQL3D. For each discharge, profiles of current and hard x-ray from simulation and measurement agree within uncertainties. An examination of the electron distribution function indicates that the LH synergy is supported by the increased width of the LH resonance plateau in the simultaneous injection case compared to independent injection.

Lower hybrid synergy↗
DOE OSTI2021

Emu v1.1

Emu is a particle-in-cell code for solving the neutrino quantum kinetic equations in 1, 2, or 3 spatial dimensions with arbitrary angular resolution in order to simulate neutrino flavor transformation in neutron star merger and core collapse supernova environments. Emu represents the neutrino distribution function as a set of particles, each of which represent a collection of neutrinos and antineutrinos with unique position and momentum. Each particle carries two density matrices to define the flavor state of the neutrinos and antineutrinos it represents. Emu includes the vacuum, matter, and neutrino self-interaction potentials. Emu calculates the self-interaction potential using PIC deposition and interpolation algorithms that efficiently compute the local number density and flux at particle locations. Emu is implemented in C++ and is based on the AMReX library for high-performance, block-structured adaptive mesh refinement. Emu is parallelized with MPI + OpenMP for CPUs and MPI + CUDA for GPUs.

Willcox, Donald↗
DOE OSTI2023

GSAS-II (GENERAL STRUCTURE AND ALAYSIS SYSTEM-II)

GSAS-II is a comprehensive data analysis and simulation software package for determining atomistic structures from single-crystal and powder diffraction data. X-ray, neutron and electron data may be used. Constant-wavelength, time-of-flight and pink-beam sources may be employed and analyses may combine different types of data. The package is oriented towards materials characterization problems, but protein crystallography, pair-distribution function analysis, small-angle scattering and reflectometry analysis are all included. GSAS-II provides tools for indexing and solving structures from powder diffraction data, in addition to structure, data and results visualization tools

TOBY, BRIANH↗
DOE OSTI2024

LTAU-FF: Loss Trajectory Analysis for Uncertainty in atomistic Force Fields

LTAU (Loss Trajectory Analysis for Uncertainty) is a technique for estimating the uncertainty in a model's predictions for any regression task by approximating the cumulative distribution function (CDF) of the model's errors over the course of training. The approximated CDF, combined with a similarity search algorithm in the model's descriptor space, can then be used to estimate the likelihood that the model's predictions on any given test point will fall below a chosen threshold. LTAU-FF (LTAU in atomistic Force Fields) is the application of LTAU specifically for use with atomistic force fields.

Vita, JoshuaA↗
DOE OSTI2024

PDF DECODER ANALYSIS CODE

SF-24-038"PDFdecoder", as a new application to explore parametrizations of parton distribution functions (PDFs) of the proton or other hadrons. The PDFs are fundamental quantities in particle physics which are necessary inputs to precise theoretical predictions for experiments at the Large Hadron Collider (LHC) and other facilities. As such, understanding how the PDFs are parametrized and associated uncertainties is a pressing need. The specific problem PDFdecoder confronts is the need of having a tractable and interpretably machine-learning (ML) framework to parametrize the PDFs and their uncertainties so as to understand how a given preferred parametrization is obtained. This problem has not been significantly addressed in the current literature. While other groups have used ML-based approaches to parametrize PDFs in the form of feed-forward neural networks, the question of tractability has not been explored in a PDF context. Our solution makes significant progress in this problem by using an array of encoder-decoder (essentially, autoencoder) architectures with varying constraints to the intermediate latent spaces based on interpretable physics. As a consequence, the trained models can be used as generative networks to produce interpretable predictions for the PDFs in a way that can be refined and studied further.

Hobbs, Timothy↗
DOE OSTI2025

NCSU-iSSRL/WEC-Mobile-HLCD

High level co-design tools for the WEC-Mobile system. These tools use public MFWAM wave forecasts to generate joint distribution functions for likely system deployment sites and calculate feasible system dimensions using a dimensional analysis.

Alam, Minhazul [Univ. of Michigan, Ann Arbor, MI (↗
DOE OSTI2023

Climate Impacts of Convective Cloud Microphysics in NCAR CAM5

Here we improved the treatments of convective cloud microphysics in the NCAR Community Atmosphere Model version 5.3 (CAM5.3) by 1) implementing new terminal velocity parameterizations for convective ice and snow particles, 2) adding graupel microphysics, 3) considering convective snow detrainment, and 4) enhancing rain initiation and generation rate in warm clouds. Furthermore, we evaluated the impacts of improved microphysics on simulated global climate, focusing on simulated cloud radiative forcing, graupel microphysics, convective cloud ice amount, and tropical precipitation. Compared to CAM5.3 with the default convective microphysics, the too-strong cloud shortwave radiative forcing due primarily to excessive convective cloud liquid is largely alleviated over the tropics and midlatitudes after rain initiation and generation rate is enhanced, in better agreement with the CERES-EBAF estimates. Geographic distributions of graupel occurrence are reasonably simulated over continents; whereas the graupel occurrence remains highly uncertain over the oceanic storm-track regions. When evaluated against the CloudSat–CALIPSO estimates, the overestimation of convective ice mass is alleviated with the improved convective ice microphysics, among which adding graupel microphysics and the accompanying increase in hydrometeor fall speed play the most important role. The probability distribution function (PDF) of rainfall intensity is sensitive to warm rain processes in convective clouds, and enhancement in warm rain production shifts the PDF toward heavier precipitation, which agrees better with the TRMM observations. Common biases of overestimating the light rain frequency and underestimating the heavy rain frequency in GCMs are mitigated.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

BSDF Data generation for daylight applications: A call for international standardization

Standardized methods for generating angle-dependent, bidirectional, solar-optical properties for complex fenestration systems do not exist, which means that energy and daylight evaluations in building performance simulations often suffer from major inaccuracies. This position paper provides an overview of state-of-the-art data-driven methods for characterizing light scattering properties of fenestration materials and blind systems (e.g. fabrics, metal slats, patterned glazing), validation via laboratory, simulation and field tests, and salient issues in support of standardization of such methods via the International Standardization Organization (ISO). The ISO standard is intended to provide the fundamental underpinnings for recently mandated daylight standards that rely on bidirectional scattering distribution function data for climate-based daylight modelling and building performance simulations.

Geisler-Moroder, D.↗
DOE OSTI2021

Powder Diffraction Crystallography : 98 Years as Plan B?

A historical essay written in 2014 to accompany a talk in the Transactions Symposium "100 Years of Crystallography" at the 2014 American Crystallographic Association Annual Meeting (talk at https://history.amercrystalassn.org/h-brian-toby-videos). Abstract: Powder diffraction patterns were used for structure solution almost contemporaneously with the Bragg’s initial single-crystal determinations. Structures have continued to be determined from powders, in the subsequent century. There has been a nearly explosive growth in this technique since the 1970’s. This was initially driven by neutron diffraction, which prompted the invention of Rietveld analysis, but synchrotrons have gone on to have at least as large an impact. There is no apparent upper limit to the scale of systems that can be structurally studied with powder diffraction; novel protein structures have been determined. When short-range ordering or disordered regions of a material are of interest, the Fourier transform of a powder diffraction pattern, the pair distribution function can be modeled.

Toby, Brian H.↗
DOE OSTI2023

Generalized Bayesian Framework for Evaluation of Integral Benchmark Experiments

A recently published generalized Bayesian optimization framework has provided a way to retract any or all of the three common assumptions underlying the conventional Generalized Linear Least Squares (GLLS) optimization method based on the concepts introduced in Ref. [2]. These assumptions are: 1. Perfection: The model used for data evaluation and the prior probability distribution function (PDF) of generalized data are perfect. 2. Normality: The prior and posterior PDF are normal. 3. Linearity: The model is linear. In this work we outline how the framework in [1] could be directly adopted for improved evaluation of nuclear criticality integral benchmark experiments (IBEs) by: 1. Removing the first assumption alone by utilizing the concept of imperfections introduced in [1] to enable evaluation in the presence of discrepancies between the data and model or of missing covariance information by a GLLS method that will be seen as a generalization of the conventional GLLS method employed by the TSURFER code, and by 2. Removing the remaining two assumptions by implementing a Markov Chain Monte Carlo method for computation of the posterior PDF in the SAMPLER code, where TSURFER and SAMPLER are the uncertainty quantification (UQ) codes for IBEs in the SCALE code system based on the GLLS and the stochastic method, respectively. The graphic in Figure 1 categorizes the methods discussed in terms of the assumptions that they employ to determine posterior PDFs.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Accelerated microwave-assisted synthesis and in situ X-ray scattering of tungsten-substituted vanadium dioxide (V 1-x W x O 2 )

Vanadium dioxide (VO 2 ) has been widely studied due to its metal-insulator phase transition at 68°C, below which it is a semiconducting monoclinic phase, P2 1 /c, and above it is a metallic tetragonal phase, P4 2 /mnm. Substituting vanadium with transition metals allows transition temperature tunability. An accelerated microwave-assisted synthesis for VO 2 and 5d tungsten substituted VO 2 presented herein decreased synthesis time by three orders of magnitude while maintaining phase purity, particle size, and transition character. Tungsten substitution amount was determined using inductively coupled plasma-optical emission spectroscopy. Differential scanning calorimetry, superconducting quantum interference device measurements, and in-situ heating and cooling experiments monitored through synchrotron x-ray diffraction (XRD) confirmed the transition temperature decreased with increased tungsten substitution. Scanning electron microscopy analyzed through the line-intercept method produced an average particle size of 3 – 5 μm. Average structure and local structure phase purity was determined through Rietveld analysis of synchrotron XRD and least-squares refinement of pair-distribution function data.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Seismic Event Characterization Using Full Moment Tensors on the Hypersphere

Moment tensor solutions provide insights into the deformation that has occurred in the source region of a seismic event and are therefore of great value in identifying different types of seismic sources, such as when monitoring for underground nuclear tests. Despite this utility, inversion of waveforms recorded by seismometers for their full seismic moment tensor is not yet routine, and development of robust methods to classify events based on this information is in its infancy. Here, we assemble an inventory of 1405 full moment tensor solutions that include explosive, earthquake, and collapse events, and investigate the use of anisotropic probability distribution functions on the 5D hypersphere to discriminate between these sources. Using a Bayesian classifier, we obtain optimal success rates of 98.4% across all events and demonstrate that modification of the prior probabilities provides a natural way to alter the balance between not missing desirable events (such as explosions) versus misclassifying large numbers of undesired events (such as earthquakes). The approach is specifically designed to progress from traditional, bipolar event screening metrics to more generalized event identification across multiple types of seismic sources. Despite current databases containing insufficient numbers of events to definitively demonstrate at present, we also find intriguing evidence of subgroupings within individual source populations on the hypersphere, for example, between chemical and nuclear explosions, raising the potential possibility of discriminating between these event types in the future.

Geosciences↗
DOE OSTI2022

Structural analyses of amorphous calcium carbonate before and after removing strontium ions from an aqueous solution

Amorphous calcium carbonate (ACC) is a precursor of crystalline calcium carbonate; hence, its structural information at the atomic level is very important for controlling the morphology of crystalline calcium carbonate. In this study, we attempted to elucidate the process of Sr extraction from aqueous solution by using ACC for the purpose of removal of radioactive Sr from the contaminated water leaked after the Fukushima Daiichi nuclear accident. Here, the pair distribution functions, g(r) obtained by X-ray and neutron diffraction (ND) measurements show that ACC has a structure partially similar to that of monohydrocalcite, suggesting that ACC is transferred to the crystalline calcium carbonate starting from its crystal nucleus. Rietveld analysis of the ND data showed that the ACC that removed Sr was crystallized to calcite. However, the Sr–O coordination analyzed using extended X-ray absorption fine structure (EXAFS) implies that the local environment of O around Sr is similar to that in crystalline calcium carbonate aragonite.

36 MATERIALS SCIENCE↗
DOE OSTI2025

High-resolution analysis of clay minerals and amorphous materials in martian analog environments

Mudrocks and mud-sized sediments (i.e., silt to clay) dominate the surface of Earth and Mars. These fine-grained sediments preserve a rich history of sedimentary processes from source to sink and shed light on ancient climates. However, both the physical and chemical nature of these materials make them difficult to fully characterize with traditional laboratory techniques. Here, we explore a cross-disciplinary and high-resolution approach using synchrotron radiation for X-ray diffraction, pair distribution function analysis, and submicrometer-scale X-ray fluorescence, combined with transmission electron microscopy, to better understand the nanostructure and composition of mud-sized sediments from a glacio-fluvial watershed in southwest Iceland. Our results demonstrate that sediments in the cold and wet climate of Iceland are more altered than previously thought, as evidenced by the identification of kaolinite and mixed-layer kaolinite-smectite. Additionally, sediments are enriched in amorphous materials and nanocrystalline phases, as determined from grain morphologies and compositions consistent with allophane, hisingerite, ferrihydrite, and halloysite. These alteration products are present as intimate mixtures that vary across depositional sites, demonstrating the dynamic nature of the secondary assemblage from source to sink. This work has implications for Mars, where, for example, basalt-sourced sedimentary rocks from Gale crater are abundant in clay minerals and amorphous materials. Finally, this work underpins the importance of using high-resolution techniques, a coordinated methodology, and developing innovative approaches for future planetary sample return missions (e.g., Mars sample return).

36 MATERIALS SCIENCE↗
DOE OSTI2022

Publishing statistical models: Getting the most out of particle physics experiments

The statistical models used to derive the results of experimental analyses are of incredible scientific value and are essential information for analysis preservation and reuse. In this paper, we make the scientific case for systematically publishing the full statistical models and discuss the technical developments that make this practical. By means of a variety of physics cases - including parton distribution functions, Higgs boson measurements, effective field theory interpretations, direct searches for new physics, heavy flavor physics, direct dark matter detection, world averages, and beyond the Standard Model global fits - we illustrate how detailed information on the statistical modelling can enhance the short- and long-term impact of experimental results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Les Houches study on inclusive jet production at NNLO+NNLL

Jet production at the LHC is a powerful probe of QCD, making it ideal for precision tests and determinations of QCD parameters such as parton distribution functions and the strong coupling constant. To make the most of the abundant jet production data collected at the LHC, precise calculations are required. While state-of-the-art calculations reach next-to-next-to-leading order (NNLO) QCD accuracy, a critical assessment of the remaining uncertainties arising from non-perturbative effects and missing higher orders remains crucial for correctly interpreting comparisons between theory and data. Scale variation is nearly always used to determine effects from missing higher orders. In this article, we reassess this method in the context of inclusive jet production by performing NNLO QCD calculations supplemented by small-jet-radius resummation through next-to-next-to-leading-logarithmic accuracy (NNLL). We find that NNLL resummation can have an appreciable impact on the scale uncertainty for inclusive jet cross sections, and, for some scale choices, can lead to sizeable shifts of the central cross section. We conclude that scale variations in fixed-order and resummed calculations can drastically underestimate the impact of higher orders for commonly used jet radius parameters, and that missing higher-order estimates obtained via scale variations should be considered unreliable. Our findings add further evidence to the importance of going beyond scale variations in jet and jet substructure calculations.

Lee, Kyle↗