Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “solutes”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

An Analytic Benchmark for Neutron Boltzmann Transport with Downscattering—Part I: Flux and Eigenvalue Solutions

Computing in the energy dimension is one of the greatest challenges confronting present-day deterministic neutron transport solvers. Accurately resolving the neutron flux as neutrons downscatter across resonances in the nuclear cross sections currently requires considerable computing power and suffers from approximation errors. Flux uncertainty resulting from the uncertainty of the resonance structure is the single-largest cause of reactivity uncertainty. Any additional reference solution for the critical neutron downscattering problem with resonance phenomena would be a boon to verification and validation of neutronics codes. This paper establishes a benchmark to verify the accuracy of neutron transport criticality solvers along the energy dimension. For the first time, the analytic solution of the flux amplitude is derived in the particular case of an infinite homogeneous medium with isotropic scattering in the center of mass and an arbitrary number of no-threshold, neutral particle reaction resonances (e.g., radiative capture, fission, and resonance scattering). Furthermore, original analytic expressions are established to quantify the discrepancy between the ψ k (E) and ψ α (E) flux amplitudes, respective solutions of the multiplication factor k, or the exponential time-evolution frequency α eigenproblems. The physical study of these relations led to analysis of their first-order relative difference near the criticality condition α=0. Finally, numerical solutions are provided to a benchmark problem constituted of the first resonance of 239 Pu, the 6.67-eV resonance of 238 U, and a scattering isotope with a flat cross section, allowing for the computational verification of the energy resolution of current neutron transport criticality codes. Through these novel results, this analytic benchmark can serve as a reference to verify the energy resolution and sensitivity analysis of neutron transport criticality calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Demonstrating Computational Equivalence Between Continuous and Discrete Adjoint Methods by Calculating Time-Dependent Adjoint Solutions with Neutron Diffusion Models

The continuous adjoint method and the discrete adjoint method are two alternative approaches used to calculate adjoint solutions for adjoint systems. The continuous adjoint method derives adjoint equations analytically from continuous forward equations and then solves the adjoint equations either analytically or numerically in a discretized form whereas the discrete adjoint method calculates the adjoint solutions directly from the discretized forward equations. With regard to the methodology development and calculation procedure, distinct differences are well recognized between the two methods. For certain reasons, both methods are exclusively preferred and commonly used by different computational communities, but limited studies clarify the connections between the two adjoint methods from either of the communities. Herein, this paper demonstrates the computational equivalence between the continuous and discrete adjoint methods by investigating time-dependent adjoint solutions to the two-group neutron diffusion model in nuclear reactor analysis problems using both methods. Adjoint solutions can be used to estimate system parameters for reactor safety analysis. Appropriate final state conditions for the adjoint systems are specified in both of the methods, and the conditions are clarified with proper physical explanations. With the help of an event-based case study on neutron diffusion models, the accuracy of the time-dependent adjoint fluxes obtained from both methods is verified, and the pros and cons of both adjoint methods are examined. More importantly, the computational equivalence of both methods is demonstrated when they are applied to multigroup neutron diffusion systems. The advantage of calculating time-dependent adjoint fluxes by directly solving time-dependent adjoint systems rather than taking steady-state approximations as in common practice is also demonstrated.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Resurgent trans-series for generalized Hastings–McLeod solutions

Here, we show that the physical Hastings–McLeod solution of the integrable Painlevé II equation generalizes in a natural way to a class of non-integrable equations, in a way that preserves many of the significant qualitative properties. The Hastings–McLeod solution of Painlevé II is an important and universal example of resurgent relations between perturbative and non-perturbative physics. We derive the trans-series structure of the generalized Hastings–McLeod solutions, demonstrating that integrability is not essential for the resurgent asymptotic properties of the solutions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Magnetic frame-dragging correction to the electromagnetic solution of a compact neutron star

ABSTRACT Neutron stars are usually modelled as spherical, rotating perfect conductors with a predominant intrinsic dipolar magnetic field anchored to their stellar crust. Due to their compactness, General Relativity corrections must be accounted for in Maxwell’s equations, leading to modified interior and exterior electromagnetic solutions. We present analytical solutions for slowly rotating magnetized neutron stars, taking into account the magnetic frame-dragging correction. For typical compactness values, i.e. Rs ∼ 0.5 [R*], we show that the new terms lead to a per cent order correction in the magnetic field orientation and strength compared to the case, with no magnetic frame-dragging correction. Also, we obtain a self-consistent redistribution of the surface azimuthal current. We verify the validity of the derived solution through two-dimensional particle-in-cell simulations of an isolated neutron star. Defining the azimuthal electric and magnetic field amplitudes during the transient phase as observables, we prove that the magnetic frame-dragging correction reduces the transient wave amplitude, as expected from the analytical solution. We show that simulations are more accurate and stable, when we include all first-order terms. The increased accuracy at lower spatiotemporal resolutions translates into a reduction in simulation runtimes.

Torres, R. (ORCID:0000000291820228)↗
DOE OSTI2023

Exact solutions for the electromagnetic fields of a flying focus

The intensity peak of a “flying” focus travels at a programmable velocity over many Rayleigh ranges while maintaining a near-constant profile. Assessing the extent to which these features can enhance laser-based applications requires an accurate description of the electromagnetic fields. Here we present exact analytical solutions to Maxwell’s equations for the electromagnetic fields of a constant-velocity flying focus, generalized for arbitrary polarization and orbital angular momentum. Furthermore, the approach combines the complex source-point method, which transforms multipole solutions into beam-like solutions, with the Lorentz invariance of Maxwell’s equations. Propagating the fields backward in space reveals the space-time profile that an optical assembly must produce to realize these fields in the laboratory. Comparisons with simpler paraxial solutions provide conditions for their reliable use when modeling a flying focus.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Testing meson portal dark sector solutions to the MiniBooNE anomaly at the Coherent CAPTAIN Mills experiment

A solution to the MiniBooNE excess invoking rare three-body decays of the charged pions and kaons to new states in the MeV mass scale was recently proposed as a dark-sector explanation. This class of solution illuminates the fact that, while the charged pions were focused in the target-mode run, their decay products were isotropically suppressed in the beam-dump-mode run in which no excess was observed. This suggests a new physics solution correlated to the mesonic sector. We investigate an extended set of phenomenological models that can explain the MiniBooNE excess as a dark sector solution, utilizing long-lived particles that might be produced in the three-body decays of the charged mesons and the two-body anomalous decays of the neutral mesons. Over a broad set of interactions with the long-lived particles, we show that these scenarios can be compatible with constraints from LSND, KARMEN, and MicroBooNE, and evaluate the sensitivity of the ongoing and future data taken by the Coherent CAPTAIN Mills experiment to a potential discovery in this parameter space. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Perturbative stability of non-Abelian electric field solutions

We consider SU(2) gauge theory with a scalar field in the fundamental representation. The model is known to contain electric field solutions sourced by the scalar field that are distinct from embedded Maxwell electric fields. We examine the perturbative stability of the solution and identify a region of parameter space where the solution is stable. In the regime where the scalar field has a negative mass squared, the solution has two branches, and we identify an instability in one of the branches.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Learning Optimal Power Flow Solutions using Linearized Models in Power Distribution Systems

Solving nonlinear optimal power flow (OPF) problem is computationally expensive, and poses scalability challenges for power distribution networks. An alternative to solving the original nonlinear OPF is the linear approximated OPF models. Although, these linear approximated OPF models are fast, the resulting solutions may result in significant optimality gap. Lately, the application of machine learning (ML) methods in successfully solving the nonlinear OPF has been reported. These methods learn and estimate the nonlinear control policies using a purely data-driven approach. In this paper, we propose an approach to complements the ML based approach to solving OPF using solutions from known linearized OPF model. Specifically, we use supervised learning to map the solutions of linear OPF to nonlinear control variables. Unlike, the traditional ML based methods for OPF that approximate the full distribution feeder model using function approximation, our approach uses a two-node approximation of radial networks. The proposed approach is validated using IEEE 123 bus test system for OPF solutions obtained using the nonlinear OPF models.

optimal power flow, power distribution systems, su↗
DOE OSTI2023

A MultiGPU Performance-Portable Solution for Array Programming Based on Kokkos

Today, multiGPU nodes are widely used in high-performance computing and data centers. However, current programming models do not provide simple, transparent, and portable support for automatically targeting multiple GPUs within a node on application areas of array programming. In this paper, we describe a new application programming interface based on the Kokkos programming model to enable array computation on multiple GPUs in a transparent and portable way across both NVIDIA and AMD GPUs. We implement different variations of this technique to accommodate the exchange of stencils (array boundaries) among different GPU memory spaces, and we provide autotuning to select the proper number of GPUs, depending on the computational cost of the operations to be computed on arrays, that is completely transparent to the programmer. We evaluate our multiGPU extension on Summit (#5 TOP500), with six NVIDIA V100 Volta GPUs per node, and Crusher that contains identical hardware/software as Frontier (#1 TOP500), with four AMD MI250X GPUs, each with 2 Graphics Compute Dies (GCDs)for a total of 8 GCDs per node. We also compare the performance of this solution against the use of MPI + Kokkos, which is the cur-rent de facto solution for multiple GPUs in Kokkos. Our evaluation shows that the new Kokkos solution provides good scalability for many GPUs and a faster and simpler solution (from a programming productivity perspective) than MPI + Kokkos.

Valero Lara, Pedro↗
DOE OSTI2023

Application of the Solute-Solvent EMF Cell to Measure Activity of NiF 2 in Molten FLiNaK

Electromotive force (emf) measurements made using a combination of solute- and solvent-based electrodes were used to determine the activity of NiF 2 in molten FLiNaK eutectic at 823 K across a concentration range of x NiF2 = 5.2 × 10 –4 –1.0 × 10 –2 . The solute emf values were measured using electrodes consisting of Ni wires immersed in FLiNaK with dissolved NiF 2 contained in graphite crucibles. The measured emf values were then converted to the FLiNaK Eutectic Potassium Electrode (FEKE) potential and used to quantify the activity of dissolved NiF 2 . This quantification was based upon comparative measurements of a reference solvent electrode consisting of a K-Bi alloy immersed in pure FLiNaK contained in a boron nitride crucible and a solute electrode. Short cell lives were characteristic of the measurements due to the corrosive nature of the fluoride salts. Quantifying the activity of NiF 2 will improve the utility of Ni 2+ /Ni reference electrodes in molten fluoride salts, which are notoriously difficult electrolytes to work with because of their reactivity. In conclusion, this work demonstrates the general nature of the solute-solvent approach as a repeatable, easily employed method for measuring the activity values of electroactive species in a variety of molten salts to improve understanding of the electroactive species behavior in these systems.

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

Denudation, solute export, landscape evolution modeling, and geographic information system data for the East River watershed, Colorado, USA (2020-2024)

This data package contains geographic information system (GIS) layers and tabular datasets associated with the study of lithologic controls on denudation, solute export, carbon-scaling relationships, and transient landscape evolution in the East River watershed near Crested Butte, Colorado, USA. The package includes GIS layers used to produce the Figure 2 map, including drainage, hillshade, lithology, sample locations, and basin polygons, together with comma-separated value (CSV) tables and matching CSV data dictionaries. One group of tables reports sample-level and catchment-level information for river-sediment samples analyzed for in situ-produced cosmogenic beryllium-10 (10Be), including sample names, outlet elevations, geographic coordinates, upstream drainage area, rock-type classes, production-rate scaling scheme, analyzed nuclide, catchment-averaged denudation rates, and associated lower and upper analytical uncertainties. Sample and catchment attributes provide the basis for comparing denudation rates across intrusive, shale, sedimentary, and mixed-lithology settings. A second group of tables reports supporting information for landscape-evolution modeling and the mapped geologic framework of the study area. Included files list parameter values and definitions for the two-phase landscape-evolution simulations, summarize full-domain model erosion fluxes and topographic metrics for different simulation configurations, provide a fixed-area carbon-model scaling table, and summarize mapped geologic units within the East River study domain, including geologic code, formation name, lithologic description, mapped area, and lithologic class grouping. Model outputs and geologic summaries support interpretation of transient landscape behavior and its relation to the mapped distribution of shale, intrusive, sedimentary, and surficial units. A third group of tables reports hydrologic and hydrochemical information used to quantify dissolved export from the watershed. Included files provide site-level values for drainage area, mean annual solute export, standard error of annual export, area-normalized solute yield, and equivalent weathering rate for five East River monitoring sites, along with metadata describing the number, sampling cadence, and date range of discharge records and partial and full total dissolved solids observations used in the solute-yield analyses. The package also contains a supplementary daily ion-load time series with daily mean discharge, discharge observation counts, dissolved concentrations, and daily loads for calcium, magnesium, sodium, potassium, chloride, sulfate, nitrate, fluoride, dissolved silica, charge-balance bicarbonate, and total dissolved solids. The package contains GIS files, comma-separated value files (.csv), CSV data dictionaries, a file-level metadata table, a package-tree text file, and a readme text file.

10Be↗
DOE OSTI2020

Benchmark Solutions for Radiation Transport in Stochastic Media with Inhomogeneous Material Statistics

Accurately solving implicit Monte Carlo (IMC) thermal photon transport problems with mixed material cells is important in realistic applications. The production IMC package at LLNL treats mixed material cells arising from ALE remap and hydrodynamics using the same approximate model. The new Imp IMC thermal photon transport package currently under development has both a material interface reconstruction (MIR) algorithm and a Levermore-Pomraning (LP) stochastic medium algorithm for treating mixed material cells. Existing stochastic medium algorithms for treating mixed material cells in IMC lack a complete theoretical basis. The IMC LP algorithm implementation has been demonstrated to reproduce published deterministic LP solutions for the particular case of spatially homogeneous material statistics. Realistic simulations will include spatially inhomogeneous material statistics (material mean chord lengths). In a previous investigation, the LP-model for transport in binary stochastic media in rod geometry was generalized to accommodate spatially varying material chord lengths, i.e., the mixing statistics were allowed to be nonhomogeneous. Analytical solutions were obtained and used to produce a verifi cation suite for the Imp IMC Levermore-Pomraning implementation for different spatial variations of the chord lengths. However, the accuracy of the LP model when the mixing statistics are nonhomogeneous has not been assessed and leaves open the question of whether local accuracy is improved or further degraded when chord lengths are not uniform. This shortcoming is rectifi ed here by developing benchmark analytic solutions for transport in binary Markovian stochastic mixtures in rod geometry with nonhomogeneous mixing statistics, using spatially varying chord lengths considered in the previous investigation based on the LP model. Methods for sampling a nonhomogeneous Poisson process (NHPP) are first described and used to construct individual realizations of the binary mixtures in rod geometry. Analytic solutions are then obtained for the forward and backward directed fluxes on a given realization, now viewed as a deterministic medium with alternating layers of the two materials with known interface locations. Finally, material averaged scalar fluxes are obtained using these sampling schemes with spatially linear and quadratic chord lengths and used to assess the accuracy of the previously obtained LP-model results.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Code and Solution Verification Assessment of the CTF Thermal Hydraulic Subchannel Code

CTF is a thermal-hydraulics subchannel code jointly developed by Oak Ridge National Laboratory and North Carolina State University. Over the past seven years, the Consortium for Advanced Simulation of Light Water Reactors (CASL) has made a significant investment in developing CTF so it can be used to model light water reactors, including nominal operating conditions, departure from nucleate boiling analysis, and transients ranging from loss of flow to reactivity insertion accidents. In addition to implementing new modeling capabilities and developing the user input and output interface, extensive work has been performed to improve the code’s quality assurance program, resulting in a development process that conforms with NQA-1 requirements. The CASL program follows the Predictive Capability Maturity Model (PCMM) approach for assessing code quality, which emphasizes performing code verification(ensuring the code converges to the correct answer) and solution verification (ensuring the code converges for the intended application). Code and solution verification are used to identify uncertainty errors introduced by numerical approximations in the code and are important for demonstrating that the model is coded without error, which is an important aspect of the Best Estimate plus Uncertainty method. This paper presents a comprehensive overview of the code and solution verification testing that has been performed on CTF. A top-down approach is taken in which the intended CTF applications are presented, followed by the code features required for their modeling. These features are then linked to the applicable code and solution verification tests that demonstrate proper functioning. Past testing efforts are summarized, and new tests are added to help close gaps in the presented test matrix. Rather than performing “one-off” exercises, these tests are added to the automated CTF regression test suite to ensure continual code quality.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Federated Architecture for Secure and Transactive Distributed Energy Resource Management Solutions (FAST-DERMS)

This document provides system-level specifications for a federated architecture for secure and transactive distributed energy resource management solutions (FAST-DERMS), presents a solution, and describes operational concepts for the proposed solution. FAST-DERMS enables the provision of reliable, resilient, and secure transmission and distribution (T&D) grid services through the scalable aggregation and near-real-time management of utility-scale and small-scale distributed energy resources (DERs). We first present the principles and objectives of FAST-DERMS. Then, after discussing important system concepts, we present the specifications for FAST-DERMS and a solution that employs a distributed and federated control methodology in which the DERs connected to a single point of common coupling with the rest of the system, such as individual substations, are optimized coordinately to provide system-level grid services. FAST-DERMS aims to aggregate and coordinate the operations of DERs to support T&D grid operations. The key optimization and control component of this FAST-DERMS reference implementation is a flexible resource scheduler (FRS) that aggregates the DERs within a substation service area. These FRSs operate at the substation level and perform constrained economic dispatch of DERs, either directly or through a transactive market or aggregator, as shown in Figure ES-1. An FRS Coordinator at the distribution system operator (DSO) level aggregates distribution substations operated by FRSs and interfaces with the transmission system operator (TSO) to provide transmission services. FAST-DERMS also allows for the integration of the FRS Coordinator with an existing distribution utility management system that could be employed by the DSO to enhance distribution grid operations.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Non-Equilibrium Plasma Interactions with Biomaterials, Biological Solutions and Tissues

Cold atmospheric pressure plasma discharges offer an abundant source of reactive oxygen and nitrogen species (RONS) at room temperature enabling unique interactions with biomaterials, biological solutions and tissues. These interactions particularly with living matter are presently an important intellectual frontier in plasma science with promising potential applications ranging from human health care to advanced biomaterial processing. Exciting case studies have been reported that illustrate the huge potential of cold atmospheric plasma technology in wound healing and cancer treatment. The interaction of plasma with conducting and dielectric biomaterials such as tissue strongly influences the plasma properties. In turn this changes the impact of the plasma on the biomaterial. Particularly in the case of living matter, liquid based solutions are ubiquitous which complicates interfacial processes. The lack of insight into the underlying mechanisms of the interaction of plasma with wounds and tumors is currently a bottleneck for the further development of the technology and gives rise to many interesting scientific questions. This project was focused on plasma properties and kinetics during plasma-biomaterial interactions. Both DC pulsed and RF driven atmospheric pressure plasma jets, extensively used by the plasma community were studied. The bio-interfaces included hydrogel as a tissue model, (saline) solutions, bacteria and virus. Plasma diagnostics used include Thomson scattering, Rayleigh scattering, Raman scattering, (two-photon absorption) laser induced fluorescence, optical emission spectroscopy, absorption spectroscopy, molecular beam mass spectrometry and fast imaging allowing to determine electron densities and temperatures, ionic species, reactive species including radicals, gas temperatures, gas composition, electric fields and solution components transferred to the gas phase.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Impact of Irradiation and Decay Time on Chemical Speciation in Uranyl Sulfate Solution

A pH-1, 200-g U/L, highly enriched uranyl sulfate solution was neutron-irradiated using an accelerator-driven neutron source (the Universal Neutron Irradiator, UNI) to fission U-235 under conditions similar to those used by SHINE to produce Mo-99 in a subcritical assembly. We completed a series of five uranyl sulfate irradiations to assess the impact of irradiation and decay time on uranium precipitation and iodine speciation in solution. The duration of each irradiation was increased in order, with the shortest irradiation taking 2.9 hours and the longest, 64.5 hours. A total of 9.5E16 U-235 nuclei were fissioned in all irradiations, and a maximum fission power density of 0.26 W/mL was reached in the final irradiation. We observed no evidence of uranyl peroxide or other precipitates, demonstrating that the addition of Fe 2+ to a final Fe concentration of 200 ppm in the target solution was successful at preventing precipitation under representative irradiation conditions. Nonvolatile I - was found to be the dominant (>70%) iodine species in solution, both immediately after irradiation and after days to weeks of decay time.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2022

Options for the Removal of Cesium and Strontium from Acidic Uranyl Sulfate Solutions from Molybdenum-99 Production

A small laboratory-scale column process aimed at removing fission products from prototypical acidic sulfate solution wastes from a prototypical Mo-99 recovery process has been operated in a radioactive materials hood at Savannah River National Laboratory. This equipment has been operated with realistic concentrations of uranium (0.4 wt % 235 U), plutonium ( 239 Pu), neptunium ( 237 Np), non-radioactive surrogates of common fission products (including Cs and Sr) and two promising absorbents AMP-PAN and CST R91290-B. The goal of this work was to demonstrate the removal and concentration of high specific activity isotopes which are expected to dominate the classification of the low-level waste from the production of 99 Mo. Both AMP-PAN and CST R9120-B are attractive absorbents for the removal of 137 Cs from an acidic sulfate waste stream. They do not provide a solution for the removal of 90 Sr from the acidic solutions, but CST should be investigated at a higher pH. Limited removal of transuranic species and almost no removal of uranium was observed in this testing. This could be important due to the relatively large amount of uranium expected to be present in such waste streams. The capacity of these absorbents for 137 Cs will likely be dependent on the Cs concentration in the solutions processed but loading in the range of 0.2 to 5 mg Cs per gram of absorbent appears to be achievable.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Interaction of Concentrated (68%) Aqueous Nitric Acid Solutions with Cheesecloth: A TGA-DSC-MS Study

The exposure of cotton cheesecloth to concentrated aqueous solutions of nitric acid reduced the onset of combustion by 30-40°C (from ~350°C to ~310°C) relative to that of both unexposed cheesecloth and cheesecloth that was rinsed immediately following exposure to the same solution. Cheesecloth exposed to concentrated (68%) nitric acid underwent three exothermic pre-combustion reactions which onset at 125-135°C, 160- 170°C and 200-210°C. These reactions give off mostly water vapor, with lesser amounts of NO and carbon dioxide. Exposure to concentrated nitric acid damages the cellulose polymers, oxidizing them, reducing them in length, and increasing their reactivity in response to thermal stress. Cheesecloth which was rinsed after exposure behaved similarly to unexposed cheesecloth, indicating that NC was not formed on exposure to the concentrated aqueous nitric acid solution. The addition of a 239Pu – 241Am spike in an aqueous nitric – hydrofluoric solution to the cheesecloth reduced the onset temperature of the first exothermic reaction to 126°C and lowered the onset temperature of combustion from ~310°C to ~275°C. When dried on to cellulose substrates in the presence of excess nitrate, metal ions may act as oxidizers, or may catalyze the oxidation of the carbonaceous material, lowering ignition temperatures. No evidence for the formation of nitrocellulose following the exposure of cheesecloth to nitric acid was observed during this study

36 MATERIALS SCIENCE↗