Search NASA⌕ Search

SEARCH · Search NASA

Results for “fault localization”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Shallow Aseismic Slip in the Delaware Basin Determined by Sentinel-1 InSAR

The Delaware Basin, Texas is currently a hot-spot of induced seismicity and ground deformation due to fluid extraction and injection associated with horizontal drilling techniques; however, the driving mechanism behind the seismicity and deformation remains under debate. Here, using vertical and east-west horizontal surface deformation measurements derived from Sentinel-1 interferometric synthetic aperture radar (InSAR), we show that the subsurface responds differently to oil and gas activity in the northern and southeastern portions of the basin. In the north, where there is little seismicity, deformation patterns display long-wavelengths and equidimensional patterns. In contrast, the southeast region hosts most of the seismicity and displays spatial deformation patterns with narrow linear features that strike parallel to the maximum principal horizontal stress and to trends in seismicity, suggesting movement along normal faults. We model a linear deformation feature using edge dislocations and show that the InSAR observations can be reproduced by slip on normal faults contained within the Delaware Mountain Group (DMG), the formation that hosts local wastewater injection and the majority of earthquakes. Our model consists of three parallel, high-angle normal faults, with two dipping toward one another in a graben structure. Slip magnitudes reach up to 25 cm and are spatially correlated with injection wells. Measured seismicity can only explain ~2% of the fault motion predicted by our fault model, suggesting that slip leading to the deformation is predominantly aseismic. We conclude that seismic and aseismic fault motion in the southeastern Delaware Basin is likely driven by wastewater injection near critically-stressed normal faults within the DMG.

58 GEOSCIENCES↗
DOE OSTI2025

EV-EVSE Fault Study: An Analysis of Thermal Events Caused by Electrical Faults during DC Charging

This report outlines multiple avenues of analysis of electrical faults associated with electric vehicles (EVs) during charging, focusing specifically on the interactions between EVs and EV supply equipment (EVSE). Key concerns include the identification and mitigation of overtemperature events that can result in fires, which are commonly initiated by localized heating of connectors, wiring, or high-impedance fault current paths.

42 - ENGINEERING↗
DOE OSTI2020

Internal structure of the San Jacinto fault zone at the Ramona Reservation, north of Anza, California, from dense array seismic data

SUMMARY We image the internal structure of the San Jacinto fault zone (SJFZ) near Anza, California, with seismic data recorded by two dense arrays (RA and RR) from ∼42 000 local and ∼180 teleseismic events occurring between 2012 and 2017. The RA linear array has short aperture (∼470 m long with 12 strong motion sensors) and recorded for the entire analysed time window, whereas the RR is a large three-component nodal array (97 geophones across a ∼2.4 km × 1.4 km area) that operated for about a month in September–October 2016. The SJFZ at the site contains three near-parallel surface traces F1, F2 and F3 from SW to NE that have accommodated several Mw > 6 earthquakes in the past 15 000 yr. Waveform changes in the fault normal direction indicate structural discontinuities that are consistent with the three fault surface traces. Relative slowness from local events and delay time analysis of teleseismic arrivals in the fault normal direction suggest a slower SW side than the NE with a core damage zone between F1 and F2. This core damage zone causes ∼0.05 s delay at stations RR26–31 in the teleseismic P arrivals compared with the SW-most station, and generates both P- and S-type fault zone trapped waves. Inversion of S trapped waves indicates the core damaged structure is ∼100 m wide, ∼4 km deep with a Q value of ∼20 and 40 per cent S-wave velocity reduction compared with bounding rocks. Fault zone head waves observed at stations SW of F3 indicate a local bimaterial interface that separates the locally faster NE block from the broad damage zone in the SW at shallow depth and merges with a deep interface that separates the regionally faster NE block from rocks to the SW with slower velocities at greater depth. The multiscale structural components observed at the site are related to the geological and earthquake rupture history at the site, and provide important information on the preferred NW propagation of earthquake ruptures on the San Jacinto fault.

Geochemistry & Geophysics↗
DOE OSTI2023

Temperature Uncertainty Modeling with Proxy Structural Data as Geostatistical Constraints for Well Siting: An Example Applied to Granite Springs Valley, NV, USA

Utilizing existing temperature and structural information around Granite Springs Valley, Nevada, we build 3D stochastic temperature models with the aim of evaluating the 3D uncertainty of temperature and choosing between candidate exploration well locations . The data used to support the modeling are measured temperatures and structural proxies from 3D geologic modeling, the latter considered "secondary" data. Two stochastic geostatistical techniques are explored for incorporating the structural proxies: cosimulation and local varying mean. With both the cosimulation and local varying mean methods, many equally likely temperature models (i.e., realizations) are produced, from which temperature probability profiles are calculated at candidate well locations. To aid in choosing between the candidate locations, two quantities summarize the temperature probabilities: Vprior and entropy. Vprior quantifies the likelihood for economic temperatures at each candidate location, whereas entropy identifies where new information has the most potential to reduce uncertainty. In general, the cosimulation realizations have smoother spatial structure, and extrapolate high temperatures at candidate locations that are located along the direction of the longest spatial correlation, which are down dip from existing temperature logs. The smooth realizations result in tight temperature probability profiles that are easier to interpret, but they have unrealistic temperature reversals in some locations because the cosimulation technique does not enforce a conductive geothermal gradient as a baseline (i.e., linearly increasing temperature with depth). The local varying mean results produce realizations with more realistic geothermal gradients, with temperatures increasing downward since a depth-temperature relationship is included. However, because they have much noisier spatial nature compared to cosimulation, it is harder to interpret the temperature probability profiles. The different local varying mean results allow the geologist to determine which proxy (e.g., dilation versus distance to fault termination) should be used given the specific geothermal system. In general, Vprior from local varying mean results identify locations that are close to high values for the structural proxies: areas with highe r probabilities for higher temperatures. The entropy results identify where uncertainty is greatest and therefore new drilling information could be most useful. Though these techniques provide useful information, even when applied to areas of sparse data, our comp arison of these two techniques demonstrates the need for new geothermal geostatistics techniques that combine the advantages of these two methods and that are tailored to the spatial uncertainty issues inherent in geothermal exploration.

3D temperature modeling↗
DOE OSTI2020

Peak Ground Velocity Spatial Variability Revealed by Dense Seismic Array in Southern California

There are an infinite number of locations where seismic waves interact with Earth's surface, yet seismologists can only record a finite set of ground motions. This requires assumptions about how known measurements or predictions can be extrapolated to nearby locations. We test to what extent extrapolation is viable using an array of 1,108 geophones installed 10–20 m apart within a 0.36 km 2 study area for 1 month atop the San Jacinto Fault in Southern California. Measuring ground motions from 38 small local earthquakes, we observe peak ground velocity (PGV) variations up to 22% across the array. We find the basin structure consistently traps energy, amplifying the seismic velocities, and that the damaged region within the fault reduces the PGVs. Comparing surface and borehole measurements (148 m deep), we find PGVs recorded at the surface are 3–10 times greater than the measurements at depth. We attribute these differences to amplification from soft sediments within the basin. We show that by leveraging information from >1,000 measurements per earthquake the ground motion velocities can be tightly constrained. These results suggest the spatial extrapolation of PGVs in regions of known faults and basins, even across a small area, should be done with caution.

58 GEOSCIENCES↗
DOE OSTI2024

Enhancing detection of topological order by local error correction

The exploration of topologically-ordered states of matter is a long-standing goal at the interface of several subfields of the physical sciences. Such states feature intriguing physical properties such as long-range entanglement, emergent gauge fields and non-local correlations, and can aid in realization of scalable fault-tolerant quantum computation. However, these same features also make creation, detection, and characterization of topologically-ordered states particularly challenging. Motivated by recent experimental demonstrations, we introduce a paradigm for quantifying topological states—locally error-corrected decoration (LED)—by combining methods of error correction with ideas of renormalization-group flow. Our approach allows for efficient and robust identification of topological order, and is applicable in the presence of incoherent noise sources, making it particularly suitable for realistic experiments. We demonstrate the power of LED using numerical simulations of the toric code under a variety of perturbations. We subsequently apply it to an experimental realization, providing new insights into a quantum spin liquid created on a Rydberg-atom simulator. Finally, we extend LED to generic topological phases, including those with non-abelian order.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Facet-Dependent Cold Welding of Au Nanorods Revealed by Liquid Cell Transmission Electron Microscopy

Cold welding of metals at the nanoscale has been demonstrated to play a significant role in bottom-up manufacturing and self-healing processes of nanostructures and nanodevices. However, the welding mechanism at the nanoscale is not well understood. In this study, a comprehensive demonstration of the cold welding process of gold nanorods with different modes is presented through in situ liquid cell transmission electron microscopy. The experimental results and molecular dynamics simulations reveal that the nanorods are welded through the facet-dependent atomic surface diffusion and rearrangement along {100} facets. The density functional theory calculations indicate that the preferred coalescence of two {100} surfaces is thermodynamically favorable. Unlike the prevalent “oriented attachment” in the nanoparticle coalescence, the misalignment of nanorod orientations and local stresses can induce grain boundaries and stacking faults in the welded interface.

42 ENGINEERING↗
DOE OSTI2020

Highly deformable Mg–Al–Ca alloy with Al 2 Ca precipitates

Magnesium (Mg) is the lightest structural metal. However, the poor formability of Mg alloys to great extent limits their applications in making structural parts. Formability is strongly correlated to both high tensile elongation and large work hardening capacity. In this work, we report a new Mg-Al-Ca alloy in which a majority of deformable Al 2 Ca precipitates form while the formation of Laves phases of Mg 17 Al 12 and Mg 2 Ca seems suppressed. Al 2 Ca precipitates imped dislocation motion, leading to large work hardening. Then, Al 2 Ca precipitates deform with dislocations and stacking faults under the enhanced flow stress, which relieve local stress concentration and improve tensile elongation. In addition, solutes Al and Ca suppress twin nucleation while promoting non-basal dislocations in Mg. Our Mg-Al-Ca alloy demonstrates one of the highest combinations of tensile elongation and work hardening capacity among existing Mg alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Toward improved urban earthquake monitoring through deep-learning-based noise suppression

Earthquake monitoring in urban settings is essential but challenging, due to the strong anthropogenic noise inherent to urban seismic recordings. Here, we develop a deep-learning-based denoising algorithm, UrbanDenoiser, to filter out urban seismological noise. UrbanDenoiser strongly suppresses noise relative to the signals, because it was trained using waveform datasets containing rich noise sources from the urban Long Beach dense array and high signal-to-noise ratio (SNR) earthquake signals from the rural San Jacinto dense array. Application to the dense array data and an earthquake sequence in an urban area shows that UrbanDenoiser can increase signal quality and recover signals at an SNR level down to ~0 dB. Earthquake location using our denoised Long Beach data does not support the presence of mantle seismicity beneath Los Angeles but suggests a fault model featuring shallow creep, intermediate locking, and localized stress concentration at the base of the seismogenic zone.

58 GEOSCIENCES↗
DOE OSTI2023

Temperature uncertainty modelling with proxy structural data as geostatistical constraints for well siting: an example applied to Granite Springs Valley, NV, USA

Utilizing existing temperature and structural geology information around Granite Springs Valley, Nevada, we build 3D stochastic temperature models with the aims of evaluating the 3D uncertainty of temperature and choosing between candidate exploration well locations. The data used to support the modelling are measured temperatures and structural proxies from 3D geologic modelling (distance to fault, distance to fault intersections and terminations, Coulomb stress change and dilation tendency), the latter considered ‘secondary’ data. Two stochastic geostatistical techniques are explored for incorporating the structural proxies: cosimulation and local varying mean. With both the cosimulation and local varying mean methods, many equally-likely temperature models (i.e. realizations) are produced, from which temperature probability profiles are calculated at candidate well locations. To aid in choosing between the candidate locations, two quantities summarize the temperature probabilities: V prior and entropy. V prior quantifies the likelihood for economic temperatures at each candidate location, whereas entropy identifies where new information has the most potential to reduce uncertainty. In general, the cosimulation realizations have smoother spatial structure, and extrapolate high temperatures at candidate locations that are located along the direction of the longest spatial correlation, which are down dip from existing temperature logs. The smooth realizations result in tight temperature probability profiles that are easier to interpret, but they have unrealistic temperature reversals in some locations because of the dipping ellipsoid shape created and that the cosimulation technique does not enforce a conductive geothermal gradient as a baseline (i.e. linearly increasing temperature with depth). The local varying mean results produce realizations with more realistic geothermal gradients, with temperatures increasing downward since a depth-temperature relationship is included. However, because they have much noisier spatial nature compared to cosimulation, it is harder to interpret the temperature probability profiles. The different local varying mean results allow the geologist to determine which proxy (e.g. dilation v. distance to fault termination) should be used given the specific geothermal system. In general, V prior from local varying mean results identify locations that are close to high values for the structural proxies: areas with higher probabilities for higher temperatures. The entropy results identify where uncertainty is greatest and therefore new drilling information could be most useful. Though these techniques provide useful information, even when applied to areas of sparse data, our comparison of these two techniques demonstrates the need for new geothermal geostatistics techniques that combine the advantages of these two methods and that are tailored to the spatial uncertainty issues inherent in geothermal exploration.

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

Advanced High-Performance Computational Modeling of the Seismic Response of High-Hazard and/or Nuclear Facilities and Critical Infrastructure at the NNSS

New methods for predicting the amplitude and variability of ground shaking from earthquakes (and explosions) are needed for seismic hazard analysis for buildings, nuclear power plants, and critical infrastructure at the NNSS. We are comparing existing 1-D and new 3-D geophysical methods for estimating the shear-wave velocity structure in the upper 30 meters of the ground surface (Vs30), which plays a major role in ground motion amplification and seismic response of buildings. We evaluate the performance of these methodologies at the U1a Complex at the NNSS and develop simple 1-D and high-resolution 3-D Vs30 models. We then emplace these high-resolution models into a background seismic velocity model. We will collaborate with Lawrence Livermore National Laboratory (LLNL) to conduct numerical modeling of the ground shaking at the NNSS using their high-performance computing technology and state-of-the-art ground motion simulation methodology. The primary work that was completed in FY 2019 was to acquire the seismic systems and familiarize staff at the NNSS with their use. We also worked on developing a collection plan with the Device Assembly Facility (DAF) at the NNSS, but due to time constraints and other ongoing projects at the DAF, we had to use U1a as a backup. We were able to coordinate the seismic survey, and we will complete the collection of seismic data in FY 2020. Additionally, during FY 2019, we completed the geologic framework model (GFM) for the U1a Complex and modeled the Yucca fault. LLNL worked with us through FY 2019 to prepare the data files for their modeling software and tested the software for reliability. In FY 2020 we will develop the end-to-end capability so that any facility could easily be modeled and the expected shaking from a local earthquake understood. The work in FY 2020 will include building fault models from the GFM and finalizing the velocity model analysis. The final simulations will be run for multiple rupture models, and final assessments will demonstrate the seismic hazard at the U1a Complex.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Dislocation nano-hydrides in nickel: Nucleation, evolution and effects on dislocation behaviors

Nano-hydrides have been predicted to precipitate at the core of edge dislocations in the Ni-H system, a mechanism that may promote hydrogen embrittlement. However, nano-hydride nucleation, growth, and effects on dislocation behavior have seldom been explored. This work combines molecular dynamics grand canonical Monte Carlo (MD-GCMC) simulations and continuum modeling to uncover a wide range of phenomena linked to dislocation nano-hydrides. Simulations reveal that nano-hydrides can be stabilized at dislocation cores with all character angles, including screw segments, due to the hydrostatic stresses around the cores of the Shockley partials. Nano-hydride nucleation takes place in these regions, and growth is dictated by the character angle θ of the perfect dislocation. The equilibrium stacking fault width d eq varies dynamically to increase the local hydrostatic stress field and facilitate the formation of the nano-hydride, forming a constriction-like feature and leading to three distinct behaviors: d eq decreases for θ > 30°, d eq remains unchanged for θ = 30°, and d eq increases for θ < 30°. Remote hydrostatic and Escaig stresses are also shown to influence the nucleation stage, implying stress concentrations such as those ahead of crack tips may facilitate nano-hydride precipitation. Moreover, we identify a new hydrogen-induced 60° dislocation reaction that emits a Shockley partial on a conjugate plane, with potential implications for twin nucleation. Testable predictions from this study are then used to reinterpret previous results from the literature. In conclusion, these findings provide a comprehensive framework to assess nano-hydride formation and evolution at dislocations in nickel and other face-centered cubic metals, with important implications to hydrogen embrittlement.

Dislocations↗
DOE OSTI2024

Dynamic Adaptive Relaying for Distribution Grids with High Inverter-Based Generation

The penetration of distributed energy resources (DERs) in distribution systems is rapidly increasing. This can lead to large variations in the operating load and the fault current seen by a relay, which can impact the reliability of overcurrent protection. Adaptive overcurrent relaying (AOCR) can address this concern by enabling the relay to adjust its pickup settings based on the system operating conditions. This paper presents a novel dynamic AOCR approach that allows the relays to locally estimate the pickup settings and dynamically modify them. The approach is based on a two-stage estimation algorithm that uses local measurements and system information to estimate the worst-case fault current and the relay pickup settings. The relay pickup settings are then dynamically modified based on the estimated fault current. The performance of the proposed approach is evaluated on multiple configurations of the EPRI J1 feeder and the IEEE 13-bus system. The results show that the proposed approach can effectively improve the reliability of overcurrent protection in distribution systems with DERs.

adaptive overcurrent protection↗
DOE OSTI2024

Intermittent Criticality Multi‐Scale Processes Leading to Large Slip Events on Rough Laboratory Faults

Abstract We discuss data of three laboratory stick‐slip experiments on Westerly Granite samples performed at elevated confining pressure and constant displacement rate on rough fracture surfaces. The experiments produced complex slip patterns including fast and slow ruptures with large and small fault slips, as well as failure events on the fault surface producing acoustic emission bursts without externally‐detectable stress drop. Preparatory processes leading to large slips were tracked with an ensemble of ten seismo‐mechanical and statistical parameters characterizing local and global damage and stress evolution, localization and clustering processes, as well as event interactions. We decompose complex spatio‐temporal trends in the lab‐quake characteristics and identify persistent effects of evolving fault roughness and damage at different length scales, and local stress evolution approaching large events. The observed trends highlight labquake localization processes on different spatial and temporal scales. The preparatory process of large slip events includes smaller events marked by confined bursts of acoustic emission activity that collectively prepare the fault surface for a system‐wide failure by conditioning the large‐scale stress field. Our results are consistent overall with an evolving process of intermittent criticality leading to large failure events, and may contribute to improved forecasting of large natural earthquakes.

Geochemistry & Geophysics↗
DOE OSTI2025

Local chemical ordering of a neutron-irradiated CrFeMnNi compositionally complex alloy

While ion-irradiation studies are a critical first step in studying compositionally complex alloys (CCAs) for nuclear applications, they do not capture all the microstructural changes occurring under the low irradiation dose rates and different particles’ scattering patterns experienced in a nuclear reactor setting. To explore these phenomena in reactor-relevant conditions for the first time in CCA, the single-phase solid-solution Cr 10 Fe 30 Mn 30 Ni 30 was neutron irradiated up to 6.61 displacements per atom at 395 and 579 °C. Irradiation-enhanced local chemical ordering (LCO) well beyond the range of short range ordering was observed, and is predicted to be the precursor to the precipitation of a coherent Ni-Mn L1 0 phase and a Cr-rich α’ phase, though TEM analysis did not indicate the presence of either in any irradiation condition. The line density of faulted dislocation loops decreased from 6.47 to 1.69 ∙ 10 15 m -2 from 3.43 to 6.61 dpa at 579 °C despite no appreciable faulted loop content in the unirradiated material. LCO is expected to increase the complexity of the energy landscape within this alloy, restricting interstitial point defect mobility and creating local regions of greater stacking fault energy. These contribute to the negative correlation between irradiation dose and faulted dislocation loop density in this study, as well as the lack of void swelling observed.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Asymmetric equilibrium core structures of pyramidal-II < c + a > dislocations in ten hexagonal-close-packed metals

The structures of pyramidal-II < c + a > dislocations, one of the most important defects in structural hexagonal-close-packed (HCP) metals, have not been fully characterized for many of the HCP metals in use today. Here, we employ ab initio informed phase-field dislocation dynamics to determine the minimum energy structure of pyramidal {1¯1¯22} < 11¯23 > dislocations in ten HCP metals, including Be, Co, Mg, Re, Ti, Zn, Cd, Hf, Y, and Zr. As input for the simulations, we calculate, using first-principles density functional theory, the {1¯1¯22} generalized stacking fault energy (GSFE) curves for all ten metals. From these calculations, it is found that magnetism in Co is necessary for achieving a local minimum in the GSFE curve. We observe in simulations that edge and screw character dislocations split into two partials separated by a low-energy intrinsic stacking fault. The splitting distance is shown to scale inversely with the local minimum energy normalized by the product of its shear modulus and Burgers vector. Interestingly, some HCP metals exhibit an asymmetric structure, with either unequal partial Burgers vectors or widths, in contrast to the symmetric configuration expected from linear elastic dislocation theory. We explain these structures by properties of the local maxima in their GSFE curves. Metals with larger degrees of elastic anisotropy result in dislocations with larger splitting distances than would be expected under the commonly used assumption of elastic isotropy. Furthermore, these findings on the sizes and asymmetry in the structures of pyramidal-II < c + a > dislocations are fundamental to understanding how these dislocations glide and interact or react with other defects when these metals are mechanically strained.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Machine-Learning-Based High-Resolution Earthquake Catalog Reveals How Complex Fault Structures Were Activated during the 2016–2017 Central Italy Sequence

The 2016–2017 central Italy seismic sequence occurred on an 80 km long normal-fault system. The sequence initiated with the Mw 6.0 Amatrice event on 24 August 2016, followed by the Mw 5.9 Visso event on 26 October and the Mw 6.5 Norcia event on 30 October. We analyze continuous data from a dense network of 139 seismic stations to build a high-precision catalog of ~900,000 earthquakes spanning a 1 yr period, based on arrival times derived using a deep-neural-network-based picker. Our catalog contains an order of magnitude more events than the catalog routinely produced by the local earthquake monitoring agency. Aftershock activity reveals the geometry of complex fault structures activated during the earthquake sequence and provides additional insights into the potential factors controlling the development of the largest events. Activated fault structures in the northern and southern regions appear complementary to faults activated during the 1997 Colfiorito and 2009 L’Aquila sequences, suggesting that earthquake triggering primarily occurs on critically stressed faults. Delineated major fault zones are relatively thick compared to estimated earthquake location uncertainties, and a large number of kilometer-long faults and diffuse seismicity were activated during the sequence. These properties might be related to fault age, roughness, and the complexity of inherited structures. The rich details resolvable in this catalog will facilitate continued investigation of this energetic and well-recorded earthquake sequence.

58 GEOSCIENCES↗
DOE OSTI2024

Chiral Kirigami for Bend‐Tolerant Reconfigurable Hologram with Continuously Variable Chirality Measures

Abstract Despite the commonality of static holograms, the holography with multiple information layers and reconfigurable grey‐scale images at communication frequencies remain a confluence of scientific challenges. One well‐known difficulty is the simultaneous modulation of phase and amplitude of electromagnetic wavefronts with a high modulation depth. A less appreciated challenge is scrambling of the information and images with hologram bending. Here, this work shows that chirality‐guided pixelation of plasmonic kirigami sheets enables tunable multiplexed holography at terahertz (THz) frequencies. The convex and concave structures with slanted Au strips exhibit gradual variations in geometries facilitating modulation of light ellipticity reaching 40 deg. Real‐time switching of 3D images of the letter “ M ” and the Mona Lisa demonstrates the possibility of complex grey‐scale information content and importance of continuously variable mirror asymmetry. Microscale chirality measures of each pixel experiences little change with bending while retaining controllable reconfigurability upon stretching, which translates to remarkable resilience of chiral holograms to bending. Simplicity of their design with local chirality measures opens the door to information technologies with fault‐tolerant THz encryption, wearable holographic devices, and new communication technologies.

36 MATERIALS SCIENCE↗