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 199 records · Page 11

Machine learning the relationship between Debye temperature and superconducting transition temperature

Recently a relationship between the Debye temperature $Θ_D$ and the superconducting transition temperature $T_c$ of conventional superconductors has been proposed [Esterlis et al., npj Quantum Mater. 3, 59 (2018)]. The relationship indicates that $T_c$ ≤ $AΘ_D$ for phonon-mediated BCS superconductors, with $A$ being a prefactor of order ~ $0.1$. In order to verify this bound, we train machine learning (ML) models with 10 330 samples in the Materials Project database to predict $Θ_D$. Here, by applying our ML models to 9860 known superconductors in the NIMS SuperCon database, we find that the conventional superconductors in the database indeed follow the proposed bound. We also perform first-principles phonon calculations for $\mathrm{H_3S}$ and $\mathrm{LaH_{10}}$ at 200 GPa. The calculation results indicate that these high-pressure hydrides essentially saturate the bound of $T_c$ versus $Θ_D$.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Non-thermal structural transformation of diamond driven by x-rays

Intense x-ray pulses can cause the non-thermal structural transformation of diamond. At the SACLA XFEL facility, pump x-ray pulses triggered this phase transition, and probe x-ray pulses produced diffraction patterns. Time delays were observed from 0 to 250 fs, and the x-ray dose varied from 0.9 to 8.0 eV/atom. The intensity of the (111), (220), and (311) diffraction peaks decreased with time, indicating a disordering of the crystal lattice. From a Debye–Waller analysis, the rms atomic displacements perpendicular to the (111) planes were observed to be significantly larger than those perpendicular to the (220) or (311) planes. At a long time delay of 33 ms, graphite (002) diffraction indicates that graphitization did occur above a threshold dose of 1.2 eV/atom. These experimental results are in qualitative agreement with XTANT+ simulations using a hybrid model based on density-functional tight-binding molecular dynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Discrete ion stochastic continuum overdamped solvent algorithm for modeling electrolytes

In this paper we develop a methodology for the mesoscale simulation of strong electrolytes. The methodology is an extension of the fluctuating immersed-boundary approach that treats a solute as discrete Lagrangian particles that interact with Eulerian hydrodynamic and electrostatic fields. In both algorithms the immersed-boundary method of Peskin is used for particle-field coupling. Hydrodynamic interactions are taken to be overdamped, with thermal noise incorporated using the fluctuating Stokes equation, including a "dry diffusion" Brownian motion to account for scales not resolved by the coarse-grained model of the solvent. Long-range electrostatic interactions are computed by solving the Poisson equation, with short-range corrections included using an immersed-boundary variant of the classical particle-particle particle-mesh technique. Also included is a short-range repulsive force based on the Weeks-Chandler-Andersen potential. This methodology is validated by comparison to Debye-Hückel theory for ion-ion pair correlation functions, and Debye-Hückel-Onsager theory for conductivity, including the Wien effect for strong electric fields. In each case, good agreement is observed, provided that hydrodynamic interactions at the typical ion-ion separation are resolved by the fluid grid.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Thermoelastic properties of tungsten at simultaneous high pressure and temperature

In this work, the compressional ( P ) and shear wave velocities ( S ) and unit cell volumes (densities) of polycrystalline tungsten (W) have been measured simultaneously up to 10.5 GPa and 1073 K using ultrasonic interferometry in conjunction with x-ray diffraction and x-radiography techniques. Thermoelastic properties of W were derived using different methods. We obtained the isothermal bulk modulus K T 0 = 310.3(1.5) GPa, its pressure derivative K ' T 0 = 4.4(3), its temperature derivative at constant pressure ( ∂ K T / ∂ T ) P = – 0.0138 ( 1 ) GPa K – 1 and at constant volume ( ∂ K T / ∂ T ) V = – 0.0050 GPa K – 1 , the thermal expansion α ( 0, T ) = 1.02(27) × 10 –5 + 7.39(3.2) × 10 –9 T (K –1 ), as well as the pressure derivative of thermal expansion ( ∂ α / ∂ P ) T = – 1.44 ( 1 ) × 10 – 7 K – 1 GP a – 1 based on the high-temperature Birch–Murnaghan equation of state (EOS), the Vinet EOS, and thermal pressure approach. Finite strain analysis allowed us to derive the elastic properties and their pressure/temperature derivatives independent of the choice of pressure scale. A least-squares fitting yielded K S 0 = 314.5(2.5) GPa, K S 0 ' = 4.45(9), (∂ K S /∂ T ) P = – 0.0076(6) GPa K –1 , G 0 = 162.4(9) GPa, G 0 ' = 1.8(1), (∂ G /∂ T ) P = – 0.0175(9) GPa K –1 , and α 298 K = 1.23 × 10 – 5 K – 1 . Fitting current data to the Mie–Grüneisen–Debye EOS with derived θ 0 = 383.4 K yielded γ 0 = 1.81 ( 6 ) and q = 0.3 . The thermoelastic parameters obtained from various approaches are consistent with one another and comparable with previous results within uncertainties. Our current study provides a complete and self-consistent dataset for the thermoelastic properties of tungsten at high P – T conditions, which is important to improve the theoretical modeling of these materials under dynamic conditions.

36 MATERIALS SCIENCE↗

Large gyro-orbit model of ion velocity distribution in plasma near a wall in a grazing-angle magnetic field

A model is presented for the ion distribution function in a plasma at a solid target with a magnetic field $\boldsymbol {B}$ inclined at a small angle, $\alpha \ll 1$ (in radians), to the target. Adiabatic electrons are assumed, requiring $\alpha \gg \sqrt {Zm_{e}/m_{i}}$ , where $m_{e}$ and $m_{i}$ are the electron and ion mass, respectively, and $Z$ is the charge state of the ion. An electric field $\boldsymbol {E}$ is present to repel electrons, and so the characteristic size of the electrostatic potential $\phi$ is set by the electron temperature $T_{e}$ , $e\phi \sim T_{e}$ , where $e$ is the proton charge. An asymptotic scale separation between the Debye length $\lambda _{D} = \sqrt {\epsilon _0 T_{{e}} / e^{2} n_{{e}} }$ , the ion sound gyro-radius $\rho _{s} = \sqrt { m_{i} ( ZT_{e} + T_{i} ) } / (ZeB)$ and the size of the collisional region $d_{c} = \alpha \lambda _{\textrm {mfp}}$ is assumed, $\lambda _{D} \ll \rho _{s} \ll d_{c}$ . Here $\epsilon _0$ is the permittivity of free space, $n_{e}$ is the electron density, $T_{i}$ is the ion temperature, $B= |\boldsymbol {B}|$ and $\lambda _{\textrm {mfp}}$ is the collisional mean free path of an ion. The form of the ion distribution function is assumed at distances $x$ from the wall such that $\rho _{s} \ll x \ll d_{c}$ , that is, collisions are not treated. A self-consistent solution of the electrostatic potential for $x \sim \rho _{s}$ is required to solve for the quasi-periodic ion trajectories and for the ion distribution function at the target. The large gyro-orbit model presented here allows to bypass the numerical solution of $\phi (x)$ and results in an analytical expression for the ion distribution function at the target. It assumes that $\tau =T_{i}/(ZT_{e})\gg 1$ , and ignores the electric force on the quasi-periodic ion trajectory until close to the target. For $\tau \gtrsim 1$ , the model provides an extremely fast approximation to energy–angle distributions of ions at the target. These can be used to make sputtering predictions.

Physics↗

Cosmological evolution of light dark photon dark matter

Light dark photons are subject to various plasma effects, such as Debye screening and resonant oscillations, which can lead to a more complex cosmological evolution than is experienced by conventional cold dark matter candidates. Maintaining a consistent history of dark photon dark matter requires ensuring that the superthermal abundance present in the early Universe (i) does not deviate significantly after the formation of the cosmic microwave background (CMB), and (ii) does not excessively leak into the Standard Model plasma after big band nucleosynthesis (BBN). We point out that the role of nonresonant absorption, which has previously been neglected in cosmological studies of this dark matter candidate, produces strong constraints on dark photon dark matter with mass as low as 10 - 22 eV . Furthermore, we show that resonant conversion of dark photons after recombination can produce excessive heating of the intergalactic medium (IGM) which is capable of prematurely reionizing hydrogen and helium, leaving a distinct imprint on both the Ly- α forest and the integrated optical depth of the CMB. Our constraints surpass existing cosmological bounds by more than 5 orders of magnitude across a wide range of dark photon masses.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Direct observations of shock-induced melting in a porous solid using time-resolved x-ray diffraction

Here in this work we provide direct evidence of shock-induced melting and associated kinetics in a porous solid (aluminum powder) using time-resolved x-ray diffraction. Unambiguous evidence of melting in 50% porous aluminum (Al) powder samples, shocked to peak pressures between ~13 - 19 GPa, was provided by the broadening of the Debye-Scherrer ring corresponding to the (111) peak. Shocked Al powder did not melt completely in any of our experiments within the durations of measurement. Incomplete (partial) melting of the powder, even after several hundreds of nanoseconds of shock loading, provides insights into thermal transport with Al powder particles under high-pressure dynamic loading. Such insights are quite valuable for developing well-constrained melting models and thermodynamic equations of state for porous Al and other porous solids relevant to planetary and materials science.

36 MATERIALS SCIENCE↗

Conceptual Design of an Electric Sail Technology Demonstration Mission Spacecraft

There is great interest in examining the outer planets of our solar system and Heliopause region (edge of Solar System) and beyond regions of interstellar space by both the Planetary and Heliophysics communities. These needs are well docu-mented in the recent National Academy of Sciences Decadal Surveys. There is significant interest in developing revolutionary propulsion techniques that will enable such Heliopause scientific missions to be completed within 10 to15 years of the launch date. One such enabling propulsion technique commonly known as Electric Sail (E-Sail) propulsion employs positively charged bare wire tethers that extend radially outward from a rotating spacecraft spinning at a rate of one revolution per hour. Around the positively charged bare-wire tethers, a Debye Sheath is created once positive voltage is applied. This sheath stands off of the bare wire tether at a sheath diameter that is proportional to the voltage in the wire coupled with the flux density of solar wind ions within the solar system (or the location of spacecraft in the solar system. The protons that are expended from the sun (solar wind) at 400 to 800 km/sec are electrostatically repelled away from these positively charged Debye sheaths and propulsive thrust is produced via the resulting momentum transfer. The amount of thrust produced is directly proportional to the total wire length. The Marshall Space Flight Center (MSFC) Electric Sail team is currently funded via a two year Phase II NASA Innovative Advanced Concepts (NIAC) awarded in July 2015. The team's current activities are: 1) Developing a Particle in Cell (PIC) numeric engineering model from the experimental data collected at MSFC's Solar Wind Facility on the interaction between simulated solar wind interaction with a charged bare wire that can be applied to a variety of missions, 2) The development of the necessary tether deployers and tethers to enable successful de-ployment of multiple, multi km length bare tethers, 3) Controllability of the space-craft via a voltage bias to steer itself through the solar system to destinations of discovery. These activities once demonstrated analytically, will require a technology demonstration mission (TDM) around the year2020 to demonstrate that all systems work together seamlessly before a Heliophysics Electrostatic Rapid Transit System (HERTS) mission could be initiated. A notional TDM spacecraft that meets the requirements of such a mission will be showcased in this paper.

Wiegmann, Bruce M.↗

Regular sensitivity computation avoiding chaotic effects in particle-in-cell plasma methods

Particle-in-cell (PIC) simulation methods are attractive for representing species distribution functions in plasmas. However, as a model, they introduce uncertain parameters, and for quantifying their prediction uncertainty it is useful to be able to assess the sensitivity of a quantity-of-interest (QoI) to these parameters. Such sensitivity information is likewise useful for optimization. However, computing sensitivity for PIC methods is challenging due to the chaotic particle dynamics, and sensitivity techniques remain underdeveloped compared to those for Eulerian discretizations. This challenge is examined from a dual particle–continuum perspective that motivates a new sensitivity discretization. Two routes to sensitivity computation are presented and compared: a direct fully-Lagrangian particle-exact approach provides sensitivities of each particle trajectory, and a new particle-pdf discretization, which is formulated from a continuum perspective but discretized by particles to take the advantages of the same type of Lagrangian particle description leveraged by PIC methods. Since the sensitivity particles in this approach are only indirectly linked to the plasma-PIC particles, they can be positioned and weighted independently for efficiency and accuracy. The corresponding numerical algorithms are presented here in mathematical detail. The advantage of the particle-pdf approach in avoiding the spurious chaotic sensitivity of the particle-exact approach is demonstrated for Debye shielding and sheath configurations. In essence, the continuum perspective makes implicit the distinctness of the particles, which circumvents the Lyapunov instability of the N -body PIC system. The cost of the particle-pdf approach is comparable to the baseline PIC simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Landau damping and steepening of interplanetary nonlinear hydromagnetic waves

According to collisionless shock theories, the thickness of a shock front should be of the order of the characteristic lengths of the plasmas (the Debye length, the proton and Larmor radii, etc.). Chao and Lepping (1974), found, however, that 30% of the observed interplanetary shocks at 1 AU have thicknesses much larger than these characteristic lengths. It is the objective of the present paper to investigate whether the competition between nonlinear steepening and Landau damping can result in a wave of finite width that does not steepen into a shock. A heuristic model of such a wave is developed and tested by the examples of two structures that are qualitatively shocklike, but thicker than expected from theory. It is found that both events are in the process of steepening and their limiting thicknesses due to Landau damping are greater than the corresponding proton Larmor radius for both structures as observed at Mariner 5 (nearer the sun than 1 AU) but are comparable to the proton Larmor radius for Explorer (near 1 AU) observations.

Barnes, A.↗

Thermal disorder and phonon softening in the ferroelectric phase transition of lead titanate

We report a molecular dynamics study of ab initio quality of the ferroelectric phase transition in crystalline PbTi⁢O3. We model anharmonicity accurately in terms of potential energy and polarization surfaces trained on density functional theory data with modern machine learning techniques. Our simulations demonstrate that the transition has a strong order-disorder character, in agreement with diffraction experiments, and provide fresh insight into the approach to equilibrium across the phase transition. We find that the emergence and disappearance of the macroscopic polarization is driven by dipolar switching at the nanometer scale. We also computed the infrared optical absorption spectra in both the ferroelectric and the paraelectric phases, finding good agreement with the experimental Raman frequencies. Often, the almost ideal displacive character of the soft mode detected by Raman scattering in the paraelectric phase has been contrasted with the order-disorder character of the transition suggested by diffraction experiments. We settle this issue by showing that the soft mode coexists with a strong Debye relaxation associated with thermal disordering of the dipoles. The Debye relaxation feature is centered at zero frequency and appears near the transition temperature in both the ferroelectric and the paraelectric phases.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Dispersion of modified fumed silica in elastomeric nanocomposites

In polymer nanocomposites, surface modification of silica aggregates can shield Coulombic interactions that inhibit agglomeration and formation of a network of agglomerates. Surface modification is usually achieved with silane coupling agents although carbon-coating during pyrolytic silica production is also possible. Pyrogenic silica with varying surface carbon contents were dispersed in styrene-butadiene (SBR) rubber to explore the impact on hierarchical dispersion, the emergence of meso-scale structures, and the rheological response. Pristine pyrogenic silica aggregates at concentrations above a critical value (related to the Debye screening length) display correlated meso-scale structures and poor filler network formation in rubber nanocomposites due to the presence of silanol groups on the surface. Here, in the present study, flame synthesized silica with sufficient surface carbon monolayers can mitigate the charge repulsion thereby impacting network structural emergence. The impact of the surface carbon on the van der Waals enthalpic attraction, a * , is determined. The van der Waals model for polymer nanocomposites is drawn through an analogy between thermal energy, k B T, and the accumulated strain, γ. The rheological response of the emergent meso-scale structures depends on the surface density of both carbon and silanol groups.

36 MATERIALS SCIENCE↗

Elastic and thermodynamic properties of cerium-doped yttrium aluminum garnets

Cerium-doped yttrium aluminum garnets (Y 3-x Ce x Al 5 O 12 , Ce:YAGs) are promising yellow light-emitting components of solid-state white light-emitting diodes. Although there have been numerous studies examining the effects of Ce concentrations on the luminescent properties of Y 3-x Ce x Al 5 O 12 , the impacts of Ce dopant on the elastic and thermodynamic properties are not well understood. In this work, we used resonant ultrasound spectroscopy (RUS) to determine the effects of Ce doping (0.025, 0.1, 1 at. %) on the elastic and thermodynamic properties of Y 3-x Ce x Al 5 O 12 . The elastic moduli calculated via the Voigt–Reuss–Hill (VRH) method demonstrated that low Ce dopant concentrations (≤0.1 at. %) induced negligible effects on the elasticity of the YAG host matrix, while a high Ce concentration (1 at. %) yielded significant softening. RUS spectral analysis and SEM images suggested that the elastic softening originated from microstructural differences induced at higher Ce dopant concentrations. Additionally, we demonstrated an increase in elastic anisotropy at higher Ce concentrations, which further elucidated the correlations between structure and elasticity of Y 3-x Ce x Al 5 O 12 . Debye temperatures ($θ_D$), heat capacities ($C_p$), and thermal conductivities ($\mathit{κ}$) were calculated for Ce:YAGs through the relations of RUS-derived parameters (sound velocities, elastic moduli) and previously determined thermal expansion coefficients. Ce:YAG was found to have a significant reduction in $θ_D$, Cp, and $\mathit{κ}$ at Ce concentrations ≥1 at. %. Lastly, extrapolation of $C_p$ and $\mathit{κ}$ to higher temperatures allowed the modeling of thermal stress experienced by Y 3-x Ce x Al 5 O 12 disks up to 1073.15 K.

36 MATERIALS SCIENCE↗

Charge Regulation Stabilizes the Formation of Ionic Liquid‐Based Amphiphilic Oligomer Droplet Interface Bilayers

Amphiphilic charged oligomers (oligodimethylsiloxane – methylimidazolium cation, ODMS-MIM (+) ), assemble into bilayers using the droplet interface bilayer (DIB) platform, possess similar size and functionality as phospholipid bilayers, but exhibit increased stability. The oligomer ionic headgroups (MIM (+) ) are covalently bound to monodisperse, short-chain (n = 13) hydrophobic tails (ODMS). These self-assemble as monolayer brushes at the oil–aqueous interface of water droplets that are influenced by both the charged cationic headgroups, and the nature of the covalently attached tails in the organic phase. Charge regulation (CR) stabilizes the formation of ordered, molecularly close-packed brush phases, which results in highly insulating, stable DIB membranes, with contributions from specific ion-pairing effects, Debye screening, and voltage-dependent electrocompressive stresses. In the oil phase, interactions between hexadecane, a good solvent for ODMS, and the hydrophobic tails result in extended waiting times for bilayer formation compared to phospholipid DIBs, for which hexadecane is a poor solvent. Close agreement between experimental values and predictions for two key parameters, the critical membrane thickness, h c , and maximal grafted headgroup density, Γ 0 , validate an electrostatic CR model consisting of adsorption and partial neutralization of counterions at a charged interface.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ion Transport in Charged Membranes: Linking Electric-Field-Driven Mechanisms to Pore Size via Perturbation Analysis

Ion-exchange membranes are a critical component in electrochemical systems. Nevertheless, the understanding and modeling of ion transport within these porous structures have been limited by particular complexity reductions, either ignoring the dimensionality of their porous network architecture or imposing geometric assumptions (i.e., overlapping double layers). Before addressing this morphology-transport gap, a framework that relates the driving forces of transport to the geometry of a single pore is required. In this work, our modeling domain consists of a two-dimensional single pore with charged walls, connecting two identical electrolyte reservoirs. Using the Poisson-Nernst-Planck equations and regular perturbation theory, we decouple the electric fields and analyze the driving forces of ion transport, specifically electromigration and induced electroosmosis within the pore. These processes are described as analytical functions of the interaction aspect ratio, ?, defined as the ratio of the pore radius to the Debye length. Using this parameter, our study (i) describes the interplay between electromigrative and electroosmotic mechanisms that set ionic conductivity, (ii) identifies a dimensionless group of intrinsic electrolyte properties that indicates the predominant driving force, and (iii) provides a qualitative, confinement-dependent perspective on selectivity in ion-conducting membranes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

MILK : a Python scripting interface to MAUD for automation of Rietveld analysis

Modern diffraction experiments ( e.g. in situ parametric studies) present scientists with many diffraction patterns to analyze. Interactive analyses via graphical user interfaces tend to slow down obtaining quantitative results such as lattice parameters and phase fractions. Furthermore, Rietveld refinement strategies ( i.e. the parameter turn-on-off sequences) tend to be instrument specific or even specific to a given dataset, such that selection of strategies can become a bottleneck for efficient data analysis. Managing multi-histogram datasets such as from multi-bank neutron diffractometers or caked 2D synchrotron data presents additional challenges due to the large number of histogram-specific parameters. To overcome these challenges in the Rietveld software Material Analysis Using Diffraction ( MAUD ), the MAUD Interface Language Kit ( MILK ) is developed along with an updated text batch interface for MAUD . The open-source software MILK is computer-platform independent and is packaged as a Python library that interfaces with MAUD . Using MILK , model selection ( e.g. various texture or peak-broadening models), Rietveld parameter manipulation and distributed parallel batch computing can be performed through a high-level Python interface. A high-level interface enables analysis workflows to be easily programmed, shared and applied to large datasets, and external tools to be integrated with MAUD . Through modification to the MAUD batch interface, plot and data exports have been improved. The resulting hierarchical folders from Rietveld refinements with MILK are compatible with Cinema: Debye–Scherrer , a tool for visualizing and inspecting the results of multi-parameter analyses of large quantities of diffraction data. In this manuscript, the combined Python scripting and visualization capability of MILK is demonstrated with a quantitative texture and phase analysis of data collected at the HIPPO neutron diffractometer.

97 MATHEMATICS AND COMPUTING↗

Dusty plasmas in the solar system

The processes that lead to charging of dust grains in a plasma are briefly reviewed. Whereas for single grains the results have been long known, the reduction of the average charge on a grain by 'Debye screening' has only recently been discovered. This reduction can be important in the Jovian ring and in the rings of Uranus. The emerging field of gravitoelectrodynamics which deals with the motion of charged grains in a planetary magnetosphere is then reviewed. Important mechanisms for distributing grains in radial distance are due to stochastic fluctuations of the grain charge and a systematic variation due to motion through plasma gradients. The electrostatic levitation model for the formation of spokes is discussed, and it is shown that the radial transport of dust contained in the spokes may be responsible for the rich radial structure in Saturn's rings. Finally, collective effects in dusty plasmas are discussed which affect various waves, such as density waves in planetary rings and low-frequency plasma waves. The possibility of charged grains forming a Coulomb lattice is briefly described.

Goertz, C. K.↗

CheMin Instrument Performance and Calibration on Mars

The CheMin (Chemistry and Mineralogy) instrument on the Mars Science Laboratory rover Curiosity uses a CCD detector and a Co-anode X-ray tube source to acquire both mineralogy (from the pattern of Co diffraction) and chemical information (from energies of fluoresced X-rays). A key component of the CheMin instrument is the ability to move grains within sample cells during analysis, providing multiple, random grain orientations that disperse diffracted X-ray photons along Debye rings rather than producing discrete Laue spots. This movement is accomplished by piezoelectric vibration of the sample cells. A cryocooler is used to maintain the CCD at a temperature at about -50 C in order to obtain energy resolution better than 250 eV, allowing discrimination of diffracted Co K X-rays from Fe K and other fluorescent X-rays. A detailed description of CheMin is provided in [1]. The CheMin flight model (FM) is mounted within the body of Curiosity and has been operating on Mars since August 6, 2012. An essentially identical sister instrument, the CheMin demonstration model (DM), is operated in a Mars environment chamber at JPL.

Vaniman, D. T.↗