Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “FACTS”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Electron kinetics in a high- Z plasmoid

The problem of the electron dynamics on a closed magnetic field line passing through a high-Z plasmoid is considered. The electron kinetic equation is integrated over bounce motion and pitch angle, reducing the independent variables to a single adiabatic invariant plus time. Integration of the full Landau self-collision operator is carried out exactly, resulting in a nonlinear integro-differential operator in the new invariant. Conservation laws and the H theorem of the integrated self-collision operator are proven. Numerical solutions of the integrated kinetic equation are obtained with a self-consistent quasineutral electric potential, given the initial condition of a cold plasmoid immersed in a hot ambient plasma. The fact that cold electrons are deeply trapped in a potential with a parabolic peak leads to exactly 3/4 the usual rate of collisional heating by the ambient plasma, independent of any other parameters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Kinetic stability of Chapman–Enskog plasmas

In this paper, we investigate the kinetic stability of classical, collisional plasma – that is, plasma in which the mean-free-path $\lambda$ of constituent particles is short compared with the length scale $L$ over which fields and bulk motions in the plasma vary macroscopically, and the collision time is short compared with the evolution time. Fluid equations are typically used to describe such plasmas, since their distribution functions are close to being Maxwellian. The small deviations from the Maxwellian distribution are calculated via the Chapman–Enskog (CE) expansion in $\lambda /L \ll 1$ , and determine macroscopic momentum and heat fluxes in the plasma. Such a calculation is only valid if the underlying CE distribution function is stable at collisionless length scales and/or time scales. We find that at sufficiently high plasma $\beta$ , the CE distribution function can be subject to numerous microinstabilities across a wide range of scales. For a particular form of the CE distribution function arising in strongly magnetised plasma ( viz. plasma in which the Larmor periods of particles are much smaller than collision times), we provide a detailed analytic characterisation of all significant microinstabilities, including peak growth rates and their associated wavenumbers. Of specific note is the discovery of several new microinstabilities, including one at sub-electron-Larmor scales (the ‘whisper instability’) whose growth rate in certain parameter regimes is large compared with other instabilities. Our approach enables us to construct the kinetic stability maps of classical, two-species collisional plasma in terms of $\lambda$ , the electron inertial scale $d_e$ and the plasma $\beta$ . This work is of general consequence in emphasising the fact that high- $\beta$ collisional plasmas can be kinetically unstable; for strongly magnetised CE plasmas, the condition for instability is $\beta \gtrsim L/\lambda$ . In this situation, the determination of transport coefficients via the standard CE approach is not valid.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

The effect of shear flow on the resistive tearing instability

We develop a new scaling theory for the resistive tearing mode instability of a current sheet with a strong shear flow across the layer. The growth rate decreases with increasing flow shear and is completely stabilized as the shear flow becomes Alfvénic: both in the constant-Ψ regime, as in previous results, but we also show that the growth rate is in fact suppressed more strongly in the nonconstant-Ψ regime. As a consequence, for sufficiently large flow shear, the maximum of the growth rate is always affected by the shear suppression, and the wavenumber at which this maximum growth rate is attained is an increasing function of the strength of the flow shear. These results may be important for the onset of reconnection in imbalanced MHD turbulence.

astrophysical plasmas↗
DOE OSTI2023

Habitat severity characteristics structure soil communities at regional and local spatial scales along the Antarctica Peninsula

Antarctic soils provide an excellent setting to test biogeographical patterns across spatial and environmental scales given their relatively simple communities and the dominance of physical factors that create strong environmental gradients. Additional urgency is given by the fact that their unique terrestrial communities are the subject of conservation efforts in a rapidly changing environment. We investigated relationships of soil community assembly and alpha and beta diversity with climatic and environmental parameters across regional and local scales in Maritime Antarctica. We sampled from a regional gradient of sites that differ in habitat severity, ranging from relatively favourable to harsher physicochemical conditions. At the regional scale, bacterial community characteristics and microarthropod abundance varied along this severity gradient, but most measures of fungal communities did not. Microarthropod and microbial communities differed in which soil and climate parameters were most influential, and the specific parameters that influenced each taxon differed across broad and fine spatial scales. This suggests that conservation efforts will need to focus on a large variety of habitat characteristics to successfully encompass diversity across taxa. Because beta diversity was the result of species turnover, conservation efforts also cannot focus on only the most biodiverse sites to effectively preserve all aspects of biodiversity.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Emerging Potassium Metal Anodes: Perspectives on Control of the Electrochemical Interfaces

Potassium metal acts as the anode in emerging potassium metal batteries (KMBs). It further serves as the counter-electrode for potassium ion battery (KIB) half-cells, with its reliable performance being critical for assessing the working electrode material. This first-of-its-kind critical review focuses on the dual challenge of controlling the potassium metal-substrate and the potassium metal-electrolyte interface so as to prevent dendrites. The discussion begins with a comparison of the physical and chemical properties of K metal anodes versus the much oft studied Li and Na metal anodes. Based on established descriptions for root causes of dendrites, the problem should be less severe for K than for Li or Na, while in fact the opposite is observed. The key reason that the K metal surface rapidly becomes dendritic in common electrolytes is its unstable solid electrolyte interphase (SEI). An unstable SEI layer is defined as being non self-passivating. No SEI is perfectly stable during cycling, and all SEI structures are heterogenous both vertically and horizontally relative to the electrolyte interface. The difference between a "stable" and an "unstable" SEI may be viewed as the relative degree to which during cycling it thickens and becomes further heterogeneous. The unstable SEI on K metal leads to a number of interrelated problems, such as low cycling Coulombic efficiency (CE), a severe impedance rise, large overpotentials, and possibly electrical shorting; all of which have been reported to occur as early as in the first ten plating/stripping cycles. Many of the traditional "interface fixes" employed for Li and Na metal anodes, such as various artificial SEIs, surface membranes, barrier layers, the secondary separators, etc. have not been attempted or optimized for the case of K. This is an important area for further exploration, with an understanding that success may come harder than with Li due to K-based SEI reactivity with both ether and ester solvents. The second critical problem with K metal anodes is that they do not thermally or electrochemically wet a standard (untreated) Cu foil current collector. Published experimental and modeling research directly highlights the weak bonding between the K atoms and a Cu surface. Existing surface treatment approaches that achieve improved K wetting are discussed, along with the general design rules for future studies. Also discussed are geometry-based methods to tune nucleation, as well dual approaches where nucleation and SEI structure are tuned through complementary schemes to achieve extended half-cell and full battery stability. We hypothesize that K metal never achieves a planar wetting morphology even at cycle one making the dendrites "baked-in". We propose that classical thin films growth models, Frank van der Merwe (F-M), Volmer Weber (V-W), and Stranski-Krastanov (S-K) can be employed to describe early stage plating behavior. It is demonstrated that island-like V-W growth is the applicable description for the natural plating behavior of K on pristine Cu. Moving forward, there are three inter-related thrusts to be pursued: First, K salt - based electrolyte formulations have to mature and become further tailored to handle the increased reactivity of a metal rather than an ion anode. Second, the K-based SEI structure needs to be further understood and ultimately tuned to be less reactive. Third, the energetics of the K metal - current collector interface must be controlled to promote planar wetting/dewetting throughout cycling.

25 ENERGY STORAGE↗
DOE OSTI2020

Oxygen-Based Anion Redox for Lithium Batteries

The importance of current Li-ion batteries (LIBs) in modern society cannot be overstated. While the energy demands of devices increase, the corresponding enhancements in energy density of battery technologies are highly sought after. Currently, many different battery concepts, such as Li-S and metal-air among many others, have been investigated. However, their practical implementation has mostly been restricted to the prototyping stage. In fact, most of these technologies require rework of existing Li-ion battery manufacturing facilities and will naturally incur resistance to change from industry. For this reason, one specifically attractive technology, anionic redox in transition metal oxides, has gained much attention in the recent years. Its ability to be directly used in already established processes and higher energy density with similar electrolyte formulation make it a key materials research direction for next generation Li-ion batteries. In regular LIBs, the redox active centers are the transition metal cation. In anion redox, both the anion (typically O) and the transition metal cation are utilized as redox centers with enormous implications for increasing energy density. This new material can be highly competitive for replacing the current LIB technologies. However, much is still unknown about its cyding mechanism. Upon activating the O redox couples, most cationic and anionic redox active materials will either evolve O 2 or undergo irreversible structural degradation with associated severe decreases in electrochemical performance. By understanding the transition from full anion redox to partial cationic and anionic redo; we hope readers can gain a deeper understanding of the topic. This Account will focus mainly on the work that was conducted by our group at Argonne National Laboratory. Herein, the phenomenon of cationic and anionic redox in a lithium-ion battery cathode will first be discussed. Our work in resonant inelastic X-ray scattering to investigate the spectroscopic features of O after delithiation has found potential "fingerprint" signals that could likely be used to identify and confirm reversible O redox if corroborated with other techniques. To follow, we will examine our work on Li-O 2 batteries. While our group and the research community have had many significant contributions and improvements to the field of Li-O 2 (such as decreasing overpotential and achieving cyclability in air environment), its practical application is still far from realization. Perhaps our most important contribution to this area is the discovery that Ir deposited on reduced graphene oxide can be used to halt the reduction of O 2 at the LiO 2 oxidation state. This not only significantly decreases the charge overpotential but also presents the important concept of oxidation-state controlled discharge. Subsequently, we will focus on our oxidation state-controlled redox-based charging of oxygen in a pure oxygen redox Li-ion battery. Future implications of this technology will be emphasized.

25 ENERGY STORAGE↗
DOE OSTI2021

Toward a First-Principles Framework for Predicting Collective Properties of Electrolytes

Conspectus Given the universal importance of electrolyte solutions, it is natural to expect that we have a nearly complete understanding of the fundamental properties of these solutions (e.g., the chemical potential) and that we can therefore explain, predict, and control the phenomena occurring in them. In fact, reality falls short of these expectations. But, recent advances in the simulation and modeling of electrolyte solutions indicate that it should soon be possible to make progress toward these goals. In this Account, we will discuss the use of first-principles interaction potentials based in quantum mechanics (QM) to enhance our understanding of electrolyte solutions. Specifically, we will focus on the use of quantum density functional theory (DFT) combined with molecular dynamics simulation (DFT-MD) as the foundation for our approach. The overarching concept is to understand and accurately reproduce the balance between local or short-ranged (SR) structural details and long-range (LR) correlations, allowing the prediction of the thermodynamics of both single ions in solution as well as the collective interactions characterized by activity/osmotic coefficients. In doing so, relevant collective motions and driving forces characterized by chemical potentials can be determined. Here, in this Account, we will make the case that understanding electrolyte solutions requires a faithful QM representation of the SR nature of the ion–ion, ion–water, and water–water interactions. However, the number of molecules that is required for collective behavior makes the direct application of high-level QM methods that contain the best SR physics untenable, making methods that balance accuracy and efficiency a practical goal. Alternatives such as continuum solvent models (CSMs) and empirically based classical molecular dynamics have been extensively employed to resolve this problem but without yet overcoming the fundamental issue of SR accuracy. We will demonstrate that accurately describing the SR interaction is imperative for predicting both intrinsic properties, namely, at infinite dilution, and collective properties of electrolyte solutions. DFT has played an important role in our understanding of condensed phase systems, e.g., bulk liquid water, the air–water interface, ions in bulk, and at the air–water interface. This approach holds huge promise to provide benchmark calculations of electrolyte solution properties that will allow for the development and improvement of more efficient methods, as well as an enhanced understanding of fundamental phenomena. However, the standard protocol using the generalized gradient approximation with van der Waals (vdW) correction requires improvement in order to achieve a high level of quantitative accuracy. Simply simulating with higher level DFT functionals may not be the best route considering the significant computational cost. Alternative methods of incorporating information from higher levels of QM should be explored; e.g., using force matching techniques on small clusters, where high level benchmark calculations are possible, to develop ideal correction terms to the DFT functional is a promising possibility. We argue that DFT with statistical mechanics is becoming an increasingly useful framework enabling the prediction of collective electrolyte properties.

Duignan, Timothy T.↗
DOE OSTI2020

Ensembles of Metastable States Govern Heterogeneous Catalysis on Dynamic Interfaces

Heterogeneous catalysis is at the heart of chemical industry. Being able to tune and design efficient catalysts for processes of interest is of an utmost importance, and for this, the molecular-level understanding of heterogeneous catalysts is the first step, and indeed a prime focus of the modern catalysis research. For a long time, a single most thermodynamically stable structure of the catalytic interface attained in reaction conditions had been envisioned as the reactive phase. However, some catalytic interfaces continue to undergo structural dynamics in the steady state, triggered by high temperatures, pressures, and binding and changing reagents. Among particularly dynamic interfaces are such widely-used catalysts as crystalline and amorphous surfaced supporting (sub-)nano metallic clusters. Recently, it became clear that this dynamic fluxionality causes the supported clusters to populate many distinct structural and stoichiometric states in catalytic conditions. Hence, the catalytic interface should be viewed as an evolving statistical ensemble of many (not one) structures. As such, every member in the ensemble contributes to the properties of the catalyst differently, and in proportion to its probability to be populated. This new notion flips the established paradigm and calls for new theory, modeling approaches, operando measurements, and updated design strategies. The statistical ensemble nature of surface-supported sub-nano cluster catalysts can be exemplified by oxide-supported and adsorbate-covered Pt, Pd, Cu, CuPd clusters, catalytic toward oxidative and non-oxidative dehydrogenation. They have access to a variety of 3D and quasi-2D shapes. The compositions of their thermal ensembles are dependent on the cluster size, leading to size-specific catalytic activities and the famous “every atom counts” phenomenon. The support and adsorbates affect catalyst structures, and state of the reacting species causes the ensemble to change in every reaction intermediate. The most stable member of the ensemble dominates the thermodynamic properties of the corresponding intermediate, whereas the kinetics can be determined by more active but less populated metastable catalyst states, and that suggests that many earlier studies might have overlooked the actual active sites. Both effects depend on the relative time-scales of catalyst restructuring and reaction dynamics. The catalyst may routinely operate off-equilibrium. Ensemble phenomena lead to surprising exceptions from established rules of catalysis, such as scaling relations, and Arrhenius behavior. Catalyst deactivation is also an ensemble property, and its extent of mitigation can be predicted through the new paradigm. These findings were enabled by advances in theory, such as global optimization and subsequent utilization of multiple local minima, and pathways sampling, as well as operando catalyst characterization. The fact that the per-site and per-species resolution is needed for the description and predicting of catalyst properties gives theory the central role in catalysis research, as most experiments provide ensemble-average information and cannot detect the crucial minority species that may be responsible for the catalytic activity.

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

A Six-Membered Concerted Mechanism for CO 2 Capture by Amines Studied under Charged Microdroplet Reaction Conditions

Carbon dioxide (CO 2 ) capture and storage represents an important technological challenge. A mechanistic understanding of interactions involved in the capture process is necessary not only for technological development but also for efficient conversion of captured CO 2 into value-added materials. Herein, we present a novel contained secondary electrospray ionization platform for studying the interactions of gaseous amines and CO 2 gas under microdroplet reaction conditions, which enables mass spectrometry (MS) characterization of CO 2 capture products and intermediates in real time. We detected [2 M + CO 2 + H] + species, which corresponds to a six-membered intermediate. DFT calculations confirmed the high stability of the protonated six-membered ring intermediate. This finding provides a plausible concerted mechanism in the microdroplet environment that excludes the involvement of thermodynamically disfavored ionic species. The carbamic acid counterpart of the final product/salt was readily characterized by tandem MS. The carbamic acid/amine salt was also isolated and characterized by Fourier transform infrared spectroscopy. By virtue of the fact that headspace vapors of amines are sampled, we were able to establish a high-throughput platform that enabled the CO 2 capture capacity of five different amines to be studied in under 2 min. The same device also enabled the absolute quantification of capture capacity.

Amines↗
DOE OSTI2021

A Catastrophic Charge Density Wave in BaFe 2 Al 9

Charge density waves (CDWs) are modulations of the electron density and the atomic lattice that develop in some crystalline materials at low temperatures. We report an unusual example of a CDW in BaFe 2 Al 9 below 100 K. In contrast to the canonical CDW phase transition, temperature-dependent physical properties of single crystals reveal a first-order phase transition. This is accompanied by a discontinuous change in the size of the crystal lattice. In fact, this large strain has catastrophic consequences for the crystals causing them to physically shatter. Single-crystal X-ray diffraction reveals superlattice peaks in the low-temperature phase signaling the development of a CDW lattice modulation. No similar low-temperature transitions are observed in BaFe 2 Al 9 . Electronic structure calculations provide one hint to the different behavior of these two compounds; the d-orbital states in the Fe compound are not completely filled. Iron compounds are renowned for their magnetism, and partly filled d-states play a key role. Finally, it is therefore surprising that BaFe 2 Al 9 develops a structural modulation at low temperature instead of magnetic order.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Understanding How Ligand Functionalization Influences CO2 and N2 Adsorption in a Sodalite Metal–Organic Framework

In this work, a detailed study is conducted to understand how ligand substitution influences the CO2 and N2 adsorption properties of two highly crystalline sodalite metal–organic frameworks (MOFs) known as Cu–BTT (BTT–3 = 1,3,5-benzenetristetrazolate) and Cu–BTTri (BTTri–3 = 1,3,5-benzenetristriazolate). The enthalpy of adsorption and observed adsorption capacities at a given pressure are significantly lower for Cu–BTTri compared to its tetrazole counterpart, Cu–BTT. In situ X-ray and neutron diffraction, which allow visualization of the CO2 and N2 binding sites on the internal surface of Cu–BTTri, provide insights into understanding the subtle differences. As expected, slightly elongated distances between the open Cu2+ sites and surface-bound CO2 in Cu–BTTri can be explained by the fact that the triazolate ligand is a better electron donor than the tetrazolate. The more pronounced Jahn–Teller effect in Cu–BTTri leads to weaker guest binding. The results of the aforementioned structural analysis were complemented by the prediction of the binding energies at each CO2 and N2 adsorption site by density functional theory calculations. In addition, variable temperature in situ diffraction measurements shed light on the fine structural changes of the framework and CO2 occupancies at different adsorption sites as a function of temperature. Finally, simulated breakthrough curves obtained for both sodalite MOFs demonstrate the materials’ potential performance in dry postcombustion CO2 capture. The simulation, which considers both framework uptake capacity and selectivity, predicts better separation performance for Cu–BTT. The information obtained in this work highlights how ligand substitution can influence adsorption properties and hence provides further insights into the material optimization for important separations.

adsorption capacities↗
DOE OSTI2022

Designing Sites in Heterogeneous Catalysis: Are We Reaching Selectivities Competitive With Those of Homogeneous Catalysts?

A critical review of different prominent nanotechnologies adapted to catalysis is provided, with focus on how they contribute to the improvement of selectivity in heterogeneous catalysis. Ways to modify catalytic sites range from the use of the reversible or irreversible adsorption of molecular modifiers to the immobilization or tethering of homogeneous catalysts and the development of well-defined catalytic sites on solid surfaces. The latter covers methods for the dispersion of single-atom sites within solid supports as well as the use of complex nanostructures, and it includes the post-modification of materials via processes such as silylation and atomic layer deposition. All these methodologies exhibit both advantages and limitations, but all offer new avenues for the design of catalysts for specific applications. Because of the high cost of most nanotechnologies and the fact that the resulting materials may exhibit limited thermal or chemical stability, they may be best aimed at improving the selective synthesis of high value-added chemicals, to be incorporated in organic synthesis schemes, but other applications are being explored as well to address problems in energy production, for instance, and to design greener chemical processes. Finally, the details of each of these approaches are discussed, and representative examples are provided. We conclude with some general remarks on the future of this field.

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

Prospects and Limitations of Predicting Fuel Ignition Properties from Low-Temperature Speciation Data

Using chemical kinetic modeling and statistical analysis, we investigate the possibility of correlating key chemical “markers”–typically small molecules–formed during very lean (φ ~ 0.001) oxidation experiments with near-stoichiometric (φ ~ 1) fuel ignition properties. One goal of this work is to evaluate the feasibility of designing a fuel-screening platform, based on small laboratory reactors that operate at low temperatures and use minimal fuel volume. Buras et al. [Combust. Flame2020,216, 472–484] have shown that convolutional neural net (CNN) fitting can be used to correlate first-stage ignition delay times (IDTs) with OH/HO2 measurements during very lean oxidation in low-T flow reactors with better than factor-of-2 accuracy. In this work, we test the limits of applying this correlation-based approach to predict the low-temperature heat release (LTHR) and total IDT, including the sensitivity of total IDT to the equivalence ratio, φ. We demonstrate that first-stage IDT can be reliably correlated with very lean oxidation measurements using compressed sensing (CS), which is simpler to implement than CNN fitting. LTHR can also be predicted via CS analysis, although the correlation quality is somewhat lower than for first-stage IDT. In contrast, the accuracy of total IDT prediction at φ = 1 is significantly lower (within a factor of 4 or worse). Furthermore, these results can be rationalized by the fact that the first-stage IDT and LTHR are primarily determined by low-temperature chemistry, whereas total IDT depends on low-, intermediate-, and high-temperature chemistry. Oxidation reactions are most important at low temperatures, and therefore, measurements of universal molecular markers of oxidation do not capture the full chemical complexity required to accurately predict the total IDT even at a single equivalence ratio. As a result, we find that φ-sensitivity of ignition delay cannot be predicted at all using solely correlation with lean low-T chemical speciation measurements.

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

Sand Migration Simulation during Gas Production from Gas Hydrate Reservoir at Kuparuk 7–11–12 site in the Prodhoe Bay Unit, Alaska

Uncontrolled sand production impedes continuous gas production from a hydrate reservoir as observed in the past field-scale gas production tests. Sand mobilization is strongly linked with sediment deformation and high pressure gradient. Sand production occurs when the mobilized sands reach the well. Throughout gas production from a hydrate reservoir, because of hydrate dissociation, deformation and pressure gradient change both in time and space. Therefore, it is necessary to consider the entire process to identify where sand is likely mobilized and how much mobilized sand could reach the well. Here, this study utilizes a coupled thermal, hydrological, chemical, and mechanical (thermo-hydro-chemo-mechanical) sand migration model to simulate a year long gas production from the hydrate reservoir at the Kuparuk 7–11–12 site in the Prodhoe Bay Unit, Alaska. It is found that sand would mainly come from the lower portion of the production zone where faster hydrate dissociation occurs. The relatively faster hydrate dissociation coupled with the fact that hydrate-bearing sediments and hydrate-free sediments have similar stiffness but different strengths causes complex stress transfer, which results in excessive shear deformation when the upper portion of the production zone starts to dissociate. This is evident not only near the well but also away from the well, leading to a large amount of sand mobilization. Another focus of the modeling study is evaluation of the effect of sand migration on gas production. The comparison of two extreme cases–with sandscreen and without sandscreen–suggests that migrated and settled solids around a sand control device can prevent pressure drop across the near-wellbore zone and lead to reduction of the gas production rate.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Comparison of the Chemical Composition of Liquids from the Pyrolysis and Hydrothermal Liquefaction of Lignocellulosic Materials

Major differences in thermal stability and hydrotreatment behavior of HTL and pyrolysis oils have been reported in the literature. However, little is known about the variations in the chemical composition of these oils that could explain such differences. Two commercial wood pyrolysis oils (Pyrovac and BTG), and their water-soluble (WS) and water-insoluble fractions (WIS) were analyzed and compared with the aqueous (WS WD-57 ) and oily (WIS WS-57 ) fractions obtained from hydrothermal liquefaction (HTL) of Douglas-fir. The samples were characterized by GC/MS, Karl Fischer titration, carbonyl content, total acid number, elemental composition, calorific value, proximate analysis, Fourier Transform Infrared Spectroscopy (FTIR), Folin-Ciocalteu (FC), and UV fluorescence. All the fractions were also analyzed by Fourier Transform Ion Cyclotron Resonance Mass Spectroscopy (FT-ICR-MS) and by Electrospray Ionization (EI). The most prevalent class of compounds in the water insoluble phases were phenols derived from lignin. Water-soluble phases contain mostly the oxygenated compounds derived from cellulose and hemicellulose and were richer in carbonyl functional groups. The water content of the resulting aqueous phases were between: 65 (WS BTG ) and 96 (WS Pyrovac ) wt. %. The bio-oil from BTG has higher water content and lower HHV, compared to Pyrovac oil. The GC/MS results of BTG oil show the presence of a more prominent acetic acid peak and higher TAN number than the Pyrovac oil. The GC/MS of Pyrovac oil showed more obvious mono-phenol peaks. The quantification of this family by Folin-Ciocalteu method confirmed higher content of monophenolic compounds compared with the BTG oil. The lower thermal stability of pyrolysis oils compared with HTL biocrudes can be partially explained by the fact that pyrolysis oils (BTG and Pyrovac) contain carbohydrates while HTL biocrude (WIS WD-57 ) doesn’t. Thus, we decided to further investigate the chemical differences between the phenolic rich fractions insoluble in water and the holocellulose derived compounds soluble in water. Even after water extraction, the acid content of the water insoluble fraction from BTG (WISBTG) was higher than the acid content of the water insoluble fraction obtained by HTL (WISWSD-57). Likewise, the acid content of the aqueous phases derived from pyrolysis oils (WS Pyrovac , WS BTG ) was also higher than for the aqueous phase obtained by HTL (WSWSD-57). This result is in part due to the use of bases in the HTL process that neutralizes the acid formed in that process. Moreover, the starting feedstock may also influence the differences between the oils. Although, the UV-Fluorescence spectra, ICR-MS and the EI analyses showed some minor differences in the molecular weight and chemical make-up of the oligomers soluble and insoluble in water from pyrolysis and HTL; the differences observed were not large enough to justify the differences in behavior between these oils reported in the literature. Our results suggest that the differences observed between HTL biocrudes and pyrolysis oils are likely partially due to the presence of holocellulose derived products in the pyrolysis oils and higher acid contents.

09 BIOMASS FUELS↗
DOE OSTI2021

Identification of Known and Novel Monomers for Poly(hydroxyurethanes) from Biobased Materials

Sourcing polymers from biobased materials is desirable for a transition to a more sustainable world. However, finding new monomers sourced from biobased molecules requires extensive experimental effort, due to the fact that most biomolecules lack the functional groups or properties needed to serve as a monomer in a particular polymeric system without chemical modification. In this work, we demonstrate a computational screening method to obtain monomers of a non-isocyanate polyurethane (NIPU), polyhydroxyurethane (PHU), from biomolecules. Our method generates candidate monomers from a list of 15 selected biobased starting molecules by performing chemical transformations that are standard to chemical catalysis and organic synthesis on these starting molecules and then screening the results for molecules with moieties that can be used as monomers for a PHU. Product molecules from two generations of reaction, a total number of 20,057, were analyzed automatically to identify molecules that have two or more epoxy or cyclic carbonate functional groups as candidates for further derivatization. These 21 candidate molecules were then subjected to another two generations of reaction, producing a set of 17,913 additional molecules. Further scrutiny of this set for desired functional groups characterizing PHUs and querying against the PubChem database was carried out, and 15 known molecules that are commercially available were identified, with the desirable features that they are predicted to be derived from biobased starting molecules in four or fewer steps. Moreover, molecule space generated from this process was primarily comprised of molecules that were not known to PubChem, indicating the framework also has the ability to generate novel candidate monomers that originate from biobased starting points.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A Simple Criterion for Feasibility of Heat Integration between Distillation Streams Based on Relative Volatilities

In a multicomponent distillation configuration, there are numerous sources and sinks of heat, and a potential way to reduce the heat duty requirement is to perform heat integration. Unfortunately, an algorithmic search of the optimal heat integration opportunities is intractable when the required temperatures of intermediate streams are computed via complex models. Instead, in this work, we introduce pressure-scaled pseudo relative volatility, a new metric to compare stream temperatures. Here, we justify the use of pseudo relative volatility by proving that this variable is a monotonically increasing function of the liquid fraction in a saturated mixture stream. Using this metric, we derive a shortcut criterion to check the feasibility of various heat integration opportunities, such as thermal coupling via heat transfer (TCH). The advantage of this approach is that it circumvents the need for explicit temperatures and instead relies on composition, component relative volatilities, and pressure—quantities that are readily available in shortcut models for optimization of distillation configurations. Leveraging this fact, we propose a new optimization framework that identifies feasible TCHs that we consider within the formulation while minimizing the total heat duty of a distillation configuration. We demonstrate, on a few examples, that our formulation can identify heat duty efficient configurations, some of which are multieffect configurations. Using this methodology, we discover configurations that are not only simpler than the fully thermally coupled (FTC) configuration but also have a much lower heat duty.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

Experimental and Theoretical Thermokinetic Analysis of the Na 3 VO 4 –CO 2 Reaction System under Different Physicochemical Conditions

Sodium vanadate (Na 3 VO 4 ) was synthesized, structural and microstructurally characterized as well as tested as possible carbon dioxide (CO 2 ) captor through thermodynamic calculations in addition to thermogravimetric dynamic and isothermal analyses. Ceramic structural characterization evidenced the formation of the Na 3 VO 4 crystal phase, while microstructural features evidenced dense agglomerates with a poor specific surface area, although synthesis temperature (600 °C) was not as high as that for other sodium ceramics. Na 3 VO 4 was investigated for the CO 2 capture process through dynamic and isothermal experiments, in the presence or absence of oxygen. All these experiments showed that Na 3 VO 4 possesses interesting CO 2 capture properties in a specific temperature range (520 and 600 °C). Isothermal products characterization allowed to elucidate the reaction pathway, implying the Na 4 V 2 O 7 formation, as possible reactive intermediate. In fact, reaction path evolution was supported with theoretical thermodynamic data. Moreover, the kinetic analysis, performed using the Avrami-Erofeev mathematical model and subsequent thermodynamic data obtention, probed the positive influence of the temperature, in both CO 2 superficial and bulk chemical sorption processes, while oxygen addition did not enhance the superficial process. Additionally, as it would be expected, kinetics were enhanced as a function of the CO 2 concentration.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗