Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “acoustical waves”

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 613 records · Page 34

Development of a Laser Ultrasonics-based Approach for Rapid Screening of High Entropy Alloys

This project utilized a laser ultrasonic technique to systematically study temperature-induced evolution of material properties in a set of interrelated binary alloys and a high entropy alloy (HEA) fabricated using arc-melting and spark plasma sintering processes. This technique involved the use of a nanosecond duration, high-intensity pulsed laser to thermo-elastically generate ultrasonic waves that propagate in the bulk of the metal alloy. Sub-nanometer-scale displacements associated with the propagating bulk ultrasonic waves were detected along the epicenter on the opposite surface of the sample using a 1 GHz bandwidth photorefractive interferometer. Phase transformations and microstructural changes were inferred from the temperature-dependent trends in the bulk acoustic velocities and features in the ultrasonic epicentral waveforms measured in the binary alloys. These inferences were then correlated with electron/optical microscopy observations and predictions using calculations of phase diagrams (CALPHAD). The results showed that the laser-generated ultrasonic pulses were strongly influenced by changes in material microstructure and could accurately track thermally driven phase transformations and detect the presence of microscale heterogeneities (grain boundaries, dendritic structures, etc.) in the set of binary alloy samples. The measurement approach was then applied to a quinary HEA sample for estimating phase transition temperature and determining microstructural heterogeneity. The rapid, non-contact and non-destructive ultrasonic testing approach demonstrated here is amenable to high throughput combinatorial investigations that can be applied to graded composition HEAs produced using advanced manufacturing methods. When paired with atomistic simulations and CALPHAD modeling, this approach can overcome the bottlenecks faced by current material characterization methods in efficiently screening the vast discovery space of HEAs that spans over a hundred million unique quinary alloy compositions.

36 MATERIALS SCIENCE↗

A Resonant Metamaterial Line Array for Ultrasound Compressive Imaging

Abstract Acoustic metamaterials have been proposed for numerous applications including subwavelength imaging, impedance matching, and lensing. Yet, their application in compressive sensing and imaging has not been fully investigated. When metamaterials are used as resonators at certain frequencies, they can generate random radiation patterns in the transmitted waves from the transducers and received waves from a target. Compressive sensing favors such randomness inasmuch as it can increase incoherence by decreasing the amount of mutual information between any two different measurements. This study aims at assessing whether the use of resonating metamaterial unit cells in a single-layered non-optimized array between a number of ultrasound transceivers and targets can improve the sensing capacity, point-spread function of the sensing array (their beam focusing ability), and imaging performance in point-like target detection. The theoretical results are promising and can open the way for more efficient metamaterial designs with the aim of enhancing ultrasound imaging with lower number of transceivers compared to the regular systems.

Acoustics↗

A simple and effective five-equation two-phase numerical model for liquid-vapor phase transition in cavitating flows

Numerical difficulties, notably the non-monotonic behavior of the Wood speed of sound and the volume fraction positivity, associated with the reduced five-equation two-phase flow model of Kapila et al. (2001) [A.K. Kapila, R. Menikoff, J.B. Bdzil, S.F. Son, D.S. Stewart, 2001. Two-phase modeling of deflagration-to-detonation transition in granular materials: reduced equations, Physics of Fluids 13(10), 3002–3024] have been resolved in the past through the introduction of a frozen speed of sound and an algebraic approach for mechanical relaxation afforded by a pressure non-equilibrium six-equation model proposed by [R. Saurel, F. Petitpas, R.A. Berry, 2009. Simple and efficient relaxation methods for interfaces separating compressible fluids, cavitating flows and shocks in multiphase mixture, J. Comput. Phys. 228, 1678–1712]. By contrast, it is explored and demonstrated in this work that these difficulties can in fact still be resolved within the numerical scheme for solving the reduced five-equation model by numerically replacing the Wood speed of sound for the estimates of wave speeds in the approximate Riemann solver HLLC with the monotonic mixture speed of sound for a transport five-equation model. For shock interface (artificial mixture separating pure or nearly pure fluids) interaction problems, with the apparent advantage of monotonic behavior of the speed of sound in the interface, the effect of the numerical replacement is also confined to the interface. Differences other than the behavior of the speed of sound within the interface in the solutions due to the replacement diminish with increasing resolution when reasonable solution can be obtained with Wood speed of sound. For cavitating/expansion problems in physical fluid mixture, it is pointed out and explained why acoustics in the numerical solutions still propagate at the Wood speed of sound (therefore consistent with the reduced five-equation model) even though in some cases a much higher speed of sound like the numerical replacement above for solving the reduced five-equation model or the frozen speed of sound for solving a six-equation model is used for the estimates of wave speeds in the HLLC scheme. A variant of the five-equation two-phase flow model by Saurel et al. (2008) [R. Saurel, F. Petitpas, R. Abgrall, 2008. Modelling phase transition in metastable liquids: application to cavitating and flashing flows, J. Fluid Mech. 607, 313–350] is then constructed for liquid-vapor phase transition in cavitating flows. The relaxation toward thermo-chemical equilibrium during phase transition is achieved by solving a simple system of algebraic equations for the equilibrium state variables for better efficiency, following Pelanti and Shyue (2014) [M. Pelanti, K.-M. Shyue, 2014. A mixture-energy-consistent six-equation two-phase numerical model for fluid with interfaces, cavitation and evaporation waves. J. Comput. Phys. 259, 331–357]. Therefore, the current model retains both the simplicity afforded by the five-equation model and the efficiency of the algebraic relaxation solver. An alternative algebraic approach for handling the non-conservative term (the so-called K∇ · u term) in the reduced five-equation model for mechanical equilibrium of a liquid-vapor mixture is also explored by enforcing the thermal equilibrium at the same time. Finally, numerical results of sample tests in both one and two dimensions in the literature as well as that in three dimensions demonstrate the effectiveness and ability of the proposed model to simulate cavitating flows. An interesting mechanism of shock generation by acoustics in water due to phase transition is then found by the numerical simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Magnetoelastic effects in the hyperhoneycomb Kitaev spin liquid

In recent years, several magnetic Mott insulators with strong spin-orbit coupling were suggested to be proximate to the Kitaev quantum spin liquid, whose one of the most exciting features is the fractionalization of spin excitations into itinerant Majorana fermions and static Z 2 fluxes. Unfortunately, the ground states of these systems cannot be easily captured by experiment, remaining featureless to conventional local probes. Here, we propose to study the signatures of fractionalized excitations by exploiting their coupling to the lattice vibrations, dubbed magnetoelastic coupling, which arises from the fact that the interaction between spins depends on the relative distance between them. We argue that the magnetoelastic coupling can lead to the distinct modification of the phonon dynamics, which can be observed by measuring renormalized phonon spectrum, the sound attenuation, and the phonon Hall viscosity. This makes the phonon dynamics a promising tool for the characterization and identification of quantum spin liquid phases. In this work, we focus on the magnetoelastic effects in the three-dimensional Kitaev model realized on the hyperhoneycomb lattice. Here, the hyperhoneycomb Kitaev spin liquid is particularly interesting since the strong Kitaev interaction was observed in the Kitaev magnet β-Li 2 IrO 3 , for which the spin-orbit entangled J eff = 1/2 moments of iridium ions form precisely the hyperhoneycomb lattice.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

The stability of expanding reactive shocks in a van der Waals fluid

Despite the extensive literature accumulated since the pioneering works of D'yakov and Kontorovich in the 1950s, the stability of steady shocks is still an open question when realistic boundary conditions are accounted. The consideration of a supporting mechanism, which is indeed a necessary condition for shock steadiness, modifies the perturbation shock dynamics in the unstable range. The Noh problem is a suitable example to form steady expanding shocks. This configuration is of great interest to the high-energy-density-physics community because of its direct application to inertial confinement fusion and astrophysics, for which the stagnation of a supersonically converging material via an accretion shock front is ubiquitous. Here in this work, we extend the generalized Noh problem, both base-flow solution and linear stability analysis, to conditions where endothermic or exothermic transformations undergo across the shock. Within the spontaneous acoustic emission conditions found for a van der Waals gas, we find that cylindrical and spherical expanding shocks become literally unstable for sufficiently high mode numbers. Counterintuitively, the effect of exothermicity or endothermicity across the shock is found to be stabilizing or destabilizing, respectively.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Phonon softening and slowing-down of charge density wave fluctuations in BaNi 2 As 2

BaNi 2 As 2 is a nonmagnetic analog of the iron pnictide superconductors, and exhibits an incommensurate charge density wave (IC-CDW) and a sizable elastoresistance. Here, in this Letter, phonons in BaNi 2 As 2 associated with the IC-CDW and uniform in-plane lattice distortions are investigated using high-resolution inelastic x-ray scattering. The in-plane transverse acoustic phonons reveal no softening at temperatures where the elastoresistance increases strongly, indicating the latter to be electronically driven. Systematic phonon measurements suggest the IC-CDW occurs in two stages upon cooling: Underdamped phonons first soften to zero energy well above the IC-CDW ordering temperature, then the resulting quasielastic IC-CDW fluctuations gradually slow down and coalesce into the static IC-CDW order. A possible origin for our observations is the IC-CDW in BaNi 2 As 2 being uniaxial, which provides an additional Ising degree of freedom favorable for disordered IC-CDW modulations, and accounts for the elastoresistance through a weak coupling to the lattice.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Airborne Infrasound Makes a Splash

Natural and anthropogenic events may create low frequency sound waves, or infrasound, that can travel for vast distances in planetary atmospheres. They permit the remote monitoring of geophysical activity over local to global scales. Most studies have utilized ground-based recorders, but it is possible to deploy acoustic sensors to altitudes of over 50 km. Such elevated platforms can capture sounds that their surface analogs cannot access. High altitude balloons and low altitude aerostats are filling this observation gap, but key environments remain out of reach of both of these. Furthermore, recent work by den Ouden, Smets et al. (2021) addressed this with a new instrumentation platform—a large seabird flying just above the ocean's surface. Their work demonstrates that, infrasound sensing using heavier-than-air platforms in windy environments is possible, which has implications both terrestrially (e.g., extending sensor networks over the oceans) and extraterrestrially (proposed or planned missions to Venus and Titan).

47 OTHER INSTRUMENTATION↗

Magnetization dynamics fingerprints of an excitonic condensate t 2 g 4 magnet

The competition between spin-orbit coupling (SOC) λ and electron-electron interaction U leads to a plethora of novel states of matter, extensively studied in the context of t 2 g 4 and t 2 g 5 materials, such as ruthenates and iridates. Excitonic magnets—the antiferromagnetic state of bounded electron-hole pairs-are prominent examples of phenomena driven by those competing energy scales. Interestingly, recent theoretical studies predicted that excitonic magnets can be found in the ground state of SOC t 2 g 4 Hubbard models. Here we present a detailed computational study of the magnetic excitations in that excitonic magnet, employing one-dimensional chains (via density matrix renormalization group) and small two-dimensional clusters (via Lanczos). Specifically, first we show that the low-energy spectrum is dominated by a dispersive (acoustic) magnonic mode, with extra features arising from the λ = 0 state in the phase diagram. Second, and more importantly, we found a novel magnetic excitation forming a high-energy optical mode with the highest intensity at wave-vector q → 0 . In the excitonic condensation regime at large U , we also have found a novel high-energy π mode composed solely of orbital excitations. Overall, these features do not appear all together in any of the neighboring states in the phase diagram and thus constitute unique fingerprints of the t 2 g 4 excitonic magnet, of importance in the analysis of neutron and resonant inelastic x-ray scattering experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Humming Trains in Seismology: An Opportune Source for Probing the Shallow Crust

Abstract Seismologists are eagerly seeking new and preferably low-cost ways to map and track changes in the complex structure of the top few kilometers of the crust. By understanding it better, they can build on what is known regarding important, practical issues. These include telling us whether imminent earthquakes and volcanic eruptions are generating telltale underground signs of hazard, about mitigation of induced seismicity such as from deep injection of wastewater, how the Earth and its atmosphere couple, and where accessible natural resources are. Passive seismic imaging usually relies on blind correlations within extended recordings of Earth’s ceaseless “hum” or coda of well-mixed, small vibrations. In this article, we propose a complementary approach. It is seismic interferometry using opportune sources—specifically ones not stationary in time and moving in a well-understood configuration. Its interpretation relies on an accurate understanding of how these sources radiate seismic waves, precise timing, careful placement of pairs of listening stations, and seismic phase differentiation (surface and body waves). Massive freight trains were only recently recognized as such a persistent, powerful cultural (human activity-caused) seismic source. One train passage may generate a tremor with an energy output of a magnitude 1 earthquake and be detectable for up to 100 km from the track. We discuss the source mechanisms of train tremors and review the basic theory on sources. Finally, we present case studies of body- and surface-wave retrieval as an aid to mineral exploration in Canada and to monitoring of a southern California fault zone. We believe noise recovery from this new signal source, together with dense data acquisition technologies such as nodes or distributed acoustic sensing, will deeply transform our ability to monitor activity in the shallow crust at sharpened resolution in time and space.

Geochemistry & Geophysics↗

Ultrasonic microphone and ultrasonic acoustic radio

This disclosure provides systems, methods, and apparatus related to an ultrasonic microphone and an ultrasonic acoustic radio. In one aspect a system includes a transmitter and a receiver. The receiver comprises a membrane. The membrane comprises a single layer or multiple layers of a two-dimensional material. The receiver is operable to receive sound waves in a frequency range, with the frequency range being the ultrasonic frequency range.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Ultrasonic microphone and ultrasonic acoustic radio

This disclosure provides systems, methods, and apparatus related to an ultrasonic microphone and an ultrasonic acoustic radio. In one aspect a system includes a transmitter and a receiver. The receiver comprises a membrane. The membrane comprises a single layer or multiple layers of a two-dimensional material. The receiver is operable to receive sound waves in a frequency range, with the frequency range being the ultrasonic frequency range.

Zhou, Qin↗

Ultrasonic microphone and ultrasonic acoustic radio

This disclosure provides systems, methods, and apparatus related to an ultrasonic microphone and an ultrasonic acoustic radio. In one aspect a system includes a transmitter and a receiver. The receiver comprises a membrane. The membrane comprises a single layer or multiple layers of a two-dimensional material. The receiver is operable to receive sound waves in a frequency range, with the frequency range being the ultrasonic frequency range.

Zhou, Qin↗

UNH TDP - Concurrent Measurements of Inflow, Power Performance, and Loads for a Grid-Synchronized Vertical Axis Cross-Flow Turbine Operating in a Tidal Estuary

This data was collected between October 12 and December 15 of 2021 at the University of New Hampshire (UNH) and Atlantic Marine Energy Center (AMEC) turbine deployment platform (TDP). This data set includes over 29 days of grid connected turbine operation during this 65 day time frame. The priority for this measurement campaign was to collect data while the turbine was electrically connected to the grid by means of a rectifier and inverter. The Fall_2021_UNH_Measurement_Timeline.png highlights when each instrument was functioning and the Fall_2021_UNH_Test_Log.jpg indicates the four main regions for analysis available from this measurement campaign. The TDP is a floating structure moored on the Portsmouth facing side of Memorial Bridge pier #2, which spans the Piscataqua River between Portsmouth, NH and Kittery, ME. The Piscataqua River connects the Great Bay Estuary to the Gulf of Maine and the river currents are dominated by tidal forcing with water velocities exceeding 2.5 m/s during spring ebb tides at this site which were previously characterized by Kaelin Chancey (Assessment Of The Localized Flow And Tidal Energy Conversion System At An Estuarine Bridge - UNH MS Thesis 2019). The turbine under test was a modified New Energy Corporation (Calgary, CA) model EVG-025 4-blade H-Darrius type vertical axis cross flow turbine that rotates in the clockwise direction with a rotor diameter of 3.2m and blade length of 1.7m. The hydro-foil profile was a NACA 0021 with a 10 inch chord length and a blade preset pitch angle of +4deg with a positive angle corresponding with the toe in direction. The standard EVG-025 has a rotor diameter of 3.4m and its rated power output is 25kW at 3 m/s. The rotor diameter was reduced to accommodate the size of the existing TDP moon-pool. This project was pursued to quantify device performance for cross flow turbines operating in a marine environment. Accurate physical models, to characterize cross flow turbine performance, require real operational data sets due to the complexity of blade fluid interactions. This data can help support model development which will help predict turbine performance when analyzing perspective project locations in the future. Instrumentation was deployed to measure; water speed/direction, electrical power output, turbine shaft speed, turbine thrust force, and platform motion. Concurrent measurements of these parameters allow for correlations (cause and affect) to be inferred, allowing for characterization of device performance over a range of operating conditions. Water currents were measured using Acoustic Doppler Current Profilers (ADCP's) and Acoustic Doppler Velocimeters (ADV's) directly upstream and downstream of the turbine for inflow, wake and turbulence measurements. Electrical power output was measured using the Voltsys rectifier and the Shark power meter. Shaft speed was calculated based on the Voltsys measurements of the permanent magnet three phase generator AC generation frequency, coupled directly to the cross flow turbine under test (i.e., no gear box). Platform motions were captured using a Yost IMU (inertial measurement unit). Turbine thrust loading was measured using a reaction arm about the turbine deployment platform spanning beam, where two bi-directional load cells were connected to the system via a pinned connection. This submission includes zipped folders for each instrument containing quality controlled (QC'd) data in daily .csv files for the relevant duration specific to each instrument, along with separate .csv file that contains the units for each variable. Some instrument daily files are quite large and can pose a challenge for a visual spreadsheet editor to open. A processing software like MATLAB or Python is recommended. Note the degree of QC varied between each instrument due to time constraints. Particular time and attention was given to perform quality control tests on the acoustic based instruments that are particularly suscep...

16 TIDAL AND WAVE POWER↗

Tanana River Test Site Model Verification Using the Marine and Hydrokinetic Toolkit (MHKiT)

The marine energy (ME) industry historically lacked a standardized data processing toolkit for common tasks such as data ingestion, quality control, and visualization. The marine and hydrokinetic toolkit (MHKiT) solved this issue by providing a public software deployment (open-source and free) toolkit for the ME industry to store and maintain commonly used functionality for wave, tidal, and river energy. This paper demonstrates an initial model verification study in MHKiT. Using Delft3D, a numerical model of the Tanana River Test Site (TRTS) at Nenana, Alaska was created. Field data from the site was collected using an Acoustic Doppler Current Profiler (ADCP) at the proposed Current Energy Converter (CEC) locations. MHKiT is used to process model simulations from Delft3D and compare them to the transect data from the ADCP measurements at TRTS. The ability to use a single tool to process simulation and field data demonstrates the ease at which the ME industry can obtain results and collaborate across specialties, reducing errors and increasing efficiency.

13 HYDRO ENERGY↗

Ultrasonic microphone and ultrasonic acoustic radio

This disclosure provides systems, methods, and apparatus related to an ultrasonic microphone and an ultrasonic acoustic radio. In one aspect a system includes a transmitter and a receiver. The receiver comprises a membrane. The membrane comprises 5 a single layer or multiple layers of a two-dimensional material. The receiver is operable to receive sound waves in a frequency range, with the frequency range being the ultrasonic frequency range

42 ENGINEERING↗

Parallel-in-Time Solution of Hyperbolic PDE Systems via Characteristic-Variable Block Preconditioning

We consider the parallel-in-time solution of both linear and nonlinear hyperbolic partial differential equation (PDE) systems in one spatial dimension. In the nonlinear setting, the discretized equations are solved with a preconditioned residual iteration based on a global linearization. The linear(ized) equation systems are approximately solved parallel-in-time using a block preconditioner applied in the characteristic variables of the underlying linear(ized) hyperbolic PDE. This change of variables is motivated by the observation that intervariable coupling between characteristic variables is weak, at least locally where spatio-temporal variations in the eigenvectors of the associated flux Jacobian are sufficiently small, while that between the original variables is not. For an ℓ-dimensional system of PDEs, applying the preconditioner consists of solving a sequence of ℓ scalar linear(ized)-advection-like problems, each associated with a different characteristic wave-speed in the underlying linear(ized) PDE. Furthermore, we approximately solve these linear advection problems using multigrid reduction-in-time (MGRIT); however, any other suitable parallel-in-time method could be used. Numerical examples are shown for the (linear) acoustics equations in heterogeneous media and for the (nonlinear) shallow water equations and Euler equations of gas dynamics with shocks and rarefactions. For many test problems, the solver converges in just a handful of iterations and with mesh-independent convergence rates.

97 MATHEMATICS AND COMPUTING↗

Numerical investigation of laser-driven shock interaction with a deformable particle

A laser-driven shock propagating through an isolated particle embedded in a plastic (CH) target was studied using the radiation-hydrodynamic code FLASH. Preliminary simulations using IONMIX equations of state (EOS) showed significant differences in the shock Hugoniot of aluminum compared to experimental data in the low-pressure regime [O(10) GPa], resulting in higher streamwise compression and deformation of an aluminum particle. Hence, a simple modification to the ideal gas EOS was developed and employed to describe the target materials and examine the particle dynamics. The evolution of the pressure field demonstrated a complex wave interaction, resulting in a highly unsteady particle drag which featured two drag minima due to shock focusing at the rear end of the particle and rarefaction stretching due to laser shut-off. Although ~30% lateral expansion and ~25% streamwise compression were observed, the aluminum particle maintained considerable integrity without significant distortion. Additional simulations examined the particle response for a range of particle densities, sizes, and acoustic impedances. The data revealed that lighter particles such as aluminum gained significant momentum, reaching up to ~96% of the shocked CH's speed, compared to ~29% for the heavier tungsten particles. Despite the differences seen in the early stage of shock interaction, particles with varying acoustic impedances ultimately reached the same peak velocity. This identified particle-to-host density ratio is an important factor in determining the inviscid terminal velocity of the particle. In addition, the modified EOS model presented in this study could be used to approximate solid materials in hydrocodes that lack material strength models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Characterization of nonlinear ultrasonic waves behavior while interacting with poor interlayer bonds in large-scale additive manufactured materials

Over the past decades, researchers have developed several nonlinear ultrasonic techniques for quality control of materials commonly used in different applications. Owing to the superior sensitivity of nonlinear ultrasound waves to small defects such as micro-cracks, their applicability in different nondestructive testing (NDT) problems has been investigated in numerous studies. These studies utilize frequency domain analysis to detect the generation of higher harmonics because of the formation of defects in the inspected medium. Frequency domain analysis based on the Fourier transform is a significant approach used in linear systems; however, it may not perform adequately on nonlinear systems. Hence, studies on nonlinear dynamics and physics consider analyzing systems' behavior in the phase-space domain. In contrast to the frequency domain analysis, which can result in information loss, analysis in the phase-space domain retains all the information regarding a system's states. Here, we investigate the nonlinearities induced by poor interlayer bonding in polymer-based additive manufactured parts in the phase-space domain. It is convenient to characterize the nonlinearity in the phase-space domain because it provides a geometrical representation of a system's states. Two types of low quality interlayer bond are considered. The first type is simulated artificially while the second type is manufactured by reducing the bond quality during the printing process. The analysis verified that the received ultrasonic signals exhibit classical nonlinear behavior in the phase-space domain while interacting with simulated poor interlayer bonds. In addition, the results showed that the behavior of ultrasonic waves is amplitude-dependent and evolves into models that have not been previously reported. Furthermore, Largest Lyapunov Exponent (LLE) is used to quantify the behavior of nonlinear ultrasonic waves while interacting with poor interlayer bonds. Using LLE, it was observed that the divergence rate of the phase-space trajectories depends on the amplitude of the excitation. This observation quantitatively proves that nonlinear behavior of ultrasound while interacting with poor interlayer bonds can be amplitude-dependent. The results of both simulated and inherent poor interlayer bond cases showed that LLE can be used as a reliable quantitative damage-sensitive feature to detect and potentially characterize weak bonds, which are difficult to detect using conventional approaches. Additionally, the reported results in the phase-space domain provide a basis for proposing a new mathematical model for ultrasonic waves interacting with poor interlayer bonds.

36 MATERIALS SCIENCE↗