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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

An Improved Microwave Semiempirical Model for the Dielectric Behavior of Moist Soils

Soil semiempirical dielectric models (SEMs) are powerful, and they are generally considered a useful hybrid of both empirical and physical models. In this work, the Wang-Schmugge dielectric model is improved to more accurately estimate the relative complex dielectric constants (CDCs) of moist soils. Instead of the Debye relaxation spectrum of liquid water located outside of the soil (i.e., free out-of-soil water) adopted in the Wang-Schmugge model, the Debye relaxation formula related to the free-water component inside the soil [i.e., free soil water (FSW)], which is correlated with the soil texture, is employed in the improved SEM. In addition, the effective conductivity loss term related to both soil texture and soil moisture is introduced to explain the ionic conductivity losses of FSW. Since the soil moisture influence is reduced at high frequencies, the effective conductivity loss term related to only the soil texture is also analyzed for 14-18 GHz. As in the Wang-Schmugge model, the relative CDC of bound soil water varies with the soil volumetric moisture content when the soil moisture is lower than the maximum bound water fraction in the new model, which takes a different approach than the Mironov mineralogy-based SEM. The proposed model obtains better fitting results than the three most widely employed SEMs. Lastly, the improved model exhibits a significantly improved accuracy with a higher correlation coefficient (R²), a closer 1:1 relationship, and a lower root-mean-square error, including in the L-band, and especially in the imaginary part of the L-band.

semiempirical dielectric model (SEM)↗
DOE OSTI2023

Resonance in radio frequency sheath admittance and enhanced impurity emission near the ion cyclotron frequency

Ion cyclotron resonance heating (ICRH) is of considerable interest among all auxiliary heating techniques, because it transfers power directly to ions, targets the high-density core, and involves the cheapest radio-frequency (RF) components. During ICRH operation, RF sheaths form on the ICRH antenna itself, nearby hardware, and far-field surfaces. These sheaths are associated with large hot-spot formation and impurity emissions. This work presents high-resolution numerical modelling of RF sheaths in nuclear fusion scenarios using hPIC2, a Debye-scale particle-in-cell code. The modelling reveals a new RF sheath phenomenon which occurs when the RF is resonant with a species ion cyclotron frequency or harmonic in an oblique magnetic field. This resonance of the RF sheath causes modifications of the energy-angle distributions of the ions impacting on the walls, with consequent increase in wall impurity emission. Due to the 1/ R scaling of the magnetic field in a tokamak, such cases are possible on divertor or vacuum chamber surfaces, depending on the geometry of the tokamak and the RF heating scenario. A simple physics interpretation using a driven damped harmonic oscillator is proposed, where the ion sheath transit time plays the role of damping the RF sheath admittance. Critically, the resonance leads to increased ion heat flux at the wall as well as increased physical sputtering, despite a lack of increase in RF rectified sheath potential.

ICRH antenna↗
DOE OSTI2021

Mechanical and Thermal Properties for Uranium and U–6Nb Alloy from First-Principles Theory

Elasticity, lattice dynamics, and thermal expansion for uranium and U–6Nb alloy (elastic moduli) are calculated from density functional theory that is extended to include orbital polarization (DFT+OP). Introducing 12.5 at.% of niobium, substitutionally, in uranium softens all the cii elastic moduli, resulting in a significantly softer shear modulus (G). Combined with a nearly invariant bulk modulus (B), the quotient B/G increases dramatically for U–6Nb, suggesting a more ductile material. Lattice dynamics from a harmonic model coupled with a DFT+OP electronic structure is applied for α uranium, and the obtained phonon density of states compares well with inelastic neutron-scattering measurements. The Debye temperature associated with the lattice dynamics falls within the range of experimentally observed Debye temperatures and it also validates our quasi-harmonic (QH) phonon model. The QH Debye–Grüneisen phonon method is combined with a DFT+OP electronic structure and used to explore the anisotropic thermal expansion in α uranium. The anomalous negative thermal expansion (contraction) of the b lattice parameter of the α-phase orthorhombic cell is relatively well reproduced from a free-energy model consisting of QH-phonon and DFT+OP electronic structure contributions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Validity condition for the local sheath impedance boundary condition and a nonlocal generalization

ICRF sheaths can cause unwanted interactions of high-power RF waves with material surfaces in magnetic fusion devices. In previous work, a local RF sheath impedance boundary condition (BC) was derived for use in ICRF codes together with a microscale (i.e., Debye or sheath width scale) model for obtaining the sheath impedance used in that BC. This local RF sheath BC matches the normal component of current and electrostatic potential across the sheath-plasma interface. Collapsing the matching conditions at the sheath-plasma interface to a BC depends on the assumption of scale separation, which can be violated when conditions along the local radius of curvature of the surface vary sufficiently rapidly. The validity condition is explored in this contribution, with special attention to the case where the magnetic field approaches being tangent to the surface. When the local sheath BC no longer applies, a non-local sheath BC is developed under the assumption of a more relaxed scale separation assumption. It is shown that the non-local sheath BC reduces to the previous local sheath BC under appropriate conditions. Furthermore, a surface-integrated sheath admittance parameter describes the 2D physics in the new BC.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Mechanical and Thermal Properties of Low-Density Al 20+x Cr 20-x Mo 20-y Ti 20 V 20+y Alloys

Refractory high-entropy alloys (RHEAs) Al 20+x Cr 20-x Mo 20-y Ti 20 V 20+y ((x, y) = (0, 0), (0, 10), and (10, 15)) were computationally studied to obtain a low density and a better mechanical property. The density functional theory (DFT) method was employed to compute the structural and mechanical properties of the alloys, based on a large unit cell model of randomly distributed elements. Debye–Grüneisen theory was used to study the thermal properties of Al 20+x Cr 20-x Mo 20-y Ti 20 V 20+y . The phase diagram calculation shows that all three RHEAs have a single body-centered cubic (BCC) structure at high temperatures ranging from 1000 K to 2000 K. The RHEA Al 30 Cr 10 Mo 5 Ti 20 V 35 has shown a low density of 5.16 g/cm 3 and a hardness of 5.56 GPa. The studied RHEAs could be potential candidates for high-temperature application materials where high hardness, ductility, and low density are required.

36 MATERIALS SCIENCE↗
NASA NTRS1985

Equilibrium thermodynamics and stochastic nonlinear acoustic fields

A crystalline solid is considered to consist of a large number of incoherent nonlinear acoustic radiation sources identified with the vibrating particles of the crystalline lattice. Randomization of the field, together with the assumption of a stochastically independent, fluctuating, radiation field at the absolue zero of temperature, leads to an expression of the temperature-dependent radiation field in terms of the zero-point field. The equation is identified with the Planck distribution formula of quantum mechanics in the linear field limit. The thermodynamic state functions are also obtained in terms of the nonlinear acoustic modal energies per unit mass and reduce to the results of the Debye-Einstein stochastic quantum oscillator model in the linear field limit.

Cantrell, J. H.↗
NASA NTRS1982

Observations of structuring in the downstream region of a large spherical model in a laboratory simulated solar wind plasma

The effects of inserting a spherical conducting model, large in comparison with the Debye length, into a free streaming high-energy 1 kV) unmagnetized hydrogen plasma are investigated in order to measure energies and compositions directly relevant to solar wind and astrophysical plasma phenomena. Holding the incident plasma parameters constant, transverse profiles of the net Langmuir probe current are plotted at various locations downstream in the model wake and are divided into three regions (the shadow, transition, and boundary). Results attributable to the use of a high-energy plasma show that enhancements in the shadow exist at downstream locations where the Mach ratio is less than one, and turbulence exists in the transition region on the shadow edges and outside in the boundary region. In addition, a small current enhancement is found in the boundary and can be attributed to the plasma/model interaction. It is concluded that many similar features observed by spacecraft downstream from planetary bodies are relatively permanent and are due to the intrinsic nature of the interaction between the solar wind plasma and the obstacle.

Intriligator, D. S.↗
DOE OSTI2024

On the force exerted on a non-spherical asymmetric dust grain from homogeneous, stationary, isotropic, non-magnetized plasma

Whereas the conventional wisdom suggests that the force between non-magnetized homogeneous, stationary, isotropic plasma, and the dust grain is only possible for the case of relative plasma–grain velocity, it is shown that stationary non-spherical asymmetric dust grain immersed in stationary, non-magnetized, isotropic plasma can experience a force caused by the grain–plasma interactions. The component of the force due to scattering of plasma particles in the limit of infinite Debye length is considered analytically. Both the particle scattering and absorption force components are modeled numerically in the limits of infinite and finite Debye length using a newly developed 2D3V Aspherical Particle-in-Cell code. The code simulates interactions of dust grain of selected non-spherical asymmetric shape with plasmas using dust shape conforming coordinates. Here, the simulations confirm the existence of the force on non-spherical asymmetric grain in stationary non-magnetized plasma and show that the plasma screening effects can lead to reversal of the force direction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Fast and accurate calculation of EXAFS Debye-Waller factors in U⁢O2 using the dynamical matrix method

Theoretical modeling of bonding dynamics in metal oxides is required for predicting their thermal conductivity, catalytic activity, and mechanical properties. A primary challenge is the scarcity of experimental methods for validating theoretical predictions of these atomic-scale dynamics. This work presents a workflow that uses experimental extended x-ray absorption fine structure (EXAFS) data collected at high temperatures to validate an interatomic force field for uranium dioxide (UO2), an important model material. The validated force field is then used to drive computationally intensive molecular dynamics (MD) simulations and as input for the much faster dynamical matrix Debye-Waller (DMDW) method. The predicted values of the Debye-Waller factors from the DMDW calculations are in good agreement with those obtained from the MD simulations, with residual pair-specific differences attributable to quantum zero-point motion at low temperatures and lattice anharmonicity at high temperatures. We further show that theoretical EXAFS spectra constructed directly from DMDW-derived Debye-Waller factors reproduce the experimental data (at relatively low temperatures) with accuracy comparable to full MD-EXAFS, providing an additional validation of the choice of the potential. This study establishes a validated, rapid computational pathway for modeling bond dynamics, naturally incorporating quantum nuclear\\\\r\\\\nstatistics absent in classical simulations, which are essential for the mechanistic understanding of complex oxide materials.

58 GEOSCIENCES↗
DOE OSTI2023

Temperature Dependence of Nuclear Quadrupole Resonance and the Observation of Metal–Ligand Covalency in Actinide Complexes: 35 Cl in Cs 2 UO 2 Cl 4

We report a study of the temperature dependence of 35 Cl nuclear quadrupole resonance (NQR) transition energies and spin–lattice relaxation times (T 1 ) for 235 U-depleted dicesium uranyl tetrachloride (Cs 2 UO 2 Cl 4 ) aimed at elucidating electronic interactions between the uranium center and atoms in the equatorial plane of the UO 2 2+ ion. The transition frequency decreases slowly with temperature below 75 K and with a more rapid linear dependence above this temperature. The spin–lattice relaxation time becomes shorter with temperature, and as temperatures increase, the T 1 decrease becomes nearly quadratic. The observed trends are reproduced by a model that assumes phonon-induced fluctuations of the electric field gradient tensor and partial electron delocalization from Cl to U. The fit of the theoretical model to experimental data allows a Debye temperature of 96 K to be estimated. Finally, the generalization of this approach to investigations of covalency in actinide–ligand bonding is examined.

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

Potentiometric Studies on Ion-Transport Selectivity in Charged Gold Nanotubes

Under ideal conditions, nanotubes with a fixed negative tube-wall charge will reject anions and transport-only cations. Because many proposed nanofluidic devices are optimized in this ideally cation-permselective state, it is important to know the experimental conditions that produce ideal responses. A parameter called Ccrit, the highest salt concentration in a contacting solution that still produces ideal cation permselectivity, is of particular importance. Pioneering potentiometric studies on gold nanotubes were interpreted using an electrostatic model that states that Ccrit should occur when the Debye length in the contacting salt solution becomes equivalent to the tube radius. Since this “double-layer overlap model” (DLOM), treats all same-charge ions as identical point charges, it predicts that all same-charged cations should produce the same Ccrit. However, the effect of cation on Ccrit in gold nanotubes was never investigated. This knowledge gap has become important because recent studies with a polymeric cation-permselective nanopore membrane showed that DLOM failed for every cation studied. To resolve this issue, we conducted potentiometric studies on the effect of salt cation on Ccrit for a 10 nm diameter gold nanotube membrane. Ccrit for all cations studied were, within experimental error, the same and identical, with values predicted by DLOM. The reason DLOM prevailed for the gold nanotubes but failed for the polymeric nanopores stems from the chemical difference between the fixed negative charges of these two membranes.

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

Pore-Scale Modeling of Electrokinetics in Geomaterials

Pore-scale finite-volume continuum models of electrokinetic processes are used to predict the Debye lengths, velocity, and potential profiles for two-dimensional arrays of circles, ellipses and squares with different orientations. The pore-scale continuum model solves the coupled Navier–Stokes, Poisson, and Nernst–Planck equations to characterize the electro-osmotic pressure and streaming potentials developed on the application of an external voltage and pressure difference, respectively. Here, this model is used to predict the macroscale permeabilities of geomaterials via the widely used Carmen–Kozeny equation and through the electrokinetic coupling coefficients. The permeability results for a two-dimensional X-ray tomography-derived sand microstructure are within the same order of magnitude as the experimentally calculated values. The effect of the particle aspect ratio and orientation on the electrokinetic coupling coefficients and subsequently the electrical and hydraulic tortuosity of the porous media has been determined. These calculations suggest a highly tortuous geomaterial can be efficient for applications like decontamination and desalination.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Superconductivity out of a non-Fermi liquid: Free energy analysis

In this paper, we present an in-depth analysis of the condensation energy E c for a superconductor in a situation when superconductivity emerges out of a non-Fermi liquid due to pairing mediated by a massless boson. This is the case for electronic-mediated pairing near a quantum-critical point in metal, for pairing in SYK-type models, and for phonon-mediated pairing in the properly defined limit, when the dressed Debye frequency vanishes. We consider a subset of these quantum-critical models, in which the pairing in a channel with a proper spatial symmetry is described by an effective 0 + 1 dimensional model with the effective dynamical interaction V(Ω m ) = g¯ γ /|Ω m | γ , where γ is model-specific (the γ model). In previous papers, we argued that the pairing in the γ model is qualitatively different from that in a Fermi liquid, and the gap equation at T = 0 has an infinite number of topologically distinct solutions, Δ n (ω m ), where an integer n, running between 0 and infinity, is the number of zeros of Δ n (ω m ) on the positive Matsubara axis. This gives rise to the set of extrema of E c at E c,n , of which E c,0 is the global minimum. Here, the spectrum E c,n is discrete for a generic γ < 2 but becomes continuous at γ = 2–0. Here, we discuss in more detail the profile of the condensation energy near each E c,n and the transformation from a discrete to a continuous spectrum at γ → 2. We also discuss the free energy and the specific heat of the γ model in the normal state.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

In situ high-temperature X-ray diffraction study of Sc-doped titanium oxide nanocrystallites

Titanium dioxide is an inexpensive wide-gap highly ionic semiconductor with striking photocatalytic capabilities in several heterogeneous photoredox reactions. A small crystal size is desirable to maximize the surface area, since photocatalytic reactions occur at the surface of a photocatalyst. Presented in this work are the synthesis and microstructural characterization of 4 at.% Sc-doped TiO 2 (4SDT) prepared by water-based co-precipitation. The crystal structure of 4SDT was examined via in situ high-temperature powder X-ray diffraction experiments from 25 to 1200°C. Rietveld analysis revealed single-phase anatase up to 875°C, while at 900°C the anatase-to-rutile phase transformation occurred and at higher temperatures additional reflections of Sc-rich phases (Sc 2 TiO 5 from 975°C and Ti 3 Sc 4 O 12 or Sc 2 O 3 at 1200°C) were observed. Debye function analysis (DFA) was applied to model the total scattering pattern directly in reciprocal space, allowing the reconstruction of Ti vacancies. Both Rietveld and DFA methods were applied to estimate the nanocrystallite size and shape with consistent growth in crystallite size with temperature: an ellipsoid shape with equatorial ~4.7 nm / axial (001) ~6.9 nm at 25°C to equatorial ~27.9 nm / axial (001) ~39.6 nm at 900°C refined by Rietveld analysis, versus a cylinder shape with D a , b = 4.3 nm and size dispersion σ ab = 1.5 nm, L c = 4.9 nm and σ c = 2.3 nm at 25°C to D a , b = 21.4 nm, σ ab = 8.3 nm, L c = 23.9 and σ c = 10.9 nm at 900°C estimated by DFA. The microstructural changes obtained by Rietveld and DFA methods were supported by high-resolution transmission electron microscopy image analysis, as well as by the less direct nitrogen sorption techniques that provide information on the size of non-agglomerated and dense particles. The Ti site-occupancy factor showed a linear increase from 0.6–0.8 at 25°C to unity at 900°C for anatase, and from ~0.7 at 900°C to unity at 1200°C for rutile, via Rietveld analysis and DFA.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Inclusive jet and hadron suppression in a multistage approach

We present a new study of jet interactions in the quark-gluon plasma created in high-energy heavy-ion collisions, using a multistage event generator within the jetscape framework. We focus on medium-induced modifications in the rate of inclusive jets and high transverse momentum (high-p T ) hadrons. Scattering-induced jet energy loss is calculated in two stages: a high virtuality stage based on the matter model, in which scattering of highly virtual partons modifies the vacuum radiation pattern, and a second stage at lower jet virtuality based on the LBT model, in which leading partons gain and lose virtuality by scattering and radiation. Coherence effects that reduce the medium-induced emission rate in the matter phase are also included. The TRENTO model is used for initial conditions, and the (2 + 1)dimensional VISHNU model is used for viscous hydrodynamic evolution. Jet interactions with the medium are modeled via 2-to-2 scattering with Debye screened potentials, in which the recoiling partons are tracked, hadronized, and included in the jet clustering. Holes left in the medium are also tracked and subtracted to conserve transverse momentum. Calculations of the nuclear modification factor (R AA ) for inclusive jets and high-p T hadrons are compared to experimental measurements at the BNL Relativistic Heavy Ion Collider (RHIC) and the CERN Large Hadron Collider (LHC). Within this framework, we find that with one extra parameter which codifies the transition between stages of jet modification—along with the typical parameters such as the coupling in the medium, the start and stop criteria, etc.—we can describe these data at all energies for central and semicentral collisions without a rescaling of the jet transport coefficient q.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

Implementation of the Crank‐Nicolson Douglas‐Gunn finite‐difference time domain with complex frequency‐shifted perfectly matched layer for modeling unbounded isotropic dispersive media in two dimensions

Abstract To model open‐domain problems with Drude, Lorentz, and Debye media, the complex frequency‐shifted perfectly matched layer (CFS‐PML) is adopted to truncate the Crank‐Nicolson Douglas‐Gunn finite‐difference time‐domain (CNDG‐FDTD) region. The auxiliary differential equation (ADE) and the bilinear Z ‐transform methods are incorporated separately into the implementations of CNDG‐CFS‐PML formulations, while the ADE, piecewise linear recursive convolution (RC), and trapezoidal RC methods are utilized to model dispersive media. The proposed formulations can not only circumvent the stability condition, but also have the advantages of the CFS‐PML in attenuating the evanescent waves and reducing the late‐time reflection. Three numerical examples have been carried out to validate these formulations. The simulation results show that the proposed CNDG‐CFS‐PML algorithm is efficient in absorbing performance and saving more computational time compared with the conventional FDTD method, which leads to extensive applicability and acts as a very good prospect.

Shi, Xue‐Yang↗
DOE OSTI2022

Anisotropic temperature-dependent elastic constants and thermal conductivities of TRISO particle coatings

Tristructural isotropic (TRISO) particles show great promise as a candidate fuel for use in several next-generation high-temperature nuclear reactor designs due to their structural integrity and fuel performance at high temperatures and burnups. Computational codes exist that can simulate TRISO fuel performance characteristics and failure probabilities under extreme conditions which require knowledge of the TRISO coatings’ thermophysical properties. The thermophysical descriptions of the TRISO particle's layers (i.e., buffer, pyrolytic carbon, and silicon carbide) currently used in fuel performance codes, however, assume that many of these properties are constant with respect to temperature or texture. In this paper, we obtain the full elastic stiffness tensors of the carbon and silicon carbide layers, which have transversely isotropic symmetry. Additionally, the calculated elastic properties of each of the coating layers are in remarkable agreement with the current models at room temperature and correct orientations. Additionally, the textured 3C-SiC layer was found to exhibit novel auxetic behavior above 1500 °C. The anisotropic temperature-dependent thermal conductivities of the carbon layers were calculated using acoustical Grüneisen-Debye theory which are in excellent agreement with current models at room temperature and correct orientations. These texture- and temperature-dependent relationships can be incorporated into the thermophysical description of TRISO particles in order to more accurately model fuel performance and failure probabilities under extreme conditions in forthcoming high-fidelity computational simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Diffraction properties of cylindrically bent KAP crystals in energy range of 2.3–7.5 keV using synchrotron radiation

Verification of physics models and computer simulations are heavily reliant upon the accuracy of experimental measurements. Calibration of instrument responses becomes an important step to achieve this goal. This paper presents systematic studies of bent potassium acid phthalate (KAP) crystals using Lawrence Berkeley National Laboratories, Advanced Light Source, beamline 9.3.1 in the energy range of 2.3 to 7.5 keV. A set of KAP crystals, gradually bent from flat up to a 50.8 mm cylindrical curvature. The measured integrated reflectivity for this set of KAP crystals shows good agreement with the X-ray Oriented Program (XOP) calculations when adjusting the Debye–Waller temperature factor and using the multilamellar model in the calculations. Significant differences in rocking curve profiles were observed between experimental measurements and theory. A forward convolution model and software code were developed to include experimental parameters, allowing the investigation of the difference between measurements and calculations. After considering the experimental parameters, good agreements were obtained for the rocking curve profiles for all bending radii with a unique set of parameters. Our results show that XOP can be a useful and reliable tool to predict performance of cylindrically bent KAP crystals in this energy range.

Wu, Ming↗