Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Slip transmission”

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

Atomistic and phase field simulations of three dimensional interactions of {$10\bar12$} twins with grain boundaries in Mg: twin transmission and dislocation emission

In polycrystals, the interaction of dislocations and twins with grain boundaries (GBs) plays a role in hardening and formability during plastic deformation. While dislocation-GB interactions are relatively well-understood, twin-GB interactions remain mostly unknown. In this work, an approach using molecular dynamics and phase-field simulations is followed to study the forward and lateral interactions between {$10\bar12$} twins and tilt grain boundaries in Mg. Molecular dynamics results show that the resolved shear stress on slip/twinning modes of the neighboring grain, not the geometric alignment, is the dominant factor in determining the outcome of the twin-GB interactions. For some lateral interaction configurations, as the misorientation angle increases, the resolved shear stress on the same {$10\bar12$} twin variant of the neighboring grain reduces while it increases for slip or I 2 stacking fault emissions or other twin modes such as {$11\bar12$} and {$10\bar11$}, explaining why twin transmission is not seen at high misorientation angles. Furthermore, lateral and forward interactions of the twin with tilt grain boundaries whose misorientation axes are normal to the coherent twin boundary show significantly different outcomes. For the forward interaction, the twin is absorbed and stacking faults are emitted when interacting for low misorientation angles (up to 30°) while the lateral interaction results in twin transmission, nucleation of a {$11\bar12$} twin, and emission of I 2 stacking faults. Finally, comparisons between twin interactions with symmetric and asymmetric tilt GBs with different GB structures show similar outcomes.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Reversing inverse Hall-Petch and direct computation of Hall-Petch coefficients

A 2D bicrystal atomistic model of dislocation transmission through a Σ11<101>{131} symmetric-tilt grain boundary reveals details of Hall-Petch breakdown in the single dislocation regime in Al, Ni, and Cu. Partly based on a previous study, this research determines the stress required for a single dislocation to be transmitted through the Σ11 boundary and finds that single dislocation transmission typically deviates from Hall-Petch behavior because the leading partial dislocation becomes trapped in the boundary and emits glissile grain boundary disconnections that cause boundary sliding deformation at stresses well below boundary transmission stresses. Thus, mechanisms of inverse Hall-Petch, namely grain boundary shear and sliding, are shown to operate in lieu of grain boundary transmission in the single dislocation regime. However, by controlling the applied stresses, inverse Hall-Petch can be reversed and typical Hall- Petch grain boundary transmission regained. This, in turn, allows the direct computation of Hall-Petch coefficients. Further, a second dislocation on the same slip plane preserves typical Hall-Petch behavior for Al and Ni and does not lead to boundary sliding events, thus implying a critical grain size for Hall-Petch breakdown based on the pileup model.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Cross-slip and easy-glide CRSS in titanium: Theoretical predictions and in-situ TEM measurements

This study investigates the mechanics of prismatic and first-order pyramidal $\langle$a$\rangle$ slip in titanium (Ti), elucidating the physics of easy-glide and cross-slip through a combination of theory and experiments. Screw-character prismatic (Pr) dislocations in Ti are of particular interest because their complex cores can be stable or unstable, leading to activation by either cross-slip or planar glide. To investigate these mechanisms, site-specific micro tensile samples were prepared using focused ion beam (FIB) milling and mounted on a push-to-pull (PTP) device for in-situ transmission electron microscopy (TEM) tensile testing. The in-situ experiments provide direct observations of the onset of dislocation motion and the precise determination of the critical resolved shear stress (CRSS) for the activated mechanisms, and its evolution with load cycling. A comprehensive theory has been developed to predict the CRSS values for easy glide, cross-slip, and multiplication of dislocations. Predicted critical stresses for pyramidal (π)-to-Pr and reverse cross-slip agree closely with the experimental measurements. The latter cross-slip stress is a factor of two higher than that of unobstructed planar slip. The model accounts for overlapping dislocation cores and employs a Wigner-Seitz based cell to evaluate misfit energies. By combining ab initio density functional theory (DFT) with anisotropic elasticity, the framework identifies minimum energy pathways for dislocation glide, which can be intermittent and zig-zag. A simplified expression utilizing (π) and Pr Schmid factor ratios is proposed for critical stress corresponding to (π)-to-Pr cross-slip transition. The results are strongly dependent on crystal orientation, underscoring non-Schmid behavior. Overall, this study explores key critical stress parameters essential for informing higher-scale simulations of plasticity in Ti.

Critical stress↗
DOE OSTI2021

Guiding light via slippery liquid-infused porous surfaces

Slippery liquid-infused porous surfaces (SLIPSs) have been explored for many applications, taking advantage of their highly non-wetting property. In this work, we explore the SLIPS as a cladding material for waveguiding. SLIPSs are prepared by infusing perfluoropolyether oil to hydrophobized nanoporous surfaces of silicon. Power loss and transmission efficiency of an HeNe laser (1.82 mW and 632.8 nm) with varying incident angles were measured through microchannels consisting of the SLIPSs as cladding layers (noil = 1.30) and water (nwater = 1.33) as a core, compared to other cladding types including a planar silicon surface and the nanoporous surfaces in hydrophilic (Wenzel state) and hydrophobic (Cassie–Baxter state) conditions with no oil infused. Agreeing with Snell's law, a total internal reflection occurs at the incident angle as high as 14° for the SLIPSs. The waveguide loss at 14° is only 1.8 dB/cm for the SLIPSs, while those for planar silicon, hydrophilic nanoporous, and hydrophobic nanoporous surfaces are 5.9, 7.4, and 4.9 dB/cm, respectively. The power transmission efficiency of the SLIPSs is independent of the porosity because the surfaces are fully covered with the oil layer, whereas those of hydrophilic and hydrophobic nanoporous surfaces, whose pores are filled with water and air, respectively, depend on the porosity. The significantly lower power loss and the insensitivity to the surface porosity are advantages of the SLIPSs over the other surfaces and can benefit in waveguiding applications such as optofluidics.

Asawa, Kaustubh (ORCID:0000000197142555)↗
DOE OSTI2020

The combined effects of hydrogen and aging condition on the deformation and fracture behavior of a precipitation-hardened nickel-base superalloy

The effect of hydrogen (H) on the deformation behavior of Monel K-500 in various isothermal heat treatment conditions (non-aged, under-aged, peak-aged, and over-aged) was assessed via uniaxial mechanical testing. H-charged and non-charged specimens were strained to failure to facilitate a comparison of ductility, fracture surface morphology, strength, and work hardening behavior. For all examined heat treatment conditions, H charging leads to a significant reduction in ductility, which is accompanied by a consistent change in fracture surface morphology from ductile microvoid coalescence to brittle intergranular fracture. While H charging led to a systematic enhancement in the yield strength of all heat treatments, the three age-hardened conditions exhibited a more than 2-fold increase relative to the non-aged heat treatment. This suggests that H modifies the dislocation–precipitate interactions, which also manifest themselves through changes in work hardening metrics related to the dislocation storage and recovery rates. In particular, the H-charged peak-aged specimen exhibited a significant increase in initial hardening (dislocation storage) rate relative to the H-charged under-aged specimen. Transmission electron microscopy of these samples revealed the onset of widespread dislocation looping in the H-charged peak-aged sample, in addition to the planar slip bands characteristic of the non-charged condition. This result suggests that hydrogen induces the particle shearing-to-looping transition at smaller particle sizes. Possible mechanistic explanations for this observed behavior are presented.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Direct imaging of short-range order and its impact on deformation in Ti-6Al

Chemical short-range order (SRO) within a nominally single-phase solid solution is known to affect the mechanical properties of alloys. While SRO has been indirectly related to deformation, direct observation of the SRO domain structure, and its effects on deformation mechanisms at the nanoscale, has remained elusive. Here, we report the direct observation of SRO in relation to deformation using energy-filtered imaging in a transmission electron microscope (TEM). The diffraction contrast is enhanced by reducing the inelastically scattered electrons, revealing subnanometer SRO-enhanced domains. The destruction of these domains by dislocation planar slip is observed after ex situ and in situ TEM mechanical testing. These results confirm the impact of SRO in Ti-Al alloys on the scale of angstroms. The direct confirmation of SRO in relationship to dislocation plasticity in metals can provide insight into how the mechanical behavior of concentrated solid solutions by the material’s thermal history.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Low-temperature failure mechanism of [001] niobium micropillars under uniaxial tension

The plasticity of body-centered cubic (bcc) metals is dependent of temperature as well as sample dimension at the micrometer scale, but the effects of cryogenic temperature on the plasticity and the related failure process in micron-sized bcc metals have not been studied under uniaxial tension. In this work, we utilized in situ cryogenic micro-tensile tests, transmission electron microscopy, and dislocation dynamic simulations to examine the plasticity and failure processes of [001]-oriented bcc niobium micropillars. Our study reveals that a strong suppression of cross-slip at low temperatures prevents dislocation multiplication and leads to a dislocation starvation state, at which no mobile dislocation exists due to the rapid annihilation of dislocations at free surfaces. New dislocations are then nucleated until stress concentration at a slip step creates a micro-crack, the propagation of which leads to catastrophic failure. As a result, this unique failure process results from the combined effects of sample dimension and temperature.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Dislocation-density evolution and pileups in bicrystalline systems

Here, a dislocation-density crystalline plasticity (DCP) framework based on total and partial dislocation densities interactions was used to investigate the behavior of Cu/Pb bicrystals with a focus on GB effects. The modeling predictions were validated with bicrystal compression micropillar experiments. A key new aspect of the modeling approach is to account for partial dislocation-densities. A GB formulation that is directly linked to GB energies was used to monitor GB transmission and blockages, such that pileups can be monitored and predicted at the GB interfaces for misorientations. The predictions indicate that pileups can form due to fully and partially blocked slip-rates and perfect and partial dislocation-densities. As the nominal strain increases from five to fifteen percent, dislocation-densities and pileups significantly increase by almost an order of magnitude. The proposed validated approach provides a microstructural scale predictive framework that accounts for a myriad of defects related to the interactions of partial and perfect dislocation densities that interact at highly misoriented GBs; it is these interactions that are critical to the formation and evolution of dislocation-density pileups that can lead to physically limiting stress accumulations in bicrystals.

36 MATERIALS SCIENCE↗
DOE OSTI2023

On the Low Risk of SSR in Type III Wind Turbines Operating in Grid-Forming Control: Preprint

We have shown in a previous work that the risk of sub-synchronous resonance (SSR) between a wind power plant with Type III wind turbines and series-compensated transmission lines is low when the wind turbines in the plant are operated in grid-forming mode, instead of the standard grid-following mode. The fundamental mechanism behind the improved damping characteristics is explained in this paper by modeling the positive and negative sequence impedances of Type III wind turbines for GFM operation mode. It is discovered that the GFM control naturally acts against the negative resistance of Type III wind turbines that results from an interaction between the proportional gain of the rotor-side converter current controller and negative slip. The developed sequence impedance models and improved damping behavior are verified using PSCAD simulations of a 2.5-MW Type III GFM wind turbine. The modeling predictions are also supported by experimentally measuring the sequence impedance response of a 2.5-MW Type III wind turbine during operation in GFL and GFM modes.

grid-following turbine↗
DOE OSTI2021

An in situ ambient and cryogenic transmission electron microscopy study of the effects of temperature on dislocation behavior in CrCoNi-based high-entropy alloys with low stacking-fault energy

Temperature is known to affect deformation mechanisms in metallic alloys. As temperature decreases, the stacking-fault energy in many face-centered cubic (fcc) alloys decreases, resulting in a change of deformation mode from dislocation slip to deformation twinning. Such an impact of temperature can be more complex in compositionally heterogeneous microstructures that exhibit, for example, local concentration fluctuation such as that in multi-principal element alloys. In this work, we compare the dislocation behavior and mechanical properties of a fcc Cr 20 Mn 10 Fe 30 Co 30 Ni 10 high-entropy alloy at ambient and liquid-nitrogen temperatures. We find that a network of stacking faults is formed by uniformly extended dislocations at ambient temperatures with low stacking-fault energy, whereas at lower temperatures, uneven dissociation of dislocations becomes significant, which results in severe dislocation pile-ups together with their pronounced entanglement. Our findings indicate that as the stacking-fault energy decreases with decreasing temperature, the heterogeneity of the distribution of elements becomes more dominant in tuning the local variation of lattice resistance. As a result, the change in dislocation behavior at low temperatures strongly affects microstructural evolution and consequently leads to significantly more pronounced work hardening.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Mechanical properties and deformation mechanisms of single crystal Mg micropillars subjected to high-strain-rate C-axis compression

Here, the mechanical properties and deformation mechanisms of single crystal magnesium under c-axis quasi-static and high-strain rate compressions are investigated through in situ scanning electron microscope (SEM) experiments and post-mortem transmission electron microscope (TEM) characterization. The findings revealed that ductility and high rates of hardening are preserved for pillars as large as 15 μm. Furthermore, rate effects result in a mild increase in flow stress with plastic deformations controlled primarily by the slip of type dislocations. Importantly and in contrast to other literature reports, plastic deformation occurs in the absence of twining. As the strain increases and plastic deformation exceeds about 4%, crystal rotation activates basal slip, <$\mathrm{a}$> type dislocations, resulting in a more rate independent flow stress. TEM observation on micropillars compressed at a strain rate of 250/s, revealed the activation of {${11}$$\bar{2}$$\bar{2}$} < $\bar{1}$$\bar{1}23$ > slip systems and high mobility of screw dislocations as major contributors to plastic strains in excess of 10% without fracture. These findings are relevant to the design of lightweight materials used in transportation systems, e.g., selection of material grain size. Moreover, the experimental data here reported provides the materials science community with a unique opportunity to validate discrete dislocation dynamics (DDD) formulations employed in multiscale design of materials.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Solute effects upon dislocation motion and recovery in Mg alloys (Final Report)

The objective the research was to develop a firmer understanding of the interactions between substitutional solute atoms and dislocations within Mg alloys. These interactions govern the absolute and relative mobilities of various types of dislocations in Mg (e.g., basal < a>, non-basal < a>, and pyramidal < c+a>). Furthermore, they have an impact on dislocation recovery processes (e.g., cross-glide, climb, rearrangement, and annihilation). Ultimately, solute-dislocation interactions strongly impact a) strain hardening, b) strain rate sensitivity, c) plastic anisotropy of textured polycrystals, d) texture evolution, e) dislocation substructures evolution, and even f) recrystallization behavior. It is broadly known that Mg alloys exhibit poor low temperature formability. Since formability is largely governed by properties (a) – (c) in this list, it is critical to develop a better understanding of solute-dislocation interactions, if one hopes to improve the situation. While the theory of static solute strengthening is well developed, especially for alloys with face centered cubic (FCC) crystal structures, there are outstanding questions related to applications to hexagonal close packed (HCP) crystal structures and dynamic strain aging (DSA) in materials of various structures. DSA has far-reaching implications for metal formability and, in the case of Mg alloys, appears to correlate with the so-called rare-earth (RE) texture which has been shown to benefit formability. The Portevin-Le Chatelier (PLC) effect, and associated negative strain rate sensitivity, occur at higher temperatures (>100℃) in Mg alloys as compared with similar Al alloys, even though they have similar melting points and solute diffusivities. Our preliminary research has shown that modern, physics-based models of DSA can be tuned to describe the behavior of Mg alloys if a rather higher activation enthalpy is assumed for cross-core diffusivity. While this partially explains the delay in DSA to higher temperatures, it is also hypothesized that this delay is due in part to the intrinsically more thermally activated (rate sensitive) nature of non-basal < a> dislocation motion which is required for macroscopic flow of Mg alloys, whereas octahedral slip in many FCC metals like aluminum is essentially athermal at room temperature. It was originally proposed to employ a combination of in-situ diffraction-based experimental characterization to validate existing theory. It was envisioned to perform in-situ transmission electron microscopy (TEM) to assess individual dislocation behavior and in-situ high-energy X-ray diffraction (HEXRD) techniques which were showing great promise for elucidating collective dislocation behavior, 2 including recovery, especially if the contributions to various forms of diffraction peak broadening (𝜂𝜂,𝜔𝜔,and 2𝜃𝜃) can be effectively integrated. Finally, it was envisioned to perform discrete dislocation dynamics (DDD) modeling approaches to aide in the interpretation of both TEM and HEXRD experiments. In the end, mechanical tests were performed on more complex Mg alloys which exhibited evidence of dynamic strain aging, and this led to the establishment of another project. Mechanical test data obtained at McMaster University served as the basis of an assessment of the applicability of Bazinski’s “stress equivalence” theory of solute strengthening to polycrystalline alloys of Mg. Although we did not succeed in applying the approach to Mg alloys, we did develop expertise with the HEXRD approach using a BCC, β-Ti alloy and demonstrated numerous new capabilities that may be applied to any polycrystalline material in collaboration with researchers at CHESS and around the world: (1) assessment of details of the elastoplastic transition (yielding) using a combination of HEXRD and full-field polycrystal plasticity modeling, (2) the first-ever experimental observation of strong stress rotation within the individual grains of a polycrystalline material, and (3) a comprehensive analysis of the grain-level dislocation density evolution based upon diffraction peak broadening along 𝜂𝜂,𝜔𝜔,and 2𝜃𝜃 directions. This final aspect allowed us to confirm that dislocations were gliding on multiple plane types and not restricted to {110} type planes, and it also provided clues as to why some grains were unloading during straining, with surprising implications for our understanding of the effects of geometrically necessary dislocations (GNDs). Finally, graduate student, Mohammed Shabana, developed a MATLAB code which confirmed the conclusions of Prof. Catalin Picu (Rensallear Polytechnic Institute, RPI) regarding the effect of solute-trapped, forest dislocations on the breaking stress of Lomer lock junctions in FCC metal alloys. He applied the same anisotropic line-tension model to a variety of dislocation junction configurations and found an inconsistency in the widely cited results of Dupuy and Fivel regarding the Hirth Lock, and he outlined an approach to extend these finding to HCP Mg alloys that we are still pursuing with discretionary fundings at UVA.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Dynamic fracture processes in hydrogen embrittled iron

Mechanistic understanding of hydrogen embrittlement is critical to the transport and storage of hydrogen as a sustainable fuel. Experimental observation of dynamic fracture processes in hydrogen environments has been challenging. Here, in-situ transmission X-ray microscopy is used to image polycrystalline iron foils that are concurrently strained and electrochemically charged with hydrogen. Using this technique, voids and microcracks down to ∼30 nm in size are identified at or near the crack tip and tracked as strain is increased. We observe both hydrogen-induced transgranular failure and hydrogen-induced intergranular failure with distinct fracture processes. We find that transgranular fracture occurs through the formation of sharp branched structures at the main crack tip that are reminiscent of slip band formation. This is followed by void growth and coalescence, and interaction with secondary cracks that lead to propagation of the main crack. Intergranular fracture is found to occur through void coalescence along grain boundaries, along with significant secondary cracking near the crack tip. The crack growth rate is greater for void-mediated intergranular failure than for slip-mediated transgranular failure. Additionally, the removal of hydrogen results in recovered ductility for a crack initially propagated under hydrogen charging.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Damage Detection and Electrical Performance Impact of Flat-Flexible Circuits

Flat-flexible circuits or ‘flex circuits’ have varied applications in electronics and in explosive detonation systems. They pose advantages over conventional coaxial cables in their ability to fit into small mechanical envelopes and routings. A typical flex circuit is composed of a laminated assembly of multiple layers of conductor and dielectric materials. Copper, because of its high electrical conductivity, is one of the primary materials used in flex circuit conductors. The conductor must maintain structural integrity and electrical continuity when the flexible circuit is bent during processing or to fit an application envelope. Under some strain conditions a cable may have degraded electrical performance even though visible damage may not be evident. Dislocation density can increase by mechanisms such as dislocation forest hardening, where gliding dislocations in primary slip planes get pinned at obstacles and form Frank-Read sources. Dislocation density increases within the conductor layers may impact electrical resistivity but would not necessarily result in visible mechanical failure. In some detonation applications, such as for building implosion for demolition, the time for signal transmission and voltage used for firing detonators can be critical. If the resistance in a single flex cable within a building demolition array is too high, the resultant demolition shape may be asymmetric, leading to unplanned damage to nearby structures and potential hazards to public safety. This study will apply non-destructive evaluation techniques to evaluate the detectability of damage accumulation in flat, flexible copper-conductor cables and discuss potential implications to their practical use.

42 ENGINEERING↗
DOE OSTI2023

On the Low Risk of SSR in Type III Wind Turbines Operating with Grid-Forming Control

We have shown in previous work that the risk of subsynchronous resonance between wind power plants with Type III wind turbines and series-compensated transmission lines is low when the wind turbines are operated in grid-forming (GFM) mode instead of the standard grid-following mode. This paper explains the fundamental mechanism behind the improved damping characteristics by modeling the positive- and negative-sequence impedances of Type III wind turbines in GFM control mode. Here it is discovered that the GFM control naturally acts against the negative resistance behavior of Type III wind turbines at subsynchronous frequencies that results from an interaction between the rotor-side converter current controller and negative slip. The developed impedance models in the sequence domain and the improved damping behavior are verified using PSCAD simulations of a 2.5-MW Type III GFM wind turbine. The modeling predictions are also supported by experimentally measuring the impedance response of a 2.5-MW Type III wind turbine in GFM mode.

17 WIND ENERGY↗
DOE OSTI2023

On the Low Risk of Subsynchronous Resonance (SSR) in Type III Wind Turbines Operating in Grid-Forming Control

We have shown in a previous work that the risk of sub-synchronous resonance (SSR) between a wind power plant with Type III wind turbines and series-compensated transmission lines is low when the wind turbines in the plant are operated in grid-forming mode, instead of the standard grid-following mode. The fundamental mechanism behind the improved damping characteristics is explained in this paper by modeling the positive and negative sequence impedances of Type III wind turbines for GFM operation mode. It is discovered that the GFM control naturally acts against the negative resistance of Type III wind turbines that results from an interaction between the proportional gain of the rotor-side converter current controller and negative slip. The developed sequence impedance models and improved damping behavior are verified using PSCAD simulations of a 2.5-MW Type III GFM wind turbine. The modeling predictions are also supported by experimentally measuring the sequence impedance response of a 2.5-MW Type III wind turbine during operation in GFL and GFM modes.

grid-following turbine↗
DOE OSTI2022

Kink mechanism in Cu/Nb nanolaminates explored by $\mathcal{in}$ $\mathcal{situ}$ pillar compression

We report Nano metallic laminates (NMLs) exhibit different failure modes depending on the loading conditions due to their mechanical anisotropies. Kinking is a typical failure mode in many NMLs compressed along a layer-parallel direction. However, a detailed description of the microstructure evolution during kink band (KB) formation and an in-depth understanding of the formation mechanisms are lacking. In this work, the KB process is investigated in Cu/Nb NMLs by in situ micro pillar compression in the scanning electron microscope (SEM) along a layer-parallel direction. Post-mortem S/TEM and transmission Kikuchi diffraction (TKD) analyses show that kink banding leads to significant microstructure changes characterized by an accumulation of geometrically necessary dislocations (GNDs) and of tilt geometrically necessary boundaries (GNBs) near KB boundaries (KBBs). The distinct microstructure evolution implies that KB formation is facilitated by the inhomogeneous microstructures resulting in constrained deformation modes. Specifically, dislocations active on slip planes nearly parallel to the interfaces make a major contribution to kink evolution after the onset of kinking. Once layer-parallel slip systems are activated, preexisting lattice dislocations and dislocations nucleating from interfaces will accumulate as GNDs near KBBs via the stochastic storage of lattice dislocations that have certain Burgers vectors. GNDs can further transform into GNBs via cross-slip and climb driven processes near the KBB. Furthermore, GNBs near KBBs can grow by incorporating more GNDs or by coalescence to accommodate the KB evolution. We further hypothesize that microstructural perturbations and their ensuing stresses can initiate KB formation in Cu/Nb NMLs.

36 MATERIALS SCIENCE↗