Search NASA⌕ Search

SEARCH · Search NASA

Results for “elastic strains”

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 289 records · Page 16

Dynamic flow stress of pure polycrystalline aluminum: Pressure-shear plate impact experiments and extension of dislocation-based modeling to large strains

The dynamic thermo-mechanical behavior of pure aluminum has attracted renewed interest lately due to experimental observations of an anomalous increase in Hugoniot Elastic Limit (HEL) at incipient plasticity and elevated temperatures in polycrystalline pure metals. In context of current dislocation-mediated plasticity models for metals, this increase in dynamic strength is indicative of a transition in the rate controlling mechanism for dislocation glide from being thermally assisted to phonon-drag restricted due to increase in phonon viscosity at elevated temperatures. Though these studies have helped to shed light on these important mechanisms operative in FCC metals at incipient plasticity, the extent to which they contribute to flow stress, particularly at larger strains, remains unclear. Here, we address these questions through a combined experimental and modeling effort focused on investigating the evolution of dynamic flow stress in polycrystalline aluminum using experimental data gathered from a series of combined pressure-and-shear plate impact (PSPI) experiments designed to reveal the flow stress of pure aluminum at strain rates ~ 10 5 /s, plastic strains of up to 40% and temperatures ranging from room to 866 K. In all cases, the flow stress of aluminum, as inferred from the measured transverse particle velocity histories at the free surface of an elastic tungsten carbide target plate, reveals saturation with increasing plastic strains at stress levels that decrease with increasing test temperatures. Numerical simulations are performed to correlate the experimentally observed temperature and strain rate dependence of flow stress at small and large plastic strains using the Austin-McDowell dislocation-mediated plasticity model parametrized to normal plate impact experiments conducted in an earlier study by Zaretsky and Kanel (2012). Extensions to the model are made to better represent dynamic behavior of pure aluminum at larger plastic strains as observed in elevated temperature split Hopkinson Pressure bar (SPHB) experiments of Samanta (1971) and Lindholm and Yeakly (1965), and the combined pressure-and-shear plate-impact experiments conducted in the present study. The main theoretical extension to the Austin-McDowell model made in this paper is the introduction of a new rate- and temperature- dependent dynamic recovery function, which can potentially allow for an effective reduction in the rate of accumulation of dislocations at large plastic strains. The numerical predictions of the revised and re-calibrated plasticity model are brought into agreement with the experimental observations, correlating sufficiently well with dynamic yield stress at incipient plasticity and the flow stress levels at larger plastic strains, plastic strain rates in the range 10 3 – 10 6 /s, and elevated temperatures up to near melt. The model, in concert with the experimental measurements, suggests that at incipient plasticity the rate governing mechanism for plastic flow is phonon-drag restricted dislocation glide, whereas, at higher magnitudes of plastic strain, it transitions to stress-assisted thermally activated glide. The transition between these two rate governing mechanisms is controlled by the evolution of dislocations throughout the deformation process.

36 MATERIALS SCIENCE↗

Study of the interplay between lower-order and higher-order energetic strain-gradient effects in polycrystal plasticity

in this report strain-gradient (SG) plasticity refers to a class of non-local theories in which gradients of plastic slip determine the storage of geometrically necessary dislocations, introducing a length-scale dependence in the mechanical behavior of crystalline materials, which is otherwise lacking in local theories. In this work, we incorporate lower-order (LO) and higher-order energetic (HOE) strain-gradient effects into a crystal plasticity fast Fourier transform (FFT)-based formulation to investigate the interplay of the length scale that each strain-gradient term introduces at the microscale, and the mechanical properties that result at the macroscale. For an applicable range of length scales, we consider two systems: a 1-D two-phase face centered cubic (FCC) laminate and a 3-D FCC polycrystal, and two uniaxial deformation modes: monotonic tension and cyclic tension–compression. We show that increases in the individual LO and HOE length scales increase the hardening rate and strength of the material, respectively. When combined, the strong LO hardening is less pronounced than the effect alone due to the lowering of the gradients due to the HOE microstress. We demonstrate that the LO and HOE hardening manifest as “isotropic” (yield surface expansion) and “kinematic” (yield surface shift) effects, respectively, consistent with their theoretical origins. We show that in cyclic loading, the Bauschinger effect emerges in both local and non-local calculations and link its origins and severity to the behavior in the strain field, slip-system rates, and the HOE microforce.

36 MATERIALS SCIENCE↗

Mapping the structural–mechanical landscape of amorphous carbon with ReaxFF molecular dynamics

We use ReaxFF molecular dynamics (MD) to investigate the relationship between structural and mechanical properties in bulk and nanostructured amorphous carbon (a-C). The liquid-quench MD method is used to generate isotropic bulk samples with mass densities ranging from 0.96 to 3.29 g/cm3. Structural analysis identifies two types of structures with distinct short- and medium-range order: lower-density sp2-dominated a-C, which is characterized by a bimodal ring-size distribution, and higher-density sp3-dominated tetrahedral amorphous carbon (ta-C), exhibiting a unimodal ring-size distribution. Stress–strain MD simulations and analysis reveal how an atomistic structure impacts elastic properties and post-yield atomic rearrangements. All stretched structures demonstrate elastic isotropy and plasticity driven by a ring-size expansion mechanism reflected in changes in ring statistics. The plastic region is substantially larger in ta-C than in a-C due to the post-yield shift from sp3 to sp2 C dominant bonding. In both a-C and ta-C, ultimate failure occurs when a reactive crack, traversed by long sp chains, forms and propagates predominantly perpendicular to the direction of the applied strain. Oxygen infiltration into the fractured region significantly reduces stress resistance, primarily through the early rupture of long sp chains. MD simulations and analysis are extended to a-C slabs, a-C nanotubes, and partially a-C nanotubes. The latter nanostructure highlights the differences between the elastically isotropic a-C walls, which develop circumferential cracking, and the crystalline walls, which tear along crystallographic directions. These results provide a strong foundation for further computational characterization of a-C materials.

Dernov, A. (ORCID:0009000004220973)↗

Examining the role of glycoside hydrolases in local rheology of Pseudomonas aeruginosa biofilms

Current research strategies in the treatment of biofilm infections have focused on dispersal, in which bacteria are made to vacate the extracellular polymeric substance (EPS) surrounding them and return to a planktonic state where antimicrobial treatments are more effective. Glycoside hydrolases (GHs), which cleave bonds in EPS polysaccharides, have been shown to promote dispersal in Pseudomonas aeruginosa biofilms. The dispersal mechanism is possibly due to GHs’ ability to directly release bacteria from the EPS, disrupt EPS’ ability to regulate the environment, or reduce overall mechanical stability. In this work, passive microrheology is used to examine the relevance of the last mechanism by exploring the effects of three GHs (α-amylase, cellulase, and xylanase) known to disperse P. aeruginosa on local biofilm viscoelasticity. Compared to control studies in wild-type strains, it is found that treatment with all three GHs results in statistically relatively less elastic and stiffer biofilms, indicating that changes to mechanical stability may be a factor in effective dispersal. Both cellulase and xylanase were observed to have the greatest impact in creating a less stiff and elastic biofilm; these GHs have been observed to be effective at dispersal in the published results. Each GH was further tested on biofilms grown with strains that produced EPS missing specific polysaccharide components. Cellulase specifically targeted Psl, which forms the major structural and mechanical backbone of the EPS, explaining its efficacy in dispersal. However, xylanase did not appear to exhibit any affinity to any polysaccharide within the EPS based on the microrheology results. Overall, these results suggest that the local microrheology of the biofilms is impacted by GHs and that may be one of the factors that is causing the ability of these therapeutics to enhance dispersal.

59 BASIC BIOLOGICAL SCIENCES↗

Measurement of the Static Nonlinear Third-Order Elastic Moduli of Rocks: Problems and Applicability

The third-order elastic (TOE) model has been used to describe the widely observed nonlinear mechanical behaviors of earth materials. In addition to linear elastic constants ( λ , μ ), three nonlinear elastic moduli ( A , B , C ) are required for isotropic rocks. Contrary to previous research on dynamic TOE moduli, this study followed the protocol to measure strain and stress under uniaxial and hydrostatic compressive tests statically, which were later used to invert for the full set of TOE moduli for four standard rock types with differing pore structures of a porous oolitic limestone, quartz-rich sandstones, and a dense crystalline basalt. The applicability of the TOE model to characterize nonlinearity depends on the fulfillment of path-independence and small-strain assumptions. Using the measured static TOE moduli, the finite element model demonstrates that the stress in the vicinity of the wellbore is more amplified than the stress in the linear elastic case, which leads to a wider zone of rock failure around the wellbore. Due to the long-standing discrepancy between static and dynamic moduli, the rarely reported full set of static TOE moduli is necessary and will benefit future research in understanding the effect of rock nonlinearity on geophysical and geomechanical applications, such as long-term safe storage of CO 2 and generating process of geohazards.

58 GEOSCIENCES↗

Distributed fiber optic strain sensing of bending deformation of a well mockup in the laboratory

Well integrity is critical to the safety and success of subsurface energy exploration and management, as leakage of fluids from subsurface reservoirs is often induced by compromised wells. This study investigates bending deformation of a mockup of an oil/gas well that is subjected to loads expected in deviated wells under reservoir compaction and assesses the feasibility of utilizing distributed fiber optic strain sensing to monitor the deformation. Here, a well mockup, which consists of outer and inner steel pipes with a cemented annulus, is tested under three-point bending loading and its strain and curvature development is monitored by Brillouin optical time domain reflectometry/analysis (BOTDR/A) as well as optical frequency domain reflectometry (OFDR). The primary objective of this research is to assess the strain sensing performance of newly fabricated fiber optic cables and to identify key cable characteristics which could improve the quality of distributed strain measurements with BOTDR/A. Results show that the tight-buffered cable is best suited for strain sensing as its maximum errors in the strain measurement were -36% and -24% against conventional sensors at the maximum elastic and plastic bending loads, respectively, whereas those of the non-tight-buffered cable were -45% and -71%, respectively. Similar trends were obtained in the bending curvature measurement. The detailed design of the tight-buffered cable is presented to elucidate key characteristics of such a cable, which will facilitate accurate distributed strain sensing in oil and gas wells.

42 ENGINEERING↗

Magnetoelastic properties of multiferroic hexagonal ErMnO 3

Here, the strength and dynamics of magnetoelastic coupling through the paramagnetic (PM) – antiferromagnetic (AFM) – ferrimagnetic (FIM) transitions in multiferroic hexagonal ErMnO 3 have been investigated by Resonant Ultrasound Spectroscopy. Elastic stiffening by up to 2% below the PM – AFM transition at 80 K arises from biquadratic coupling between strain and the magnetic order parameter with relaxation times longer than ~10 -6 s for the response of spins to changes in strain. In contrast with YMnO 3 , the PM – AFM transition in ErMnO 3 is accompanied by a peak in acoustic loss immediately below the Néel point which is interpreted in terms of strain relaxation accompanying ordering of spins of Er 3+ at 4b sites. Changes in the magnetic ordering scheme at the AFM – FIM transition near 3 K are accompanied by elastic softening of ~0.03 %. During poling of the low temperature ferrimagnetic structure round magnetic hysteresis loops, small changes in elastic stiffness which arise due to the contribution of piezomagnetic and/or piezoelectric moduli are detected. Contributions of piezoelectric moduli to acoustic resonance frequencies also permit changes in the configuration of ferroelectric domains to be detected in response both to cycling through this transition and to application of a magnetic field. A peak in acoustic loss in the vicinity of 250 K is attributed to strain-mediated pinning/freezing of some aspect of the domain microstructure with an activation energy of ~0.25–0.3 eV. A return to the original elastic properties on heating to temperatures above ~250 K is interpreted in terms of backswitching of domains to the configuration they had at the start. These observations confirm the existence of subtle variations in magnetoelastic coupling behaviour relating to both the magnetic order parameters and magnetic domain structures.

36 MATERIALS SCIENCE↗

Multi-stage load partitioning in additively manufactured Al6061+TiC nanocomposite characterized by in-situ neutron diffraction

Metal matrix composites (MMCs), combining metal matrix and ceramic particles, exhibit high mechanical strength and stiffness compared to conventional metallic materials. During fusion-based additive manufacturing (AM), MMCs undergo a rapid heating and cooling thermal history, which affects the interfacial bonding, thermal misfit stress and mechanical load transfer behavior between the metal matrix and particles, thus impacting the bulk mechanical performance. Here, we investigate the deformation dynamics and phase-specific load transfer in fusion-based AMed Al6061 +TiC MMCs through in-situ neutron diffraction. By calculating the phase-specific lattice strains, a multi-stage load transfer and deformation behavior during uniaxial compression is revealed: (1) elastic deformation in Stage I for both phases; (2) sudden stress rebalance between two phases in Stage II, evidenced by a sudden decrease in Al lattice strain and increase in TiC; (3) active load carrying by both phases in Stage III, with both phases experiencing an increase in lattice strains. Further, the deformation mechanism of each stage is deducted by correlating the evolutions of the interphase stresses and peaks’ broadening and intensity of the Al matrix. The local plastic deformation of the Al matrix near the phase interface is triggered and leads to stress rebalancing in Stage II. The global plastic deformation subsequently propagates throughout the Al matrix in Stage III, and the composite is further hardened through both dislocation multiplication and load transfer. The findings offer valuable insights into deformation and load-sharing behavior in AMed MMCs.

36 MATERIALS SCIENCE↗

Tensile testing and fracture mechanism analysis of polyvinyl alcohol nanofibrous webs

Abstract A tensile properties testing study was conducted to understand the influence of thickness, cross‐head speed (speed of testing), gauge length (GL; specimen test length), and sample shape on important tensile properties of polyvinyl alcohol (PVA) nanofiber webs. The effects of each testing parameter on load at break, extension at break, Young's modulus, and tensile stress–strain curve of PVA nanofiber webs are analyzed. The Welch two sample t ‐tests show the significant difference among tested data. Using interaction plots, two‐way analysis of variance, and margin mean plots, the interaction effects among testing parameters have been analyzed. Of all the factors, cross‐head speed, the interaction among GL, and sample thickness (GL: Thickness) and the interaction among GL, testing speed and sample thickness (GL: Speed: Thickness) have significant influence on the tensile properties of PVA nanofiber webs. Moreover, the hypothesized model of mechanism of tensile strain–stress curve of PVA nanofiber webs has been proposed. Based on the model, the tensile strain–stress curve can be split into three stages: linear elastic, partial break up, and complete breakage. This study will provide a better understanding of tensile testing parameters' effects and their interaction effects on the tensile properties of nanowebs.

Lou, Lihua↗

Seismic Visibility of Eclogite in the Earth’s Upper Mantle—Implications From High Pressure‐Temperature Single‐Crystal Elastic Properties of Omphacite

Abstract Identifying and locating the geochemical and geophysical heterogeneities in the Earth’s interior is one of the most important and challenging tasks for the deep Earth scientists. Subducted oceanic crust metamorphizes into the dense eclogite in the upper mantle and is considered as a major cause of geochemical and geophysical heterogeneities in the deep Earth. In order to detect eclogitic materials inside the Earth, precise measurements of the high pressure‐temperature single‐crystal elasticity of major minerals in eclogite are thus exceedingly important. Omphacite, a Na,Al‐bearing clinopyroxene, constitutes up to 75 vol% of eclogite. In the present study, we performed the first high pressure‐temperature single‐crystal elasticity measurements of omphacite using Brillouin spectroscopy. Utilizing the finite‐strain approach, we obtained the following thermoelastic parameters for omphacite: K S0 ’ = 4.5(1), G 0 ’ = 1.53(5), ∂K S0 /∂ T = −0.029(5) GPa/K, ∂G 0 /∂ T = −0.013(5) GPa/K, with K S0 = 123(3) GPa, G 0 = 74(2) GPa, and ρ 0 = 3.34(1) g/cm 3 . We found that the seismic velocities of undeformed eclogite are similar to pyrolite at the depths of 200–300 and 410–500 km, thus eclogite is seismically invisible at these depths. Combined with the lattice‐preferred orientations of the omphacite in naturally deformed eclogites, we also modeled seismic anisotropy of eclogite at various pressure‐temperature conditions. A 10 km thick subducted eclogitic crust can result in ∼0.2 s shear wave splitting in the Earth’s upper mantle.

Hao, Ming↗

Geodetic Evidence for Distributed Shear Below the Brittle Crust of the Walker Lane, Western United States

Abstract Models of active deformation of the Earth's crust are predominantly represented with dislocations having a downdip continuation into the lower crust, where the fault slips continuously. This model predicts surface strain accumulation concentrated near the fault during the interseismic period. In an alternative model, faults do not extend beneath the elastic portion of the crust and are accompanied by a wide zone of distributed shear underneath, predicting a more constant strain rate lacking concentrations at the faults. We use high‐precision GPS data collected across the northern and central Walker Lane, USA— a region of complex faulting near the western edge of the Basin and Range Province to evaluate which model is appropriate. Despite the existence of dense continuous and semi‐continuous geodetic networks that have been surveyed for ∼20 years, the horizontal velocities reveal no evidence of localized strain accumulation across the fault surface expressions. Instead, deformation within the Walker Lane is uniformly linear, suggesting that the surface deformation reflects distributed shear within the ductile crust rather than focused deformation at faults. This suggests no downdip extension of the faults below the seismogenic layer. The shear zone is 172 ± 6 km wide in the northernmost Walker Lane narrowing to 116 ± 4 km in the central Walker Lane. The total velocity budget across the shear zone is 7.2 ± 0.1 mm/yr in the north, increasing to 10.1 ± 0.1 mm/yr in the central Walker Lane. We conclude that assuming the presence of lower crustal dislocations when estimating geodetic faults slip rates may be inappropriate.

Geochemistry & Geophysics↗

Modeling the contributions to acoustic nonlinearity from complex dislocation networks using 3D dislocation dynamics

Nonlinear ultrasonic parameters are highly sensitive to microstructural features that affect macroscale material behavior, providing a nondestructive means to characterize their evolution. Although dislocations are known to be a strong source of acoustic nonlinearity, establishing quantitative links between the acoustic nonlinearity parameter (β), measured via Second Harmonic Generation, and dislocation morphology—such as dislocation length and density—remains an open challenge. This work advances the numerical modeling of dislocation–β relationships using 3D dislocation dynamics (DD) simulations in two approaches: a “static” method computing strain and stress fields from dislocation configurations in the absence of external loading, and a “quasi-static” method to estimate β from the curvature of dislocation lines under applied load. First, the static method is combined with finite element analysis to investigate a recent assertion that heterogeneous initial strain fields can induce higher harmonic generation in a linear elastic medium; the present results do not corroborate this outcome. Then, the quasi-static method is applied to multiple-dislocation scenarios through parametric studies, revealing behaviors not predicted by analytical models, such as the competing interactions of edge and screw dislocations and the significant influence of applied stress on β. Finally, the simulations are used to model SHG experimental results and validate the hypothesis that β can decrease during plastic deformation, despite increasing dislocation density. As the DD code used here is open-source, it provides a practical platform for future investigation into microstructure–β relationships important to the interpretation of SHG results.

Materials science↗

Magnetoelastic standing waves induced in UO 2 by microsecond magnetic field pulses

Significance UO 2 is of particular interest due to its use as a nuclear fuel, its unusually poor thermal conductivity, and its uniquely strong magnetoelastic properties characterized by a linear coupling between strain and magnetization. We clearly observe mechanical standing waves and elastic resonances in a single crystal, induced by microsecond magnetic field pulses. The phase and amplitude of the mechanical resonances are modulated by the piezomagnetic switching behavior, offering a sensitive probe for otherwise hard to detect magnetic dynamics in UO 2 . We highlight the importance of magnetic domains in the attenuation of the mechanical resonances and the robustness of the antiferromagnetic state in UO 2 against externally applied magnetic fields to 150 T.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Non-invasive Pipe Pressure Monitoring for Safeguards

Our project explores a new solution to complement current safeguards toolbox by independently, and noninvasively, determining the pressure inside uranium hexafluoride (UF 6 )-carrying Gas Centrifuge Enrichment Plant (GCEP) pipes. Proof-of-principle, non-invasive pipe pressure monitoring tests were conducted through measurements of hoop and axial strain using an optical-interference dilatometry technique based on fiber Bragg gratings (FBGs) inscribed in the core of telecommunication-type 125 µm diameter optical fibers. The fibers were affixed to commercial stainless steel and aluminum tubes and pipes of the type expected to be found at GCEP’s, each measuring four inches in diameter. When the internal air pressure in the pipes changes their external diameter and length follow, and the pitch in the FBG shrinks and/or expands accordingly. In turn, these changes result in a shift of the light wavelength reflected by the FBGs. Several FBGs physically attached to the pipe, as well as control FBGs not attached, were interrogated during the process using state-of-the-art commercial interrogators that operate on swept-wavelength laser technology. Strains measured were contrasted against calculations based on the known elastic and mechanical properties of the pipes. Good agreement is found, and strategies are proposed for further improvement of the experimental resolution as well as future implementation of this novel approach for GCEP safeguards.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Interim Development of new Class B Code Case with Variable Design Lifetime and enhanced design rules to guard against cyclic structural failure mode

This report summarizes updates on the ongoing development of the new American Society of Mechanical Engineers Section III, Division 5, Class B rules to address the potential failure modes of Class B components. An initial design-by-analysis approach and supporting design rules are presented. The elastic-perfectly plastic (EPP) analysis approach has been adopted for primary stress limit and ratcheting check to prevent the component failure against the structural failure mode against primary load and strain accumulation due to cyclic loads. The creep-fatigue damage assessment uses explicitly defined elastic follow-up calculated from the stress concentration region of given component. Damage fraction calculation uses a novel coupled approach to capture the influence of the elastic follow-up and provide adequate conservatism for Class B components. The proposed Class B rules does not use the stress classification approach and uses the combined loads to assess component.

36 MATERIALS SCIENCE↗

Mechanical Properties of Nominally Ductile and Brittle Zones within the Caney Shale Formation

We present a case study from the Caney Shale, which is an emerging hydrocarbon-bearing formation in Oklahoma. Open hole logs lead to a nominal distinction of five zones in an exploration well, two of which are designated “ductile” and three of which are designated as less ductile “reservoir” zones. Subsequently, a detailed mechanical characterization was carried out on core plugs taken from each of these zones. Results show nominally “ductile” zones are weaker in terms of UCS, tensile strength, fracture toughness, and confined compressive strength. However, by a variety of brittleness indices, it is found the nominal “ductile” regions do not behave as expected. They are weaker, but not consistently more ductile. Furthermore, some brittleness indices lead to contradictory conclusions. For example, zones which would be considered the most brittle according to the proportion of nonlinear stress-strain behavior prior to failure and/or by indices based on elastic properties would be considered the least brittle according to their internal friction angle.

Benge, Margaret↗

Enhancing fatigue life by ductile-transformable multicomponent B2 precipitates in a high-entropy alloy

Catastrophic accidents caused by fatigue failures often occur in engineering structures. Thus, a fundamental understanding of cyclic-deformation and fatigue-failure mechanisms is critical for the development of fatigue-resistant structural materials. Here we report a high-entropy alloy with enhanced fatigue life by ductile-transformable multicomponent B2 precipitates. Its cyclic-deformation mechanisms are revealed by real-time in-situ neutron diffraction, transmission-electron microscopy, crystal-plasticity modeling, and Monte-Carlo simulations. Multiple cyclic-deformation mechanisms, including dislocation slips, precipitation strengthening, deformation twinning, and reversible martensitic phase transformation, are observed in the studied high-entropy alloy. Its improved fatigue performance at low strain amplitudes, i.e., the high fatigue-crack-initiation resistance, is attributed to the high elasticity, plastic deformability, and martensitic transformation of the B2-strengthening phase. This study shows that fatigue-resistant alloys can be developed by incorporating strengthening ductile-transformable multicomponent intermetallic phases.

36 MATERIALS SCIENCE↗

Dynamic Strength of Iron at High Pressures and Strain Rates

Accurate modeling of meteorite impacts, and deformation of planetary cores require characterization of the flow strength and in-elasticity of iron in its different phases. In this letter, we investigate the flow strength of both the ambient α-phase and high-pressure ϵ-phase of iron at strain-rates of 1 × 10 5 s –1 and pressures up to 42 GPa using high pressure - pressure shear plate impact experiments. We report the strength of the ϵ-iron to be significantly higher than α-phase but consequently one order smaller than the previously reported dynamic strength at high pressures. As a result, the complete stress-strain response of the ϵ-phase is reported for the first time.

36 MATERIALS SCIENCE↗