Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamic profiling”

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 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
DOE OSTI2012

Investigation of Surface Treatments to Improve the Friction and Wear of Titanium Alloys for Diesel Engine Components

This final report summarizes experimental and analytical work performed under an agreement between the U.S. Department of Energy, Office of Energy Efficiency and Renewable Energy, Office of Transportation Technologies, and UT-Battelle LLC. The project was directed by Jerry Gibbs, of the U.S. Department of Energy’s Propulsion Materials Program, with management by D. P. Stinton and J. A. Haynes of ORNL. Participants included Peter J. Blau (Principal Investigator), Kevin M. Cooley (senior technician), Melanie J. Kirkham (materials scientist) of the Materials Science and Technology Division or ORNL, and Dinesh G. Bansal, a post doctoral fellow employed by Oak Ridge Associated Universities (ORAU) and who, at the time of this writing, is an engineer with Cummins, Inc. This report covers a three-year effort that involved two stages. In the first stage, and after a review of the literature and discussions with surface treatment experts, a series of candidate alloys and surface treatments for titanium alloy (Ti-6Al-4V) was selected for initial screening. After pre-screening using an ASTM standard test method, the more promising surface treatments were tested in Phase 2 using a variable loading apparatus that was designed and built to simulate the changing load patterns in a typical connecting rod bearing. Information on load profiles from the literature was supplemented with the help of T.C. Chen and Howard Savage of Cummins, Inc. Considering the dynamic and evolving nature of materials technology, this report presents a snapshot of commercial and experimental bearing surface technologies for titanium alloys that were available during the period of this work. Undoubtedly, further improvements in surface engineering methods for titanium will evolve.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2024

Enabling "Fuel" Switching through On-Demand Wind Control

This work demonstrates that classical shear-flow stability theory can be successfully applied to modify wind turbine wakes and also explains the success of several emerging, empirically-arrived control methods (i.e., dynamic induction and helix control). Linear stability theory predictions are extended to include the effects of non-axisymmetric inflow profiles, such as wind shear, which is shown to not strongly affect the primary forcing frequency. The predictions, as well as idealized large-eddy simulations using actuator-line representation of the turbine blades, agree that the n = 0 and ±1 modes have faster initial growth rates than higher-order modes, suggesting the lower-order modes are more appropriate for wake control. Exciting the lower-order modes with periodic pitching of the blades produces higher entrainment into the wake and consequently faster wake recovery.

17 WIND ENERGY↗
DOE OSTI2023

Design and Characterization of a Lens-Coupled System for Dynamic X-Ray Diffraction

X-ray diffraction (XRD) is a necessary technique for understanding states of materials under static and dynamic loading conditions. The higher-pressure Equation of State (EOS) of many materials can only be explored via shock or ramp compression at temperatures and pressures of interest. While static XRD work has yielded EOS measurements in the 100 - 200 GPa regime, dynamic X-ray diffraction (DXRD) can explore EOS phases in the TPa regime, which closely resembles inner-core planetary conditions. DXRD hinges on the ability to measure the exact phase or phase change of a material while under dynamic loading conditions. Macroscopic diagnostic systems (e.g. velocimetry and pyrometry) can infer a phase change but not identify the specific phase entered by a material. While microscopic (atomic-level) diagnostic systems (e.g. DXRD) have been designed and implemented in Department of Energy’s (DOE) National Laboratories complex, the unique nature of Sandia National Laboratories’ Pulsed Power Facility (Z Machine) prohibits the use of such devices. The destructive nature of Z experiments presents a challenge to data capture and retrieval. Furthermore there are electromagnetic interference, X-ray background, and mechanical constraints to consider. Thus, a multi-part X-ray diagnostic for use on the Z Machine and Z-Beamlet Laser system has been designed and analyzed. Portions of this new DYnamic SCintillator Optic (DYSCO) have been built, tested and fielded. A data analysis software has been written. Finally, the radiance profile of the DYSCO’s scintillator has been characterized through experiments performed at the University of Arizona.

36 MATERIALS SCIENCE↗
DOE OSTI2016

THE FIRST SPECTRUM OF THE COLDEST BROWN DWARF

The recently discovered brown dwarf WISE 0855 presents the first opportunity to directly study an object outside the solar system that is nearly as cold as our own gas giant planets. However, the traditional methodology for characterizing brown dwarfs—near-infrared spectroscopy—is not currently feasible, as WISE 0855 is too cold and faint. To characterize this frozen extrasolar world we obtained a 4.5–5.2 μ m spectrum, the same bandpass long used to study Jupiter’s deep thermal emission. Our spectrum reveals the presence of atmospheric water vapor and clouds, with an absorption profile that is strikingly similar to Jupiter’s. The spectrum quality is high enough to allow for the investigation of dynamical and chemical processes that have long been studied in Jupiter’s atmosphere, but now on an extrasolar world.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

Simulation of Hollow Electron Beam Collimation in the Fermilab Tevatron Collider

The concept of augmenting the conventional collimation system of high-energy storage rings with a hollow electron beam was successfully demonstrated in experiments at the Tevatron. A reliable numerical model is required for understanding particle dynamics in the presence of a hollow beam collimator. Several models were developed to describe imperfections of the electron beam profile and alignment. The features of the imperfections are estimated from electron beam profile measurements. Numerical simulations of halo removal rates are compared with experimental data taken at the Tevatron.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Experimental Investigation of Vapor Formation in Liquid CO2 Flow Through a Converging-Diverging Nozzle

Carbon dioxide is an attractive working fluid for many cycles, including for pumped thermal energy storage (PTES). A challenge with some proposed sCO2 PTES cycles is the operation of sCO2 machinery outside the typical bounds of experience, with local phase change from the liquid state being particularly unknown. Presently, there is insufficient data in the literature regarding multiphase CO2 to adequately design a multiphase-tolerant turbine, so generation of foundational data is required. This experimental study investigates the flow characteristics of sub-sonic liquid CO2 undergoing expansion and phase change in a converging-diverging nozzle. The nozzle is instrumented to measure static pressure, unsteady pressure, temperature, and density. The static pressure transducers are located at 27 axial locations to accurately characterize the pressure profile in the nozzle. High-accuracy RTDs are located at the entrance and exit of the nozzle, and three dynamic pressure transducers are strategically located to capture any unsteady phenomena. During testing, values of mass flow and nozzle inlet pressure are swept to vary the pressure drop and fluid properties. The measured total pressure drop in the nozzle is compared to a homogenous model and the Lockhart-Martinelli correlation method, with the latter predicting loss quite closely. The resulting data set is valuable for validating multiphase numerical models in a simple geometry before implementation of these models in turbomachinery design.

25 ENERGY STORAGE↗
DOE OSTI2014

Tropical temperature trends in Atmospheric General Circulation Model simulations and the impact of uncertainties in observed SSTs

Abstract The comparison of trends in various climate indices in observations and models is of fundamental importance for judging the credibility of climate projections. Tropical tropospheric temperature trends have attracted particular attention as this comparison may suggest a model deficiency. One can think of this problem as composed of two parts: one focused on tropical surface temperature trends and the associated issues related to forcing, feedbacks, and ocean heat uptake and a second part focusing on connections between surface and tropospheric temperatures and the vertical profile of trends in temperature. Here we focus on the atmospheric component of the problem. We show that two ensembles of Geophysical Fluid Dynamics Laboratory HiRAM model runs (similar results are shown for National Center for Atmospheric Research's CAM4 model) with different commonly used prescribed sea surface temperatures (SSTs), namely, the HadISST1 and “Hurrell” data sets, have a difference in upper tropical tropospheric temperature trends (∼0.1 K/decade at 300 hPa for the period 1984–2008) that is about a factor 3 larger than expected from moist adiabatic scaling of the tropical average SST trend difference. We show that this surprisingly large discrepancy in temperature trends is a consequence of SST trend differences being largest in regions of deep convection. Further, trends, and the degree of agreement with observations, not only depend on SST data set and the particular atmospheric temperature data set but also on the period chosen for comparison. Due to the large impact on atmospheric temperatures, these systematic uncertainties in SSTs need to be resolved before the fidelity of climate models' tropical temperature trend profiles can be assessed.

Flannaghan, T. J.↗
DOE OSTI2018

Dynamic Power Distribution System Management With a Locally Connected Communication Network

Coordinated optimization and control of distribution-level assets can enable a reliable and optimal integration of massive amount of distributed energy resources (DERs) and facilitate distribution system management (DSM). Accordingly, the objective is to coordinate the power injection at the DERs to maintain certain quantities across the network, e.g., voltage magnitude, line flows, or line losses, to be close to a desired profile. By and large, the performance of the DSM algorithms has been challenged by two factors: i) the possibly non-strongly connected communication network over DERs that hinders the coordination; ii) the dynamics of the real system caused by the DERs with heterogeneous capabilities, time-varying operating conditions, and real-time measurement mismatches. In this paper, we investigate the modeling and algorithm design and analysis with the consideration of these two factors. In particular, a game theoretic characterization is first proposed to account for a locally connected communication network over DERs, along with the analysis of the existence and uniqueness of the Nash equilibrium (NE) therein. To achieve the equilibrium in a distributed fashion, a projected-gradient-based asynchronous DSM algorithm is then advocated. The algorithm performance, including the convergence speed and the tracking error, is analytically guaranteed under the dynamic setting. In conclusion, extensive numerical tests on both synthetic and realistic cases corroborate the analytical results derived.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2016

THE HANLE AND ZEEMAN POLARIZATION SIGNALS OF THE SOLAR Ca II 8542 Å LINE

We highlight the main results of a three-dimensional (3D) multilevel radiative transfer investigation about the solar disk-center polarization of the Ca ii 8542 Å line. First, through the use of a 3D model of the solar atmosphere, we investigate the linear polarization that occurs due to the atomic level polarization produced by the absorption and scattering of anisotropic radiation, taking into account the symmetry-breaking effects caused by its thermal, dynamic, and magnetic structure. Second, we study the contribution of the Zeeman effect to the linear and circular polarization. Finally, we show examples of the Stokes profiles produced by the joint action of the atomic level polarization and the Hanle and Zeeman effects. We find that the Zeeman effect tends to dominate the linear polarization signals only in the localized patches of opposite magnetic polarity, where the magnetic field is relatively strong and slightly inclined; outside such very localized patches, the linear polarization is often dominated by the contribution of atomic level polarization. We demonstrate that a correct modeling of this last contribution requires taking into account the symmetry-breaking effects caused by the thermal, dynamic, and magnetic structure of the solar atmosphere, and that in the 3D model used the Hanle effect in forward-scattering geometry (disk-center observation) mainly reduces the polarization corresponding to the zero-field case. We emphasize that, in general, a reliable modeling of the linear polarization in the Ca ii 8542 Å line requires taking into account the joint action of atomic level polarization and the Hanle and Zeeman effects.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Multiphasic size-dependent growth dynamics of nanoparticle ensembles

Colloidal nanoparticles are of great interest in modern science and industry. However, the thermodynamic mechanism and dynamics of nanoparticle growth have yet to be understood. Addressing these issues, we tracked hundreds of in-situ growth trajectories of a nanoparticle ensemble using liquid-phase TEM and found that the nanoparticle growth, including coalescence, exhibits nanoparticle size-dependent multiphasic dynamics, unexplainable by current theories. Motivated by this finding, we developed a model and theory for an ensemble of growing nanoparticles, providing a unified, quantitative understanding of the time-dependent mean and fluctuation of nanoparticle size and size-dependent growth rate profiles across various nanoparticle systems and experimental conditions. Our work reveals that the chemical potential in a small nanoparticle strongly deviates from the Gibbs-Thomson equation, shedding light on how it governs the size-dependent growth dynamics of nanoparticles.

Kim, Ji-Hyun↗
DOE OSTI2015

Kiloparsec-scale simulations of star formation in disk galaxies. III. Structure and dynamics of filaments and clumps in giant molecular clouds

We present hydrodynamic simulations of self-gravitating dense gas in a galactic disk, exploring scales ranging from 1 kpc down to ∼0.1 pc. Our primary goal is to understand how dense filaments form in giant molecular clouds (GMCs). These structures, often observed as infrared dark clouds (IRDCs) in the Galactic plane, are thought to be the precursors to massive stars and star clusters, so their formation may be the rate-limiting step controlling global star formation rates in galactic systems as described by the Kennicutt–Schmidt relation. Our study follows on from Van Loo et al., which carried out simulations to 0.5 pc resolution and examined global aspects of the formation of dense gas clumps and the resulting star formation rate. Here, using our higher resolution, we examine the detailed structural, kinematic, and dynamical properties of dense filaments and clumps, including mass surface density (Σ) probability distribution functions, filament mass per unit length and its dispersion, lateral Σ profiles, filament fragmentation, filament velocity gradients and infall, and degree of filament and clump virialization. Where possible, these properties are compared to observations of IRDCs. By many metrics, especially too large mass fractions of high Σ>1 g cm{sup −2} material, too high mass per unit length dispersion due to dense clump formation, too high velocity gradients, and too high velocity dispersion for a given mass per unit length, the simulated filaments differ from observed IRDCs. We thus conclude that IRDCs do not form from global fast collapse of GMCs. Rather, we expect that IRDC formation and collapse are slowed significantly by the influence of dynamically important magnetic fields, which may thus play a crucial role in regulating galactic star formation rates.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

A profile monitor for proton radiography experiments at the Los Alamos Neutron Science Center

The Proton Radiography (pRad) facility at the Los Alamos Neutron Science Center utilizes pulses of protons delivered by the 800 MeV linear accelerator to produce a series of radiographic images to study the dynamic behavior of materials under extreme conditions. Radiographs taken with an empty field of view, or beam pictures, are used to normalize transmission. However, because the center of the proton beam shifts between pulses, an in situ method for measuring beam position is required to normalize images for beam movement to perform absolute radiography. The beam profile monitor described here uses an array of scintillating fibers positioned in the beam path to produce light proportional to beam intensity across the beam cross section. This light is detected using fast photodiodes and a digital oscilloscope, providing a response time of several nanoseconds—suitable for measuring the 50-ns proton pulses used in pRad. The profile monitor achieves a measured position precision of 40 μm and an intensity precision of 0.7%, allowing for beam movement corrections to be applied to images, thereby improving data accuracy and image quality.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2019

DynamicME: dynamic simulation and refinement of integrated models of metabolism and protein expression

Background Genome-scale models of metabolism and macromolecular expression (ME models) enable systems-level computation of proteome allocation coupled to metabolic phenotype. Results We develop DynamicME, an algorithm enabling time-course simulation of cell metabolism and protein expression. DynamicME correctly predicted the substrate utilization hierarchy on a mixed carbon substrate medium. We also found good agreement between predicted and measured time-course expression profiles. ME models involve considerably more parameters than metabolic models (M models). We thus generate an ensemble of models (each model having its rate constants perturbed), and then analyze the models by identifying archetypal time-course metabolite concentration profiles. Furthermore, we use a metaheuristic optimization method to calibrate ME model parameters using time-course measurements such as from a (fed-) batch culture. Finally, we show that constraints on protein concentration dynamics (“inertia”) alter the metabolic response to environmental fluctuations, including increased substrate-level phosphorylation and lowered oxidative phosphorylation. Conclusions Overall, DynamicME provides a novel method for understanding proteome allocation and metabolism under complex and transient environments, and to utilize time-course cell culture data for model-based interpretation or model refinement.

Mathematical & Computational Biology↗
DOE OSTI2015

Development of an x-ray Talbot-Lau moiré deflectometer for fast density profile measurements of dense plasmas generated by beam-target interactions

We report that for the first time an x-ray Talbot-Lau moiré deflectometer is being developed that will use a flash tube source and fast detector for dynamic density gradient measurements. In Talbot-Lau moiré deflectometry, an x-ray grating makes an image of itself on a second grating (the Talbot effect) to produce a moiré pattern on a detector. The test object is placed between these gratings, with variations in index of refraction changing the pattern. A third grating in front of an incoherent x-ray source produces an array of coherent sources. With a 150 kV x-ray flash tube as the source, the gratings are placed in a glancing angle setup for performance at ~60 keV. The detector is a gated CCD with a fast scintillator for x-ray conversion. This diagnostic, designed for the Dual-Axis Radiographic Hydrodynamic Test facility (DARHT) at Los Alamos National Laboratory, measures the density profile of dense plasma plumes ejected from beam-target interactions. DARHT has two high-current, pulsed, inductive linear electron accelerators with bremsstrahlung targets at the end of each beam line to create 2-D radiographic images of hydrodynamic tests. One multi-pulse accelerator has up to four beam pulses striking the same target within 2 μs. Computer simulations that model target evolution and ejected material between pulses are used to design these targets for optimal radiographic performance; the x-ray deflectometer will directly measure density gradients in the ejected plumes and provide the first experimental constraints to these models. During the first year, currently underway, the diagnostic systems are being designed. In year two, the flash tube and fast detector will be deployed at DARHT for radiographic imaging while the deflectometer is built and tested on the bench with a continuous source. In year three, the fast deflectometer will be installed on DARHT and density measurements will be performed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Cusps and cores in the presence of galactic bulges

In this paper, we study how the presence of bulge formation in galaxies influence their inner density profile, by means of an extended version of the Del Popolo (2009) semi-analytical model. As in Del Popolo (2009), the model takes into account the effect of baryons adiabatic contraction, ordered and random angular momentum, dynamical friction, and adds to the previous the effect of gas cooling, star formation, supernova feedback, and reionization. Our model shows that dwarf galaxies are bulgeless, in agreement with observations showing that the large majority of them has no stellar bulges, and are characterized by a flat profile well described by a Burkert profile. We then studied the effect of a bulge, added to the cored DM halo, on the density profile. In the case of a galaxy having a mass 10{sup 11}M{sub ⊙} the inner density profile has a slope α ≅ 0.65, for a bulge of 4.5 × 10{sup 9}M{sub ⊙} , while if bulge formation is not considered, the slope would be α ≅ 0.55. If the bulge is larger, 6.5 × 10{sup 9}M{sub ⊙} the slope is α ≅ 0.7. In the case of a larger galaxy with 10{sup 12}M{sub ⊙} the slope is α ≅ 0.85, while in absence of bulge it is α ≅ 0.75. We finally study how the inner slope α changes with the bulge mass, and we find a correlation among the two quantities. The result shows that bulge formation has an important role in shaping the inner DM density profile in agreement with Inoue and Saitoh (2011) result. The result implies that going from Sc to SO Hubble type the slope is slightly steepening due to the bulge formation, and due to the fact that early type galaxies have larger bulges.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

RADIO MONITORING OF THE TIDAL DISRUPTION EVENT SWIFT J164449.3+573451. I. JET ENERGETICS AND THE PRISTINE PARSEC-SCALE ENVIRONMENT OF A SUPERMASSIVE BLACK HOLE

We present continued radio observations of the tidal disruption event Swift J164449.3+573451 extending to {delta}t Almost-Equal-To 216 days after discovery. The data were obtained with the EVLA, AMI Large Array, CARMA, the SMA, and the VLBA+Effelsberg as part of a long-term program to monitor the expansion and energy scale of the relativistic outflow, and to trace the parsec-scale environment around a previously dormant supermassive black hole (SMBH). The new observations reveal a significant change in the radio evolution starting at {delta}t Almost-Equal-To 1 month, with a brightening at all frequencies that requires an increase in the energy by about an order of magnitude, and an overall density profile around the SMBH of {rho}{proportional_to}r{sup -3/2} (0.1-1.2 pc) with a significant flattening at r Almost-Equal-To 0.4-0.6 pc. The increase in energy cannot be explained with continuous injection from an L{proportional_to}t{sup -5/3} tail, which is observed in the X-rays. Instead, we conclude that the relativistic jet was launched with a wide range of Lorentz factors, obeying E(> {Gamma}{sub j}){proportional_to}{Gamma}{sup -2.5}{sub j}. The similar ratios of duration to dynamical timescale for Sw 1644+57 and gamma-ray bursts (GRBs) suggest that this result may be applicable to GRB jets as well. The radial density profile may be indicative of Bondi accretion, with the inferred flattening at r {approx} 0.5 pc in good agreement with the Bondi radius for a {approx}few Multiplication-Sign 10{sup 6} M{sub Sun} black hole. The density at {approx}0.5 pc is about a factor of 30 times lower than inferred for the Milky Way Galactic Center, potentially due to a smaller number of mass-shedding massive stars. From our latest observations ({delta}t Almost-Equal-To 216 days) we find that the jet energy is E{sub j,iso} Almost-Equal-To 5 Multiplication-Sign 10{sup 53} erg (E{sub j} Almost-Equal-To 2.4 Multiplication-Sign 10{sup 51} erg for {theta}{sub j} = 0.1), the radius is r Almost-Equal-To 1.2 pc, the Lorentz factor is {Gamma}{sub j} Almost-Equal-To 2.2, the ambient density is n Almost-Equal-To 0.2 cm{sup -3}, and the projected angular size is r{sub proj} Almost-Equal-To 25 {mu}as, below the resolution of the VLBA+Effelsberg. Assuming no future changes in the observed evolution and a final integrated total energy of E{sub j} Almost-Equal-To 10{sup 52} erg, we predict that the radio emission from Sw 1644+57 should be detectable with the EVLA for several decades and will be resolvable with very long baseline interferometry in a few years.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

THE OUTSKIRTS OF GLOBULAR CLUSTERS AS MODIFIED GRAVITY PROBES

In the context of theories of gravity modified to account for the observed dynamics of galactic systems without the need to invoke the existence of dark matter, a prediction often appears regarding low-acceleration systems: wherever a falls below a{sub 0}, one should expect a transition from the classical to the modified gravity regime. This modified gravity regime will be characterized by equilibrium velocities that become independent of distance and that scale with the fourth root of the total baryonic mass, V{sup 4}{proportional_to}M. The two conditions above are the well-known flat rotation curves and Tully-Fisher relations of the galactic regime. Recently, however, a similar phenomenology has been hinted at, at the outskirts of Galactic globular clusters, precisely in the region where a < a{sub 0}. Radial profiles of the projected velocity dispersion have been observed to stop decreasing along Keplerian expectations and to level off at constant values beyond the radii where a < a{sub 0}. We have constructed gravitational equilibrium dynamical models for a number of globular clusters for which the above gravitational anomaly has been reported, using a modified Newtonian force law that yields equilibrium velocities equivalent to modified Newtonian dynamics. We find models having an inner Newtonian region and an outer modified gravity regime, which reproduce all observational constraints, surface brightness profiles, total masses, and line-of-sight velocity dispersion profiles, can be easily constructed. Through the use of detailed single stellar population models tuned individually to each of the globular clusters in question, we derive estimates of the total masses for these systems. Interestingly, we find that the asymptotic values of the velocity dispersion profiles are consistent with scaling with the fourth root of the total masses, as expected under modified gravity scenarios.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Spall fracture in additive manufactured tantalum

In this work, we present a series of experiments on the response of additive manufactured (AM) tantalum to dynamic loading, specifically the spall strength. Rectangular plates of AM tantalum were produced, with subsequent characterization revealing a highly anisotropic microstructure. Samples were taken from these plates to investigate the effect of anisotropy on the spall strength: the resistance to high strain-rate tensile damage. A conventional, wrought tantalum sample, possessing an equiaxed microstructure, was also tested to serve as a control. Additionally, shock loading was performed via light gas-gun flyer-plate impact experiments, with laser velocimetry on the rear of the samples to record the shock wave profiles and soft-recovery techniques to allow post-mortem analysis. In general, the AM samples were found to have a higher Hugoniot elastic limit, the dynamic yield strength under shock loading, while having a reduced spall strength, when compared to the wrought tantalum samples.

36 MATERIALS SCIENCE↗