Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 OSTI2018

Measurement and modeling of the multi-wavelength optical properties ofuncoated flame-generated soot

Abstract. Optical properties of flame-generated black carbon (BC) containing soot particles were quantified at multiple wavelengths for particles produced using two different flames, a methane diffusion flame and an ethylene premixed flame. Measurements were made for: (i) nascent soot particles, (ii) thermally denuded nascent particles, and (iii) particles that were coated then thermally denuded, leading to collapse of the initially lacy, fractal-like morphology. The measured mass absorption coefficients (MAC) depended on soot maturity and generation, but were similar between flames for similar conditions. For mature soot, here corresponding to particles with volume-equivalent diameters > ~ 160 nm, the MAC and absorption Angstrom exponent values were independent of particle collapse while the single scatter albedo increased. The MAC values for these larger particles were also size-independent. Effective, theory-specific complex refractive index ( RI ) values are derived from the observations with two widely-used methods: Lorenz-Mie theory and the Rayleigh-Debye-Gans (RDG) approximation. Mie theory systematically under-predicts the observed absorption cross-sections at all wavelengths for larger particles (with x > 0.9) independent of the complex RI used, while RDG provides good agreement. Importantly, this implies that the use of Mie theory within air quality and climate models, as is common, likely leads to under-predictions in the absorption by BC, with the extent of under-prediction depending on the assumed BC size distribution and complex RI used. We suggest that it is more appropriate to assume a constant, size-independent (but wavelength-specific) MAC to represent absorption by uncoated BC particles within models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Determination of the electron-lattice coupling strength of copper with ultrafast MeV electron diffraction

Electron-lattice coupling strength governs the energy transfer between electrons and the lattice and is important for understanding the material behavior under highly non-equilibrium conditions. We report the results of employing time-resolved electron diffraction at MeV energies to directly study the electron-lattice coupling strength in 40-nm-thick polycrystalline copper excited by femtosecond optical lasers. The temporal evolution of lattice temperature at various pump fluence conditions were obtained from the measurements of the Debye-Waller decay of multiple diffraction peaks. We observed the temperature dependence of the electron-lattice relaxation time which is a result of the temperature dependence of electron heat capacity. Comparison with two-temperature model simulations reveals an electron-lattice coupling strength of (0.9 ± 0.1) × 10 17 W/m 3 /K for copper.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Capacitive response of biological membranes

We present a minimal model to analyze the capacitive response of a biological membrane subjected to a step voltage via blocking electrodes. Through a perturbative analysis of the underlying electrolyte transport equations, we show that the leading-order relaxation of the transmembrane potential is governed by a capacitive timescale, τ_{C}=λ_{D}L/D(2+Γδ^{M}/L/4+Γδ^{M}/λ_{D}), where λ_{D} is the Debye screening length, L is the electrolyte width, Γ is the ratio of the permittivity of the electrolyte to the membrane, δ^{M} is the membrane thickness, and D is the ionic diffusivity. This timescale is considerably shorter than the traditional RC timescale λ_{D}L/D for a bare electrolyte due to the membrane's low permittivity and finite thickness. Beyond the linear regime, however, salt diffusion in the bulk electrolyte drives a secondary, nonlinear relaxation process of the transmembrane potential over a longer timescale τ_{L}=L^{2}/4π^{2}D. A simple equivalent-circuit model accurately captures the linear behavior, and the perturbation expansion remains applicable across the entire range of observed physiological transmembrane potentials. Together, these findings underscore the importance of the faster capacitive timescale and nonlinear effects on the bulk diffusion timescale in determining transmembrane potential dynamics for a range of biological systems.

Farhadi, Jafar↗
NASA NTRS1994

Computing Plasma Interactions Of A Spacecraft

NASCAP/LEO (NASA Spacecraft Charging Analyzer Program for Low Earth Orbit) implements collection of mathematical models and algorithms designed to study electrostatic interaction between cold, dense plasma and spacecraft surfaces. Computes variety of electrostatic, plasma, and flow effects. Appropriate for conditions in which temperature of plasma small in comparison with spacecraft-generated potentials and Debye screening length short in comparison with dimensions of spacecraft. Related NASCAP/GEO (LEW-12973) code (NASA Charging Analyzer Program for Geosynchronous Orbit, denoted as NASCAP) appropriate for conditions which spacecraft differential potentials result from interactions with hot plasma and Debye screening length larger than dimensions of spacecraft. Object-definition portion of NASCAP/GEO code provided as part of package. NASCAP/LEO written in FORTRAN 77 and C language.

Mandell, M. J.↗
NASA NTRS1994

A two-dimensional particle-in-cell model of a dusty plasma

Dusty plasmas are present in comets, in the ring systems of the outer planets, and in the interstellar medium. A two-dimensional particle-in-cell (PIC) model of a dusty plasma is presented in this paper. The PIC code is best suited for modeling the plasma-dust interaction for large grains, with diameters of the order of a centimeter. We have modeled the charging process for an individual dust grain and the associated potential pattern in the surrounding plasma. We have also considered the case of a large number of grains in a plasma, with intergrain separations of the order of the Debye length, and have shown that the plasma becomes depleted and the charge on a dust grain is reduced, as other workers in this field have predicted (cf. C. K. Goertz, 1989). We examine the electron and ion distribution functions in the vicinity of a charged grain and demonstrate that the ions near a grain have clearly been accelerated by the electrostatic potential.

Young, B.↗
NASA NTRS2020

Plasma-environment effects on K lines of astrophysical interest III. IPs, K thresholds, radiative rates, and Auger widths in Fe ix – Fe xvi

Aims. In the context of black-hole accretion disks, we aim to compute the plasma-environment e ects on the atomic parameters used to model the decay of K-vacancy states in moderately charged iron ions, namely Fe ix – Fe xvi. Methods. We used the fully relativistic multiconfiguration Dirac–Fock method approximating the plasma electron–nucleus and electron–electron screenings with a time-averaged Debye–Hückel potential. Results. We report modified ionization potentials, K-threshold energies, wavelengths, radiative emission rates, and Auger widths for plasmas characterized by electron temperatures and densities in the ranges 10^(5)-10^(7) K and 10^(18)-10^(22)/cu. cm. Conclusions. This study confirms that the high-resolution X-ray spectrometers onboard the future XRISM and Athena space missions will be capable of detecting the lowering of the K edges of these ions due to the extreme plasma conditions occurring in accretion disks around compact objects.

J. Deprince↗
DOE OSTI2022

Modeling electrokinetic flows with the discrete ion stochastic continuum overdamped solvent algorithm

In this article we develop an algorithm for the efficient simulation of electrolytes in the presence of physical boundaries. In previous work the discrete ion stochastic continuum overdamped solvent (DISCOS) algorithm was derived for triply periodic domains, and was validated through ion-ion pair correlation functions and Debye-Hückel-Onsager theory for conductivity, including the Wien effect for strong electric fields. In extending this approach to include an accurate treatment of physical boundaries we must address several important issues. First, the modifications to the spreading and interpolation operators necessary to incorporate interactions of the ions with the boundary are described. Next we discuss the modifications to the electrostatic solver to handle the influence of charges near either a fixed potential or dielectric boundary. An additional short-ranged potential is also introduced to represent interaction of the ions with a solid wall. Finally, the dry diffusion term is modified to account for the reduced mobility of ions near a boundary, which introduces an additional stochastic drift correction. Herein, several validation tests are presented confirming the correct equilibrium distribution of ions in a channel. Additionally, the methodology is demonstrated using electro-osmosis and induced-charge electro-osmosis, with comparison made to theory and other numerical methods. Notably, the DISCOS approach achieves greater accuracy than a continuum electrostatic simulation method. We also examine the effect of under-resolving hydrodynamic effects using a “dry diffusion” approach, and find that considerable computational speedup can be achieved with a negligible impact on accuracy.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Predicting Excitation Energies of Twisted Intramolecular Charge-Transfer States with the Time-Dependent Density Functional Theory: Comparison with Experimental Measurements in the Gas Phase and Solvents Ranging from Hexanes to Acetonitrile

Electronically excited states characterized by intramolecular charge transfer play an essential role in many biological processes and optical devices. The ability to make quantitative ab initio predictions of the relative energetics involved is a challenging yet desirable goal, especially for large molecules in solution. In this work, we present a data set of 61 experimental measurements of absorption and emission processes, both in the gas phase and in solvents representing a broad range of polarities, which involve intramolecular charge transfer mediated by a nonzero, "twisted" dihedral angle between one or more donor and acceptor subunits. Among a variety of density functionals investigated within the framework of linear-response theory, the "optimally tuned" LRC-ωPBE functional, which utilizes a system-specific yet nonempirical procedure to specify the range-separation parameter, emerges as the preferred choice. For the entire set of excitation energies, involving changes in dipole moment ranging from 4 to >20 Debye, the mean signed and absolute errors are 0.02 and 0.18 eV, respectively (compared, e.g., to -0.30 and 0.30 for PBE0, 0.44 and 0.47 for LRC-ωPBEh, 0.83 and 0.83 for ωB97X-V). We analyze the performance of polarizable continuum solvation models available in Q-Chem that partition the solvent response into fast and slow time scales, and clear trends emerge when measurements corresponding to the four small 4-(dimethylamino)benzonitrile (DMABN)-like molecules and a charged species are excluded. We make the case that the large errors found only for small molecules in the gas phase and weak solvents cannot be expected to improve via the optimal tuning procedure, which enforces a condition that is exact only in the well-separated donor-acceptor limit, and present empirical evidence implicating the outsized importance for small donor-acceptor systems of relaxation effects that cannot be accounted for by the linear-response time-dependent density functional theory within the adiabatic approximation. Finally, we demonstrate the utility of the optimally tuned density functional approach by targeting the charge-transfer states of a large biomimetic model system for light-harvesting structures in Photosystem II.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2012

Plasma Wake Simulations and Object Charging in a Shadowed Lunar Crater During a Solar Storm

Within a permanently shadowed lunar crater the horizontal flow of solar wind is obstructed by upstream topography, forming a plasma wake that electrostatically diverts ions toward the crater floor and generates a surface potential that can reach kilovolts. In the present work kinetic plasma simulations are employed to investigate the morphology of a lunar crater wake during passage of a solar storm. Results are cast in terms of leading dimensionless ratios including the ion Mach number, ratio of crater depth to plasma Debye length, peak secondary electron yield, and electron temperature vs. electron impact energy at peak secondary yield. This small set of ratios allows generalization to a much wider range of scenarios. The kinetic simulation results are fed forward into an equivalent-circuit model of a roving astronaut. In very low-plasma-current environments triboelectric charging of the astronaut suit becomes effectively perpetual, representing a critical engineering concern for roving within shadowed lunar regions. Finally, simulated ion fluxes are used to explore sputtering and implantation processes within an idealized crater. It is suggested that the physics of plasma mini-wakes formed in the vicinity of permanently shadowed topography may play a critical role in modulating the enigmatic spatial distribution of volatiles at the lunar poles.

Zimmerman, Michael I.↗
NASA NTRS2020

Plasma environment effects on K lines of astrophysical interest IV. IPs, K thresholds, radiative rates, and Auger widths in Fe ii–Fe viii

Aims. Within the framework of compact-object accretion disks, we calculate plasma environment e ects on the atomic structure and decay parameters used in the modeling of K lines in lowly charged iron ions, namely Fe ii–Fe viii. Methods. For this study, we used the fully relativistic multiconfiguration Dirac–Fock method approximating the plasma electron–nucleus and electron-electron screenings with a time-averaged Debye-Hückel potential. Results. We report modified ionization potentials, K-threshold energies, wavelengths, radiative emission rates, and Auger widths for plasmas characterized by electron temperatures and densities in the ranges 10(exp 5)-10(exp 7) K and 10(exp 18)-10(exp 22) per cu.cm. In addition, we propose two universal fitting formulae to predict the IP and K-threshold lowerings in any elemental ion. Conclusions. We conclude that the high-resolution X-ray spectrometers onboard the future XRISM and ATHENA space missions will be able to detect the lowering of the K edges of these Fe ions due to the extreme plasma conditions occurring in the accretion disks around compact objects.

J. Deprince↗
DOE OSTI2018

Measurement and modeling of the multiwavelength optical properties of uncoated flame-generated soot

Optical properties of flame-generated black carbon (BC) containing soot particles were quantified at multiple wavelengths for particles produced using two different flames: a methane diffusion flame and an ethylene premixed flame. Measurements were made for (i) nascent soot particles,(ii) thermally denuded nascent particles, and (iii) particles that were coated and then thermally denuded, leading to the collapse of the initiallylacy, fractal-like morphology. The measured mass absorption coefficients(MACs) depended on soot maturity and generation but were similar between flames for similar conditions. For mature soot, here corresponding to particles with volume-equivalent diameters > ~ 160 nm, the MAC and absorption Ångström exponent (AAE) values were independent of particle collapse while the single-scatter albedo increased. The MAC values for these larger particles were also size-independent. The mean MAC value at 532 nm for larger particles was 9.1±1.1 m 2 g -1 , about 17 % higher than that recommended by Bond and Bergstrom (2006), and the AAE was close to unity. Effective, theory-specific complex refractive index(RI) values are derived from the observations with two widely used methods:Lorenz–Mie theory and the Rayleigh–Debye–Gans (RDG) approximation. Mietheory systematically underpredicts the observed absorption cross sections at all wavelengths for larger particles (with x >0.9 ) independent of the complex RI used, while RDG provides good agreement. (The dimensionless size parameter x = π d p / λ , where d p is particle diameter and λ is wavelength.) Importantly, this implies that the use of Mie theory within air quality and climate models, as is common, likely leads to under predictions in the absorption by BC, with the extent of underprediction depending on the assumed BC size distribution and complex RI used. We suggest that it is more appropriate to assume a constant,size-independent (but wavelength-specific) MAC to represent absorption by uncoated BC particles within models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Type-II superconductivity at 9K in Pb–Bi alloy

In the present work, we report the synthesis of Pb–Bi alloy with enhanced Tc of up to 9K, which is higher than that of Pb. The alloy is synthesized via a solid-state reaction route in the vacuum-encapsulated quartz tube at 700°C in an automated furnace. The synthesized sample is characterized by X-ray Diffraction(XRD) and Energy dispersive X-ray analysis(EDAX) for its phase purity and elemental composition. Rietveld refinement of XRD reveals that the end product is a majority hexagonal Pb 7 Bi 3 , with minor rhombohedral Bi. The electronic transport measurement shows metallic behavior with the Debye temperature of 108K and a superconductivity transition temperature (T c ) below 9K, which is the maximum to date for any reported Pb–Bi alloy, Pb or Bi at ambient pressure. Partial substitution of Bi at the Pb site may modify the free density of electronic states within the BCS model to attain the optimum T c , which is higher by around 2K from the reported T c of Pb. The superconductor phase diagram derived from magneto-transport measurements reveals that the synthesized alloy is a conventional superconductor with an upper critical field (H c2 ) of 3.9 T, which lies well within the Pauli paramagnetic limit. The magnetization measurements carried out following ZFC(Zero Field Cool) protocols infer that the synthesized alloy is a bulk superconductor below 9K. The isothermal M-H(Magnetization vs. Field) measurements performed below T c establish it as a type-II superconductor. Furthermore, the specific heat capacity measurements show that the Pb–Bi alloy is a strongly coupled bulk superconductor below around 9K with possibly two superconducting gaps.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Lectures on statistical mechanics

Presented here is a transcription of the lecture notes from Professor Allan N. Kaufman’s graduate statistical mechanics course Physics 212A and 212B at the University of California Berkeley from the 1972–1973 academic year. 212A addressed equilibrium statistical mechanics with topics: fundamentals (micro-canonical and sub-canonical ensembles, adiabatic law and action conservation, fluctuations, pressure, and virial theorem), classical fluids and other systems (equation of state, deviations from ideality, virial coefficients and van der Waals potential, canonical ensemble and partition function, quasistatic evolution, grand-canonical ensemble and partition function, chemical potential, simple model of a phase transition, quantum virial expansion, numerical simulation of equations of state, and phase transition), chemical equilibrium (systems with multiple species and chemical reactions, law of mass action, Saha equation, chemical equilibrium including ionization and excited states), and long-range interactions (including Coulomb, dipole, and gravitational interactions, Debye–Hückel theory, and shielding). 212B addressed nonequilibrium statistical mechanics with topics: fundamentals (definitions: realizations, moments, characteristic function, and discrete variables), Brownian motion (Langevin equation, fluctuation–dissipation theorem, spatial diffusion, Boltzmann’s H-theorem), Liouville and Klimontovich equations, Landau equation (derivation, elaboration, and H-theorem, and irreversibility), Markov processes and Fokker–Planck equation (derivations of the Fokker–Planck equation and a master equation), linear response and transport theory (linear Boltzmann equation, linear response theory of Kubo and Mori, relation of entropy production to electrical conductivity, transport relations and coefficients, normal mode solutions of the transport equations, sketch of a generalized Langevin equation method for transport theory), and an introduction to nonequilibrium quantum statistical mechanics.

plasma dynamics↗
DOE OSTI2018

Accounting for Debye sheath expansion for proud Langmuir probes in magnetic confinement fusion plasmas

A Child-Langmuir law-based method for accounting for Debye sheath expansion while fitting the current-voltage I-V characteristic of proud Langmuir probes (electrodes that extend into the volume of the plasma) is described. For Langmuir probes of a typical size used in tokamak plasmas, these new estimates of electron temperature and ion saturation current density values decreased by up to 60% compared to methods that did not account for sheath expansion. Changes to the collection area are modeled using the Child-Langmuir law and effective expansion perimeter l p , and the model is thus referred to as the “perimeter sheath expansion method.” l p is determined solely from electrode geometry, so the method may be employed without prior measurement of the magnitude of the sheath expansion effects for a given Langmuir probe and can be used for electrodes of different geometries. This method correctly predicts the non-saturating ΔI/ΔV slope for cold, low-density plasmas where sheath-expansion effects are strong, as well as for hot plasmas where ΔI/ΔV ~ 0, though it is shown that the sheath can still significantly affect the collection area in these hot conditions. The perimeter sheath expansion method has several advantages compared to methods where the non-saturating current is fitted: (1) It is more resilient to scatter in the I-V characteristics observed in turbulent plasmas. (2) It is able to separate the contributions to the ΔI/ΔV slope from sheath expansion to that of the high energy electron tail in high Te conditions. (3) It calculates the change in the collection area due to the Debye sheath for conditions where ΔI/ΔV ~ 0 and for V = V f .

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

A gyrokinetic simulation model for 2D equilibrium potential in the scrape-off layer of a field-reversed configuration

The equilibrium potential structure in the scrape-off layer (SOL) of the field-reversed configuration (FRC) can be affected by the penetration of edge biasing applied at the divertor ends. The primary focus of the paper is to establish a formulation that accurately captures both parallel and radial variations of the two-dimensional (2D) potential in SOL. The formulation mainly describes a quasi-neutral plasma with a logical sheath boundary. A full-f gyrokinetic ion model and a massless electron model are implemented in the GTC-X code to solve for the self-consistent equilibrium potential, given fixed radial potential profiles at the boundaries. The first essential point of this 2D model lies in its ability to couple radial and parallel dynamics stemming from resistive currents and drag force on ions. The model successfully recovers the fluid force balance and continuity equations. These collisional effects on 2D potential mainly appear through the density profile changes, modifying the potential through electron pressure gradient. This means an accurate prescription of electron density and temperature profiles is important in predicting the potential structure in the FRC SOL. The Debye sheath potential and the potential profiles applied at the boundaries can be additional factors contributing to the 2D variations in SOL. This comprehensive full-f scheme holds promise for future investigations into turbulent transport in the presence of the self-consistent 2D potential together with the non-Maxwellian distributions and open boundary conditions in the FRC SOL.

Physics↗
DOE OSTI2018

Anharmonic stabilization and lattice heat transport in rocksalt β -GeTe

Peierls-Boltzmann transport equation, coupled with third-order anharmonic lattice dynamics calculations, has been widely used to model lattice thermal conductivity (κ l ) in bulk crystals. However, its application to materials with structural phase transition at relatively high temperature is fundamentally challenged by the presence of lattice instabilities (imaginary phonon modes). Additionally, its accuracy suffers from the absence of higher-than-third-order phonon scattering processes, which are important near/above the Debye temperature. In this letter, we present an effective scheme that combines temperature-induced anharmonic phonon renormalization and four-phonon scattering to resolve these two theoretical challenges. Here, we apply this scheme to investigate the lattice dynamics and thermal transport properties of GeTe, which undergoes a second-order ferroelectric phase transition from rhombohedral α-GeTe to rocksalt β-GeTe at about 700 K. Our results on the high-temperature phase β-GeTe at 800 K confirm the stabilization of β-GeTe by temperature effects. We find that considering only three-phonon scattering leads to significantly overestimated κ l of 3.8 W/mK at 800 K, whereas including four-phonon scattering reduces κ l to 1.7 W/mK, a value comparable with experiments. To explore the possibility to further suppress κ l , we show that alloying β-GeTe with heavy cations such as Pb and Bi can effectively reduce κ l to about 1.0 W/mK, whereas particle size needs to be around 10 nm through nanostructuring to achieve a comparable reduction in κ l .

36 MATERIALS SCIENCE↗
NASA NTRS1981

Additional extensions to the NASCAP computer code, volume 1

Extensions and revisions to a computer code that comprehensively analyzes problems of spacecraft charging (NASCAP) are documented. Using a fully three dimensional approach, it can accurately predict spacecraft potentials under a variety of conditions. Among the extensions are a multiple electron/ion gun test tank capability, and the ability to model anisotropic and time dependent space environments. Also documented are a greatly extended MATCHG program and the preliminary version of NASCAP/LEO. The interactive MATCHG code was developed into an extremely powerful tool for the study of material-environment interactions. The NASCAP/LEO, a three dimensional code to study current collection under conditions of high voltages and short Debye lengths, was distributed for preliminary testing.

Mandell, M. J.↗
DOE OSTI2021

A Theoretical Model for Computing Freezing Point Depression of Lithium-Ion Battery Electrolytes

Reliable prediction of freezing point depression in liquid electrolytes will accelerate the development of improved Li-ion batteries which can operate in low temperature environments. In this work we establish a computational methodology to calculate activity coefficients and liquidus lines for battery-relevant liquid electrolytes. Electronic structure methods are used in conjuction with classical molecular dynamics simulations and theoretical expressions for Born solvation energy, ion-atmosphere effects from Debye-Hückel theory and solvent entropic effects. The framework uses no a priori knowledge beyond neat solvent properties and the concentration of salt. LiPF 6 in propylene carbonate (PC), LiPF 6 in dimethyl carbonate (DMC) and LiClO 4 in DMC are investigated up to 1 molal with accuracy better than 3 °C when compared to experimental freezing point measurements. Finally, we find that the difference in freezing point depression between the propylene carbonate-based electrolyte and the dimethyl carbonate electrolytes originates from the difference in the solvent dielectric constant.

25 ENERGY STORAGE↗