Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffusion coefficient”

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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Modeling the Size Distribution and Chemical Composition of Secondary Organic Aerosols during the Reactive Uptake of Isoprene-Derived Epoxydiols under Low-Humidity Condition

Reactive uptake of isoprene epoxydiols (IEPOX), which are isoprene oxidation products, onto acidic sulfate aerosols is recognized to be an important mechanism for the formation of isoprene-derived secondary organic aerosol (SOA). While a mechanistic understanding of IEPOX-SOA formation exists, several processes affecting their formation remain uncertain. Evaluating mechanistic IEPOX-SOA models with controlled laboratory experiments under longer atmospherically relevant timescales is critical. Here, we implement our latest understanding of IEPOX-SOA formation within a box model to simulate the measured reactive uptake of IEPOX on polydisperse ammonium bisulfate seed aerosols within an environmental Teflon chamber. The model is evaluated with single-particle measurements of size distribution, volume, density, and composition of aerosols due to IEPOX-SOA formation at timescales of hours. Additionally, we find that the model can simulate the growth of particles due to IEPOX multiphase chemistry, as reflected in increases of the mean particle size and volume concentrations, and a shift of number size distribution to larger sizes. The model also predicts the observed evolution of particle number mean diameter and total volume concentrations at the end of the experiment. We show that in addition to self-limiting effects of IEPOX-SOA coatings, accounting for the molar balance between inorganic and organic sulfate and mass accommodation coefficient of IEPOX are important parameters governing modeling of IEPOX-SOA formation. Thus, models which do not account for the molar sulfate balance and/or diffusion limitations within IEPOX-SOA coatings are likely to predict too high IEPOX-SOA formation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Integrating Ultra-Coarse-Grained Protein Models into Accessible Workflows for Multiscale Molecular Dynamics

To capture protein conformational transitions using molecular dynamics (MD), several simulation resolutions covering different spatial and temporal scales are typically needed. All-atom (AA) simulations provide fine resolution, but are computationally infeasible for large systems over longer durations. Coarse-grained (CG) and ultra-coarse-grained (UCG) models have a lower resolution and computational cost while still being able to conserve essential protein features. Prior work on a Multiscale Machinelearned Modeling Infrastructure (MuMMI) combined both AA and CG simulations to study RAS-RAF protein interactions, leveraging CG models for longer time scales and using AA to investigate unusual conformations in greater detail. However, MuMMI is still resource-intensive, and this study aims to maximize exploration of the protein conformational space while reducing computational cost. In this paper, we build on prior work that integrates UCG models based on heterogeneous elastic network modeling (hENM) into the MuMMI workflow. We demonstrate that UCG models enable accurate sampling of protein conformations, focusing on simulating RAS-RAF protein interactions. Using higher-resolution CG Martini simulation data, we can automatically refine intramolecular interactions in UCG models. We present a scalable Python package that uses fluctuations observed in higher-resolution CG Martini simulations to estimate bond coefficients of the UCG model. We built novel machine learning-based backmapping methods to recover more detailed CG Martini structures from UCG structures, using diffusion models to learn the mapping between scales. Finally, we present UCG-mini-MuMMI, an accessible and less compute-intensive version of MuMMI as a resource for the scientific community. Incorporating UCG models into MD studies is applicable to a broad range of systems and proteins, and our study offers insights into the advantages and limitations of these methods.

Chemical structure↗
DOE OSTI2022

Study of quasi-collisional effects in laboratory and astrophysical plasmas

High-amplitude turbulence excited in plasmas at small-scales can leave imprint on the radiation produced by (accelerated) plasma particles -- mostly electrons in electron-ion plasmas and both electrons and positrons in lepton pair plasmas. Furthermore, turbulence is known to introduce "effective collisional" effects -- anomalous resistivity, dissipation, diffusion, etc. -- in the otherwise collisionless plasmas. These effective collisions affect radiation transfer and modify optical and magneto-optic effects -- the transmittance and reflectivity coefficients, Faraday rotation, etc. Note that the effective collisions considered in this project are not resonant wave-particle interactions, but much less conventional randomization of particles’ paths. Therefore, we colloquially refer them to as "quasi-collisions". The study of quasi-collisions on transport, radiative and optical properties of plasmas is the main focus of this project.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Distinguishing Surface and Bulk Reactivity: Concentration-Dependent Kinetics of Iodide Oxidation by Ozone in Microdroplets

Iodine oxidation reactions play an important role in environmental, biological, and industrial contexts. The multiphase reaction between aqueous iodide and ozone is of particular interest due to its prevalence in the marine atmosphere and unique reactivity at the air–water interface. Here, we explore the concentration dependence of the I – + O 3 reaction in levitated microdroplets under both acidic and basic conditions. To interpret the experimental kinetics, molecular simulations are used to benchmark a kinetic model, which enables insight into the reactivity of the interface, the nanometer-scale subsurface region, and the bulk interior of the droplet. For all experiments, a kinetic description of gas- and liquid-phase diffusion is critical to interpreting the results. We find that the surface dominates the iodide oxidation kinetics under concentrated and acidic conditions, with the reactive uptake coefficient approaching an upper limit of 10 –2 at pH 3. In contrast, reactions in the subsurface dominate under more dilute and alkaline conditions, with inhibition of the surface reaction at pH 12 and an uptake coefficient that is 10× smaller. In conclusion, the origin of a changing surface mechanism with pH is explored and compared to previous ozone-dependent measurements.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

TOUGHREACT V4.12

TOUGHREACT V4.12 is a 3-D reactive transport code that handles to non-isothermal, multiphase/multicomponent fluid flow, heat transport, aqueous and gaseous species advection-diffusion, and equilibrium/kinetic water-gas-rock-biological reactions. TOUGHREACT V4.12 is based on TOUGHREACT V3.3, with several new features including: an option to compute activity coefficients using the Pitzer ion interaction model; implementation of ECO2N V2.0 for temperatures up to 300 C; options to simulate forward and reverse osmosis, including 2D Poiseuille flow; addition of EOS8 and revisions for EOS5 and EWASG with H2(g); addition of gas species decay, various new output formats and time step control options; many bug fixes.

Spycher, Nicolas↗
DOE OSTI2025

Mesoscale Linear Elastic Modeling and Homogenization of Marine Energy Composites

The design of fiber-reinforced composite (FRC)-based components for marine energy applications necessitates a fundamental understanding of material properties and the resulting geometry to predict long-term performance. In this work, we present a modeling workflow to predict linear elastic and diffusive bulk properties at the mesoscale for an idealized geometry based on knowledge of fiber and resin properties. A parametric study was performed to identify the key model input parameters that influence bulk properties. Furthermore, we demonstrate how bulk properties can be leveraged in high-fidelity image-based simulations, where imperfections in tow geometry and voids captured during X-ray computed tomography imaging are explicitly represented within the simulation. Bulk properties of interest include moduli, Poisson’s ratios, hygroscopic swelling, diffusivity, and moisture uptake, which are key parameters for characterizing FRC performance within marine environments. Modeling predictions agreed well with experimental data, except for estimating swelling coefficients, likely due to crack accumulation as a function of moisture uptake. The mesoscale modeling workflow ultimately highlights a versatile framework for understanding the influence of material and geometric properties, which can be leveraged to rapidly assess new FRC-based components.

computational mechanics↗
DOE OSTI2025

Explosive Soot Challenge (Final Report)

This project assembled a broad ensemble of modeling and experimentation tools to study the morphological and optical properties of detonation soots in explosive fireballs. A gram-scale hemispherical high explosive was studied in a low-pressure controlled environment using in-situ experimentation with diffusely illuminated visible absorption spectroscopy, particle sizing through light scattering techniques, and post-test collections with subsequent morphological analysis. Hydrocode modeling was performed to replicate the detonation flow observations, and subsequent aerosol kinetics models provided particle size distributions and extinction coefficients from the hydrocode results. Experimentally observed soot morphologies agreed with expectation from the literature - a bimodal distribution was found, brought upon by the particles growing to a size where their inertia and fluid wakes are non-negligible. The aerosol kinetics model did not replicate the observed bimodal size distribution for lack of a coagulation kernel to represent the behavior. To recover particulate optical properties, a spectrally resolved absorption spectroscopy method termed Spectral diffuse back-illuminated extinction imaging (SBI-EI) was developed and implemented on two explosive types. Inverting the absorption spectra using a Kramers-Kronig consistent method yielded the complex index of refraction for the soots produced by the explosives. This method resulted in an unrealistic index of refraction for one of the two explosives, and this is suspected to be due to the model neglecting scattering brought upon by the large particle sizes observed. In addition to the core work, three additional studies were performed in parallel. These investigated the impact of scattering on diffuse absorption spectroscopy, studied how soots oxidate and sublimate in a well-controlled shock tube, and laid the theoretical groundwork for a new collision kernel to replicate the bimodal size distribution from the observations. Summaries of these efforts are included at the end of this report.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2025

Microstructure characteristics of LPBF&HIP fabricated graded composite transition joint between ferritic steel and austenitic stainless steel

Graded composite transition joints (GCTJs) offer a promising alternative to conventional dissimilar metal welds (DMWs) by enabling smooth compositional and microstructural transitions. However, GCTJs fabricated solely through additive manufacturing (AM) face challenges such as heat accumulation, complex parameter control, and elemental segregation. In this study, we propose a novel approach that relies on AM to design a spatially graded structure in one alloy and then employs hot isostatic pressing (HIP) as a diffusion bonding method to join it with a second alloy. Here, this method combines the flexibility of AM with the powder net-shaping advantage of HIP. Specifically, a series of closely packed austenitic stainless steel 304 conical structures were printed using laser powder bed fusion (LPBF) and then combined with ferritic steel P91 powder via HIP. By using electron backscatter diffraction (EBSD), electron probe microanalysis (EPMA), and transmission electron microscopy (TEM) techniques, the microstructure characteristics of the GCTJ of 304&P91, especially the interdiffusion zone (IDZ), have been systematically investigated. The microstructure at the interface transitions from austenite-ferrite (A+F) to austenite-martensite-ferrite (A+M+F), and finally to martensite-ferrite (M+F) due to diffusion. Additionally, the diffusion width between 304 and P91 increases with the volume fraction of P91. This unique design also ensures a gradual transition in both hardness and thermal expansion coefficient from 304 to P91, thereby enabling a smooth gradient in functional properties. Overall, this study proposes a novel approach for fabricating GCTJs and contributes to advancing design concepts in the field of dissimilar metal joining.

Additively manufacturing (AM)↗
DOE OSTI2023

Irradiation Effects on Stability of δ-UZr2 phase in U-50 wt% Zr Alloy

U-50wt%Zr is a candidate metallic nuclear fuel with potential application in light water reactors due to its excellent thermal properties and high radiation tolerance. The Zr-rich UZr fuels possess greater swelling resistance and fission gas release characteristics compared with U-rich UZr fuels. In this current study, the δ-phase U-50wt%Zr is proton irradiated at room temperature to 1 displacement per atom (dpa) to provide insights on phase stability under irradiation conditions. High resolution characterization of Transmission Electron Microscopy (TEM) and Atom Probe Tomography (APT) characterization techniques are used to elucidate microstructural changes due to irradiation. TEM and APT results show highly oriented bcc β-Zr-rich platelet precipitates nucleating adjacent to α-U phases inside the UZr 2 matrix. Formation of this platelet morphology is characteristic of Widmänstatten structure which can be attributed to a variety of factors such as differences in thermal expansion coefficient between the phases, grain size, alloy composition, and cooling rate. The phases present are distinctly different than those observed through in situ annealing, but irradiation accelerates diffusion and phase separation kinetics. These microstructural changes in U-50wt%Zr are different from those achieved by pure thermodynamic or high temperature heavy ion irradiation experiments. Finally, our work, together with previous ones, highlight the necessity to study the U-Zr phase diagram under non-equilibrium thermodynamics conditions to support this material's deployment as a viable nuclear fuel form.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Coupled Interfacial and Bulk Kinetics Govern the Timescales of Multiphase Ozonolysis Reactions

Chemical transformations in aerosols impact the lifetime of particle phase species, the fate of atmospheric pollutants, and both climate- and health-relevant aerosol properties. Timescales for multiphase reactions of ozone in atmospheric aqueous phases are governed by coupled kinetic processes between the gas phase, the particle interface, and its bulk, which respond dynamically to reactive consumption of O 3 . However, models of atmospheric aerosol reactivity often do not account for the coupled nature of multiphase processes. Herein, to examine these dynamics, we use new and prior experimental observations of aqueous droplet reaction kinetics, including three systems with a range of surface affinities and ozonolysis rate coefficients ( trans -aconitic acid (C 6 H 6 O 6 ), maleic acid (C 4 H 4 O 4 ), and sodium nitrite (NaNO 2 )). Using literature rate coefficients and thermodynamic properties, we constrain a simple two-compartment stochastic kinetic model which resolves the interface from the particle bulk and represents O 3 partitioning, diffusion, and reaction as a coupled kinetic system. Our kinetic model accurately predicts decay kinetics across all three systems, demonstrating that both the thermodynamic properties of O 3 and the coupled kinetic and diffusion processes are key to making accurate predictions. An enhanced concentration of adsorbed O 3 , compared to gas and bulk phases is rapidly maintained and remains constant even as O 3 is consumed by reaction. Multiphase systems dynamically seek to achieve equilibrium in response to reactive O 3 loss, but this is hampered at solute concentrations relevant to aqueous aerosol by the rate of O 3 arrival in the bulk by diffusion. As a result, bulk-phase O 3 becomes depleted from its Henry's law solubility. This bulk-phase O 3 depletion limits reaction timescales for relatively slow-reacting organic solutes with low interfacial affinity (i.e., trans-aconitic and maleic acids, with k rxn ≈ 10 3 -10 4 M -1 s -1 ), which is in contrast to fast-reacting solutes with higher surface affinity (i.e., nitrite, with k rxn ≈ 10 5 M -1 s -1 ) where surface reactions strongly impact the observed decay kinetics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Self-Trapped-Exciton Radiative Recombination in β–Ga 2 O 3 : Impact of Two Concurrent Nonradiative Auger Processes

The peculiarities of radiative and nonradiative processes associated with self-trapped intrinsic eXcitons in the excited β-Ga 2 O 3 crystals are studied via time-resolved techniques of induced absorption, transient grating, and photoluminescence (PL) at room temperature. The excitation above the bandgap is produced by laser pulses with linear light polarization parallel and orthogonal in the (–201) and (001) planes. We elucidate that the nonradiative recombination rate occurring in the eXciton prevails over its radiative emission rate in a wide range of free carrier concentration composed of excited and equilibrium electrons. Hence, the nonradiative recombination has no effect on the strong anisotropy and the shape of the eXciton emission band. However, we find out that the conventional ABC model of electron effective lifetime is insufficient for explanation of the excitation dependences. Inclusion of two nonradiative Auger mechanisms in a modified ABC formula provides excellent agreement of these dependences. We conclude that the trap-assisted Auger process is in proportion to the free electron density with coefficient B = 1.1 × 10 –11 cm 3 /s and appears at low/intermediate excitation, while the triple-particle Auger process is in proportion to Δn 2 with coefficient C = 8 × 10 –30 cm 6 /s and appears at high excitation conditions. The transition between two Auger mechanisms is accompanied by a rise of the eXciton diffusivity in preferred crystallographic directions where the radiative PL intensity is maximal. The diffusion length LD in these directions can reach values ~300 nm, but, at high excitations, LD becomes limited by Auger lifetimes. These findings pave the way for the implementation of self-trapped eXcitons into specific optoelectronic devices.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Dependence of the impurity transport on the dominant turbulent regime in ELM-y H-mode discharges on the DIII-D tokamak

Laser blow-off injections of aluminum and tungsten have been performed on the DIII-D tokamak to investigate the variation of impurity transport in a set of dedicated ion and electron heating scans with a fixed value of the external torque. The particle transport is quantified via the Bayesian inference method, which, constrained by a combination of a charge exchange recombination spectroscopy, soft X-ray measurements, and VUV spectroscopy provides a detailed uncertainty quantification of transport coefficients. Contrasting discharge phases with a dominant electron and ion heating reveal 3-fold drop in the impurity confinement time and order of magnitude increase in midradius impurity diffusion, when additional electron heating is applied. Further, the calculated stationary aluminum density profiles reverse from peaked in electron heated to hollow in the ion heated case, following a similar trend as electron and carbon density. Comparable values of a core diffusion have been observed for W and Al ions, while differences in the propagation dynamics of these impurities are attributed to pedestal and edge transport. Modeling of the core transport with non-linear gyrokinetics code CGYRO, significantly underpredicts the magnitude of the variation in Al transport. Diffusion increases 3-times steeper with additional electron heat flux and 10-times lower diffusion is observed in ion heated case than predicted by the modeling. Here, the CGYRO model quantitatively matches the increase of the Al diffusion when approaching the linear threshold for the transition from the ion temperature gradient (ITG) to trapped electron mode (TEM).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Experimental inference of neutral and impurity transport in Alcator C-Mod using high-resolution x-ray and ultra-violet spectra

We present experimental inferences of cross-field impurity transport coefficients for AlcatorC-Mod plasmas without edge-localized modes, leveraging a novel forward model for the entire Ca Kα spectrum, including satellite lines within the spectral range, to compare to high-resolution x-ray imaging crystal spectroscopy (XICS). These measurements are complemented by extreme ultra-violet (EUV) spectroscopy that constrains transport closer to the edge. Using new atomic data sets for both XICS and EUV analysis has enabled consideration of line ratios across both spectral ranges and has increased the accuracy of inferred transport coefficients. Inclusion of charge exchange between edge thermal neutrals and impurities is shown to be extremely important in C-Mod pedestals. Atomic D neutral densities from experimental D Lyα measurements at the midplane are compared to SOLPS-ITER simulations, finding good agreement. Bayesian inferences of impurity transport coefficients are presented for L-, EDA H-, and I-mode discharges, making use of the Aurora package for forward modeling and combining our spectroscopic constraints. Experimentally inferred diffusion profiles are found to match turbulent transport models at midradius within uncertainties, using both quasilinear gyro-fluid TGLF SAT-1 and nonlinear ion-scale gyrokinetic CGYRO simulations. Furthermore, significant discrepancies in convection are observed in some cases, suggesting difficulties in predictions of flat or hollow impurity profiles.

neutrals↗
DOE OSTI2022

Synthesis, Crystal Structure, and Physical Properties of BaSnS 2

Phase-pure BaSnS 2 , with space group P2 1 /c, is synthesized, and the structural and physical properties are investigated. We report thermal properties and optical measurements are reported for the first time. The Debye temperature and Sommerfeld coefficient are obtained from temperature-dependent heat capacity measurements, the latter indicating that BaSnS 2 is an electrical insulator. A direct bandgap of 2.4 eV is obtained from diffuse reflectance and photoluminescence spectroscopy. The findings herein lay the foundation for understanding the physical properties of this material and are part of a continuing effort to investigate previously unexplored ternary chalcogenides.

semiconductors↗
DOE OSTI2023

Synthesis of high-entropy carbides from multi-metal polymer precursors

Three compositions of high-entropy carbides, (TiHfVNbTa)C, (TiZrHfNbTa)C and (TiZrNbTaW)C were synthesized via a modified Pechini process, in which citric acid served as both a cation chelating agent and a carbon source. Through pyrolysis and spark plasma sintering, single phase high-entropy carbides were formed from the homogeneous precursors at a relatively low temperature of 1800 °C. The dispersion of cations in the polymer precursor facilitated shorter diffusion distances in polymer-derived materials, and thus compositional homogeneity was significantly improved relative to materials produced by a solid-state method, as quantified by a defined coefficient of variation applied to energy-dispersive X-ray spectroscopy. Finally, a finer microstructure in polymer-derived materials results in improved fracture toughness with a K IC value of 4.29 MPa·m 1/2 achieved for (TiHfVNbTa)C.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Prediction of divertor heat flux width for ITER pre-fusion power operation using BOUT++ transport code

Prediction of divertor heat flux width is performed for the first and the second pre-fusion power operation (PFPO) phases specified in the new ITER research plan using BOUT++ transport code (Li et al 2018 Comput. Phys. Commun.228 69–82). Here the initial plasma profiles inside the separatrix are taken from CORSICA scenario studies. Transport coefficients in transport code are calculated by inverting the plasma profiles inside the separatrix and are assumed to be constants in the scrape-off-layer. An anomalous thermal diffusivity scan is performed with E × B and magnetic drifts. The results in two scenarios identifying two distinct regimes: a drift-dominant regime when diffusivity is smaller than the respective critical thermal diffusivity χ c and a turbulence-dominant regime when diffusivity is larger than it. The Goldston heuristic drift model and the ITPA multi-machine experimental scaling yield a lower limit of the width λ q . From transport simulations, we obtain the critical thermal diffusivity χ c = 0.5 m 2 s –1 for the PFPO-1 scenario with toroidal magnetic field B = 1.77 T and plasma current I p = 5 MA, and χ c = 0.3 m 2 s –1 for the PFPO-2 scenario with toroidal magnetic field B = 2.65 T and plasma current I p = 7.5 MA. Separatrix temperature and collisionality also have a significant impact on the heat flux width in the drift-dominant regime. The investigation clearly yields a scaling for critical thermal diffusivity ${\chi }_{\text{c}}\propto {A}^{1/2}/(Z{\left(1+Z\right)}^{\frac{1}{2}}{B}_{\text{p}}^{2})$ using ITER scenarios with fixed safety factor q 95 , major radius R, aspect ratio R/a, and the separatrix temperature T sep , and establishes the connection with CFETR and C-Mod discharges. This scaling implies that for a given tokamak device with q 95 , R, R/a, and T sep fixed, a reduction of poloidal magnetic field by a factor of 3 leads to a 9 times higher critical value of thermal diffusivity χ c , possibly yielding a transition from turbulence- to drift-dominant regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Phase-field model of freeze casting

Directional solidification of water-based solutions has emerged as a versatile technique for templating hierarchical porous materials. However, the underlying mechanisms of pattern formation remain incompletely understood. In this work, we present a detailed derivation and analysis of a quantitative phase-field model for simulating this nonequilibrium process. The phase-field model extends the thin-interface formulation of dilute binary alloy solidification with antitrapping to incorporate the highly anisotropic energetic and kinetic properties of the partially faceted ice-water interface. This interface is faceted in the basal plane normal to the ⟨0001⟩ directions and atomically rough in other directions within the basal plane. On the basal plane, the model reproduces a linear or nonlinear relationship between the interface growth rate and the kinetic undercooling that can be linked to experimental measurements. In both cases, spontaneous parity breaking of the solidification front is observed when the preferred growth direction is aligned with the temperature gradient. This phenomenon leads to the formation of partially faceted ice lamellae that drift laterally in one of the ⟨0001⟩ directions. Here, we demonstrate that the drifting velocity of the ice lamellae is controlled by the kinetics on the basal plane and converges as the thickness of the diffuse solid-liquid interface decreases. Furthermore, we examine the effect of the form of the kinetic anisotropy, which is chosen here such that the inverse of the kinetic coefficient varies linearly from a finite value in the ⟨0001⟩ directions to zero in all other directions within the basal plane, consistent with the assumption that the interface grows in local thermodynamic equilibrium in this plane. Our results indicate that the drifting velocity of ice lamellae is not affected by the slope of this linear relation, and the radius and undercooling at the tip of an ice lamella converge at relatively small slope values. Consequently, the phase-field simulations remain quantitative with computationally tractable choices of both the interface thickness and the slope assumed in the form of the kinetic anisotropy.

Materials science↗