Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “magnetic 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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Dynamic heterogeneity and nanophase separation in rubber-toughened amine-cured highly cross-linked polymer networks

Solid state nuclear magnetic resonance (NMR) spectroscopy and small-to wide-angle X-ray scattering (SWAXS) methods were used to characterize the heterogeneous dynamics and polymer domain structure in rubber modified thermoset materials containing the diglycidyl ether of bisphenol A (DGEBA) epoxy resin and a mixture of Jeffamine reactive rubber and 4,4-diaminodicyclohexylmethane (PACM) amine curing agent. The polymer chain dynamics and morphologies as a function of the PACM/Jeffamine ratio were determined. Using dipolar-filtered NMR experiments, the resulting networks are shown to be composed of mobile and rigid regions that are separated on nanometer length scales, along with a dynamically immobilized interface region. Proton NMR spin diffusion experiments measured the dimensions of the mobile phase to range between 9 and 66 nm and varied with the relative PACM concentration. Solid state 13 C magic angle spinning NMR experiments show that the highly mobile phase is composed entirely of the dynamically flexible polyether chains of the Jeffamine rubber, the immobilized interface region is a mixture of DGEBA, PACM, and the Jeffamine rubber, with the PACM cross-linked to DGEBA predominantly residing in the rigid phase. The SWAXS results showed compositional nanophase separation spanning the 11–77 nm range. These measurements of the nanoscale compositional and dynamic heterogeneity provide molecular level insight into the very broad and controllable glass transition temperature distributions observed for these highly cross-linked polymer networks.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Machine Learning assisted optimization and parameter space exploration dataset of spin ice Hamiltonian

This repository contains both simulated and experimental structure factor data for the data challenge involving the inverse scattering problem. The simulated data were generated during a machine-learning-assisted optimization routine described in ref[1]. The experimental structure factor was measured on a rare-earth oxide, Dy2Ti2O7 using diffuse neutron scattering from time-of-flight techniques on the CORELLI instrument at the Spallation Neutron Source, Oak Ridge National Laboratory. A Metropolis Monte Carlo code implemented to run in a High-performance computing setting was used to calculated simulated structure factors for the spin-ice Hamiltonian at 680 mK, which is the same temperature as for the experimental data. The total size of all the files in this repository is 5.12 GB. A detailed description of the files is given below. ExperimentalData_630mK.dat – A linearized version of 3-dimensional experimental data of size 61×81×21. This data was processed to remove an estimation of non-magnetic background, including nuclear scattering signal and instrumentation background. Parameters.dat – 6700 samples were evaluated over the 4-dimensional parameter space (J_1, J_2, J_3 and J_(3^' )). There is an additional parameter, D in the spin Hamiltonian to mimic the dipolar interaction between magnetic ions. However, this parameter, D was fixed to a value determined by prior work. This file contains five columns for the parameters J_1, J_2, J_3, J_(3^' ) and D respectively. 3D_Simulation_Data.dat – The simulated structure factor, S(Q) data are included in this file. Each raw contains a linearized array of 3D volumes of S(Q) calculated for the parameter set given in the corresponding row of the file Parameters.dat. The size of the volume data was matched to the experimental data. Qx(h,-h,0).dat, Qy(k,k,-2k).dat, Qz(l,l,l).dat – These files contain the h, k, and l values along with the reciprocal vectors [h,-h,0], [k,k,-2k] and [l,l,l] respectively.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Performance and Economic Analysis of Organosolv Softwood and Herbaceous Lignins to Activated Carbons as Electrode Materials in Supercapacitors

In this work, yellow pine (YP, softwood) and switchgrass (SG, grass) lignins were extracted as a coproduct of an organosolv γ-valerolactone (GVL) biorefinery that also produces biofuels and furfural. The extracted lignins were converted to carbon precursors for synthesizing porous activated carbon electrodes for high energy-density supercapacitors. This research details the impact of lignin composition on the derived porous structures and electrochemical properties of activated carbons. Lignin precursors with various syringyl (S) to guaiacyl (G) contents were characterized using 31 P nuclear magnetic resonance (NMR) and two-dimensional 1 H‒ 13 C NMR. A two-step activation process, using steam and carbon dioxide as the activating agents, enabled the formation of porous carbons structures with high surface area. The capacitive behavior of supercapacitors was systematically characterized by cyclic voltammetry, charge-discharge cycling, and electrochemical impedance spectroscopy. The specific capacitance of YP and SG capacitors reached 367 and 221 F g −1 , respectively. Both types of capacitors demonstrated remarkably stable capacitance (capacitance retention >90%) along with excellent Coulombic efficiency (>99%) over 10,000 cycles. Compared to SG electrode, the better electrochemical performance achieved with YP electrodes was mainly due to shorter diffusion path, improved ionic mobility, and increased active surface area. The inexpensive lignin-based porous electrodes synthesized in this work can be used for various electrochemical devices for improved performance, decreased cost, and enhanced durability. This work also demonstrates that the selection of feedstock and appropriate processing conditions can tailor the structure of carbon composites for targeted applications. Techno-economic analysis indicates that YP and SG activated carbons can be produced at a minimum selling price of $8,493 and $6,670 per ton, respectively, which is competitive with the commercially available supercapacitor-grade activated carbons.

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

Observations and Simulations of Dropout Events and Flux Decays in October 2013: Comparing MEO Equatorial With LEO Polar Orbit

We compare ESA PROBA-V observations of electron flux at LEO with those from the NASA Van Allen Probes mostly at MEO for October 2013. Dropouts are visible at all energy during four storms from both satellites. Equatorially trapped electron fluxes are higher than at LEO by 10 2 (<1 MeV) to 10 5 (>2.5 MeV). We observe a quite isotropic structure of the outer belt during quiet times, contrary to the inner belt, and pitch angle dependence of high energy injection. We find a very good overlap of the outer belt at MEO and LEO at ~0.5 MeV. We use test-particle simulations of the energetic electrons trapped in the terrestrial magnetic field to study the outer radiation belt electron flux changes during geomagnetic storms. We show that the Dst (Disturbance storm time) effect during the main phase of a geomagnetic storm results in a betatron mechanism causing outward radial drift and a deceleration of the electrons. This outward drift motion is energy independent, pitch angle-dependent, and represents a significant distance (~1 L-shell at L = 5 for moderate storms). At fixed L-shell, this causes a decay of the LEO precipitating flux (adiabatic outward motion), followed by a return to the normal state (adiabatic inward motion) during main and recovery phases. Dst effect, associated with magnetopause shadowing and radial diffusion can explain the main characteristics of outer radiation belt electron dropouts in October 2013. We also use Fokker-Planck simulations with event-driven diffusion coefficients at high temporal resolution, to distinguish instantaneous loss from the gradual scattering that depopulates the slot region and the outer belt after storms. Simulations reproduce the slot formation and the gradual loss in the outer belt. The typical energy dependence of these losses leads to the absence of scattering for relativistic and ultra-relativistic electrons in the outer belt, oppositely to dropouts.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Optical Detection of Long Electron Spin Transport Lengths in a Monolayer Semiconductor

Using a spatially resolved optical pump-probe experiment, we measure the lateral transport of spin-valley polarized electrons over very long distances (tens of micrometers) in a single WSe 2 monolayer. By locally pumping the Fermi sea of 2D electrons to a high degree of spin-valley polarization (up to 75%) using circularly polarized light, the lateral diffusion of the electron polarization can be mapped out via the photoluminescence induced by a spatially separated and linearly polarized probe laser. In this work, up to 25% spin-valley polarization is observed at pump-probe separations up to 20 μm. Characteristic spin-valley diffusion lengths of 18 ± 3 μm are revealed at low temperatures. The dependence on temperature, pump helicity, pump intensity, and electron density highlight the key roles played by spin relaxation time and pumping efficiency on polarized electron transport in monolayer semiconductors possessing spin-valley locking.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Advances in physics and applications of 3D magnetic perturbations on the J-TEXT tokamak

In the last two years, three major technical improvements have been made on J-TEXT in supporting of the expanded operation regions and diagnostic capabilities. (1) The successful commission of the 105 GHz/500 kW/1 s electron cyclotron resonance heating (ECRH) system increasing the core electron temperature from 0.9 keV up to around 1.5 keV. (2) The poloidal divertor configuration with an X-point in the high-field side has been achieved. In particular, the 400 kW electron cyclotron wave has also been successfully injected into the diverted plasma. (3) A 256-channel electron cyclotron emission imaging diagnostic system and two sets of four-channel Doppler backscattering diagnostics have been successfully developed on J-TEXT, allowing detailed measurement of the electron temperature and density fluctuations for turbulence and MHD research. The locked mode (LM), especially the 2/1 LM, is one of the biggest threats to the plasma operation. Both the thresholds of 2/1 and 3/1 LM are observed to vary non-monotonically on electron density. The electrode biasing was applied successfully to unlock the LM from either a rotating or static resonant magnetic perturbation (RMP) field. In the presence of 2/1 LM, three kinds of standing wave (SW) structures have been observed to share a similar connection to the island structure, i.e. the nodes of the SWs locate around the O- or X-points of the 2/1 island. The control and mitigation of disruption is essential to the safe operation of ITER, and it has been systematically studied by applying a RMP field, massive gas injection (MGI) and shattered pellet injection on J-TEXT. When the RMP-induced 2/1 LM is larger than a critical width, the MGI shutdown process can be significantly influenced. If the phase difference between the O-point of LM and the MGI valve is +90° (or –90°), the penetration depth and the assimilation of impurities can be enhanced (or suppressed) during the pre-thermal quench (TQ) phase and result in a faster (or slower) TQ. A secondary MGI can also suppress the runaway electron (RE) generation, if the additional high-Z impurity gas arrives at the plasma edge before TQ. When the secondary MGI has been applied after the formation of the RE current plateau, the RE current can be dissipated, and the dissipation rate increases with the injected impurity quantity but saturates with a maximum of 28 MA s –1 . The non-local transport is experimentally observed in the ion transport channel. The electron thermal diffusivity significantly increases with the ECRH power. Theoretical work shows that significant intrinsic current can be driven by electromagnetic turbulence, and the robust formation mechanism of the E × B staircase is identified from the Hasegawa–Wakatani system.

divertor↗
DOE OSTI2022

Gyrokinetic simulations of momentum flux parasitic to free-energy transfer

Ion Landau damping interacts with a portion of the E×B drift to cause a non-diffusive outward flux of co-current toroidal angular momentum. Quantitative evaluation of this momentum flux requires nonlinear simulations to determine fL, the fraction of fluctuation free energy that passes through ion Landau damping, in fully developed turbulence. Nonlinear gyrokinetic simulations with the GKW code confirm the presence of the systematic symmetry-breaking momentum flux. For simulations with adiabatic electrons, fL scales inversely with the ion temperature gradient, because only the ion curvature drift can transfer free energy to the electrostatic potential. Although kinetic electrons should in principle relax this restriction, the ion Landau damping measured in collisionless kinetic-electron simulations remained at low levels comparable with ion-curvature-drift transfer, except when magnetic shear was strong. A set of simulations scanning the electron pitch-angle scattering rate showed only a weak variation of fL with the electron collisionality. However, collisional-electron simulations with electron temperature greater than ion temperature unambiguously showed electron-curvature-drift transfer supporting ion Landau damping, leading to a corresponding enhancement of the symmetry-breaking momentum flux.

intrinsic rotation↗
DOE OSTI2023

Coherences of Photoinduced Electron Spin Qubit Pair States in Photosystem I

Here, this publication presents the first comprehensive experimental study of electron spin coherences in photosynthetic reaction center proteins, specifically focusing on photosystem I (PSI). The ultrafast electron transfer in PSI generates spin-correlated radical pairs (SCRPs), which are entangled spin pairs formed in well-defined spin states (Bell states). Since their discovery in our group in the 1980s, SCRPs have been extensively used to enhance our understanding of structure–function relationships in photosynthetic proteins. More recently, SCRPs have been utilized as tools for quantum sensing. Electron spin decoherence poses a significant challenge in realizing practical applications of electron spin qubits, particularly the creation of quantum entanglement between multiple electron spins. This work is focused on the systematic characterization of decoherence in SCRPs of PSI. These decoherence times were measured as electron spin echo decay times, termed phase memory times (T M ), at various temperatures. Decoherence was recorded on both transient SCRP states P 700 + A 1 – and thermalized states. Our study reveals that T M exhibits minimal dependence on the biological species, biochemical treatment, and paramagnetic species. The analysis indicates that nuclear spin diffusion and instantaneous diffusion mechanisms alone cannot explain the observed decoherence. As a plausible explanation we discuss the assumption that the low-temperature dynamics of methyl groups in the protein surrounding the unpaired electron spin centers is the main factor governing the loss of the spin coherence in PSI.

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

Hydrodynamics of quantum spin liquids

Quantum spin liquids are topological states of matter that arise in frustrated quantum magnets at low temperatures. At low energies, such states exhibit emergent gauge fields and fractionalized quasiparticles and can also possess enhanced global symmetries compared to their parent microscopic Hamiltonians. Here we study the consequences of this emergent gauge and symmetry structure for the hydrodynamics of quantum spin liquids. Specifically, we analyze two cases, the U(1) spin liquid with a Fermi surface and the SU(4)-symmetric “algebraic” spin liquid. We show that the emergent degrees of freedom in the spin liquid phase lead to a variety of additional hydrodynamic modes compared to the high-temperature paramagnetic phase. We identify a hydrodynamic regime for the internal U(1) gauge field common to both states, characterized by slow diffusion of the internal transverse photon.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Unconventional short-range structural fluctuations in cuprate superconductors

Abstract The interplay between structural and electronic degrees of freedom in complex materials is the subject of extensive debate in physics and materials science. Particularly interesting questions pertain to the nature and extent of pre-transitional short-range order in diverse systems ranging from shape-memory alloys to unconventional superconductors, and how this microstructure affects macroscopic properties. Here we use neutron and X-ray diffuse scattering to uncover universal structural fluctuations in La 2-x Sr x CuO 4 and Tl 2 Ba 2 CuO 6+δ , two cuprate superconductors with distinct point disorder effects and with optimal superconducting transition temperatures that differ by more than a factor of two. The fluctuations are present in wide doping and temperature ranges, including compositions that maintain high average structural symmetry, and they exhibit unusual, yet simple scaling behaviour. The scaling regime is robust and universal, similar to the well-known critical fluctuations close to second-order phase transitions, but with a distinctly different physical origin. We relate this behaviour to pre-transitional phenomena in a broad class of systems with structural and magnetic transitions, and propose an explanation based on rare structural fluctuations caused by intrinsic nanoscale inhomogeneity. We also uncover parallels with superconducting fluctuations, which indicates that the underlying inhomogeneity plays an important role in cuprate physics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

LA-R43S6-L1 Postshot Report [Slides]

The Ranchero LA-R43S6-L1 test (formerly referred to as LA-43S-L1. The R specifically designates it is a Ranchero device and the 6 indicates a 6” high explosive charge diameter in anticipation of using larger diameters in the future.) was conducted on 3/8/17 with good success. The experiment was the first test of a “Swooped” Ranchero flux compression generator (FCG), and the load for the test was an aluminum imploding liner. In addition to the swoop, the stator was anodized to provide an insulating layer instead of one or more layers of polyethylene as had been the case for all earlier Ranchero tests. Using the Firing point TA-39-88 capacitor bank, 3.5 MA was delivered to the FCG for the initial magnetic field, and 36.4 MA were delivered to the load having an initial inductance of 2.5 nH. Implosion speeds over 1.1 cm/us were recorded, which exceeds previous LANL solid liner experimental results. The test was diagnosed with 18 PDV channels, 12 of which recorded liner implosion data, and the other 6 of which tracked the FCG armature during flux compression. These armature expansion data recorded during flux compression are the first of their kind, and it was previously unknown, due to unknown effects of the SF 6 in the flux compression volume, whether or not the expansion could be recorded using the PDV technique. Post shot 2D MHD calculations employed a complete external circuit model, and an improved flux diffusion model not available prior to the test. Results from these simulations show very good agreement with the experimental result. Four PDV channels recorded the liner implosion viewed radially outward from the center of the liner. All four of these probes recorded implosion velocities greater that 1 cm/µs, with the cylindrical center displaced by ~2 mm from actual center. Four PDV channels looked at +10° angles from the CMU wall, and four looked at -10° angles. A glide plane interaction is shown moving toward the center of the liner in MHD calculations, and PDV probes provide confirmation. The high explosive (HE) in the LA-R43S6-L1 test was PBX 9501 which had to be glued together in many pieces, and the test was preceded by a camera test to verify that the PBX 9501 could be assembled with acceptable tolerance in glue joints to prevent severing the armature during expansion. The camera test verified an acceptable armature performance, and results are given in complete detail in a post shot report LA-UR-19-20124, and summarized here. This report provides complete detail of the considerable body of data obtained on the test and in the post shot calculations. It is prepared in Power Point for ease of preparation and future review. Shot documentation available in the LANL on-line library are included as references, and non-referenceable documents will be cited and stored in an LA- R43S6-L1 post shot folder stored on an M-6 shared drive. Fabrication drawings are also maintained in an M-6 shared drive and paper files are available from M-6 personnel. This report will be maintained as a PowerPoint document on the shared drive, as well, since included movies will play and graphs from Excel files retain information in the PowerPoint versions. PDF versions are required for clearance and are not to large to send by e-mail.

42 ENGINEERING↗
DOE OSTI2025

Spatial profile of argon (1s 5 ) metastables in an electron beam generated plasma

Electron beams with an applied magnetic field generate a secondary cold plasma with a selective chemical composition, featuring low-energy ions and metastable species in the discharge periphery, ideal for low-damage plasma treatment of material substrates. In this work, we studied the plasma generated by an e-beam using a 4 kV voltage in a pure argon gas environment under a magnetic field of 150 G and in the pressure range of 25–90 mTorr. We measured the absolute spatial density profile of argon (1 s 5 ) metastables in an electron beam generated plasma by laser-induced fluorescence and found it to be of the order of 10 16 m −3 . The electron temperature and the electron density measured by a Langmuir probe were of the order of 10 16 m −3 and less than an eV respectively. Electron-impact quenching was identified as a significant loss mechanism for the Ar(1s 5 ) state, leading to the saturation of the metastable density at higher pressures. Outside the primary ionization region, the spatial distribution of argon metastables followed a linear diffusion profile, indicating negligible additional production in those regions.

EEDF↗
DOE OSTI2023

Dark Field Proton Radiography (Los Alamos LDRD Report)

Proton radiography at LANL (pRad) takes advantage of the high flux, high proton energy, and flexible timing structure of LANSCE's 800-MeV proton beam to visualize dense, dynamic systems under extreme strain. With recent needs identified for the development of the W93, the pRad team has been actively searching for ways to increase the facility's radiographic capabilities. Recently, we discovered that utilizing two magnetic lenses (one upstream and one downstream of the object), and specific combinations of collimator choices within the two lenses, that the signal-to-noise ratio was improved by more than a factor of two. However, we did not fully understand the mechanism, and found that results did not match forward modeling predictions. This work specifically sought to implement an engineering solution to solve this problem, through the development of an actuated collimator that reduced the time required to swap collimator settings from a typical 4-8 hours, to less than 30 seconds. This allowed, in combination with a unique multi-material, multi-resolution stepwedge, to acquire an unprecedently enormous dataset, with in total 24 collimation combinations, and multiple diffuser settings and pulse patterns, equivalent to almost a full run cycle's worth of data in a few days. This massive dataset is currently being harnessed to make a fully-predictive, ultra-fast Python-based model of the system, to dial in specific settings for any future experiments that fully optimize that radiographic sensitivity to the material of interest.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Simulation of plasma transport in MPS-LD linear plasma device by using BOUT++

Here a linear plasma device (LPD) module based on the 2D transport module under the BOUT++ framework is developed in this paper to simulate plasma transport in the LPD. The LPD module includes three parts, i.e. magnetic field calculation, simulation mesh generation and plasma transport setup. The magnetic field is calculated with the circular current loop using the location and current information of each coil. The mesh generation code can produce a simulation mesh for LPD by employing the magnetic field. The plasma transport model is based on the reduced Braginskii equations, which consist of the continuity equation, momentum equation, and energy equation. The fluid neutral model is applied for neutral particle simulation. Deuterium (D) and helium (He) ions and atoms can be simulated by the model. The first attempt to simulate plasma transport in the new LPD, called multiple plasma simulation linear plasma device (MPS-LD), is presented using the developed model. The D plasma transport in the MPS-LD is simulated and benchmarked against the two-point model, showing the validation of the BOUT++ simulation. The effects of radial transport on the heat load to the target are studied, which illustrates the significant impact of the radial diffusivities $D_⟂$ and $\mathcal{X}_{⟂i,e}$ on the plasma. Moreover, the He impurity injection process during the discharge is studied with emphasis on plasma–impurity interactions. The simulation results show that He injection can reduce the plasma energy load to the target significantly, and the efficiency depends on the He source density and injection velocity. The present work provides an alternative and flexible simulation tool for plasma transport in LPDs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

NMR spectroscopy of coin cell batteries with metal casings

Battery cells with metal casings are commonly considered incompatible with nuclear magnetic resonance (NMR) spectroscopy because the oscillating radio-frequency magnetic fields (“rf fields”) responsible for excitation and detection of NMR active nuclei do not penetrate metals. Here, we show that rf fields can still efficiently penetrate nonmetallic layers of coin cells with metal casings provided “B 1 damming” configurations are avoided. With this understanding, we demonstrate noninvasive high-field in situ 7 Li and 19 F NMR of coin cells with metal casings using a traditional external NMR coil. This includes the first NMR measurements of an unmodified commercial off-the-shelf rechargeable battery in operando, from which we detect, resolve, and separate 7 Li NMR signals from elemental Li, anodic β-LiAl, and cathodic Li x MnO 2 compounds. Real-time changes of β-LiAl lithium diffusion rates and variable β-LiAl 7 Li NMR Knight shifts are observed and tied to electrochemically driven changes of the β-LiAl defect structure.

25 ENERGY STORAGE↗
DOE OSTI2024

Thermodynamic and Kinetic Modulation of Methylammonium Lead Bromide Crystallization Revealed by In Situ Monitoring

Hybrid organic–inorganic perovskite (HOIP) crystals are promising optoelectronic materials, but little is known about either the thermodynamic and kinetic controls on crystal growth or the underlying growth mechanism(s). Herein, we use fluid cell atomic force microscopy (AFM) and solution nuclear magnetic resonance (NMR) spectroscopy to investigate growth of the model HOIP crystal CH 3 NH 3 PbBr 3 (MAPbBr 3 ) and to determine how formic acid (HCOOH) modulates the thermodynamics and kinetics of growth. The results show that growth of MAPbBr 3 in dimethylformamide (DMF) proceeds through the classical pathway by the spreading of atomic crystal steps generated at screw dislocations on the {100} surface. Temperature dependent step velocity measurements demonstrate that with increasing concentration, HCOOH decreases both the solubility of MAPbBr 3 and the kinetic coefficient (b) of step movement. 1 H-NMR measurements indicate that HCOOH increases the lifetime of the methylammonium (MA + ) ions and promotes the association of MAPbBr 3 , thus tuning the solubility of the perovskite. HCOOH also alters the molecular tumbling motion and bulk diffusion of the MA + ions, possibly via H-bonding. Further, our findings establish a direct correlation between the mesoscale crystal growth kinetics and the molecular-scale interactions between organic additives and constituent ions, providing unprecedented insights for developing predictive syntheses of HOIP crystals with defined size, crystal habit and shape, and defect distribution.

36 MATERIALS SCIENCE↗