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 469 records · Page 26

Enhanced mobility of cations and anions in the redox state: The polaronium mechanism

Hugely enhanced slow-ion diffusivity has been widely observed under extreme redox conditions and for unclear reasons. Aided by first-principles calculations on model systems of ZrO 2 , CeO 2 , BaTiO 3 and Li 4/3 Mn 2/3 O 2 , here we successfully explained the intriguing phenomenon by a polaronium mechanism. We found a polaronium, defined as a transitory complex of a polaronic electron or hole and a migrating counterion, becomes highly mobile when the counterion comes from a d 0 or f 0 cation (e.g., Zr 4+ , Ti 4+ and Ce 4 ) or a p 6 anion (e.g., O 2- ) in the host compound. Upon a redox reaction, the complex attains a d 1 /f 1 or p 5 configuration, which spontaneously forms because it is favored by an electron-phonon interaction (manifest as the Jahn-Teller effect in high symmetry systems) that enables local relaxation and lowers the system energy. Here our calculations found such interaction reaching its peak at the saddle point where the local environment is softest, so soft that it allows a reorientation of the anisotropic d/f/p orbital to minimize the electron repulsion locally. Since the complex may dissolve after a successful ion-migration event, the redox electron/hole can be recycled to form another free-radical-like polaronium elsewhere, thereby enhancing ion migration repeatedly. The proposed polaronium mechanism, which also operates in ceramics doped with mixed-valence cations, is most relevant under dynamic and extreme thermal/field/irradiation conditions where extra electrons/holes are abundantly generated by non-equilibrium redox reactions. For such operations, some with emerging applications, our diffusion-enhancing mechanism may provide new theoretical insight to help understand their material/microstructure stability and performance.

36 MATERIALS SCIENCE↗

Tracking Dendritic Growth in Hydrogen-Based Hematite Reduction via Computer Vision

The reduction of hematite to metallic iron using hydrogen (H2) as a reducing agent presents a promising pathway for decarbonizing steel production. In this study, we employ a combination of in situ confocal scanning laser microscopy (CSLM) and advanced computer vision techniques to quantitatively analyze dendritic growth of ferrite during H2-based reduction of iron oxide at high temperatures. A workflow integrating Watershed Image Segmentation (WIS) and Lucas-Kanade Optical Flow (LKOF) is developed to extract both global and local kinetic information from time-resolved micrograph sequences. H2 reduction experiments conducted at 1400 degrees C and 1500 degrees C demonstrate a clear correlation between temperature and reduction rate, as evidenced by accuracy of fitted Johnson-Mehl-Avrami-Kolmogorov (JMAK) parameters. Optical flow analysis further elucidates the anisotropic and branched nature of dendritic growth, providing spatially resolved velocity fields that correlate well with global transformation kinetics. The proposed methodology demonstrates strong agreement with experimental measurements and literature values, offering a robust framework for automated image-based analysis to study kinetics through microstructural evolution in the reduction of iron ore, and likely other reaction-diffusion phenomena.

08 HYDROGEN↗

Turbulence-driven anisotropic electron tail generation during magnetic reconnection

Magnetic reconnection (MR) plays an important role in particle transport, energization, and acceleration in space, astrophysical, and laboratory plasmas. In the Madison Symmetric Torus reversed field pinch, discrete MR events release large amounts of energy from the equilibrium magnetic field, a fraction of which is transferred to electrons and ions. Previous experiments revealed an anisotropic electron tail that favors the perpendicular direction and is symmetric in the parallel. New profile measurements of x-ray emission show that the tail distribution is localized near the magnetic axis, consistent modeling of the bremsstrahlung emission. The tail appears first near the magnetic axis and then spreads radially, and the dynamics in the anisotropy and diffusion are discussed. The data presented imply that the electron tail formation likely results from a turbulent wave-particle interaction and provides evidence that high energy electrons are escaping the core-localized region through pitch angle scattering into the parallel direction, followed by stochastic parallel transport to the plasma edge. Finally, new measurements also show a strong correlation between high energy x-ray measurements and tearing mode dynamics, suggesting that the coupling between core and edge tearing modes is essential for energetic electron tail formation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Effect of Radiation on Chromospheric Magnetic Reconnection: Reactive and Collisional Multi-fluid Simulations

We study magnetic reconnection under chromospheric conditions in five different ionization levels from 0.5% to 50% using a self-consistent two-fluid (ions + neutrals) model that accounts for compressibility, collisional effects, chemical inequilibrium, and anisotropic heat conduction. Results with and without radiation are compared, using two models for the radiative losses: an optically thin radiation loss function, and an approximation of the radiative losses of a plasma with photospheric abundances. The results without radiation show that reconnection occurs faster for the weakly ionized cases as a result of the effect of ambipolar diffusion and fast recombination. The tearing mode instability appears earlier in the low ionized cases and grows rapidly. We find that radiative losses have a stronger effect than was found in previous results as the cooling changes the plasma pressure and the concentration of ions inside the current sheet. This affects the ambipolar diffusion and the chemical equilibrium, resulting in thin current sheets and enhanced reconnection. The results quantify this complex nonlinear interaction by showing that a strong cooling produces faster reconnections than have been found in models without radiation. The results accounting for radiation show timescales and outflows comparable to spicules and chromospheric jets.

79 ASTRONOMY AND ASTROPHYSICS↗

Atomic to molecular hydrogen transition in interstellar clouds

The photodestruction of molecular hydrogen is reexamined in the context of recent observational data. Improvements in the theory of radiation transfer are made by using a more accurate self-shielding integral, treating anisotropic scattering by dust grains, and including photodissociation from excited rotational states. The total column density of hydrogen required to obtain 10% conversion decreases with a parameter, epsilon, the ratio of the formation rate to the destruction rate, as epsilon to the -m power, with m in the range 1.4-1.6. Satisfactory agreement with the data for thicker diffuse clouds is obtained for epsilon of about 0.00006.

Federman, S. R.↗

Twin Proliferation and Prolongation under Kinetic Control: Pd–Au Janus Icosahedra versus Pd@Au Core–Shell Starfishes

Heterogeneous bimetallic nanocrystals featuring explicit spatial configurations and abundant twin defects can simultaneously enable geometric and ligand effects to enhance catalytic and photonic applications. Herein, we report two growth patterns of Au atoms on penta-twinned Pd decahedra, involving twin proliferation to generate asymmetric Pd–Au Janus icosahedra and twin elongation to produce anisotropic Pd@Au core–shell starfishes, respectively. Mechanistic analysis indicates that the injection rate determines the lower-limit number (n low ) of Au(III) ions in the steady state and thus controls the growth pattern. When n low ≤ 5.5, the kinetic rate is slow enough to initiate asymmetrical one-side growth but fast enough to outpace surface diffusion; Au tetrahedral subunits are successively proliferated along the axial $\langle$110$\rangle$ direction of Pd decahedra to form Pd–Au Janus icosahedra. Composed of five Pd and 15 Au tetrahedral subunits, such a heterogeneous icosahedron supports high (2.2 GPa) tensile strain and high strain difference up to +21.9%. In contrast, when n low > 5.5, the fast reduction kinetics promotes symmetric growth with inadequate surface diffusion. As such, Au atoms are laterally deposited along five high-indexed $\langle$211$\rangle$ ridges of Pd decahedra to generate concave Pd@Au core–shell starfishes with tunable sizes (28–40 nm), twin elongation ratios (33.82–162.08%), and lattice expansion ratios (8.82–20.10%).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Role of Magnetic Reconnection in MHD Turbulence

The current understanding of MHD turbulence envisions turbulent eddies which are anisotropic in all three directions. In the plane perpendicular to the local mean magnetic field, this implies that such eddies become current-sheet-like structures at small scales. We analyze the role of magnetic reconnection in these structures and conclude that reconnection becomes important at a scale ... This scale is larger than the scale set by the resistive diffusion of eddies, therefore implying a fundamentally different route to energy dissipation than that predicted by the Kolmogorov-like phenomenology. In particular, our analysis predicts the existence of the sub-inertial, reconnection interval of MHD turbulence ...

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Annual Research Briefs

This report contains the 1997 annual progress reports of the research fellows and students supported by the Center for Turbulence Research (CTR). Titles include: Invariant modeling in large-eddy simulation of turbulence; Validation of large-eddy simulation in a plain asymmetric diffuser; Progress in large-eddy simulation of trailing-edge turbulence and aeronautics; Resolution requirements in large-eddy simulations of shear flows; A general theory of discrete filtering for LES in complex geometry; On the use of discrete filters for large eddy simulation; Wall models in large eddy simulation of separated flow; Perspectives for ensemble average LES; Anisotropic grid-based formulas for subgrid-scale models; Some modeling requirements for wall models in large eddy simulation; Numerical simulation of 3D turbulent boundary layers using the V2F model; Accurate modeling of impinging jet heat transfer; Application of turbulence models to high-lift airfoils; Advances in structure-based turbulence modeling; Incorporating realistic chemistry into direct numerical simulations of turbulent non-premixed combustion; Effects of small-scale structure on turbulent mixing; Turbulent premixed combustion in the laminar flamelet and the thin reaction zone regime; Large eddy simulation of combustion instabilities in turbulent premixed burners; On the generation of vorticity at a free-surface; Active control of turbulent channel flow; A generalized framework for robust control in fluid mechanics; Combined immersed-boundary/B-spline methods for simulations of flow in complex geometries; and DNS of shock boundary-layer interaction - preliminary results for compression ramp flow.

Spinks, Debra↗

Verification of Unstructured Grid Adaptation Components

Adaptive unstructured grid techniques have made limited impact on production analysis workflows where the control of discretization error is critical to obtaining reliable simulation results. Recent progress has matured a number of independent implementations of flow solvers, error estimation methods, and anisotropic grid adaptation mechanics. Known differences and previously unknown differences in grid adaptation components and their integrated processes are identified here for study. Unstructured grid adaptation tools are verified using analytic functions and the Code Comparison Principle. Three analytic functions with different smoothness properties are adapted to show the impact of smoothness on implementation differences. A scalar advection-diffusion problem with an analytic solution that models a boundary layer is adapted to test individual grid adaptation components. Laminar flow over a delta wing and turbulent flow over an ONERA M6 wing are verified with multiple, independent grid adaptation procedures to show consistent convergence to fine-grid forces and a moment. The scalar problems illustrate known differences in a grid adaptation component implementation and a previously unknown interaction between components. The wing adaptation cases in the current study document a clear improvement to existing grid adaptation procedures. The stage is set for the infusion of verified grid adaptation into production fluid flow simulations.

Park, Michael A.↗

Verification of Unstructured Grid Adaptation Components

Adaptive unstructured grid techniques have made limited impact on production analysis workflows where the control of discretization error is critical to obtaining reliable simulation results. Recent progress has matured a number of independent implementations of flow solvers, error estimation methods, and anisotropic grid adaptation mechanics. Known differences and previously unknown differences in grid adaptation components and their integrated processes are identified here for study. Unstructured grid adaptation tools are verified using analytic functions and the Code Comparison Principle. Three analytic functions with different smoothness properties are adapted to show the impact of smoothness on implementation differences. A scalar advection-diffusion problem with an analytic solution that models a boundary layer is adapted to test individual grid adaptation components. The scalar problems illustrate known differences in a grid adaptation component implementation and a previously unknown interaction between components. Laminar flow over a delta wing is verified with multiple, independent grid adaptation procedures to show consistent convergence to fine-grid forces and pitching moment.

Park, Michael A.↗

Signatures of low-dimensional magnetism and short-range spin fluctuations in Co-based trirutiles

Features of low dimensional magnetism resulting from a square-net arrangement of Co atoms in trirutile CoTa2O6 is studied in the present work by means of density functional theory and is compared with the experimental results of specific heat and neutron diffraction. The small total energy differences between the ferromagnetic (FM) and antiferromagnetic (AFM) configuration of CoTa2O6 shows that competing magnetic ground states exist, with the possibility of transition from FM to AFM phase at low temperature. Our calculation further suggests the semi-conducting behavior for CoTa2O6 with a band gap of ~0.41 eV. The calculated magnetic anisotropy energy is ~2.5 meV with its easy axis along the [100] (in-plane) direction. Studying the evolution of magnetism in Co1-xMgxTa2O6 (x = 0, 0.1, 0.3, 0.5, 0.7 and 1). it is found that the sharp AFM transition exhibited by CoTa2O6 at TN = 6.2 K in it’s heat capacity vanishes with Mg-dilution, indicating the obvious effect of weakening the superexchange pathways of Co. The current specific heat study reveals the robust nature of TN for CoTa2O6 in applied magnetic fields. Clear indication of short-range magnetism is obtained from the magnetic entropy, however, diffuse components are absent in neutron diffraction data. At TN , CoTa2O6 enters a long-range ordered magnetic state which can be described using a propagation vector, ( 1 1 0). Upon Mg-dilution at x =0.1, the long-range ordered magnetism is destroyed. The present results should motivate an investigation of magnetic excitations in this low-dimensional anisotropic magnet.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

First Protein Crystallization Experiments on The International Space Station: Sweet Success in Space With Thaumatin

We determined the room temperature 1.2 A structure of thaumatin using a crystal grown in the first protein crystallization experiment conducted aboard the International Space Station (ISS). The crystals were grown in the Enhanced Gaseous Nitrogen Dewar (EGN) developed by Alexander McPherson and co-workers. EGN transports frozen solutions contained in tygon tubing in a liquid nitrogen Dewar to ISS where the tubes then thaw. Batch, free interface diffusion (FID), or vapor diffusion crystallization occurs after thawing. EGN was flown to the ISS on STS-106 on September 8, 2000. This was a "risk mitigation" flight that tested EGN performance and the process of conducting experiments on ISS. We focused on how to map a hanging drop crystallization recipe to the EGN FID method. Thaumatin was chosen as the test system. Three series of crystallization recipes were set-up. Each series tested different volume ratios of protein-rich solution to precipitant-rich solution. The series differed from each other by fixing either the protein concentration or the amount of protein in the solutions. Upon return of the samples to Earth on October 24 by STS-92, bubbles that spanned the diameter of the tubing were observed in all tubes. Such bubbles interrupt liquid-liquid diffusion and force vapor diffusion equilibration to occur instead. Nonetheless, crystals grew in 9 of 30 tubes. Many large crystals were grown, the largest being 2.0 x 1.1 x 1.0 cubic mm. The largest crystal was used to collect data at room temperature on beamline 7-1 of the Stanford Synchrotron Radiation Source to a maximum resolution of 1.2 A. The structure was refined anisotropically using SHELX with a data to parameter ratio of 4.5 to give an R(sub factor) of 15.8% (R(sub free) = 18.2%) for ail reflections without generated hydrogens. This refinement is proceeding. Comparisons of this 1.2 A microgravity structure to previous reports of the thaumatin structure at 1.75 A and to ground control crystals will be presented.

Kundrot, Craig E.↗

Fourier-like Thermal Relaxation of Nanoscale Explosive Hot Spots

Hot spots are local regions of high temperature that are widely considered to govern explosive initiation. Hot spot dynamics rests on a delicate balance between heat generation due to chemical reactions and heat loss through thermal conduction, making accurate determinations of the conductivity under extreme conditions a key component of predictive explosive models. We develop here an approach to directly determine the thermal transport properties of explosive hot spots with realistic initial structures through a combination of molecular dynamics (MD) and diffusive heat equation (HEq) modeling. Effective thermal conductivity values are determined by fitting HEq models to MD predictions of long timescale hot spot relaxation. The approach is applied to model hot spots in the molecular crystalline explosive 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) for a range of shock strengths and two limiting cases for impact orientation. Isotropic and anisotropic HEq models yield similar results, despite TATB exhibiting some of the largest and most anisotropic thermal conductivity values for explosive near normal conditions. The conductivity is found to be a strong function of density, which parametrically captures dependence on temperature, pressure, and material state. The associated root-mean-square errors of the fitted HEq models are approximately 5% of MD predicted final equilibrium temperatures. Here, the conductivity values determined here for TATB hot spots are considerably larger than those used in a prior hot spot criticality study, which may significantly impact predictions for critical hot spot sizes. The approach provides a convenient foundation for determining the effective thermal conductivity for hot spot problems in other explosives and directly yields information on reasonable approximations that might be taken in higher-level models for those materials.

36 MATERIALS SCIENCE↗

Electrochemical Degradation of Pt–Ni Nanocatalysts: An Identical Location Aberration-Corrected Scanning Transmission Electron Microscopy Study

The evolution of Pt–Ni nanoparticles supported on amorphous carbon is investigated before and after electrochemical potential cycling (0.6-1.1 V), using aberration-corrected scanning transmission electron microscopy (STEM) and energy dispersive X-ray spectroscopy (EDS). During voltage cycling and due to the dissolution of nanoparticles, single ions/atoms and ionic/atomic clusters emerge and diffuse across the carbon support toward larger nanoparticles, where they redeposit. We observe that the preferred locations for the dissolution are the steps and corners of the nanoparticles. On the other hand, the redeposition process happens often on {111} type planes. In addition, contrary to the conventional view, where larger particles grow isotropically from smaller ones, this research work shows that anisotropic growth of smaller particles occurs during potential cycling. The reason for this behavior seems to be related to the fact that smaller particles with thicker Pt-rich shells trigger the nucleation and deposition of Pt.

catalysts↗

Modeling Coupled Chemo-Mechanical Behavior of Randomly Oriented NMC811 Polycrystalline Li-Ion Battery Cathodes

This paper develops a three-dimensional, transient, chemo-mechanical model that predicts the performance of single secondary particle Li-ion battery cathodes. The secondary particles are composed of numerous (approximately 60) randomly oriented single-crystal primary particles. The model incorporates concentration-dependent and anisotropic material properties. As much as possible, electrochemical, transport, and structural properties for crystalline NMC811 (Li x Ni 0.8 Mn 0.1 Co 0.1 O 2 ) are taken from prior publications. Weak Van der Waals bonding between primary particles is modeled empirically using a spring analogy, which enables local primary-particle separations (disintegration) and subsequent reattachments. The model fully couples Li diffusion and the mechanical response. Results include predictions of local Li-concentrations and stresses. High stresses are found near grain boundaries, especially when the lattice orientations are greatly misaligned. Particle separations are characterized in terms of a damage parameter. The model is used to predict the effects of design and operating conditions, including charge/discharge rates, cycling scenarios, and particle sizes.

25 ENERGY STORAGE↗

Non-random walk diffusion enhances the sink strength of semicoherent interfaces

Clean, safe and economical nuclear energy requires new materials capable of withstanding severe radiation damage. One strategy of imparting radiation resistance to solids is to incorporate into them a high density of solid-phase interfaces capable of absorbing and annihilating radiation-induced defects. Here we show that elastic interactions between point defects and semicoherent interfaces lead to a marked enhancement in interface sink strength. Our conclusions stem from simulations that integrate first principles, object kinetic Monte Carlo and anisotropic elasticity calculations. Surprisingly, the enhancement in sink strength is not due primarily to increased thermodynamic driving forces, but rather to reduced defect migration barriers, which induce a preferential drift of defects towards interfaces. The sink strength enhancement is highly sensitive to the detailed character of interfacial stresses, suggesting that ‘super-sink’ interfaces may be designed by optimizing interface stress fields. Lastly, such interfaces may be used to create materials with unprecedented resistance to radiation-induced damage.

36 MATERIALS SCIENCE↗

Interface defects in SiC power MOSFETs - An electrically detected magnetic resonance study based on spin dependent recombination

This study presents electrically detected magnetic resonance (EDMR) measurements on a silicon carbide (SiC) MOSFET having the structure of a double-diffused silicon MOSFET (DMOS). The resonance pattern of a SiC DMOS was measured by monitoring the change of the recombination current between the source/body and the drain. The amplitude of the response has a maximum when the device is biased in depletion due to the equal concentrations of electrons and holes at the interface resulting in the most efficient recombination. The measured anisotropic g-tensor has axial symmetry with g{sub ∥} = 2.0051(4) (B ‖ c-axis), and g{sub ⊥} = 2.0029(4) (B⊥ c-axis) and the pattern shows several hyperfine (HF) peaks. We tentatively identify the observed defect as a silicon vacancy located directly at the interface.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗