Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “vortices”

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 307 records · Page 17

Multi-year assessment of the impact of ship-borne radiosonde observations on polar WRF forecasts in the Arctic

Abstract To compensate for the lack of conventional observations over the Arctic Ocean, ship-borne radiosonde observations have been regularly carried out during summer Arctic expeditions and the observed data have been broadcast via the global telecommunication system since 2017. With these data obtained over the data-sparse Arctic Ocean, observing system experiments were carried out using a polar-optimized version of the Weather Research and Forecasting (WRF) model and the WRF Data Assimilation (WRFDA) system to investigate their effects on analyses and forecasts over the Arctic. The results of verification against reanalysis data reveal: (1) DA effects on analyses and forecasts; (2) the reason for the year-to-year variability of DA effects; and (3) the possible role of upper-level potential vorticity in delayed DA effects. The overall assimilation effects of the extra data on the analyses and forecasts over the Arctic are positive. Initially, the DA effects are the most apparent in the temperature variables in the middle/lower troposphere, which spread to the wind variables in the upper troposphere. The effects decrease with time but reappear after approximately 120 h, even in the 240-h forecasts. The effects on forecasts vary depending on the proximity of the radiosonde observation locations to the high synoptic variability. The upper-level potential vorticity is known to play an important role in the development of Arctic cyclones, and it is suggested as a possible explanation for the delayed DA effects after about 120 h.

Geology↗

TEAMER: Electrically Engaged Undulation (EEL) System

The Electrically Engaged UnduLation (EEL) system is a buoyancy-driven submersible device for powering oceanographic instruments. Physically, EEL is a slender body whose flexible spine is made up of energy units interconnected by uniaxial hinges. Each unit consists of a pair of piezoelectric elements that converts the bending stress into electrical current to a battery charging circuit. An outer plastic skin forms a seal against water and allows for flexibility at hinge locations. At the top is a bluff body with electronics that holds a ballast for buoyancy adjustment. The bluff body is also responsible for creating fluid instabilities in its wake. When gliding through the water (mode 2), the spine will flex in response to the alternating vortices that shed from the head. This "lock-in" phenomenon occurs when the frequency at which vortices shed resonates with the EEL natural frequency, during which the efficient gaits were found in species of sea snake, eels, and fish. For active propulsion, a single motor can be placed at the first segment and provide the oscillatory input for propulsion similar to a dolphin's kick. Such anguilliform swimming is both efficient and nearly silent compared to a spinning propeller. Ultimately, mimicking bio-locomotion provides a viable path to a drag-reduced, self-propelled energy harvesting system for ocean monitoring.

16 TIDAL AND WAVE POWER↗

On the Effect of Beating during Nonlinear Frequency Chirping

Spectroscopic analyses of energetic particle (EP) driven bursts of MHD fluctuations in magnetically confined plasmas often exhibit chirps that occur simultaneously in groups of two or more. While the superposition of oscillations at multiple frequencies necessarily causes beating in the signal acquired by a localized external probe, self-consistent hybrid simulations of chirping EP modes in a JT-60U tokamak plasma have demonstrated the possibility of global beating, where the mode's electromagnetic field vanishes globally between beats and reappears with opposite phase [Bierwage et al., Nucl. Fusion 57, 016036 (2017)]. This implies that there can be a single coherent field mode that oscillates at multiple frequencies simultaneously when it is resonantly driven by multiple density waves in EP phase space. Conversely, this means that the EP density waves are mutually coupled and interfere with each other via the jointly driven field, a mechanism ignored in some theories of chirping. In this thesis-style treatise, we study the role of field pulsations in general and beating in particular using the Hamiltonian guiding center orbit-following code ORBIT with a reduced wave-particle interaction model in realistic geometry. Beating is found to drive the evolution of EP phase space structures. Additionally, a key mechanism is the pulsation of effective phase space islands combined with the alternation of their effective O- and X-points due to phase jumps between each beat. Observations: (1) Beating causes density wave fronts to advance radially in a pulsed manner and the resulting chirps become staircase-like. (2) The pulsations facilitate convective transfer of material between neighboring layers of phase space density waves. On the one hand, this may inhibit the early detachment of solitary phase space vortices. On the other hand, it facilitates the accumulation of hole and clump fragments into larger structures. (3) Long-range chirping is observed when massive holes or clumps detach and drift away from the turbulent belt around the seed resonance. It is remarkable that the detached vortices remain robust and,on average, maintain their concentric nested layers while being visibly perturbed by the field's continued beating.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Initial Atmospheric Transport of Particles

A new version of the Initial Atmospheric Transport (IAT) model has been developed. This report fully characterizes the new model. IAT is designed to model heat release into the atmosphere that contains contaminated particles. The heat release forms a buoyant puff or plume. Buoyancy is quickly dissipated in 10's to 100's of seconds through drag and entrainment of ambient air. The final location of contaminant particles in the atmosphere after dissipation of abnormal heat is an important input to long term transportation codes such as the National Oceanic and Atmospheric Administration's (NOAA) Hybrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT) code. Ultimately the settling of contaminants after long term transport provides a basis for estimating health consequences. In comparison to the old version, the new model includes explicit simulation of particles, newly formulated equations of motion for a rising puff, and a new plume model. The plume model releases many instances of a single puff along a single trajectory in a way that balances energy released over a longer period than would be appropriate for a single puff. The puff model has additions for vorticity, turbulence, virtual mass, drag, and entrainment. The particle model leads to direct interfacing between IAT and the long-term transport code HYSPLIT. Previously, particle locations were given a pattern based on IAT's puff trajectory with no particle simulation. Now, particles are explicitly simulated, and 3-D spatial distributions of particles are output based on final positions. Ballistic particles are thrown outward at high velocity from explosive events and tend to hit the ground before the end of the simulation. Non-ballistic particles tend to get pulled into the rising puff s ring of vorticity with some escaping to ambient conditions as buoyant energy dissipates. A new model verification and validation effort was made for the new plume model. The Wallops Island Solid Rocket Propellant (SRP) fire plume was extracted to a 3-D transient dataset. After calibration, the model performs better than previous single puff attempts. Even so, error levels in the process used to estimate location and ambient weather conditions approach 30%. This leads to the conclusion that unplanned photogrammetry evaluations of the IAT model are useful but not sufficiently accurate to quantify model parameters. Like previous attempts, the vertical motion entrainment parameters come out much too high in comparison to values for other studies in the literature. It is suspected that a new model of shear entrainment would alleviate this problem. Unfortunately, IAT already has more calibration parameters than the data comparisons thus far can resolve with statistical certainty. This underscores the need to verify and validate the IAT model using detailed CFD studies that allow ideal statistical comparison with the IAT model for hundreds of cases. Even so, validation by comparison to carefully planned experiments (i.e. with acceptable error levels) is also necessary because the dynamics of entrainment are at a sub-grid level for CFD. Though better validation and verification are most important, many new enhancements are still envisioned for future IAT development. These include particle rain out due to condensation and agglomeration of water, direct connection to the North American Mesoscale (NAM) weather dataset, dynamic weather conditions, including ambient particle entrainment into the puff, and detrainment at the puff to ambient mixing interface. In general, including detailed aerosol dynamics is the greatest coming challenge which will require an entire new sub-model in IAT. This is very important if contaminant particles change in such a way that they increase or decrease in their potential harm to the environment.

54 ENVIRONMENTAL SCIENCES↗

Non-Equilibrium relaxation, aging scaling, and critical depinning dynamics of skyrmions in disordered magnetic films. Final Report

Technical applications of type-II superconductors, especially high-Tc compounds, in external magnetic fields require an effective flux pinning mechanism to reduce Ohmic losses due to flux creep and flow. Naturally abundant point defects, dislocations, and twin boundaries as well as artificially introduced disorder may serve as effective pinning centers for magnetic vortices. The properties of interacting flux lines subject to strong thermal fluctuations and point-like or extended pinning centers have therefore been a major research focus in condensed matter physics, and considerable progress has been achieved in the characterization and mathematical description of the equilibrium phases and non-equilibrium steady states of vortex matter. Detailed studies of the rich phase diagrams that reveal novel glassy states as well as the associated complex dynamical phenomena have in turn enriched the theoretical understanding of condensed matter systems, elucidating specifically the kinetics of elastic objects driven through disordered media and the important role of topological defects such as vortices and dislocations in various materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Using Taylor Green Vortex to Assess Modeling Turbulence Transition

The Taylor Green Vortex problem is an interesting problem to study the transition from laminar flow to fully developed turbulence. The problem consists initially of large-scale vortices that interact and eventually form smaller vortices that create a cascade of energy from larger scales to smaller scales. This document describes a) when fully developed turbulence occurs, and b) the scale separation in the inertial range. In addition, the resolution needed to capture perturbations in transitional flow is also proposed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Taylor-Couette Flow Modeling (Rev.1)

The Taylor-Couette Flow is a form of viscous flow in which fluid flows between two cylindrical walls with a prescribed rotation. A diagram is shown on the right (Figure 1) where Ω is the rotational velocity, H is the distance between the two cylinders, R is the radius the outer cylinder, r is the radius of the inner cylinder and L is the length of the cylinders. When the rotating inner cylinder reaches a critical velocity, the flow of the liquid becomes unstable and creates vortices in between the layers. In this experiment, the Taylor-Couette cylinders were modeled so that the vortices would be created after reaching a critical velocity. After the initial experiment, the velocity of the inner cylinder was decreased and increased to look at how the results changed.

42 ENGINEERING↗

Atmospheric teleconnection processes linking winter air stagnation and haze extremes in China with regional Arctic sea ice decline

Abstract. Recent studies suggested significant impacts of boreal cryosphere changes on wintertime air stagnation and haze pollution extremes in China. However, the underlying mechanisms of such a teleconnection relationship remains unclear. Here we use the Whole Atmosphere Community Climate Model (WACCM) to investigate dynamic processes leading to atmospheric circulation and air stagnation responses to Arctic sea ice changes. We conduct four climate sensitivity experiments by perturbing sea ice concentrations (SIC) and corresponding sea surface temperature (SST) in autumn and early winter over the whole Arctic and three subregions in the climate model. The results indicate distinct responses in circulation patterns and regional ventilation to the region-specific Arctic changes, with the largest increase of both the probability (by 132 %) and the intensity (by 30 %) of monthly air stagnation extremes being found in the experiment driven by SIC and SST changes over the Pacific sector of the Arctic (the East Siberian and Chukchi seas). The increased air stagnation extremes are mainly driven by an amplified planetary-scale atmospheric teleconnection pattern that resembles the negative phase of the Eurasian (EU) pattern. Dynamical diagnostics suggest that convergence of transient eddy forcing in the vicinity of Scandinavia in winter is largely responsible for the amplification of the teleconnection pattern. Transient eddy vorticity fluxes dominate the transient eddy forcing and produce a barotropic anticyclonic anomaly near Scandinavia and wave train propagation across Eurasia to the downstream regions in East Asia. The piecewise potential vorticity inversion analysis reveals that this long-range atmospheric teleconnection of Arctic origin takes place primarily via the middle and upper troposphere. The anomalous ridge over East Asia in the middle and upper troposphere worsens regional ventilation conditions by weakening monsoon northwesterlies and enhancing temperature inversions near the surface, leading to more and stronger air stagnation and pollution extremes over eastern China in winter. Ensemble projections based on state-of-the-art climate models in the Coupled Model Intercomparison Project Phase 6 (CMIP6) corroborate this teleconnection relationship between high-latitude environmental changes and midlatitude weather extremes, though the tendency and magnitude vary considerably among each participating model.

54 ENVIRONMENTAL SCIENCES↗

Control-Oriented Model for Secondary Effects of Wake Steering

This paper presents a model to incorporate the secondary effects of wake steering in large arrays of turbines. Previous models have focused on the aerodynamic interaction of wake steering between two turbines. The model proposed in this paper builds on these models to include yaw-induced wake recovery and secondary steering seen in large arrays of turbines when wake steering is performed. Turbines operating in yaw-misaligned conditions generate counter-rotating vortices that entrain momentum and contribute to the deformation and deflection of the wake at downstream turbines. Rows of turbines can compound the effects of wake steering that benefit turbines far downstream. This model quantifies these effects and demonstrates that wake steering has greater potential to increase the performance of a wind farm due to these counter-rotating vortices especially for large rows of turbines. This is validated using numerous large-eddy simulations for three-turbine, five-turbine, and wind farm scenarios.

17 WIND ENERGY↗

Measurements of Magnetic Field Penetration of Materials for Superconducting Radiofrequency Cavities

Superconducting Radio Frequency (SRF) cavities used in particle accelerators are typically formed from or coated with superconducting materials. Currently high purity niobium is the material of choice for SRF cavities which have been optimized to operate near their theoretical field limits. This brings about the need for significant R&D efforts to develop next generation superconducting materials which could outperform Nb and keep up with the demands of new accelerator facilities. To achieve high quality factors and accelerating gradients, the cavity material should be able to remain in the superconducting Meissner state under high RF magnetic field without penetration of quantized magnetic vortices through the cavity wall. Therefore, the magnetic field at which vortices penetrate in a supercon- ductor is one of the key parameters of merit of SRF cavities. Techniques to measure the onset of magnetic field penetration on thin film samples need to be developed to mitigate the issues with the conventional magnetometry measurements which are strongly influenced by the film orientation and shape and edge effects. The applied magnetic field is also needed to be parallel to the one side of the superconductor to resemble the magnetic field profile at the surface of the SRF cavities operating at fundamental accelerating mode. In this work we report the development of an experimental setup called Magnetic Field Penetration (MFP) magnetometer to measure the field of full flux penetration through bulk, thin films and multilayered superconductors. Our system combines a small superconducting solenoid which can generate the magnetic field up to 500 mT at the sample surface and three Hall probes to detect the full flux penetration through the planner superconductor with 2 inch diameter. This setup was used to study alternative materials which could potentially outperform niobium, as well as SIS multilayer coatings on niobium.

Senevirathne, Iresha↗

Numerical Analysis of Novel Plate Type Heat Exchanger with Oval-Twisted Channels

The novel heat exchanger (HX) designs implementing geometrical and surface modification combinations are expected to perform better than traditional HX technologies. Applied enhancement techniques seek to achieve (1) higher overall heat transfer performance, (2) increased compactness, and (3) simplified or comparable manufacturability. One innovative enhancement technique is to generate swirling flow vortices induced by channel or tube twisting. The turbulence generated by the twisted cross-section greatly enhances the heat transfer rate with minimal increases in pressure drop. Plate-type HXs and Printed Circuit Heat Exchangers (PCHXs) are compact designs that achieve high heat transfer rates per unit volume by utilizing several small channels, which maximizes the heat transfer surface area between the hot and cold fluids. Currently, advanced manufacturing technologies enable the design and fabrication of compact-type units with complicated channel geometries to achieve the highest performance and meet the compactness criteria of innovative HX technology. The proposed HX design concept combines the compactness of plate-type HXs and twisted channels, which provide additional turbulence and flow swirl enhancement. The plate-type oval-twisted HX (PTOTHX) is a crossflow configuration, with 16 short channels on one side (for hot fluid) and 8 long channels on the perpendicular side (for cold fluid). The inlet plenums have flow guide vanes to redirect flow and produce uniformity across the various flow paths. The compact size and purportedly improved heat transfer performance of the PTOTHX investigated herein prove its viability in various applications. Some notable potential nuclear applications of the PTOTHX include reactor core, spent fuel cooling, and residual heat dissipation. To establish a reference case, circular channels (PTCHX) are also considered in the present study for comparison with (PTOTHX). This paper aims to outline the numerical analysis procedures for determining the viability of the PTOTHX by comparing its heat transfer performance with the PTCHX units. The computational study used STAR-CCM+, a commercial computational fluid dynamics (CFD) code. Sensitivity analysis and model selection studies are conducted to determine the appropriate mesh density and turbulence model to provide the reported results. Numerical analyses comparing the Nusselt number (Nu) of the PTOTHX design with a comparable HX unit, including a circular cross-section and no twisting (PTCHX), show an overall heat transfer performance increase of 29-55% for balanced flow and 29-59% for imbalanced flow. Oval-cross-sectional twisted channels induce swirling flow vortices, enhancing the working fluid's convective heat transfer capabilities.

25 ENERGY STORAGE↗

Braiding of Majorana zero modes in vortex cores

Here, we demonstrate the successful simulation of $\sqrt{𝑍-}$, $\sqrt{𝑋-}$, and 𝑋-quantum gates using Majorana zero modes (MZMs) that emerge in magnetic vortices located in topological superconductors. We compute the transition probabilities and geometric phase differences accounting for the full many-body dynamics and show that qubit states can be read out by fusing the vortex core MZMs and measuring the resulting charge density. We visualize the gate processes using the time- and energy-dependent nonequilibrium local density of states. Our results demonstrate the feasibility of employing vortex core MZMs for the realization of fault-tolerant topological quantum computing.

Majorana bound states↗

Critique on “A Revisit of Navier-Stokes Equation”

In a recently published work (Sheng, 2020), Sheng identified inconsistencies in the concepts and principles leading to the Navier–Stokes equations and proposed a new asymmetric “friction tensor” to resolve these inconsistencies. We show that this friction tensor can be interpreted as the non-isotropic part of a vorticity-influenced viscous stress tensor following the analysis in Berdahl and Strang (1986). The asymmetric stress tensor however violates the principles of mechanical equilibrium in fluid flows that are free of external or body couples and is therefore non-physical in classical continuum hydrodynamics. We refute the arguments in favour of the asymmetric stress tensor in Sheng (2020) whilst also demonstrating that the friction tensor leads to the incorrect form of the Navier–Stokes equations in non-Cartesian co-ordinates and its ability to recover the correct governing equations in Cartesian co-ordinates is fortuitous.

asymmetric↗

Numerical simulation of evaporating wavy falling liquid films in laminar gas streams

Numerical simulations are performed to investigate the interfacial heat and mass transfer from evaporating wavy falling liquid films in interaction with laminar gas streams. The OpenFOAM solver has been used to conduct the simulations where the liquid-gas interface is resolved using the Volume of Fluid method of solving for three phases, i.e., liquid, vapor, and air. The configuration considered is a falling liquid (water) film on a heated vertical plate with a confined laminar moist air (gas) flow that is either (a) co-current or (b) counter-current to the downward liquid flow. The evaporation at the liquid-gas interface is driven by the interfacial gradient of the vapor mass fraction. Interfacial waves are triggered using a monochromatic forcing disturbance that leads to sinusoidal or solitary waves forming at the liquid-gas interface under respective forcing frequencies. Further, the numerical model is validated with the available experimental data. The results show nearly a 15% enhancement in time-averaged Sherwood number (Sh) due to film waviness (sinusoidal or solitary) at the lower volumetric gas flow rate, Q g = +50 (co-current) and Q g = -50 (counter-current). This enhancement in the Sh for both the waves further increases by 11% with Q g = +800 and 196% with Q g = -800. A closer examination of the mass transfer process over a wave demonstrates that with Q g = +50, the concentration of the gas side streamlines at the trough locations of the wave leads to higher values of Sh at these locations. However, with Q g = +800, although the overall Sh increases, vortices appear at the wave trough locations, leading to a corresponding decrease in the local Sh values. Correlations are proposed for predicting Sh under co-current and counter-current gas flow effects.

36 MATERIALS SCIENCE↗

Hyperon global polarization in isobar Ru+Ru and Zr+Zr collisions at $\sqrt{s_{NN}}$ = 200Ge V

The polarization of Λ, $\overline{Λ}$, $Ξ$ – , and $\overline{Ξ}$ + hyperons along the angular momentum of the system has been measured in isobar collisions of Ru+Ru and Zr+Zr at $\sqrt{s_{NN}}$ = 200 GeV with the STAR detector at RHIC. The polarization dependence on collision centrality exhibits an increasing trend in more peripheral collisions. Λ and $\overline{Λ}$ polarization dependence on the transverse momentum and pseudorapidity have been investigated, but no significant dependence was observed. The polarizations of Λ and $\overline{Λ}$ are found to be consistent with each other, indicating little contribution of the spin-magnetic coupling to the measured polarization. Comparison to previously measured polarization in Au+Au collisions show no obvious system size dependence. The results are qualitatively consistent with hydrodynamic calculations including contributions from shear-induced polarization and thermal vorticity. For the first time in heavy-ion collisions, the dependence of the global polarization on the hyperon’s emission azimuthal angle relative to the second-order event plane has been measured, indicating stronger polarization for the in-plane emitted hyperons at the level of 2.4 σ significance in 20–50 % centrality. The $Ξ$ hyperon polarization measurements via polarization transfer analysis yield finite positive values with 2.9 σ significance in 20–50 % centrality, slightly larger compared to the inclusive Λ polarization.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Pacific decadal oscillation causes fewer near-equatorial cyclones in the North Indian Ocean

Tropical cyclones do not form easily near the equator but can intensify rapidly, leaving little time for preparation. We investigate the number of near-equatorial (originating between 5°N and 11°N) tropical cyclones over the north Indian Ocean during post-monsoon season (October to December) over the past 60 years. The study reveals a marked 43% decline in the number of such cyclones in recent decades (1981–2010) compared to earlier (1951–1980). Here, we show this decline in tropical cyclone frequency is primarily due to the weakened low-level vorticity modulated by the Pacific Decadal Oscillation (PDO) and increased vertical wind shear. In the presence of low-latitude basin-wide warming and a favorable phase of the PDO, both the intensity and frequency of such cyclones are expected to increase. Such dramatic and unique changes in tropical cyclonic activity due to the interplay between natural variability and climate change call for appropriate planning and mitigation strategies.

54 ENVIRONMENTAL SCIENCES↗

Vortex phases and glassy dynamics in the highly anisotropic superconductor HgBa 2 CuO 4+δ

We present an extensive study of vortex dynamics in a high-quality single crystal of HgBa 2 CuO 4+ δ , a highly anisotropic superconductor that is a model system for studying the effects of anisotropy. From magnetization M measurements over a wide range of temperatures T and fields H , we construct a detailed vortex phase diagram. We find that the temperature-dependent vortex penetration field H p ( T ), second magnetization peak H smp ( T ), and irreversibility field H irr ( T ) all decay exponentially at low temperatures and exhibit an abrupt change in behavior at high temperatures T / T c >~0.5. By measuring the rates of thermally activated vortex motion (creep) S ( T , H ) = | d ln M ( T , H )/ d ln t |, we reveal glassy behavior involving collective creep of bundles of 2D pancake vortices as well as temperature- and time-tuned crossovers from elastic (collective) dynamics to plastic flow. Based on the creep results, we show that the second magnetization peak coincides with the elastic-to-plastic crossover at low T , yet the mechanism changes at higher temperatures.

36 MATERIALS SCIENCE↗

Designing high-performance superconductors with nanoparticle inclusions: Comparisons to strong pinning theory

One of the most promising routes for achieving high critical currents in superconductors is to incorporate dispersed, non-superconducting nanoparticles to control the dissipative motion of vortices. However, these inclusions reduce the overall superconducting volume and can strain the interlaying superconducting matrix, which can detrimentally reduce T c . Consequently, an optimal balance must be achieved between the nanoparticle density n p and size d. Determining this balance requires garnering a better understanding of vortex–nanoparticle interactions, described by strong pinning theory. Here, we map the dependence of the critical current on nanoparticle size and density in (Y 0.77 , Gd 0.23 )Ba 2 Cu 3 O 7-δ films in magnetic fields of up to 35 T and compare the trends to recent results from time-dependent Ginzburg–Landau simulations. We identify consistency between the field-dependent critical current J c (B) and expectations from strong pinning theory. Specifically, we find that J c ∝ B -α , where α decreases from 0.66 to 0.2 with increasing density of nanoparticles and increases roughly linearly with nanoparticle size d/ξ (normalized to the coherence length). At high fields, the critical current decays faster (~B -1 ), suggesting that each nanoparticle has captured a vortex. When nanoparticles capture more than one vortex, a small, high-field peak is expected in J c (B). Due to a spread in defect sizes, this novel peak effect remains unresolved here. Finally, we reveal that the dependence of the vortex creep rate S on nanoparticle size and density roughly mirrors that of α, and we compare our results to low-T nonlinearities in S(T) that are predicted by strong pinning theory.

36 MATERIALS SCIENCE↗