Search NASA⌕ Search

SEARCH · Search NASA

Results for “Debye Model”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Quantitative Analysis of the Morphology of {101} and {001} Faceted Anatase TiO 2 Nanocrystals and Its Implication on Photocatalytic Activity

The atomistic structure and morphology (shape and size) of nanomaterials have strong influences on their physical and chemical properties. However, many characterization techniques focus exclusively on one length-scale regime or another when developing quantitative morphology/structural models. In this article, we demonstrate that powder X-ray diffraction and neutron pair distribution function (PDF) can be used to obtain accurate average morphology and atomistic structure of {001} and {101} faceted anatase TiO 2 nanocrystals based on differential evolution refinements using Debye scattering equation calculations. It is also demonstrated that the morphology polydispersity of TiO 2 nanocrystals can be effectively obtained from the diffraction data via a numerical refinement routine. The morphology refinement results are in good agreement with those from transmission electron microscopy and the modeling of small angle neutron scattering data. This method is successfully used to quantify the facet-specified photocatalytic hydrogen evolution activity of anatase TiO 2 nanocrystals with different {001} to {101} ratios. It is found that the sample with an intermediate amount of both {001} and {101} facets shows the best photocatalytic hydrogen evolution reaction (HER) activity. It is expected that the simultaneous structure and morphology refinement technique can be generally used to study the relationship between morphology and functionality of nanomaterials.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Expansion-driven Weibel instability in magnetohydrodynamics: Linear theory for static systems

The Weibel dispersion relation is obtained in the magnetohydrodynamic (MHD) approximation by including a tensor expression for temperature. MHD gives an upper cutoff wavenumber for Weibel growth identical to kinetic theory if the electron drift velocity is included in the evolution of the electron temperature. MHD overestimates growth rates compared to kinetic theory and gives maximum growth at larger wavenumbers, but it only leads to divergent results as temperature anisotropy tends to infinity. Thermal conduction in MHD lowers the growth rates and shifts maximum growth to smaller wavenumbers, but flux-limited or nonlocal thermal conduction is found to significantly limit these effects. The results are used as a first step in evaluating the potential of MHD to simulate expansion-driven Weibel instability, which occurs due to the temperature tensor only cooling in the direction of expansion. Finally, the small temperature anisotropy found in expansion-driven Weibel instability means that MHD could be an adequate model because growth rates are much less than the electron plasma frequency, and wavelengths are much greater than the Debye length.

Davies, J. R. [Univ. of Rochester, NY (United Stat↗
DOE OSTI2025

Influence of plasma screening on high-density inverse bremsstrahlung absorption

A spherical-implosion platform diagnosed with the “beamlets” scattered-light detector provides high sensitivity to the impact of plasma screening on inverse bremsstrahlung absorption. Contrary to the more restrictive screening length suggested previously [D. Turnbull et al., Phys. Rev. Lett. 130, 145103 (2023); D. Turnbull et al., Phys. Plasmas 31, 063304 (2024)], the beamlets data indicate that the electron-only Debye length is the relevant screening length for high-density inverse bremsstrahlung absorption. Furthermore, using the updated absorption model, we simulate the OMEGA direct-drive inertial confinement fusion implosion database and show that bang times are well reproduced without any ad hoc multipliers.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2014

Dependence of Lunar Surface Charging on Solar Wind Plasma Conditions and Solar Irradiation

The surface of the Moon is electrically charged by exposure to solar radiation on its dayside, as well as by the continuous flux of charged particles from the various plasma environments that surround it. An electric potential develops between the lunar surface and ambient plasma, which manifests itself in a near-surface plasma sheath with a scale height of order the Debye length. This study investigates surface charging on the lunar dayside and near-terminator regions in the solar wind, for which the dominant current sources are usually from the pohotoemission of electrons, J(sub p), and the collection of plasma electrons J(sub e) and ions J(sub i). These currents are dependent on the following six parameters: plasma concentration n(sub 0), electron temperature T(sub e), ion temperature T(sub i), bulk flow velocity V, photoemission current at normal incidence J(sub P0), and photo electron temperature T(sub p). Using a numerical model, derived from a set of eleven basic assumptions, the influence of these six parameters on surface charging - characterized by the equilibrium surface potential, Debye length, and surface electric field - is investigated as a function of solar zenith angle. Overall, T(sub e) is the most important parameter, especially near the terminator, while J(sub P0) and T(sub p) dominate over most of the dayside.

Photoemission↗
DOE OSTI2022

Modeling multiple scattering transient of an ultrashort laser pulse by spherical particles

The multiple scattering of an ultrashort laser pulse by a turbid dispersive medium (namely a cloud of bubbles in water) is investigated by means of Monte Carlo simulations. The theory of Gouesbet and Gréhan (2000) is used to derive an energetic model of the scattering transient. It is shown that the spreading and extinction of the pulse can be decoupled from the transient of scattering, which allows to describe each phenomenon individually. The transient of scattering is modeled with the Lorenz-Mie Theory and thus is also valid for a relative refractive index lower than one, contrary to the Debye series expansion which does not converge close to the critical angle. This is made possible after the introduction of a new physical object, the Scattering Impulse Response Function (SIRF) which allows to detect the different modes of scattering transient, in time and direction. Here the present approach is more generic, as it enables to simulate clouds of air bubbles in water, which was not possible previously. Two different approaches are proposed within the Monte Carlo framework. The first is a pure Monte Carlo approach where the delay due to the scattering is randomly drawn at each event, while the second is based on the transport of the whole scattering signal. They are both embedded in the Monte Carlo code Scatter3D. Both models produce equivalent trends and are validated against published numerical results. They are then applied to the multiple scattering of ultra short pulse by a cloud of bubble in water in the forward direction. The pulse spread due to the propagation in water is computed for a wide range of traveled distances and pulse durations, and the optimal pulse duration is given to minimize the pulse spread at a given distance. The main result is that the scattered photons exit the turbid medium earlier than the ballistic photons and produce a double peak related to the refraction in the bubble. This demonstrates the possibility to develop new diagnostics to characterize dynamic bubbly flows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Numerical-heating effects in atmospheric pressure streamer discharges simulated with a PIC code

Artificial heating in plasma simulations is a well-known phenomenon which occurs when, among other things, the Debye length is poorly resolved by the simulation mesh. Here, in this work, the degree to which numerical-heating occurs during a simulation of a nanosecond atmospheric pressure streamer discharge is examined. The streamer is simulated using a two-dimensional finite-element, particle-in-cell code Empire, which uses direct simulation Monte Carlo for binary particle interactions. Initially, an estimate of the numerical-heating rate applied to Empire is performed using a simple plasma model. Second, a positive atmospheric pressure streamer discharge simulation is performed to study the effects of numerical heating on plasma density, electron temperature, and streamer velocity. The nominal Debye length is approximately 1 μm and the amount of numerical heating introduced in the simulation is varied by using mesh sizes ranging from 2 μm to 20 μm. A measurable numerical heating quantity is proposed that can be used to estimate the appropriate element size and quantify the numerical-heating that can be expected over the simulation time for an atmospheric pressure streamer. In conclusion while Δx/λ D violations can be an issue it is not likely to be an issue with streamer discharges that are temporally short and occur in environments where collision frequencies are high. This result validates the rationale of grid size choices for a large amount of previously published works where Δx/λ D violation was not clearly addressed. Primary finding of this work is that numerical heating is of minor concern for plasma simulations where electron–neutral collisions are numerous such that multiple collisions can occur within a single plasma period.

Nikic, Dejan [University of New Mexico, Albuquerqu↗
NASA NTRS1977

Modeling of spacecraft charging

Three types of modeling of spacecraft charging are discussed: statistical models, parametric models, and physical models. Local time dependence of circuit upset for DoD and communication satellites, and electron current to a sphere with an assumed Debye potential distribution are presented. Four regions were involved in spacecraft charging: (1) undisturbed plasma, (2) plasma sheath region, (3) spacecraft surface, and (4) spacecraft equivalent circuit.

Whipple, E. C., Jr.↗
DOE OSTI2019

Analysis of Water Coupling in Inelastic Neutron Spectra of Uranyl Fluoride

Inelastic neutron scattering (INS) is uniquely sensitive to hydrogen due to its comparatively large thermal neutron scattering cross-section (82 b). Consequently, the inclusion of water in real samples presents significant challenges to INS data analysis due directly to the scattering strength of hydrogen. Here, we investigate uranyl fluoride (UO 2 F 2 ) with inelastic neutron scattering. UO 2 F 2 is the hydrolysis product of uranium hexafluoride (UF 6 ), and is a hygroscopic, uranyl-ion containing particulate. Raman spectral signatures are commonly used for inferential understanding of the chemical environment for the uranyl ion in UO 2 F 2 , but no direct measurement of the influence of absorbed water molecules on the overall lattice dynamics has been performed until now. To deconvolute the influence of waters on the observed INS spectra, we use density functional theory with full spectral modeling to separate lattice motion from water coupling. In particular, we present a careful and novel analysis of the Q-dependent Debye–Waller factor, allowing us to separate spectral contributions by mass, which reveals preferential water coupling to the uranyl stretching vibrations. Coupled with the detailed partial phonon densities of states calculated via DFT, we infer the probable adsorption locations of interlayer waters. We explain that a common spectral feature in Raman spectra of uranyl fluoride originates from the interaction of water molecules with the uranyl ion based on this analysis. In conclusion, the Debye–Waller analysis is applicable to all INS spectra and could be used to identify light element contributions in other systems.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
NASA NTRS1972

RF sheath and admittance characteristics of a spherical plasma probe.

Development of a radio-frequency sheath model for a spherical probe in a collisionless plasma. The method of solution is based on the quasi-static approximation and the electrostatic probe theory of Bernstein and Rabinowitz (1959). The resistive part of the admittance is ascribed to the sheath transit-time collisionless dissipation mechanism suggested by Mayer (1963) and developed by Gould (1964). Expressions are obtained for the effective sheath thickness and the equivalent resistance of the transit-time dissipation. The sheath model and, hence, the admittance are completely determined in terms of the bias potential, the probe radius, the plasma frequency, and the Debye length - i.e., there are no adjustable parameters in the proposed theory which are to be determined by experiment. The results obtained agree favorably with Cohen and Bekefi's (1971) experimental data on the conductance resonant frequency and the width of the conductance peak.

Kan, J. R.↗
DOE OSTI2026

IRMA

IRMA (In)elastic Representation of Materials As S(α,β) evaluations IRMA turns one phonon model into three outputs that usually require three separate tool chains: an evaluated nuclear-data file, predicted neutron-scattering spectra, and scattering kernels for Monte Carlo transport. The three outputs draw on a single, consistent description of the material, so the evaluation, the spectroscopy that can validate it, and the transport that uses it always agree about the physics. Nuclear data. IRMA writes ENDF-6 File 7 thermal scattering evaluations on automatically constructed (α, β) grids. This part reimplements and generalizes NJOY's LEAPR: the classic kernels reproduce freshly generated NJOY2016 tapes digit for digit and published reference tapes to about 1e-4, and the generalized paths add the exact coherent one-phonon term, anisotropic Debye-Waller tensors, coherent elastic for arbitrary crystals, and a per-species partition for polyatomic materials. The tapes feed NJOY, AMPX, FUDGE, and every transport code downstream of them. Neutron spectroscopy. The irma.spectra forward model projects the same physics onto an instrument's kinematics and resolution: INS spectra for VISION and generic indirect geometries, and 2-D S(Q,E) powder maps for direct-geometry spectrometers, from a phonopy model or straight from a phonon DOS. It can be used to predict a proposed measurement before beam time; in analysis, it supplies the calculated single-scattering counterpart of a measured spectrum, from the same material description the evaluation was built from. Monte Carlo transport. The irma.ncrystal exporter writes per-temperature scattering kernels for the companion NCrystal plugin, so McStas, OpenMC, and other NCrystal-aware codes sample the same physics. The exported kernels carry the per-site anisotropic Debye-Waller tensors, keeping directional coherent-elastic physics that NCrystal's standard scalar treatment does not represent. With the same physics inside a transport code, an entire beamline becomes a virtual experiment: IRMA's end-to-end validation ran a custom McStas implementation of the ARCS spectrometer, assembled from the existing McVine and McStas models, against measured data. From a bare crystal structure. The irma mlip front end builds the phonon model itself: a structure file and a choice of potential are enough. Nine pretrained machine-learned interatomic potentials are supported, on a laptop CPU, with no first-principles calculation; an approximate phonon model for a new material costs minutes, not a DFT campaign, and the build emits prefilled inputs for all three outputs. The result is a good starting point rather than a finished evaluation: survey-quality physics with every parameter exposed for review. A converged atomistic calculation enters the same way, as a phonopy model, when higher fidelity is needed.

Ramic, Kemal [Oak Ridge National Laboratory (ORNL)↗
NASA NTRS2013

Athermal Energy Loss from X-rays Deposited in Thin Superconducting Films on Solid Substrates

When energy is deposited in a thin-film cryogenic detector, such as from the absorption of an X-ray, an important feature that determines the energy resolution is the amount of athermal energy that can be lost to the heat bath prior to the elementary excitation systems coming into thermal equilibrium. This form of energy loss will be position-dependent and therefore can limit the detector energy resolution. An understanding of the physical processes that occur when elementary excitations are generated in metal films on dielectric substrates is important for the design and optimization of a number of different types of low temperature detector. We have measured the total energy loss in one relatively simple geometry that allows us to study these processes and compare measurements with calculation based upon a model for the various di.erent processes. We have modeled the athermal phonon energy loss in this device by finding an evolving phonon distribution function that solves the system of kinetic equations for the interacting system of electrons and phonons. Using measurements of device parameters such as the Debye energy and the thermal di.usivity we have calculated the expected energy loss from this detector geometry, and also the position-dependent variation of this loss. We have also calculated the predicted impact on measured spectral line-shapes, and shown that they agree well with measurements. In addition, we have tested this model by using it to predict the performance of a number of other types of detector with di.erent geometries, where good agreement is also found.

Kozorezov, Alexander G.↗
DOE OSTI2015

A first-principles study of cementite (Fe{sub 3}C) and its alloyed counterparts: Elastic constants, elastic anisotropies, and isotropic elastic moduli

A comprehensive computational study of elastic properties of cementite (Fe{sub 3}C) and its alloyed counterparts (M{sub 3}C (M = Al, Co, Cr, Cu, Fe, Hf, Mn, Mo, Nb, Ni, Si, Ta, Ti, V, W, Zr, Cr{sub 2}FeC and CrFe{sub 2}C) having the crystal structure of Fe{sub 3}C is carried out employing electronic density-functional theory (DFT), all-electron PAW pseudopotentials and the generalized gradient approximation for the exchange-correlation energy (GGA). Specifically, as a part of our systematic study of cohesive properties of solids and in the spirit of materials genome, following properties are calculated: (i) single-crystal elastic constants, C{sub ij}, of above M{sub 3}Cs; (ii) anisotropies of bulk, Young’s and shear moduli, and Poisson’s ratio based on calculated C{sub ij}s, demonstrating their extreme anisotropies; (iii) isotropic (polycrystalline) elastic moduli (bulk, shear, Young’s moduli and Poisson’s ratio) of M{sub 3}Cs by homogenization of calculated C{sub ij}s; and (iv) acoustic Debye temperature, θ{sub D}, of M{sub 3}Cs based on calculated C{sub ij}s. We provide a critical appraisal of available data of polycrystalline elastic properties of alloyed cementite. Calculated single crystal properties may be incorporated in anisotropic constitutive models to develop and test microstructure-processing-property-performance links in multi-phase materials where cementite is a constituent phase.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Accuracy of the explicit energy-conserving particle-in-cell method for under-resolved simulations of capacitively coupled plasma discharges

The traditional explicit electrostatic momentum-conserving particle-in-cell algorithm requires strict resolution of the electron Debye length to deliver numerical stability and accuracy. The explicit electrostatic energy-conserving particle-in-cell algorithm alleviates this constraint with minimal modification to the traditional algorithm, retaining its simplicity, ease of parallelization, and acceleration on modern supercomputing architectures. In this article, we apply the algorithm to model a one-dimensional radio frequency capacitively coupled plasma discharge relevant to industrial applications. The energy-conserving approach closely matches the results from the momentum-conserving algorithm and retains accuracy even for cell sizes up to 8 times the electron Debye length. For even larger cells, the algorithm loses accuracy due to poor resolution of steep gradients within the radio frequency sheath. Accuracy can be recovered by adopting a non-uniform grid, which resolves the sheath and allows for cell sizes up to 32 times the electron Debye length in the quasi-neutral bulk of the discharge. The effect is an up to 8 times reduction in the number of required simulation cells, an improvement that can compound in higher-dimensional simulations. We therefore consider the explicit energy-conserving algorithm as a promising approach to significantly reduce the computational cost of full-scale device simulations and a pathway to delivering kinetic simulation capabilities of use to industry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1978

Differential photoelectric charging of nonconducting surfaces in space

The photoelectric charging caused by an infinitely long strip of sunlight across a nonconducting plane is studied by use of a model which contains an electrical cutoff radius, and the results of numerical calculations are presented. The model simulates charging of a sunlit area with dimensions equal to the strip's width, exposed to a plasma with a comparatively large Debye length. Uniform potential is quickly established on a uniformly sunlit strip as a result of charge redistribution by low-energy photoelectrons. The results are in accord with a theoretical surface conductivity derived for photoelectron sheaths above highly charged sunlit areas. The surface potential, which drops sharply across the sunlight-shadow boundary, is discussed.

Pelizzari, M. A.↗
DOE OSTI2015

Proton Transport in Imidazoles: Unraveling the Role of Supramolecular Structure

The impact of supramolecular hydrogen bonded networks on dynamics and charge transport in 2-ethyl-4-methylimidazole (2E4MIm), a model proton-conducting system, is investigated by broadband dielectric spectroscopy, depolarized dynamic light scattering, viscometry, and calorimetry. It is observed that the slow, Debye-like relaxation reflecting the supramolecular structure in neat 2E4MIm is eliminated upon the addition of minute amounts of levulinic acid. This is attributed to the dissociation of imidazole molecules and the breaking down of hydrogen-bonded chains, which leads to a 10-fold enhancement of ionic conductivity.

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

A generalized Boltzmann kinetic theory for strongly magnetized plasmas with application to friction

Coulomb collisions in plasmas are typically modeled using the Boltzmann collision operator, or its variants, which apply to weakly magnetized plasmas in which the typical gyroradius of particles significantly exceeds the Debye length. Conversely, O’Neil has developed a kinetic theory to treat plasmas that are so strongly magnetized that the typical gyroradius of particles is much smaller than the distance of the closest approach in a binary collision. Here, we develop a generalized collision operator that applies across the full range of magnetization strength. Since there is no closed-form solution for the scattering cross section when plasma is strongly magnetized, the input to the collision operator is obtained by numerically calculating particle trajectories. To demonstrate novel physics associated with strong magnetization, it is used to compute the friction force on a massive test charge. In addition to the traditional stopping power component, this is found to exhibit a transverse component that is perpendicular to both the velocity and Lorentz force vectors in the strongly magnetized regime, as was predicted recently using linear response theory. As a result, good agreement is found between the collision theory and linear response theory in the regime in which both apply, but the new collision theory also applies to stronger magnetization strength regimes than the linear response theory is expected to apply in.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Intrinsic properties and strengthening mechanism of monocrystalline Ni-containing ternary concentrated solid solutions

Ternary single-phase concentrated solid solution alloys (SP-CSAs), so-called "medium entropy alloys", not only possess notable mechanical and physical properties but also form a model system linking the relatively simple binary alloys to the complex high entropy alloys. Our knowledge of their intrinsic properties is vital to understand the material behavior and to prompt future applications. To this end, three model alloys NiCoFe, NiCoCr, and NiFe-20Cr have been selected and grown as single crystals. We measured their elastic constants using an ultrasonic method, and several key materials properties, such as shear modulus, bulk modulus, elastic anisotropy, and Debye temperatures have been derived. Furthermore, nanoindentation tests have been performed on these three alloys together with Ni, NiCo and NiFe on their (100) surface, to investigate the strengthening mechanisms. NiCoCr has the highest hardness, NiFe, NiCoFe and NiFe-20Cr share a similar hardness that is apparently lower than NiCoCr; NiCo has the lowest hardness in the alloys, which is similar to elemental Ni. The Labusch-type solid solution model has been applied to interpret the nanoindentation data, with two approaches used to calculate the lattice mismatch. Finally, by adopting an interatomic spacing matrix method, the Labusch model can reasonably predict the hardening effects for the whole set of materials.

36 MATERIALS SCIENCE↗