Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Phase transformations”

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 1,099 records · Page 61

In situ Synchrotron X‐ray Metrology Boosted by Automated Data Analysis for Real‐time Monitoring of Cathode Calcination

Abstract Synchrotron X‐ray‐based in situ metrology is advantageous for monitoring the synthesis of battery materials, offering high throughput, high spatial and temporal resolution, and chemical sensitivity. However, the rapid generation of massive data poses a challenge to on‐site, on‐the‐fly analysis needed for real‐time process monitoring. Here, a weighted lagged cross‐correlation (WLCC) similarity approach is presented for automated data analysis, which merges with in situ synchrotron X‐ray diffraction metrology to monitor the calcination process of the archetypal nickel‐based cathode, LiNiO 2 . The WLCC approach, incorporating variables that account for peak shifts and width changes associated with structural transformations, enables rapid extraction of phase progression within 10 seconds from tens of diffraction patterns. Details are captured, from initial precursors to intermediates and the final layered LiNiO 2 , providing information for agile on‐site adjustments during experiments and complementing post hoc diffraction analysis by offering insights into early‐stage phase nucleation and growth. Expanding this data‐powered platform paves the way for real time calcination process monitoring and control, which is pivotal to quality control in battery cathode manufacturing.

36 MATERIALS SCIENCE↗

Computational Methodology to Simulate Pyrometallurgical Processes in a Secondary Lead Furnace

Pyrometallurgical recovery of nonferrous metals involves a combination of thermally intensive transformations during exothermic gas-phase reactions, endothermic decomposition of solid charge, and melting of simpler solids. In the recovery of secondary lead, simultaneous thermal effects in a reverberatory-style furnace cause a melt pool to accumulate at the bottom, with lighter solids (slag) floating above and gaseous products from decomposition of the charge diffusing through the gas–slag interface. Species from oxy-fuel combustion of natural gas, species profiles from smelting reactions, and the formation of a melt pool consisting primarily of lead are simulated via a time-averaged formulation. Predictions of outflow are compared with preset inflow profiles to ensure conservation of mass. Thermal profiles for solid, liquid, and gas phases are presented by species. A novel method is implemented to model the latent heat of fusion using a heterogeneous chemical reaction. The simulation is conducted in Simcenter STAR-CCM+ v. 16.02.009-R8.

Rao, Vivek↗

Ni5Ga3 Catalysts for CO2 Reduction to Methanol: Exploring the Role of Ga Surface Oxidation/reductions on Catalytic Activity

Ad-Ni5Ga3/SiO2 catalyst, which is highly active and stable for thermal CO2 hydrogenation to methanol, was investigated to understand its surface dynamics during reaction conditions. The catalyst was prepared, tested and characterized using a multitude of techniques, including ex-situ XRD, TEM, H2-TPR, CO chemisorption and in-situ ETEM and APXPS. The formation of Ga2O3 upon air exposure was observed and the consequences on catalytic activity were investigated. Results revealed that upon air exposure Ga migrates from the subsurface region to the surface of the nanoparticles forming a Ga-oxide shell surrounding a metallic core. The oxide shell can be reduced completely at high temperatures and the temperature of the reduction activation treatment plays a crucial role on the catalytic activity. This air exposure-reduction procedure results in a redistribution of the Ga species in the surface/subsurface region, promoted by a higher Ga affinity for oxygen. XRD analysis of the post-reaction samples provide evidence for the partial transformation of a Ni5Ga3 intermetallic phase into a Ni-rich Ni3Ga phase. APXPS data indicate that the Ni-rich phase forms in the subsurface region. A Ga-rich phase exists on the surface, where the Ni/Ga ratio reaches 2/3 upon activation in hydrogen. At low reduction temperatures an amorphous Ga2O3 shell remains, indicating a promotion effect on methanol synthesis activity.

Ni-Ga, intermetallic, methanol synthesis, in-situ ↗

Transformative high entropy alloy conquers the strength-ductility paradigm by massive interface strengthening

Recent metastable alloy designs have demonstrated simultaneous attainment of high ultimate tensile strength (UTS) and ductility in high entropy alloys but with low yield strength. In this work, we present new strategy for improving the work hardenability and yield strength (YS) together in Fe 38.5 Mn 20 Co 20 Cr 15 Si 5 Cu 1.5 high entropy alloy (Cu-HEA). Drastic increase in the YS (1.5 GPa) is attributed to the formation of γ/ε, ε/ε (twin type) and ε/ε (plate type) interfaces in the microstructure due to extreme grain refinement whereas high UTS-ductility synergy (2.2 GPa, 15%) is attained by dynamic Hall-Petch hardening across these interfaces (i.e. massive interface strengthening) and transformation induced plasticity in γ phase. Thus, this harmonious combination of YS and UTS-ductility synergy in Cu-HEA outperform all structural materials till date. Therefore, deformation-induced massive interface strengthening is a new, yet cost-effective pathway for synergizing the benefits of advanced steels and high entropy alloys together in a material by conventional processing route.

36 MATERIALS SCIENCE↗

Phonon Spectroscopy in Antimony and Tellurium Oxides

α-Sb 2 O 3 (senarmontite), β-Sb 2 O 3 (valentinite), and α-TeO 2 (paratellurite) are compounds with pronounced stereochemically active Sb and Te lone pairs. The vibrational and lattice properties of each have been previously studied but often lead to incomplete or unreliable results due to modes being inactive in infrared or Raman spectroscopy. Here in this paper, we present a study of the relationship between bonding and lattice dynamics of these compounds. Mössbauer spectroscopy is used to study the structure of Sb in α-Sb 2 O 3 and β-Sb 2 O 3 , whereas the vibrational modes of Sb and Te for each oxide are investigated using nuclear inelastic scattering, and further information on O vibrational modes is obtained using inelastic neutron scattering. Additionally, vibrational frequencies obtained by density functional theory (DFT) calculations are compared with experimental results in order to assess the validity of the utilized functional. Good agreement was found between DFT-calculated and experimental density of phonon states with a 7% scaling factor. The Sb–O–Sb wagging mode of α-Sb 2 O 3 whose frequency was not clear in most previous studies is experimentally observed for the first time at ~340 cm -1 . Softer lattice vibrational modes occur in orthorhombic β-Sb 2 O 3 compared to cubic α-Sb 2 O 3 , indicating that the antimony bonds are weakened upon transforming from the molecular α phase to the layer-chained β structure. The resulting vibrational entropy increase of 0.45 ± 0.1 kB/Sb 2 O 3 at 880 K accounts for about half of the α–β transition entropy. The comparison of experimental and theoretical approaches presented here provides a detailed picture of the lattice dynamics in these oxides beyond the zone center and shows that the accuracy of DFT is sufficient for future calculations of similar material structures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Skyrmions Near Defects

We study the impact of an exchange-reducing defect on a skyrmion for a thin film of finite thickness. Attraction of the skyrmion to a defect is demonstrated in a lattice model by computing the micromagnetic energy accounting for the exchange, Dzyaloshinskii-Moriya interaction, magnetic anisotropy, and dipole-dipole coupling. The spiraling dynamics of the skyrmion towards the defect is illustrated by solving numerically the full Landau-Lifshitz-Gilbert equations on a lattice and, independently, the Thiele equation, with the two methods in agreement with each other. We find that the skyrmion deforms when approaching the defect. Here, deformation is small in the parameter space far from the phase boundary that determines stability of skyrmions. It increases dramatically near the phase boundary, leading to the transformation of the skyrmion by the defect into a snake-like magnetic domain.

36 MATERIALS SCIENCE↗

Modeling the Effects of Microwave Heating on Densification in Chemical Vapor Infiltration

Microwave heating has great potential to accelerate the synthesis of ceramic matrix composites (CMCs) by a process called Chemical Vapor Infiltration (CVI). In CVI, reactive gases ingress a porous preform and undergo chemical transformation to deposit solid ceramic phase within the pores at high temperature thus, densifying the preform. However, the competing effects of chemical kinetics and gas transport are known to result in non-uniform densification as the outer surfaces of the preform experience faster depositions compared to the core. Achieving spatial temperature control plays a key role improving the densification quality. Microwave heating can potentially create temperature inversion such that the core of the preform is hotter than the outer surface and subsequently, lead to improved densification while keeping the manufacturing times and costs low. In the present work, a computational modeling strategy has been developed that accounts for the key physical phenomena responsible for densification of porous preforms using microwave heating. Specifically, a chemical kinetics model has been formulated for Silicon Carbide (SiC) deposition from MTS/H2 precursor. The model is implemented in a pore-resolved reactive transport solver, Quilt. The CVI simulations are performed for several conditions. Initially, parameterized temperature control is imposed to identify optimum conditions for good densification quality at a fraction of processing time. Further, simulations of microwave heating of porous SiC preforms are performed using OpenFOAM. A strategy to achieve and enhance temperature inversion is identified. The resulting temperature profiles are used in the pore-resolved densification simulations to analyze the densification behavior. It is observed that the temperature inversion achieved by microwave heating leads to significant improvements in densification quality and at the same time, keeps the manufacturing time low.

36 MATERIALS SCIENCE↗

Characterization of Organosulfates and Organonitrates in Vertically Resolved Aerosols over the Southern Great Plains Field Campaign Report

Aerosol is the largest individual source of uncertainty in assessing the Earth’s radiative balance due to poorly understood aerosol processes and aerosol-cloud interactions. Thus, to constrain aerosol’s climate impacts, it is critical to understand aerosol compositions at high elevations. Organic material accounts for about 20-80% of total aerosol mass in the atmosphere and is known to play a key role in aerosol processes and cloud formation. In prior research, characterization of aerosol molecular composition has been mostly carried out at ground level. In these ground-level measurements, total speciated organosulfates and organonitrates may contribute significantly to organic aerosol mass on the order of 10-40%, most of which are from abundant biogenic volatile organic compounds such as isoprene and monoterpenes. Organosulfates and organonitrates have also been shown to be present in substantial amounts in cloud water, suggesting that they can be transported and/or directly formed at high elevations. On the other hand, the less common and more expensive aircraft studies have performed high-altitude measurements, but were often limited in the number of different altitudes sampled and in sample collection under consistent conditions. Therefore, the vertical distributions of key aerosol components such as organosulfates and organonitrates are not well understood and have not been extensively studied. Nevertheless, it is important to understand how organosulfates and organonitrates in aerosols are distributed, transformed, participate in aqueous-phase processes, and interact with high-altitude clouds.

54 ENVIRONMENTAL SCIENCES↗

Between Harmonic Crystal and Glass: Solids with Dimpled Potential-Energy Surfaces Having Multiple Local Energy Minima

Solids with dimpled potential-energy surfaces are ubiquitous in nature and, typically, exhibit structural (elastic or phonon) instabilities. Dimpled potentials are not harmonic; thus, the conventional quasiharmonic approximation at finite temperatures fails to describe anharmonic vibrations in such solids. At sufficiently high temperatures, their crystal structure is stabilized by entropy; in this phase, a diffraction pattern of a periodic crystal is combined with vibrational properties of a phonon glass. As temperature is lowered, the solid undergoes a symmetry-breaking transition and transforms into a lower-symmetry phase with lower lattice entropy. Here, we identify specific features in the potential-energy surface that lead to such polymorphic behavior; we establish reliable estimates for the relative energies and temperatures associated with the anharmonic vibrations and the solid–solid symmetry-breaking phase transitions. We show that computational phonon methods can be applied to address anharmonic vibrations in a polymorphic solid at fixed temperature. To illustrate the ubiquity of this class of materials, we present a range of examples (elemental metals, a shape-memory alloy, and a layered charge-density-wave system); we show that our theoretical predictions compare well with known experimental data.

36 MATERIALS SCIENCE↗

Repurposing Offshore Infrastructure for Clean Energy (ROICE) vs. Decommissioning – Commercial Considerations

Abstract The Repurposing Offshore Infrastructure for Clean Energy (ROICE) Program, a collaboration of the energy industry and University of Houston, proposes extending the life of up to 1,500 oil and gas platforms in the Outer Continental Shelf of the Gulf of Mexico, USA. Rather than decommissioning or converting them to underwater reefs at the end of their oil and gas production phase, the platforms could be transformed for 10–20+ years of renewable energy development, such as green hydrogen production or carbon dioxide injection and storage, resulting in significant economic, environmental, and social benefits. ROICE has already published papers on the technical and regulatory considerations for such repurposing projects. This paper focuses on the commercial considerations needed to ensure ROICE projects are economically viable and sustainable for all involved. There are potentially many entities that can come together to progress a repurposing project, including oil and gas asset owners, operators, investors, developers, contractors, manufacturers, and regulators. This paper looks at possible combinations of these entities in a ROICE project team and the various commercial agreements that will be needed to ensure a mutually successful outcome. The paper proposes adapting existing industry agreements and templates to suit, including asset transfer agreements (ATA), asset purchase agreements (APA) and joint purchasing agreements (JOA). An ATA or APA can be drafted to address the obligations and issues involved in the sale or transfer of an existing oil and gas facility to the developer, operator, non-operating interests, and financial investors in a ROICE project, while a JOA can be drafted to govern the rights, obligations, and financial requirements of the parties involved in the development and operation.

Legge, G. [Endeavor Management, USA]↗

Burning conditions and transportation pathways determine biomass-burning aerosol properties in the Ascension Island marine boundary layer

African biomass-burning aerosol (BBA) in the southeast Atlantic Ocean (SEA) marine boundary layer (MBL) is an important contributor to Earth’s radiation budget yet its representation remains poorly constrained in regional and global climate models. Data from the Layered Atlantic Smoke Interactions with Clouds (LASIC) field campaign on Ascension Island (-7.95° N, -14.36° E) detail how fire source regions (burning conditions and fuel type), transport pathways, and longer-term chemical processing affect the chemical, microphysical, and optical properties of the BBA in the remote MBL between June and September of 2017. Ten individual plume events characterize the seasonal evolution of BBA characteristics. Inefficient burning conditions, determined by the mass ratio of refractory black carbon to above-background carbon monoxide (rBC:ΔCO), enhance organic- and sulfate-rich aerosol concentrations in June–July. In contrast, the heart of the burning season exhibited higher rBC:ΔCO values indicative of efficient burning conditions, correlating with more rBC-enriched BBA. Toward the end of the burning season, a mix of burning conditions results in increased variation of the BBA properties. The BBA transit to Ascension Island was predominantly through slow-moving pathways in the MBL and lower free troposphere (FT), facilitating prolonged chemical transformations through heterogeneous and aqueous phase processes. Heterogeneous oxidation can persist for up to 10 days, resulting in a considerable decrease in organic aerosol (OA) mass. OA to rBC mass ratios (OA:rBC) in the MBL between 2 and 5 contrast to higher values of 5 to 15 observed in the nearby FT. Conversely, early-season aqueous-phase processes primarily contributed to aerosol oxidation and some aerosol production, but not appreciable aerosol removal. These two chemical processes yield more light-absorbing BBA in the MBL than in the FT and explain the notably low scattering albedo at 530 nm (SSA 530 ) values (< 0.80) at Ascension Island. This study establishes a robust correlation between SSA 530 and OA:rBC across both MBL and FT, underscoring the dependency of optical properties on chemical composition. These findings highlight how the interplay between chemical composition and atmospheric processing can be improved in global and regional climate models. Questions remain on the mixing of aerosols with different pathway histories, and on what accounts for the doubling of the mass absorption coefficient in the boundary layer.

54 ENVIRONMENTAL SCIENCES↗

Systematic and statistical uncertainties of the hilbert-transform based high-precision FID frequency extraction method

Pulsed nuclear magnetic resonance (NMR) is widely used in high-precision magnetic field measurements. The absolute value of the magnetic field is determined from the precession frequency of nuclear magnetic moments. The Hilbert transform is widely used to extract the phase function from the observed free induction decay (FID) signal and then its frequency. In this paper, a detailed implementation of a Hilbert-transform based FID frequency extraction method is described. How artifacts and noise level in the FID signal affect the extracted phase function are derived analytically. A method of mitigating the artifacts in the extracted phase function of an FID is discussed. Correlations between noises of the phase function samples are studied for different noise spectra. We discovered that the error covariance matrix for the extracted phase function is nearly singular and improper for constructing the $\chi^2$ used in the fitting routine. A down-sampling method for fixing the singular covariance matrix has been developed, so that the minimum $\chi^2$-fit yields properly the statistical uncertainty of the extracted frequency. Other practical methods of obtaining the statistical uncertainty are also discussed.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

The short-term comprehensive impact of the phase-out of global coal combustion on air pollution and climate change

With the continuous intensification of global warming, the reduction and ultimate phase-out of coal combustion is an inevitable trend in the future global energy transformation. This study comprehensively analyzed the impact of phasing out coal combustion on global emissions and concentrations of air pollutants, radiative fluxes, meteorology and climate using Community Earth System Model 2 (CESM2). The results indicate that after the global phase-out of coal combustion, there is a marked decrease in the concentrations of sulfur dioxide (SO 2 ), nitrogen oxides (NO x ) and fine particulate matter (PM 2.5 ), with some regions experiencing a reduction of exceeding 50%. There is no significant change in global ozone (O 3 ) concentration. There are decreasing AOD and positive radiative fluxes globally in the short term, though the cloud contributes minor negative radiative fluxes. The global air temperature may increase by approximately (0.02 ± 0.15) °C on average with regional and seasonal variations, and the precipitation may potentially increase by approximately (2.7 ± 40.6) mm yr −1 globally and over 20% in equatorial regions in the short term. But combined with the decreasing trend of cloud water content in the Northern Hemisphere, it indicates a potential increase in the extremity of precipitation events. In conclusion, this study provides references for global control of air pollution, mitigation strategies of climate change, and transformation of energy structures under the objective of “carbon neutrality”, such as focusing on the negative climate impacts of exacerbating regional warming and increasing extreme precipitation resulting from the rapid reduction of aerosols in the short term.

01 COAL, LIGNITE, AND PEAT↗

Black Start of Unbalanced Microgrids Harmonizing Single- and Three-Phase Grid-Forming Inverters: Preprint

As power systems are transforming with increasing penetrations of inverter-based resources (IBRs), system restoration using IBRs has drawn attention. Using distributed grid-forming (GFM) assets located near critical loads, either three-phase or single-phase, to establish microgrid voltages in the absence of a bulk grid, a distribution system could obtain high system survivability. For swift and secure recovery of a critical load in a single-phase lateral, local single-phase GFM inverters can form a microgrid, and then it can be combined with a neighbouring grid with the inverters remaining in GFM mode for voltage and frequency regulation until the bulk grid comes online. It leads to dynamic interoperation of single-phase GFM inverters with three-phase ones in the black start process. This paper studies the novel approach with electromagnetic transient (EMT) simulations. To evaluate the potential and the technical challenges of the heterogeneous IBR-driven black start, three-phase and single-phase GFM inverter models are developed, including negative-sequence control for voltage balance and a phase-by-phase current limiter (three-phase) and current magnitude limiter (single-phase). To examine dynamic aspects of the black-start process, the EMT simulation also models transformer and motor dynamics emulating their inrush and startup behavior as well as network dynamics. Involvement of grid-following assets to facilitate the black-start process is also modeled. By allowing multiple GFM inverters to collectively black start without leader-follower coordination, regardless of phases, a system can achieve extreme resilience. An inverter-driven black start of a heavily unbalanced 2-MVA distribution feeder using 1 three-phase and 3 single-phase GFM inverters is demonstrated. The simulation shows the heterogeneous system can maintain stability with the single-phase GFM dynamics coupled with the three-phase one.

grid-forming inverter↗

Open Source Fault-tolerant Grid Frequency Measurement for Solar Inverters

The Discrete Fourier transform (DFT) based measurement algorithms are one of the most common measurement algorithms for grid parameter estimation such as rms, phase angle, frequency. Over the past few years, many DFT based algorithms have been developed to enhance its measurement accuracy under steady-state and/or dynamic grid conditions. For example, an adaptive band-pass filter utilizing exponential modulation filter has been proposed to reduce measurement errors at the presence of large frequency deviations. Measurement accuracy of different algorithms including FIR filter, extended Kalman filtering (EKF), and enhanced DFT method have been compared in detail under different grid conditions. Two artificial signals that have 90-degree phase difference were constructed by the Clarke transformation to address the frequency spectrum leakage of DFT. A multi-module approach was developed to enhance both steady-state and dynamic measurement accuracies, in which each module was developed to eliminate some specific errors. Besides DFT-based measurement algorithms, some signal model-based algorithms have been developed to further improve the accuracy under dynamic conditions. However, a key drawback of the state-of-the-art algorithms is that they cannot perform measurements accurately during system transient faults. In the Blue Cut Fire event, there was a phase angle jump of about 26 degrees in the voltage waveform during the transient fault. The phase angle jump fault will cause waveform discontinuity, and these algorithms will fail to provide reliable measurements during this period because they typically assume the waveform to be measured is continuous, no matter what method (DFT, PLL, EKF, FIR, or Taylor WLS) is used for estimation. In fact, the measurement errors during the system transient faults like phase-jump is not required in the IEEE Standard. As a result, although a measurement instrument can pass the strict IEEE Standard, it could still be the source of the problem in the future if we have similar system transient faults, which could happen again. Therefore, developing the fault-tolerant measurement technology is the key to solve the problem.

14 SOLAR ENERGY↗

Proton beam quality enhancement by spectral phase control of a PW-class laser system

We report on experimental investigations of proton acceleration from solid foils irradiated with PW-class laser-pulses, where highest proton cut-off energies were achieved for temporal pulse parameters that varied significantly from those of an ideally Fourier transform limited (FTL) pulse. Controlled spectral phase modulation of the driver laser by means of an acousto-optic programmable dispersive filter enabled us to manipulate the temporal shape of the last picoseconds around the main pulse and to study the effect on proton acceleration from thin foil targets. The results show that applying positive third order dispersion values to short pulses is favourable for proton acceleration and can lead to maximum energies of 70 MeV in target normal direction at 18 J laser energy for thin plastic foils, significantly enhancing the maximum energy compared to ideally compressed FTL pulses. The paper further proves the robustness and applicability of this enhancement effect for the use of different target materials and thicknesses as well as laser energy and temporal intensity contrast settings. We demonstrate that application relevant proton beam quality was reliably achieved over many months of operation with appropriate control of spectral phase and temporal contrast conditions using a state-of-the-art high-repetition rate PW laser system.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Thermodynamic modeling of the Ti-Nb-Sn system and its preliminary applications in aging precipitation and solidification simulations

Abstract Gasping the thermodynamic and phase diagram of the materials being studied is essential for further guiding the development process, like composition design, heat treatment optimization and so on. Nb and Sn are common alloying elements in Ti-based biomedical alloys. In this work, the Ti-Nb-Sn system has been thermodynamically assessed based on experimental phase equilibria. A self-consistent thermodynamic description for the Ti-Nb-Sn system including two ternary compounds Ti 3 Nb 3 Sn 2 and Ti 3 NbSn has been obtained. Based on the comparison between the calculated and experimental phase equilibria, the reliability of the presently obtained thermodynamic description has been verified. The influence of Nb and Sn contents on the precipitation temperature and amount of hcp α -Ti phase in bcc β -(Ti, Nb) has been investigated based on thermodynamic calculations. Results show that the increase of Nb or Sn content in Ti-based alloys would significantly reduce the precipitation amount and transformation temperature of hcp α -Ti phase in bcc β -(Ti, Nb) and the influence of Nb content is more effective than that of Sn. The solidified phases in several as-cast alloys have been thermodynamically predicted based on Scheil solidification simulations. The presently developed thermodynamic description of the Ti-Nb-Sn system can serve as a sub-ternary one in the thermodynamic database for multicomponent Ti-based biomedical alloys.

Peng, Yingbiao↗

Order-to-disorder transition due to entropy in layered and 2D carbides

In compositionally complex materials, there is controversy on the effect of enthalpy versus entropy on the structure and short-range ordering in so-called high-entropy materials. Here, to help address this controversy, we synthesized and probed 40 M 4 AlC 3 layered carbide phases containing two to nine metals and found that short-range ordering from enthalpy was present until the entropy increased enough to achieve complete disordering of the transition metals in their atomic planes. We transformed all of these layered carbide phases into two-dimensional (2D) sheets and showed the effects of the order versus disorder on their surface properties and electronic behavior. This study suggests the key effect that the competition between enthalpy and entropy has on short-range order in multicompositional materials.

Wyatt, Brian C. [Purdue Univ., West Lafayette, IN ↗