Search NASA⌕ Search

SEARCH · Search NASA

Results for “Joule heating”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Three-dimensional feedback processes in current-driven metal

Using three-dimensional (3D) magnetohydrodynamic simulations, we study how a pit on a metal surface evolves when driven by intense electrical current density j. Redistribution of j around the pit initiates a feedback loop: j both reacts to and alters the electrical conductivity σ, through Joule heating and hydrodynamic expansion, so that j and σ are constantly in flux. Thus, the pit transforms into larger striation and filament structures predicted by the electrothermal instability theory. Both structures are important in applications of current-driven metal: Here, the striation constitutes a density perturbation that can seed the magneto-Rayleigh-Taylor instability, while the filament provides a more rapid path to plasma formation, through 3D j redistribution. Simulations predict distinctive self-emission patterns, thus allowing for experimental observation and comparison.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Pulsed studies of intervalley transfer in Al 0.35 In 0.65 As: A paradigm for valley photovoltaics

Al 0.35 In 0.65 As is a direct semiconductor whose multivalley conduction-band structure has been proposed for use in so-called valley photovoltaics. In such hot-carrier solar cells, energetic (hot) photocarriers are stored in satellite valleys away from the Γ point, allowing them to be extracted prior to thermalization and to thereby increase power-conversion efficiency. While prior theoretical work has highlighted the potential of Al 0.35 In 0.65 As—a widely used barrier material in electronic and optoelectronic devices, for use in valley photovoltaics—surprisingly little is known about its electrical properties, especially how these are impacted by the application of high fields. In this work, we therefore undertake a detailed characterization of the electrical properties of Te-doped ($\textit{n}$-type) Al 0.35 In 0.65 As, over wide ranges of temperature (3–400 K) and electric field (<50 kV/cm). Using pulsed measurements to suppress the influence of Joule heating, we reveal the presence of clear negative-differential conductance in the current-voltage characteristics of the films, suggestive of the intervalley transfer of hot electrons. This conclusion is supported by the results of ensemble Monte Carlo simulations of the hot-carrier action, which confirm the connection of the observed negative-differential conductance to hot-electron transfer from the conduction-band (Γ) minimum, to the side valleys at the $\textit{L}$ point. Further, the quantitative features of the experimentally determined velocity-field curves are found to be in good agreement with the results of these calculations, providing further confidence in the role of the implied intervalley transfer mechanism. Overall, these results confirm the excellent potential of Al 0.35 In 0.65 As for use as the absorber material in hot-carrier solar cell technology.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Low Offset and Noise in High Biased GaN 2DEG Hall-Effect Plates Investigated With Infrared Microscopy

This article presents GaN two-dimensional electron gas (2DEG) Hall plates with low residual offset and noise at 3 V input bias. We studied devices made from three consecutive fabrication generations through current spinning offset measurements in a zero-field chamber. When operated above 1 V, the first-generation devices charted high residual offsets >1 mT. We reduced these residual offsets by three orders of magnitude in later device generations. Two experiments were performed to confirm the improvements. First, the GaN 2DEG Hall-effect plates were measured with infrared microscopy during current spinning to investigate current crowding and Joule heating. The low offset devices had minimal thermal gradients while the high offset devices had large thermal variations during current spinning. Second, zero-field current spinning was used to compare GaN 2DEG Hall-effect plates against typical commercial silicon Hall-effect plates. The Si devices had residual offsets >10μT biased at 3 V and AlGaN/GaN devices had residual offset <4 μT at 3 V, and an InAlN/GaN devices operated at 2 V had a residual offset <4 μT . Thus, for the first time, GaN 2DEG Hall-effect plates were shown to operate with high bias with low offsets and noise, which enables high fidelity sensing applications in extreme environments.

2DEG↗
DOE OSTI2022

Combined Normal and Disruption, Electromagnetic Transient, Thermal, and Structural Analysis of COMPASS Upgrade

We report that COMPASS Upgrade is a replacement for the COMPASS device at the Institute of Plasma Physics (IPP) of the Czech Academy of Sciences. It is a cryogenic copper machine with a major radius of 0.89 m, with 5 T at the plasma centerline and 2-MA plasma current. It is intended to develop ITER relevant plasma densities and high power fluxes in the divertor region. The entire vacuum vessel is planned to reach 500 °C. The COMPASS-U TF coil is a bitter plate-like design with radial conductor plates that extend the width of the TF inner leg. Currents redistribute radially based on resistive and inductive effects. This necessitates an electromagnetic (EM) transient analysis of the coil current distribution. For copper coils like COMPASS Upgrade, the solution must be coupled with a thermal solution to properly include the temperature-dependent resistive effects. The EM transient solution includes the TF magnetic loads, and a stress pass on the results produces in-plane loads. The EM solution chosen requires inclusion of all the conducting and nonconducting regions surrounding the TF coil. With the TF current distribution solved, the addition of the poloidal coils and resulting background fields allows the determination of the out-of-plane loads. Integration of the Joule heating produced from the EM solution produces temperature distributions throughout the shot, which can be read in to a structural pass to include thermal stresses in the evaluation. This is true of both TF and PF temperature prediction with appropriate packing fractions. Currents in the PF coils, with computed background fields, produce Lorentz loads. Time transients of the PF currents produce startup eddy currents in the structures included in the model; inclusion of the vessel and other passive structures allows these to be included in a stress pass. The EM transient solution includes all the elements needed for a disruption analysis with prescribed motions and quench of the plasma. From a model originally intended to simulate current diffusion in the TF coil, the analysis can be refined to be used to predict the performance of most of the components of the tokamak or provide boundary conditions for more detailed submodels. In this article, the use of the current diffusion model for bounding assessments of the COMPASS Upgrade TF, PF, vessel, external structure, and disruption analyses will be presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Effect of water vapor and thermal history on nuclear waste feed conversion to glass

Water affects the glass melting process by interacting with the foam layer at the glass melt surface and by influencing the batch conversion reactions. Water vapor from the nuclear waste slurry feed maintains a high water vapor pressure in the melter atmosphere. To investigate to what extent water vapor affects the vitrification of nuclear waste in joule-heated, cold-top melters, a series of feed expansion experiments were performed under humid and dry atmospheres using samples of low-activity waste (LAW) melter feed simulants. Melting of feed pellets in the presence of water vapor slightly decreased the temperature of primary foam onset, but did not significantly affect the feed volume expansion by foaming or the foam-collapse temperature. Sets of feed expansion experiments and evolved gas analyses were also performed to check the effect of thermal history on LAW feed samples tested as direct slurry, loose powder, and slow- and fast-dried pellets.

Marcial, Jose↗
DOE OSTI2020

High Voltage Surge Protection by a Varistor-Filled Air Gap

The electric discharge across a varistor granule filled air gap under a fast-rising voltage pulse was investigated for surge protection applications. The effects of temperature and pressure on the arc and the electrical conduction were analyzed by the characteristic changes in voltage waveforms triggered by a fast-rising high voltage pulse. In addition to the gap size, experimental results show that competing mechanisms among arc conduction, conduction through the varistor granule network, thermionic emission from Joule heating at granule-to-granule contact points, and the magnitude of the switching voltage dictate the maximum surge protection voltage for the filled air gap. Experimental evidence indicated that accumulated degradation was created at small contact points between varistor granules by repetitive assaults from longer duration, high voltage pulses. The uniqueness of using varistor over other dielectric granules in an air gap for surge protection is identified and discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Impact of non‐framework cation mixing on the structure and crystallization behavior of model high‐level waste glasses

Abstract Spinel crystallization is known to be detrimental to the operation of Joule heated ceramic melters during the vitrification of iron‐rich high‐level nuclear wastes (HLW) into borosilicate glasses. The literature on this subject focuses on tackling the problem by developing empirical constraints to design compositions, which limit the fraction of spinels formed in the melter or by developing empirical models to predict the settling behavior of spinels in the melter as a function of the glass composition. While these empirical models can predict the behavior of most of the compositions, they are not failsafe as there are always some compositions, whose behavior is beyond the predictive ability of these models. This can lead to undesirable situations during the vitrification of the nuclear waste, and therefore an in‐depth investigation of the chemo‐structural descriptors controlling the crystallization behavior in these glasses is warranted. Accordingly, the present study aims to understand the impact of non‐framework cation mixing (i.e., Li + /Na + and Ca 2+ /Na + ) on the structure (through Raman spectroscopy and Mössbauer spectroscopy) and crystallization behavior (through XRD, SEM‐EDS, and vibrating sample magnetometry) of iron‐rich model HLW glasses in the system: (mol.%) x M y O–(25− x ) Na 2 O–9.12 B 2 O 3 –6.4 Al 2 O 3 –51.25 SiO 2 –7.22 Fe 2 O 3 –0.38 MnO–0.08 Cr 2 O 3 –0.55NiO (M y O = Li 2 O or CaO).

Balasubramanya, Nikhila↗
DOE OSTI2021

Modeling Nanosecond-Pulsed Spark Discharge and Flame Kernel Evolution

Dilute combustion, either using exhaust gas recirculation or with excess air, is considered a promising strategy to improve the thermal efficiency of internal combustion engines. However, the dilute air-fuel mixture, especially under intensified turbulence and high-pressure conditions, poses significant challenges for ignitability and combustion stability, which may limit the attainable efficiency benefits. In-depth knowledge of the flame kernel evolution to stabilize ignition and combustion in a challenging environment is crucial for effective engine development and optimization. To date, a comprehensive understanding of ignition processes that result in the development of fully predictive ignition models usable by the automotive industry does not yet exist. Spark-ignition consists of a wide range of physics that includes electrical discharge, plasma evolution, joule-heating of gas, and flame kernel initiation and growth into a self-sustainable flame. In this study, an advanced approach is proposed to model spark-ignition energy deposition and flame kernel growth. To decouple the flame kernel growth from the electrical discharge, a nanosecond-pulsed high-voltage discharge is used to trigger spark-ignition in an optically accessible small ignition test vessel with a quiescent mixture of air and methane. Initial conditions for the flame kernel, including its thermodynamic state and species composition, are derived from a plasma-chemical equilibrium calculation. The geometric shape and dimension of the kernel are characterized using a multi-dimensional thermal plasma solver. Here, the proposed modeling approach is evaluated using a high-fidelity computational fluid dynamics procedure to compare the simulated flame kernel evolution against flame boundaries from companion Schlieren images.

42 ENGINEERING↗
DOE OSTI2022

Multi-principal elemental intermetallic nanoparticles synthesized via a disorder-to-order transition

Nanoscale multi-principal element intermetallics (MPEIs) may provide a broad and tunable compositional space of active, high–surface area materials with potential applications such as catalysis and magnetics. However, MPEI nanoparticles are challenging to fabricate because of the tendency of the particles to grow/agglomerate or phase-separated during annealing. Here, we demonstrate a disorder-to-order phase transition approach that enables the synthesis of ultrasmall (4 to 5 nm) and stable MPEI nanoparticles (up to eight elements). We apply just 5 min of Joule heating to promote the phase transition of the nanoparticles into L1 0 intermetallic structure, which is then preserved by rapidly cooling. This disorder-to-order transition results in phase-stable nanoscale MPEIs with compositions (e.g., PtPdAuFeCoNiCuSn), which have not been previously attained by traditional synthetic methods. This synthesis strategy offers a new paradigm for developing previously unexplored MPEI nanoparticles by accessing a nanoscale-size regime and novel compositions with potentially broad applications.

25 ENERGY STORAGE↗
DOE OSTI2021

Spatiotemporal characterization of the field-induced insulator-to-metal transition

Watching a metal filament grow Resistive switching is a process in which the electrical resistance of a sample changes abruptly in response to a voltage pulse, often by orders of magnitude. This process is at the heart of many neuromorphic computing approaches but visualizing it in both space and time is tricky. del Valle et al . monitored the resistive switching in three different vanadium oxide compounds by measuring time- and space-resolved optical reflectivity (see the Perspective by Hilgenkamp and Gao). A characteristic conducting filament was quickly nucleated on the inhomogeneities in the sample and then propagated due to Joule heating. —JS

Science & Technology - Other Topics↗
DOE OSTI2023

A New Analytical Expression for Estimating the Adiabatic Temperature Rise in Lithium-Ion Batteries During High-Power Pulses

A coupled, thermal-electrochemical model is used to explain why the Joule heating assumption (i.e., I 2 R) does not provide a good representation of the temperature rise during high power pulses in lithium-ion batteries, even in cases where the reversible heat generation can be neglected. The poor agreement occurs because the internal resistance changes during the pulse due to the opposing effects of mass transport in the electrolyte, which raises the resistance, and heat generation (temperature rise), which lowers the resistance. These insights are used to propose a new analytical expression for predicting the temperature rise during adiabatic pulses with limited experimental and physical data. The expression accounts for thermal effects using a Taylor series expansion of an Arrhenius-type equation. It accounts for transport effects using an approximate solution to the one-dimensional diffusion equation. The new expression is shown to accurately estimate the simulated, adiabatic temperature rise across a range of loadings (1 to 4 mAh cm -2 ) and C-rates (1C to 10C) for cells containing LiNi 0.5 Mn 0.3 Co 0.2 O 2 positive electrodes and graphite negative electrodes. In conclusion, it is also shown to accurately estimate the experimental, adiabatic temperature rise measured for cells with LiFePO 4 positive electrodes and mesocarbon negative electrodes with minimal changes to the fitted parameters.

25 ENERGY STORAGE↗
DOE OSTI2020

Coupling between interfacial charge and mechanical deformation at high temperatures in ceramics. Final report

Basic research, under this project, has led to new fundamental discoveries in the effect of electric field (and temperature) on diffusion kinetics, electroluminescence and electronic conductivity in ceramic materials. For example, yttria stabilized zirconia sinters in just seconds well below 1000°C, while, at the same time, exhibiting a highly non-linear rise in conductivity, and intense electroluminescence. The phenomenon has been shown to occur in nearly all ceramics whether they are ionic or electronic conductors, or insulators. The rise in conductivity leads to Joule heating but in-situ measurements with a platinum standard at APS and Brookhaven (NSLS-II) have shown that the temperatures are much below what would have been needed to sinter in a few seconds. We have explored the hypothesis that huge concentrations of Frenkel pairs of vacancies and interstitials are produced. While these point defects accelerate sintering (by the interstitials being transported to the pores and the vacancies to the grain boundaries), they may also produce electron-hole pairs which can impart high electrical conductivity or recombine to emit photons. Many measurements, including in-situ measurements of lattice expansion are consistent with this hypothesis. However, the energy barrier for generating Frenkels in zirconia for example is much larger than can be expected at temperature and electrical fields in these experiments. Most recent work suggests that short wavelength phonons play a critical role in the genesis of the “flash” phenomenon. For example, the Debye temperature has been shown to be the lower bound for initiating flash. Molecular dynamics simulations have shown that large concentrations of Frenkel pairs can form under the proliferation of “optical” phonons at the edge of the Brillouin zone, that is with a wavelength at or below the lattice parameter. This project has led to the development of new experimental techniques as well as a software hardware interface for process and control and data analysis in real time, with a time scale of less than one second. This approach spells a new paradigm in materials science research. It potential impact on high rate manufacturing is self evident.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Vitrification of High-Cr Glass in Research-Scale Melter

This test report describes the experimental results from a small-scale test using the research-scale melter (RSM) at Pacific Northwest National Laboratory to demonstrate processing of a high-Cr simulated feed stream, designated HLW-HCr-16. The RSM is a small, joule-heated melter capable of processing melter feed continuously. The melter is equipped with Inconel ® 693 electrodes, Monofrax ® K-3 refractory, and an Inconel 690 pour spout. An electric kiln surrounds the melter body and minimizes heat loss from the melter body during operation. The RSM is equipped with an offgas treatment system that employs quenching, wet scrubbing, and high-efficiency mist elimination. The glass-discharge section is heated to facilitate pouring of the glass. The melter is fitted with a melt cavity that is ~25 cm (10 in.) in diameter with a nominal glass depth of 8.9 cm (3.5 in.). The melter was operated with a target glass temperature of 1150°C and target plenum temperature between 550°C and 700°C for this test. The air bubbling rate was 4.2 L/min. Overall, during the continuous operation of the melter for ~ 103 hours, ~ 141 kg of glass was produced. At the conclusion of the test, the melter and exhaust lines were visually inspected for particulate deposition and corrosion. Entrained material had adhered to the underside of the melter lid and to the exhaust piping. Enrichments in elements such as Cl, F, B, K, Li, P, Na, and S were measured in these deposits through inductively coupled plasma–optical emission spectroscopy and X-ray fluorescence analysis. When the melter electrodes and air bubbler tube were removed from the glass in the RSM, the electrodes appeared discolored, but no significant loss of metal was observed. The processing of a high-Cr simulant, HLW-HCr-16, in the RSM produced glass at an average rate of 1.36 kg/h, equaling a melter-surface-area normalized glass generation rate of 654 kg/day/m 2 . The resulting glass met the toxicity characteristic leaching procedure requirement. Test results of crystallinity, electrical conductivity, and viscosity showed good processing properties of this high-Cr high-level waste glass. RSM offgas was also sampled and analyzed at periodic intervals during steady-state operating conditions. The total decontamination factor averaged by four sampling periods was 134. The concentrations of CO and NO x in emissions were 237 to 422 and 69 to 94 parts per million by volume, respectively.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Development of Quantum Computing Interconnect Based on Aerosol Jet Printing and Electrochemical Deposition of Rhenium

The electrodeposition of rhenium on to a metal seed layer on flexible substrates is presented as a means to creating superconducting flexible cable connectors in an enabling plug-and-play approach for quantum computing. Cryogenic quantum electronics are currently connected using masses of stainless-steel coaxial cables that are bulky, rigid - both in form and design - and lead to significant joule heating, thermal noise, and cross talk. Here, we present an unprecedented approach to integrate an aerosol jet printed (AJP) metal seed layer with rhenium electrodeposition on a flexible substrate in the advancement of superconducting interconnect technologies. Silver and gold were printed using the ‘Nanojet’ aerosol jet printer on Kapton films. Adhesion of gold was found to be far superior to that of silver and adhesion on roughened Kapton surpassed that of its smooth counterpart. Electrodeposition of rhenium was successful on both silver and gold and an amorphous Re film was confirmed by XRD. Results for both materials are presented however due to the poor adhesion of silver to Kapton it was ruled out as a viable candidate. Composite materials were characterized by profilometry, EDS, XRD and FIBSEM. Electrical measurements of the composite at ambient temperature showed a critical temperature (T c ), where the resistance drops to 0, of 5.8 K, well above 4.2 K, the temperature of liquid helium. Stress-strain tests of the Ag-Re and Au-Re composites on roughened and smooth Kapton were completed. Cyclic flexure testing (200 cycles) to 1.25% strain showed smooth Kapton samples reach a stress of ~16 MPa, while Kapton roughened with sandpaper, reaches ~20MPa of stress for the same 1.25% strain.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Eddy current studies for the beam screen of the Electron-Ion Collider hadron storage ring superconducting magnets

A copper-clad stainless steel beam screen will be installed in each of the Relativistic Heavy Ion Collider (RHIC) superconducting (SC) magnets used for making the Electron-Ion Collider (EIC) Hadron Storage Ring (HSR). Eddy currents are induced within conductors to oppose a changing magnetic flux. Eddy currents will appear on the EIC HSR beam screens during ramp up and ramp down, after a magnet quenches and when reversing the polarity of the γ transition jump quadrupoles during γ transition crossing. This study evaluates the magnitude of the eddy currents, the temperature increase due to Joule heating and the stress at the beam screens from the eddy-current induced force.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Compositional Variation Tests on DuraMelter 100 with LAW Sub-Envelope C1 Feed (LAWC22) Glass in Support of the LAW Pilot Melter (Final Report)

The Vitreous State Laboratory at The Catholic University of America (VSL) is developing and testing glass formulations for RPP-WTP waste envelopes to provide data to meet the RPP-WTP contract requirements and to support system design activities. That work is based upon small-scale batch melts (“crucible melts”) using waste envelope simulants. Selected formulations have also been tested in small-scale continuously-fed joule-heated melters (DM10 and DM100 systems) and, ultimately, in the LAW Pilot Melter. Such melter tests provide information on key process factors such as feed processing behavior, dynamic effects during processing, sulfate incorporation, processing rates, off-gas amounts and compositions, foaming control, etc., that cannot be reliably obtained from crucible melts. This sequential scale-up approach in the vitrification testing program ensures that maximum benefit is obtained from the more costly pilot-scale tests and that the most effective use is made of that resource. The principal objective of the work described in this report was to collect the necessary small-scale melter test data with a LAW Sub-Envelope C1 waste simulant in order to support the next series of tests with the same melter feed on the LAW Pilot Melter. This work was conducted under a corresponding Test Specification and Test Plan.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Exploration of signal processing methods for superconducting magnet and quench data

Quenching is the phenomenon of a superconducting magnetic material carrying current transitioning into a regular conducting material. This may cause severe and irreparable damage to the superconductor due to Joule heating. The Magnet Department at Fermi National Accelerator Laboratory (FNAL) has acquired experimental data through quench antenna arrays that are recorded when the quench is detected. These data are in terms of voltage signals that are sampled at 100kHz for several minutes. There are multiple channels and each channel provides a data set of more than 20 million observations, while there is one channel, called the trigger channel which shows the time when quench is detected. Despite some advancements that were made including machine learning, data complexity still shadows the progress. In this work, we studied a multi-resolution analysis of the quench antenna data through the Haar wavelet transform. In particular, we applied the maximally overlapped discrete w avelet transform (MODWT) of a suitable level L to the given data and then projected it onto the wavelet basis. This decomposes a given signal (Original data) $x ϵ \mathbb{R}^N$ into $L + 1$ subspaces of $\mathbb{R}^N$. One of the subspaces called the approximation, captures the trend of the signal, and the others, called the details, capture the fluctuations at different frequency bands. This decomposition provides a clear trend of the data at a suitable level and also various activities (spikes) are seen in the details of the decomposition at every level. These spikes might reveal some information about the quench under investigation but in any case, give information about magnet behavior. Also, this decomposition is seen to be very useful in removing noise present in the data due to the source or mechanism of the experiment.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

The Influence of Carbon Fiber Composite Specimen Design Parameters on Artificial Lightning Strike Current Dissipation and Material Thermal Damage

Previous artificial lightning strike direct effect research has examined a broad range of specimen design parameters. No works have studied how such specimen design parameters and electrical boundary conditions impact the dissipation of electric current flow through individual plies. This article assesses the influence of carbon fiber composite specimen design parameters (design parameters = specimen size, shape, and stacking sequence) and electrical boundary conditions on the dissipation of current and the spread of damage resulting from Joule heating. Thermal-electric finite element (FE) modelling is used and laboratory scale (<1 m long) and aircraft scale (>1 m long) models are generated in which laminated ply current dissipation is predicted, considering a fixed artificial lightning current waveform. The simulation results establish a positive correlation between the current exiting the specimen from a given ply and the amount of thermal damage in that ply. The results also establish that the distance to ground, from the strike location to the zero potential boundary conditions (ground), is the controlling factor which dictates the electric current dissipation in each ply. Significantly, this distance to ground is dependent on each of the specimen shape, dimensions, stacking sequence, and location of ground boundary conditions. Therefore, it is not possible to decouple current dissipation and damage from specimen design and boundary condition setup. However, it is possible to define a specimen size for a given specimen shape, stacking sequence, and waveform which limit the influence of specimen dimensions on the resulting current distribution and damage. For a rectangular specimen design which appears in literature multiple times, as 100 × 150 mm and with a stacking sequence of [45/0/-45/90] 4s , a specimen design of greater than 300 × 200 mm is required to limit the influence of specimen dimensions on current distribution and damage.

36 MATERIALS SCIENCE↗