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 109 records · Page 6

Buffered Coordination Modulation as a Means of Controlling Crystal Morphology and Molecular Diffusion in an Anisotropic Metal-Organic Framework

We know significant advances have been made in the synthesis of chemically selective environments within metal-organic frameworks, yet materials development and industrial implementation have been hindered by the inability to predictively control crystallite size and shape. One common strategy to control crystal growth is the inclusion of coordination modulators, which are molecular species designed to compete with the linker for metal coordination during synthesis. However, these modulators can simultaneously alter the pH of the reaction solution, an effect that can also significantly influence crystal morphology. Herein, noncoordinating buffers are used to independently control reaction pH during metal-organic framework synthesis, enabling direct interrogation of the role of the coordinating species on crystal growth. We demonstrate the efficacy of this strategy in the synthesis of low-dispersity single-crystals of the framework Co 2 (dobdc) (dobdc 4- = 2,5-dioxido-1,4-benzenedicarboxylate) in a pH 7-buffered solution using cobalt(II) acetate as the metal source. Density functional theory calculations reveal that acetate competitively binds to Co during crystallization, and by using a series of cobalt(II) salts with carboxylate anions of varying coordination strength, it is possible to control crystal growth along the c -direction. Finally, we use zero length column chromatography to show that crystal morphology has a direct impact on guest diffusional path length for the industrially important hydrocarbon m -xylene. Together, these results provide molecular-level insight into the use of modulators in governing crystallite morphology and a powerful strategy for the control of molecular diffusion rates within metal-organic frameworks.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Orientational Pair Correlations and Local Structure of Benzonitrile from Molecular Dynamics Simulations with Comparisons to Experiments [plus supplemental information]

Here, we present an experimentally parametrized molecular dynamics study of single-molecule and collective orientational relaxation in neat benzonitrile through the analysis of the reorientational anisotropy and polarizability anisotropy time-correlation function (PA-TCF). The simulations show that the PA-TCF is dominated by collective reorientation after 20 ps. Collective reorientation is found to be slower than single-molecule reorientation by a factor of 1.67, consistent with recent experiments. The simulations provide direct evidence of local antiparallel benzonitrile configurations. These structures, which have been the center of some debate, are responsible for the slower rate of collective vs. single-molecule reorientation in the liquid. Further structural analysis indicates that significant Coulombic interactions between the nitrile group and hydrogen atoms on adjacent molecules play a role in the formation of the antiparallel structures. The single-molecule dynamics reflected in the anisotropy are complex and consist of a ballistic regime, restricted angular diffusion, and spatially anisotropic free diffusion. The principal components of the rotational diffusion tensor are independently obtained and shown to reproduce the free diffusion regime of the anisotropy for each principal axis according to the predictions of a previous theory.

14 SOLAR ENERGY↗

Enhanced oxygen vacancy diffusion in Ta{sub 2}O{sub 5} resistive memory devices due to infinitely adaptive crystal structure

Metal oxide resistive memory devices based on Ta{sub 2}O{sub 5} have demonstrated high switching speed, long endurance, and low set voltage. However, the physical origin of this improved performance is still unclear. Ta{sub 2}O{sub 5} is an important archetype of a class of materials that possess an adaptive crystal structure that can respond easily to the presence of defects. Using first principles nudged elastic band calculations, we show that this adaptive crystal structure leads to low energy barriers for in-plane diffusion of oxygen vacancies in λ phase Ta{sub 2}O{sub 5}. Identified diffusion paths are associated with collective motion of neighboring atoms. The overall vacancy diffusion is anisotropic with higher diffusion barriers found for oxygen vacancy movement between Ta-O planes. Coupled with the fact that oxygen vacancy formation energy in Ta{sub 2}O{sub 5} is relatively small, our calculated low diffusion barriers can help explain the low set voltage in Ta{sub 2}O{sub 5} based resistive memory devices. Our work shows that other oxides with adaptive crystal structures could serve as potential candidates for resistive random access memory devices. We also discuss some general characteristics for ideal resistive RAM oxides that could be used in future computational material searches.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A signature of anisotropic cosmic-ray transport in the gamma-ray sky

A crucial process in Galactic cosmic-ray (CR) transport is the spatial diffusion due to the interaction with the interstellar turbulent magnetic field. Usually, CR diffusion is assumed to be uniform and isotropic all across the Galaxy. However, this picture is clearly inaccurate: several data-driven and theoretical arguments, as well as dedicated numerical simulations, show that diffusion exhibits highly anisotropic properties with respect to the direction of a background (ordered) magnetic field (i.e., parallel or perpendicular to it). In this paper we focus on a recently discovered anomaly in the hadronic CR spectrum inferred by the Fermi-LAT gamma-ray data at different positions in the Galaxy, i.e. the progressive hardening of the proton slope at low Galactocentric radii. We propose the idea that this feature can be interpreted as a signature of anisotropic diffusion in the complex Galactic magnetic field: in particular, the harder slope in the inner Galaxy is due, in our scenario, to the parallel diffusive escape along the poloidal component of the large-scale, regular, magnetic field. We implement this idea in a numerical framework, based on the DRAGON code, and perform detailed numerical tests on the accuracy of our setup. We discuss how the effect proposed depends on the relevant free parameters involved. Based on low-energy extrapolation of the few focused numerical simulations aimed at determining the scalings of the anisotropic diffusion coefficients, we finally present a set of plausible models that reproduce the behavior of the CR proton slopes inferred by gamma-ray data.

79 ASTRONOMY AND ASTROPHYSICS↗

Calculations of Diffuser Flows with an Anisotropic K-Epsilon Model

A newly developed anisotropic K-epsilon model is applied to calculate three axisymmetric diffuser flows with or without separation. The new model uses a quadratic stress-strain relation and satisfies the realizability conditions, i.e., it ensures both the positivity of the turbulent normal stresses and the Schwarz' inequality between any fluctuating velocities. Calculations are carried out with a finite-element method. A second-order accurate, bounded convection scheme and sufficiently fine grids are used to ensure numerical credibility of the solutions. The standard K-epsilon model is also used in order to highlight the performance of the new model. Comparison with the experimental data shows that the anisotropic K-epsilon model performs consistently better than does the standard K-epsilon model in all of the three test cases.

Zhu, J.↗

STOCHASTIC PARTICLE ACCELERATION IN TURBULENCE GENERATED BY MAGNETOROTATIONAL INSTABILITY

We investigate stochastic particle acceleration in accretion flows. It is believed that magnetorotational instability (MRI) generates turbulence inside accretion flows and that cosmic rays (CRs) are accelerated by the turbulence. We calculate equations of motion for CRs in the turbulent fields generated by MRI with the shearing box approximation and without back reaction to the field. Our results show that the CRs randomly gain or lose their energy through interaction with the turbulent fields. The CRs diffuse in the configuration space anisotropically: the diffusion coefficient in the direction of the unperturbed flow is about 20 times higher than the Bohm coefficient, while those in the other directions are only a few times higher than the Bohm. The momentum distribution is isotropic and its evolution can be described by the diffusion equation in momentum space where the diffusion coefficient is a power-law function of the CR momentum. We show that the shear acceleration works efficiently for energetic particles. We also cautiously note that in the shearing box approximation, particles that cross the simulation box many times along the radial direction undergo unphysical runaway acceleration by the Lorentz transformation, which needs to be taken into account with special care.

79 ASTRONOMY AND ASTROPHYSICS↗

Model atmospheres and radiation of magnetic neutron stars. I - The fully ionized case

Model neutron star atmospheres are calculated for typical cooling stars with a strong magnetic field and effective temperatures of 10 exp 5 to 10 exp 6 K. The effect of anisotropic photon diffusion in two normal modes are examined under the assumption that the opacity is due solely to the bremsstrahlung and Thomson scattering processes under conditions of LTE that are expected to prevail at the temperatures and densities obtained. The main aspects of anisotropic photon diffusion, and an original procedure for calculating model atmospheres and emitted spectra are discussed. Representative calculated spectra are given, and it is found that the hard spectral excess characterizing the nonmagnetic case, while still present, becomes less prominent in the presence of magnetic fields in the range of 10 exp 11 to 10 exp 13 G.

Shibanov, Iu. A.↗

Single molecule characterization of anomalous transport in a thin, anisotropic film

Here, the diffusion of small, charged molecules incorporated in an anisotropic polyelectrolyte multilayers (PEMs) was tracked in three dimensions by combining single-molecule fluorescence localization (to characterize lateral diffusion) with Förster resonance energy transfer (FRET) between diffusing molecules and the supporting surface (to measure diffusion in the surface-normal direction). Analysis of the surface-normal diffusion required model-based statistical analysis to account for the inherently noisy FRET signal. Combining these distinct single-molecule methods, which are inherently sensitive to different length scales, permitted simultaneous characterization of severely anisotropic diffusion, which was more than three orders of magnitude slower in the surface-normal direction. We hypothesize that the anomalously slow surface-normal diffusion was related to the periodic distribution of charge in the PEM, which created electrostatic barriers. The motion was strongly subdiffusive, with anomalous temporal scaling exponents in lateral and normal directions, suggesting a connection to the transient, random fractal conformation of polymer chains in the film’s matrix.

36 MATERIALS SCIENCE↗

Heat conduction in an irregular magnetic field. Part 2. Heat transport as a measure of the effective non-integrable volume

Given the large anisotropy of transport processes in magnetized plasmas, the magnetic field structure can strongly impact heat diffusion: magnetic surfaces and cantori form barriers to transport while chaotic layers and island structures can degrade confinement. When a small but non-zero amount of perpendicular diffusion is included, the structure of the magnetic field becomes less important, allowing pressure gradients to be supported across chaotic regions and island chains. We introduce a metric for the effective volume over which the local parallel diffusion dominates based on the solution to the anisotropic heat diffusion equation. To validate this metric, we consider model fields with a single island chain and a strongly chaotic layer for which analytic predictions of the relative parallel and perpendicular transport can be made. We also analyse critically chaotic fields produced from different sets of perturbations, highlighting the impact of the mode number spectrum on the heat transport. Our results indicate that this metric coincides with the effective volume of non-integrability in the limit $\kappa _{\perp } \rightarrow 0$ , where $\kappa_{\perp}$ is the perpendicular diffusion coefficient. We propose that this metric be used to assess the impact of non-integrability on the heat transport in stellarator equilibria.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Investigating the Effects of Anisotropic Mass Transport on Dendrite Growth in High Energy Density Lithium Batteries

Dendrite formation on the electrode surface of high energy density lithium (Li) batteries causes safety problems and limits their applications. Suppressing dendrite growth could significantly improve Li battery performance. Dendrite growth and morphology is a function of the mixing in the electrolyte near the anode interface. Most research into dendrites in batteries focuses on dendrite formation in isotropic electrolytes (i.e., electrolytes with isotropic diffusion coefficient). In this work, an anisotropic diffusion reaction model is developed to study the anisotropic mixing effect on dendrite growth in Li batteries. The model uses a Lagrangian particle-based method to model dendrite growth in an anisotropic electrolyte solution. The model is verified by comparing the numerical simulation results with analytical solutions, and its accuracy is shown to be better than previous particle-based anisotropic diffusion models. Several parametric studies of dendrite growth in an anisotropic electrolyte are performed and the results demonstrate the effects of anisotropic transport on dendrite growth and morphology, and show the possible advantages of anisotropic electrolytes for dendrite suppression.

SPH, high energy density lithium batteries, precip↗

An Efficient High-Order Solver for Diffusion Equations with Strong Anisotropy on Non-Anisotropy-Aligned Meshes

This paper concerns numerical solution of the diffusion equation with strong anisotropy on meshes not aligned with the anisotropic vector field. In order to resolve the numerical pollution for simulations on a non-anisotropy-aligned mesh and reduce the associated high computational cost we propose an effective preconditioner, extending our previous work. Similar to the anisotropy-aligned mesh case, we apply the auxiliary space preconditioning framework to design a preconditioner where a continuous finite element space is used as the auxiliary space for the discontinuous finite element space. The key component is an effective line smoother that can mitigate the high-frequency errors perpendicular to the magnetic field. We design a graph-based approach to find such a line smoother that is approximately perpendicular to the vector fields when the mesh does not align with the anisotropy. Finally, numerical experiments for several benchmark problems are presented, demonstrating the effectiveness and robustness of the proposed preconditioner when applied to Krylov iterative methods.

97 MATHEMATICS AND COMPUTING↗

Techniques for improved statistical convergence in quantification of eddy diffusivity moments

While recent approaches, such as the macroscopic forcing method (MFM) or Green's function-based approaches, can be used to compute Reynolds-averaged Navier-Stokes closure operators using forced direct numerical simulations, MFM can also be used to directly compute moments of the effective nonlocal and anisotropic eddy diffusivities. The low-order spatial and temporal moments contain limited information about the eddy diffusivity but are often sufficient for quantification and modeling of nonlocal and anisotropic effects. However, when using MFM to compute eddy diffusivity moments, the statistical convergence can be slow for higher-order moments. In this work, we demonstrate that using the same direct numerical simulation (DNS) for all forced MFM simulations improves statistical convergence of the eddy diffusivity moments. We present its implementation in conjunction with a decomposition method that handles the MFM forcing semianalytically and allows for consistent boundary condition treatment, which we develop for both scalar and momentum transport. We demonstrate that for a two-dimensional Rayleigh-Taylor instability case study, using the same DNS for all forced MFM simulations results in convergence with 𝒪⁡(100) simulations rather than 𝒪⁡(1000) simulations. In conclusion, we then demonstrate the impacts of improved convergence on the quantification of the eddy diffusivity.

general physics↗

Density functional theory calculations of self- and Xe diffusion in U 3 Si 2

Uranium silicide, U 3 Si 2 , has been proposed as an advanced nuclear fuel to be used in light water reactors (LWRs). Development of this alternative to the predominant current fuel, UO 2 , is motivated by enhanced accident tolerance as a result of higher thermal conductivity as well as improved fuel cycle economics through increased uranium density. In order to accurately model the fuel performance of U 3 Si 2 , the diffusion rate of point defects, which is related to self-diffusion, and of fission gas atoms must be determined. DFT calculations are used to predict the U and Si point defect concentrations, the corresponding self-diffusivities, the preferred Xe trap site and the Xe diffusivity. Effects of irradiation are not considered. A low defect formation energy and a high entropy for Si interstitials give rise to Si-rich non-stoichiometry at elevated temperatures. Both U and Si self-diffusion and Xe diffusion are anisotropic as a consequence of the tetragonal crystal structure of U 3 Si 2 . Si diffusion occurs by interstitial mechanisms in both the a-b plane and along the c axis, while the U c axis diffusion rate is controlled by a vacancy mechanism. Interstitial diffusion of U is very fast in the a-b plane of the U 3 Si 2 crystal structure. Xe atoms prefer to occupy U vacancy trap sites. The highest Xe diffusion rate occurs by a vacancy mechanism in both the a-b plane and along the c axis. The diffusion rate is similar in the a-b plane and along the c axis. U and Si self-diffusion and Xe diffusion are all faster in U3Si2 than intrinsic U and Xe diffusion in conventional UO 2 nuclear fuel.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Effects of surface vacancy-adatom pair formation on PFC tungsten surface morphological response

Abstract We report a simulation study of the effect of He-irradiation-induced surface vacancy-adatom pair formation on the surface morphological evolution of plasma-facing component (PFC) tungsten and examine a number of factors that impact such evolution. Our analysis is based on self-consistent dynamical simulations according to an atomistically-informed, continuum-scale surface evolution model that has been developed following a hierarchical multiscale modeling strategy and can access the spatiotemporal scales of relevance to fuzz formation. The model accounts for the flux of surface adatoms generated as a result of the surface vacancy-adatom pair formation effect upon He implantation, which contributes to the anisotropic growth of surface nanostructural features due to the different rates of adatom diffusion along and across step edges of islands on the tungsten surface. We have carried out atomic-scale computations of optimal diffusion pathways along and across island step edges on the W(110) surface and calculated Ehrlich–Schwoebel (ES) barriers in adatom diffusion along and across such step edges. This aspect of surface adatom diffusion contributes to anisotropic surface atomic fluxes, terrace and step diffusive currents, and has been incorporated into our PFC surface evolution model, which predicts the formation of preferentially aligned nanoridge stripe patterns on the PFC surface. We establish that these anisotropic diffusive currents accelerate nanotendril growth on the PFC surface and the onset of surface nanostructure pattern formation. We also explore systematically the dependence of the PFC surface morphological response on the surface temperature and He ion incident flux, characterize in detail the resulting surface topographies and growth kinetics, and compare the predicted surface morphologies with experimental observations. Our simulation predictions for the emerging surface nanostructure patterns under certain plasma exposure conditions are consistent with experimental findings in the literature.

36 MATERIALS SCIENCE↗

NON-AXISYMMETRIC PERPENDICULAR DIFFUSION OF CHARGED PARTICLES AND THEIR TRANSPORT ACROSS TANGENTIAL MAGNETIC DISCONTINUITIES

We investigate the transport of charged particles across magnetic discontinuities, focusing specifically on stream interfaces associated with co-rotating interaction regions in the solar wind. We argue that the magnetic field fluctuations perpendicular to the magnetic discontinuity, and usually also perpendicular to the mean magnetic field, are strongly damped in the vicinity of such a magnetic structure, leading to anisotropic perpendicular diffusion. Assuming that perpendicular diffusion arises from drifts in a turbulent magnetic field, we adopt a simplified approach to derive the relevant perpendicular diffusion coefficient. This approach, which we believe gives the correct principal dependences as expected from more elaborate calculations, allows us to investigate transport in different turbulent geometries, such as longitudinal compressional turbulence that may be present near the heliopause. Although highly dependent on the (possibly anisotropic) perpendicular length scales and turbulence levels, we generally find perpendicular diffusion to be strongly damped at magnetic discontinuities, which may in turn provide an explanation for the large particle gradients associated with these structures.

79 ASTRONOMY AND ASTROPHYSICS↗

Modulation of structure and dynamics of water under alternating electric field and the role of hydrogen bonding

Using Molecular Dynamics simulations, we investigate the effect of alternating (AC) electric field on static and dynamic properties of water. Here, the central question we address is how hydrogen bonds respond to perpetual field-induced dipole reorientations. We assess structural perturbations of water network and changes of hydrogen bond dynamics in a range of alternating electric field strengths and frequencies using a non-polarizable water model, SPC/E, and two distinct polarizable models: SWM4-NDP and BK3. We confirm that AC field causes only moderate structural perturbations. Dynamic properties, including the rates of bond breaking, switching of hydrogen-bonding partners, and diffusion, accelerate with the strength of AC fields. All models reveal a nonmonotonic frequency dependence with fastest dynamics at frequencies around 200 GHz where the period of the field oscillation is commensurate with the average time it takes a typical proton to switch from one acceptor to another. Higher frequencies result in smaller amplitudes of angle oscillations and in reduced probability to complete the switch to another acceptor before the field reversal restores the original configuration. As frequency increases, these effects gradually weaken the influence of the field on the kinetics of hydrogen bonding and the associated rates of translational and rotational diffusion in water.

74 ATOMIC AND MOLECULAR PHYSICS↗

Migration energy barriers and diffusion anisotropy of point defects on tungsten surfaces

We report tungsten is one of the most promising candidates for plasma-facing materials in future fusion devices, owing to its high performance under extreme irradiation conditions. However, irradiation-induced surface morphology varies significantly, depending on the irradiation type, fluence, and flux. Therefore, it is critical to examine the dynamics of point defects on tungsten (W) surfaces to understand how irradiation affects surface morphology. In this study, we employ the Self-Evolving Atomistic Kinetic Monte Carlo (SEAKMC) method to search for potential migration paths of point defects on W (1 0 0), (1 1 0) and (1 1 1) surfaces. The first-principles calculations are then used to accurately determine the migration energy barriers. The obtained paths and barriers are incorporated into a kinetic Monte Carlo (KMC) model to determine trajectories and diffusivities, which are described by diffusion tensors to demonstrate their anisotropic features. Multiple diffusion mechanisms have been identified on different surfaces, with various anisotropy factors. Particularly, point defects on the W (1 1 0) surface have the highest diffusivities, with anisotropy factors independent of temperatures. In comparison, the anisotropy of the W (1 0 0) surface decreases as temperature increases, while the W (1 1 1) surface is isotropic for point defects. This study provides insights into defect transport properties on different surfaces, which are essential for understanding the early stages of irradiation-induced microstructural evolutions and surface morphology of tungsten.

36 MATERIALS SCIENCE↗

Ultra‐Narrow Phosphorene Nanoribbons Produced by Facile Electrochemical Process

Abstract Phosphorene nanoribbons (PNRs) have inspired strong research interests to explore their exciting properties that are associated with the unique two‐dimensional (2D) structure of phosphorene as well as the additional quantum confinement of the nanoribbon morphology, providing new materials strategy for electronic and optoelectronic applications. Despite several important properties of PNRs, the production of these structures with narrow widths is still a great challenge. Here, a facile and straightforward approach to synthesize PNRs via an electrochemical process that utilize the anisotropic Na + diffusion barrier in black phosphorus (BP) along the [001] zigzag direction against the [100] armchair direction, is reported. The produced PNRs display widths of good uniformity (10.3 ± 3.8 nm) observed by high‐resolution transmission electron microscopy, and the suppressed B 2g vibrational mode from Raman spectroscopy results. More interestingly, when used in field‐effect transistors, synthesized bundles exhibit the n‐type behavior, which is dramatically different from bulk BP flakes which are p‐type. This work provides insights into a new synthesis approach of PNRs with confined widths, paving the way toward the development of phosphorene and other highly anisotropic nanoribbon materials for high‐quality electronic applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗