Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Ni”

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 829 records · Page 46

Joining AA7099 to Ni-Cr-Mo Steel Using Single Pass Friction Stir Dovetailing and AA6061 Butter Layer

A novel single pass friction stir dovetailing (FSD) was developed to join thick plate AA7099 to Ni-Cr-Mo steel in a lap configuration through the use of an AA6061 butter layer. The butter layer approach can be shown to significantly reduce interfacial failures by inhibiting the formation of brittle Fe–Al intermetallic compounds (IMCs) which are typically formed during conventional friction stir welding (FSW) of AA7099 directly to steel. This study evaluates the strength and microstructural evolution of the aluminum/steel joint interface between single pass FSD of AA7099/AA6061(butter)/steel fabricated using three different FSD tool configurations. The butter layer was achieved by placement of a free-floating strip of AA6061 within the dovetail groove that was pre-machined in the steel substrate. Friction stir welding was then used to simultaneously join the thick plate AA7099, the AA6061 strip and the AA7099 in a single weld pass. To optimize the joining process, three different pin geometries (full thread, half smooth/half thread, no thread) were evaluated; it was demonstrated that the half smooth/half thread tool produced the strongest weld joint due to the absence of a Zn rich layer at the Al/steel interface. The efficacy of this unique processing technique was tested using lap shear tensile testing as well as high resolution transmission electron microscopy (TEM). A discussion of the new FSD technique, joint configurations, and process parameters are provided along with joint microstructural analyses and mechanical performance.

Olszta, Matthew J.↗

Development of a heterogeneous nanostructure through abnormal recrystallization of a nanotwinned Ni superalloy

This work explores the development of a heterogeneous nanostructured material through leveraging abnormal recrystallization, which is a prominent phenomenon in coarse-grained Ni-based superalloys. Additionally, through synthesis of a sputtered Inconel 725 film with a heterogeneous distribution of stored energy and subsequent aging treatments at 730°C, a unique combination of grain sizes and morphologies was observed throughout the thickness of the material. Three distinct domains are formed in the aged microstructure, where abnormally large grains are observed in-between a nanocrystalline and a nanotwinned region. In order to investigate the transitions towards a heterogeneous structure, crystallographic orientation and elemental mapping at interval aging times up to 8 h revealed the microstructural evolution and precipitation behavior. From the experimental observations and the detailed analysis of this study, the current methodology can be utilized to further expand the design space of current heterogeneous nanostructured materials.

36 MATERIALS SCIENCE↗

Design of a V–Ti–Ni alloy with superelastic nano-precipitates

Stress-induced martensitic transformations enable metastable alloys to exhibit enhanced strain hardening capacity, leading to improved formability and toughness. As is well-known from transformation-induced plasticity (TRIP) steels, however, the resulting martensite can limit ductility and fatigue life due to its intrinsic brittleness. In this work, we explore an alloy design strategy that utilizes stress-induced martensitic transformations but does not retain the martensite phase. This strategy is based on the introduction of superelastic nano-precipitates, which exhibit reverse transformation after initial stress-induced forward transformation. To this end, utilizing ab-initio simulations and thermodynamic calculations we designed and produced a V 45 Ti 30 Ni 25 (at%) alloy. In this alloy, TiNi is present as nano-precipitates uniformly distributed within a ductile V-rich base-centered cubic (bcc) β matrix, as well as being present as a larger matrix phase. We characterized the microstructure of the produced alloy using various scanning electron microscopy (SEM) and transmission electron microscopy (TEM) methods. Furthermore, the bulk mechanical properties of the alloy are demonstrated through tensile tests, and the reversible transformation in each of the TiNi morphologies were confirmed by in-situ TEM micro-pillar compression experiments, in-situ high-energy diffraction synchrotron cyclic tensile tests, indentation experiments, and differential scanning calorimetry experiments. The observed transformation pathways and variables impacting phase stability are critically discussed.

59 BASIC BIOLOGICAL SCIENCES↗

Plastic deformation of single crystals of the equiatomic Cr-Mn-Fe-Co-Ni high-entropy alloy in tension and compression from 10 K to 1273 K

In this report the plastic deformation behavior of single crystals of the quinary, equiatomic Cr-Mn-Fe-Co-Ni high-entropy alloy (HEA) with the face-centered cubic structure has been investigated in tension and compression as a function of crystal orientation and temperature from 10 K to 1373 K. The critical resolved shear stress (CRSS) for {111}<110> slip at room temperature is 42-45 MPa. It does not depend much on crystal orientation (i.e., the Schmid law holds true) and the sense (tension vs. compression) of the applied load. The CRSS for {111}<110> slip increases with the decrease in temperature, without showing any significant inertia effects at cryogenic temperatures below 77 K. Extrapolation from the measured yield stresses down to 10 K yields a CRSS value at 0 K of 168 MPa. At cryogenic temperatures, the measured strain-rate sensitivity of flow stress is consistent with a very small activation volume. The concept of stress equivalence holds true both for the temperature dependence of CRSS and the stress dependence of activation volume, indicating that solid-solution hardening is the major strengthening mechanism. Deformation twinning occurs at 77 K but not at room temperature, resulting in higher tensile elongation to failure at 77 K than at room temperature. Deformation twinning at 77 K occurs at a shear stress of 378 MPa on conjugate ($\bar1\bar11$) planes in the form of Lüders deformation after large plastic strain (about 85%) achieved by the stage I (easy glide) and stage II (linear work-hardening) deformation.

36 MATERIALS SCIENCE↗

An additively manufactured near-eutectic Al-Ce-Ni-Mn-Zr alloy with high creep resistance

A new additively manufactured (AM) Al-7.5Ce-4.5Ni-0.4Mn-0.7Zr (wt.%) near-eutectic alloy is reported, which shows unprecedented creep resistance up to 400 °C (a homologous temperature of 0.72). The eutectic solidification microstructure comprises ~ 27 vol% of coarsening-resistant second phase network with an ultrafine (<100 nm) inter-phase spacing. Both Mn and Zr contribute to creep resistance of the alloy. Small amount of Mn addition promotes selection of coarsening resistant phases without compromising the alloy processability. Zr not only improves hot-tearing resistance, but further enhances the second phase coarsening resistance resulting in improved creep resistance. Neutron diffraction performed during creep deformation reveals that the underlying mechanism for creep resistance in this alloy is impedance to dislocation motion stemming from the ultrafine eutectic solidification microstructure, whereas load transfer strengthening becomes less effective as the creep temperature increases. The second phase forms a continuous network in the as-fabricated condition, which is maintained during long-term creep at 300 °C. However, this network is fragmented into fine dispersoids at higher temperatures. It is proposed that the rate-limiting deformation mechanism at 300–400 °C is (i) dislocation climb for the alloy with fragmented second phase dispersoids and (ii) Orowan looping for the alloy with a continuous second phase network. In conclusion, the present design of an AM-processable multicomponent eutectic alloy with high creep resistance can be applied to other metallic systems exhibiting eutectic reactions, with expected extreme creep resistance.

36 MATERIALS SCIENCE↗

Influence of microstructural heterogeneities on small-scale mechanical properties of an additively manufactured Al-Ce-Ni-Mn alloy

Laser powder bed fusion-based additive manufacturing (AM) is a promising method to fabricate creep-resistant Al-rare earth alloys (Al-Ce-Ni-Mn) with stable microstructures at up to 400°C. However, creep testing of these alloys at high temperatures shows that void coalescence and failure initiation occurs along the melt pool boundaries in the microstructure. Hence it is crucial to understand how the local mechanical behavior of the melt pool boundaries would influence the global properties of the AM produced alloy. In this study, in situ nanoindentation conducted at room temperature and 300°C revealed a reduced hardness at the melt pool boundaries. Similarly, micro-pillar compression showed a slight decline in yield strength at these boundaries, indicating that they are the weak spots in the microstructure. Such multimodal local mechanical property studies are necessary for understanding the influence of melt pool boundaries on the bulk response of fusion-based AM alloys.

36 MATERIALS SCIENCE↗

Synergy of tensile strength and high cycle fatigue properties in a novel additively manufactured Al-Ni-Ti-Zr alloy with a heterogeneous microstructure

Alloy design strategies in additive manufacturing (AM) to achieve grain refinement and terminal eutectic solidification have been introduced to engineer Al alloys having microstructural hierarchy and heterogeneity. Such alloy design strategies enable crack-free builds with an expanded AM processing window and pushed the strength limit in Al alloys. However, fatigue performance of Al alloys made by AM is restricted by the presence of process induced defects and its stochasticity. In this work, tensile and high cycle fatigue (HCF) behavior of a novel Al-Ni-Ti-Zr alloy with a heterogeneous microstructure is studied in the as-built condition, supplemented by detailed microstructural and mechanical characterization. Excellent strength-ductility synergy of 342 MPa and 16% failure strain achieved in the alloy was associated with the microstructural attributes that pertain to the novel alloy. Additionally, the alloy showed excellent HCF performance with a fatigue endurance limit to ultimate tensile strength ratio of 0.29 in flexural fatigue mode. The study revealed the existence of multiple crack retardation mechanisms and favorable crack propagation pathways through the fine-grained regions which enabled good fatigue performance to the alloy. Further, a probabilistic model has been used to estimate the fatigue life of the alloy as a function of the stochastic microstructure by utilizing the statistical distribution of pores, solid-state inclusions, and grains in the AM Al alloy. Finally, the model parametric trends are consistent with the experimental observations.

36 MATERIALS SCIENCE↗

Creep ductility limiting mechanisms in an additively manufactured Al-Ce-Ni-Mn-Zr alloy

Tensile creep response and cavitation damage evolution in an additively manufactured Al-7.5Ce-4.5Ni-0.4Mn-0.7Zr (wt%) alloy with peak-aging and overaging treatments were investigated in the 300–400 ºC range. Microstructural heterogeneity and its response to heat treatment and subsequent creep deformation were studied to understand the interplay between cavity formation, creep lifetime and ductility. Increasing the applied stress activated the nucleation of more cavities, an experimental observation that is well described using the vacancy accumulation model. Cavities nucleated prematurely due to localized plasticity in the denuded zones that formed at/near melt-pool or grain boundaries. Microstructure/deformation heterogeneity with consequent evolution of stress triaxiality, especially at lower stresses, causes accelerated cavitation, thus producing low creep ductility (∼ 0.2–2.4 %), compared to (∼12–21 %) ductility of the alloy measured by regular tensile tests at equivalent temperatures. A constrained diffusional cavity growth mechanism with continuous cavity nucleation during creep is established as the dominant mechanism, implying that cavitation involves vacancy diffusion, yet its growth rate is dictated by the minimum creep rate. In conclusion, the ductility-limiting creep and cavitation mechanisms discussed here provide new insight into the creep behavior of 3D-printed metallic alloys.

Additive manufacturing↗

Ab initio simulations on the pure Cr lattice stability at 0K: Verification with the Fe-Cr and Ni-Cr binary systems

Significant discrepancies have been observed and discussed on the lattice stability of Cr between the predictions from the ab initio calculations and the CALPHAD approach. In the current work, we carefully examined the possible structures for pure Cr and reviewed the history back from how Kaufman originally determined the Gibbs energy of FCC-Cr in the 1970s. The reliability of Cr lattice stability derived by the CALPHAD and ab initio approaches was systematically discussed. It is concluded that the Cr lattice stability based on the CALPHAD approach has large uncertainty. Meanwhile, we cannot claim that the ab initio H FCC-Cr is error-free as FCC-Cr is an unstable phase under ambient conditions. The present work shows that the ab initio HFCC-Cr can be a viable scientific approach. As both approaches have their limitations, the present work propose to integrate the ab initio results into the CALPHAD platform for the development of the next generation CALPHAD database. As a result, the Fe-Cr and Ni-Cr binary systems were chosen as two case studies demonstrating the capability to adopt the ab initio Cr lattice stability directly into the current CALPHAD database framework.

36 MATERIALS SCIENCE↗

Unlocking naphtha from polyolefins using Ni-based hydrocracking catalysts

Naphtha (C 5 -C 12 alkanes) is a platform feedstock in the petroleum industry and an ideal target for reutilizing plastic waste at scale. Hydrocracking has emerged as a promising technology for deconstructing polyolefins to naphtha. However, contemporary catalysts rely on Pt to achieve high rates, while earth-abundant metals (EAM) perform poorly. Herein, we develop high-performance Ni/BEA catalysts for polyolefin hydrocracking, achieving complete low-density polyethylene (LDPE) deconstruction with 80 % maximum naphtha yield within 12 h at 250 °C, 60 bar H 2 , and a catalyst-to-polymer ratio of 1:100. These catalysts are versatile, accommodating virgin resins and commercial polymer products, and are directly reusable and regenerable. Comparative analysis highlights that 5Ni/BEA(25) exhibits the highest naphtha productivity of 11.4 $\huge{(}$$\frac{g_{Naphtha}}{g_{cat}{•h}}$$\huge{)}$, surpassing previously reported Pt- and EAM-based catalysts by 2.5–6.2x. Technoeconomic and life-cycle analyses reveal naphtha from LDPE hydrocracking is economically and environmentally competitive to fossil-fuel-derived naphtha. In conclusion, these advancements pave the way for the industrial implementation of earth-abundant Ni-based catalysts in polyolefin hydrocracking, followed by retrofitting steam crackers to address plastics waste at scale.

Circular Economy↗

Applicability of the Delplot method for the determination of catalytic reaction sequences: Hydrodeoxygenation of $\gamma$-valerolactone on Ni 2 P/MCM-41

Contact-time measurements based on sites are used to investigate the reaction sequence of the hydrodeoxygenation of γ-valerolactone (GVL) on Ni 2 P/MCM-41. Here the study is carried out at 300 °C and 0.5 MPa with site times of 0.29–6.1 s. A reaction network is proposed in which pentanoic acid and 2-methyltetrahydrofuran (2-MTHF) are formed in parallel initially and pentanal, butane and pentane are formed subsequently. A first-rank Delplot analysis (selectivity) gives a good fit, but higher-rank Delplots do not. The lack of applicability is explained by the magnitude of the rate constants, in which the rate constant for one of the later steps (the consumption of pentanal) is substantially larger than those for the preceding steps. For this reason, pentanal demonstrated the characteristics of a primary product and consequently affected the discrimination of higher rank products. A simulation of a model sequence with assumed rate constants confirms the reason why the Delplot method does not work for this reaction. Since the rate constants cannot be determined a priori to be of close magnitude, the Delplot method cannot be generally used to ascertain reaction sequences.

42 ENGINEERING↗

Quantifying the effects of artificial aging on the ignition and self-propagating reactions of Ni(V)/Al multilayers

Bimetallic, reactive multilayers are uniformly structured materials composed of alternating nanoscale layers that may be ignited to produce self-propagating intermetallic-formation reactions. When reactive multilayers age, there is a change in local composition and loss of stored chemical energy due to mass transport and rearrangement at the interfaces. Here, to quantify the long-term reliability of commercial Ni(V)/Al multilayers, the effects of accelerated aging on both ignition sensitivity and self-propagating reactions have been examined. Thermally aged samples were characterized using transmission electron microscopy, differential scanning calorimetry, and laser ignition combined with high-speed videography. The analytically quantifiable nature of both continuous wave laser ignition and reactive wave propagation enabled calculations of Arrhenius rate constants for as-received and various heat-treated multilayers. With heat treatment, there is a change in the intermixed thickness and interfacial chemistry that decreases the activation energy for point ignition but increases it for self-propagating reactions. This finding implies that with increased thermal aging, the reaction becomes easier to facilitate in the solid state but harder in the liquid phase.

Energetic material↗

Fractal coatings of Ni and NiYSZ for high-temperature corrosion mitigation in solar salt

This paper introduces a novel approach to corrosion mitigation using electrodeposited, fractal Ni and NiYSZ coatings. Highly adherent, durable multiscale coatings with fractal dimensions above 1.80 were fabricated. Corrosion characteristics of coated plain and etched surfaces of In800H, SS310, SS316 and SS347, in comparison to Ha230, were studied in molten 60% NaNO 3 + 40% KNO 3 at 600 °C. The corrosion rate of coated Fe-based alloys reduced by up to 68% compared to uncoated alloys and was on par with or lower than the corrosion rate of Ha230. Furthermore, the corrosion-resistant coatings provide for use of low-cost alloys in high temperature corrosive environments.

36 MATERIALS SCIENCE↗

Simultaneous Single Crystal Growth and Segregation of Ni-Rich Cathode Enabled by Nanoscale Phase Separation for Advanced Lithium-Ion Batteries

Here, a novel nanoscale phase separation process has been discovered to promote the growth and segregation of single-crystal LiNi 0.8 Mn 0.1 Co 0.1 O 2 (NMC811). This process occurs directly during high-temperature calcination without significant agglomeration. The key lies in converting transition metal hydroxide (TM(OH) 2 ) precursors with well-controlled morphology into transition metal oxide (TMO) intermediates before reacting them with lithium salt. The nanoscale redistribution of Ni in TMO, resulting from the concurrent formation of spinel and rock salt phases, helps to deagglomerate the clusters of later-formed NMC811 crystals. The as-prepared single-crystal NMC811 is further validated in a 2Ah pouch cell, demonstrating 1,000 stable cycles. The fundamentally new reaction mechanism of single-crystal growth and segregation provides a new direction for large-scale synthesis of a broad range of single crystals for advanced energy storage.

25 ENERGY STORAGE↗

Transport of coexisting Ni-Cu sulfide liquid and silicate melt in partially molten peridotite

Transport of coexisting sulfide and silicate melts in partially molten peridotite contributes to the redistribution of chalcophile elements within the upper mantle as well as the genesis of magmatic Ni-Cu sulfide deposits, but has not been investigated systematically. Using laboratory experiments, theoretical calculations, and X-ray synchrotron microtomography, this study documents the topology and considers controls on the extraction of two immiscible liquids during partial melting of mantle peridotite. Under hydrostatic conditions, the measured dihedral angle at silicate melt-mineral-mineral junctions is 13.7–21.3°. Silicate melt is distributed along grain edges forming incompletely interconnected melt channels that disconnected by some dead ends at melt fractions ~7–9 vol%. Application of theoretically predicted permeability (k ~10 -14 -10 -16 m 2 ) permits estimation of the extraction velocity of silicate melt of 0.7-11.1 μm/day within a single interconnected melt channel. In the absence of silicate melt, isolated sulfide droplets (3.77 vol%) show a sulfide-olivine-olivine dihedral angle of 91.5–101.3°. However, in the presence of silicate melt, sulfide droplets (average size ~2.53±2.14μm, 1σ) are partially surrounded by silicate melt and stranded in triple junctions or melt pockets due to the limitation of the smallest dimension (0.3μm) of melt channels. Thus, the extraction of sulfide liquid is highly restricted by these dead ends and the smallest dimension of melt channels during porous flow of silicate melt. In contrast, during large-strain shear deformation (strain ~1.6–2.5), initially stranded sulfide droplets were elongated and extracted with silicate melt into liquid-rich sheets with a length of several hundred microns, constantly oriented at 14.3±4.5°to the shear plane and antithetic to the shear direction. The angle is lower than that (18-30°) of those sheets containing sulfide liquid only, indicating that silicate melt dominates liquid-rich sheets. Driven by stress, silicate melt-dominated liquid-rich sheets open the appropriately oriented grain boundaries between silicate minerals, thus providing an efficient pathway for the extraction of sulfide liquid during deformation. When such sheet-like channels remain open, sulfide droplets (> millimeter-scale) can be potentially mobile through high strain domains of the upper mantle, contributing to the addition of chalcophile elements and the fertilization of the lithospheric mantle.

high P-T experiment↗

Activity model for 36 elements in Fe-Ni-Si-S-C liquids with application to terrestrial planet accretion and mantle geochemistry: New data for Ru, Re, Pt, Os, Ti, Nb, and Ta

Understanding siderophile element partitioning between metal and silicate melts under diverse conditions can be used to place important constraints on the materials and conditions of planetary accretion and core formation, as well as post core formation processes. However, the effects of Si on the partitioning and activity coefficients for these elements are not well known, despite Si likely being one of the dominant light elements in Earth’s core. To address this gap in understanding, we have undertaken a systematic study of the highly siderophile elements Re, Pt, Os, and Ru, and the refractory lithophile elements Nb, Ta and Ti at 1600 °C and 1 GPa, to derive epsilon interaction parameters for these elements in FeSi metallic liquids. Positive epsilon interaction parameters were measured for Nb, Ta, Ti, Ru, Re, Pt, and Os, indicating that dissolved Si in Fe liquids causes a decrease in their metal/silicate partition coefficients (or ‘silicophobic’ behavior). Furthermore, ε$_{Re, Os, or Ru}^{Si}$ > ε$_{Re, Os, or Ru}^{S}$ which means Si causes a larger decrease in D(metal/silicate) than S, and the chalcophile behavior expected from some elements will be completely masked by the presence of Si in a metallic liquid. The new parameters are used to update an activity model that now includes 36 siderophile elements in Fe-Ni-Si-S-C liquids (27 trace elements considered here). Systematic assessment of these 27 elements shows which have the strongest affinity for Si, C, and S, and also how activity coefficients for these elements would vary during accretion and core formation in Earth, Mars, and Mercury of widely differing fO 2 and core compositional conditions. The activity model is combined with new partitioning expressions for Mo, W, Cr, Re Ru, Pt, and Os and applied to aspects of post core formation mantle geochemistry of Earth, Mars, and Mercury. Our updated expressions show that the BSE Mo/W ratio can easily be achieved with metal/silicate partitioning during growth of the Earth, whereas Re, Os and Ru become lower than and highly fractionated compared with BSE values during core formation and accretion, and thus nearly 99% of their BSE abundances are likely contributed by late accretion. Ru isotopes should be a very good indicator of the source material for the late accretion. The high Pt/Os and Re/Os developed in a deepening magma ocean during the growth of the Earth, indicates 186 Os and 187 Os isotopes could be coupled if this ancient material remained isolated and subsequently became entrained in mantle plumes and measured in surficial lavas. The extent to which this occurred will be limited by the low Os content of this ancient material, thus requiring mixing as a major component in plume sources. Martian mantle Hf/W ratio stays low during accretion and core formation modelling, suggesting that W isotope anomalies are more likely due to solid/liquid silicate fractionation than to core formation. Finally, Ti contents measured by MESSENGER at Mercury’s surface can be explained by segregation of either a metallic core (IW-6 to -8) or metallic core + sulfide (IW-4 to -7.5) followed by mantle melting.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Interplay between hydrogen, temperature, and character angle on dissociated dislocation energies in Fe–Ni–Cr austenitic stainless steels

Dislocation energy has an important role in the mechanical performance of structural metals. While dislocation energies cannot be fully obtained from continuum theories due to the contribution of the dislocation core, they have been calculated via atomistic simulations in elemental metals. However, constraints on the local atomic environments have prevented the use of such approaches in systems that incorporate alloying or interstitial solutes. In this work, we develop robust molecular dynamics methods to resolve these issues through a geometric construction of dislocation dipoles and the calculation of time-averaged energies. Furthermore, we apply these methods to calculate dislocation energies (including core energies) in an Fe 70 Ni 11 Cr 19 austenitic steel at a variety of character angles, hydrogen concentrations, temperatures, and dipole spacings. The resulting highly converged energies show an excellent agreement with continuum expressions. Overall, hydrogen concentrations up to 1.0 % do not have a significant effect on the elastic parameters and dislocation energy. The methods and insights derived in this work have the potential to facilitate the calculation of dislocation energies in a wide range of systems, and to guide our understanding of hydrogen embrittlement.

Alloyed systems↗

Solid-solid phase transition via the liquid in a Pd 43 Cu 27 Ni 10 P 20 bulk metallic glass under conventional conditions

The microscopic mechanism of solid-solid phase transitions is a long-standing fundamental issue in materials science. Here we directly report the experimental evidence of the existence of an intermediate liquid during the solid-solid phase transition in a bulk Pd 43 Cu 27 Ni 10 P 20 alloy by in situ high temperature high energy X-ray diffraction, combined with the differential scanning calorimetry at slow heating rates down to 5 K/min. We elucidate that this intermediate liquid is attributed to the melting of the low-melting-point crystalline phases, rather than induced by the elastic stress energy along the interface proposed in literatures. This could also be a general feature in multicomponent alloys if a low-melting-point crystalline phase is formed together with other crystalline phases during heating. Furthermore, all results obtained here are scientifically sound and provide new insight into the phase transition theory, which will trigger more studies on the microscopic mechanism of solid-solid phase transition via an intermediate liquid.

36 MATERIALS SCIENCE↗