Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computational pyrolysis”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Multi-Scale Lithofacies and Chemostratigraphic Analysis of Two Middle Devonian Marcellus Shale Wells in Northern West Virginia, USA

The Marcellus Shale Energy and Environmental Laboratory (MSEEL) consists of two project areas within the dry gas producing region of the Marcellus shale play in Monongalia County, West Virginia. MSEEL is a collaborative field project led by West Virginia University, with Northeast Natural Energy LLC, several industrial partners, and sponsored by the US Department of Energy National Energy Technology Laboratory. The study areas are drilled approximately 8.5 miles apart to better understand the vertical and lateral changes in stratigraphy over a short distance. Two vertical pilot wells, MIP-3H and Boggess 17H were drilled in the fall of 2015 and spring of 2019, respectively. Core was recovered from the MIP-3H (API: 47-061-01707-00-00) 112 feet (34m) between depths of 7445 to 7557 feet, and from the Boggess 17H (API: 47-061-01812-00-00) 139 feet (42m) between depths of 7908 and 8012 ft. A full suite of triple combo (gamma ray, neutron, density logs), image logs, and advanced logging tools were run in both wells and calibrated to core analysis. Core analysis includes medical computed tomography (CT) scans, mineralogy and chemostratigraphy determined from handheld X-Ray fluorescence (hhXRF) and X-Ray powder diffraction (XRD) measurements, and determination of total organic content (TOC). Lithofacies were determined at core-scale using traditional core description techniques and medical CT-scan images. Log-scale facies are based on mineralogy and TOC data and developed using petrophysical logging data calibrated to core data (XRD and pyrolysis data). Chemostratigraphic analysis utilized hhXRF data to determine the major and trace element trends in the cores. In the two wells six shale lithofacies were recognized at the core and log scale. Both wells show organic-rich facies (TOC > 6.5%) primarily in the middle and lower Marcellus, with a slight decrease in thickness of this interval in the Boggess 17H. This interval is interpreted as an increase in paleo-productivity (increased Ni, Zn, and V), decreased sedimentation (decreased detrital proxies), and anoxic to euxinic conditions (increased Mo and chalcophile elements). Paleo-redox conditions in both wells are dynamic throughout deposition transitioning between euxinic/anoxic to dysoxic/oxic. This trend is seen through elemental proxies and calcite/pyrite concretion distributions.

03 NATURAL GAS↗
DOE OSTI2021

Predicting Thermal Excursions During in situ Oxidative Regeneration of Packed Bed Catalytic Fast Pyrolysis Catalyst

Ex situ catalytic fast pyrolysis (CFP) uses a secondary reactor to upgrade biomass pyrolysis vapors to stabilized CFP oils with reduced oxygen content. In one configuration, the secondary reactor is operated as a packed-bed swing reactor system which allows coke-deactivated beds to be decarbonized in situ while other beds remain online for vapor upgrading. In situ decarbonization must be done carefully to avoid irreversible deactivation and/or physical degradation of catalyst pellets. Given that packed bed reactors are well known to have poor heat transfer characteristics, this is a critical issue impacting scaleability and commercial viability of the technology. To predict thermal excursions during regeneration, finite element computational models have been built to assist in scaling up oxidative decarbonization of a Pt/TiO2 CFP catalyst (0.5 mm spheres) from a bench scale packed bed with 100 g of catalyst to a pilot scale packed bed with 2 kg of catalyst and internal cooling tubes. Based on transient measurements of outlet temperature and effluent CO2 concentration, and using an assumed coke profile and activation energy, this paper demonstrates that specific combinations of effective thermal conductivity and wall heat transfer coefficient can fit bench scale oxidative regeneration data equally well. For the upscaled 2 kg bed, four bench-scale “best fit” parameter pairs give different predictions for location and magnitude of thermal excursions, with the maximum computed bed temperature gradients ranging from 30 °C cm-1 to as high as 3000 °C cm-1. The larger the fraction of heat removal by conduction through the cooling tubes, the greater the differences between the parameter pairs. The modelling results presented in this paper cast doubt on the industrial viability of the proposed combination of catalyst, bed and regeneration process, and point to the need for alternate reactor designs. However, there is considerable uncertainly in some of the key model parameters. The reliability of model predictions can be increased by adding more temperature measurements at key bed locations, testing additional variations in process conditions, performing careful bed dissections to determine the true coke profile, and perhaps most importantly, directly measuring the effective thermal conductivity of the catalyst pellets.

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

Association Kinetics for Perfluorinated n -Alkyl Radicals

Radical-radical reaction channels are important in the pyrolysis and oxidation chemistry of perfluoroalkyl substances (PFAS). In particular, unimolecular dissociation reactions within unbranched n-perfluoroalkyl chains, and their corresponding reverse barrierless association reactions, are expected to be significant contributors to the gas-phase thermal decomposition of families of species such as perfluorinated carboxylic acids and perfluorinated sulfonic acids. Unfortunately, experimental data for these reactions are scarce and uncertain. Furthermore, obtaining reliable theoretical predictions for such reactions is a laborious and computationally intensive task. Here, in this work, the chemical kinetics of the various association/decomposition reactions producing/decomposing the C 2 -C 4 series of unbranched n-perfluoroalkanes (C 2 F 6 , C 3 F 8 , and C 4 F 10 ) are examined using state-of-the-art ab initio transition-state-theory-based master-equation calculations. The variable-reaction-coordinate transition-state theory (VRC-TST) formalism is employed in computing the microcanonical and canonical rates for the association reactions. Reaction thermochemistry is obtained via composite quantum chemistry calculations and the laddering of error-canceling reaction schemes via a connectivity-based hierarchy approach employing ANL1/ANL0-style reference energies. Lennard-Jones collision model parameters for the considered systems were estimated by a direct dynamics approach, and collisional energy transfer parameters were obtained from analogies to systems of similar size and heavy-atom connectivity. A one-dimensional master equation approach was used to convert the microcanonical rate coefficients from the VRC-TST analysis into temperature- and pressure-dependent rate constants for the association reactions and the reverse dissociation reactions. The data are reported in standardized formats for usage in comprehensive chemical kinetic models for PFAS thermal destruction.

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

Spray Pyrolysis-Aerosol Deposition for the Production of Thick Yttria-Stabilized Zirconia Coatings

Aerosol deposition (AD) is a coating technique wherein particles are impacted onto a target substrate at reduced pressures, and supersonic particle impact velocities lead to coating consolidation. The limiting step in AD application is often not supersonic deposition operation, but aerosolization of powder particles with the proper size distribution; the translational impact velocity is strongly size-dependent. It is demonstrated that by directly synthesizing particles in the gas phase, size-controlled ceramic particles can be injected into AD systems. This in situ formation step obviates the need for particle aerosolization. Ultrasonic spray pyrolysis (USP) is applied to produce yttria-stabilized zirconia (YSZ), and USP is directly coupled with AD to produce consolidated, thick, YSZ coatings on metal substrates. USP-AD yields YSZ coatings on stainless steel and aluminum substrates with porosities <0.20, which grow to thicknesses beyond 100 μm. Aerodynamic particle spectrometry and electron microscopy reveal that the depositing particles are 200 nm–1.2 μm in diameter, though each particle is composed of nanocrystalline YSZ. Supporting computational fluid dynamics calculations demonstrate that the YSZ particle impact speeds are above 300 m s -1 . Thermal conductivity measurements demonstrate that USP-AD coatings have conductivities consistent with those produced from high-temperature processes.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A comprehensive experimental and kinetic modeling study of di-isobutylene isomers: Part 2

A wide variety of high temperature experimental data obtained in this study complement the data on the oxidation of the two di-isobutylene isomers presented in Part I and offers a basis for an extensive validation of the kinetic model developed in this study. Due to the increasing importance of unimolecular decomposition reactions in high-temperature combustion, we have investigated the di-isobutylene isomers in high dilution utilizing a pyrolysis microflow reactor and detected radical intermediates and stable products using vacuum ultraviolet (VUV) synchrotron radiation and photoelectron photoion coincidence (PEPICO) spectroscopy. Additional speciation data at oxidative conditions were also recorded utilizing a plug flow reactor at atmospheric pressure in the temperature range 725-1150 K at equivalence ratios of 1.0 and 3.0 and at residence times of 0.35 s and 0.22 s, respectively. Combustion products were analyzed using gas chromatography (GC) and mass spectrometry (MS). Ignition delay time measurements for di-isobutylene were performed at pressures of 15 and 30 bar at equivalence ratios of 0.5, 1.0, and 2.0 diluted in 'air' in the temperature range 900-1400 K using a high-pressure shock-tube facility. New measurements of the laminar burning velocities of di-isobutylene/air flames are also presented. The experiments were performed using the heat flux method at atmospheric pressure and initial temperatures of 298-358 K. Moreover, data consistency was assessed with the help of analysis of the temperature and pressure dependencies of laminar burning velocity measurements, which was interpreted using an empirical power-law expression. Electronic structure calculations were performed to compute the energy barriers to the formation of many of the product species formed. The predictions of the present mechanism were found to be in adequate agreement with the wide variety of experimental measurements performed.

09 BIOMASS FUELS↗
DOE OSTI2024

Understanding and controlling water-organic co-transport in amorphous microporous materials

The transport of molecules in microporous material is a significant and active area of research in separation applications. The movement of vapor/liquid molecules in a high-loading condition in microporous spaces is especially challenging to interrogate. This project aims to understand and control the transport of complex water and organic solvent mixtures in varying structures of microporous carbon molecular sieve (CMS) and activated carbon membranes. Modeling of transport behavior in such carbonaceous samples using computational modeling is difficult due to their amorphous structure. Therefore, this research intends to study the transport mechanism of water-organic mixtures by experimentally revealing fundamental transport properties, such as guest sorption amounts as well as diffusion and permeation rates. Considering the different size and guest-host affinity of water and organic solvent molecules, the objective of this research is to understand the structural conditions within the ultramicropores and micropores that generate different types of molecular transport and selection mechanisms within such complex systems. Carbonaceous materials are developed using tailored pyrolysis techniques to pyrolyze polymeric precursors, including polyvinylidene chloride (PVDC),polyvinylidene fluoride (PVDF), polymer of intrinsic microporosity (PIM) -1, and fully aromatic polyamide. The initial stage of the research will focus on the microscopic diffusion and sorption studies of pure component water, xylene isomers, and n,n- dimethylformamide within various carbons. Additional structural investigation on CMS, such as gas physisorption and neutron scattering studies, will be executed to gain deeper insight into these structure-transport relationships. Scanning electron microscopy and X-ray photoelectron microscopy are also used to further characterize these materials. The research will progress to probe mutual diffusion of complex water-organic mixtures within various CMS microstructures. Competitive sorption and permeation studies of the binary mixtures in different microstructural CMS will be followed. The obtained binary mixture transport parameters were used to predict binary water-organic mixtures transport.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Design and additive manufacturing of optimized electrodes for energy storage applications

Supercapacitors exhibit fast charging/discharging ability and have attracted considerable attention within the automotive, aerospace, and telecommunication industries. Porous carbons, prized for their high electrical conductivity and high surface area, have been attractive candidates for supercapacitor electrodes. Moving to thick electrodes is one strategy to further increase energy density due to a higher volume fraction of active material. However, thick electrodes suffer from sluggish charged species transport, which is why thin electrodes are currently favored. In this work, we investigate the use of computational optimization and additive manufacturing to design and fabricate thick porous electrodes with improved performance. Electrode performance was maximized by designing their morphologies via topology optimization and printing by projection micro stereolithography (PμSL) using commercial resin (PR48). The PR48 resin was then pyrolyzed (PR48-P) to create the final conductive electrode. The optimized PR48-P electrodes exhibited 99% improvement in capacitance compared to control electrodes printed with cubic lattice morphologies. To further improve performance, we formulated a resin combining graphene oxide (GO) and trimethylolpropane triacrylate (TMPTA). Electrodes printed with 3 wt% GO in TMPTA exhibited improved capacitance retention after pyrolysis compared to the PR48-P electrodes. Finally, this work demonstrates the benefits of using topology optimization to design electrodes and material development to improve functional properties of 3D printable electrodes.

25 ENERGY STORAGE↗
DOE OSTI2023

Atomistic insight into the effects of electrostatic fields on hydrocarbon reaction kinetics

Reactive Molecular Dynamics (MD) and Density Functional Theory (DFT) computations are performed to provide insight into the effects of external electrostatic fields on hydrocarbon reaction kinetics. By comparing the results from MD and DFT, the suitability of the MD method in modeling electrodynamics is first assessed. Results show that the electric field-induced polarization predicted by the MD charge equilibration method is in good agreement with various DFT charge partitioning schemes. Then, the effects of oriented external electric fields on the transition pathways of non-redox reactions are investigated. Results on the minimum energy path suggest that electric fields can cause catalysis or inhibition of oxidation reactions, whereas pyrolysis reactions are not affected due to the weaker electronegativity of the hydrogen and carbon atoms. MD simulations of isolated reactions show that the reaction kinetics is also affected by applied external Lorentz forces and interatomic Coulomb forces since they can increase or decrease the energy of collision depending on the molecular conformation. In addition, electric fields can affect the kinetics of polar species and force them to align in the direction of field lines. These effects are attributed to energy transfer via intermolecular collisions and stabilization under the external Lorentz force. The kinetics of apolar species is not significantly affected mainly due to the weak induced dipole moment even under strong electric fields. The dynamics and reaction rates of species are studied by means of large-scale combustion simulations of n-dodecane and oxygen mixtures. Results show that under strong electric fields, the fuel, oxidizer, and most product molecules experience translational and rotational acceleration mainly due to close charge transfer along with a reduction in their vibrational energy due to stabilization. This study will serve as a basis to improve the current methods used in MD and to develop novel methodologies for the modeling of macroscale reacting flows under external electrostatic fields.

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

Analysis of Infrastructures for Processing Plastic Waste using Pyrolysis-Based Chemical Upcycling Pathways

Modern mechanical recycling infrastructure for plastic is capable of processing only a small subset of waste plastics, reinforcing the need for parallel disposal methods such as landfilling and incineration. Emerging pyrolysis-based chemical technologies can "upcycle" plastic waste into high-value polymer and chemical products and process a broader range of waste plastics. In this work, we study the economic and environmental benefits of deploying an upcycling infrastructure in the continental United States for producing low-density polyethylene (LDPE) and polypropylene (PP) from post-consumer mixed plastic waste. Our analysis aims to determine the market size that the infrastructure can create, the degree of circularity that it can achieve, the prices for waste and derived products it can propagate, and the environmental benefits of diverting plastic waste from landfill and incineration facilities it can produce. We apply a computational framework that integrates techno-economic analysis, life cycle assessment, and value chain optimization. Our results demonstrate that the infrastructure generates an economy of nearly 20 billion USD and positive prices for plastic waste, opening opportunities for compensation to residents who provide plastic waste. Our analysis also indicates that the infrastructure can achieve a plastic-to-plastic degree of circularity of 34% and remains viable under various external factors (including technology efficiencies, capital investment budgets, and polymer market values). Finally, we present significant environmental benefits of upcycling over alternative landfill and incineration waste disposal methods, and comment on ongoing work expanding our modeling methodology to other chemical upcycling pathway case studies, including hydroformylation of specific plastics to chemicals.

Interdisciplinary↗
DOE OSTI2022

A Conical Intersection Influences the Ground State Rearrangement of Fulvene to Benzene

The rearrangement of fulvene to benzene is believed to play an important role in the formation of soot during hydrocarbon combustion. Previous work has identified two possible mechanisms for the rearrangement-a unimolecular path and a hydrogen-atom-assisted, bimolecular path. Computational results to date have suggested that the unimolecular mechanism faces a barrier of about 74 kcal/mol, which makes it unable to compete with the bimolecular mechanism under typical combustion conditions. This computed barrier is about 10 kcal/mol higher than the experimental value, which is an unusually large discrepancy for modern electronic structure theory. In the present work, we have reinvestigated the unimolecular mechanism computationally, and we have found a second transition state that is approximately 10 kcal/mol lower in energy than the previously identified one and, therefore, in excellent agreement with the experimental value. The existence of two transition states for the same rearrangement arises because there is a conical intersection between the two lowest singlet states which occurs in the vicinity of the reaction coordinates. The two possible paths around the cone on the lower adiabatic surface give rise to the two distinct saddle points. The lower barrier for the unimolecular mechanism now makes it competitive with the bimolecular one, according to our calculations. In support of this conclusion, we have reanalyzed some previous experimental results on anisole pyrolysis, which leads to benzene as a significant product and have shown that the unimolecular and bimolecular mechanisms for fulvene ? benzene must be occurring competitively in that system. Finally, we have identified that similar conical intersections arise during the isomerizations of benzofulvene and isobenzofulvene to naphthalene.

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

Determining the Role of Synthesis and Reaction Conditions on the Exposed Surface and Shape of ..beta..-Mo2C Catalysts

The production of drop-in quality biofuels typically can be achieved through catalytic fast pyrolysis (CFP) of biomass, followed by hydrodeoxygenation (HDO) of the vapor-phase product to reduce the oxygen content and stabilize the bio-oil product. To favor selectivity toward desired deoxygenation products, bifunctional catalysts with both acidic and hydrogenation sites are often employed. Due to both its oxophilicity and its ability to promote hydrogenation reactions, molybdenum carbide catalysts have shown promise for promoting the HDO of CFP vapors. However, the surface structure and composition of carbide catalyst particles are not well understood, precluding the establishment of fundamental structure-function relationships needed for the rational development of molybdedum carbide for HDO, as well as other reactions. For orthorhombic ..beta..-Mo2C, the (100) facet is largely studied for computational surface-reaction analyses, possibly due to its similarities with the close-packed (111) facet of noble-metal catalysts. However, it has previously been postulated through ab initio thermodynamic studies that the (111), (110), and (011) facets are most stable, though there is no consensus. Furthermore, while a majority of previous surface stability analyses of ..beta..-Mo2C focus on bulk terminations, the highly carburizing synthesis conditions and the use of ..beta..-Mo2C in carbon-containing reaction applications could introduce non-stoichiometric coverages of surface carbon. As such, this study uses ab initio thermodynamics calculations to assess the relative surface stability of the low-index facets of ..beta..-Mo2C, including consideration of varying the coverage of surface carbon. From this information, Wulff constructions are utilized to predict the synthesis and reaction condition-dependent shape and surface carbon composition of ..beta..-Mo2C catalyst particles.

BIOMASS FUELS,INORGANIC, ORGANIC, PHYSICAL, AND AN↗
DOE OSTI2021

Recent Advancements in Catalytic Fast Pyrolysis for the Production of Fuels and Chemicals from Biomass

Catalytic fast pyrolysis (CFP) is a promising biorefinery technology platform for generating both fuel blendstocks and chemical co-products. The versatility of this platform is afforded by (1) the coupling of a thermochemical deconstruction step to a catalytic step that enables tuning of the product slate and (2) the diversity of accessible products across gaseous, liquid, and solid streams. However, there are a number of challenges impeding commercial viability for this CFP platform, including high process costs caused by low carbon efficiency to bio-oil due to coking and light gas formation, the quality and properties of the bio-oil dictate processability for downstream operations, and the heterogeneity of the bio-oil limits quality of resulting fuel blendstocks. This presentation will discuss our efforts to address these challenges through a combined experimental-computational approach to catalyst design and process development, with a specific focus on (1) modifying product composition through catalyst design and (2) generating olefinic and oxygenated chemical co-products.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2021

Atomic-Level Features for Kinetic Monte Carlo Models of Complex Chemistry from Molecular Dynamics Simulations

The high computational cost of evaluating atomic interactions recently motivated the development of computationally inexpensive kinetic models, which can be parametrized from MD simulations of complex chemistry of thousands of species or other processes and accelerate the prediction of the chemical evolution by up to four order of magnitude. Such models go beyond the commonly employed potential energy surface fitting methods in that they are aimed purely at describing kinetic effects. So far, such kinetic models utilize molecular descriptions of reactions and have been constrained to only reproduce molecules previously observed in MD simulations. Therefore, these descriptions fail to predict the reactivity of unobserved molecules, for example in the case of large molecules or solids. In this work, we propose a new approach for the extraction of reaction mechanisms and reaction rates from MD simulations, namely the use of atomic-level features. Using the complex chemical network of hydrocarbon pyrolysis as example, it is demonstrated that kinetic models built using atomic features are able to explore chemical reaction pathways never observed in the MD simulations used to parametrize them, a critical feature to describe rare events. Atomic-level features are shown to construct reaction mechanisms and estimate reaction rates of unknown molecular species from elementary atomic events. Through comparisons of the model ability to extrapolate to longer simulation timescales and different chemical compositions than the ones used for parameterization, it is demonstrated that kinetic models employing atomic features retain the same level of accuracy and transferability as the use of features based on molecular species, while being more compact and parametrized with less data. We also find that atomic features can better describe the formation of large molecules enabling the simultaneous description of small molecules and condensed phases.

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

Theoretical studies of chemical reactions related to the formation and growth of polycyclic aromatic hydrocarbons (PAH) and molecular properties of their key intermediates (Final Progress Report)

The formation mechanisms of polycyclic aromatic hydrocarbons, (PAHs) – organic molecules carrying fused benzene rings – are of great interest to scientists and engineers due to their importance in combustion chemistry and astrochemistry. On Earth, PAHs are largely produced in incomplete combustion of fossil fuel and are considered as critical precursors to unwanted soot particles leading to combustion inefficiency and causing air pollution along with detrimental health effects. Simple PAH molecules initially formed in the gas phase, are further involved in a build-up process in combustion flames leading to larger PAH, bowl-shaped nanostructures, fullerenes, and solid-phase species including carbonaceous dust, graphene particles, and soot. In deep space, PAH and their derivatives are potential key intermediates and nucleation sites leading eventually to carbonaceous nanoparticles (“interstellar grains”). Therefore, the understanding of the key processes in the synthesis of PAHs along with their precursors and their degradation mechanisms in combustion systems and in interstellar, circumstellar, and planetary atmospheric environments will provide critical insights into how complex aromatic structures, carbonaceous nanoparticles, and fullerenes are formed and destroyed. Achieving this understanding is an important step in the development of the efficient combustion processes and of the ecofriendly devices with reduced environmental pollution as well as technological strategies for the production of hydrogen and solid carbon through thermal or plasma-assisted pyrolysis of natural gas and biomass. Also, the understanding of the key processes of PAH and soot growth will help in our comprehension of chemical evolution in the universe. Detailed information on the mechanisms and reliable rate constants of the key elementary chemical reactions involved in PAH formation and destruction processes and in inception of soot particles is often missing, with the main deficiencies being the absence of temperature- and pressure-dependent rate constants for the broad range of conditions occurring in various terrestrial and interstellar processes and the lack of data on the reaction products and their branching ratios. Complementary to experimental studies, these gaps in knowledge can be filled by using quantum chemical calculations of reaction potential energy surfaces providing us with accurate energies of reaction products, intermediates, and transition states, revealing the reaction mechanism, and giving the molecular properties required to compute rate constants for relevant reaction steps and product branching ratios using the RRKM-Master Equation (ME) method. Molecular dynamics (MD) simulations can be used in cases when a reaction rate cannot be properly described by statistical theories. During the terminal renewal project period we employed these ab initio/RRKM-ME and MD approaches to complete our studies on several key reactions relevant to the formation/growth of PAH and inception of soot particles including (1) the reaction mechanism and kinetics of the resonance stabilized fulvenallenyl radical with propargyl and C 3 H 4 isomers; (2) the reaction mechanism and kinetics for the C + indene and C 2 + styrene reactions producing naphthyl or azulenyl radicals in low-temperature environments; (3) the MD study of non-equilibrium dimerization of acepyrene and coronene and its radical. The information derived from our theoretical calculations contributed to a better fundamental understanding of the reaction mechanisms and provide missing critical kinetic data to improve combustion models of hydrocarbon fuels and astrochemical models of the growth of carbonaceous molecules and particles in cold molecular clouds, circumstellar envelopes, and planetary atmospheres.

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

Fulcrum House – Towards a Sustainable Composite Architecture

In the last decades, the accessibility to digital fabrication technologies positioned composites as an affordable, structurally efficient and low-maintenance alternative to traditional mineral-based and high-energy materials. Of special interest are carbon nanotubes (CNT) and carbon foam (CFoam) from methane and coal pyrolysis respectively, as an environmentally friendly substitute to glass and carbon fiber. This paper presents the design of a principled single-family house, as a first case study to identify the constraints of such a material system, including technology gap, techno-economic and life cycle assessment, as well as building code requirements. In turn, the investigation proposes a suitable design-to-fabrication workflow, which potentially could bring CNT and CFoam materials to the construction of buildings.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2025

Machine learning for reparameterization of multi-scale closures

Scientific machine learning (ML) is becoming increasingly useful in learning closure models for multi-scale physics problems; however, many ML approaches require a vast array of training data and can struggle with generalization and interpretability. Here, rather than learning an entire closure operator, we adopt an existing reduced-dimension model of the microphysics and learn an optimal re-parameterization of the solver. We demonstrate two approaches for training the reduced dimension closure model (1) an a priori method that optimizes the closure parameterization and the neural network parameters separately and (2) an a posteriori method that simultaneously optimizes both. Using the simulation of biomass pyrolysis as a motivating example, we show that the a posteriori method achieves better target losses and is less dependent on training dataset size for generalizability. We then demonstrate the impact that implementing this reparameterization has at the macroscale, showing improved predictive performance with no modification to the underlying macroscale solvers.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Rapid data-driven model reduction of nonlinear dynamical systems including chemical reaction networks using ℓ 1 -regularization

We develop a new data-driven paradigm for efficient model reduction of a broad class of nonlinear dynamical systems. Our model reduction method directly enables the interpretation of key components of the dynamical system, unlike traditional projection-based model reduction methods that focus on reducing computational complexity more than interpretability. Our method is not application specific and is simple to implement on nonlinear dynamical systems arising from a variety of different fields. It requires minimal parameterization using a single parameter to trade-off between model complexity and estimation error. We use a data-driven paradigm to formulate model reduction as an efficient convex optimization problem that scales polynomially in the original size of the complex system, enabling systems with as many as thousands of components to be reduced in a matter of minutes.

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

Simulation of Transport Phenomena in Bubble Column Reactors

The production of carbon and hydrogen from methane is a viable pathway for monetizing the vast reserves of US natural gas. Liquid (metals and/or salts) bubble column reactors operating at temperatures greater than 1000° C allow the conversion of methane into hydrogen without CO 2 emission. The low-density carbon floats on top of the liquid surface, allowing an easier removal, compared to other reactor concepts for methane pyrolysis. In this current report, we present results from simulations of transport phenomena in a liquid bubble column reactor. The goal is to investigate and understand parameters that are crucial for bubble column reactor design and scale up. The report presents the formulation and numerical methods used in the multiphase direct numerical simulation software, Quilt. This software was validated for the bubbly flow using numerical and physics based test problems. Quilt was then used to perform simulations of the bubble column reactor as a parametric study, varying the bubble injection parameters. The simulations were visualized for qualitative analysis and followed up with a quantitative analysis. The effect of injection parameters on various quantitative measures of the bubble residence time, interfacial area and motion through the column are presented.

03 NATURAL GAS↗