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 127 records · Page 7

Flowing Plasma Interaction with an Electric Sail Tether Element

Electric sails are a relatively new concept for providing high speed propellant-less propulsion. Employing multiple tethers biased to high positive voltage levels (kV), electric sails are designed to gain momentum from the solar wind by repelling solar wind protons. To maximize the area of the sail that interacts with the solar wind, electric sails rely on the formation of a large plasma sheath around each small diameter tether. Motivated by interest in advancing the development of electric sails, a set of laboratory tests has been conducted to study the interaction of a drifting plasma with a sheath formed around a small diameter tether element biased at positive voltages. The laboratory test setup was created with Debye length scaling in mind to offer a path to extrapolate (via modeling) to full scale electric sail missions. Using an instrument known as a Differential Ion Flux Probe (DIFP) the interaction between a positively biased tether element and a drifting plasma has been measured for several scenarios. Clear evidence of the tether element sheath deflecting ions has been obtained. Maps of the flow angle downstream from the tether element have been made and they show the influence of the plasma sheath. Finally, electron current collection measurements have been made for a wide range of plasma conditions and tether element bias voltages. The electron collection data will have an impact on electric sail power requirements, as high voltage power supplies and electron guns will have to be sized to accommodate the electron currents collected by each tether.

Schneider, Todd↗

Verification of an energy-conserving semi-implicit electrostatic particle-in-cell scheme for modeling high-density plasma at scale

A verification study of a semi-implicit energy-conserving electrostatic particle-in-cell algorithm is presented. The algorithm relaxes the time-step and mesh-size constraints that require resolution of the plasma period and Debye length associated with traditional explicit momentum-conserving particle-in-cell algorithms. Physical implications and applicability of using the semi-implicit scheme for modeling high-density plasmas are discussed. Where possible, numerical results are compared against analytical solutions. The simulation results indicate that the algorithm is stable at time steps larger than twice the inverse plasma frequency and cell sizes larger than the Debye length. It is found that the algorithm gives adequate results, provided that the distribution function and the spatiotemporal scales dictating the physics of the problem are resolved. As such, the algorithm may provide a robust method for kinetic modeling of high-density plasmas at scale.

Cyclotron resonance↗

Efficient simulation of optical spectra via machine learning and physical decomposition of environmental effects

Simulations of optical spectra can provide key insights to aid experimental interpretation of electronic excitation phenomena. For chromophores in the condensed phase, these spectra, which incorporate the coupling between electronic excitation and molecular and solvent nuclear motions, can be simulated using excitation energies obtained from molecular dynamics simulations of the chromophore and solvent. Here, we present a hybrid scheme that exploits machine learning and physically informed spectral densities to show that as few as 25 ground and excited state energetic gradient calculations can be used to construct models that accurately predict environment-influenced vibronic coupling in optical spectra. We demonstrate our approach for the green fluorescent protein chromophore in water and the cresyl violet chromophore in methanol. We show that our hybrid approach, employing a machine learning model for the high-frequency spectral density and an ab initio parameterized Debye spectral density for the low-frequency, results in a systematic improvement of the optical absorption lineshape, leading to a simple machine learning scheme that can be used for the simulation of spectral densities and optical spectra.

Snider, Andrew [Univ. of California, Merced, CA (U↗

Two-dimensional potential double layers and discrete auroras

This paper is concerned with the formation of the acceleration region for electrons which produce the visible auroral arc and with the formation of the inverted V precipitation region. The former is embedded in the latter, and both are associated with field-aligned current sheets carried by plasma sheet electrons. It is shown that an electron current sheet driven from the plasma sheet into the ionosphere leads to the formation of a two-dimensional potential double layer. For a current sheet of a thickness less than the proton gyrodiameter solutions are obtained in which the field-aligned potential drop is distributed over a length much greater than the Debye length. For a current sheet of a thickness much greater than the proton gyrodiameter solutions are obtained in which the potential drop is confined to a distance on the order of the Debye length. The electric field in the two-dimensional double-layer model is the zeroth-order field inherent to the current sheet configuration, in contrast to those models in which the electric field is attributed to the first-order field due to current instabilities or turbulences. The maximum potential in the two-dimensional double-layer models is on the order of the thermal energy of plasma sheet protons, which ranges from 1 to 10 keV.

Kan, J. R.↗

The electrostatics of a dusty plasma

The potential distribution in a plasma containing dust grains were derived where the Debye length can be larger or smaller than the average intergrain spacing. Three models were treated for the grain-plasma system, with the assumption that the system of dust and plasma is charge-neutral: a permeable grain model, an impermeable grain model, and a capacitor model that does not require the nearest neighbor approximation of the other two models. A gauge-invariant form of Poisson's equation was used which is linearized about the average potential in the system. The charging currents to a grain are functions of the difference between the grain potential and this average potential. Expressions were obtained for the equilibrium potential of the grain and for the gauge-invariant capacitance between the grain and the plasma. The charge on a grain is determined by the product of this capacitance and the grain-plasma potential difference.

Whipple, E. C.↗

Determination of a Two-Phase Structure of Nanocrystals: GaN and SiC

The properties of nano-crystalline materials are critically dependent on the structure of the constituent grains. Experimental conditions necessary to perform structural analysis of nanocrystalline materials as a two-phase core-surface shell system are discussed. It is shown, that a standard X-ray diffraction measurements and analysis are insufficient and may lead to incorrect conclusions as to the real structure of the materials. A new method of evaluation of powder diffraction data based on the analysis of the shift of the Bragg reflections from their perfect-lattice positions was developed. "Apparent lattice parameters" quantity, alp, was introduced and calculated from the actual positions of each individual Bragg reflection. The alp values plotted versus diffraction vector (Q) show characteristic features that are used for evaluation of the experimental results. The study was based on modeling of nano-grains and simulations of theoretical intensity profiles using the Debye functions. The method was applied to the analysis of synchrotron X-ray diffraction data of GaN and SiC nanocrystals. A presence of strained surface shell and a considerable internal pressure (GaN) in the nanoparticles was concluded.

Palosz, W.↗

Vibrational, Elastic Properties and Sound Velocities of MgTi{sub 2}O{sub 4} Spinel

Zone-centre phonon frequencies, elastic properties and sound velocities along high-symmetry directions and Debye temperature of MgTi{sub 2}O{sub 4} spinel is calculated using rigid ion model. The interatomic interactions up to third nearest neighbours are also calculated. The main results of the present paper reveals that the tetrahedral bonding is stronger than the octahedral bonding. Calculated results are compared with previously available experimental and theoretical data. It is observed that the obtained results are in good agreement with the available results in literature.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

H2OI95: A Stand-Alone Fortran Code for Evaluating the IAPWS-95 Equation-of-State Model for Water (Rev. 1)

H2OI95 is a stand-alone Fortran code for evaluating the IAPWS-95 equation-of-state model for the thermodynamic properties of water. It also evaluates the corresponding thermochemical properties of water consistent with the CODATA recommendations. The IAPWS-95 model is based on a model equation for the dimensionless Helmholtz energy for which the primary variables are the inverse reduced temperature (τ = T cr /T ) and reduced density (δ = ρ/ρ cr ). Here T is the absolute temperature (K), ρ is density (kg/m³), and the subscript “ cr ” refers to the critical point of water (647.096 K and 22.064 MPa pressure in this model, for which ρ cr is 322 kg/m³). The code solves four basic types of problems, distinguished by the specified inputs: 1. Temperature (K) and density (ρ kg/m³) or reduced density (δ) 2. Temperature (K) and pressure (MPa). 3. Temperature (K) on the saturation (liquid-vapor equilibrium) curve 4. Pressure (MPa) on the saturation curve. Each type of problem is run using a corresponding input (text) file. All but the first type of problem require iteration. For example, to solve for desired temperature and pressure, the reduced density must be adjusted to give the desired pressure. Iteration is accomplished using the Newton-Raphson method, though the secant method is also used in solving the fourth type of problem. For the last three types of problems, H2OI95 has been used to conduct numerical studies of convergence and the problem of multiple numerical solutions, some of which are not physically valid. Obtaining valid results depends mainly on appropriate choice of starting values for the density. The default values used in H2OI95 appear to consistently lead to generally desired results. With modification (not addressed here), H2OI95 can be used to support SUPCRT92 and similar codes that compute chemical thermodynamic properties of species and reactions over a wide range of temperature and pressure (273.16-1273K and 0-1000 MPa). Version 1.1 of this software calculates an extended set of water properties, including some derivable solely from IAPWS-95, such as the isobaric thermal expansivity and the isothermal compressibility, and still other properties based on supplementary models sanctioned by IAPWS, including the dynamic viscosity, thermal conductivity, surface tension, and static dielectric constant, and various other properties (notably Debye-Hückel parameters) derivable from the dielectric constant and its partial derivatives with respect to temperature and pressure. A coding error in Version 1.0 was identified while testing the implementation of the viscosity model. This error is described in Appendix C. It mainly affects calculation of “higher order” properties such as compressibility and heat capacity near the critical point.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Model for Dielectric Constant of Seawater based on L-band Measurements with Conductivity by Definition

This article reports an improvement in the model for the dielectric constant of seawater used to fit laboratory measurements of the dielectric constant at the L-band. The new model (dielectric constant as a function of salinity, temperature, and frequency) is based on the response of a polar molecule proposed by Debye and fits the same measurement as reported in earlier work but uses a functional form for conductivity, σ(S,T) , that is given by the definition of salinity. The new version of this model fits the data well and has the advantages that the relaxation time constant is allowed to be a function of temperature and salinity and is well behaved when extrapolated to high salinities.

L-band↗

Size- and Temperature-Dependent Suppression of Phonon Thermal Conductivity in Carbon Nanotube Thermoelectric Films

Heat transport in nanoscale carbon materials such as carbon nanotubes and graphene is normally dominated by phonons. Here, measurements of in-plane thermal conductivity, electrical conductivity, and thermopower are presented from 77–350 K on two films with thickness <100 nm formed from semiconducting single-walled carbon nanotubes. These measurements are made with silicon–nitride membrane thermal isolation platforms. The two films, formed from disordered networks of tubes with differing tube and bundle size, have very different thermal conductivity. One film matches a simple model of heat conduction assuming constant phonon velocity and mean free path, and 3D Debye heat capacity with a Debye temperature of 770 K. The second film shows a more complicated temperature dependence, with a dramatic drop in a relatively narrow window near 200 K where phonon contributions to thermal conductivity essentially vanish. This causes a corresponding large increase in thermoelectric figure-of-merit at the same temperature. A better understanding of this behavior can allow significant improvement in thermoelectric efficiency of these low-cost earth-abundant, organic electronic materials. Finally, heat and charge conductivity near room temperature is also presented as a function of doping, which provides further information on the interaction of dopant molecules and phonon transport in disordered nanotube films.

36 MATERIALS SCIENCE↗

Quantifying Dynamic Signal Spread in Real-Time High-Energy X-ray Diffraction

Abstract Measured intensity in high-energy monochromatic X-ray diffraction (HEXD) experiments provides information regarding the microstructure of the crystalline material under study. The location of intensity on an areal detector is determined by the lattice spacing and orientation of crystals so that changes in the heterogeneity of these quantities are reflected in the spreading of diffraction peaks over time. High temporal resolution of such dynamics can now be experimentally observed using technologies such as the mixed-mode pixel array detector (MM-PAD) which facilitates in situ dynamic HEXD experiments to study plasticity and its underlying mechanisms. In this paper, we define and demonstrate a feature computed directly from such diffraction time series data quantifying signal spread in a manner that is correlated with plastic deformation of the sample. A distinguishing characteristic of the analysis is the capability to describe the evolution from the distinct diffraction peaks of an undeformed alloy sample through to the non-uniform Debye–Scherrer rings developed upon significant plastic deformation. We build on our previous work modeling data using an overcomplete dictionary by treating temporal measurements jointly to improve signal spread recovery. We demonstrate our approach in simulations and on experimental HEXD measurements captured using the MM-PAD. Our method for characterizing the temporal evolution of signal spread is shown to provide an informative means of data analysis that adds to the capabilities of existing methods. Our work draws on ideas from convolutional sparse coding and requires solving a coupled convex optimization problem based on the alternating direction method of multipliers.

Engineering↗

Higher Ion Selectivity with Lower Energy Usage Promoted by Electro-osmotic Flow in the Transport through Conical Nanopores

The tradeoff between selectivity and throughput presents fundamental challenges to improve desalination and charge storage, salinity gradient-based energy harvesting, memory device/circuit development, and so forth. The well-known ion-current rectification and the recently resolved time-dependent transport hysteresis in conical nanopores or asymmetric nanointerfaces offer new opportunities for the selective transport of matter. This report shows that electro-osmotic flow (EOF) is an overlooked factor that increases ion selectivity while maintaining enhanced transport throughput in rectified nanoscale electrokinetic transport. The increased selectivity originates primarily from the suppression of anion flux by the fluid flow in the opposite direction under the applied electrical field. By solving the Poisson and Nernst–Planck (PNP) equations without and with coupled Navier–Stokes (PNP–NS) equations, the EOF effects on cation and anion transport are unequivocally revealed in asymmetric nanopipettes. The flux of cations and anions as well as the transference number and flow velocity are elucidated using the models and boundary conditions validated by previous experiments. A dimensionless parameter, radius over the Debye length, reveals optimal ion selectivity and energy cost at intermediate ion concentrations and nanopore sizes, up to hundreds of millimolars and tens of nanometers under time-dependent potential stimulus. Furthermore, the fundamental insights into EOF at nanointerfaces suggest new routes/strategies for better separation, analysis, and energy applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effect of double layers on magnetosphere-ionosphere coupling

The dynamic aspects of auroral current structures are reviewed with emphasis on consequences for models of microscopic turbulence (MT). A number of models of MT are introduced into a large-scale model of Alfven wave propagation to determine the effect of various models on the overall structure of auroral currents. The effect of a double layer (DL) electric field which scales with the plasma temperature and the Debye length is compared with the effect of anomalous resistivity due to electrostatic ion cyclotron turbulence in which the electric field scales with the magnetic field strength. It is shown that the DL model is less diffusive than the resistive model, indicating the possibility of narrow intense current structures.

Lysak, Robert L.↗

Fluctuating hydrodynamics of electrolytes at electroneutral scales

At mesoscopic scales electrolyte solutions are modeled by the fluctuating generalized Poisson-Nernst-Planck (PNP) equations. However, at length and time scales larger than the Debye scales, electrolytes are effectively electroneutral, and the charged-fluid PNP equations become too stiff to solve numerically. Here we formulate the isothermal incompressible equations of fluctuating hydrodynamics for reactive multispecies mixtures involving charged species in the electroneutral limit, and design a numerical algorithm to solve these equations. Our model does not assume a dilute electrolyte solution but rather treats all species on an equal footing, accounting for cross-diffusion and non-ideality using Maxwell-Stefan theory. By enforcing local electroneutrality as a constraint, we obtain an elliptic equation for the electric potential that replaces the Poisson equation in the fluctuating PNP equations. We develop a second-order midpoint predictor-corrector algorithm to solve either the charged-fluid or electroneutral equations with only a change of the elliptic solver. We use the electroneutral algorithm to study a gravitational fingering instability, triggered by thermal fluctuations, at an interface where an acid and base react to neutralize each other. Our results demonstrate that, because the four ions diffuse with very different coefficients, one must treat each ion as an individual species, and cannot treat the acid, base, and salt as neutral species. This emphasizes the differences between electrodiffusion and classical Fickian diffusion, even at electroneutral scales.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effect of double layers on magnetosphere-ionosphere coupling

The Earth's auroral zone contains dynamic processes occurring on scales from the length of an auroral zone field line which characterizes Alfven wave propagation to the scale of microscopic processes which occur over a few Debye lengths. These processes interact in a time-dependent fashion since the current carried by the Alfven waves can excite microscopic turbulence which can in turn provide dissipation of the Alfven wave energy. This review will first describe the dynamic aspects of auroral current structures with emphasis on consequences for models of microscopic turbulence. A number of models of microscopic turbulence will be introduced into a large-scale model of Alfven wave propagation to determine the effect of various models on the overall structure of auroral currents. In particular, the effects of a double layer electric field which scales with the plasma temperature and Debye length is compared with the effect of anomalous resistivity due to electrostatic ion cyclotron turbulence in which the electric field scales with the magnetic field strength. It is found that the double layer model is less diffusive than in the resistive model leading to the possibility of narrow, intense current structures.

Lysak, Robert L.↗

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.↗

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↗