Search NASA⌕ Search

SEARCH · Search NASA

Results for “High loading rate”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Using Co-Simulation to Model Interconnect-Scale Power Systems from Loads to Generators

Co-simulation is a modeling technique that allows analysts to combine simulation tools and their corresponding models to exchange data during run-time, allowing the creation of larger and more complex models across heterogeneous domains. HELICS is a co-simulation platform developed over the past six years that has been shown to be effective for these multi-domain analysis. Recently, a HELICS-based analysis was completed where the ERCOT electrical interconnect in the United States was modeled in high detail from bulk power system generation to individual customer loads. This model was used to evaluate a flat-rate and transactive energy tariff with integrated wholesale and retail real-time and day-ahead energy markets. This modeling allows detailed analysis showing how the operations of the power system under these tariffs impact all actors in the power system, from individual customers to bulk power system operators.

co-simulation, HELICS, transactive energy system, ↗
DOE OSTI2025

Laser ablation of high-loading Li-ion battery electrodes improves accessible capacity and cycle life for Behind-the-Meter Storage

Adoption of Behind-the-Meter Storage (BTMS) requires design of batteries that enable high safety, long cycle life, and low cost at the system level. Pairing Li 4 Ti 5 O 12 (LTO) with LiMn 2 O 4 (LMO) achieves targets related to safety and cycle life, but these materials' low energy densities contribute to higher cost at the system scale. Increasing electrode loading is a simple approach to improve energy density, but comes with a trade-off in electrode utilization due to long, tortuous Li + diffusion pathways. Here, laser ablation is used to microstructure (pattern) high-loading electrodes to enhance electrode performance through improved Li + diffusion pathways. Four cell types, comprising combinations of standard or patterned anode and cathode, were prepared to evaluate the effects of laser ablation at each electrode. A rate test shows that patterning electrodes enhances active material utilization at ≳1C rates. Patterning the cathode yields the most benefit, as cells with a patterned cathode demonstrate a ~20% higher accessible capacity than those without at 1.4C. Additionally, 1C capacity retention of cells with patterned cathode (91% through 3000 cycles) is significantly improved over cells with only the anode patterned (64%) and non-patterned electrodes (50%). Characterization of post-mortem cells before and after refreshing their electrolyte suggests that 1C capacity retention is improved by mitigation of electrode "dry-out". We hypothesize that the microstructure acts as a reservoir of additional electrolyte, or a path for gas to escape, so that active material remains wetted throughout long-term cycling, and/or the microstructure may reduce localized, gas-forming overpotentials in the high-loading electrode.

25 ENERGY STORAGE↗
DOE OSTI2021

Analysis of Residential Time-of-Use Utility Rate Structures and Economic Implications for Thermal Energy Storage

Thermal energy storage (TES) is an increasingly popular tool to level out the daily electrical demand and add stability to the electrical grid as more intermittent renewable energy sources are installed. TES systems can locally decouple high thermal loads from the operation of a heat pump or reduce the electrical energy demand of the heat pump by providing a more favorable temperature gradient. In addition, many policy makers and utility providers have introduced time-of-use (TOU) rate schedules for residential customers to better reflect the price of electricity generation and demand for specific times. TOU rate schedules price grid-provided electricity differently throughout the day depending on the region’s climate, time of year, and electrical production. Large differences between on-peak and off-peak electrical prices may create an economic advantage for a residential customer to install a TES system. In this work, the economic and energy savings are calculated for a modeled 2,400 square foot residential building with water/ice-based TES using a TOU rate structure. The weather data is from Fresno County, CA, ASHRAE climate zone 3B, and a representative residential TOU utility rate structure from Pacific Gas and Electric (PG&E) was used. The results showed that total energy consumption could be reduced by 14.5% with an 87.5% reduction in on-peak energy usage when the TES is installed. The cost of operating this system for space cooling was reduced by nearly 20% using the sample utility rate plan.

Sultan, Sara↗
DOE OSTI2021

An investigation into the effects of cyclic strain rate on the high cycle and very high cycle fatigue behaviors of wrought and additively manufactured Inconel 718

Additive Manufacturing (AM) has increasingly been used to fabricate parts in aerospace applications, which may require service lives beyond ten-million cycles due to the imposed high loading frequencies. Understanding the very high cycle fatigue (VHCF) behavior of these additive manufactured (AM) parts is an important step towards their design and qualification processes. This study focuses on the high cycle fatigue (HCF) and VHCF behaviors of both wrought and laser beam-powder bed fusion (LB-PBF) fabricated Inconel 718 in machined/polished surface condition, emphasizing on the influence of test frequency (i.e., cyclic strain rate). Uniaxial, fully-reversed force- and stress-controlled fatigue tests were conducted utilizing a servo-hydraulic and an ultrasonic test system operating at 5 Hz and 20 kHz, respectively, on wrought as well as LB-PBF vertically and diagonally built specimens. Fatigue cracks in the majority of the specimens were found to initiate from intra-granular slip bands near or at the surface, which gives rise to strong anisotropy in fatigue resistance in LB-PBF specimens due to the presence of columnar grains along the build directions. Longer fatigue lives were obtained at 20 kHz, which was ascribed to possibly lower-than-intended stresses applied in the ultrasonic tests. The corrected stress-life fatigue data at 20 kHz were found to converge to the one obtained from conventional testing at 5 Hz, implying no effect of cyclic strain rate on the fatigue behavior of Inconel 718 regardless of the fabrication process. The findings of this work confirm the use of ultrasonic fatigue testing to expedite generation of AM materials data to keep up with the current demand; however, the applied stress may need to be corrected.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Nucleation rate controlled grain boundary and lattice creep

Nucleation versus diffusion rate-limited bicrystal and single crystal creep exhibit different scaling dependencies that enable the mechanisms to be isolated when measured as a function of sample size. It has recently been suggested that nucleation rate-limited kinetic models generally describe the non-Newtonian portion of the creep response well, but more direct evidence is required. This work analyzes the grain boundary creep response of UO 2 , a pyrochlore high entropy oxide, silver, and palladium, along with the lattice creep of silver using small-scale in situ loading in the transmission electron microscope. At small sizes, each system exhibits scale dependence associated with nucleation rate-limited kinetics. Fits of the data produce activation volumes on the order of a few Burgers vectors cubed with positive temperature coefficients as expected for nucleation kinetics. The activation enthalpies fall in the range of about 0.4 eV to 1.7 eV, being lower for the metals and higher for the oxides.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Progressing Analysis of Variable Electric Rates (PAVER) Study

The Progressing Analysis of Variable Electric Rates (PAVER) study analyzed the impact of a range of time-varying electric rates on the performance of a regional electric grid and the resulting costs for participating and non-participating customers. This analysis leveraged and extended the work of PNNL’s Distribution System Operator with Transactive (DSO+T) study. Five different rate designs were included: a flat volumetric energy charge, a typical Time of Use (TOU) rate, a dynamic energy (DE) rate (based on wholesale locational marginal prices), a dynamic energy and capacity (DE+C) rate, and, finally, a Block and Swing (B&S) rate that billed customers based on their average load profile at constant pricing, but used the DE+C dynamic price for load deviations from their average profile. These rates were analyzed in a large-scale co-simulation of an entire regional grid with a customer population representative of the current state. A large fraction (80%) of residential and commercial customers were assumed to participate in these time-varying rates with automatically controlled HVAC, water heaters, electric vehicles, and batteries. This study assumed no industrial sector participation. The DE and DE+C rates saw system peak loads reduced by 6-7%, while the large participation in the TOU rate case saw a significant rebound effect and a resulting peak load increase of >5%. The impacts to the annual and peak system demand impacted system wholesale prices and the overall grid operating costs. This cost structure determined the revenue needed to be collected from customers by each rate design. Participating customers on the DE and DE+C rates (located in one of the modeled DSOs) saw reductions in average annual electricity bills of 11-17% with average increases in monthly bill variation of no more than 13%. At such high participation levels, TOU customers saw 10% higher average annual bills (due to system-wide rebound effects) and average increased monthly bill variation of 16%. Residential owners of large flexible loads (such as electric vehicles) saw larger bill savings (17-20%) when on a fully dynamic rate. The presence of on-site generation (such as rooftop solar) did not appear to appreciably change customer outcomes. Customers on the Block and Swing rate did see 6% lower monthly bill variation (as intended) than the flat rate case, but at the expense of appreciable bill savings, which were only 3%, comparable to the savings seen by non-participants. Given this finding we recommend that additional research be conducted into how best various bill protection mechanisms can balance minimizing customer bill variation with providing financial incentives commensurate with the flexibility customers provide. We also recommend that customer outcomes be explored across a range of regions using current actual customer and system cost data.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2021

Anisotropic thermal-conductivity degradation in the phase-field method accounting for crack directionality

Dynamic loads applied to metals may lead to brittle or ductile fracture depending on the imposed strain rates, material properties, and specimen geometry. With an increase in the applied velocity, brittle to ductile failure mode transition may also be observed. Furthermore, at high strain rates, ductile fracture may be preceded by shear bands, which are narrow bands of intense plastic deformation, typically accompanied by a significant rise in temperature. Reliable models are needed to predict the response of metals subject to dynamic loads. In particular, capturing the interplay between heat conduction and crack propagation using the phase-field fracture method is still an open research field. To accurately capture the heat transfer physics across crack surfaces, damage models degrading thermal-conductivity are necessary. While isotropic thermal-conductivity degradation models were proposed, they may lead to errors as they are indifferent to crack directionality. In our current work, a new anisotropic approach is proposed in which thermal-conductivity, which depends on the phase-field gradient, is degraded solely across the crack. Here, it is shown that this approach improves the near-field approximation of temperature and heat flux compared with isotropic degradation, when taking the discontinuous crack solutions as reference. Additionally, the anisotropic degradation approach is implemented in a unified dynamic fracture model and a couple of examples demonstrate the viability of this approach.

42 ENGINEERING↗
DOE OSTI2026

Experimental characterization and analysis of phase change material-based thermal energy storage system for refrigerated display case

Refrigerated display cases that are used to store and exhibit food products in supermarkets and retail spaces consume a significant portion of the buildings’ total electricity. Importantly, the refrigeration-related energy cost and demand charges are greatly affected by the time-of-use electricity pricing and demand rates, which are at their maximum during peak hours, typically when refrigeration energy consumption is also high. Using energy storage to shift the refrigeration load from peak to off-peak hours can greatly reduce the operational costs in the supermarket. This study demonstrates a phase change material-based thermal energy storage (TES) system, specifically designed in stackable units, that can be integrated with an open vertical refrigerated display case. We perform numerical and experimental characterization that includes finite-difference modeling for the TES, prototype fabrication, and laboratory evaluation, followed by a preliminary system-level analysis to predict the impact of TES on the refrigerated case performance, energy use, and energy cost. The results show that the dedicated latent TES for refrigerated cases can provide a specific energy of 50.4 Wh/kg and a specific power of 15.5 W/kg. The TES can be charged during 12 h of the off-peak period and discharged at various rates during 4 to 6 h of the peak period, thereby shifting the refrigeration load from the peak to the off-peak period. Consequently, annual cost savings up to 19% can be achieved, depending on the thermal load, the summer/winter peak electricity pricing, and the transition temperature of the phase change material used.

25 ENERGY STORAGE↗
DOE OSTI2020

Hydrogen production from lignocellulosic hydrolysate in an up-scaled microbial electrolysis cell with stacked bio-electrodes

Hydrogen production from renewable resources via microbial electrolysis cells (MECs) is a promising approach for sustainable energy production. Yet high hydrogen yield from real feedstocks has not been demonstrated in up-scaled MECs. In this study, a 10-L single chamber MEC with a high electrode surface area to volume ratio (66 m 2 /m 3 ) was constructed and electroactive cathodic biofilms were enriched for hydrogen evolution reaction. A high hydrogen yield of 91% was achieved using lignocellulosic hydrolysate with a hydrogen production rate of 0.71 L/L/D at an organic loading rate of 0.4 g/D. The anodic and cathodic microbial communities, with Enterococcus spp. as the known electroactive bacteria, were capable of achieving current densities of 13.7 A/m 2 and 16.5 A/m 2 , respectively. Finally, a machine learning algorithm was used to investigate the correlation between community data and electrochemical performance, and the critical genera on determining current density were identified.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Decomposition characteristics of C4F7N-based SF6-alternative gas mixtures

C4F7N [2,3,3,3-tetrafluoro-2-(trifluoromethyl)propanenitrile]/CO2 gas mixtures are being developed as an eco-friendly electrical insulator to replace SF6, the most potent greenhouse industrial gaseous dielectric. However, recent studies have reported complicated and often conflicting decomposition pathways for C4F7N/CO2 gas mixtures, which has raised concerns. In this work, the decomposition characteristics of C4F7N/CO2 gas mixtures were studied comprehensively by both designed computations and experiments. Computations were performed starting from fundamental propositions of C4F7N/CO2 decompositions, which were further experimentally verified by pyrolysis, long-term thermal aging with/without catalytic materials (industrial-grade molecular sieves 4A), and electrical decomposition by spark discharge. The results of both computations and experiments suggest that in an ideal thermal decomposition, C4F7N is likely to decompose into C2F6 and small fluoronitriles first at high temperatures. The generation of C3F6 and C2N2 from C4F7N thermal decomposition at lower temperatures appears because of the catalytic effect of incompatible materials, for example, the industrial-grade molecular sieves 4A that we tested. The electron impact dissociation of C4F7N plays an important role in C4F7N electrical decomposition, leading to additional formation of distinctive small molecules of CF4 and C2N2 of low concentrations. It was pointed out based on a real arcing test in a load disconnector that the decomposition of C4F7N gas mixtures in real applications will be at a much moderate and manageable rate than what was obtained from the highly accelerated laboratory tests presented in this work. The signatures of decomposition products extracted in this study provide invaluable guidance for developing decomposition-based diagnosis and fixation of decomposition byproducts toward SF6-free power grids.

Physics↗
DOE OSTI2025

Micron Size NaCrO 2 Particles Enable High‐loading Dry‐processed Electrode for Sodium Ion Batteries

Dry-process fabrication using fibrillatable binder is emerging as a promising method to produce high-loading electrodes for energy storage applications, favored by its cost-efficiency and eco-friendliness. While previous studies have demonstrated the advantages of dry process over the traditional slurry method, there remains a gap in understanding how the particle size of active materials influences the mechanical and electrochemical performance of dry electrodes. In this study, four different particle size NaCrO 2 materials (Average size, S-NCO: 0.6 µm, M1-NCO 1.5 µm, M2-NCO: 4.4 µm, and L-NCO: 9.9 µm) are synthesized to investigate the effect of particle size on dry-processed high-loading electrodes. The findings reveal that the larger micron-sized (>4.4 µm) NCO dry films exhibit significantly improved tensile strength and electrochemical performance, primarily ascribed to the low film porosity, abundant inter-particle connection by the binder, comprehensive carbon coverage, and efficient percolation of the conductive pathway. Notably, a full cell incorporated with a high loading (5.2 mAh cm −2 ) and high active material ratio (96.5 wt.%) L-NCO film electrode demonstrates promising cycling stability and rate capability. Furthermore, these results provide valuable insights regarding the design and fabrication of dry-processed electrodes for future energy storage applications.

25 ENERGY STORAGE↗
DOE OSTI2023

Modeling Electrochemical and Rheological Characteristics of Suspension-Based Electrodes for Redox Flow Cells

Flowable suspension-based electrodes (FSEs) have gained attention in recent years, as the integration of solid materials into electrochemical flow cells can offer improved performance and flexible operation. However, under conditions that engender favorable electrochemical properties (e.g., high particle loading, high conductivity, high surface area), FSEs can exhibit non-Newtonian characteristics that impose large pumping losses and flow-dependent transport rates. These multifaceted trade-offs motivate the use of models to broadly explore scaling relationships and better understand design rules for electrochemical devices. To this end, we present a one-dimensional model, integrating porous electrode theory with FSE rheology as well as flow-dependent electron and mass transport under pressure-driven flow. We study FSE behavior as a function of material properties and operating conditions, identifying key dimensionless groups that describe the underlying physical processes. We assess flow cell performance by quantifying electrode polarization and relative pumping losses, establishing generalized property-performance relationships for FSEs. Importantly, we expound relevant operating regimes—based on a subset of dimensionless groups—that inform practical operating envelopes, ultimately helping to guide FSE and cell engineering for electrochemical systems.

25 ENERGY STORAGE↗
DOE OSTI2026

Measurement Adequacy for Monitoring Data Center Oscillations

Artificial intelligence (AI) training data centers with periodic load profiles can induce sustained grid oscillations across a wide frequency range, making accurate monitoring essential for reliable power system operation. This report evaluates the adequacy of existing measurement systems for monitoring such oscillations, focusing on phasor measurement units (PMUs) and point-on-wave (POW) measurement systems. The analysis shows that while PMUs are highly effective for monitoring low-frequency electromechanical oscillations, they have inherent limitations in accurately representing higher-frequency oscillations due to constraints imposed by reporting rates and the bandwidth of phasor estimation filters. Even when configured with higher reporting rates, the filtering inherent in the phasor estimation process can significantly attenuate oscillation magnitudes, potentially leading to underestimation of oscillatory behavior. This has important implications for compliance and performance monitoring of large loads. To address the limitations associated with PMU-based monitoring, the report examines the use of high-resolution POW measurements, which can capture oscillations across a broader frequency range. However, continuous POW monitoring introduces practical challenges related to large data volumes, communication bandwidth, and real-time data processing. For this reason, the report also discusses emerging approaches that use POW measurements as a complementary capability alongside PMUs to improve observability of oscillations from large data center loads.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

The Global Fingerprint of Modern Ice-Mass Loss on 3-D Crustal Motion

Crustal motion generated by rapid ice-mass loss from Earth's glaciers and ice sheets has previously been considered in Global Navigational Satellite System (GNSS) analyses and numerical models across regions of ice retreat. However, the fingerprint of ice-mass loss is not limited to glaciated areas, but is characterized by a global pattern of 3-D crustal deformation. We compute “far-field” vertical and horizontal deformation rates that occurred in response to early 21st century mass flux from the Greenland and Antarctic Ice Sheets, global glaciers and ice caps, and associated ocean loading. We demonstrate that mass changes in the Greenland Ice Sheet and high latitude glacier systems each generated average crustal motion of 0.1–0.4 mm/yr across much of the Northern Hemisphere, with significant year-to-year variability in magnitude and direction. Horizontal motions associated with ice-mass loss exceed vertical rates in many far-field areas, and both should be considered in future analysis of GNSS measurements.

58 GEOSCIENCES↗
DOE OSTI2026

Characterizing peak electricity demand for U.S. households: an assessment of end-use loads and demand factors

Understanding household peak electricity demand is critical to evaluate the technical need for electrical infrastructure upgrades. This study characterizes peak loads for existing and new equipment using metered data from a convenience sample of 11,940 U.S. dwellings from four sources, including 911 from two sources with end-use metering. After standardized data cleaning and labeling, we derived descriptive statistics for key metrics, such as maximum demand and demand factors, and developed predictive models relating 60- to 15-min demand for the National Electrical Code (NEC). Mean 15-min maximum demand was 9.7 kW (median 9.0 kW; IQR 7.0–11.5 kW, 95% CI 9.6–9.8 kW), indicating spare capacity in 98% of homes with hypothetical 100 A panels. Maximum demand increased with floor area and number of high-demand loads. Dwelling maximum demand was driven by higher-power, longer-duration heating appliances and vehicle charging, while most user-operated appliances contributed little. Demand factors are used to account for how most devices contribute less than their rated power to maximum demand. Existing load mean demand factors (28%; median 10%; IQR 0–58%; CI 28–29%) were higher than those for new loads (21%; median 7%; IQR 0–35%; CI 20–21%), because new loads changed the timing and magnitude of maximum demand. New high-demand loads had higher than average demand factors (40–60%). Whole dwelling demand factors support the NEC's 40% assumption, but they challenge its conservative 100% treatment of new HVAC. We propose a data-driven 50% demand factor for new equipment, which would align with metered data, improve affordability, and modernize electrical codes.

Appliances↗
DOE OSTI2024

High‐Rate Polymeric Redox in MXene‐Based Superlattice‐Like Heterostructure for Ammonium Ion Storage

Abstract Achieving both high redox activity and rapid ion transport is a critical and pervasive challenge in electrochemical energy storage applications. This challenge is significantly magnified when using large‐sized charge carriers, such as the sustainable ammonium ion (NH 4 + ). A self‐assembled MXene/n‐type conjugated polyelectrolyte (CPE) superlattice‐like heterostructure that enables redox‐active, fast, and reversible ammonium storage is reported. The superlattice‐like structure persists as the CPE:MXene ratio increases, accompanied by a linear increase in the interlayer spacing of MXene flakes and a greater overlap of CPEs. Concurrently, the redox activity per unit of CPE unexpectedly intensifies, a phenomenon that can be explained by the enhanced de‐solvation of ammonium due to the increased volume of 3 Å‐sized pores, as indicated by molecular dynamic simulations. At the maximum CPE mass loading (MXene:CPE ratio = 2:1), the heterostructure demonstrates the strongest polymeric redox activity with a high ammonium storage capacity of 126.1 C g −1 and a superior rate capability at 10 A g −1 . This work unveils an effective strategy for designing tunable superlattice‐like heterostructures to enhance redox activity and achieve rapid charge transfer for ions beyond lithium.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Lithium promotion of Pt/m-ZrO 2 catalysts for low temperature water-gas shift

Low temperature water-gas shift (LTS) is an important reaction occurring in a fuel processor for producing and purifying hydrogen. Platinum supported on m-ZrO 2 belongs to a family of catalysts consisting of metal nanoparticles and an active partially reducible oxide, with the catalysis proposed to occur at the boundary between metal particles and the support. In this investigation, increasing the loading of lithium dopant increased the LTS rate up to 0.54 wt % lithium, where conversion was 2.4 times that of the unpromoted catalyst at 260°C. Further increases in lithium loading up to 1.5 wt % decreased the rate, although it remained higher than that of the unpromoted catalyst. Infrared spectroscopy and CO 2 temperature programmed desorption experiments showed three effects with increasing lithium loading: (1) lithium promoter weakened the C–H bond of formate, the proposed rate limiting step of the interfacial surface formate mechanism; (2) high levels of lithium suppressed the platinum site capacity required for hydrogen transfer; and (3) high levels of lithium increased catalyst basicity. Aspects (2) and (3) tended to inhibit desorption of product CO 2 , an acidic molecule the removal of which is metal-catalyzed. XANES and XPS experiments revealed that electron transfer to enrich Pt nanoparticles is unlikely the root cause of C–H bond weakening in formate. However, other electronic effects (e.g., electrostatic effects or molecular rearrangement due to enhanced basicity) were not ruled out.

08 HYDROGEN↗
DOE OSTI2021

Anisotropic strength behavior of single-crystal TATB

High-rate strength behavior plays an important role in the shock initiation of high explosives, with plastic deformation serving to localize heat into hot spots and as a mechanochemical means to enhance reactivity. Recent simulations predict that detonation-like shocks produce highly reactive nanoscale shear bands in the layered crystalline explosive TATB (1,3,5-triamino-2,4,6-trinitrobenzene), but the thresholds leading to this response are poorly understood. We utilize molecular dynamics to simulate the high-rate compressive stress–strain response of TATB, with a focus on understanding flow behavior. Here, the dependence of strength on pressure and loading axis (crystal orientation) is explored. The deformation mechanisms fall broadly into two categories, with compression along crystal layers activating a buckling/twinning mode and compression normal to the layers producing nanoscale shear bands. Despite the complexity of the underlying mechanisms, the crystal exhibits relatively straightforward stress–strain curves. Most of the crystal orientations studied show rapid strain softening following the onset of yielding, which settles to a steady flow state. Trajectories are analyzed using five metrics for local states and structural order, but most of these metrics yield similar distributions for these deformation mechanisms. On the other hand, a recently proposed measure of intramolecular strain energy is found to most cleanly distinguish between these mechanisms, while also providing a plausible connection with mechanochemically accelerated decomposition kinetics. Localization of intramolecular strain energy is found to depend strongly on crystal orientation and pressure.

36 MATERIALS SCIENCE↗