Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Phase slips”

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 253 records · Page 14

The performance of suppressed carrier receivers in a pulsed RFI environment

This paper investigates the tracking performance of a Costas loop in the presence of single and multiple in-band pulsed RFI signals and receiver additive noise. It is shown that the loop's tracking performance is very sensitive to the pulse durations of the RFI signals, as well as their repetition rates and peak powers relative to those of the desired signal. The presence of in-band RFI can affect the loop's tracking performance in basically three ways. 1) The loop tracks the composite phase of the RFI instead of the phase of the desired signal when the effective RFI power is greater than that of the signal. 2) When the effective RFI power is not sufficiently large to capture the loop, there still will exist static phase error biases which cause the loop to lock at a phase between that of the desired signal and the composite RFI. 3) There is an increase in the rms phase jitter over the nominal design. These latter two effects degrade the bit error probability performance of the receiver. More significantly, if the nominal loop bandwidth is not designed properly to account for the increase in jitter phenomenon, frequency cycle slipping may result.

Simon, M. K.↗

Progress of High Efficiency Centrifugal Compressor Simulations Using TURBO

Three-dimensional, time-accurate, and phase-lagged computational fluid dynamics (CFD) simulations of the High Efficiency Centrifugal Compressor (HECC) stage were generated using the TURBO solver. Changes to the TURBO Parallel Version 4 source code were made in order to properly model the no-slip boundary condition along the spinning hub region for centrifugal impellers. A startup procedure was developed to generate a converged flow field in TURBO. This procedure initialized computations on a coarsened mesh generated by the Turbomachinery Gridding System (TGS) and relied on a method of systematically increasing wheel speed and backpressure. Baseline design-speed TURBO results generally overpredicted total pressure ratio, adiabatic efficiency, and the choking flow rate of the HECC stage as compared with the design-intent CFD results of Code Leo. Including diffuser fillet geometry in the TURBO computation resulted in a 0.6 percent reduction in the choking flow rate and led to a better match with design-intent CFD. Diffuser fillets reduced annulus cross-sectional area but also reduced corner separation, and thus blockage, in the diffuser passage. It was found that the TURBO computations are somewhat insensitive to inlet total pressure changing from the TURBO default inlet pressure of 14.7 pounds per square inch (101.35 kilopascals) down to 11.0 pounds per square inch (75.83 kilopascals), the inlet pressure of the component test. Off-design tip clearance was modeled in TURBO in two computations: one in which the blade tip geometry was trimmed by 12 mils (0.3048 millimeters), and another in which the hub flow path was moved to reflect a 12-mil axial shift in the impeller hub, creating a step at the hub. The one-dimensional results of these two computations indicate non-negligible differences between the two modeling approaches.

turbomachinery↗

Evaluation of Nuclear Spent Fuel Disposal in Clay-Bearing Rock - Process Model Development and Experimental Studies

The DOE R&D program under the Spent Fuel Waste Science Technology (SFWST) campaign has made key progress in modeling and experimental approaches towards the characterization of chemical and physical phenomena that could impact the long-term safety assessment of heat-generating nuclear waste disposition in deep clay/shale/argillaceous rock. International collaboration activities such as heater tests and postmortem analysis of samples recovered from these have elucidated key information regarding changes in the engineered barrier system (EBS) material exposed to years of thermal loads. Chemical and structural analyses of sampled bentonite material from such tests has as well as experiments conducted on these are key to the characterization of thermal effects affecting bentonite clay barrier performance and the extent of sacrificial zones in the EBS during the thermal period. Thermal, hydrologic, and chemical data collected from heater tests and laboratory experiments has been used in the development, validation, and calibration of THMC simulators to model near-field coupled processes. This information leads to the development of simulation approaches (e.g., continuum vs. discrete) to tackle issues related to flow and transport at various scales of the host-rock and EBS design concept. Consideration of direct disposal of large capacity dual-purpose canisters (DPCs) as part of the back-end SNF waste disposition strategy has generated interest in improving our understanding of the effects of elevated temperatures on the EBS design. This is particularly important for backfilled repository concepts where temperature plays a key role in the EBS behavior and long-term performance. This report describes multiple R&D efforts on disposal in argillaceous geologic media through development and application of coupled THMC process models, experimental studies on clay/metal/cement barrier and host-rock (argillite) material interactions, molecular dynamic (MD) simulations of water transport during (swelling) clay dehydration, first-principles studies of metaschoepite (UO 2 corrosion product) stability, and advances in thermodynamic plus surface complexation database development. Drift-scale URL experiments provides key data for testing hydrological-chemical (HC) model involving strong couplings of fluid mixing and barrier material chemical interactions. The THM modeling focuses on heater test experiments in argillite rock and gas migration in bentonite as part of international collaboration activities at underground research laboratories (URLs). In addition, field testing at an URL involves in situ analysis of fault slip behavior and fault permeability. Pore-scale modeling of gas bubble migration is also being investigated within the gas migration modeling effort. Interaction experiments on bentonite samples from heater test under ambient and elevated temperatures permit the evaluation of ion exchange, phase stability, and mineral transformation changes that could impact clay swelling. Advances in the development, testing, and implementation of a spent nuclear fuel (SNF) degradation model coupled with canister corrosion focus on the effects of hydrogen gas generation and its integration with Geologic Disposal Safety Assessment (GDSA). GDSA integration activities includes evaluation of groundwater chemistries in shale formations.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Mechanisms for Microseismicity Occurrence Due to CO 2 Injection at Decatur, Illinois: A Coupled Multiphase Flow and Geomechanics Perspective

Here, we numerically investigate the mechanisms that resulted in induced seismicity occurrence associated with CO 2 injection at the Illinois Basin–Decatur Project (IBDP). We build a geologically consistent model that honors key stratigraphic horizons and 3D fault surfaces interpreted using surface seismic data and microseismicity locations. We populate our model with reservoir and geomechanical properties estimated using well-log and core data. We then performed coupled multiphase flow and geomechanics modeling to investigate the impact of CO 2 injection on fault stability using the Coulomb failure criteria. We calibrate our flow model using measured reservoir pressure during the CO 2 injection phase. Our model results show that pore-pressure diffusion along faults connecting the injection interval to the basement is essential to explain the destabilization of the regions where microseismicity occurred, and that poroelastic stresses alone would result in stabilization of those regions. Slip tendency analysis indicates that, due to their orientations with respect to the maximum horizontal stress direction, the faults where the microseismicity occurred were very close to failure prior to injection. These model results highlight the importance of accurate subsurface fault characterization for CO 2 sequestration operations.

58 GEOSCIENCES↗

Molecular Dynamics Simulation of a Multi-Walled Carbon Nanotube Based Gear

We used molecular dynamics to investigate the properties of a multi-walled carbon nanotube based gear. Previous work computationally suggested that molecular gears fashioned from (14,0) single-walled carbon nanotubes operate well at 50-100 gigahertz. The gears were formed from nanotubes with teeth added via a benzyne reaction known to occur with C60. A modified, parallelized version of Brenner's potential was used to model interatomic forces within each molecule. A Leonard-Jones 6-12 potential was used for forces between molecules. The gear in this study was based on the smallest multi-walled nanotube supported by some experimental evidence. Each gear was a (52,0) nanotube surrounding a (37,10) nanotube with approximate 20.4 and 16,8 A radii respectively. These sizes were chosen to be consistent with inter-tube spacing observed by and were slightly larger than graphite inter-layer spacings. The benzyne teeth were attached via 2+4 cycloaddition to exterior of the (52,0) tube. 2+4 bonds were used rather than the 2+2 bonds observed by Hoke since 2+4 bonds are preferred by naphthalene and quantum calculations by Jaffe suggest that 2+4 bonds are preferred on carbon nanotubes of sufficient diameter. One gear was 'powered' by forcing the atoms near the end of the outside buckytube to rotate to simulate a motor. A second gear was allowed to rotate by keeping the atoms near the end of its outside buckytube on a cylinder. The ends of both gears were constrained to stay in an approximately constant position relative to each other, simulating a casing, to insure that the gear teeth meshed. The stiff meshing aromatic gear teeth transferred angular momentum from the powered gear to the driven gear. The simulation was performed in a vacuum and with a software thermostat. Preliminary results suggest that the powered gear had trouble turning the driven gear without slip. The larger radius and greater mass of these gears relative to the (14,0) gears previously studied requires a smaller rotation rate and multiple rows of teeth to avoid excessive force on the gear teeth resulting, in slip and failure of the driven gear to turn. We hope that studies such as these will eventually lead to synthesis of components that can be assembled into atomically precise fullerene machines. These machines, in turn, may someday be used in machine-phase fullerene materials with remarkable properties.

Han, Jie↗

Characterization of local parameters of inclined upward air–water bubbly two-phase flows in a round pipe

The present work seeks to develop a detailed experimental database and characterize the local effects of inclination on horizontal and upward bubbly air–water two-phase flows. Experimental conditions are selected to investigate bubbly flows at a variety of orientations while reducing the influence of flow regime transition. Detailed time-averaged data is collected using state-of-the-art four-sensor conductivity probes for measurements across the flow area for five flow conditions at five inclination orientations in a 25.4 mm inner diameter pipe. Using the newly established database, effects of inclination on bubbly flows are analyzed, as are changes in liquid and gas flow rate at different orientations. The observations include (1) bubble migration and agglomeration at the top pipe wall as orientation changes from vertical to horizontal; (2) reducing the angle from vertical reduces the local gas velocity, and correspondingly reduces the global slip from greater than unity at 90° to less than unity for 60° and lower; (3) changes in the bubble size distribution are primarily at the interface between the bubble cluster and the liquid region of the pipe, not coalescence within the bubble cluster. In conclusion, physical mechanisms for these observations are proposed, and the present trends at different orientations support previous experiments in horizontal two-phase flows.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Fluctuations and Pinning for Individually Manipulated Skyrmions

We numerically examine the dynamics of individually dragged skyrmions interacting simultaneously with an array of other skyrmions and quenched disorder. For drives just above depinning, we observe a broadband noise signal with a 1/ f characteristic, while at higher drives, narrowband or white noise appears. Even in the absence of quenched disorder, the threshold force that must be applied to translate the driven skyrmion is finite due to elastic interactions with other skyrmions. The depinning threshold increases as the strength of the quenched disorder is raised. Above the depinning force, the skyrmion moves faster in the presence of quenched disorder than in a disorder-free system since the pinning sites prevent other skyrmions from being dragged along with the driven skyrmion. For strong pinning, we find a stick-slip motion of the driven skyrmion which produces a telegraph noise signature. The depinning threshold increases monotonically with skyrmion density in the absence of quenched disorder, but when pinning is present, the depinning threshold changes nonmonotonically with skyrmion density, and there are reentrant pinned phases due to a competition between pinning induced by the quenched disorder and that produced by the elastic interactions of the skyrmion lattice.

36 MATERIALS SCIENCE↗

Tectonics of the Outer Planet Satellites

Tectonic features on the satellites of the outer planets range from the familiar, such as clearly recognizable graben on many satellites, to the bizarre, such as the ubiquitous double ridges on Europa, the twisting sets of ridges on Triton, or the isolated giant mountains rising from Io's surface. All of the large and middle-sized outer planet satellites except Io are dominated by water ice near their surfaces. Though ice is a brittle material at the cold temperatures found in the outer solar system, the amount of energy it takes to bring it close to its melting point is lower than for a rocky body. Therefore, some unique features of icy satellite tectonics may be influenced by a near-surface ductile layer beneath the brittle surface material, and several of the icy satellites may possess subsurface oceans. Sources of stress to drive tectonism are commonly dominated by the tides that deform these satellites as they orbit their primary giant planets. On several satellites, the observed tectonic features may be the result of changes in their tidal figures, or motions of their solid surfaces with respect to their tidal figures. Other driving mechanisms for tectonics include volume changes due to ice or water phase changes in the interior, thermoelastic stress, deformation of the surface above rising diapirs of warm ice, and motion of subsurface material toward large impact basins as they fill in and relax. Most satellites exhibit evidence for extensional deformation, and some exhibit strike-slip faulting, whereas contractional tectonism appears to be rare. Io s surface is unique, exhibiting huge isolated mountains that may be blocks of crust tilting and foundering into the rapidly emptying interior as the surface is constantly buried by deposits from hyperactive volcanoes. Of the satellites, diminutive Enceladus is spectacularly active; its south polar terrain is a site of young tectonism, copious heat flow, and tall plumes venting into space. Europa's surface is pervasively tectonized, covered with a diverse array of exotic and incompletely understood tectonic features. The paucity of impact craters on Europa suggests that its tectonic activity is ongoing. Geysers on Triton show that some degree of current activity, while tectonic features that cross sparsely cratered terrain indicate that it may also be tectonically active. Ganymede and Miranda both exhibit ancient terrains that have been pulled apart by normal faulting. On Ganymede these faults form a global network, while they are confined to regional provinces on Miranda. Ariel, Dione, Tethys, Rhea, and Titania all have systems of normal faults cutting across their surfaces, though the rifting is less pronounced than it is on Ganymede and Miranda. Iapetus exhibits a giant equatorial ridge that has defied simple explanation. The rest of the large and middle-sized satellites show very little evidence for tectonic features on their surfaces, though the exploration of Titan's surface has just begun.

McKinnon, W. B.↗

Fifth-degree elastic energy for predictive continuum stress–strain relations and elastic instabilities under large strain and complex loading in silicon

Materials under complex loading develop large strains and often phase transformation via an elastic instability, as observed in both simple and complex systems. Here, we represent a material (exemplified for Si I) under large Lagrangian strains within a continuum description by a 5th-order elastic energy found by minimizing error relative to density functional theory (DFT) results. The Cauchy stress—Lagrangian strain curves for arbitrary complex loadings are in excellent correspondence with DFT results, including the elastic instability driving the Si I → II phase transformation (PT) and the shear instabilities. PT conditions for Si I → II under action of cubic axial stresses are linear in Cauchy stresses in agreement with DFT predictions. Such continuum elastic energy permits study of elastic instabilities and orientational dependence leading to different PTs, slip, twinning, or fracture, providing a fundamental basis for continuum physics simulations of crystal behavior under extreme loading.

36 MATERIALS SCIENCE↗

Phase field modeling of dislocations and obstacles in InSb

We present a phase-field dislocation dynamics (PFDD) model informed by first-principle calculations to elucidate the competitive dislocation nucleation and propagation between the glide and shuffle sets in InSb diamond cubic crystal. The calculations are directly informed with generalized stacking fault energy curves on the (111) slip plane for both the “glide set,” with the smaller interplanar spacing, and the “shuffle set,” with the larger interplanar spacing. The formulation also includes elastic anisotropy and the gradient term associated with the dislocation core. The PFDD calculations show that under no stress the equilibrium structure of screw glide set dislocations dissociates into Shockley partials, while those of the shuffle set dislocations do not dissociate, remaining compact. The calculated dislocation core widths of these InSb dislocations agree well with the measured values for other semiconductor materials, such as Si and GaN. We find that a shuffle set dislocation emits from a dislocation source at an applied stress about three times smaller than that needed to emit leading and trailing partials successively on the glide set plane. Once the partial dislocations in the glide set are emitted, they propagate faster than the shuffle set perfect dislocation at the same stress level.

36 MATERIALS SCIENCE↗

Fifth-degree elastic energy for predictive continuum stress-strain relations and elastic instabilities under large strain and complex loading in silicon

Materials under complex loading develop large strains and often phase transformation via an elastic instability, as observed in both simple and complex systems. Here, we represent a material (exemplified for Si I) under large Lagrangian strains within a continuum description by a 5th-order elastic energy found by minimizing error relative to density functional theory (DFT) results. The Cauchy stress—Lagrangian strain curves for arbitrary complex loadings are in excellent correspondence with DFT results, including the elastic instability driving the Si I → II phase transformation (PT) and the shear instabilities. PT conditions for Si I → II under action of cubic axial stresses are linear in Cauchy stresses in agreement with DFT predictions. Such continuum elastic energy permits study of elastic instabilities and orientational dependence leading to different PTs, slip, twinning, or fracture, providing a fundamental basis for continuum physics simulations of crystal behavior under extreme loading.

Levitas, Valery↗

Enhanced work hardening from oxygen-stabilized ω precipitates in an aged metastable β Ti-Nb alloy

High levels of oxygen in solid solution in Ti alloys are considered detrimental to mechanical properties because of embrittlement concerns. In metastable β titanium alloys, the formation of isothermal ω precipitates is also known to cause severe embrittlement and ductility reduction. However, oxygen has been shown to partition to the ω phase during ageing, and this partitioning behavior may potentially impact ω’s mechanical contribution. Using micropillar compression, we compared the deformation behavior of Ti-20Nb (at. %) with oxygen-stabilized ω precipitates to the behavior of oxygen-free specimens. The oxygen-stabilized microstructures showed increased compressive yield strength and enhanced work hardening behavior compared to oxygen-free specimens. In the absence of oxygen, the compressed pillars showed slip band formation and catastrophic failure, and transmission electron microscopy imaging revealed that ω precipitates were sheared within the continuous deformation channels resulting in slip localization. In contrast, oxygen-stabilized ω precipitates were harder to shear and the formation of continuous deformation channels was suppressed during compression, leading to improved work hardening behavior up to 15% strain. Furthermore, this counter-intuitive role of oxygen may offer design strategies to address the significant embrittlement and loss of ductility observed for ω-strengthened β Ti alloys without oxygen and avenues to expand the use of β Ti alloys.

36 MATERIALS SCIENCE↗

Mechanochemical α to β phase transition of U 3 O 8

Recognized as the most stable phase of uranium oxide, α-U 3 O 8 is widely used throughout the nuclear fuel cycle for safe storage and transportation. There is, however, a lack of understanding of the thermodynamic relationship between α-U 3 O 8 and its polymorph, β-U 3 O 8 . This research contributes new knowledge regarding the kinetics and mechanism of a previously described mechanochemical phase transition between the α- and β-U 3 O 8 polymorphs to improve the understanding of their thermodynamic relationship. In this work, tumble milling of α-U 3 O 8 using different milling media shows an ingrowth of the β-phase over time, with an observed correlation between media density, percent conversion and lattice strain. Anisotropic peak broadening observed within collected X-ray diffraction powder patterns suggests the preservation of uranium polyhedral layered sheets throughout the milling process. Preservation of the sheets implies a shear-induced slip mechanism occurring through in-plane shifting of uranium polyhedra along lattice planes perpendicular to the polyhedral sheets.

Roach, Jordan M. [Oak Ridge National Laboratory (O↗

Journey to Time-Variable Moment Tensors through Inversion of Acoustic and Seismoacoustic Data

We explore the capability of acoustic and seismoacoustic datasets to directly resolve a complex, time-variable source consisting of a buried mechanism, represented as a moment tensor, and a spall mechanism, represented as a vertical force at the surface. Traditionally, each component of a resolved moment tensor assumes one underlying source time function, which likely fails to capture the full evolution of a dynamic source, such as an explosion followed by slip on near-source joints or development of spallation. Specifically, we expand previous work to resolve a time-variable moment tensor using single-modality and joint-modality inversion frameworks through analysis of infrasound and seismoacoustic data recorded as part of the Source Physics Experiment Phase II: Dry Alluvium Geology (DAG). We investigate the impact of including signals from seismic-to-air coupling that are local to each infrasound sensor in comparison to mainly atmosphere-propagating acoustic signals, which occur from coupling of the wavefield from the subsurface to the atmosphere directly above the source. Additionally, we assess the ability of our inversion algorithm to fit observed infrasound data using a variety of time-variable source mechanisms. First, we consider the buried moment tensor source alone, which assumes that the determined Green’s functions incorporate effects from spallation or that the impact from spallation is minimal. Second, we examine the estimated buried moment tensor and vertical surface spallation as terms that must both be resolved in the inversion. Third, we assess the ability for an estimated vertical surface spallation source to fit the acoustic data on its own. Finally, we compare results from the joint inversion of both seismic geophone and infrasound acoustic data for the buried-only source compared to buried and spallation sources. Our results are a preliminary investigation into the applications of the inversion technique to recorded datasets and show the technique has limited capabilities using acoustic data alone. Instead, this method shows promise for seismic and seismoacoustic datasets to resolve the time-variable mechanisms of a buried source.

47 OTHER INSTRUMENTATION↗

Carrier and symbol synchronization system performance study

Results pertinent to predicting the performance of convolutionally encoded binary phase-shift keyed communication links were presented. The details of the development are provided in four sections. These sections are concerned with developing the bit error probability performance degradations due to PN despreading by a time-shared delay locked loop, the Costas demodulation process, symbol synchronization effects and cycle slipping phenomena in the Costas loop. In addition, Costas cycle slipping probabilities are studied as functions of Doppler count time and signal-to-noise conditions. The effect of cycle slipping in the symbol synchronizer is also studied as a function of channel Doppler and other frequency uncertainties.

Lindsey, W. C.↗

The Plunge Phase of Friction Stir Welding

The many advantages of Friction Stir Welding have led to a relatively rapid acceptance in the often conservative welding community. Because the process is so different from traditional fusion welding, with which most investigators are most familiar, there remain many aspects of FSW for which there is no clear consensus. For example, the well known onion rings seen in transverse sections have been variously interpreted as grain size variations, variation in density of second phase particles and parts of the carousel of material rotating with the pin that have been shed from the carousel. Using Orientation Imaging Microscopy, Schneider has recently noted that the onion rings have a different orientation (and hence etch differently) than the surrounding material, and this orientation is consistent with slip plane orientations at the edge of the carousel. Likewise, the forces and torque exerted by the FSW tool on the work piece largely remain unaccounted for. Although these forces are routinely measured by investigators with commercial instrumented welders, they are rarely reported or even qualitatively analyzed. This paper will introduce a model based on a carousel or disk of material that rotates with the tool to estimate the torque and plunge force required to plunge a tool into the work piece. A stationary tool is modeled rather than the moving tool because effects such as thermal transients and metallurgical changes in the sample (primarily aging in aluminum) can be more easily accounted for. It is believed, however, that with some modifications the model should be applicable to a moving tool also.

McClure, John C.↗

Observations of the Polarized Solar Corona During the Annular Eclipse of 14 October 2023

We present results of a dual eclipse expedition to observe the solar corona from two sites during the annular solar eclipse of 14 October 2023 using a novel coronagraph designed to be accessible for amateurs and students to build and deploy. The coronagraph (CATEcor) builds on the standardized eclipse observing equipment developed for the Citizen CATE 2024 experiment. The observing sites were selected for likelihood of clear observations, for historic relevance (near the Climax site in the Colorado Rocky Mountains), and for centrality to the annular eclipse path (atop Sandia Peak above Albuquerque, New Mexico). The novel portion of CATEcor is an external occulter assembly that slips over the front of a conventional dioptric telescope, forming a shaded-truss externally occulted coronagraph. CATEcor is specifically designed to be easily constructed in a garage or “makerspace” environment. We successfully observed some bright features in the solar corona to an altitude of approximately 2.25 R during the annular phases of the eclipse. Future improvements to the design, in progress now, will reduce both stray light and image artifacts; our objective is to develop a design that can be operated successfully by amateur astronomers at sufficient altitude even without the darkened skies of a partial or annular eclipse.

Solar instrumentation↗

A spall and diffraction study of nanosecond pressure release across the iron ε-α phase boundary

The extreme response of polycrystalline iron at high pressures and high strain rates is revealed by means of high-power laser pulses. The compression portion of the pulse coupled with x-ray diffraction identifies the expected body-centered cubic (α) to hexagonal close packed (ε) displacive transformation. Upon release, observation shows that the complete reverse transformation takes approximately 8 ns and that the structure returns to its initial microstructural configuration, in a reversible transformation path. This is in good agreement with molecular dynamics (MD) simulations which predict an inverse dependence between transformation time and strain rate. The grain size is reduced from μm to nm range during compression and begins increasing back to the original grain size on decompression. The kinetics of the transition is dictated by heterogenous nucleation as it follows the Johnson-Mehl-Avrami-Kolmogorov equation with the appropriate time exponent of ~1. This is confirmed by MD simulations which also identify profuse twinning and dislocation generation. Here, the tensile pulse generated upon reflection at the free surface is captured by time-resolved free surface velocity measurements from which a peak tensile stress of 7 GPa is obtained, in stark contrast with its quasi-static value of ~200 MPa. At these strain rates, the strength of grain interiors, which is determined by twinning and slip exceeds the strength of the boundaries, and failure initiates preferentially in the latter.

36 MATERIALS SCIENCE↗