Search NASA⌕ Search

SEARCH · Search NASA

Results for “BOLTZMANN DISTRIBUTION”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

A velocity space hybridization-based Boltzmann equation solver

In the present research, a new method for simulation of rarefied gas flows is proposed, a velocity-space hybrid of both a DSMC representation of particles and a discrete velocity quasi-particle representation of the distribution function. The hybridization scheme is discussed in detail, and is numerically verified for two test-cases: the BKW relaxation problem and a stationary Maxwellian distribution. It is demonstrated that such a velocity-space hybridization can provide computational benefits when compared to a pure discrete velocity method or pure DSMC approach, while retaining some of the more attractive properties of discrete velocity methods. Further possible improvements to the velocity-space hybrid approach are discussed.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1981

Direct simulation Monte Carlo technique for modeling of the environment in the vicinity of the Space Shuttle Orbiter

A program for analyzing the flowfield parameters in the neighborhood of the Space Shuttle Orbiter has been developed. The program uses the direct simulation Monte Carlo method, which is a completely probabilistic Monte Carlo technique capable of analyzing 3-dimensional steady or unsteady flow with prescribed internal and external boundary conditions. The freestream flux densities incident on the external flowfield boundaries are calculated from the drifting Maxwellian gas properties of the freestream. The flux entering the flowfield from the internal boundary is calculated from the outgassing flux density distribution over the Shuttle external surface and the prescribed discrete source fluxes. This technique produces a numerical flowfield solution which is the probabilistic equivalent of a complete solution of the time-dependent, 3-dimensional Boltzmann equation. Flowfield results are presented for the following configurations: (1) Shuttle angles of attack of 0 and 90 deg, (2) freestream density 10 to the 9th - 10 to the 11th per cu cm, (3) Shuttle outgassing flux density, and (4) operation of the aft downfiring vernier Reaction Control System (RCS) engine. Results are presented for column density distribution of outgassed and engine species and for the flux density of outgassed and engine species incident on the Shuttle bay and the windshield.

Melfi, L. T., Jr.↗
DOE OSTI2020

Fast collective neutrino oscillations inside the neutrino sphere in core-collapse supernovae

Neutrinos are believed to have a key role in the explosion mechanism of core-collapse supernovae as they carry most of the energy released by the gravitational collapse of a massive star. If their flavor is converted fast inside the neutrino sphere, the supernova explosion may be influenced. This paper is reporting the results of the extended work of our previous paper. Here, we perform a thorough survey of the electron lepton number (ELN) crossing in one of our self-consistent, realistic Boltzmann simulations in two spatial dimensions under axisymmetry for the existence of the crossings between $ν_e$ and $\bar{ν}_e$ angular distributions, or the ELN crossing. We report for the first time the positive detections deep inside the core of the massive star in the vicinity of neutrino sphere at $\textit{r}$ ≈ 16 – 21 km . We find that low values of the electron fraction $Y_e$ produced by convective motions together with the appearance of light elements are critically important to give rise to the ELN crossing by enhancing the chemical potential difference between proton and neutron, and hence by mitigating the Fermi-degeneracy of $ν_e$. Since the region of positive detection are sustained and, in fact, expanding with time, it may have an impact on the explosion of core-collapse supernovae, observational neutrino astronomy, and nucleosynthesis of heavy nuclei.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Benchmark of the KGMf with a coupled Boltzmann equation solver

The Kinetic Global Model framework (KGMf) is an open-source general-purpose global model (spatially averaged) simulation code developed to explore the reaction kinetics and pathways in plasma discharge systems. It contains species continuity and electron energy balance equations, with a time-dependent evaluated electron energy distribution function (EEDF) for electron impact reactions. The EEDF is utilized to determine the rate coefficients for electron impact reactions, which can have profound impact on the temporal evolutions of plasma parameters. Previously, the EEDF was commonly assumed as Maxwellian or an analytical function of the effective electron temperature. In this work, the KGMf is coupled with a Boltzmann equation (BE) solver to self-consistently compute the EEDF. The EEDF evolution frequency is determined based on relative changes of the reduced electric field. The KGMf is benchmarked with the ZDPlasKin code based on high-pressure low-temperature argon plasma discharge cases. Additionally, the temporal evolutions of reduced electric field, electron temperature, EEDF, reaction rates, and species densities, are obtained and compared under different discharge conditions, showing good agreement between the KGMf and the ZDPlasKin simulations. The application of the KGMf for predicting breakdown times in high power microwave discharges is also presented, which shows qualitative agreement with particle-in-cell simulations. The KGMf can be further applied for more complicated plasma discharge systems, where the reaction kinetics are more intricate (e.g., plasma-assisted combustion systems).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Temperature-dependent phonon lifetimes and thermal conductivity of silicon by inelastic neutron scattering and ab initio calculations

In this paper, inelastic neutron scattering on a single crystal of silicon was performed at temperatures from 100 to 1500 K. These experimental data were reduced to obtain phonon spectral intensity at all wave vectors $\overset{→}{Q}$ and frequencies ω in the first Brillouin zone. Thermal broadenings of the phonon peaks were obtained by fitting and by calculating with an iterative ab initio method that uses thermal atom displacements on an ensemble of superlattices. Agreement between the calculated and experimental broadenings was good, with possible discrepancies at the highest temperatures. Distributions of phonon widths versus phonon energy had similar shapes for computation and experiment. These distributions grew with temperature but maintained similar shapes. Parameters from the ab initio calculations were used to obtain the thermal conductivity from the Boltzmann transport equation, which was in good agreement with experimental data. Despite the high group velocities of longitudinal acoustic phonons, their shorter lifetimes reduced their contribution to the thermal conductivity, which was dominated by transverse acoustic modes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2002

Force Evaluation in the Lattice Boltzmann Method Involving Curved Geometry

The present work investigates two approaches for force evaluation in the lattice Boltzmann equation: the momentum- exchange method and the stress-integration method on the surface of a body. The boundary condition for the particle distribution functions on curved geometries is handled with second order accuracy based on our recent works. The stress-integration method is computationally laborious for two-dimensional flows and in general difficult to implement for three-dimensional flows, while the momentum-exchange method is reliable, accurate, and easy to implement for both two-dimensional and three-dimensional flows. Several test cases are selected to evaluate the present methods, including: (i) two-dimensional pressure-driven channel flow; (ii) two-dimensional uniform flow past a column of cylinders; (iii) two-dimensional flow past a cylinder asymmetrically placed in a channel (with vortex shedding); (iv) three-dimensional pressure-driven flow in a circular pipe; and (v) three-dimensional flow past a sphere. The drag evaluated by using the momentum-exchange method agrees well with the exact or other published results.

Mei, Renwei↗
DOE OSTI2020

Interplaying mechanisms behind inclusive jet R and extraction of jet energy loss distributions

The observed inclusive jet suppression in heavy-ion collisions at LHC has a very weak p T dependence over a large range of p T =50–1000 GeV and is almost independent of the colliding energy, though the initial energy density of the bulk medium has increased from √s=2.76 to 5.02 TeV by about 20%. This interesting phenomenon is investigated in the linear Boltzmann transport (LBT) model for jet propagation in an event-by-event 3+1D hydro background. We show that the p T dependence of jet R AA is determined by the initial spectrum in p+p collisions and p T dependence of jet energy loss. Furthermore, jet energy loss distributions for inclusive jet and γ-jet at both LHC energies are extracted directly from experimental data through the state-of-art Bayesian analysis. Additionally, the averaged jet energy loss has a weak p T dependence and the scaled jet energy loss distributions have a large width, both of which are consistent with the LBT simulations and indicate that jet quenching is caused by only a few out-of-cone jet medium scatterings.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Velocity-Space Hybridization of Direct Simulation Monte Carlo and a Quasi-Particle Boltzmann Solver

This paper presents a new method for modeling rarefied gas flows based on hybridization of direct simulation Monte Carlo (DSMC) and discrete velocity method (DVM)-based quasi-particle representations of the velocity distribution function. It is aimed at improving the resolution of the tails of the distribution function (compared with DSMC) and computational efficiency (compared with DVM). Details of the method, such as the collision algorithm and the particle merging scheme, are discussed. The hybrid approach is applied to the study of noise in a Maxwellian distribution, computation of electron-impact ionization rate coefficient, as well as numerical simulation of a supersonic Couette flow. The hybrid-based solver is compared with pure DSMC and DVM approaches in terms of accuracy, computational speed, and memory use. It is shown that such a hybrid approach can provide a lower computational cost than a pure DVM approach, while being able to retain accuracy in modeling high-velocity tails of the distribution function. For problems where trace species have a significant impact on the flow physics, the proposed method is shown to be capable of providing better computational efficiency and accuracy compared with standard fixed-weight DSMC.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1972

Perturbation of the Maxwellian distribution in chemically reacting gas mixtures

Chemically reacting gas mixtures are considered for the case in which inelastic collisions occur significantly less frequently than elastic ones. The first approximation of the general solution of the Boltzmann equation is developed and is shown to contain an additional scaler term. The concentration, the mean mass velocity, and the temperature of the mixture are determined and the rate of reaction in the first approximation is formulated. Tables are presented which give reaction rate values obtained by both the first and second approximation for different activation energies. The tables show that the first approximation is sufficient for practical purposes for calculating the reaction rate.

Alekseyev, B. V.↗
DOE OSTI2021

Imprints of high-momentum nucleons in nuclei on hard photons from heavy-ion collisions near the Fermi energy

The short-range correlation (SRC) induced by the tensor force in the isosinglet neutron-proton interaction channel leads to a high-momentum tail (HMT) in the single-nucleon momentum distributions n(k) in nuclei. Owing to the remaining uncertainties about the tensor force, the shape of the nucleon HMT may be significantly different from the dilute interacting Fermi gas model prediction n(k) ~ 1/k 4 similar to the HMT in cold atoms near the unitary limit. Within an isospin- and momentum-dependent Boltzmann-Uehling-Uhlenbeck transport model incorporating approximately the nucleon HMT, we investigate hard photon emissions in 14 N + 12 C and 48 Ca + 124 Sn reactions at beam energies around the Fermi energy. Imprints of different shapes of the HMT on the energy spectrum, angular distribution and transverse momentum spectrum of hard photons are studied. Furthermore, while the angular distribution does not carry any information about the shape of the nucleon HMT, the energy spectra and especially the mid-rapidity transverse momentum spectra of hard photons are found to bare strong imprints of the shapes of nucleon HMTs in the two colliding nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Traces of nonequilibrium effects, initial condition, bulk dynamics, and elementary collisions in the charm observables

Heavy quarks produced in relativistic heavy-ion collisions are known to be sensitive probes of the hot and dense QCD matter they traverse. In this paper we study how their dynamics is affected by the nature of the bulk evolution of the QCD matter, the initial condition of the system, and the treatment of elementary interactions between heavy quarks and the surrounding medium. For the same initial condition and the same quark-gluon plasma (QGP) expansion scenario we discuss the consequences of the assumption of a local equilibrium by comparing the consequences for the nuclear modification factor $R_{AA}$ and the elliptic flows of charm quarks, scrutinizing the different components of the final distribution of charm quarks. For this purpose we employ the parton-hadron-string dynamics (PHSD) model, which is an off-shell microscopic transport approach, as well as the linearized-Boltzmann (LB) scheme obtained by coarse graining the PHSD bulk and assuming local equilibrium for the interactions of the charm quarks with the bulk. The $R_{AA}$ of charm quarks stemming from the later LB approach is also compared to a genuine fluid dynamics evolution initiated by the coarse grained PHSD, which allows us to further assess the consequences of reducing the full $\textit{n}$-body dynamics. We then proceed to a systematic comparison of PHSD (in its LB approximation) with MC@HQ, another transport model for heavy flavors which also relies on the LB approach. In particular, we investigate the consequences for the nuclear modification factor of charm quarks if we vary separately the initial heavy quark distribution function in matter, the expansion dynamics of the QGP, and the elementary interactions of heavy quarks of these models. We find that the results for both models vary significantly depending on the details of the calculation. However, both models achieve very similar predictions for key heavy quark observables for certain combinations of initial condition, bulk evolution, and interactions. In this work, we conclude that this ambiguity limits our ability to determine the different properties of the system based on the current set of observables.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Bayesian inference of in-medium baryon-baryon scattering cross sections from HADES proton flow data

Within a Bayesian statistical framework using a Gaussian Process emulator for an isospin-dependent Boltzmann-Uehling-Uhlenbeck (IBUU) transport model simulator of heavy-ion reactions at intermediate energies, we infer from the HADES proton flow data the posterior probability distribution functions of in-medium baryon-baryon scattering cross section modification factor X with respect to free-space and the corresponding incompressibility K of nuclear matter as well as their correlation function. In conclusion, the mean value of X is found to be $X$ = $1.32_ {+0.28} ^ {-0.40} $ at 68% confidence level assuming the nuclear incompressibility K will not exceed 400 MeV, providing circumstantial evidence for enhanced baryon-braryon scattering cross sections in hot and dense nuclear matter.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Moment extraction using an unfolding protocol without binning

Deconvolving (“unfolding”) detector distortions is a critical step in the comparison of cross-section measurements with theoretical predictions in particle and nuclear physics. However, most existing approaches require histogram binning while many theoretical predictions are at the level of statistical moments. We develop a new approach to directly unfold distribution moments as a function of another observable without having to first discretize the data. Our moment unfolding technique uses machine learning and is inspired by Boltzmann weight factors and generative adversarial networks (GANs). We demonstrate the performance of this approach using jet substructure measurements in collider physics. With this illustrative example, we find that our moment unfolding protocol is more precise than bin-based approaches and is as or more precise than completely unbinned methods.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS2023

Quantum Image Denoising: A Framework via Boltzmann Machines, QUBO, and Quantum Annealing

We investigate a framework for binary image denoising via restricted Boltzmann machines (RBMs) that introduces a denoising objective in quadratic unconstrained binary optimization (QUBO) form and is well-suited for quantum annealing. The denoising objective is attained by balancing the distribution learned by a trained RBM with a penalty term for derivations from the noisy image. We derive the statistically optimal choice of the penalty parameter assuming the target distribution has been well-approximated, and further suggest an empirically supported modification to make the method robust to that idealistic assumption. We also show under additional assumptions that the denoised images attained by our method are, in expectation, strictly closer to the noise-free images than the noisy images are. While we frame the model as an image denoising model, it can be applied to any binary data. As the QUBO formulation is well-suited for implementation on quantum annealers, we test the model on a D-Wave Advantage machine, and also test on data too large for current quantum annealers by approximating QUBO solutions through classical heuristics.

restricted Boltzmann machine↗
DOE OSTI2022

Training a Quantum Annealing Based Restricted Boltzmann Machine on Cybersecurity Data

A restricted Boltzmann machine (RBM) is a generative model that could be used in effectively balancing a cybersecurity dataset because the synthetic data a RBM generates follows the probability distribution of the training data. RBM training can be performed using contrastive divergence (CD) and quantum annealing (QA). QA-based RBM training is fundamentally different from CD and requires samples from a quantum computer. We present a real-world application that uses a quantum computer. Specifically, we train a RBM using QA for cybersecurity applications. The D-Wave 2000Q has been used to implement QA. RBMs are trained on the ISCX data, which is a benchmark dataset for cybersecurity. For comparison, RBMs are also trained using CD. CD is a commonly used method for RBM training. Our analysis of the ISCX data shows that the dataset is imbalanced. We present two different schemes to balance the training dataset before feeding it to a classifier. The first scheme is based on the undersampling of benign instances. The imbalanced training dataset is divided into five sub-datasets that are trained separately. A majority voting is then performed to get the result. Our results show the majority vote increases the classification accuracy up from 90.24% to 95.68%, in the case of CD. For the case of QA, the classification accuracy increases from 74.14% to 80.04%. In the second scheme, a RBM is used to generate synthetic data to balance the training dataset. We show that both QA and CD-trained RBM can be used to generate useful synthetic data. Balanced training data is used to evaluate several classifiers. Among the classifiers investigated, K-Nearest Neighbor (KNN) and Neural Network (NN) perform better than other classifiers. They both show an accuracy of 93%. Our results show a proof-of-concept that a QA-based RBM can be trained on a 64-bit binary dataset. The illustrative example suggests the possibility to migrate many practical classification problems to QA-based techniques. Further, we show that synthetic data generated from a RBM can be used to balance the original dataset.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Improved pore network models to simulate single-phase flow in porous media by coupling with lattice Boltzmann method

In this paper, different pore network models to simulate single-phase flow in porous media are built and their accuracy is evaluated. In addition to the conventional pore network model (CPNM) which consists of regular pore bodies and throat bonds, three improved pore network models (IPNMs) are developed allowing to better describing the real pore and throat geometry. The first improved pore network model (IPNM1) replaces the regular throat bond with a throat bond showing the real throat cross section. The second improvement (IPNM2) uses a series of sub-throat bonds with varying cross sections to better describe the real throat geometry, which is firstly proposed in this paper. The third model (IPNM3) extracts the real pore-throat-pore geometry without simplification. The conductance of fluid flow through these more realistic throat bonds is calculated by the lattice Boltzmann method (LBM). The accuracy and computational efficiency of the different pore network models are evaluated taking the LBM simulation over the whole porous medium as reference solution. The global permeability and detailed pressure distributions in the pores for the different pore network models are validated. The results show that the accuracy of the pore network model increases from CPNM to IPNM3, but at the expense of increasing computational cost. This study suggests that IPNM3 can replace a whole-domain LBM simulation with similar accuracy but much lower computational cost. As a first-order approximation the newly proposed IPNM2 is suggested as good compromise between accuracy and computational cost.

42 ENGINEERING↗
DOE OSTI2020

Two-Phase Fluid Flow Properties of Rough Fractures With Heterogeneous Wettability: Analysis With Lattice Boltzmann Simulations

Fractures are conduits for fluid flow in low-permeability geological formations. Multiphase flow properties of fractures are important in natural processes and in engineering applications such as the evaluation of the sealing capacity of caprocks and productivity of hydrocarbon-bearing tight rocks. Investigations of flow and transport through fractures typically focus on the effects of fracture geometric and mechanical factors such as aperture, roughness, and compressibility. The wettability of the fracture surfaces and its influence on microscale interfacial phenomena and macroscale effective transport properties are seldom studied. In this study, we investigated the effect of heterogeneous wetting properties on the displacement of water by supercritical CO2 through a series of lattice Boltzmann method simulations. The results show the evolution of the CO2 plume within a fracture is controlled by both the roughness of the aperture field and the wetting distribution. We combined these factors into a capillary pressure map that can be related to the macroscopic flow behavior of the fracture. We observed that heterogeneous wetting distributions promote the residual trapping of water where lower capillary pressures allowed for isolated water pockets in higher capillary pressure zones. Analysis of fracture unsteady relative permeability shows the effect of wetting on permeability evolution and provides support for the viscous-coupling relative permeability model. Finally, analysis of the steady-state relative permeability and saturation demonstrates a strong correlation between permeability and the standard deviation of the capillary pressure field. Thus, characterizing the distribution of wetting properties of fractures is crucial to understanding multiphase fracture flow and transport properties.

58 GEOSCIENCES↗
NASA NTRS1973

Equations for the kinetic modeling of supersonically flowing electrically excited lasers

The equations for the kinetic modeling of a supersonically flowing electrically excited laser system are presented. The work focuses on the use of diatomic gases, in particular carbon monoxide mixtures. The equations presented include the vibrational rate equation which describes the vibrational population distribution, the electron, ion and electronic level rate equations, the gasdynamic equations for an ionized gas in the presence of an applied electric field, and the free electron Boltzmann equation including flow and gradient coupling terms. The model developed accounts for vibration-vibration collisions, vibration-translation collisions, electron-molecule inelastic excitation and superelastic de-excitation collisions, charge particle collisions, ionization and three body recombination collisions, elastic collisions, and radiative decay, all of which take place in such a system. A simplified form of the free electron Boltzmann equation is developed and discussed with emphasis placed on its coupling with the supersonic flow. A brief description of a possible solution procedure for the set of coupled equations is then discussed.

Lind, R. C.↗