Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “process cooling”

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 91 records · Page 5

Advanced CO 2 Capture Solvent Systems for Dynamic Power Generation: Quarterly Research Performance Progress Report, QR4 (Q4FY24)

We developed an integrated Computational Fluid Dynamics (CFD) model to simulate the multi-physics coupled cooling process of mixed gas by cold water within a Direct Contact Cooler (DCC) equipped with a rotating packing bed (RPB). The model captures the interactions between fluid dynamics, heat transfer, mass transport, and phase transitions, while accounting for key operational variables such as RPB rotational speed and the mass flow rates of both liquid and gas. The CFD model has been validated using experimental data, specifically by comparing predicted outflow gas and liquid temperatures to measured results. Our findings demonstrate the significant effects of RPB rotational speed and mass flow rates on cooling performance, providing valuable insights for optimizing DCC efficiency in industrial applications.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Understanding formation mechanisms of intermetallic compounds in dissimilar Al/steel joint processed by resistance spot welding

Here, this paper confirmed the formation mechanism of intermetallic compounds (IMCs) in Al/steel resistance spot welds with transmission electron microscopy, electron backscatter diffraction, nanoindentation and thermodynamic calculations. In particular, the formation of AlFe with BCC_B2 structure, which is not commonly seen in welds, was identified. The formation mechanism of IMCs at the high welding energy region is described as follows. Al 13 Fe 4 first nucleates from the Al side, followed by Al 5 Fe 2 growth with Fe atoms accumulating at the Al 13 Fe 4 grain boundaries. Then, Al 5 Fe 2 grains grow continuously to coarse columnar grains, and AlFe forms at the interface of Al 5 Fe 2 grains and the ferrite phase. Lastly, the needle-like Al 13 Fe 4 forms in the cooling process. At the middle welding energy region, only equiaxed Al 5 Fe 2 and small Al 13 Fe 4 grains are formed at the interface because of the lower diffusion rates of Al and Fe, hence postponing the growth of Al 5 Fe 2 and Al 13 Fe 4 . At the low welding energy region, only sporadic Al 5 Fe 2 and Al 13 Fe 4 grains are formed surrounded by the Al phase.

42 ENGINEERING↗

Thermodynamic Spectrum of Solar Flares Based on SDO/EVE Observations: Techniques and First Results

The Solar Dynamics Observatory (SDO)/EUV Variability Experiment (EVE) provides rich information on the thermodynamic processes of solar activities, particularly on solar flares. Here, we develop a method to construct thermodynamic spectrum (TDS) charts based on the EVE spectral lines. This tool could potentially be useful for extreme ultraviolet (EUV) astronomy to learn about the eruptive activities on distant astronomical objects. Through several cases, we illustrate what we can learn from the TDS charts. Furthermore, we apply the TDS method to 74 flares equal to or greater than the M5.0 class, and reach the following statistical results. First, EUV peaks are always behind the soft X-ray (SXR) peaks and stronger flares tend to have faster cooling rates. There is a power-law correlation between the peak delay times and the cooling rates, suggesting a coherent cooling process of flares from SXR to EUV emissions. Second, there are two distinct temperature drift patterns, called Type I and Type II. For Type I flares, the enhanced emission drifts from high to low temperature like a quadrilateral, whereas for Type II flares the drift pattern looks like a triangle. Statistical analysis suggests that Type II flares are more impulsive than Type I flares. Third, for late-phase flares, the peak intensity ratio of the late phase to the main phase is roughly correlated with the flare class, and the flares with a strong late phase are all confined. We believe that the re-deposition of the energy carried by a flux rope, which unsuccessfully erupts out, into thermal emissions is responsible for the strong late phase found in a confined flare. Furthermore, we show the signatures of the flare thermodynamic process in the chromosphere and transition region in the TDS charts. These results provide new clues to advance our understanding of the thermodynamic processes of solar flares and associated solar eruptions, e.g., coronal mass ejections.

Wang, Yuming↗

Enriched immersed boundary method (EIBM) for interface-coupled multi-physics and applications to convective conjugate heat transfer

It is challenging to develop numerical methods that simultaneously maintain accuracy of resolving boundary conditions and mesh flexibility to handle the interface in interface-coupled multi-physics problems involving large property discontinuity. On the one hand, boundary-fitted methods possess high accuracy in capturing the interfacial phenomenon but involve complicated volumetric mesh generation, mesh-motion, and even re-meshing procedures. On the other hand, immersed boundary methods (IBM) provide mesh flexibility but sometimes suffer from inferior interface representations, leading to poor enforcement of boundary conditions. This paper presents an enriched immersed boundary method (EIBM) to overcome this challenge and demonstrates its efficacy in conjugate heat transfer, a representative example in interface-coupled multi-physics systems that have implications for many industrial processes. Further, the core technique of the method is to enhance IBM’s accuracy of the fluid–solid interface by enriching the degrees of freedom of the cut elements to enforce temperature and flux compatibilities and resolve all the physical unknowns on the background mesh to simplify the volumetric mesh generation. We implement the EIBM under the framework of a variational multiscale formulation for coupled Navier–Stokes and thermodynamics equations. The enriched DoFs enable better enforcement of temperature and flux compatibilities with large conductivity ratios across the fluid–solid interface. At the same time, the immersed nature of the proposed method still attains mesh flexibility. The EIBM’s accuracy is thoroughly evaluated through a set of examples, ranging from benchmark problems with analytical solutions to real-world cooling processes of a moving metallic structure with complex geometry.

42 ENGINEERING↗

A fundamental study of liquid spray evaporation at a heated surface

The influence of feedwater temperature and the influence of pulsing nozzle flow on the heat transfer characteristics in the spray evaporative cooling process were investigated. A commercial water circulation heater was installed in the feedwater line between the supply reservoir and the nozzle. A pass through the heater provided water at temperatures up to 50 C at the nozzle; higher temperatures were obtained by multiple heater passes. A thermocouple installed in the water line just upstream of the nozzle provided a measure of the feedwater temperature. Initial tests were run with belljar exhaust orifice number two which provided a pressure in a range about 7mm Hg abs. The feedwater temperature was heated to a temperature of 45 C. Results indicate that for a given system operating pressure, there is an upper limit on the feedwater temperature above which the spray will flash to vapor before hitting the surface to be cooled.

Wierum, F. A.↗

Rate of Para-Orthohydrogen Conversion in Cryogenic Vortex Tube

Raising efficiency of the cooling process for cryogenic hydrogen and minimizing hydrogen boil-off during storage and transportation of liquid hydrogen are critically important factors for widening implementation of hydrogen-based clean energy systems. One novel approach involves utilization of relatively simple catalyzed vortex tubes, where hydrogen in the primarily para- nuclear spin isomer form is converted into the ortho- form while consuming a significant amount of heat. To design such systems effectively, better understanding is required about catalyst-assisted para-orthohydrogen conversion rates in high-speed vortical flows. In this study, a computational fluid dynamics simulation has been set up to model an experimental system with cryogenic hydrogen vortex tubes. Mesh-verification and a validation study has been conducted first for a non-catalyzed tube. Then, the rate coefficient of the para-orthohydrogen conversion of cryogenic hydrogen has been determined by matching numerical results with experimental data available for a catalyzed vortex tube. Furthermore, the presented information can help design and optimize cryogenic hydrogen cooling within vortex tubes.

08 HYDROGEN↗

Separate sensible and latent cooling: Carnot limits and system taxonomy

Energy efficiency and thermal comfort can be improved by independent control of temperature and humidity. In recent work, an emerging class of cooling systems known as separate sensible and latent cooling (SSLC) systems have demonstrated 14 to 47% energy savings, while simultaneously providing better comfort control than conventional air conditioning. In this work, SSLC is defined as using two or more cooling processes with unequal sensible heat ratio (SHR) to control temperature and humidity independently and dynamically. Fundamental efficiency limits of SSLC systems are developed for the first time, demonstrating that the Carnot efficiency for SSLC is higher than that of conventional systems. It is revealed that SSLC has particularly high potential for moderate outdoor temperatures, high sensible building loads, and drier indoor conditions. Additionally, this work proposes a classification scheme for SSLC systems, revealing an enormous number of feasible cycle permutations. This provides a framework for the important research task of developing SSLC systems that will save energy, be practical to construct, and be controllable.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Single-Nucleon Transfer on Unstable 59 Cu Probes the NiCu Cycle in Astrophysical X-Ray Bursts

Recent models of the rapid proton (𝑟⁢𝑝) capture process indicate that a competition between the 59 Cu⁢(𝑝,𝛾) ⁢60 Zn and 59 Cu⁢(𝑝,𝛼) ⁢56 Ni reactions may result in the formation of a nickel-copper (NiCu) cycle that traps the flux of material between 56 Ni and 60 Zn . Here, in this work, we report the identification of 15 proton-unbound levels in 60 Zn , populated via 59 Cu⁢(𝑑,𝑛) transfer, which govern the rate of the 59 Cu⁢(𝑝,𝛾)⁢ 60 Zn reaction in XRBs. Precise excitation energies for levels in 60 Zn were obtained from observed 𝛾 decays, and spectroscopic factors were determined from angle-integrated cross sections. Incorporating these results into stellar-model calculations, we find that with experimentally constrained uncertainties a NiCu cycle in XRBs is indeed possible, though we limit its branching strength to less than 38%. While modest, such a branching has significant impact on the light curve, motivating further studies of the relevant rates. Our calculations also indicate that a significant NiCu cycle leads to an increase in the amount of odd-𝐴 nuclei in the burst ashes, which may affect Urca cooling processes in neutron star crusts.

direct reactions↗

Microwave synthesis of single-phase nanoparticles made of multi-principal element alloys

In this work, metal nanoparticles of multi-principal element alloys (MPEA) with a single crystalline phase have been synthesized by flash heating/cooling of nanosized metals encapsulated in micelle vesicles dispersed in an oil phase (e.g., cyclohexane). Flash heating is realized by selective absorption of a microwave pulse in metals to rapidly heat metals into uniform melts. The oil phase barely absorbs microwave and maintains the low temperature, which can rapidly quench the high-temperature metal melts to enable the flash cooling process. The precursor ions of four metals, including Au, Pt, Pd, and Cu, can be simultaneously reduced by hydrazine in the aqueous solution encapsulated in the micelle vesicles. The resulting metals efficiently absorb microwave energy to locally reach a temperature high enough to melt themselves into a uniform mixture. The duration of microwave pulse is crucial to ensure the reduced metals mix uniformly, while the temperature of oil phase is still low to rapidly quench the metals and freeze the single-phase crystalline lattices in alloy nanoparticles. The microwave-enabled flash heating/cooling provides a new method to synthesize single-phase MPEA nanoparticles of many metal combinations when the appropriate water-in-oil micelle systems and the appropriate reduction reactions of metal precursors are available.

36 MATERIALS SCIENCE↗

Water desalination through the dewpoint evaporative system

This paper presents the first numerical analysis of a novel Dew Point Desalination unit (DPD), which utilizes the phenomenon of the evaporative cooling process to minimize energy consumption. The proposed solution is based on a modular structure, which allows it to obtain any required water production capacity. The analyses for various operational parameters and various climatic conditions are based on the numerical simulations conducted using a validated mathematical model. The system performance was described with different factors, including specific energy consumption (SEC) and daily water production rate respected to 1 m^3*s^-1 of passing air stream (DWP). It is shown that the DPD system can operate in almost every climate zone, and the desalinated water production consumes less than 1.5 kWh*m^-3 of electric energy It is noted the best performance of the DPD system is achieved in the semi-arid, desert and subtropical climate zones. In these climate zones, the 2-stage DPD system allows the production of desalinated water with an electric energy consumption SEC lower than 0.5 kWh*m^-3. The additional benefit of the DPD system is its modularity. The presented DPD system stands as a module, that can be multiplied as needed to achieve the desired size of the desalination plant.

Desalination↗

A reduced-order modeling of a tubular solar reactor for long duration thermochemical energy storage

The storage of solar energy in a solid form, referred to as a “solar fuel”, can be achieved through a process known as endothermic solar thermochemistry. This process transforms the absorbed solar energy into a stable and retrievable form that can be stored for extended periods of time. This paper presents a low–order heat transfer model of a counter–current tubular falling bed reactor designed to produce thermally reduced magnesium manganese oxide pellets for long duration thermochemical energy storage. The energy required for the endothermic reduction was supplied by concentrated solar energy or renewable electricity via indirect heating of the gas and solid reactants flowing in a ceramic tube. The counter-current gas flow enhances the mixing of the solid particles with the heat recuperation zone, allowing the gas and particles to enter and exit the tubular reactor close to room temperature. Further, the reactor was vertically oriented and was heated circumferentially by an adjustable level heat flux along a finite segment of its length. The temperature distribution of the reactor in response to transient changes along the tube was modeled by considering conduction, convection, and radiation heat transfer. Governing equations for the heat transfer model were solved by discretizing the reactor tube into a finite number of control volumes and using an energy balance for the heat exchange between the reactor wall, gas, and particles within the control volume. The energy absorbed during this endothermic reaction was modeled numerically by fitting the data of the chemical conversion rate with the corresponding temperature of particles in the heating zone. The numerical model has been experimentally validated using a reactor prototype made of a 121.92 cm alumina tube heated by a 7kW electric tube–furnace. The alumina tube receives magnesium manganese oxide pellets of 3.66±0.516 mm in diameter from the top, and a counter–current gas flow from the bottom. The reactor wall temperature was monitored by six thermocouples installed along the reactor tube length. The experimental procedure was numerically simulated, and the temperature variation along the reactor tube was compared with a matrix of experimental runs for a range of particles mass flowrates (0.75–1.25g/s) and corresponding gas flowrates (36–65 SLPM). The reactor system was heated gradually from room temperature to a steady state temperature of 1673K, and then cooled down to room temperature. The heating and cooling processes were simulated, and the numerical and experimental results were compared throughout processes. The numerical model showed similar trends to the experimental results, with an error of 0.69 to 7.9% for the particle inlet and 0.7 to 7.9% for the gas inlet during steady-state operation. The proposed numerical model can be implemented as a simplified physical model to design a feedback control system to regulate reactor temperature.

14 SOLAR ENERGY↗

Thermal Unequilibrium of PdSn Intermetallic Nanocatalysts: From In Situ Tailored Synthesis to Unexpected Hydrogenation Selectivity

Abstract Effective control on chemoselectivity in the catalytic hydrogenation of C=O over C=C bonds is uncommon with Pd‐based catalysts because of the favored adsorption of C=C bonds on Pd surface. Here we report a unique orthorhombic PdSn intermetallic phase with unprecedented chemoselectivity toward C=O hydrogenation. We observed the formation and metastability of this PdSn phase in situ. During a natural cooling process, the PdSn nanoparticles readily revert to the favored Pd 3 Sn 2 phase. Instead, using a thermal quenching method, we prepared a pure‐phase PdSn nanocatalyst. PdSn shows an >96 % selectivity toward hydrogenating C=O bonds of various α,β‐unsaturated aldehydes, highest in reported Pd‐based catalysts. Further study suggests that efficient quenching prevents the reversion from PdSn‐ to Pd 3 Sn 2 ‐structured surface, the key to the desired catalytic performance. Density functional theory calculations and analysis of reaction kinetics provide an explanation for the observed high selectivity.

Chen, Minda↗

Thermal Unequilibrium of PdSn Intermetallic Nanocatalysts: From In Situ Tailored Synthesis to Unexpected Hydrogenation Selectivity

Effective control on chemoselectivity in the catalytic hydrogenation of C=O over C=C bonds is uncommon with Pd-based catalysts because of the favored adsorption of C=C bonds on Pd surface. Here we report a unique orthorhombic PdSn intermetallic phase with unprecedented chemoselectivity toward C=O hydrogenation. We observed the formation and metastability of this PdSn phase in situ. During a natural cooling process, the PdSn nanoparticles readily revert to the favored Pd 3 Sn 2 phase. Instead, using a thermal quenching method, we prepared a pure-phase PdSn nanocatalyst. Furthermore, PdSn shows an >96% selectivity toward hydrogenating C=O bonds of various α,β-unsaturated aldehydes, highest in reported Pd-based catalysts. Further study suggests that efficient quenching prevents the reversion from PdSn- to Pd 3 Sn 2 -structured surface, the key to the desired catalytic performance. Density functional theory calculations and analysis of reaction kinetics provide an explanation for the observed high selectivity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Crystal plasticity model of residual stress in additive manufacturing using the element elimination and reactivation method

Selective laser melting is receiving increasing interest as an additive manufacturing technique. Residual stresses induced by the large temperature gradients and inhomogeneous cooling process can favour the generation of cracks. In this work, a crystal plasticity finite element model is developed to simulate the formation of residual stresses and to understand the correlation between plastic deformation, grain orientation and residual stresses in the additive manufacturing process. The temperature profile and grain structure from thermal-fluid flow and grain growth simulations are implemented into the crystal plasticity model. An element elimination and reactivation method is proposed to model the melting and solidification and to reinitialize state variables, such as the plastic deformation, in the reactivated elements. The accuracy of this method is judged against previous method based on the stiffness degradation of liquid regions by comparing the plastic deformation as a function of time induced by thermal stresses. The method is used to investigate residual stresses parallel and perpendicular to the laser scan direction, and the correlation with the maximum Schmid factor of the grains along those directions. The magnitude of the residual stress can be predicted as a function of the depth, grain orientation and position with respect to the molten pool. The simulation results are directly comparable to X-ray diffraction experiments and stress–strain curves.

316 stainless steel↗

Air quality and comfort constrained energy efficient operation of multi-zone buildings

Maintaining indoor air quality (IAQ) through effective ventilation is essential for the well-being and productivity of building occupants. Control strategies aimed at improving the efficiency of heating, ventilation and air conditioning (HVAC) systems must jointly determine ventilation and heating and cooling processes. Here, in this paper, we study the problem of minimizing the energy consumption of the HVAC system in a multi-zone building, while meeting thermal comfort and IAQ requirements. We first perform a steady state analysis of the zonal carbon dioxide (CO 2 ) concentration and the temperature dynamics. The resulting expressions are convex in the zonal mass flow rates and zonal temperatures. Guided by the steady state solutions for meeting the thermal comfort constraints, we develop two control policies for improving the energy efficiency of building HVAC systems while jointly satisfying indoor temperature and IAQ constraints. We compare the performance of our proposed approaches with those of multiple baseline approaches which implement separate regimes for controlling zonal temperature and IAQ for a typical work-day in a multi-zone campus building. We have evaluated the performance of our proposed approaches under varying levels of flexibility in zonal temperatures. We have shown that zonal temperature flexibility can result in energy savings up to 32% (for the same control strategies) as compared to the case where no such flexibility is permitted. Our proposed approaches were seen to offer potential savings of nearly 29% compared to the baseline.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Development of experimental and computational frameworks to predict subcooled flow boiling in the LANL Isotope Production Facility

Cooling is crucial to maintain the integrity of target systems in isotope production facilities. At Los Alamos National Laboratory (LANL)’s Isotope Production Facility (IPF), multiple encapsulated targets are stacked and irradiated in tandem with a 100 MeV, ~250μA proton beam. To facilitate effective heat removal, these stacked targets are separated and cooled via a series of water channels. At these beam currents, this high-energy proton beam heats the target system, likely initiating subcooled flow boiling in the cooling channels. However, in-beam monitoring of the IPF target system is not possible due to the extreme radiation environment, and the necessarily significant shielding. To better understand high-power target performance, we developed ex-situ experimental and computational frameworks to predict the behavior of subcooled flow boiling at IPF. Subcooled flow boiling experiments on Inconel 625 samples under IPF conditions (2 bar pressure, 10 GPM flow rate (i.e., 2249 kg/m 2 /s), 85 K subcooling) revealed that IPF's average operating power is at the early stage of boiling with a heat transfer coefficient of 48,000 W/m 2 /s. The proposed modeling framework enables us to predict a complete boiling curve, i.e., single-phase heat transfer, onset of nucleate boiling, two-phase heat transfer, and critical heat flux (CHF), with specification of input boiling parameters up to intermediate heat flux levels. The estimated CHF under IPF conditions is 5.2 MW/m 2 . Experimental data under reduced conditions (2 bar pressure, 1.5 GPM flow rate (i.e., 337 kg/m 2 /s), 45 K subcooling) served as validation cases for the computational modeling. This computational model can be further extended to more complicated systems replicating the real IPF configuration, for instance, to study void distribution as a function of the incident proton beam profile and coolant velocity profile of multiple cooling channels. Finally, the proposed experimental and computational frameworks provide a means to better understand cooling systems in the isotope production facilities at different accelerators, where in-beam monitoring of the cooling process is not available.

07 ISOTOPE AND RADIATION SOURCES↗

Spark plasma sintering (SPS) densified U 3 Si 2 pellets: Microstructure control and enhanced mechanical and oxidation properties

Dense U 3 Si 2 pellets with controlled grain structure and enhanced thermal-mechanical and oxidation properties are synthesized with spark plasma sintering (SPS). Microstructure and phase composition of the SPS densified pellets are characterized systematically using SEM, EDS, and XRD. Thermal-mechanical properties and oxidation behavior of the sintered silicide fuel pellets are analyzed by laser flash, indentation, and dynamic thermogravimetric analysis. Dense U 3 Si 2 pellets are consolidated by combining high energy ball milling and rapid sintering by SPS, and the microstructure structures are controlled from micron-sized (~5.7 μm grain size) for conventional silicide to a nanocrystalline matrix with an average grain size of ~280 nm. A dominant phase of distorted U 3 Si 2 was identified with lattice expansion due to residual thermal stress upon SPS consolidation and rapid cooling processes. Both micron-sized and nano-sized pellets show exceptional thermal transport properties, consistent with monolithic silicides reported in literature. The SPS-densified pellets possess simultaneously high hardness and fracture toughness. The SPS-densified silicide pellets also demonstrate exceptional oxidation performance with extended onset oxidation temperature above 500 °C and reduced oxidation kinetics, particularly for nano-sized pellets. A strong strain effect was proposed in which compressive stress in nano-sized pellets enhances the oxidation resistance of silicide fuels, as evidenced by the degradation of oxidation performance upon strain relaxation by isothermal annealing. The correlation among the sintering process – microstructure control – physical properties and fuel behavior is established. A new concept of strain engineering is proposed further properties optimization, enabling the development of potential oxidation and corrosion-resistant silicides with extended performance, the key technological challenge of U 3 Si 2 as the leading concept of accident tolerant fuels.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Depicting the roles of CuO secondary phase and heat treatment in driving the magnetic and magnetocaloric features of Pr 2/3 Sr 1/3 MnO 3 manganite

In this work, we report a detailed experimental study regarding the impact of copper oxide (CuO) secondary phase and heat treatments on the structural, magnetic, and magnetocaloric properties of the near-room temperature Pr 2/3 Sr 1/3 MnO 3 (PSMO) magnetic refrigerant. Our investigations are carried out by using structural and microstructural analyses, alongside magnetization measurements. Here, the analysis of X-ray diffraction data of PSMO(95%)-CuO(5%) (PSMO-CuO) samples shows the coexistence of both CuO and PSMO phases. Further, the microstructural analysis of PSMO-CuO reveals that the addition of CuO significantly enhances the grains size. On the other hand, the added secondary phase markedly reduces the Curie temperature (T c ) from 293 K for PSMO to about 273 K for PSMO-CuO composite while increasing the magnetocaloric effect. This decrease in T c is associated with a significant change from 162° to about 156° in the Mn-OII-Mn bond angle respectively. Moreover, performed investigations regarding the role of heat treatments unveil that the observed changes in the structural and magnetic features are mainly driven by the secondary phase that modify grains size and double-exchange interactions in the PSMO compound. Interestingly, our findings demonstrate that the Curie temperature of the PSMO and accordingly its magnetocaloric effect can be tailored by adding small amounts of CuO without need to substitution on cation sites. In the light of obtained results, a multilayered refrigerant composed of PSMO and PSMO-CuO is proposed to cover the magnetic cooling temperature range close to room-temperature. The resulting entropy change remains practically constant between 273 K and 293 K. Such a behavior is highly appreciated from a practical point of view, particularly in cases where the cooling process is carried out by using the AMR and Ericsson cycles.

36 MATERIALS SCIENCE↗