Search NASA⌕ Search

SEARCH · Search NASA

Results for “conditional approach”

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 91 records · Page 5

Micro Baselines for Operational Technology Environments

Critical infrastructure stakeholders need to baseline their networks to understand expected communications. Top-down approaches to baselining rely on observables that are generally available but lack properties upon which traditional statistical tools depend. We propose to construct micro-baselines: signatures within operational networks based on observables associated with specific events. Such observables are informed by precursor analysis reports of historical cyber attacks on operational environments developed by Cybersecurity for Operational Technology Environments (CyOTE). Baseline measurements depend upon context beyond the cyber domain. An energy plant's baseline running in the summer may statistically differ from a similar facility in a colder region. Domain knowledge must be integrated to apply general micro-baselining algorithms to a facility-specific context. Therefore, we propose to explore the feasibility of transferring micro baselining algorithms across different facilities. Facilities that implement the same processes in different geographic locations will be compared relative to observable measurements used in micro-baselining for comparable events. One evaluation approach would condition or augment dynamic observables measured within a facility network testbed with additional observables derived from geographic context or infrastructure dependencies such as those provided by the All-Hazards Analysis tool.

97 MATHEMATICS AND COMPUTING↗

Surrogates for Valve-Controlled Pipe Flow: Accelerating Nuclear Reactor Design

Neural surrogate models are developed to replace expensive steady-state RANS CFD simulations for valve-controlled pipe flow in nuclear reactor design. Using parametric CFD data generated with MOOSE Pronghorn across a range of valve geometry and flow conditions, three approaches are compared: a POD-based reduced-order model, a structured UNet on a cylindrical grid, and unstructured models (DeepONet and BiStride MeshGraphNet) on nondimensionalized point clouds. POD achieves the highest accuracy (99%) with fast inference but requires storing all solution snapshots, while the DeepONet and BSMS-GNN both achieve ~89% accuracy at sub-second inference, with the BSMS-GNN offering superior geometric generalizability. These surrogates enable rapid ranking of candidate valve designs and can warm-start CFD solvers to accelerate convergence, supporting agentic design iteration on the Prometheus platform.

42 - ENGINEERING↗

First Phase Consensus Roadmap for Development of Condition-Based Cable Reliability Assurance

The objective of this work was to develop a first phase consensus roadmap for condition-based qualification (CBQ) of electrical cables. With CBQ, qualification of Class 1E electrical cables moves from a time-based approach to a condition-based approach, which is anticipated to be safer in terms of reliability and conservatism, and more cost effective in the long run. However, due to barriers, the CBQ approach has not yet been adopted by U.S. nuclear power plants (NPPs). Based upon a review of current work evaluating CBQ, the limitation of available condition monitoring technology seems to be the largest barrier. The importance of condition monitoring, or more specifically selecting appropriate condition indicators, during CBQ cannot be understated. However, selecting appropriate condition indicators is challenged by techniques that are destructive and only evaluate cable degradation locally. Further, arguably, no one identified condition indicator fully establishes cable condition. Thus, additional work is necessary to evaluate potential condition indicators towards CBQ. In addition to the requirements of IEC/IEEE Std. 60780-323, ideal condition indicators should include a) both destructive and non-destructive approaches, b) both local and global measurements, c) real-time (i.e., online) monitoring that trends with degradation, d) enable correlation with qualified levels of degradation, and e) be established within a repository of condition indicators with applicable materials and/or components and their acceptance criteria. Additional work is needed in development of technology and methodology prior to adoption of CBQ, especially for extending qualified life of installed components. Education and early experience by the industry and regulators will be required for this change in approach as an alternative to re-analysis. A series of workshops that bring together stakeholders to identify and address gaps will be needed. The longstanding cooperative working group of cable researchers from the U.S. Department of Energy, the Electric Power Research Institute, and the Nuclear Regulatory Commission forms a valuable starting point for development of a consensus roadmap to condition-based qualification approach as a viable options for qualification of cable systems in U.S. light water reactors.

42 ENGINEERING↗

Polyethylene Upcycling to Liquid Alkanes in Molten Salts under Neat and External Hydrogen Source-Free Conditions

Development of facile approaches to convert plastic waste into liquid fuels under neat conditions is highly desired but challenging, particularly without noble metal catalysts and an external hydrogen source. Herein, highly efficient and selective polyethylene-to-gasoline oil (branched C 6 –C 12 alkanes) conversion was achieved under mild conditions (<170 °C) using commercially available AlCl 3 -containing molten salts as reaction media and to provide catalytic sites (no extra solvents, additives, or hydrogen feeding). The high catalytic efficiency and selectivity was ensured by the abundant active Al sites with strong Lewis acidity (comparable to the Al type in acidic zeolite) and highly ionic nature of the molten salts to stabilize the carbenium intermediates. Dynamic genesis of the Al sites was elucidated via time-resolved Al K-edge soft X-ray and 27 Al NMR, confirming the tricoordinated Al 3+ as active sites and its coordination with the as-generated alkene/aromatic intermediates. Further, the carbenium formation and polyethylene chain variation was illustrated by inelastic neutron scattering (INS) and an isotope-labeling experiment. Theoretical simulations further demonstrated the successive hydride abstraction, β-scission, isomerization, and internal hydrogen transfer reaction pathway with AlCl 3 as active sites. This facile catalytic system can further achieve the conversion of robust, densely assembled, and high molecular weight plastic model compounds to liquid alkane products in the diesel range.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Local practically safe extremum seeking with assignable rate of attractivity to the safe set

We present Assignably Safe Extremum Seeking (ASfES), an algorithm designed to minimize a measured, static objective function while maintaining a measured, static metric of safety (a control barrier function or CBF) to be positive in a practical sense. We ensure that for trajectories with safe initial conditions, the violation of safety can be made arbitrarily small through appropriately chosen design constants. We also guarantee an assignable “attractivity” rate: from unsafe initial conditions, the trajectories approach the safe set, in the sense of the measured CBF, at a rate no slower than a user-assigned rate. Similarly, from safe initial conditions, the trajectories approach the unsafe set, in the sense of the CBF, no faster than the assigned attractivity rate. The feature of assignable attractivity is not present in the semiglobal version of safe extremum seeking, where the semiglobality of convergence is achieved by slowing the adaptation. We also demonstrate local convergence of the parameter to a neighborhood of the minimum of a quadratic objective function constrained to the safe set with a linear CBF. The ASfES algorithm and analysis are multivariable, but we also extend the algorithm to a Newton-Based ASfES scheme which we show is only useful in the scalar case. The proven properties of the designs are illustrated through simulation examples.

42 ENGINEERING↗

Water heating energy use reductions from EPA WaterSense lavatory plumbing fittings

Hot water savings from water-efficient lavatory fittings lead to reductions in water heating energy consumption, and ultimately to decreases in carbon emissions. This paper characterizes existing and proposed approaches used to estimate hot water savings and carbon emissions reductions stemming from the U.S. Environmental Protection Agency's WaterSense program. Also described are refinements that improve the accuracy of residential hot water use percentage estimates of lavatory fittings. The authors conclude that (1) hot water percentages for showers and faucets calculated using up-to-date, publicly available national data are consistent with those found by regional studies and household-level models of water use; (2) the accuracy of heating energy savings estimates attributable to WaterSense-labeled lavatory products, as well as associated emissions reductions, can be refined by modifying the existing energy factor/uniform energy factor (EF/UEF)-based estimation approach with available data. The refined approach accounts for more nuanced conditions than the EF/UED-based approach but depends on data not always available; and (3) the approaches described and intermediate outputs can be generalized for other water conservation programs or estimating purposes.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Stable parallel training of Wasserstein conditional generative adversarial neural networks

In this work, we propose a stable, parallel approach to train Wasserstein conditional generative adversarial neural networks (W-CGANs) under the constraint of a fixed computational budget. Differently from previous distributed GANs training techniques, our approach avoids inter-process communications, reduces the risk of mode collapse and enhances scalability by using multiple generators, each one of them concurrently trained on a single data label. The use of the Wasserstein metric also reduces the risk of cycling by stabilizing the training of each generator. We illustrate the approach on the CIFAR10, CIFAR100, and ImageNet1k datasets, three standard benchmark image datasets, maintaining the original resolution of the images for each dataset. Performance is assessed in terms of scalability and final accuracy within a limited fixed computational time and computational resources. To measure accuracy, we use the inception score, the Fréchet inception distance, and image quality. An improvement in inception score and Fréchet inception distance is shown in comparison to previous results obtained by performing the parallel approach on deep convolutional conditional generative adversarial neural networks as well as an improvement of image quality of the new images created by the GANs approach. Weak scaling is attained on both datasets using up to 2000 NVIDIA V100 GPUs on the OLCF supercomputer Summit.

97 MATHEMATICS AND COMPUTING↗

WEC fault modelling and condition monitoring: A graph-theoretic approach

The nature of wave resources usually requires wave energy converter (WEC) components to handle peak loads (i.e., torques, forces, and powers) that are many times greater than their average loads, accelerating equipment degradation. Moreover, due to their isolated nature and harsh operating environment, WEC systems are projected to possess high operations and maintenance (O&M) cost, i.e., around 27% of their leveled cost of energy. As such, developing techniques to mitigate these costs through the application of condition monitoring and fault tolerant control will significantly impact the economic feasibility of grid connected WEC power. Toward this goal, models of faulty components are developed in the open source modeling platform, WEC-Sim, to estimate the performance and measurable states of a WEC operating with likely device and sensor failures. Two types of faulty component models are then applied to a point absorber WEC model with basic controller damping and spring forces. Resulting changes in device behavior are recorded as a benchmark, and a graph-theoretic approach is proposed for fault detection and identification utilizing multivariate time series. Simulation results demonstrate that these faults can greatly affect the WEC performance, and that the proposed method can effectively detect and classify different types of faults.

16 TIDAL AND WAVE POWER↗

High-precision inversion of dynamic radiography using hydrodynamic features

While radiography is routinely used to probe complex, evolving density fields in research areas ranging from materials science to shock physics to inertial confinement fusion and other national security applications, complications resulting from noise, scatter, complex beam dynamics, etc. prevent current methods of reconstructing density from being accurate enough to identify the underlying physics with sufficient confidence. In this work, we show that using only features that are robustly identifiable in radiographs and combining them with the underlying hydrodynamic equations of motion using a machine learning approach of a conditional generative adversarial network (cGAN) provides a new and effective approach to determine density fields from a dynamic sequence of radiographs. In particular, we demonstrate the ability of this method to outperform a traditional, direct radiograph to density reconstruction in the presence of scatter, even when relatively small amounts of scatter are present. Our experiments on synthetic data show that the approach can produce high quality, robust reconstructions. We also show that the distance (in feature space) between a testing radiograph and the training set can serve as a diagnostic of the accuracy of the reconstruction.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Systems and methods for randomized, packet-based power management of conditionally-controlled loads and bi-directional distributed energy storage systems

The present disclosure provides a distributed and anonymous approach to demand response of an electricity system. The approach conceptualizes energy consumption and production of distributed-energy resources (DERs) via discrete energy packets that are coordinated by a cyber computing entity that grants or denies energy packet requests from the DERs. The approach leverages a condition of a DER, which is particularly useful for (1) thermostatically-controlled loads, (2) non-thermostatic conditionally-controlled loads, and (3) bi-directional distributed energy storage systems. In a first aspect of the present approach, each DER independently requests the authority to switch on for a fixed amount of time (i.e., packet duration). The coordinator determines whether to grant or deny each request based electric grid and/or energy or power market conditions. In a second aspect, bi-directional DERs, such as distributed-energy storage systems (DESSs) are further able to request to supply energy to the grid.

Frolik, Jeff↗

Physical discovery in representation learning via conditioning on prior knowledge

Recent advances in electron, scanning probe, optical, and chemical imaging and spectroscopy yield bespoke data sets containing the information of structure and functionality of complex systems. In many cases, the resulting data sets are underpinned by low-dimensional simple representations encoding the factors of variability within the data. The representation learning methods seek to discover these factors of variability, ideally further connecting them with relevant physical mechanisms. However, generally, the task of identifying the latent variables corresponding to actual physical mechanisms is extremely complex. Here, we present an empirical study of an approach based on conditioning the data on the known (continuous) physical parameters and systematically compare it with the previously introduced approach based on the invariant variational autoencoders. The conditional variational autoencoder (cVAE) approach does not rely on the existence of the invariant transforms and hence allows for much greater flexibility and applicability. Interestingly, cVAE allows for limited extrapolation outside of the original domain of the conditional variable. However, this extrapolation is limited compared to the cases when true physical mechanisms are known, and the physical factor of variability can be disentangled in full. We further show that introducing the known conditioning results in the simplification of the latent distribution if the conditioning vector is correlated with the factor of variability in the data, thus allowing us to separate relevant physical factors. We initially demonstrate this approach using 1D and 2D examples on a synthetic data set and then extend it to the analysis of experimental data on ferroelectric domain dynamics visualized via piezoresponse force microscopy.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Optimizing the fuel efficiency of an opposed piston engine for electric power generation

This paper investigates the optimal crankshaft motion for an opposed piston (OP) engine in a novel hybrid architecture to maximize fuel efficiency. The OP engine was selected for this work due to its inherent thermodynamic benefits and the balanced nature of the engine which can achieve downsizing through reducing the number of cylinders rather than the individual cylinder volume. The typical geartrain required on an OP engine was exchanged for two electric motors, reducing friction loss and decoupling the crankshafts. Using the motors to control the crankshaft motion profiles, this architecture introduces capabilities to dynamically vary compression ratio, combustion volume, and scavenging dynamics. To leverage these opportunities, an optimization scheme was developed utilizing nonlinear optimization of a 0-D model to compute the crankshaft motion profile that maximizes the work generated by the system. This optimization was then iteratively coupled with a high fidelity model which supplies the cylinder flow boundary conditions. This iterative approach reduces the model complexity used in the optimal control problem (OCP) while capturing the gas exchange dynamics critical to the 2-stroke cycle of the OP engine. By using the rate of change of motor torque as the input to the OCP, the torque fluctuation in a single cycle can be limited to ensure tracking feasibility. The results show crankshaft velocity slows during the compression stroke and conversely accelerates during the expansion stroke, reducing the peak motor torque required for control and thus reducing the motor losses. The extended residence time at top dead center, however, leads to an increase in heat transfer, illustrating the trade-off between the work extraction efficiency and the indicated engine efficiency.

Engineering↗

Role of Bifunctional Ru/Acid Catalysts in the Selective Hydrocracking of Polyethylene and Polypropylene Waste to Liquid Hydrocarbons

Hydrogenolysis of C–C bonds over Ru-based catalysts has emerged as a deconstruction strategy to convert single-use polyolefin waste to liquid alkanes at relatively mild conditions, but this approach exhibits limitations, including methane formation resulting from terminal C–C bond scission. In this study, a variety of catalysts were investigated for the reductive deconstruction of polyethylene (PE) and polypropylene (PP) to identify supports that promote nonterminal C–C bond scission. We found that Ru nanoparticles supported on Brønsted-acidic zeolites with the faujasite (FAU) and Beta (BEA) topologies were highly active for the cleavage of C–C bonds in PE and PP, exhibiting improved liquid yields and suppressed methane formation. For the deconstruction of PE, supporting ruthenium nanoparticles (5 wt %) on FAU increased the yields of liquid alkanes to 67% compared to 33% over an inert silica support (5 wt % Ru/SiO 2 ) at 200 °C, 16 h, under 30 bar of H2. A dramatic selectivity enhancement toward liquid hydrocarbons was also observed for PP over Ru/FAU and Ru/BEA compared to Ru/SiO 2 . To understand the origin of this selectivity improvement, a combination of ex situ and operando characterization techniques were used to reveal that both catalyst structure and acidity play key roles in PE and PP conversion. Operando X-ray absorption spectroscopy studies with model polyolefins over Ru-supported catalysts with varying acidity levels revealed that the local chemical environment of Ru [0] during the reaction is consistent across multiple acidic supports, although the onset of reduction during synthesis of the nanoparticles varies across different supports. These results, combined with reactivity data, demonstrate the importance of the acid-noble metal cooperativity in promoting selective C–C bond scission toward liquid alkanes that shifts the mechanism from hydrogenolysis to ideal hydrocracking.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Estimation of the Surface Fluxes for Heat and Momentum in Unstable Conditions with Machine Learning and Similarity Approaches for the LAFE Data Set

Abstract Measurements of three flux towers operated during the land atmosphere feedback experiment (LAFE) are used to investigate relationships between surface fluxes and variables of the land–atmosphere system. We study these relations by means of two machine learning (ML) techniques: multilayer perceptrons (MLP) and extreme gradient boosting (XGB). We compare their flux derivation performance with Monin–Obukhov similarity theory (MOST) and a similarity relationship using the bulk Richardson number (BRN). The ML approaches outperform MOST and BRN. Best agreement with the observations is achieved for the friction velocity. For the sensible heat flux and even more so for the latent heat flux, MOST and BRN deviate from the observations while MLP and XGB yield more accurate predictions. Using MOST and BRN for latent heat flux, the root mean square errors (RMSE) are 107 Wm $$^{-2}$$ - 2 and 121 Wm $$^{-2}$$ - 2 , respectively, as well as the intercepts of the regression lines are $$\approx 110$$ ≈ 110 Wm $$^{-2}$$ - 2 . For the ML methods, the RMSEs reduce to 31 Wm $$^{-2}$$ - 2 for MLP and 33 Wm $$^{-2}$$ - 2 for XGB as well as the intercepts to just 4 Wm $$^{-2}$$ - 2 for MLP and $$-1$$ - 1 Wm $$^{-2}$$ - 2 for XGB with slopes of the regression lines close to 1, respectively. These results indicate significant deficiencies of MOST and BRN, particularly for the derivation of the latent heat flux. In fact, in contrast to the established theories, feature importance weighting demonstrates that the ML methods base their improved derivations on net radiation, the incoming and outgoing shortwave radiations, the air temperature gradient, and the available water contents, but not on the water vapor gradient. The results imply that further studies of surface fluxes and other turbulent variables with ML techniques provide great promise for deriving advanced flux parameterizations and their implementation in land–atmosphere system models.

54 ENVIRONMENTAL SCIENCES↗

Scanning Electrochemical Microscopy for Kinetic Investigations in Viscous Deep Eutectic Solvents: Identifying Practical Approach Curves and Deviations from Electron Transfer Models

Determining heterogeneous electrochemical electron transfer (ET) kinetics in electrolytes with a wide range of physical properties is of great interest for achieving high-performance redox flow batteries. Among such electrolytes, concentrated hydrogen-bonded electrolytes (CoHBEs), including deep eutectic solvents (DESs), have recently garnered significant attention. Unfortunately, traditional Tafel analysis using macroelectrodes often encounters issues with mass transfer limitations in CoHBEs with high viscosities, thereby restricting kinetic analysis to a narrow potential window. Here, in this work, we introduce a methodology for evaluating ET kinetics in viscous DES using the scanning electrochemical microscopy (SECM). We first determined practical solutions to SECM tip positioning in ethaline DES, which yield pseudopositive feedback responses. Lattice Boltzmann method (LBM) simulations helped us rationalize the impact of the fluid and concentration fields, as well as tip geometry, tip approach velocity v, and the solvent viscosity ηs, on the shape of the approach curves. In addition to successfully recreating approach curves over a variety of conditions, we found that approaching a conductor ensured a practical point where the normalized tip response (Ni T = 2) converged at L = 0.7 within ∼10% error regardless of tip velocity. With positioning capabilities at hand, we investigated the kinetics of Fe 3+ /Fe 2+ redox couple in aqueous and the ethaline media. The experimental kinetic results were interpreted using the Butler–Volmer (BV) and Marcus–Hush–Chidsey (MHC) models. For ethaline, a nonideal kinetic behavior was observed, potentially attributed to solvent dynamics within DESs or to the interplay of chloride anions in the charge transfer process.

electrodes↗

The importance of capturing power system operational details in resource adequacy assessments

Traditional methods for assessing the resource adequacy (RA) of a power system are becoming obsolete due to emerging trends such as the increasing deployment of variable renewable energy and storage. Consequently, analysts are recommending that RA be assessed using a Monte Carlo simulation approach that models chronological power system operations over many instances of possible operating conditions. However, this approach is necessarily more complex and computationally demanding, which is an obstacle to real-world implementation. Here, in this study, we investigate which operational details of power systems are important to capture in order to accurately evaluate a system's RA, versus details that add complexity but do not meaningfully affect RA results. To do so, we develop a probabilistic RA assessment framework by adapting an existing production cost model and apply it to a case study based on the IEEE Reliability Test System. Our results indicate that multi-year data, storage dispatch, and transmission limits are key details to incorporate. Accurate RA results can be obtained using non-economic dispatch strategies as long as they are coordinated with detailed operational strategies. We also demonstrate how popular expectation-based RA metrics can mask important differences in the characteristics of loss of load events.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Deep Vadose Zone Monitoring Test Bed (FY23 Status Report)

A subsurface air injection at the Hanford Site’s Deep Vadose Zone Monitoring Test Bed was completed to realize a change in subsurface hydrologic conditions in accordance with a soil desiccation remedy. The injection mimicked a previous injection at the site that was performed in accordance with a vadose zone treatability study. Unlike the previous test which relied on electrical methods only, the change in hydrologic conditions during the recent test was also monitored using cross-hole seismic sensing methods to assess the ability of the seismic methods to evaluate changes in moisture conditions of desiccated sediment at the field scale. Data from in situ neutron probes indicates a reduction of soil moisture in the vicinity of the injection well due to the air injection. Similar changes were observed in the time-lapse electrical resistivity and seismic data, which indicates a loss of soil moisture over time. Tomographic inversions of the time-lapse geophysical data illustrate the 2D and 3D features of the soil moisture distribution over time. Time-lapse electrical resistivity tomography (ERT) results show a reduction in the electrical conductivity of the subsurface in the vicinity of the injection well, with most changes occurring within the screened interval. Similar patterns are observed in the seismic tomography results, with both methods illustrating two lobe-shaped features of reduced soil moisture. Use of seismic and ERT technologies in tandem takes advantage of two complementary geophysical monitoring technologies, providing increased sensitivity to specific hydrologic conditions. The multiphysics approach, therefore, has the potential to improve the ability to estimate subsurface moisture conditions from sensor-based and remotely sensed geophysical data that will ultimately improve the ability of remediation contractors to evaluate remedy performance.

58 GEOSCIENCES↗

Microstructure and CPO evolution of dynamically recrystallized olivine during complex deformation conditions: a full-field numerical modeling approach

The rheological properties of mantle rocks are strongly dependent on their crystallographic preferred orientation (CPO). Olivine CPO, defined by the orientation of seismically fast [100] axes parallel to flow direction, is also thought to be a dominant contributor to seismic anisotropy in the Earth's upper mantle. However, the amount of deformation needed to overprint a new CPO on a pre-existing fabric and the impact of the inherited CPOs on the transient microstructure evolution, remain unknown. This study employs a full-field numerical approach (VPFFT-ELLE) to explore the dynamic recrystallization and microstructural evolution of olivine polycrystalline aggregates under complex deformation conditions. We test four combinations of successive pure shear and simple shear boundary conditions. Findings indicate that inherited CPOs influence subsequent deformation in a manner dependent on the kinematic relationship between successive stages. In all cases, a minor strain increment (ε ∼0.3–0.6) is sufficient to erase the previous microstructure and CPO. However, when deformation conditions change dramatically (e.g., stretching direction changes orthogonally), the intensity of the new CPO developed is significantly lower and strain distribution are specially altered. During a transient strain stage, pre-existing microstructures undergo extensive reworking, especially when deformation conditions are changed dramatically, such as switching from simple shear to a pure shear condition with a parallel shortening direction relative to the stretching direction. We estimate the significance of these results in interpreting observations of seismic velocity anisotropy, concluding that P-wave seismic anisotropy is significantly and positively correlated with the evolution of olivine CPO with deformation history. This research underscores the transient nature of microstructural rearrangement in olivine aggregates and the necessity for caution in interpreting seismic anisotropy in regions with complex deformation histories, as inherited CPOs can influence current fabric development and induce deviation to the present deformation conditions.

36 MATERIALS SCIENCE↗