Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mesoscale processes”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Analysis of a new multispecies tumor growth model coupling 3D phase-fields with a 1D vascular network

In this study, we present and analyze a mathematical model for tumor growth incorporating ECM erosion, interstitial flow, and the effect of vascular flow and nutrient transport. The model is of phase-field or diffused-interface type in which multiple phases of cell species and other constituents are separated by smooth evolving interfaces. The model involves a mesoscale version of Darcy’s law to capture the flow mechanism in the tissue matrix. Modeling flow and transport processes in the vasculature supplying the healthy and cancerous tissue, one-dimensional (1D) equations are considered. Since the models governing the transport and flow processes are defined together with cell species models on a three-dimensional (3D) domain, we obtain a 3D–1D coupled model.

3D–1D coupled blood flow models↗
DOE OSTI2020

Mesoscale informed parameter estimation through machine learning: A case-study in fracture modeling

Scale bridging is a critical need in computational sciences, where the modeling community has developed accurate physics models from first principles, of processes at lower length and time scales that influence the behavior at the higher scales of interest. However, it is not computationally feasible to incorporate all of the lower length scale physics directly into upscaled models. This is an area where machine learning has shown promise in building emulators of the lower length scale models, which incur a mere fraction of the computational cost of the original higher fidelity models. We demonstrate the use of machine learning using an example in materials science estimating continuum scale parameters by emulating, with uncertainties, complicated mesoscale physics. Additionally, we describe a new framework to emulate the fine scale physics, especially in the presence of microstructures, using machine learning, and showcase its usefulness by providing an example from modeling fracture propagation. Our approach can be thought of as a data-driven dimension reduction technique that yields probabilistic emulators. Our results show well-calibrated predictions for the quantities of interests in a low-strain simulation of fracture propagation at the mesoscale level. Furthermore, on average, we achieve ~10% relative errors on time-varying quantities like total damage and maximum stresses. Successfully replicating mesoscale scale physics within the continuum models is a crucial step towards predictive capability in multi-scale problems.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Global Precipitation Experiment—A New World Climate Research Programme Lighthouse Activity

The future state of the global water cycle and the prediction of freshwater availability for humans around the world remain among the challenges of climate research and are relevant to several United Nations Sustainable Development Goals. The Global Precipitation Experiment (GPEX) takes on the challenge of improving the prediction of precipitation quantity, phase, timing, and intensity, characteristics that are products of a complex integrated system. It will achieve this by leveraging existing World Climate Research Programme (WCRP) activities and community capabilities in satellite, surface-based, and airborne observations, modeling, and experimental research and by conducting new and focused activities. It was launched in October 2023 as a WCRP Lighthouse Activity. Here, we present an overview of the GPEX Science Plan that articulates the primary science questions related to precipitation measurements, process understanding, model performance and improvements, and plans for capacity development. The central phase of GPEX is the WCRP Years of Precipitation for 2–3 years with coordinated global field campaigns focusing on different storm types (atmospheric rivers, mesoscale convective systems, monsoons, and tropical cyclones, among others) over different regions and seasons. Activities are planned over the three phases (before, during, and after the Years of Precipitation) spanning a decade. These include gridded data evaluation and development, advanced modeling, enhanced understanding of processes critical to precipitation, multiscale prediction of precipitation events across scales, and capacity development. These activities will be further developed as part of the GPEX Implementation Plan.

Climate change↗
DOE OSTI2022

Simulating Extreme Precipitation in the United States in the Energy Exascale Earth System Model: Investigating the Importance of Representing Convective Intensity versus Dynamic Structure (Final Report)

Weather events that produce extreme precipitation are associated with severe flooding and winds that result in thousands of deaths and billions of dollars in damages annually in the United States (U.S.). This project aims to improve understanding of the small- and large-scale processes that govern these events and improve our ability to project changes in these extreme events under the influences of natural variability and human activities. We focus on two promising directions in the development of the U.S. Department of Energy’s Energy Exascale Earth System Model (E3SM) to investigate the tradeoffs between resolving the convective-scale processes that control the intensity versus the intermediate-scale processes that control the dynamic structure of the most prominent extreme precipitation events that impact the U.S. throughout the year (i.e., mesoscale convective systems, tropical cyclones, and extratropical cyclone).

58 GEOSCIENCES↗
DOE OSTI2023

Improving the representation of isopycnal mixing in E3SM (Final Report)

The DOE-supported Energy Exascale Earth System Model (E3SM) is a major attempt by the Department of Energy to develop a new Earth System Model using unstructured grids, which allow for the model to concentrate resolution where it is most needed and avoid numerical artifacts experienced by previous Earth System Models (ESMs) using regular grids in which the lateral resolution becomes extremely fine as the resolution shrinks at the poles. This effort required rewriting many of the algorithms previously developed for regular grids. One algorithm that was problematic in the first version of E3SM was the representation of mixing due to turbulent ocean eddies with scales smaller than the model grid. These mesoscale eddies are the primary way by which tracers are stirred along density surfaces in the ocean interior. In particular, the process of isopycnal mixing exchanges fresher, more oxygenated waters from polar regions with saltier, nutrient-rich and oxygen poor waters from tropical regions. Previous work in our group has shown that this mechanism is vitally important for bringing oxygen into poorly ventilated tropical regions (Gnanadesikan et al., 2012, 2013) and plays a significant role in determining the uptake of anthropogenic carbon dioxide (Gnanadesikan et al., 2015). However, in E3SMv1 this process was turned off, as turning it on caused the model to become unstable. Additionally, the rate of isopycnal mixing is directly proportional to a mixing coefficient $A_{Redi}$ whose value varies from less than 400 m 2 /s to 2000 m 2 /s across contemporary climate models. The main thrusts of this proposal therefore were to: 1.) Improve the numerics of isopycnal mixing in E3SM; 2.) Identify ways in which climate models are sensitive to the isopycnal mixing coefficient; and 3.) Develop better representations of the isopycnal mixing coefficient.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Hidden domain boundary dynamics toward crystalline perfection

A central paradigm of nonequilibrium physics concerns the dynamics of heterogeneity and disorder, impacting processes ranging from the behavior of glasses to the emergent functionality of active matter. Understanding these complex mesoscopic systems requires probing the microscopic trajectories associated with irreversible processes, the role of fluctuations and entropy growth, and the timescales on which nonequilibrium responses are ultimately maintained. Approaches that illuminate these processes in model systems may enable a more general understanding of other heterogeneous nonequilibrium phenomena, and potentially define ultimate speed and energy cost limits for information processing technologies. Here, we apply ultrafast single-shot X-ray photon correlation spectroscopy to resolve the nonequilibrium, heterogeneous, and irreversible mesoscale dynamics during a light-induced phase transition in a (PbTiO 3 ) 16 /(SrTiO 3 ) 16 superlattice. Such ferroelectric superlattice systems are a useful platform to study phase transitions and topological dynamics due to their high degree of tunability. This provides an approach for capturing the nucleation of the light-induced phase, the formation of transient mesoscale defects at the boundaries of the nuclei, and the eventual annihilation of these defects, even in systems with complex polarization topologies. We identify a nonequilibrium correlation response spanning >10 orders of magnitude in timescales, with multistep behavior similar to the plateaus observed in supercooled liquids and glasses. We further show how the observed time-dependent long-time correlations can be understood in terms of stochastic and non-Markovian dynamics of domain walls, encoded in waiting-time distributions with power-law tails. This work defines possibilities for probing the nonequilibrium and correlated dynamics of disordered and heterogeneous media.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Mechanism of Mesoscale Woodpile Development via Photoelectrochemical Deposition of Se–Te

A combination of experiments and optical modeling provided insight into the mechanism of mesoscale woodpile formation in response to an orthogonal shift in polarization during photoelectrochemical deposition of Se–Te. Cathodic deposition of semiconducting Se–Te using spatially uniform, linearly polarized illumination produced arrays of lamellae that were aligned parallel to the optical E-field oscillation. Continued deposition in conjunction with an orthogonal shift in the polarization direction then produced aligned bridging features that spanned the void space between, and were orthogonal to, the preexisting lamellae. The height and pitch, respectively, in each layer of the woodpile were a function of the charge density and illumination wavelength during deposition. A Monte Carlo model, in which material addition was scaled by the absorption magnitude obtained from electromagnetic simulations, produced morphologies that were nominally identical to those observed experimentally. Here, the formation of mesoscale woodpiles is consistent with a mechanism that involves a series of spontaneously initiated, concerted light–matter interactions during the photoelectrochemical deposition process.

absorption↗
DOE OSTI2022

Learning two-phase microstructure evolution using neural operators and autoencoder architectures

Abstract Phase-field modeling is an effective but computationally expensive method for capturing the mesoscale morphological and microstructure evolution in materials. Hence, fast and generalizable surrogate models are needed to alleviate the cost of computationally taxing processes such as in optimization and design of materials. The intrinsic discontinuous nature of the physical phenomena incurred by the presence of sharp phase boundaries makes the training of the surrogate model cumbersome. We develop a framework that integrates a convolutional autoencoder architecture with a deep neural operator (DeepONet) to learn the dynamic evolution of a two-phase mixture and accelerate time-to-solution in predicting the microstructure evolution. We utilize the convolutional autoencoder to provide a compact representation of the microstructure data in a low-dimensional latent space. After DeepONet is trained in the latent space, it can be used to replace the high-fidelity phase-field numerical solver in interpolation tasks or to accelerate the numerical solver in extrapolation tasks.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Tunability and Long-Range Enhancement of Resonance Energy Transfer Facilitated by Plasmonic Nanorods

Resonance energy transfer (RET) between molecules or quantum dots is an important process in many energy-related applications. Different environmental structures have been studied and demonstrated to be able to enhance the RET rates between a nearby donor–acceptor pair. In particular, cylindrical silver nanorods and nanowires have shown an extraordinary ability to transfer energy along their longitudinal axes over large distances. However, the detailed mechanism of such transfer and the effects on the molecular RET process are as yet elusive. In this study, we use the recently developed computational tool based on the plasmon-coupled RET method to systematically study the effects of nanorods with different dimensions on the RET rates. We find that highly frequency-dependent coupling factor (CF) spectra, whose amplitudes determine RET rates, can be obtained due to the localized surface plasmon polariton modes of the rods with nanoscale dimensions. Simple phenomenological models can be derived for the wavelengths of CF peaks in relation to the length and width of the nanorods, providing easy tunability for enhancing the RET rate in specific wavelength ranges. When coupled to longer rods with mesoscale lengths, exponential decay of the CF over long donor–acceptor distances with a small decay constant is observed, leading to the possibility of long-range RET processes. Furthermore, drumhead resonance modes emerge on the flat ends of the rod when the rod’s diameter reaches 300 nm, resulting in extra enhancement to RET rate compared to certain thinner rods. Furthermore, these findings shed new light on the mechanism of plasmonic enhancement with silver nanorods and establish design principles for how to optimally utilize these structures to manipulate RET processes for various applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Multiscale and Machine Learning Modeling for Process-informed Microstructure Prediction in Additively Manufactured Materials using MALAMUTE

The Advanced Materials and Manufacturing Technologies (AMMT) program under the Department of Energy Office of Nuclear Energy aims to develop and qualify additively manufactured materials for nuclear applications. One key challenge to this is the microstructural variability observed in the additively manufactured products and their impact on the properties and performance of the material in extreme environments. AMMT is using a combination of high-throughput experimental and modeling techniques to accelerate qualification. Conventionally, in-situ and ex-situ characterizations and testing are performed to correlate different aspects of the additive manufacturing process to the final product and its performance. However, adopting a trial-and-error approach to experimentally evaluate the vast range of process parameters required to capture microstructural variability is cost-prohibitive. Modeling and simulation provide a comparatively inexpensive way to understand and correlate the microstructural evolution to the processing conditions. The modeling and simulation work-packages within the AMMT program aims to use physics-based and machine learning models to develop a digital twin for additive manufacturing that can correlate the process conditions to the final product and establish a process-structure-property-performance (PSPP) correlation. The melting and subsequent solidification that occurs during the additive process is a complex phenomenon that requires multiscale multiphysics analysis. This work package focuses on understanding the role of process variabilities on the unique microstructural characteristics of additively manufactured materials. Microstructural features at the subgrain level, such as compositional micro-heterogeneity and dislocation cells, are of particular interest here since they can influence the creep properties and radiation performance. Idaho National Laboratory’s Multiphysics Object-Oriented Simulation Environment (MOOSE), specifically the MOOSE Application Library for Advanced Manufacturing UTilitiEs (MALAMUTE) software, provides an ideal platform for developing the multiphysics multiscale model to explore the intricacies of the microstructural evolution during the AM processes within a single framework. Furthermore, given that such full-fidelity simulations can be computationally intensive, reduced order models are necessary to explore the PSPP space for additively manufactured materials in an efficient, reliable, and cost-effective way. This work focuses on capturing the microstructural variabilities at the subgrain level that are often missing in the part-scale models. In fiscal year 2025, we significantly advanced upon our work in the last fiscal year, in terms of the predictive capabilities of the physics-based and ML models, by adding the capabilities to capture subgrain-level micro-segregation during solidification using phase-field model and to predict the time-dependent dynamics of the AM process through the MOGPAR model. The alloy solidification model in MOOSE incorporates the thermodynamic properties and free energy relevant to 316 stainless steel. The model demonstrates the Cr and Ni segregation that occurs during solidification, including that the rate of solidification. The microstructural evolution model is connected to the process conditions via the surrogate model developed in this work. This enables predictions of the final microstructure in conjunctions with the manufacturing process. This work supports AMMT's rapid qualification goals by laying the foundation for an efficient and cost-effective model establishing the PSPP correlation for AM. The generated microstructures and predicted micro-segregation can be used by other work packages under AMMT to evaluate the properties and environmental response of the material at the mesoscale. Thus, this work helps to identify the key microstructural features at the subgrain level that are significant in property and performance predictions of additively manufactured components. This work will also provide inputs to the large-scale process variability models to reevaluate and validate assumptions and simplifications made in the part-scale models. Furthermore, through active learning this work can help identify the data need from both modeling and experimental sides for development of a robust digital twin for additive manufacturing and accelerate the AMMT's qualification efforts.

36 - MATERIALS SCIENCE↗
DOE OSTI2024

Moisture Recycling through Pumping by Mesoscale Convective Systems

Abstract Moisture recycling, the contribution of local evapotranspiration (ET) to precipitation, has been studied using bulk models assuming a well-mixed atmosphere. The latter is inconsistent with a climatologically stratified atmosphere that slants across latitudes. Reconciling the two views requires an understanding of overturning associated with different weather systems. In this study, we aim to better understand moisture recycling associated with mesoscale convective systems (MCSs). Using a convection-permitting WRF simulation equipped with water vapor tracers (WRF-WVT), we tag moisture from terrestrial ET in the U.S. Southern Great Plains during May 2015, when more than 20 MCS events occurred and produced significant precipitation and flooding. Water budget analysis reveals that approximately 76% of terrestrial ET is advected away from the region while the remaining 24% of terrestrial ET is “pumped” upward within the region, accounting for 12% of precipitation. Moisture recycling peaks during early night hours (1800–2400 LT) due to the mixing of the daytime accumulated ET by active convection. By focusing on five “diurnally driven” MCSs with less large-scale circulation influence than other MCSs during the same period, we find an upright pumping of terrestrial ET at the MCS initiation and development stages, which diverges into two branches during the MCS mature and decaying stages. One branch in the upper level advects the ET-sourced moisture downstream, while the other branch in the mid-to-upper level contributes to the trailing precipitation upstream. Overall, our analysis depicts a pumping mechanism associated with MCSs that mixes local ET vertically, highlighting its specific contributions to enhancing convective precipitation processes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Microstructure based modeling of friction stir welded joint of dissimilar metals using crystal plasticity

We report the friction stir welding (FSW) process shows promising results in joining dissimilar metals which are otherwise almost impossible to join using traditional welding techniques. Being a new technique, the deformation and the failure mechanism of the joints made by the FSW process needs to be investigated. In this work, a joint between AZ31 Mg alloy and DP590 steel is modeled using phenomenological crystal plasticity formulation on the mesoscale in the form of a representative volume element (RVE). The interface of the two materials is modeled using a cohesive zone model. A parametric study has been performed to understand the effect of grain size and interface fracture toughness as well as strength on the mechanical performance of the joint. It was found that the grain size of AZ31 Mg alloy, as well as DP590 steel, has little effect on the overall joint performance. On the other hand, interface fracture toughness and strength have a significant impact on the mechanical properties of the joint

36 MATERIALS SCIENCE↗
DOE OSTI2023

Effects from Time Dependence of Ice Nucleus Activity for Contrasting Cloud Types

Abstract The role of time-dependent freezing of ice nucleating particles (INPs) is evaluated with the “Aerosol–Cloud” (AC) model in 1) deep convection observed over Oklahoma during the Midlatitude Continental Convective Cloud Experiment (MC3E), 2) orographic clouds observed over North California during the Atmospheric Radiation Measurement (ARM) Cloud Aerosol Precipitation Experiment (ACAPEX), and 3) supercooled, stratiform clouds over the United Kingdom, observed during the Aerosol Properties, Processes And Influences on the Earth’s climate (APPRAISE) campaign. AC uses the dynamical core of the WRF Model and has hybrid bin–bulk microphysics and a 3D mesoscale domain. AC is validated against coincident aircraft, ground-based, and satellite observations for all three cases. Filtered concentrations of ice (>0.1–0.2 mm) agree with those observed at all sampled levels. AC predicts the INP activity of various types of aerosol particles with an empirical parameterization (EP), which follows a singular approach (no time dependence). Here, the EP is modified to represent time-dependent INP activity by a purely empirical approach, using our published laboratory observations of time-dependent INP activity. In all simulated clouds, the inclusion of time dependence increases the predicted INP activity of mineral dust particles by 0.5–1 order of magnitude. However, there is little impact on the cloud glaciation because the total ice is mostly (80%–90%) from secondary ice production (SIP) at levels warmer than about −36°C. The Hallett–Mossop process and fragmentation in ice–ice collisions together initiate about 70% of the total ice, whereas fragmentation during both raindrop freezing and sublimation contributes <10%. Overall, total ice concentrations and SIP are unaffected by time-dependent INP activity. In the simulated APPRAISE case, the main causes of persistence of long-lived clouds and precipitation are predicted to be SIP in weak embedded convection and reactivation following recirculation of dust particles in supercooled layer cloud.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Robotic Planning under Uncertainty in Spatiotemporal Environments in Expeditionary Science

In the expeditionary sciences, spatiotemporally varying environments -- hydrothermal plumes, algal blooms, lava flows, or animal migrations -- are ubiquitous. Mobile robots are uniquely well-suited to study these dynamic, mesoscale natural environments. We formalize expeditionary science as a sequential decision-making problem, modeled using the language of partially-observable Markov decision processes (POMDPs). Solving the expeditionary science POMDP under real-world constraints requires efficient probabilistic modeling and decision-making in problems with complex dynamics and observational models. Previous work in informative path planning, adaptive sampling, and experimental design have shown compelling results, largely in static environments, using data-driven models and information-based rewards. However, these methodologies do not trivially extend to expeditionary science in spatiotemporal environments: they generally do not make use of scientific knowledge such as equations of state dynamics, they focus on information gathering as opposed to scientific task execution, and they make use of decision-making approaches that scale poorly to large, continuous problems with long planning horizons and real-time operational constraints. In this work, we discuss these and other challenges related to probabilistic modeling and decision-making in expeditionary science, and present some of our preliminary work that addresses these gaps. We ground our results in a real expeditionary science deployment of an autonomous underwater vehicle (AUV) in the deep ocean for hydrothermal vent discovery and characterization. Our concluding thoughts highlight remaining work to be done, and the challenges that merit consideration by the reinforcement learning and decision-making community.

Preston, Victoria↗
DOE OSTI2021

Influence of the distribution of expansive sites in aggregates on microscopic damage caused by alkali-silica reaction: Insights into the mechanical origin of expansion

The origin of damage in concrete due to the alkali-silica reaction (ASR) is attributed to the expansion site in the aggregate. To investigate the cracking process of the aggregate during ASR and its consequences on concrete damage, the effect of the distribution of the expansive sites in the aggregate on ASR expansion and the crack patterns must be evaluated. Thus, in this study, a mesoscale discrete model was applied to ASR modeling to represent the propagation of cracking during ASR accurately. The distribution of the expansive sites in the aggregate was based on the gel pocket and reaction rim models, which are two ASR mechanisms reported in the literature. These two expansion models highlight the different crack patterns obtained based on the aggregate characteristics. Further, the expansion cracking processes determined based on the gel pocket and reaction rim models are consistent with the evolution of cracking with the expansion level.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Surface engineering to enhance heat generation and joint strength in dissimilar materials AZ31 and DP590 ultrasonic welding

A multiscale simulation approach was developed and employed to optimize the sheet surface conditions for higher interfacial temperature and joint strength in ultrasonic welding of magnesium alloy AZ31 and dual-phase steel DP590. First, a mesoscale model was used to study the relationship between friction coefficient and surface roughness, which can be modified by various engineering methods. Then a macroscopic process model was employed to study the effects of surface roughness on heat generation, indicating that a temperature increase can be achieved with rougher surfaces on two sides of both DP590 and AZ31 sheets. Samples prepared by sanding and filing, as well as grinding, were first characterized for surface roughness and then welded under ultrasonic vibration. An infrared camera was used to measure temperatures in situ for model validation. Overall, lap shear test results for the welded joint showed that the joint strength can be improved by 10~25% using filing and round grinding methods as a result of the enhanced heat generation and mechanical interlocking on the interface.

42 ENGINEERING↗
DOE OSTI2021

A Lagrangian perspective on tropical anvil cloud lifecycle in present and future climate

The evolution of tropical anvil clouds from their origin in deep convective cores to their slow decay determines the climatic effects of clouds in tropical convective regions. Despite the relevance of anvil clouds for climate and responses of clouds to global warming, processes dominating their evolution are not well understood. Currently available observational data reveal instantaneous snapshots of anvil cloud properties, but cannot provide a process-based perspective on anvil evolution. We therefore conduct simulations with the high resolution version of the Exascale Earth System Model in which we track mesoscale convective systems over the Tropical Western Pacific and compute trajectories that follow air parcels detrained from peaks of convective activity. With this approach we gain new insight into the anvil cloud evolution both in present day and future climate. Comparison with geostationary satellite data shows that the model is able to simulate maritime mesoscale convective systems reasonably well. Trajectory results indicate that anvil cloud lifetime is about 15 hours with no significant change in a warmer climate. The anvil ice mixing ratio is larger in a warmer climate due to a larger source of ice by detrainment and larger depositional growth leading to a more negative net cloud radiative effect along detrained trajectories. However, the increases in sources are counteracted by increases in sinks of ice, particularly snow formation and sedimentation. Furthermore, we find that the mean anvil cloud feedback along trajectories is positive and consistent with results from more traditional cloud feedback calculation methods.

climate, weather, mesocale cellular convection↗
DOE OSTI2021

Selective dopant segregation modulates mesoscale reaction kinetics in layered transition metal oxide

Incorporation of foreign elements into the cathode material is broadly adopted by both academia and industry to improve the battery performance. The lack of an in-depth understanding for the underlying mechanism, however, makes it a largely try-and-error process with unsatisfactory efficiency and effectiveness. This is particularly true for the electrochemical reaction kinetics that is heterogeneous over a broad range of length scales and is determined collectively by the cathode’s electronic structure, lattice configuration, and micro-morphology. In this work, we unveiled a facet-dependent dopant segregation effect in Zr-modified single-crystal LiNi 0.6 Co 0.2 Mn 0.2 O 2 cathode. By forming kinetically favored corners on the cathode particles, the presence of a trace amount of Zr critically modulates the mesoscale reaction kinetics. Our findings suggest that a delicately controlled dopant distribution is a viable strategy for designing the next-generation battery cathode with superior structural and chemical robustness.

36 MATERIALS SCIENCE↗