Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Elastic energy”

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 19 records

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↗

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↗

Absorptive corrections to the electromagnetic form factor in high-energy elastic proton-proton scattering

Recently, it was noted that absorptive corrections to the electromagnetic form factor in high-energy proton-proton scattering are important for the theoretical interpretation of the p ↑ p and p ↑ A analyzing power A N ( t ) measurements with the Hydrogen Jet Target polarimeter (HJET) at RHIC. Here, a concise expression for the absorptive correction was derived within the eikonal approach. The resulting analysis reveals a systematic bias, nearly independent of the beam energy, in the experimental determination of the real-to-imaginary ratio ρ when absorption effects are overlooked in the data analysis. Quantification of this bias, as ρ meas = ρ + ( 0.036 ± 0.016 ) bias , was achieved using a Regge fit applied to available proton-proton measurements of ρ meas ( s ) and σ tot meas ( s ) . Considering the potential impact of such an effect on the experimentally determined A N ( t ) , one may enhance consistency between the HJET and STAR measurements of the hadronic spin-flip amplitude. While the sign of the bias in the value of ρ aligns with the anticipated effective increase in the proton charge radius in p p scattering due to absorption, it amplifies the observed discrepancy between σ tot meas and ρ meas values at s = 13 TeV as measured in the TOTEM experiment. Evaluation (using published TOTEM data) of the measured proton-proton d σ / d t dependence on the absorptive corrections indicated that possible soft photon corrections to the hadronic amplitude slope may be essential for such data analysis. Published by the American Physical Society 2024

43 PARTICLE ACCELERATORS↗

Form Factors and Two-Photon Exchange in High-Energy Elastic Electron-Proton Scattering

We report new precision measurements of the elastic electron-proton scattering cross section for momentum transfer squared (Q$^2$) up to 15.75 (GeV/c)$^2$. These data allow for improved extraction of the proton magnetic form factor at high Q$^2$ and nearly double the Q$^2$ range of direct longitudinal/transverse separated cross sections. A comparison of our results to polarization measurements establishes the presence of hard two-photon exchange in the $e$-$p$ elastic scattering cross section at > 95 % confidence level for Q$^2$ up to 8 (GeV/c)$^2$.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Ferri-ionic coupling in Cu In P 2 S 6 nanoflakes: Polarization states and controllable negative capacitance

We consider nanoflakes of van der Waals ferrielectric Cu In P 2 S 6 covered by an ionic surface charge and reveal the appearance of polar states with a relatively large polarization of approximately 5 µC/ cm 2 and stored free charge of around 10 µC/ cm 2 , which can mimic “midgap” states associated with a surface-field-induced transfer of Cu and/or In ions in the van der Waals gap. The change in the ionic screening degree and mismatch strains induce a broad range of the transitions between paraelectric phase, antiferroelectric, ferrielectric, and ferri-ionic states in Cu In P 2 S 6 nanoflakes. The states’ stability and/or metastability is determined by the minimum of the system free energy consisting of electrostatic energy, elastic energy, and a Landau-type four-well potential of the ferrielectric dipole polarization. The possibility of governing the transitions by strain and ionic screening can be useful for controlling the tunneling barrier in thin-film devices based on Cu In P 2 S 6 nanoflakes. Additionally, we predict that the Cu In P 2 S 6 nanoflakes reveal features of the controllable negative capacitance effect, which make them attractive for advanced electronic devices, such as nanocapacitors and gate oxide nanomaterials with reduced heat generation. Published by the American Physical Society 2024

Physics↗

The effects of introducing elasticity using different interpolation schemes to the grand potential phase field model

Introducing elastic energy in the phase field method has been shown to influence interfacial energy, depending on the elastic interpolation scheme. This study investigates the impact of the elastic energy when using a grand potential-based phase field method, comparing the result of Khachaturyan’s strain interpolation scheme (KHS) and Voight-Taylor’s elastic energy interpolation scheme (VTS). The KHS model leads to a decrease in the interfacial energy, while the VTS model leads to an increase. The change in interfacial energy is greater with the VTS model than the KHS model, which suggests that the KHS model is more appropriate to limit the artificial impact of the elastic energy on the interfacial energy. When the contribution at the interface is not negligible, it is shown that both the microstructure evolution kinetics and the equilibrium microstructure can be influenced by the choice of the elastic scheme being used. In addition, this paper shows that the grand potential model might not be appropriate when the system requires the introduction of a composition-dependent term in the elastic energy contribution. This limitation is due to the need for an explicit and invertible relation between the total potential and the composition.

36 MATERIALS SCIENCE↗

Comparison of excess free energy at an interface according to the applied interpolation scheme for elasticity: A phase-field method

Phase-field modeling is an effective simulation technique for modeling microstructure evolution of elastically anisotropic systems. To introduce the elastic energy contribution in a phase field model, an interpolation scheme is used to define the mechanical properties within the phases and across the continuous interface. Several existing interpolation schemes introduce a potential excess elastic energy at the interface, which undesirable effect on microstructure evolution needs to be evaluated. In this study, we focused on three interpolation schemes including Khachaturyan’ scheme (KHS), Voigt–Taylor’s scheme (VTS), and Steinbach–Apel’s scheme (SAS). Comparisons of these schemes’ performances were performed in three configuration types using the MOOSE (Multiphysics Object-Oriented Simulation Environment) framework: bi-crystal, isotropic particle-matrix and anisotropic particle-matrix. The contribution of excess elastic energy on the interface energy as a function of interface width and the computational time to steady-state were evaluated in these three configurations. SAS introduces the lowest excess elastic energy contribution and the VTS has the biggest contribution amongst the considered schemes. Moreover, when modeling precipitation in an anisotropic elastic material, the SAS approach seems to predict more physical convex shapes during growth, making it preferable to KHS and VTS. Finally, as currently implemented, SAS requires the largest computational time and KHS requires the smallest time to reach steady-state amongst the considered schemes.

36 MATERIALS SCIENCE↗

Tuning energy barriers by doping 2D group-IV monochalcogenides

Structural degeneracies underpin the ferroic behavior of anisotropic next-generation two-dimensional materials and lead to peculiar two-dimensional structural transformations under external fields, charge doping, and/or temperature. The most direct indicator of the ease of these transformations is an elastic energy barrier, defined as the energy difference between the (degenerate) structural ground state unit cell and a unit cell with an increased structural symmetry. Proximity of a two-dimensional material to a bulk substrate can affect the magnitude of the critical fields and/or temperature at which these transformations occur, with the first effect being a relative charge transfer, which could trigger a structural quantum phase transition. With this physical picture in mind, we report the effect of modest charge doping (within -0.2 and +0.2 electrons per unit cell) on the elastic energy barrier J s of ferroelastic black phosphorene and nine ferroelectric/ferroelastic monochalcogenide monolayers. J s is the energy needed to create a Pnm2 1 →P4/nmm two-dimensional structural transformation, and it is sensitive to the orbital character of the electronic charge added or removed. Similar to the effect on the elastic energy barrier of ferroelastic SnO monolayers, group-IV monochalcogenide monolayers show a tunable elastic energy barrier for similar amounts of doping, and a decrease (increase) of J s can be engineered under a modest hole (electron) doping of not more than one-tenth of an electron or a hole per atom.

36 MATERIALS SCIENCE↗

Effect of recycled fibers and shredded intermediates variation on the mechanical properties and energy absorption of fiber‐reinforced composite panels

Abstract The global composite industry generates large quantities of waste which mostly ends up in landfills due to a lack of established end‐use applications for multiple waste streams. The scrap from end‐of‐life (EoL) includes manufacturing waste such as dry chopped fiber tows, loose fibers, shredded fibers from fabric textile operations, cured/semi‐cured prepregs, and fully cured composite structure waste from aircraft, automobiles, wind blades, boats, and pressure vessels. In this work, different composite waste streams were reduced to shredded intermediates, followed by simple blending, and subjected to wet compression molding to produce composite panels. The panels/plaques were tested for mechanical properties (flexure and impact), fiber‐matrix wet‐out, and property bounds. It was found that wet‐compression molding was a viable and scale‐able approach to produce recycled panels from EoL composites shredded scrap. Furthermore, full‐scale size panels for use in truck bodies and intermodal shipping container flooring were manufactured and their impact resistance was tested using a drop weight impact test. They were tested both for high‐ and low‐velocity load. In the case of high‐velocity load, the average impact load was 14,673 N; the average absorbed energy was 101.6 J; the average elastic energy was 11.7 J and the impact resistance was 1065 J/m. In the case of the low‐velocity drop weight impact test, it was found that the average impact load was 7877.358 N; the average absorbed energy was 9.718 J; the average elastic energy was 9.14 J, and the impact resistance was 184.1 J/m. The shredded composite was shown to be a candidate material for the manufacture of truck bodies and intermodal containers’ flooring panels. Highlights By using shredded intermediates from different composite waste streams, it is possible to manufacture composite panels. Wet–compression process is an appropriate technique for manufacturing recycled fiber composite panels. Regardless of the source of scrap, the mechanical properties of the produced composite panels were improved. The recycling process technology can be transformed to commercial scale to produce full‐size transportation flooring panels. A product pathway is established in consideration of lower cost and improved recyclability.

Vaidya, Uday↗

Nanocrystal Geometry Governs Phase Transformation Pathways in Palladium Hydride

Pathways and structural dynamics of phase transformations impact performance of materials in energy and information storage technologies. Palladium hydride (PdH x ) nanocrystals are an ideal model system for studying solute-induced phase transformations, where elastic energy from lattice mismatch between α-PdH x and β-PdH x phases is often considered a key to determining the transformation pathways. α/β-PdH x interfacial elastic energy is affected by the confined geometry of a nanocrystal. However, how nanocrystal geometry influences phase transformation pathways is largely unknown. Using in situ liquid phase transmission electron microscopy, we directly visualize hydrogenation in Pd nanocrystals with two geometries, a nanocube and a hexagonal nanoplate. Both follow similar sequences of an initially curved nucleus, interface flattening, and reverse-stage nucleation; however, their evolving α/β-PdH x interfaces exhibit geometry-dependent crystallographic alignments. In nanocubes, {100}-aligned configurations conform to static elastic energy ordering, representing a pathway that maintains a local mechanical equilibrium, whereas nanoplates display both {110}- and {211}-aligned interfaces. Theoretical simulations show that geometry determines the accessibility of alternative phase transformation pathways as the system is driven far from equilibrium during hydrogenation. In conclusion, these findings identify geometry as a fundamental parameter for directing phase transformation pathways, offering design principles for accessing atypical configurations and improving properties of intercalation-based devices.

Lee, Daewon [Lawrence Berkeley National Laboratory↗

Novel ceramic capacitors with ultrahigh energy density and efficiency (Final Technical Report)

Antiferroelectric ceramics are a special class of material that have shown great potential as the dielectric in electrical capacitors due to their high energy- and power-density. During each charge-discharge cycle, the ceramic undergoes transformation to a ferroelectric phase and resumes its antiferroelectric phase. The hysteresis associated with the transitions leads to a mediocre energy efficiency and service lifetime of antiferroelectric capacitors and, hence, their almost absence in commercial products. Under the support of this research project, we first formulated a universal lattice-compatibility theory that included electrostatic polarization energy along with elastic energy and thermal energy to understand the origin of the hysteresis in antiferroelectric oxides. Guided by this compatibility theory, we conducted high-throughput density functional theory (DFT) calculations to assess chemical modifiers and their effect on crystal structures of 400+ PbZrO 3 -based compositions. Down-selected compositions were experimentally validated for their suppressed hysteresis and higher energy efficiency. The verified low-hysteresis compositions were then expanded to an antiferroelectric ceramic library with nearly 500 new compositions (more than 1,500 samples) using high-throughput experiments involving ceramic synthesis and property screening. The large quantity of data generated (theory and experimental) in these tasks were processed by machine-learning techniques and identified trends were fed to the next iteration. In the end, we successfully discovered four compositions with near-zero hysteresis, yielding a world-record energy efficiency of 98.2% at an energy density of 3.0 J/cm 3 . Furthermore, our antiferroelectric ceramic capacitor reaches 79.5 million charge-discharge cycles lifetime, a factor of 80 enhancement over previous antiferroelectric ceramics with large hysteresis. These research accomplishments have not only met the milestones set in the SOPO, but also led to two patent filings, three journal publications (one of them was in Advanced Materials, impact factor 29.4), and nine oral presentations at various venues. Through the course of the project, three postdocs, four Ph.D. students, and one M.S. student were trained. In short, our project established a new methodology in searching next-generation functional ceramics on the fundamental side and discovered several high-efficiency antiferroelectric compositions for capacitors on the applied side. Once fabricated into the multilayer form for commercial applications, these ceramic capacitors can potentially enable the high temperature high power density DC-link capacitors that are critical for the next generation inverters in electric vehicles. The project also significantly contributed to the nation’s workforce development in the STEM fields.

36 MATERIALS SCIENCE↗

Clarifying the definition of ‘transonic’ screw dislocations

We report a number of recent Molecular Dynamics (MD) simulations have demonstrated that screw dislocations in face centred cubic (fcc) metals can achieve stable steady state motion above the lowest shear wave speed (v shear ) which is parallel to their direction of motion (often referred to as transonic motion). This is in direct contrast to classical continuum analyses which predict a divergence in the elastic energy of the host material at a crystal geometry dependent ‘critical’ velocity v crit . Within this work, we first demonstrate through analytic analyses that the elastic energy of the host material diverges at a dislocation velocity (v crit ) which is greater than v shear , i.e. v crit >v shear . We argue that it is this latter derived velocity (v crit ) which separates ‘subsonic’ and ‘supersonic’ regimes of dislocation motion in the analytic solution. In addition to our analyses, we also present a comprehensive suite of MD simulation results of steady state screw dislocation motion for a range of stresses and several cubic metals at both cryogenic and room temperatures. At room temperature, both our independent MD simulations and the earlier works find stable screw dislocation motion only below our derived v crit . Nonetheless, in real-world polycrystalline materials v crit cannot be interpreted as a hard limit for subsonic dislocation motion. In fact, at very low temperatures our MD simulations of Cu at 10 Kelvin confirm a recent claim in the literature that true ‘supersonic’ screw dislocations with dislocation velocities v>v crit are possible at very low temperatures.

74 ATOMIC AND MOLECULAR PHYSICS↗

Velocity map imaging of inelastic and elastic low energy electron scattering in organic nanoparticles

Electron transport is of fundamental importance, and has application in a variety of fields. Different scattering mechanisms affect electron transport in solids. It is important to comprehensively understand these mechanisms and their scattering cross-sections to predict electron transport properties. Whereas electron transport is well understood for high kinetic energy (KE) electrons, there are inconsistencies in the low KE regime. In this study, velocity map imaging soft X-ray photoelectron spectroscopy is applied to unsupported organic nanoparticles to extract experimental values of inelastic and elastic mean free paths. The obtained data is used to calculate corresponding scattering cross-sections. The data demonstrates a decrease of the Inelastic Mean Free Path and increase of the Elastic Mean Free Path with increasing electron KE between 10-50 eV.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Low-Energy Electron Elastic Total Cross Sections for Ho, Er, Tm, Yb, Lu, and Hf Atoms

The robust Regge-pole methodology wherein is fully embedded the essential electron-electron correlation effects and the vital core polarization interaction has been used to explore negative ion formation in the large lanthanide Ho, Er, Tm, Yb, Lu, and Hf atoms through the electron elastic total cross sections (TCSs) calculations. These TCSs are characterized generally by dramatically sharp resonances manifesting ground, metastable, and excited negative ion formation during the collisions, Ramsauer-Townsend minima, and shape resonances. The novelty and generality of the Regge-pole approach is in the extraction of the negative ion binding energies (BEs) of complex heavy systems from the calculated electron TCSs. The extracted anionic BEs from the ground state TCSs for Ho, Er, Tm, Yb, Lu, and Hf atoms are 3.51 eV, 3.53 eV, 3.36 eV, 3.49 eV, 4.09 eV and 1.68 eV, respectively. The TCSs are presented and the extracted from the ground; metastable and excited anionic states BEs are compared with the available measured and/or calculated electron affinities. We conclude with a remark on the existing inconsistencies in the meaning of the electron affinity among the various measurements and/or calculations in the investigated atoms and make a recommendation to resolve the ambiguity.

Felfli, Zineb↗

SISGR: Chemomechanics of Far-From-Equilibrium Interfaces (COFFEI)

Portable, reliable, and deployable devices for energy storage and conversion require fundamental changes in design of solid-state composites comprising ceramics and metals. These materials comprise the electrodes and electrolytes of next-generation solid-oxide fuel cells and all solid-state batteries, forming solid-state functional composites. The advent of solid-state batteries – which replace liquid electrolytes with solid electrolytes capable of lithium ion transport for reliable energy storage in portable batteries – and the increased demand space for all solid-state fuel cells capable of oxygen reduction at intermediate temperatures remain important challenges for improved material stability and decreased system cost. However, little is understood about three fundamental facets of materials that enable such solid-state energy applications. First, how do such materials deform, fracture, or delaminate under operando conditions? Second, how does such mechanical deformation limit or facilitate electronic and ionic transport within and across such material interfaces? Third, how we can predictably design interface-rich composites to engineer both structural and electrochemical stability? This COFFEI Group comprised expertise from Materials Science & Engineering and Nuclear Science & Engineering to integrate unique in situ experiments, simulations, and fabricated interfaces that address these fundamental questions in solid-state interfaces of nanoscale composites that will guide solid-state electrochemistry, transport kinetics, and mechanical deformation for nonstoichiometric materials that enable such applications. In particular, we built on COFFEI’s understanding of chemomechanical coupling among defect concentrations, ionic transport, electron transport, and stored elastic energy that is particularly acute in the far-from-equilibrium conditions typical of energy device applications. By tailoring our focus to solid-state interfaces, we addressed these important issues by (a) developing and applying advanced in situ and ex situ characterization tools to characterize model materials and interfaces synthesized with molecular-level control, under both laboratory-controlled and extreme environments representative of energy device operation; and (b) employing computational modeling and simulation frameworks to predict transport mechanisms, reactivity and stability of these model materials and interfaces under significant chemical strains typical of energy device operation. Recent progress provided insights to additional materials systems and electrochemomechanical fatigue and fracture that were not fully envisioned when the program was initiated. Specifically, in the final three years of COFFEI we pursued two integrated thrusts, with complementary focus. Thrust I focused on failure-resistant electrochemomechanical composites, while Thrust II focused on strain-modulated conductivity and reactivity across interfaces. In contrast to our initial COFFEI focus, these thrusts concentrated wholly on solid-state material interfacial interactions and included greater integration of multiscale visualization including in situ electron microscopy of strained structures/interactions and mesoscale simulations. Successful development of functionally superior and long-lived battery and fuel cell systems and stress adaptable oxides requires a deeper, fundamental understanding of the coupling among the historically important subfields of solid-state electrochemistry, transport kinetics, and mechanical deformation for nonstoichiometric metal oxide electrodes. In this program, the understanding and the application of chemomechanical coupling of defect concentrations, ionic transport, electro-catalytic activity and stored elastic energy, particularly acute in the far-from-equilibrium conditions typical of energy device applications, are being refined and implications for device operation clarified, including for miniaturized solid-state batteries and fuel cells.

36 MATERIALS SCIENCE↗

Low-Energy Electron Elastic Collisions with Actinide Atoms Am, Cm, Bk, Es, No and Lr: Negative-Ion Formation

The rigorous Regge-pole method is used to investigate negative-ion formation in actinide atoms through electron elastic total cross sections (TCSs) calculation. The TCSs are found to be characterized generally by negative-ion formations, shape resonances and Ramsauer-Townsend(R-T) minima, and they exhibit both atomic and fullerene molecular behavior near the threshold. Additionally, a polarization-induced metastable cross section with a deep R-T minimum is identified near the threshold in the Am, Cm and Bk TCSs, which flips over to a shape resonance appearing very close to the threshold in the TCSs for Es, No and Lr. We attribute these new manifestations to size effects and orbital collapse significantly impacting the polarization interaction. From the TCSs unambiguous and reliable ground, metastable and excited states negative-ion binding energies (BEs) for Am−, Cm−, Bk−, Es−, No− and Lr− anions formed during the collisions are extracted and compared with existing electron affinities (EAs) of the atoms. The novelty of the Regge-pole approach is in the extraction of the negative-ion BEs from the TCSs. We conclude that the existing theoretical EAs of the actinide atoms and the recently measured EA of Th correspond to excited anionic BEs.

Msezane, Alfred Z. (ORCID:0000000212485661)↗

Thermodynamic analysis of hydrogel swelling in aqueous sodium chloride solutions

Here, a comprehensive thermodynamic model is presented for swelling and salt partitioning of poly (n-isopropyl acrylamide) hydrogels in aqueous sodium chloride (NaCl) solutions at 298 K. Together with a Helmholtz energy expression for network elastic energy, a modified electrolyte Nonradom Two-Liquid activity coefficient model is used to represent both the long-range electrostatic interactions and the short-range van der Waals interactions present in the gel phase. The model parameters include one network parameter per hydrogel system plus two binary parameters per interaction pair in the polymer–solvent-salt systems. With semiquantitave agreement with experimental data for both hydrogel swelling and salt partitioning up to salt saturation, the analysis suggests the ion hydration between NaCl and water is the root cause for hydrogel deswelling in aqueous NaCl solution at high NaCl concentrations, i.e., NaCl weight fraction >0.03. On the other hand, the strong repulsive interaction between NaCl and hydrogel polymer is the dominant factor for salt partitioning.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Nanometer-Scale Acoustic Wave Packets Generated by Stochastic Core-Level Photoionization Events

We demonstrate that the absorption of femtosecond hard x-ray pulses excites quasispherical, high-amplitude, and high-wave-vector coherent acoustic phonon wave packets using an all hard-x-ray pump-probe scattering experiment. The time- and momentum-resolved diffuse scattering signal is consistent with an ensemble of 3D strain wave packets induced by the rapid electron cascade dynamics following photoionization at uncorrelated excitation centers. We quantify key parameters of this process, including the localization size of the stress field and the photon energy conversion efficiency into elastic energy. The parameters are determined by the photoelectron and Auger electron cascade dynamics, as well as the electron-phonon interaction. In particular, we obtain the localization size of the observed strain wave packet to be 1.5 and 2.5 nm for bulk SrTiO 3 and KTaO 3 single crystals, respectively. The results provide crucial information on the mechanism of x-ray energy deposition into matter and shed light on the shortest collective length scales accessible to coherent acoustic phonon generation using x-ray excitation. Published by the American Physical Society 2024

36 MATERIALS SCIENCE↗