Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Chemical initiation”

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 415 records · Page 23

Effects of non-thermal termolecular reactions on wedge-induced oblique detonation waves

The shock-induced combustion ramjet (Shcramjet) based on oblique detonation waves (ODWs) is among the promising choices for hypersonic propulsion systems. An understanding of the ignition, propagation, and stability of ODWs is critical to harnessing their propulsive potential. In such high speed reacting flows, there is a high probability of occurrence of non-thermal reactions due to the presence of non-trivial amounts of highly reactive radicals including H, O and OH apart from O 2 as demonstrated recently [M. P. Burke, S. J. Klippenstein, Nat. Chem. 9 (2017) 1078–1082, Y. Tao, A. W. Jasper, Y. Georgievskii, S. J. Klippenstein, R. Sivaramakrishnan, Proc. Combust. Inst. 38 (2021) 515–522]. The present work focuses on examining the initiation, propagation and structure of oblique detonation waves in stoichiometric H 2 -air mixtures through numerical simulations with and without non-thermal reactivity on a two-dimensional adaptive grid. Non-thermal reactions were included in the macroscopic kinetic model as chemically termolecular reactions facilitated by the H + OH radical-radical recombination and the H + O 2 radical-molecule association reactions. Since, the non-thermal reactions result in a corresponding decrease in the reaction fluxes of the incipient recombination/association reactions, an additional simulation was performed by applying corrections to the respective incipient recombination/association rate constants using the methodology demonstrated by Tao et al. [Y. Tao, A. W. Jasper, Y. Georgievskii, S. J. Klippenstein, R. Sivaramakrishnan, Proc. Combust. Inst. 38 (2021) 515–522]. Results show that, under ODWE relevant conditions, non-thermal reactivity fundamentally alters the induction length, intensity as well as the structure of the ODW. Specifically, it is found that non-thermal reactivity leads to a noticeable reduction in initiation length and a simultaneous increase in instantaneous peak heat release rate and the degree of unsteadiness of the ODW. Finally, statistical analysis of key thermo-chemical variables is performed to elucidate the important species as well as reactions responsible for the observed variations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Chemical preintercalation synthesis approach for the formation of new layered tungsten oxides

Tungsten oxide, WO 3 ·nH 2 O, is a unique layered oxide material that offers enhanced performance in electrochromic and energy storage applications. Herein, we report the formation of a new, never previously synthesized, Na-containing layered tungsten oxide phase, Na 0.20 WO 3 ·0.81H 2 O, using a chemical preintercalation approach. The structure and composition of this novel phase were investigated via microscopy, spectroscopy, and diffraction methods. Electrochemical cycling of Na 0.20 WO 3 ·0.81H 2 O electrodes revealed initial discharge capacities of 37.43 mAh g -1 , 480.8 mAh g -1 , and 253.2 mAh g -1 in aqueous H 2 SO 4 cells (potential window of - 0.2–0.8 V vs. Ag/AgCl), non-aqueous Li-ion cells (potential window of 0.1–4.0 V vs. Li/Li + ), and non-aqueous Na-ion cells (potential window of 0.1–4.0 V vs. Na/Na + ), respectively. Additionally, a reversible, pressure-induced color change from pale yellow to dark brown/black was observed for the Na 0.20 WO 3 ·0.81H 2 O sample when it was placed under pressures of 1000 psi or higher. Our results demonstrate the viability of chemical preintercalation synthesis approach to produce new oxide phases with interesting functional properties.

36 MATERIALS SCIENCE↗

Stress-Alteration Enhancement of the Reactivity of Aluminum Nanoparticles in the Catalytic Decomposition of exo -Tetrahydrodicyclopentadiene (JP-10)

High energy density aluminum nanoparticles (AlNP) upon thermal annealing followed by super-quenching results in elevated stress levels in the metallic core and reduced surface energy at the core-shell interface. Isomer-selective vacuum ultraviolet based photoionization mass spectrometry coupled to a high temperature chemical microreactor reveals that these stress-altered (SA) AlNPs exhibit distinctive temperature-dependent reactivities towards catalytic decomposition of the hydrocarbon jet fuel exo-tetrahydrodicyclopentadiene (JP-10, C 10 H 16 ) compared to untreated (UN) AlNPs. SA-AlNPs show a delayed initiation of the decomposition for JP-10 by 200 K relative to the UN-AlNPs, however, the full decomposition is achieved at 100 K lower temperature. Furthermore, there is less oxygenated products that are generated from the alumina surface-induced heterogeneous oxidation process, and a larger fraction of closed and open-shell hydrocarbons. Furthermore, chemical insight bridging the reactivity order of SA-AlNPs at low and high temperatures, simultaneously, is obtained via a detailed examination of the product branching ratios obtained in this study.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Fluorine-Free Superhydrophobic Coatings: Rapid Fabrication and Highly Efficient Oil/Water Separation

In this work, a rapid method is demonstrated to obtain a fluorine-free superhydrophobic/superoleophilic coating by a simple two-step method: dip-coating and oven curing. The chemical structure of the coating is based on the crosslinking reaction of an alkyl methacrylate, a dimethacrylate (crosslinker) and a silane (adhesive), using AIBN as initiator, toluene as solvent, and silica nanoparticles to enhance surface roughness. Chemorheology results show that the coating is fully cured even in 20 min at 100 °C, exhibiting a water contact angle (CA) of 162 ± 2°, sliding angle (SA) of 4 ± 1°, nanometrical structures throughout the surface, and excellent adhesion properties for the mesh screen. The coating exhibits outstanding oil/water separation efficiency (96–99%) for seven different types of oil (gasoline, diesel, petroleum ether, hexane, toluene, chloroform, and dichloromethane), and in addition presents high recyclability. Based on nonisothermal TGA, the activation energy of degradation, calculated using the Kissinger and Ozawa models, is 113.5 and 114 kJ mol -1 , respectively. Finally, the coating is thermally aged at 150 °C: the CA exhibits a smooth decrease up to 129 ± 1° at 120 min, and FTIR analysis shows structural changes possibly related to the generation of thermal oxidation products.

36 MATERIALS SCIENCE↗

The microscopic mechanisms of high temperature oxidation of Haynes 282

Nickel-based superalloy finds widespread applications in aerospace and extreme environments. They are known for their high temperature oxidation resistance under extended periods. However, the oxide formation and evolution which sets the stages for the later parabolic oxidation kinetics is not fully understood. This paper aims to provide new insights into the transient stage oxidation mechanism and kinetics of a typical Ni-based superalloy (Haynes 282). While tracking the chemical and microstructural evolution under micrometer scale at 800°C, we show that the oxide scale and its grain boundary species change significantly during the initial stages of oxidation and can have a profound impact on the oxidation kinetics. Cr 2 O 3 forms initially at the grain boundary along with minor amount of Al 2 O 3 . Then, the grain boundary region is enriched with copious amounts TiO 2 while Cr migrates away from the grain boundary. A noticeable change in the isothermal oxidation kinetics observed likely results from a mechanistic change from uniform surface oxidation to preferential outward diffusion of Ti 4+ ions through the grain boundary. Through first-principles calculations combined with energy dispersive spectroscopy, we confirmed the preferential outward diffusion of Ti 4+ ions as the diffusion barrier is lower for Ti than Cr along the grain boundaries. In conclusion, these findings highlight the critical role of early-stage grain boundary oxidation dynamics in dictating the long-term oxidation resistance of Ni-based superalloys and provide a foundation for future strategies to enhance their performance in extreme environments.

DFT↗

Tracking the Oxidation of Silicon Anodes Using Cryo-EELS upon Battery Cycling

Silicon is a high-capacity material for the anode of a rechargeable lithium-ion battery. One of the fundamental challenges for using Si in anodes is capacity fading, which has been revealed to be partially associated with the interfacial instability between the Si and liquid electrolyte due to the large volume swing of Si upon charging and discharging. Smart nanoscale design concepts, either pre-synthesized or formed in-situ, have led to the mitigation of the detrimental factors associated with the volume swing of Si. However, it has never been clear as how the chemical state of Si evolves and contributes to the capacity fading upon battery cycling. Here, we use cryo-electron energy loss spectroscopy (EELS) to directly monitor, at sub-nanometer scale, the chemical evolution of Si upon battery cycling. We discover that during the cycling process, Si particles are progressively oxidized to form SiO 2 , which is initiated from the particle surface and gradually penetrates toward the interior of the particle, directly contributing to the capacity fading. Possible mechanisms of Si oxidation are postulated. We further show how the cycling stability can be improved by an electrolyte additive to form an effective passivation layer, representatively, even a small concentration of fluoroethylene carbonate (FEC) causes the formation of an LiF layer on the Si nanoparticle surface that prevents Si oxidation and improves cycling stability. Furthermore, the present work unveils Si oxidation as a previously unrecognized factor that contributes to capacity fading, therefore providing insight into the design of anodes with Si based materials.

25 ENERGY STORAGE↗

Experimental demonstration of necessary conditions for X-ray induced synthesis of cesium superoxide

Absorption of sufficiently energetic X-ray photons by a molecular system results in a cascade of ultrafast electronic relaxation processes which leads to a distortion and dissociation of its molecular structure. Here, in this work, we demonstrate that only decomposition of powdered cesium oxalate monohydrate induced by monochromatic X-ray irradiation under high pressure leads to the formation of cesium superoxide. Whereas, for an unhydrated form of cesium oxalate subjected to the same extreme conditions, only degradation of the electron density distribution is observed. Moreover, the corresponding model of X-ray induced electronic relaxation cascades with an emphasis on water molecules’ critical role is proposed. Our experimental results suggest that the presence of water molecules in initially solid-state systems (i.e. additional electronic relaxation channels) together with applied high pressure (reduced interatomic/intermolecular distance) could potentially be a universal criteria for chemical and structural synthesis of novel compounds via X-ray induced photochemistry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Solubility of An(OH) 3 am and An(OH) 3 cr in the presence of Citrate and EDTA

The solubility of An(III), which pertains to Pu(III), Am(III), and Cm(III), in brine is an important contributor in WIPP (Waste Isolation Pilot Plant) performance assessment (PA) models to calculate the potential release of transuranic elements in the near-field environment of the WIPP repository. Americium, although present in much smaller quantities than plutonium (~ 336 kg on year 2033 is potentially WIPP-bound primarily as Am-241), has a significant impact on release due to its higher activity during the earlier times of repository history (~26% initially, decreasing to 17% and ~0% at 1000 and 10,000 years after emplacement) (SOTERM, 2019). Lanthanide analogues such as Nd 3+ possess physical and chemical characteristics that allow them to be used to examine the chemical behavior of trivalent actinides (Lucchini, et al., 2007).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

COMPUTATIONAL MODELING OF IGNITION AND PREMIXED FLAME PROPAGATION INITIATED BY A PRE-CHAMBER TURBULENT JET

Addressing the pressing need for reduced carbon emissions, Turbulent Jet Ignition (TJI) emerges as a promising technology for ultra-lean combustion, offering enhanced thermal efficiencies and minimized cyclic variability in spark-ignited engines. To facilitate rapid testing and integration of this technology, a robust computational modeling framework is crucial. This study delves into the predictive capabilities of computational models for main-chamber ignition and premixed flame propagation using a single-cycle TJI rig measured by Biswas et al. (Applied Thermal Engineering, vol 106, 2016). Employing an open-source compressible flow simulation solver with Large Eddy Simulation (LES) for turbulence modeling, the investigation integrates the conventional Laminar Finite Rate Chemistry (LFRC) model alongside the transported Probability Density Method (PDF) for turbulence-chemistry interaction. A fully-consistent Eulerian Monte-Carlo Fields (EMCF) method is utilized to approximate the transported PDF, while Interaction by Exchange with Mean is employed to close micro-mixing terms in stochastic differential equations. A reduced chemical reaction mechanism with 21 species and 84 reactions (DRM-19) is used for solving chemical kinetics, and a double Gaussian energy deposition model is used to approximate the spark ignition in the pre-chamber. An unstructured O-grid mesh with 0.3 million cells in the prechamber and 1 million cells in the main chamber is employed. Results are divided into two phases: pre-chamber initialization and full TJI simulations. Validation of the predicted pre-chamber flame propagation and the lean ignition in the main-chamber is carried out by using available experimental data. Under quiescent conditions, both the LFRC and transported PDF methods largely underestimate the flame speed and subsequent pressure growth in the pre-chamber. A linear momentum forcing technique is applied to investigate the impact of initial turbulence in the pre-chamber, demonstrating a notable influence on flame propagation. Fine-tuning of the forcing coefficient reproduces the sudden pressure growth observed in the experiment. The experimentally validated pre-chamber simulation serves as the initial condition for the full TJI simulations. It is found that the LFRC model fails to predict lean-ignition in the main-chamber, resulting in a misfiring event. Incorporation of turbulence-chemistry interaction using the transported PDF method substantially improves the prediction of the ignition event in the main-chamber, achieving fair qualitative agreement and quantitative validation of combustion parameters within ±10% of the reported experimental data. The rich simulation results consisting of a full set of statistical description of the thermo-chemical states enable us to gain deep insights into the ignition mechanisms in the main chamber, which is limited when done experimentally. A novel dual ignition phenomenon is revealed in the TJI rig for the first time. Initially, a primary ignition kernel is formed at a downstream location which eventually detaches from the main jet. As the jet momentum decreases, a secondary ignition event follows, this time at a more upstream location which eventually combines with the primary ignition kernel to form a single connected flame front. Investigation of these ignition sequences in chemical composition space reveal distinct differences between the two. The primary ignition event in the main-chamber is followed by a large concentration of active radicals from the pre-chamber jet, accelerating the chain-branching steps, characterizing what has been referred to as flame ignition. In contrast, the secondary ignition occurs in the absence of active radicals in the pre-chamber jet, hence characterized as jet ignition. Further analysis of the effect of pre-chamber jet characteristics on lean ignition in the main-chamber is conducted by setting up cases with different initial pressure ratios (por) between the two chambers, a non-dimensional parameter, ranging from 1.2 to 3.2. As the initial pressure ratio increases, jet momentum increases, with dual ignition observed in cases above por= 2.2. Case with por= 3.2 lead to misfiring. The effect of ignition sequence on global combustion characteristics of TJI is analyzed. Dual ignition events lead to non-monotonicity in combustion characteristics such as global reaction progress variable, flame penetration, and global heat release rate. In dual ignition events, although the rate of fuel consumption and global heat release rate is initially lower, the secondary ignition leads to a sudden increase in flame surface area, resulting in a sudden jump and promoting the overall performance of the TJI system.

42 ENGINEERING↗

Reduction of Ultrafine Hematite for a Spouting Bed Chemical Looping Process - Presentation

Chemical looping combustion research has been a topic of great interest since the mid-1990s, although the concept can be traced much further back. It is a combustion technology for fuel utilization that uses a two-reactor system to produce a carbon free effluent stream from the air reactor and a carbon dioxide rich stream from the fuel reactor with oxygen shuttled between the reactors with an oxygen carrier. Chemical looping combustion has been studied extensively because of this inherent separation of the carbon dioxide combustion product from the atmospheric nitrogen over that past 25 years with hundreds of research papers and tens of books. The oxygen carrier is, in general, a solid such as naturally occurring hematite or an engineered carrier like a copper and iron oxide supported on alumina. System technical and economic research activities conducted at the National Energy Technology Laboratory have shown that the costs associated with the carrier, initial and replacement, are a significant factor in the economic viability of the technology. To overcome this economic challenge, the work reported on in this paper investigates the use of a waste iron oxide as the chemical looping oxygen carrier. This waste material is quite unique as it has a particle size of 0.7 um. Given the particle size, the reactions are quite fast. However, the development of a contactor reactor to mix the fuel with the carrier is a fluidization challenge as the material is a Geldart group C material and highly cohesive. Normally, Geldart group C materials are fluidized with either a vibrating bed. This technique is acceptable for low temperature pharmaceuticals production but does not lend itself to much larger reactors used for power production or operation at temperatures in the order of 1100 to 1200 C. Therefore, a hydrodynamic solution was needed to break up the cohesive structures to replace the vibration. In this instance, it was hypothesized that a spouting bed of large Geldart group B or Geldart group D particles could provide the necessary energy to break up the cohesive structures of the Geldart group C material. This paper presents the reduction performance for this novel carrier material in a laboratory spouting bed reactor.

Breault, Ronald↗

Epsilon Metal: A Waste Form for Noble Metals from Used Nuclear Fuel

Epsilon metal (e-metal) is the metallic phase that forms as inclusions at the grain boundaries in the UO2 fuel during reactor operation. This metal is composed of Pd, Mo, Rh, Ru, and Tc. These metallic inclusions are insoluble in strong acid and remnants of these metallic inclusions have been found in the UO2 matrix that remains from the natural reactors in Gabon that were active 1.8 billion years ago, therefore e-metal should be an excellent waste form for the immobilization of the long-lived isotopes 107Pd (6.5 ? 106 a) and 99Tc (2.13 ? 105 a), with 99Tc being the isotope of interest for repository performance. Therefore, the chemical durability of this potential waste form is assessed in this study. Typically, corrosion rates for metallic materials are measured electrochemically because they are quick, inexpensive, and can reveal the mechanism by which a metal corrodes, at least initially. However, in a repository the waste form would be subjected to slowly flowing water without an applied electrical potential over long time periods. Therefore, the corrosion rates of e-metal specimens were measured with both electrochemical tests and the single-pass flow-through test (SPFT). Potentiodynamic and potentiostatic polarization results suggest that a thin passive film exists on the alloy surface, which seems to be responsible for its high corrosion resistance. Additionally, X-Ray photoelectron spectroscopic results suggest that Pd oxides are significantly enriched in the passive film Results from the SPFT show that the dissolution rates were weakly dependent on pH. Only Mo and Re were found in solution and were used for the calculation of the dissolution rates. In general, the electrochemically determined corrosion rates agree reasonably well with the initial dissolution rate measured with the SPFT test, but they are about one or two orders of magnitude higher than the steady state rates. The causes for this discrepancy are discussed.

epsilon-metal, nuclear waste form, corrosion, sing↗

Solid-to-plasma transition of polystyrene induced by a nanosecond laser pulse within the context of inertial confinement fusion

Laser Direct Drive (LDD) Inertial Confinement Fusion (ICF) involves irradiating a spherical target of thermonuclear fuel coated with an ablator, usually made of polystyrene. Laser energy absorption near the target surface leads to matter ablation, hydrodynamic shocks, and ultimately capsule implosion. The conservation of spherical symmetry is crucial for implosion efficiency, yet spatial modulations in laser intensity can induce nonuniformities, causing the laser imprint phenomenon. Understanding laser imprint, especially considering the initial solid state, is essential for advancing LDD ICF. A first microscopic model of solid-to-plasma transition was built in 2019, accounting for laser absorption in the solid state with a band structure based ionization model. This model has been improved to include chemical fragmentation and a more accurate description of electron collision frequency in various matter states. The latest development involves assessing the model reliability by comparing theoretical predictions with experimental observations. Despite the success of this approach, questions remain, leading to further investigations and observations under different irradiation conditions. This work presents an experiment with a nanosecond pulse, taking into account hydrodynamic effects, and measures transmission dynamics over the entire laser beam area to observe two-dimensional effects. As a result, the objective is to adapt the theoretical model, couple it with a hydrodynamic code, and observe additional effects related to the initial solid state.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Far-from-equilibrium search for the QCD critical point

Initial conditions for relativistic heavy-ion collisions may be far from equilibrium (i.e., there are large initial contributions from the shear stress tensor and bulk pressure) but it is expected that on very short time scales the dynamics converge to a universal attractor that defines hydrodynamic behavior. Thus far, studies of this nature have only considered an idealized situation at LHC energies (high temperatures T and vanishing baryon chemical potential µ B = 0) but, in this work, we investigate for the first time how far-from-equilibrium effects may influence experimentally driven searches for the Quantum Chromodynamic critical point at RHIC. We find that the path to the critical point is heavily influenced by far from equilibrium initial conditions where viscous effects lead to dramatically different {T, µ B } trajectories through the QCD phase diagram. We compare hydrodynamic equations of motion with shear and bulk coupled together at finite µB for both DNMR and phenomenological Israel-Stewart equations of motion and discuss their influence on potential attractors at finite µ B and their corresponding {T, µ B } trajectories.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Towards a heat- and mass-balanced kinetic model of TATB decomposition

We report TATB (1,3,5-triamino-2,4,6-trinitrobenzene) was thermally degraded by two small-scale analytical methods – simultaneous differential scanning calorimetry and thermogravimetric analysis (SDT) and a hot-stage microscope with Fourier Transform Infrared (FTIR) analysis capabilities. SDT used ramped heating, isothermal soaking, and thermal pretreatment at various conditions. The heat flow and mass loss were monitored during various treatment conditions to derive chemical decomposition kinetics and Arrhenius parameters. FTIR experiments used isothermal heating, and changes were monitored spectroscopically. Solid samples generated at specific conditions were collected from both methods and were analyzed by DMSO extraction followed by chemical speciation by optical and mass spectrometric methods. Characterization provided the following reaction insights: TATB decreases in a sigmoidal pattern in isothermally heated samples. Other soluble products gradually increase in concentration and then abruptly decline in concentration during the second exotherm, such as diamino-dinitro-benzofurazan and amino-nitro-benzodifurazan. FTIR showed gradual changes in the amino and nitro functionality, shifting positions and decreasing intensity for the first 40 min. Then the solid gradually appeared more like an amorphous C with N incorporated, similar to previous studies on thermally degraded TATB-type materials. Extracted residues (DMSO-soluble components removed) examined by FTIR showed an abrupt change in chemical composition between 40- and 45-min isothermal treatment, indicating early forming solids are different than later forming residues. A reliable mass- and energy-balanced global reaction network must include at least two autocatalytic reactions, either in parallel or series, and at least one must have an explicit initiation reaction having a low activation energy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The Reactivity of Polyethylene Microplastics in Water under Low Oxygen Conditions Using Radiation Chemistry

Polyethylene (PE) is an intensely utilized polymer, which has consequently led to it becoming a common environmental contaminant. PE and other plastic waste are known to be highly persistent in surface waters; however, chemical and physical changes do take place over time, dependent mostly on highly variable natural conditions, such as oxygen (O2) availability. Gamma radiation was used to generate reactive oxygen species, namely hydroxyl radicals, in initially aerated aqueous solutions to simulate the natural weathering of microplastics in waters where there are fluctuations and often depletions in dissolved O2. The headspace of the irradiated PE-containing solutions was probed for the formation of degradation products using solid-phase microextraction (SPME) fibers in combination with gas chromatography mass spectrometry (GCMS). The major species detected were n-dodecane, with trace levels of tridecane, 2-dodecanone, and hexadecane, which were believed to be predominately adsorbed in the PE microplastics in excess of their aqueous solubility limits. Surface characterization by Raman spectroscopy and light and dark field microscopy indicated no change in the chemical composition of the irradiated PE microplastics under low O2 to anaerobic conditions. However, morphological changes were observed, indicating radical combination reactions.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Molecular and Microstructural Bases for Understanding Microplastic Origin, Transport, and Fate

The formation of micro and nanoplastics from parent plastics is one of the most challenging ecological threats since plastic pollution reaches even the most desolate places on earth and is expected to double by 2030. The goal of this project is to understand physico-chemical properties on microplastics formation, transport, and fate by analyzing the microplastics found in environmentally collected samples at SRS, supported by lab-based experiments. An initial assessment found and identified microplastics and micronizing plastics in Fourmile Branch and the Savannah River. Fourmile Branch also carried trace-level radioactive contaminants associated with the Site’s legacy missions, which may incorporate into entrained microplastics.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Regional Feedstock Partnership Biomass Quality Assessment Final Report

The United States (U.S.) Department of Energy (DOE) developed the Billion-Ton Vision to enable production of one-billion tons of sustainable, reliable biomass for the bioenergy industry by 2030 (Perlack et al., 2005). The Sun Grant Regional Feedstock Partnership (RFP) was organized to fill information gaps and validate biomass yield assumptions related to the Billion-Ton Study (Owens, 2018; Owens, Karlen, and Lacey, 2016). Along with the more than 130 scientific publications generated from these studies, yield and sustainability data from the RFP field trials not only validated the Billion-Ton estimates, but were critical in developing both the U.S. Billion-Ton Update report in 2011 and the 2016 Billion-Ton Report (DOE, 2011; 2016). The intention of this biomass quality assessment report is to build on these initial successes from the RFP field trials by focusing on variability in biomass quality data necessary to evaluate conversion performance. This report contains a summary of chemical quality results from samples collected as part of the RFP field trials. This report focuses on assessment of the impact of experimental agronomic designs on biomass properties followed by analyses of the impact of environmental and production variables on biomass properties. Datasets include species and other genetic variables, fertilizer treatments, harvest information, and yield, as well as other publicly available data such as precipitation, temperature, soil properties, and drought. The key outcomes from this chemical quality focused assessment have included: • Complete evaluation of the impacts of agronomic designs, genetics, and environmental conditions on chemical properties for Miscanthus, switchgrass, sorghum, energycane, mixed perennial grasses, and shrub willow short-rotation feedstocks • Over 30 peer review publications and technical reports focused on variability in quality data • Development of spatial and temporal environmental quality prediction maps for Miscanthus and switchgrass feedstocks allowing for comprehensive evaluation of variability in feedstock chemical quality across U.S. regions and over multiple harvest years

09 BIOMASS FUELS↗

HfZr_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys hafnium-zirconium (Hf-Zr). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Hf and Zr. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,100 randomized atomic structures over 31 chemical compositions.

36 MATERIALS SCIENCE↗