Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Fractional operators”

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 361 records · Page 20

Rheology and engine performance of very low sulfur fuel oil blended with 10% fast pyrolysis and hydrothermal liquefaction oils in a 2-stroke crosshead engine

The performance and emissions for a downscaled single-cylinder 2-stroke crosshead engine were determined for a very low sulfur fuel oil (VLSFO) when blended with 10 wt.% fast pyrolysis (FP) or hydrothermal liquefaction (HTL) bio-intermediates. The FP and HTL oils were derived from biomass and were observed to contain lower molecular weight (MW) hydrocarbons than neat VLSFO (which was evaluated as a baseline comparison). The addition of either biofuel reduced the overall viscosity of the VLSFO. Aging tests at 50, 90, and 120°C showed that the dynamic viscosity of VLSFO increased with exposure time up to two weeks. Similar trends were observed for the FP and HTL blends, but a pronounced spike in viscosity occurred for these fuels during the early period of exposure. None of the viscosity increases exceeded the operational limits of fuel system pumps. Engine performance studies were conducted under low, medium and high load operational settings. The relative performance of the test fuels was highly dependent on operating condition. In general, the engine results for the three test fuels were similar, but modest improvements in brake thermal efficiency and brake specific fuel consumption were observed, which may be attributed to the heightened reactivity of low molecular weight fraction of the FP and HTL oils.

09 BIOMASS FUELS↗

Metaproteomics reveals enzymatic strategies deployed by anaerobic microbiomes to maintain lignocellulose deconstruction at high solids

Economically viable production of cellulosic biofuels requires operation at high solids loadings—on the order of 15 wt%. To this end we characterize Nature’s ability to deconstruct and utilize mid-season switchgrass at increasing solid loadings using an anaerobic methanogenic microbiome. This community exhibits undiminished fractional carbohydrate solubilization at loadings ranging from 30 g/L to 150 g/L. Metaproteomic interrogation reveals marked increases in the abundance of specific carbohydrate-active enzyme classes. Significant enrichment of auxiliary activity family 6 enzymes at higher solids suggests a role for Fenton chemistry. Stress-response proteins accompanying these reactions are similarly upregulated at higher solids, as are β-glucosidases, xylosidases, carbohydrate-debranching, and pectin-acting enzymes—all of which indicate that removal of deconstruction inhibitors is important for observed undiminished solubilization. Our work provides insights into the mechanisms by which natural microbiomes effectively deconstruct and utilize lignocellulose at high solids loadings, informing the future development of defined cultures for efficient bioconversion.

09 BIOMASS FUELS↗

Impact of non‐framework cation mixing on the structure and crystallization behavior of model high‐level waste glasses

Abstract Spinel crystallization is known to be detrimental to the operation of Joule heated ceramic melters during the vitrification of iron‐rich high‐level nuclear wastes (HLW) into borosilicate glasses. The literature on this subject focuses on tackling the problem by developing empirical constraints to design compositions, which limit the fraction of spinels formed in the melter or by developing empirical models to predict the settling behavior of spinels in the melter as a function of the glass composition. While these empirical models can predict the behavior of most of the compositions, they are not failsafe as there are always some compositions, whose behavior is beyond the predictive ability of these models. This can lead to undesirable situations during the vitrification of the nuclear waste, and therefore an in‐depth investigation of the chemo‐structural descriptors controlling the crystallization behavior in these glasses is warranted. Accordingly, the present study aims to understand the impact of non‐framework cation mixing (i.e., Li + /Na + and Ca 2+ /Na + ) on the structure (through Raman spectroscopy and Mössbauer spectroscopy) and crystallization behavior (through XRD, SEM‐EDS, and vibrating sample magnetometry) of iron‐rich model HLW glasses in the system: (mol.%) x M y O–(25− x ) Na 2 O–9.12 B 2 O 3 –6.4 Al 2 O 3 –51.25 SiO 2 –7.22 Fe 2 O 3 –0.38 MnO–0.08 Cr 2 O 3 –0.55NiO (M y O = Li 2 O or CaO).

Balasubramanya, Nikhila↗

On the Flow of a Cement Suspension: The Effects of Nano-Silica and Fly Ash Particles

Additives such as nano-silica and fly ash are widely used in cement and concrete materials to improve the rheology of fresh cement and concrete and the performance of hardened materials and increase the sustainability of the cement and concrete industry by reducing the usage of Portland cement. Therefore, it is important to study the effect of these additives on the rheological behavior of fresh cement. In this paper, we study the pulsating Poiseuille flow of fresh cement in a horizontal pipe by considering two different additives and when they are combined (nano-silica, fly ash, combined nano-silica, and fly ash). To model the fresh cement suspension, we used a modified form of the power-law model to demonstrate the dependency of the cement viscosity on the shear rate and volume fraction of cement and the additive particles. The convection–diffusion equation was used to solve for the volume fraction. After solving the equations in the dimensionless forms, we conducted a parametric study to analyze the effects of nano-silica, fly ash, and combined nano-silica and fly ash additives on the velocity and volume fraction profiles of the cement suspension. According to the parametric study presented here, larger nano-silica content results in lower centerline velocity of the cement suspension and larger non-uniformity of the volume fraction. Compared to nano-silica, fly ash exhibits an opposite effect on the velocity. Larger fly ash content results in higher centerline velocity, while the effect of the fly ash on the volume fraction is not obvious. For cement suspension containing combined nano-silica and fly ash additives, nano-silica plays a dominant role in the flow behavior of the suspension. The findings of the study can help the design and operation of the pulsating flow of fresh cement mortars and concrete in the 3D printing industry.

36 MATERIALS SCIENCE↗

NIRP Core Software Suite v.5.1

SAND2023-06610O The Nuclear Incident Response Program (NIRP) Core Software Suite is a set of code that supports multiple applications. It includes miscellaneous base code for data objects, mathematic equations, and user interface components that are used for developing large applications, including those used for data management. Among the types of applications: • Analyst Manager is used for managing contact information, which is also often included in reports and identifies the owners of calculations. • Radionuclide Viewer provides access to the Dose Coefficient File Package (DCFPAK) radiological data, and it complements the Mixture Manager tool. • Mixture Manager creates and manages radionuclides mixtures that are commonly used in other applications. • High Explosive Manager manages explosives and properties. • Chart Viewer is used for data charts such as meteorology charts and other charts specific to data needs. • Unit Converter converts numeric values to/from different units. This package includes object-oriented base code such as data objects describing the date, time, and location of an event; weather conditions; data types that support predefined ranges of values; a unit library, and others. The package also includes mathematic equations such as determining the distance between two points; erf functions; factorial functions; incidence fractions of probit curves; and others. Other features are customized user interface components for text fields; unit selectors (combo boxes); menus; panels; bug-tracking; and automated checks for software updates. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Collins, Ian↗

The impact of fuel and injection strategy on combustion characteristics, emissions and efficiency in gasoline compression ignition operation

Gasoline compression ignition in diesel engines has been proposed in order to meet increasingly stringent emission regulations without sacrificing efficiency. In this study, a six-cylinder heavy-duty diesel engine was operated in a mixing controlled gasoline compression ignition mode to investigate the influence of fuels and injection strategies on the combustion characteristics, emissions, and thermal efficiencies. Fuels, including ethanol (E), isobutanol (IB), and diisobutylene (DIB), were blended with a gasoline fuel to form E10, E30, IB30, and DIB30 based on volumetric fraction. These four blends along with gasoline formed the five test fuels. With these fuels, three injections strategies were investigated, including late pilot injection, early pilot injection, and port fuel injection/direct injection. The impact of moderate exhaust gas recirculation on nitrogen oxides and soot emissions was examined to determine the most promising fuel/injection strategy for emissions reduction. In addition, first and second law analyses were performed to provide insights into the efficiency, loss, and exergy destruction of the various gasoline fuel blends at low and medium load conditions. Overall, the emission output, thermal efficiency, and combustion performances of the five fuels were found to be similar and their differences are modest under most test conditions. E30 with the port/direct fuel injection strategy obtained the peak brake thermal efficiency (46.9%) and gross indicated thermal efficiency (52.2%) at 14 bar. At 7 bar, the peak brake thermal efficiency (44.2%) and gross indicated thermal efficiency (52.5%) also belonged to E30 with the late pilot injection strategy. The second law analysis revealed that exergy destruction was mostly influenced by the fuel type and not by the injection strategy for the 14 bar operating conditions. However, the opposite trend was observed for the 7 bar case. In addition, port fuel injection in combination with oxygenated fuel blends was found to be helpful to suppress soot emissions for operation with exhaust gas recirculation.

co-optima fuels↗

Going off topics to demix quark and gluon jets in α S extractions

Quantum chromodynamics is the theory of the strong interaction between quarks and gluons; the coupling strength of the interaction, α S , is the least precisely-known of all interactions in nature. An extraction of the strong coupling from the radiation pattern within jets would provide a complementary approach to conventional extractions from jet production rates and hadronic event shapes, and would be a key achievement of jet substructure at the Large Hadron Collider (LHC). Presently, the relative fraction of quark and gluon jets in a sample is the limiting factor in such extractions, as this fraction is degenerate with the value of αS for the most well-understood observables. To overcome this limitation, we apply recently proposed techniques to statistically demix multiple mixtures of jets and obtain purified quark and gluon distributions based on an operational definiton. We illustrate that studying quark and gluon jet substructure separately can significantly improve the sensitivity of such extractions of the strong coupling. We also discuss how using machine learning techniques or infrared- and collinear-unsafe information can improve the demixing performance without the loss of theoretical control. While theoretical research is required to connect the extract topics with the quark and gluon objects in cross section calculations, our study illustrates the potential of demixing to reduce the dominant uncertainty for the α S extraction from jet substructure at the LHC.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

An Atom Probe Tomography Study of the Through Wall Attenuation Effect on Cu-rich Precipitate Formation in a Reactor Pressure Vessel Steel

In this study, high-Cu weld material harvested from an ex service reactor pressure vessel (RPV) steel from Unit 1 of the decommissioned Zion Nuclear Generating Station has been characterized using atom probe tomography. Samples taken from 4 different positions through the thicknesses of the pressure vessel wall from the water-side to the air-side were characterized, along with an unirradiated baseline material. In the baseline material, no precipitates were found and the Cu was observed to be fully in solid solution; however, scanning transmission electron microscopy combined with energy dispersive spectroscopy (STEM-EDS) revealed the presence of ϵ -Cu that form during processing of the material and results in the concomitant decrease of matrix Cu. Following irradiation, a high number density of nano-scale Cu-rich precipitates (CRPs) were observed, uniformly distributed throughout the matrix. The Cu content within the CRPs was found to be ~30-35 at.% regardless of location in the wall. No statistically significant variation in the compositions, mean radius, number density, or volume fraction as a function of location within the wall was observed. The measured matrix Cu level excluding CRPs contribution was found to be ~90 appm higher than the solubility limit suggesting that further nucleation and growth of the CRPs under continued operations would have occurred. These results clearly demonstrate that the neutron energy attenuation has no significant effect on the precipitation kinetics of CRPs regardless of location in the wall in high-Cu RPV steels under irradiation.

36 MATERIALS SCIENCE↗

Ion-Specific Effects on PuO 2 Nanoparticle Aggregation and Dissolution in Concentrated Electrolytes

Hydrolytic PuO 2 nanoparticles (NPs) are a dominant aqueous Pu-bearing phase in high ionic strength nuclear wastes, yet their reactivity in nonideal brines remains poorly constrained. We quantify how electrolyte identity and concentration control PuO 2 NP aggregation and ligand-assisted dissolution in acidic, high salinity solutions (NaCl, NaNO 3 , NaClO 4 , Na 2 SO 4 , Na 2 C 2 O 4 up to 5 M). A multitechnique workflow combining liquid scintillation counting (operationally defined aqueous [Pu]), scattering/electrokinetic measurements (aggregate size and zeta potential), and spectroscopy (UV–vis, XPS, Raman) resolves electrolyte-dependent partitioning between colloidal and molecular Pu species. Weakly coordinating anions (ClO 4 – , Cl – , NO 3 – ) largely preserve the (aggregated) nanoparticulate fraction but generate distinct dissolved species at high concentration, i.e., Pu(IV)–nitrato complexes in NaNO 3 and Pu(VI)–chloro complexes in NaCl. In contrast, stronger ligands substantially perturb PuO 2 NP stability: sulfate promotes partial dissolution to Pu(IV)–sulfate complexes at low concentration but reduces aqueous [Pu] at higher sulfate levels via secondary Pu(IV) sulfate formation, whereas oxalate drives strong dissolution to aqueous Pu–oxalate complexes. Aged NPs show similar trends with reduced aqueous fractions and more dominant aggregation mechanisms. In conclusion, these spectroscopically constrained speciation data provide a foundation for incorporating PuO 2 NP reactivity into thermodynamic and reactive transport models for high salinity waste and brine environments.

Aggregation↗

Search for h c → π + π − J / ψ via ψ ( 3686 ) → π 0 h c

Using ( 2712.4 ± 14.3 ) × 10 6 ψ ( 3686 ) events collected with the BESIII detector operating at the BEPCII collider, we search for the hadronic transition h c → π + π − J / ψ via ψ ( 3686 ) → π 0 h c . No significant signal is observed. We set the most stringent upper limits to date on the branching fractions B ( ψ ( 3686 ) → π 0 h c ) × B ( h c → π + π − J / ψ ) and B ( h c → π + π − J / ψ ) at the 90% confidence level, which are determined to be 6.7 × 10 − 7 and 9.4 × 10 − 4 , respectively. Published by the American Physical Society 2024

Ablikim, M.↗

Opportunities for Solar Industrial Process Heat in the United States

Energy used in the production and processing of materials and products represents a significant fraction of the overall energy footprint of the industrial sector. Solar technologies may provide suitable alternative to existing combustion technologies, but their relevance for industries in the United States has been largely unexplored at the necessary levels of spatial, temporal, and operational detail. To provide the first assessment of the technical opportunities for solar technologies to meet the demand for process heat in the United States, we first develop estimates of energy use by process temperature at the county and hourly resolutions for all manufacturing industries. We then estimate the opportunities for seven solar technology packages—comprised of solar thermal and photovoltaic-connected electrotechnologies—to meet process heat demand given solar resources and available land area by county.

14 SOLAR ENERGY↗

Shading and Energy Impacts for Farmer-Friendly Agrivoltaics Arrays

When designing agrivoltaic solar arrays, farmers’ natural row tilling directions might be reasonable choices for array alignment, especially in cases where the tilling direction helps prevent soil erosion or has other benefits for the farmer’s operations. We examined the light field and energy production effects when installing arrays that have azimuth values that deviate from the typical ideal for three important array types: fixed-tilt, single-axis tracking, and vertical bifacial arrays. On a fractional annual energy production basis, though, all three array types lose very little output when rotated away from their ideal. This offers farmers and solar developers a wider range of system design to possibly align solar array rows with the farmer’s preferred tilling directions.

14 SOLAR ENERGY↗

Observation of the decays 𝐵 + → Σ 𝑐 ⁢(2455) ++ $\bar{Ξ}^-_c$ and 𝐵 0 → Σ 𝑐 ⁢(2455) 0⁢ $\bar{Ξ}^0_𝑐$

We report the first observation of the two-body baryonic decays 𝐵 + → Σ 𝑐 ⁢(2455) ++ $\bar{Ξ}^-_c$ and 𝐵 0 → Σ 𝑐 ⁢(2455) 0⁢ $\bar{Ξ}^0_𝑐$ with significances of 7.3⁢𝜎 and 6.2⁢𝜎, respectively, including statistical and systematic uncertainties. The branching fractions are measured to be ℬ⁡(𝐵 + → Σ 𝑐 ⁢(2455) ++ $\bar{Ξ}^-_c$) = (5.74 ± 1.11 ± 0.4⁢2$^{+2.47}_{−1.53}$) ×10 −4 and ℬ⁡(𝐵 0 → Σ 𝑐 ⁢(2455) 0⁢ $\bar{Ξ}^0_𝑐$) = (4.83 ± 1.12 ± 0.3⁢7$^{+0.72}_{−0.60}$) ×10 −4 . The first and second uncertainties are statistical and systematic, respectively, while the third ones arise from the absolute branching fractions of $\bar{Ξ}^-_c$ or $\bar{Ξ}^0_𝑐$ decays. The data samples used for this analysis have integrated luminosities of 711 fb −1 and 365 fb −1 , and were collected at the ϒ⁡(4⁢𝑆) resonance by the Belle and Belle II detectors operating at the KEKB and SuperKEKB asymmetric-energy 𝑒 + ⁢𝑒 − colliders, respectively.

branching fraction↗

A sub-5 ps jitter time-to-digital converter ASIC with back-gate delay tuning in 22 nm CMOS

We present the design of a Time-to-Digital Converter (TDC) ASIC together with performance characterization results at cryogenic temperature (6-8 K), and at room temperature using emulated Low-Gain Avalanche Detector (LGAD) signals. The TDC design uses a two-step architecture with a ring-oscillator based counter and a Vernier delay line fine TDC. The TDC is implemented in an FD-SOI process, and back-gate tuning is used not only to correct threshold variation due to cryogenic operation, but also as a novel way to tune TDC delay elements with very little overhead. Using this technique, we demonstrate a design with very low power (0.5 mW) and area (0.003 mm 2 ). We present test results using the Fermilab Constant Fraction Discriminator (FCFD) ASIC to produce discriminated signals from an internal charge injection mechanism that mimics the waveform and signal amplitude of minimum ionizing particles impinging on LGAD sensors. We characterize the time precision of the full system and verify that the TDC ASIC contributes a negligible amount to the total system time precision, fully consistent with the expected jitter contribution of less than 5 ps. In conclusion, the results presented here demonstrate the utility of our TDC for applications in physics and quantum communications.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Barium stars as tracers of s -process nucleosynthesis in AGB stars

Barium (Ba) stars help to verify asymptotic giant branch (AGB) star nucleosynthesis models since they experienced pollution from an AGB binary companion and thus their spectra carry the signatures of the slow neutron capture process (s process). For a large number (180) of Ba stars, we searched for AGB stellar models that match the observed abundance patterns. We aim to uncover any systematic deviations of the sample abundances from the predictions of the nucleosynthesis models. We employed three machine learning algorithms as classifiers: a Random Forest method, developed for this work, and the two classifiers used in our previous study. Compared to that work, we also expanded our observational sample with 11 Ba stars available in the supersolar metallicity range. We studied the statistical behaviour of the different s-process elements in the observational sample to investigate if the AGB models systematically under- or overpredict the abundances observed in the Ba stars and show the results in the form of violin plots of the residuals between spectroscopic abundances and model predictions. We inspected the correlations between the observed [Fe/H], the s-process elemental abundances, and the residuals. We employed the [Zr/Fe] and [Nb/Fe] abundances as a thermometer to constrain the operational temperature that rules the production of these elements in the sample stars, assuming a steady-state s process. We also investigated the mass distribution of the identified polluter AGB stars and the behaviour of the δ parameter, which describes the fraction of accreted AGB material relative to the Ba star envelope. We find a significant trend in the residuals that implies an underproduction of the elements just after the first s-process peak (Nb, Mo, and Ru) in the models relative to the observations. This may originate from a neutron-capture process (e.g. the intermediate neutron-capture process, i process) not yet included in the AGB models of metallicity from solar to roughly 1/5 solar, corresponding to the range of the Ba stars. Correlations are found between the residuals of these peculiar elements, suggesting a common origin for the deviations from the models. In addition, there is a weak metallicity dependence of the residuals of these elements. The s-process temperatures derived with the [Zr/Fe] – [Nb/Fe] thermometer have an unrealistic value for the majority of our stars. The most likely explanation is that at least a fraction of these elements are not produced in a steady-state s process, and instead may be due to processes not included in the AGB models. The mass distribution of the identified models confirms that our sample of Ba stars was polluted by low-mass AGB stars (< 4 M ⊙ ). Most of the matching AGB models require low accreted mass, but a few systems with high accreted mass are needed to explain the observations.

79 ASTRONOMY AND ASTROPHYSICS↗

Divertor detachment and heat exhaust mitigation control in KSTAR with tungsten divertor

KSTAR has recently undergone an upgrade to use a new tungsten divertor to run experiments in ITER-relevant scenarios. Even with a high melting point of tungsten, it is important to control the heat flux impinging on tungsten divertor targets to minimize sputtering and contamination of the core plasma. Heat flux on the divertor is often controlled by increasing the degree of detachment of scrape-off layer plasma from the target plates. In this work, we have demonstrated successful divertor detachment and heat exhaust dissipation control experiments using two different methods. The first method uses attachment fraction as a control variable which is estimated using ion saturation current measurements from embedded Langmuir probes in the divertor. The second method uses a novel machine-learning-based surrogate model of 2D UEDGE simulation database, DivControlNN. We demonstrated running inference operation of DivControlNN in realtime to estimate heat flux at the divertor and use it as the control variable in a feedback loop with impurity gas flow. We present interesting insights from these experiments including a systematic approach to tuning controllers and discuss future improvements in the control infrastructure and control variables for future burning plasma experiments.

KSTAR tungsten divertor operations↗

Simulation and Optimization of Volatile Fatty Acid Upgrading Strategies for Sustainable Transportation Fuel Production

The multicomponent nature of and variability of biomass makes chemical upgrading via a single process stream to a single end use infeasible for most feedstocks. A more promising approach is to identify upgrading strategies (encompassing bioprocessing, catalysis, and separations) which valorize varied biomass fractions to distinct products to which each is best suited. A robust procedure to carry out this goal should incorporate comparison of various upgrading procedures as well as suitability of products to a slate of identified end uses. We applied this strategy to one biomass-derived feedstock, volatile fatty acids (VFAs) derived from wet waste arrested anaerobic digestion, by developing a computer program, VFA Upgrading to Liquid Transportation fUels Refinery Estimation (VULTURE) which evaluates VFA catalytic upgrading to liquid transportation fuels. VULTURE considers multiple separations, catalysis (ketonization, hydrogenation), and fuel application options, generating hundreds of candidate scenarios for a given VFA stream, then selects several promising strategies that optimize bio-content of products with properties best suited for target fuel types. We find that VFAs are upgraded most efficiently when separate light alcohol (C3-6) and heavy hydrocarbon or alcohol (C7-13) fractions are targeted to create gasoline and heavy-duty (diesel or jet) fuels or fuel blends. Surrogate property testing of VULTURE-recommended fuels reveals that most predictive models employed are robust, while rigorous process simulation shows that the simple unit operation assumptions used in VULTURE are largely accurate, especially for heavy-duty fuel synthesis. Techno-economic and life-cycle analyses show that VFA-derived fuels are profitable and have dramatically (=57%) lower carbon intensities than fossil analogs.

BIOMASS FUELS↗

Quantifying the effect of hydride microstructure on zirconium alloys embrittlement using image analysis

Since hydride microstructure can significantly influence Zr alloy nuclear fuel cladding’s ductility, a new metric has been developed to quantify hydride microstructure in 2D micrographs and relate it to crack propagation. As cladding failure usually results from a hoop stress, this new metric, called the Radial Hydride Continuous Path (RHCP), is based on quantifying the continuity of brittle hydride particles along the radial direction of the cladding tube. Compared to previous metrics, this approach more closely relates to the propensity of a crack to propagate radially through the cladding tube thickness. The RHCP takes into account hydride length, orientation, and connectivity to choose the optimal path for crack propagation through the cladding thickness. The RHCP can therefore be more closely linked to hydride embrittlement of the Zr alloy material, thus creating a relationship between material structure, properties, and performance. The new definition, along with previously proposed metrics such as the Radial Hydride Fraction (RHF), the Hydride Continuity Coefficient (HCC), and the Radial Hydride Continuity Factor (RHCF), have been implemented and automated in MATLAB. These codes are made available with this publication as supplementary materials. These metrics were verified by comparing their predictions of hydride morphology against expected values in simple cases, and the implementation of the new metric was validated by comparing its predictions with manual measurements of hydride microstructure performed on ImageJ. The RHCP was also validated against experimental measurements of fracture behavior and it was shown to correlate with cladding failure better than previous metrics. The information provided by these metrics will help accurately assess cladding integrity during operation, transportation, and storage.

36 MATERIALS SCIENCE↗