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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Critical role of electron heat flux on Bohm criterion

Bohm criterion, originally derived for an isothermal-electron and cold-ion plasma, is often used as a rule of thumb for more general plasmas. Here, we establish a more precise determination of the Bohm criterion that are quantitatively useful for understanding and modeling collisional plasmas that still have collisional mean-free-path much greater than plasma Debye length. Specifically, it is shown that electron heat flux, rather than the isothermal electron assumption, is what sets the Bohm speed to be $\sqrt{k_B(T_e||+3T_i||)/m_i}$ with T e,i∥ the electron and ion parallel temperature at the sheath entrance and m i the ion mass.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2004

Diffraction Studies of the Atomic Vibrations of Bulk and Surface Atoms in the Reciprocal and Real Spaces of Nanocrystalline SiC

To describe and evaluate the vibrational properties of nanoparticles it is necessary to distinguish between the surface and the core of the particles. Theoretical calculations show that vibrational density of states of the inner atoms of nanograins is similar to bulk material but shifted to higher energies which can be explained by the fact that the gain core is stressed (hardened) due to the presence of internal pressure. Theoretical calculations also show that there is a difference between vibrational properties of a crystal lattice of the grain interior in isolated particles and in a dense (sintered) nanocrystalline material. This is probably due to a coupling of the modes inside the grains via the grain boundaries in dense nanocrystalline bodies. We examined strains present in the surface shell based on examination of diamond and Sic nanocrystals in reciprocal (Bragg-type scattering) and real (PDF analysis) space analysis of neutron diffraction data. Recently we examined the atomic thermal motions in nanocrystalline Sic based on the assumption of a simple Einstein model for uncorrelated atomic notions. According to this model, the Bragg intensity is attenuated as a function of scattering angle by the Debye-Waller factor. Based on this assumption overall temperature factors were determined from the Wilson plots.

Stelmakh, S.↗
DOE OSTI2022

Accurate quantification of lattice temperature dynamics from ultrafast electron diffraction of single-crystal films using dynamical scattering simulations

In ultrafast electron diffraction (UED) experiments, accurate retrieval of time-resolved structural parameters, such as atomic coordinates and thermal displacement parameters, requires an accurate scattering model. Unfortunately, kinematical models are often inaccurate even for relativistic electron probes, especially for dense, oriented single crystals where strong channeling and multiple scattering effects are present. This article introduces and demonstrates dynamical scattering models tailored for quantitative analysis of UED experiments performed on single-crystal films. As a case study, we examine ultrafast laser heating of single-crystal gold films. Comparison of kinematical and dynamical models reveals the strong effects of dynamical scattering within nm-scale films and their dependence on sample topography and probe kinetic energy. Applying to UED experiments on an 11 nm thick film using 750 keV electron probe pulses, the dynamical models provide a tenfold improvement over a comparable kinematical model in matching the measured UED patterns. Also, the retrieved lattice temperature rise is in very good agreement with predictions based on previously measured optical constants of gold, whereas fitting the Debye–Waller factor retrieves values that are more than three times lower. Altogether, these results show the importance of a dynamical scattering theory for quantitative analysis of UED and demonstrate models that can be practically applied to single-crystal materials and heterostructures.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Molecular dynamics and continuum analyses of the electrokinetic zeta potential in nanostructured slit channels

This work presents a theoretical and numerical study of electrokinetic flow and the zeta potential for the case of slit channels with nanoscale surface features of dimensions comparable to the Debye length, by employing molecular dynamics simulations and continuum-level analyses. Additionally, a simple analytical model for considering the average effect of such surface features producing nanoscale roughness is proposed by employing matched asymptotic solutions for the charge density and fluid flow field, and matching conditions that satisfy electroneutrality and the Onsager reciprocal relation between the electroosmotic flow rate and streaming current. The proposed analytical model quantitatively accounts for results from molecular dynamics simulations that consider the presence of ion solvation shells and surface hydration layers. Our analysis indicates that a simultaneous knowledge of the electroosmotic and pressure-driven flow rate or streaming current can be instrumental to unambiguously determine the zeta potential in the presence of nanoscale surface features.

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

Computational exploration of biomedical HfNbTaTiZr and Hf 0.5 Nb 0.5 Ta 0.5 Ti 1.5 Zr refractory high-entropy alloys

Refractory high entropy alloys (RHEAs) have been proven to be a potential candidate in the biomedical field due to their balanced mechanical properties and biocompatible composition. Recent experimental findings show that RHEAs like HfNbTaTiZr and Hf 0.5 Nb 0.5 Ta 0.5 Ti 1.5 Zr have good mechanical properties such as high polarization and wear resistance than others which establish them as potential materials for biomedical application. In this work, we performed first-principles density functional theory calculations on the mechanical and thermal properties of HfNbTaTiZr and Hf 0.5 Nb 0.5 Ta 0.5 Ti 1.5 Zr. The predicted lattice constant, density, Young's modulus, and Vickers hardness are consistent with the available experimental report, which verifies the accuracy of the applied model. The thermal coefficient of linear expansion of both RHEAs has been investigated by utilizing the Debye theory. The present methods could be applied to study other future RHEAs on exploration of their physical properties.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Direct implicit and explicit energy-conserving particle-in-cell methods for modeling of capacitively coupled plasma devices

Achieving large-scale kinetic modeling is a crucial task for the development and optimization of modern plasma devices. With the trend of decreasing pressure in applications, such as plasma etching, kinetic simulations are necessary to self-consistently capture the particle dynamics. The standard, explicit, electrostatic, momentum-conserving particle-in-cell method suffers from restrictive stability constraints on spatial cell size and temporal time step, requiring resolution of the electron Debye length and electron plasma period, respectively. This results in a very high computational cost, making the technique prohibitive for large volume device modeling. We investigate the direct implicit algorithm and the explicit energy conserving algorithm as alternatives to the standard approach, both of which can reduce computational cost with a minimal (or controllable) impact on results. These algorithms are implemented into the well-tested EDIPIC-2D and LTP-PIC codes, and their performance is evaluated via 2D capacitively coupled plasma discharge simulations. The investigation reveals that both approaches enable the utilization of cell sizes larger than the Debye length, resulting in a reduced runtime, while incurring only minor inaccuracies in plasma parameters. The direct implicit method also allows for time steps larger than the electron plasma period; however, care must be taken to avoid numerical heating or cooling. It is demonstrated that by appropriately adjusting the ratio of cell size to time step, it is possible to mitigate this effect to an acceptable level.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Low Mach number fluctuating hydrodynamics for electrolytes

Here, we formulate and study computationally the low Mach number fluctuating hydrodynamic equations for electrolyte solutions. We are also interested in studying transport in mixtures of charged species at the mesoscale, down to scales below the Debye length, where thermal fluctuations have a significant impact on the dynamics. Continuing our previous work on fluctuating hydrodynamics of multicomponent mixtures of incompressible isothermal miscible liquids (A. Donev, et al., Physics of Fluids, 27, 3, 2015), we now include the effect of charged species using a quasielectrostatic approximation. Localized charges create an electric field, which in turn provides additional forcing in the mass and momentum equations. Our low Mach number formulation eliminates sound waves from the fully compressible formulation and leads to a more computationally efficient quasi-incompressible formulation. Furthermore, we demonstrate our ability to model saltwater (NaCl) solutions in both equilibrium and nonequilibrium settings. We show that our algorithm is second-order in the deterministic setting, and for length scales much greater than the Debye length gives results consistent with an electroneutral/ambipolar approximation. In the stochastic setting, our model captures the predicted dynamics of equilibrium and nonequilibrium fluctuations. We also identify and model an instability that appears when diffusive mixing occurs in the presence of an applied electric field.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1985

Wakes and differential charging of large bodies in low Earth orbit

Highlights of earlier results using the Inside-Out WAKE code on wake structures of LEO spacecraft are reviewed. For conducting bodies of radius large compared with the Debye length, a high Mach number wake develops a negative potential well. Quasineutrality is violated in the very near wake region, and the wake is relatively empty for a distance downstream of about one half of a Mach number of radii. There is also a suggestion of a core of high density along the axis. A comparison of rigorous numerical solutions with in situ wake data from the AE-C satellite suggests that the so called neutral approximation for ions (straight line trajectories, independent of fields) may be a reasonable approximation except near the center of the near wake. This approximation is adopted for very large bodies. Work concerned with the wake point potential of very large nonconducting bodies such as the shuttle orbiter is described. Using a cylindrical model for bodies of this size or larger in LEO (body radius up to 10 to the 5th power Debye lengths), approximate solutions are presented based on the neutral approximation (but with rigorous trajectory calculations for surface current balance). There is a negative potential well if the body is conducting, and no well if the body is nonconducting. In the latter case the wake surface itself becomes highly negative. The wake point potential is governed by the ion drift energy.

Parker, L. W.↗
DOE OSTI2016

Monte Carlo simulations of ionization potential depression in dense plasmas

A particle-particle grand canonical Monte Carlo model with Coulomb pair potential interaction was used to simulate modification of ionization potentials by electrostatic microfields. The Barnes-Hut tree algorithm [J. Barnes and P. Hut, Nature 324, 446 (1986)] was used to speed up calculations of electric potential. Atomic levels were approximated to be independent of the microfields as was assumed in the original paper by Ecker and Kröll [Phys. Fluids 6, 62 (1963)]; however, the available levels were limited by the corresponding mean inter-particle distance. The code was tested on hydrogen and dense aluminum plasmas. The amount of depression was up to 50% higher in the Debye-Hückel regime for hydrogen plasmas, in the high density limit, reasonable agreement was found with the Ecker-Kröll model for hydrogen plasmas and with the Stewart-Pyatt model [J. Stewart and K. Pyatt, Jr., Astrophys. J. 144, 1203 (1966)] for aluminum plasmas. Our 3D code is an improvement over the spherically symmetric simplifications of the Ecker-Kröll and Stewart-Pyatt models and is also not limited to high atomic numbers as is the underlying Thomas-Fermi model used in the Stewart-Pyatt model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Multiple-Scale Physics During Magnetic Reconnection

Magnetic reconnection is a key fundamental process in magnetized plasmas wherein the global magnetic topology is modified and stored energy is transferred from fields to particles. Reconnection is an inherently local process, and mechanisms to couple global-scale dynamics are not well understood. This dissertation explores two different mechanisms for cross-scale coupling during magnetic reconnection. As one example, we theoretically examine reconnection in a collisionless plasma using particle-in-cell simulations and demonstrate that large scale reconnection physics can couple to and drive microscopic instabilities, even in two-dimensional systems if significant scale separation exists between the Debye length and the electron skin depth. The physics underlying these instabilities is explained using simple theoretical models, and their potential connection to existing discrepancies between laboratory experiments and numerical simulations is explored. In three-dimensional systems, these instabilities are shown to generate anomalous resistivity that balances a substantial fraction of the electric field. In contrast, we also use experiments to investigate cross-scale couplings during reconnection in a collisional plasma. A leading candidate for coupling global and local scales is the hierarchical breakdown of elongated, reconnecting current sheets into numerous smaller current sheets -– the plasmoid instability. In the Magnetic Reconnection Experiment (MRX), recent hardware improvements have extended the accessible parameter space allowing for the study of long-lived, elongated current sheets. Moreover, by using Argon, reproducible and collisional plasmas are produced, which allow for a detailed statistical study of collisional reconnection. As a result, we have conclusively measured the onset of sub-ion-scale plasmoids during resistive, anti-parallel reconnection for the first time. The current sheet thickness is intermediate between ion and electron kinetic scales such that the plasma is in the Hall-MHD regime. Surprisingly, plasmoids are observed at Lundquist numbers < 100 well below theoretical predictions (> 10,000). The number of plasmoids scales with both Lundquist number and current sheet aspect ratio. The Hall quadrupolar fields are shown to suppress plasmoids. Finally, plasmoids are shown to couple local and global physics by enhancing the reconnection rate. These results are compared with prior studies of tearing and plasmoid instability, and implications for astrophysical plasmas, laboratory experiments, and theoretical studies of reconnection are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Nanocrystals in Molten Salts and Ionic Liquids: Experimental Observation of Ionic Correlations Extending beyond the Debye Length

The nature of the interface between the solute and the solvent in a colloidal solution has attracted attention for a long time. For example, the surface of colloidal nanocrystals (NCs) is specially created to impart colloidal stability in a variety of polar and nonpolar solvents. This work focuses on a special type of colloids where the solvent is a molten inorganic salt or organic ionic liquid. The stability of such colloids is difficult to rationalize because solvents with high density of mobile charges efficiently screen the electrostatic double-layer repulsion, and purely ionic molten salts represent an extreme case where the Debye length is only similar to 1 angstrom. We introduce a detailed investigation of NC dispersions in molten salts and ionic liquids using small-angle X-ray scattering (SAXS), atomic pair distribution function (PDF) analysis and molecular dynamics (MD) simulations. Our SAXS analysis confirms that a wide variety of NCs (Pt, CdSe/CdS, InP, InAs, ZrO 2 ) can be uniformly dispersed in molten salts like AlCl 3 /NaCl/KCl (AlCl 3 /AlCl 4 - ) and NaSCN/KSCN and in ionic liquids like 1-butyl-3-methylimidazolium halides (BMIM + X - , where X = Cl, Br, I). By using a combination of PDF analysis and molecular modeling, we demonstrate that the NC surface induces a solvent restructuring with electrostatic correlations extending an order of magnitude beyond the Debye screening length. These strong oscillatory ion-ion correlations, which are not accounted by the traditional mechanisms of steric and electrostatic stabilization of colloids, offer additional insight into solvent-solute interactions and enable apparently "impossible" colloidal stabilization in highly ionized media.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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.↗