Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “anisotropic diffusion”

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 415 records · Page 23

Monte-Carlo modeling of phonon thermal transport using DFT-based anisotropic dispersion relations over the full Brillouin zone

In this work, we present a Monte Carlo (MC) approach to solve the phonon Boltzmann transport equation (BTE) in which the anisotropic phonon dispersion relations over the full Brillouin zone (BZ) are used. In this approach, the discretization of the BZ used to compute the phonon relaxation time places constraints on the direction of scattered phonons in the real-space simulation domain. The phonon dispersion and phonon relaxation times are calculated using the density functional theory (DFT) approach. The modified MC approach is validated by a close examination of its ability to simulate phonon transport in both the ballistic and diffusive regimes for multiple materials including GaAs, InAs, ThO 2 , and α-U. In doing so, the phonon thermal conductivities from 100 K to 1000 K are calculated and compared with traditional non-transport solution of the phonon BTE. It is found that the phonon thermal conductivities of α-U and ThO 2 obtained from MC simulations using isotropic dispersion are larger than the values obtained using anisotropic phonon dispersion relations over the full BZ. The effect of phonon-defect scattering on the thermal conductivity of ThO 2 is also studied as an application of the current MC approach and found to agree with previously computed values in the literature. The MC solver developed here has been parallelized as a step to demonstrate its potential to solving computationally intensive phonon thermal transport problems at the mesoscale.

36 MATERIALS SCIENCE↗

Constructing stable interface layer for boosting high-voltage cycling performance of single-crystal Ni-rich cathodes

Increasing demand for electric vehicles (EVs) worldwide and the requirements for environmental protection have greatly promoted the development of low-Co, Ni-rich layered cathodes due to their high energy density, reasonable cost, and less environmental pollution. Nevertheless, serious performance degradation and safety concerns resulting from structural/interfacial instability under high operating voltages (≥4.3 V) have greatly hindered its commercialization. Herein, we propose a feasible surface modification strategy by introducing a multifunctional Li 2 SiO 3 (LSO) coating layer onto the surface of Co less single-crystalline LiNi 0.63 Co 0.07 Mn 0.3 O 2 (NCM) to solve the above challenges. Additionally, the uniform Li 2 SiO 3 coating layer has a unique three-dimensional (3D) ion diffusion channel, which greatly promotes the transmission of Li + and alleviates interfacial stress accumulation. Equally important, the electrochemically inert Li 2 SiO 3 coating layer has a stable structure framework with strong Si-O bonds, which inhibit anisotropic particle volume expansion and the occurrence of side reactions during long-term cycling. As a consequence, LSO modified NCM exhibits greatly improved electrochemical properties even at harsh testing conditions and achieves capacity retention of 79.1% after 200 cycles at 25°C within 2.95-4.5 V in half cells. Furthermore, it provides an outstanding capacity retention of 94.7% after 300 cycles at 25°C within 2.95-4.4 V when tested in pouch full cells.

25 ENERGY STORAGE↗

Simultaneous in Situ X-ray Scattering and Infrared Imaging of Polymer Extrusion in Additive Manufacturing

In situ wide-angle X-ray scattering together with infrared imaging was performed during three-dimensional material extrusion printing and correlated with the development of the crystalline structure and subsequent thermomechanical properties. Identical samples were printed with nozzle motion either along the short axis or the long axis. The short axis mode had higher thermal retention, which resulted in later onset of crystal structure. The longer time spent at temperatures between the glass transition and the melting point produced samples with higher degree of crystallinity but also significantly increased brittleness. The tracer diffusion coefficient D ( T ), together with its temperature dependence, was measured using neutron reflectivity, and the total interdiffusion length between filaments was then calculated using D ( T ) for each temperature point, as determined by the measured thermal profiles. This allowed us to define the time/temperature plane that yielded the minimum diffusion length ΔL that provides mechanical integrity of the printed features ( ΔL less than the radius of gyration of the poly(l-lactide)). The model was probed by printing structures at four nozzle temperatures and measuring the time dependence of the thermal profiles at filaments in the horizontal and vertical positions. In conclusion, the data indicated that the thermal retention was anisotropic, where higher values were obtained in the horizontal plane. Mechanical measurements indicated large differential increases in the torsional strength, corresponding to the direction with increased thermal retention.

36 MATERIALS SCIENCE↗

High-Order Hyperbolic Residual-Distribution Schemes on Arbitrary Triangular Grids

In this paper, we construct high-order hyperbolic residual-distribution schemes for general advection-diffusion problems on arbitrary triangular grids. We demonstrate that the second-order accuracy of the hyperbolic schemes can be greatly improved by requiring the scheme to preserve exact quadratic solutions. We also show that the improved second-order scheme can be easily extended to third-order by further requiring the exactness for cubic solutions. We construct these schemes based on the LDA and the SUPG methodology formulated in the framework of the residual-distribution method. For both second- and third-order-schemes, we construct a fully implicit solver by the exact residual Jacobian of the second-order scheme, and demonstrate rapid convergence of 10-15 iterations to reduce the residuals by 10 orders of magnitude. We demonstrate also that these schemes can be constructed based on a separate treatment of the advective and diffusive terms, which paves the way for the construction of hyperbolic residual-distribution schemes for the compressible Navier-Stokes equations. Numerical results show that these schemes produce exceptionally accurate and smooth solution gradients on highly skewed and anisotropic triangular grids, including curved boundary problems, using linear elements. We also present Fourier analysis performed on the constructed linear system and show that an under-relaxation parameter is needed for stabilization of Gauss-Seidel relaxation.

Mazaheri, Alireza↗

Nonlinear diffusion acceleration for the multigroup transport equation discretized with S{sub N} and continuous FEM with rattlesnake

Nonlinear diffusion acceleration (NDA) can improve the performance of a neutron transport solver significantly especially for the multigroup eigenvalue problems. The high-order transport equation and the transport-corrected low-order diffusion equation form a nonlinear system in NDA, which can be solved via a Picard iteration. The consistency of the correction of the low-order equation is important to ensure the stabilization and effectiveness of the iteration. It also makes the low-order equation preserve the scalar flux of the high-order equation. In this paper, the consistent correction for a particular discretization scheme, self-adjoint angular flux (SAAF) formulation with discrete ordinates method (S{sub N}) and continuous finite element method (CFEM) is proposed for the multigroup neutron transport equation. Equations with the anisotropic scatterings and a void treatment are included. The Picard iteration with this scheme has been implemented and tested with RattleS{sub N}ake, a MOOSE-based application at INL. Convergence results are presented. (authors)

97 MATHEMATICS AND COMPUTING↗

Probing the Dipole-Bound State in the 9-Phenanthrolate Anion by Photodetachment Spectroscopy, Resonant Two-Photon Photoelectron Imaging, and Resonant Photoelectron Spectroscopy

Valence-bound anions with a dipolar core can support dipole-bound states (DBSs) below the electron detachment threshold. The highly diffuse DBS observed is usually of σ symmetry with an s-like orbital. Recently, a π-type DBS was observed experimentally in the 9-anthrolate anion (9AT – ) and it was shown to be stabilized due to the large anisotropic polarizability of the 9AT core. To confirm the general existence of π-DBS and its structural dependence, here we report an investigation of the 9-phenanthrolate anion (9PT – ), which has a different structure and lower symmetry than 9AT – . Photodetachment spectroscopy revealed a DBS 257 cm -1 below the detachment threshold of 9PT – at 19,627 cm -1 (2.4334 eV). Resonant two-photon photoelectron imaging indeed showed a π symmetry for the DBS. Similar to that observed in 9AT – , the π-DBS in 9PT – is also stabilized by the anisotropic polarizability of the 9PT core and accessed via non-adiabatic population transfer from the initially populated σ-DBS. Photodetachment spectroscopy unveiled nine above-threshold vibrational resonances of the DBS, resulting in nine highly non-Franck-Condon resonant photoelectron spectra by tuning the detachment laser to the vibrational resonances. In conclusion, the combination of photodetachment spectroscopy and resonant photoelectron spectroscopy allowed frequencies for nine vibrational modes of the 9-phenathroxy radical to be measured, including the six lowest frequency bending modes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Synthesis of Armchair Graphene Nanoribbons on Germanium-on-Silicon

The synthesis of graphene nanoribbons on complementary metal– oxide–semiconductor-compatible substrates is a significant challenge hindering their integration into commercial semiconductor electronics. Here, the bottom-up synthesis of armchair graphene nanoribbons on epilayers of Ge on Si(001) via chemical vapor deposition is demonstrated. The synthesis leverages the previous discovery that graphene crystal growth can be driven with an extreme shape anisotropy on Ge(001) surfaces to directly form nanoribbons. However, compared to nanoribbon synthesis on Ge(001), synthesis on Ge/Si(001) is complicated by the possibility of Si diffusion to the Ge surface and the presence of threading dislocations. Herein, we demonstrate that similar to Ge(001), armchair nanoribbons with faceted edges and sub-10 nm widths can indeed be grown on Ge/Si(001). The nanoribbon synthesis proceeds faster than the diffusion of Si to the Ge surface if a sufficiently thick Ge epilayer is used (e.g., 3 μm), thereby avoiding competing reactions that form SiC. Moreover, threading dislocations in the Ge epilayer are not observed to affect the nucleation or anisotropic growth kinetics of the nanoribbons. Furthermore, these results demonstrate that previous successes in graphene nanoribbon synthesis on Ge(001) can be extended to Si(001) by using epilayers of Ge, motivating the future exploration of this promising platform for hybrid semiconducting graphene/Si electronics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

New mechanistic insight into the oxygen reduction reaction on Ruddlesden–Popper cathodes for intermediate-temperature solid oxide fuel cells

Ruddlesden–Popper (R–P) phase materials have been investigated widely as cathode candidates for IT-SOFCs. However, widespread application of R–P phase cathodes demands further improvement in electrode activity whose progress is hindered by the limited information in the oxygen reduction reaction (ORR). The ORR mechanism for the R–P phase is therefore investigated in this paper using (LaSr) 2 NiO 4±δ as an example. Accurate characterization of the surface oxygen exchange process is realized by developing thin and dense polycrystalline LSNO layers via a versatile spray-modified pressing method we invented before to avoid perceptible bulk diffusion contribution, surface enrichment and geometry complication. Hence, the governing factors of the ORR are identified as oxygen adsorption and incorporation based on the findings in reaction orders from electrochemical impedance spectroscopy (EIS), stoichiometry-related chemical capacitance and intrinsic anisotropic properties.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Layered structures of assembled imine-linked macrocycles and two-dimensional covalent organic frameworks give rise to prolonged exciton lifetimes

Ordered organic materials and assemblies have great potential to be tailored to have desirable properties for optoelectronic applications, such as long exciton lifetime and high directional exciton mobility. Framework materials, such as two-dimensional covalent organic frameworks (2D COFs), as well as their truncated macrocyclic analogues, are versatile platforms to organize functional aromatic systems into designed assemblies and robust materials. In this study we investigate the exciton dynamics in a 2D COF, its corresponding hexagonal macrocycle, and extended nanotubes comprised of stacked macrocycles. The excitonic behavior of these three systems provide an understanding of excitonic processes that occur in the plane of the covalently bonded 2D macromolecules and between layers of the nanotubes and 2D COF. The nanotube and analogous 2D COF exhibit longer excited-state lifetimes (~100 ps) compared to the individual, solvated macrocycles (<0.5 ps). These differences are attributed to the internal conversion facilitated by the internal motions of the imine linkages which are significantly reduced in the assembled macrocycles in the nanotube and 2D COF sheets in the layered structures. The exciton diffusion processes in the assembled nanotubes and 2D COF systems were characterized by the autocorrelations of the transition dipole moment of the excitons, giving the depolarization time constants for both systems to be ~1 ps. This work also reveals the anisotropic exciton dynamics related to the in-plane and inter-plane structural factors in these systems. These studies provide guidance for the design of future COF materials, where the longer excited state lifetimes imparted by assembly are beneficial for optoelectronic applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Lattice modification in KTiOPO{sub 4} by hydrogen and helium sequentially implantation in submicrometer depth

We investigated lattice modification and its physical mechanism in H and He co-implanted, z-cut potassium titanyl phosphate (KTiOPO{sub 4}). The samples were implanted with 110 keV H and 190 keV He, both to a fluence of 4 × 10{sup 16 }cm{sup −2}, at room temperature. Rutherford backscattering/channeling, high-resolution x-ray diffraction, and transmission electron microscopy were used to examine the implantation-induced structural changes and strain. Experimental and simulated x-ray diffraction results show that the strain in the implanted KTiOPO{sub 4} crystal is caused by interstitial atoms. The strain and stress are anisotropic and depend on the crystal's orientation. Transmission electron microscopy studies indicate that ion implantation produces many dislocations in the as-implanted samples. Annealing can induce ion aggregation to form nanobubbles, but plastic deformation and ion out-diffusion prevent the KTiOPO{sub 4} surface from blistering.

36 MATERIALS SCIENCE↗

Theory of plasma contactors in ground-based experiments and low Earth orbit

Previous theoretical work on plasma contactors as current collectors has fallen into two categories: collisionless double layer theory (describing space charge limited contactor clouds) and collisional quasineutral theory. Ground based experiments at low current are well explained by double layer theory, but this theory does not scale well to power generation by electrodynamic tethers in space, since very high anode potentials are needed to draw a substantial ambient electron current across the magnetic field in the absence of collisions (or effective collisions due to turbulence). Isotropic quasineutral models of contactor clouds, extending over a region where the effective collision frequency upsilon sub e exceeds the electron cyclotron frequency omega sub ce, have low anode potentials, but would collect very little ambient electron current, much less than the emitted ion current. A new model is presented, for an anisotropic contactor cloud oriented along the magnetic field, with upsilon sub e less than omega sub ce. The electron motion along the magnetic field is nearly collisionless, forming double layers in that direction, while across the magnetic field the electrons diffuse collisionally and the potential profile is determined by quasineutrality. Using a simplified expression for upsilon sub e due to ion acoustic turbulence, an analytic solution has been found for this model, which should be applicable to current collection in space. The anode potential is low and the collected ambient electron current can be several times the emitted ion current.

Gerver, M. J.↗

Combined energy and pitch angle diffusion of pickup ions at Comet Halley

It is well known that cometary pickup ions, e.g., H2O(+), OH(+), O(+), CO(+), H(+), initially form a ring-beam distribution in the solar wind reference frame, which is highly unstable to the growth of MHD waves (such as ion-cyclotron waves). The low-frequency magnetic fluctuations (or waves), which were observed upstream of Comet Halley, cannot only pitch-angle-scatter the pickup ions so that the distribution becomes at least partially isotropized, but also stochastically accelerate the ions, resulting in the energetic ion populations observed in the vicinity of Comet Halley. Here, numerical solutions of the quasi-linear diffusion equation were used to investigate the cometary ion pickup process at Comet Halley. Both pitch angle and energy diffusion are taken into account. Many quasi-linear models of cometary pickup ions exist which involve one type of diffusion or the other but not both types at once. It is found that the pitch angle scattering occurs faster than the energy diffusion, as expected. Moreover, the results demonstrate that the distribution of accelerated energetic ions is more isotropic than that of ions which have just been picked up. In fact, the ion distribution function on the initial pickup shell is quite anisotropic, even close to the Comet-Halley bow shock.

Ye, Gang↗

Dynamics of a reconnection-driven runaway ion tail in a reversed field pinch plasma

While reconnection-driven ion heating is common in laboratory and astrophysical plasmas, the underlying mechanisms for converting magnetic to kinetic energy remain not fully understood. Reversed field pinch discharges are often characterized by rapid ion heating during impulsive reconnection, generating an ion distribution with an enhanced bulk temperature, mainly perpendicular to magnetic field. In the Madison Symmetric Torus, a subset of discharges with the strongest reconnection events develop a very anisotropic, high energy tail parallel to magnetic field in addition to bulk perpendicular heating, which produces a fusion neutron flux orders of magnitude higher than that expected from a Maxwellian distribution. Here, we demonstrate that two factors in addition to a perpendicular bulk heating mechanism must be considered to explain this distribution. First, ion runaway can occur in the strong parallel-to-B electric field induced by a rapid equilibrium change triggered by reconnection-based relaxation; this effect is particularly strong on perpendicularly heated ions which experience a reduced frictional drag relative to bulk ions. Second, the confinement of ions varies dramatically as a function of velocity. Whereas thermal ions are governed by stochastic diffusion along tearing-altered field lines (and radial diffusion increases with parallel speed), sufficiently energetic ions are well confined, only weakly affected by a stochastic magnetic field. High energy ions traveling mainly in the direction of toroidal plasma current are nearly classically confined, while counter-propagating ions experience an intermediate confinement, greater than that of thermal ions but significantly less than classical expectations. The details of ion confinement tend to reinforce the asymmetric drive of the parallel electric field, resulting in a very asymmetric, anisotropic distribution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Connecting Collisionless Landau Fluid Closures to Collisional Plasma Physics Models

In this work, collisionless Landau fluid models are combined with collisional fluid models using a Padé approximation that is accurate in both collisionality limits. The collisionless models capture Landau damping through a nonlocal integro-differential sub-diffusion operator with ballistic characteristics. Collisional extensions of Landau fluid models are derived by analyzing the higher order moment equations which combine a Landau closure with collisional friction forces. The model derived here evolves fluid moments for density, parallel velocity, and anisotropic pressure and includes the frictional heat flux, the parallel thermal force and anisotropic electrical conductivity. Since anisotropies must vanish in the collisional limit, a simple closure can be derived if the friction force neglects pitch-angle scattering for the closure moments themselves. The resulting plasma physics model is potentially quite useful for applications in magnetic fusion and astrophysics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Thermal anisotropy enhanced by phonon size effects in nanoporous materials

While thermal anisotropy is a desirable materials property for many applications, including transverse thermoelectrics and thermal management in electronic devices, it remains elusive in practical natural compounds. In this work, we show how nanoporous materials with anisotropic pore lattices can be used as a platform for inducing strong heat transport directionality in isotropic materials. Using density functional theory and the phonon Boltzmann transport equation, we calculate the phonon-size effects and thermal conductivity of nanoporous silicon with different anisotropic pore lattices. Here, our calculations predict a strong directionality in the thermal conductivity, dictated by the difference in the pore-pore distances, i.e., the phonon bottleneck, along the two Cartesian axes. Using Fourier’s law, we also compute the diffusive heat transport for the same geometries obtaining significantly smaller anisotropy, revealing the crucial role of phonon-size effects in tuning thermal transport directionality. Besides enhancing our understanding of nanoscale heat transport, our results demonstrate the promise of nanoporous materials for modulating anisotropy in thermal conductivity.

36 MATERIALS SCIENCE↗

Reanalysis of the Apollo cosmic gamma-ray spectrum in the 0.3 to 10 MeV energy region

Additional data obtained from the Apollo 16 and 17 missions, together with collateral calculations on background radiation effects, have made possible an improved subtraction of unwanted backgrounds from the diffuse cosmic gamma ray data previously reported from Apollo 15. As a result, the 1 to 10 MeV spectrum is lowered significantly and connects smoothly with recent data at other energies. The inflection reported previously is much less pronounced and has no more than 1.5 sigma significance. Sky occultation by the Apollo 16 spacecraft shows the bulk of the 0.3 to 1 MeV radiation to be diffuse. The analysis of spurious backgrounds points to important improvements for future experiments designed for this spectral region. A light-weight satellite design can give a fourfold improvement in the signal to noise for such a measurement. Use of an anisotropic central crystal, which spins quickly compared with possible time variations in detector background, would enable sensitive limits to be set on galactic plane and point source contributions.

Dyer, C. S.↗

Interplay between out-of-plane anisotropic L1{sub 1}-type CoPt and in-plane anisotropic NiFe layers in CoPt/NiFe exchange springs

Films of L1{sub 1}-type CoPt/NiFe exchange springs were grown with different NiFe (Permalloy) layer thickness (t{sub NiFe} = 0–10 nm). X-ray diffraction analysis reveals that the characteristic peak position of NiFe(111) is not affected by the CoPt-layer—confirming the absence of any inter-diffusion between the CoPt and NiFe layers. Magnetic studies indicate that the magnetization orientation of NiFe layer can be tuned through varying t{sub NiFe} and the perpendicular magnetic anisotropy of L1{sub 1}-type CoPt/NiFe films cannot sustain for t{sub NiFe} larger than 3.0 nm due to the existence of exchange interaction at the interface of L1{sub 1}-CoPt and NiFe layers. Magnetic force microscopy analysis on the as-grown samples shows the changes in morphology from maze-like domains with good contrast to hazy domains when t{sub NiFe} ≥ 3.0 nm. The three-dimensional micro-magnetic simulation results demonstrate that the magnetization orientation in NiFe layer is not uniform, which continuously increases from the interface to the top of NiFe layer. Furthermore, the tilt angle of the topmost NiFe layers can be changed over a very wide range from a small number to about 75° by varying t{sub NiFe} from 1 to 10 nm. It is worth noting that there is an abrupt change in the magnetization direction at the interface, for all the t{sub NiFe} investigated. The results of present study demonstrate that the tunable tilted exchange springs can be realized with L1{sub 1}-type CoPt/NiFe bilayers for future applications in three-axis magnetic sensors or advanced spintronic devices demanding inclined magnetic anisotropy.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Formation of a conical distribution and intense ion heating in the presence of hydrogen cyclotron waves

In the considered investigation, it is assumed that the field aligned currents are responsible for producing electrostatic harmonic cyclotron waves (EHC). Using a one-dimensional simulation model in which the electron velocity distribution is maintained by a constant injection of the initial distribution, it is shown that, in contrast to earlier initial value simulations, EHC waves grow to a large amplitude, resulting in the formation of an anisotropic ion velocity distribution. Both the heating rate and the anisotropy are in reasonable agreement with the quasi-linear theory, taking into account the cyclotron resonance. The results show that the saturation is due to the combined effects of wave induced diffusion in an electron velocity space and the heating of ions perpendicularly. Both these effects reduce the growth rate.

Okuda, H.↗