Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “softening”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Further Development of the Tamped Richtmyer-Meshkov Instability Method and Application to Molybdenum Dynamic Strength Calibration and Tabulation

The high pressure and high strain rate dynamic strength of Mo is experimentally and computationally investigated in the 3–20 GPa stress, 50–600 C temperature, and 10 5 –10 6 /s strain rate regimes using a modified tamped Richtmyer-Meshkov instability (RMI) method. Modifications to the tamped RMI method include a method to determine loading states during strain, a new strength calibration function based on interface shape, and a robust uncertainty quantification method. These modifications improve fidelity of the tamped RMI method, allowing evaluation of the compensating effects of pressure hardening, strain rate hardening, strain hardening, and thermal softening. The new calibration function based on interface shape is not limited to sinusoidal corrugations and could be applied to additional interface shapes. Plate impact experiments are performed at Argonne National Laboratory’s Advanced Photon Source’s Dynamic Compression Sector operated by Washington State University (DCS), driving a planar shock front through a corrugated Mo-D 2 O or Mo-C 8 F 18 interface, forcing the corrugation to significantly deform. The extent of interfacial deformation, RMI growth, is experimentally observed using X-ray phase contrast imaging at the DCS. RMI jet lengths and jet shapes are extracted from the experimental radiographs, then used to calibrate numerical simulations performed with the Sandia National Laboratories (SNL) hydrocode CTH. Mo yield strength, Y, as a function of shock pressure, P, strain rate, $\dot{\varepsilon }$, accumulated strain, ϵ, relative volumetric compression, RD , and temperature, T , is determined for each impact experiment and presented. The calibrated Mo yield strength values range 1.2–1.8 GPa, with strength generally decreasing as the impact stress increases. This trend is likely caused by thermal softening or strain localization. The tabular yield strength versus loading condition data presented in this paper can be used to fit complex strength models.

Voorhees, T. J. [Sandia National Lab. (SNL-CA), Li↗
DOE OSTI2021

Hall–Petch strengthening limit through partially active segregation in nanocrystalline Ag-Cu alloys

The breakdown from grain-size strengthening to softening mechanisms is generally well understood for high-purity nanocrystalline materials when the mean grain size decreases to the nanometer range. In nanocrystalline alloys, however, the stabilization of nanosized grains by grain-boundary solute segregation complicates the above mechanisms. Moreover, current segregation models have little predictive power for determining the optimal solute content that maximizes Hall-Petch strengthening effects. In this article, using large-scale hybrid Monte-Carlo/molecular dynamic simulations, we present a systematic study of the Hall–Petch breakdown in Cu-segregated Ag alloys with grain sizes ranging from 8 nm to 59 nm, where three concentration-dependent regimes of plasticity are described: (1) Classical segregation strengthening behavior at low solute contents, (2) shear band-induced softening at high solute contents, and (3) a previously unknown, but extended plateau of maximum strengths for intermediate solute contents from 4 to 15 at.%, which we term as nanocrystalline Sterling silver. We find that flow strengths in nanocrystalline Sterling alloys naturally exhibit a zero-slope limit at the smallest grain sizes that is well below the ideal Hall–Petch strengthening trend. This phenomenon results from partially active grain-boundary segregation that acts to influence interfacial plasticity in some, but not all, grain boundary regions. Our findings amplify the atomic nature of solute segregation and interaction at grain boundaries and its complex roles on grain boundary-mediated plasticity mechanisms in nanocrystalline alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Extreme hardness at high temperature with a lightweight additively manufactured multi-principal element superalloy

Materials are needed that can tolerate increasingly harsh environments, especially ones that retain high strength at extreme temperatures. Higher melting temperature alloys, like those consisting primarily of refractory elements, can greatly increase the efficiency of turbomachinery used in grid electricity production worldwide. Existing alloys, including Ni- and Co-based superalloys, used in components like turbine blades, bearings, and seals, remain a performance limiting factor due to their propensity, despite extensive optimization efforts, for softening and diffusion-driven elongation at temperatures often well above half their melting point. To address this critical materials challenge, we present results from integrating additive manufacturing and alloy design to guide significant improvements in performance via traditionally difficult-to-manufacture refractory alloys. We present an example of a multi-principal element alloy (MPEA), consisting of five refractory elements and aluminum, that exhibited high hardness and specific strength surpassing other known alloys, including superalloys. The alloy shows negligible softening up to 800°C and consists of four compositionally distinct phases, in distinction to previous work on MPEAs. Density functional theory calculations reveal a thermodynamic explanation for the observed temperature-independent hardness and favorability for the formation of this multiplicity of phases.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A B-spline based gradient-enhanced micropolar implicit material point method for large localized inelastic deformations

The quasi-brittle response of cohesive-frictional materials in numerical simulations is commonly represented by softening plasticity or continuum damage models, either individually or in combination. However, classical models, particularly when coupled with non-associated plasticity, often suffer from ill-posedness and a lack of objectivity in numerical simulations. Moreover, the performance of the finite element method significantly degrades in simulations involving finite strains when mesh distortion reaches excessive levels. This represents a challenge for modeling cohesive-frictional materials, given their tendency to experience strongly localized deformations, such as those occurring during shear band dominated failure. Hence, accurate modeling of the response of cohesive-frictional solids is a demanding task. To address these challenges, we present an extension of the material point method (MPM) for the unified gradient-enhanced micropolar continuum, aiming at the analysis of finite localized inelastic deformations in cohesive-frictional materials. The generalized gradient-enhanced micropolar continuum formulation is employed to tackle challenges related to localization and softening material behavior, while the MPM addresses issues arising from excessive deformations. The method utilizes a B-spline formulation for the rigid background mesh to mitigate the well-known cell crossing errors of the MPM. To demonstrate the performance of the method, 2D and 3D numerical studies on localized failure in sandstone in plane strain compression and triaxial extension tests are presented. A comparison with finite element results confirms the suitability of the formulation. Moreover, an efficient numerical implementation of the formulation is presented, and it is demonstrated that the additional MPM specific overhead is negligible.

B-spline↗
DOE OSTI2020

Effect of grain structure and strain rate on dynamic recrystallization and deformation behavior: A phase field-crystal plasticity model

Herein, a phase field (PF)-crystal plasticity (CP) model has been developed to study the effects of initial grain structures, recrystallized grain orientation, and strain rates on dynamic recrystallization and deformation behavior under stress states typical of the Shear Assisted Processing and Extrusion (ShAPE) process. The model takes the geometrically necessary dislocation density as the nucleation criterion for recrystallization and minimizes the deformation energy as the driving force of recrystallized grain growth. Using this model, the effects of grain orientations and strain rates on dynamic recrystallization and deformation behavior were simulated. The results demonstrate that (1) a polycrystalline structure with the texture observed in ShAPE process has the lowest yield stress under the ShAPE stress conditions, (2) the recrystallized grains with the texture observed in ShAPE process largely soften the materials, and (3) higher shear strain rate or rotation speed results in larger magnitude of material softening.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Co-treating flue gas desulfurized effluent and produced water enables novel waste management and recovery of critical minerals

Herein this study reports a novel approach of resource recovery from co-managing two geographically co-located and chemically complementary wastewaters using a pilot-scale treatment process. Designed to treat flue gas desulfurized (FGD) effluent from combustion powerplants and produced water (PW) from energy industries, the process consists of soda-ash softening, nanofiltration (NF), and reverse osmosis (RO). Recovered products are barite, calcite, and low-salinity water. Using field-collected waters, the results show that softening at pH 8.5 produces calcite (yield: 30 kg/m 3 treated water), a chemical used as SO ₂(g) scrubbers. NF treatment under an applied pressure of 3.5 MPa yields a permeate stream laden with monovalent ions (water recovery 60%) and a concentrate stream with a sulfate concentration 1.8 times of the feedwater concentration. Mixing the NF concentrate and PW at a volumetric ratio of 1.0 precipitates a high-density barite material (4.1 g/cm 3 , yield: ~7.5 kg/m 3 mixture) – a critical mineral commonly used as a weighting agent in drilling. The RO treatment recovers >64% water as the permeate, which can be readily used as cooling make-up water at the powerplants. The RO concentrate stream can be further processed in a thermal evaporative system for additional water recovery and brine production.

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

Anomalous elastic behavior of tantalum at high pressures: Experimental and theoretical studies

Elastic wave velocities of polycrystalline tantalum have been measured up to 13.6 GPa at room temperature using ultrasonic interferometry technique in a multi-anvil apparatus. The bulk and shear moduli, as well as their pressure derivatives are obtained using Eulerian 3rd order finite strain equations, yielding K S0 = 193.9(29) GPa, G 0 = 69.3(10) GPa, K S0' = 3.18(5) and G 0' = 0.88(1). Here the Steinberg-Guinan yield stress model is examined using the current shear modulus and its pressure derivative, where a slightly lower A 0 = 0.0127 GPa –1 compared to previous studies is obtained. First-principles calculations on the elasticity of tantalum have also been performed up to 200 GPa. Anomalous softening of the elastic shear constant C 44 and shear wave velocity V S is observed between 80– 180 GPa, whereas the anisotropy shows a softening-stiffening behavior with increasing pressure and reaches a minimum at around 150 GPa. Such anomalies might commonly originate from the nesting of Fermi surface for Group VB metals.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Microstructural evolution and intermetallic formation in Zn-3Mg (wt%) powder mixture processed by high-pressure torsion

Severe plastic deformation (SPD) techniques have been used extensively over the past 40 years for producing strong metals and alloys. High-pressure torsion (HPT) is one of the most promising SPD techniques for achieving high strength through nanoscale grain refinement and phase transformation. In this research, a mixture of pure zinc (Zn) and magnesium (Mg) powders, Zn-3Mg (wt%), was HPT-processed under a pressure of 6 GPa for 1, 5, 10, 20, and 30 turns at room temperature to achieve a high strength biodegradable material. In order to understand the effects of pre-consolidation on the resulting microstructure and hardness, HPT processing was performed on loose powders placed in the die and also on a pre-compacted powder mixture and the resulting HPT disks were characterized by X-ray diffraction, scanning electron microscopy, atom probe tomography, and Vickers microhardness. In both cases, the HPT disk microstructures contained nanoscale grains, and stable and metastable strain-induced intermetallics, but an unusual softening appeared at large shear strains. Grain size, grain morphology, and the formation of different intermetallics were analyzed to explain the unusual hardness distribution, and it was found that an inverse Hall-Petch relationship between hardness and grain size exists. It is suggested that thermally-activated phenomena such as grain boundary sliding contributed to the strain-induced softening of this nano-structured biomaterial due to its low melting point. Finally, the current results are compared with those for HPT-processed cast alloys and hybrids of the same composition.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Peak microhardness recovery via T6 heat treatment with extremely soft annealed Al7075 feedstock in additive friction stir deposition

With standard T6 heat treatments, precipitate-hardened alloys such as Al7075 and Al6061 fabricated using Additive Friction Stir Deposition (AFSD) method restore to noticeablely different peak mechanical properties. Previous research observed similar results but did not provide a viable solution. This study proposes a new processing strategy, namely enhanced plastic deformation by subjecting the Al7075 feedstock to annealed softening treatments. Specifically, by annealing the standard T651-Al7075 feedstock and deposition substrate (179 HV), their Microhardness were first degraded to its lowest state (80 HV), which is called a Soft-to-Soft system. For the Soft-to-Soft system subjected to AFSD, after undergoing appropriate T6 heat treatment, the microhardness of the AFSD 7075 alloy was restored to its peak value (183.1 ± 1.8 HV). The results indicate that this softening strategy can effectively reduce the non-strengthening phase in the AFSD matrix, thereby increasing the strengthening phase (η′ phase). This addresses the issue of traditional AFSD-prepared aluminum alloys being unable to achieve peak performance through T6 heat treatment.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Effects of irradiation spectrum on the microstructural and mechanical properties of bulk metallic glasses

Zr 52.5 Cu 17.9 Ni 14.6 Al 10 Ti 5 (BAM-11) and Cu60Zr 20 Hf 10 Ti 10 bulk metallic glass (BMG) specimens were ion irradiated to characterize and compare their irradiation induced microstructural and mechanical property evolution. For the ion irradiations, samples of each BMG were exposed to 9 MeV Ni 3+ and 5.5 MeV C + ions to a midrange (~1.5 μm depth) dose of 0.5 displacements per atom (dpa) at temperatures ranging from room temperature to 360 °C in order to examine potential effects associated with low (5.5 MeV C) and high (9 MeV Ni) average primary knock on energies. BAM-11 BMG samples were also irradiated by neutrons (E > 0.1 MeV) at ~70 °C to a fluence of 1.4 × 10 20 n/cm 2 (dose of 0.1 dpa). Importantly, no pronounced irradiation spectrum effects were observed for the nanoindentation mechanical properties of the two BMGs following ion irradiation. However, significant softening was observed in the BAM-11 BMG sample irradiated by neutrons, while annealing of the unirradiated samples led to a marked increase in hardening. Finally, results of the nanoindentation experiments indicate that softening is caused by irradiation induced creation of free volume defects, while the hardening is caused by their annihilation.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Strengthening mechanisms, hardening/softening behavior, and microstructure evolution in an LPSO magnesium alloy at elevated temperatures

Here, rapid degradation of the mechanical properties of conventional Mg alloys with temperature precludes their wide application in industry. The LPSO (long period stacking ordered) Mg alloys, on the other hand, show considerable potential to achieve excellent high-temperature mechanical performance. In this study, a Mg-5.1Y-2.4Zn (wt.%) extrusion alloy with the microstructure comprising of α-Mg matrix and LPSO phase was prepared. It is shown that the yield strength of the alloy at 200 °C retained ~83% of the yield strength of 198 MPa at room temperature (RT). Importantly, the ductility also improved, owing to the increased non-basal slip activities at elevated temperatures. In contrast to RT and 200 °C, the alloy exhibited softening tensile behavior with an apparent loss of yield strength at 300 °C. In-situ lattice strain analysis revealed that load transfer from α-Mg to LPSO remains effective at temperatures up to 200 °C, governing the high yield strength of the alloy. The post-mortem analysis by means of electron microscopy confirmed that the LPSO phase is structurally stable at 200 °C. In particular, the high strain hardenability is rationalized by the sandwiched LPSO structures together with stacking faults serving as effective obstacles to the motion of non-basal dislocations. The load transfer effect disappears at the temperature of 300 °C, resulting in the fast deterioration of strength of the alloy. The cracking of the LPSO phase along with the recrystallizing of α-Mg grains were found to be the main causes for the softening behavior of the alloy at 300 °C.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Size effect in stainless steel thin wires under tension

Miniaturization of modern devices drives the exploration into how materials behave in small dimensions. However, the mechanical behaviors of micro-scale specimens under varied loading conditions are still challenging to characterize and predict. In this work, specimen size-dependent weakening effect was experimentally revealed in 316L stainless steel thin wires under tension with wire diameter decreasing from 200 μm down to 61 μm. The underlying mechanisms were primarily associated with surface grain softening and strain localization, which were largely dependent on the number of grains across wire diameter. Through an analog composite model with surface and interior grains that follows the Hall-Petch relation, the extent of surface grain softening was found to coincide with the flow stress reduction due to lack of boundary strengthening. Further using elasto-plastic self-consistent polycrystal modeling, the stress-strain behaviors of interior and surface grains were well described in analogy to those of bulk materials with similar-size grains and extremely large grains, respectively. Subsequently, the analog composite model well simulated the specimen size-dependent weakening effect under tension. This article provided more insights into the size effect exploration and proposed a practical model for describing size-dependent mechanical behaviors of small-scale materials to be applied to miniaturized devices.

316L stainless steel↗
DOE OSTI2025

Hydrogen influences thermal activation parameters for dislocation glide during low cycle fatigue of 316L stainless steel

Measurements of activation areas are used to investigate the effect of hydrogen on the kinetics of dislocation glide during cyclic deformation in cold-worked 316L stainless steel. Non-charged and hydrogen-precharged (H-precharged) specimens were tested in low cycle fatigue (LCF) under plastic strain control. A series of plastic strain rate changes was performed periodically at the peak true plastic strain from the first cycle to half-life, and at various plastic strain values around stable hysteresis loops near half-life to determine the operational activation area, Δ a ∗. Both material conditions demonstrate a rapid increase in Δ a ∗ during the initial rapid softening followed by a region of approximately constant values coinciding with a reduced rate of softening. Near half-life, hydrogen reduces Δ a ∗ at a given true stress due to its effect on the activation distance and obstacle spacing. The magnitudes of Δ a ∗ reveal that bypassing solutes, cutting forest dislocations, and initiating cross slip are important mechanisms of thermally activated dislocation glide at all amplitudes, except hydrogen suppresses cross slip at the lowest plastic strain amplitudes. These results are supported by electron microscopy characterization of deformed microstructures. A Haasen plot analysis indicates that forest dislocations control the kinetics of deformation in both material conditions. It also reveals the presence of athermal obstacles in both non-charged and H-precharged conditions, likely dense dislocation tangles and cell walls. Additionally, the effect of hydrogen on microstructure evolution (by reducing the propensity for cross slip) leads to a dependence of athermal stress on plastic strain amplitude.

Activation area↗
DOE OSTI2023

Micro-beam bending of FCC bicrystals: A comparison between defect dynamics simulations and experiments

To understand the role of the grain boundary (GB) in plasticity at small scale, a concurrently coupled mesoscale plasticity model was developed to simulate micro-bending of bicrystalline micron-sized beams. By coupling dislocation dynamics (DD) with a finite element model (FEM), a novel defect dynamics model provides the means to investigate intricate interactions between dislocations and GBs under various loading conditions. Our simulations of micro-bending agree well with corresponding micro-bending experiments, and they show that mechanical response of bicrystals could have not only hardening but also softening depending on the characters of the GB. In addition, changing the location of the GB in the microbeams results in different mechanical responses; GBs located at the neutral plane show softening compared to single crystals, while inclined GBs located halfway along the length of the beam show little effect. In conclusion, simulation results could provide a clear picture on detailed dislocation-GB interactions, and quantitative resolved shear stress analysis supplemented by dislocation density distribution is used to analyze the mechanical response of bicrystalline samples.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Multi-cycle reload analysis of a long cycle gas-cooled fast modular reactor

There is currently significant interest in deploying HALEU-fueled fast reactors, including the General Atomics (GA) Fast Modular Reactor (FMR). Such reactors can achieve very long fuel cycles, but with multi-batch loading will take decades to reach equilibrium. This motivates design and analysis of both the initial core and multi-cycle reload, which is typically performed using fast-running, deterministic fast reactor codes such as the Argonne Reactor Computation (ARC) codes. In this paper, multicycle reload of the GA FMR is analyzed using the ARC codes. The GA FMR utilizes 19.75 % enriched fuel in a 16 year cycle with a three-batch strategy, with twice-burned fuel placed on the core periphery. The GA FMR has a softened neutron spectrum due to reflecting elements in the core, so the neutronic solution is first benchmarked against the OpenMC Monte Carlo code. Discrepancy on k eff is 400–600 pcm, likely due to the softened neutron spectrum, heterogeneous fuel assembly design and central reflector. However, the rms discrepancy on the assembly power distribution is only 0.6 %, despite the presence of the central reflector. A reload strategy is devised for the first three cycles of such a reactor, ultimately spanning the first 45–48 years of its operation. The fresh core uses 19.75 %, 19.25 % and 16.75 % enriched fuel in place of fresh, once-burned and twice-burned and is then subsequently refueled with only 19.75 % enriched fuel. The cycle length is varied over 3 cycles of operation to balance fuel utilization and reactor availability, specifically with use of an extended 18-year Cycle 1, followed by a shortened 11-year Cycle 2. Cycle 3 is close to the target 16-year length. Finally, placing twice burned assemblies next to the GA FMR central reflector can reduce power peaking by 3 %, at the expense of slightly reducing the cycle length.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Shear localization in metallic materials at high strain rates

Three factors govern adiabatic shear localization: strain hardening (or softening), strain-rate hardening, and thermal softening. It is typically associated with large shear strains (>1), high strain rates (10 3 –10 7 s –1 ), and high temperatures (40–100% of melting point), all of which happen within narrow regions with widths of about 1–200 μm. It is often an undesirable phenomenon, leading to failure, but there are situations where it is desirable, e. g., the generation of machining chips. Here, we review the development of both theoretical and experimental achievements, from the initiation of shear bands to their propagation with emphasis on three aspects: novel experimental techniques, novel materials, and nano/microstructural effects. The principal characteristics of adiabatic shear bands in metallic materials at the nano- and micro-scale are described. Bands that were formerly identified as transformed actually consist of nanocrystalline/ultrafine grains. These grains result from the breakup of the microstructure by a rotational recrystallization process. The evolution of the microstructure inside shear bands and their interactions for hcp, bcc, and fcc alloys, high-entropy alloys, nanocrystalline alloys, and metallic glasses are analyzed mechanistically. In conclusion, the gaps in the field and opportunities for future research are identified. Modern experimental characterization and computational techniques enable a more profound and predictive understanding of adiabatic shear localization and its avoidance in advanced materials.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Understanding the origin of softness in structurally tailored and engineered macromolecular (STEM) gels: A DPD study

Structurally tailored and engineered macromolecular (STEM) gels are polymer networks containing a primary network (the “STEM-0” gel) and side chains that are anchored to this network (forming the “STEM-1” gel). The ability to vary the features of the side chains provides a robust means of tailoring the macroscopic properties of the material. We use dissipative particle dynamics (DPD) simulations to determine the mechanical response of STEM-1 gels to uniaxial compression for various values of the side chain length, n sc . While the Young’s modulus of the material is significantly decreased with increasing n sc (at a fixed density of grafted side chains), above a certain saturation value, increasing n sc does not lead to any further softening of the sample. Using the simulations, we calculate the relevant stresses, which are related to the Young’s modulus and the number of entanglements in the network. We show that the backbone chains become more spread out upon addition of sidechains, leading to a decrease in the physical entanglements between backbone segments. This observation accounts for the decrease of stress in the backbone network and the softening of the networks through the addition of sidechains. For long chains, however, the number density of physical crosslinks between side chains is independent of side chain length; this observation explains the observed stress saturation in the STEM-1 gels. Our approach allows us to correlate the molecular architecture of the gels to the resultant macroscopic mechanical behavior and provide guidelines for fabricating STEM gels with well-defined mechanical properties, which allow the materials to be used for a variety of applications.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A phase-field model of quasi-brittle fracture for pressurized cracks: Application to UO 2 high-burnup microstructure fragmentation

Here, in this paper, we present a phase-field model of quasi-brittle fracture with pressurized cracks, with dedicated applications for polycrystalline materials. The model is formulated as a minimization problem within the variational framework. The external work done by pressure on the crack surfaces is included in the objective function. Several careful modeling choices lead to a regularization-length-independent critical strength. The model is constructed to give a softening response with an underlying linear traction-separation law. The pressure-dependent softening response and the regularization of the prescribed pressure are demonstrated with a (quasi) one-dimensional numerical analysis. In a two-dimensional numerical analysis under plane strain assumptions, the critical stress corresponding to crack propagation (as predicted by our quasi-brittle fracture model) is compared with linear elastic fracture mechanics (LEFM) analytical solutions. Our model is further utilized to simulate fission-gas-induced fragmentation of the UO 2 high-burnup structure (HBS). Simulation results show that pressurized bubbles can cause crack nucleation and propagation, and that the bubble size and the surrounding external pressure affect the critical pressure corresponding to crack nucleation. Simulations of a partial HBS at different recrystallization stages show that different grain structures (due to recrystallization) also influence crack paths and fragmentation morphology.

36 MATERIALS SCIENCE↗