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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Recommended DSMC Collision Model Parameters for Planetary Entry

Hybrid simulation frameworks combining Computational Fluid Dynamics (CFD) and Direct Simulation Monte Carlo (DSMC) are frequently employed to efficiently perform high-fidelity simulations of environments consisting of both continuum and rarified flow. DSMC is a stochastic, particle-based method which solves the fundamental Boltzmann equation and is therefore necessary for high-Knudsen flow where continuum-based assumptions governing CFD break down. However, the DSMC methodology is generally computationally inefficient to model the continuum regime. In a CFD/DSMC hybrid approach, obtaining an accurate, high-fidelity solution hinges on the consistent treatment of transport properties and the thermo-chemical models employed within the two solvers. In principle, in regions where CFD and DSMC are both employed, the same gas mixture under the same conditions should have the same properties, regardless of simulation type. Observed differences should be due to non-equilibrium processes, rather than differences in physical models. The goal of this work is to provide a comprehensive DSMC transport database encompassing the vast majority of known gas species encountered during Earth or other planetary atmospheric entry. This goal is largely possible due to recently performed ab initio quantum chemistry calculations. Combined with other high-fidelity data, the planned database will consist of collision integral data for over 200 neutral and ionized species and over 20000 binary collisions. From these collision integrals, Nelder-Mead optimization is used to compute collision-specific Variable Soft Sphere (VSS) collision model parameters, fitted from 300 K to 20000 K. Initial comparisons of transport properties of relevant equilibrium gas mixtures show great agreement between CFD and DSMC-derived results. The completed database can be readily applied to model binary collisions of any gas mixture containing the included species over the specified temperature range, making it a valuable tool for future planetary probe modeling efforts. An example is shown below in Fig. 1. Equilibrium mixture transport properties for a 35-species mixture composed originally of 10% air and 90% pyrolysis species of a carbon-phenolic ablator material [4] are computed using both fitted VSS Parameters and the original CFD collision integral values. Deviations in computed properties between the two approaches are less than 5% for the entire temperature range.

M. R. Gosma↗
NASA NTRSDate not supplied

Analysis of Mars 2020 Entry with Coupled Material Response and CFD

When computing the response of an ablating thermal protection system during entry, the aerothermal environment and material response are traditionally computed independently, with the introduction of a blowing correction term in the material response model to account for the outgassing of char and pyrolysis gases. This study presents an approach, where the pyrolysis blowing gases, calculated within the PATO material response code [1], are integrated into the DPLR hypersonic CFD code [2] via a blowing boundary condition. This methodology employs an iterative process, whereby the blown pyrolysis gas products from PATO are incorporated into DPLR, refining surface heating predictions. The NEQAIR code is used to compute radiative heating [3]. The method is applied in computing 3D material response of the Mars 2020 entry. When compared with uncoupled material response, the coupled simulations show a lower surface heat flux initially and a higher heat flux at peak heating as shown in Fig. 1. This is in agreement with previous work using a sphere case with environments from MSL [4].

Heat Transfer↗
NASA NTRSDate not supplied

Analysis of Mars 2020 Entry with Coupled Material Response and CFD

When computing the response of an ablating thermal protection system during entry, the aerothermal environment and material response are traditionally computed independently, with the introduction of a blowing correction term in the material response model to account for the outgassing of char and pyrolysis gases. This study presents an approach, where the pyrolysis blowing gases, calculated within the PATO material response code, are integrated into the DPLR hypersonic CFD code via a blowing boundary condition. This methodology employs an iterative process, whereby the blown pyrolysis gas products from PATO are incorporated into DPLR, refining surface heating predictions. The NEQAIR code is used to compute radiative heating. The method is applied in computing 3D material response of the Mars 2020 entry. When compared with uncoupled material response, the coupled simulations show a lower surface heat flux initially and a higher heat flux at peak heating as shown in Fig. 1. This is in agreement with previous work using a sphere case with environments from MSL.

Heat Transfer↗
DOE OSTI2022

Catalysts by pyrolysis: Direct observation of transformations during re-pyrolysis of transition metal-nitrogen-carbon materials leading to state-of-the-art platinum group metal-free electrocatalyst

Transition metal-nitrogen-carbon (M-N-C) materials have been the focus of scientists’ efforts to address the rising need for earth-abundant materials solutions for energy technology and decarbonization of the economy. They are viewed as one of the most promising candidates to replace platinum group metal (PGM) catalysts in the fuel cell and energy conversion fields, including the application of oxygen reduction reaction, carbon dioxide reduction reaction, and nitrogen reduction reaction. In the effort to improve M-N-C materials properties and achieve atomic dispersity of the transition metal in the carbonaceous matrix, a re-pyrolysis process has been proposed. This secondary heat treatment process of already obtained primary pyrolysis-derived M-N-C materials has been widely reported to substantially improve the electrochemical performance and operational stability of the catalysts. Here, we report a systematic investigation of this process used on samples of templated M-N-C catalysts to obtain state-of-the-art catalysts via in situ heating X-ray photoelectron spectroscopy (XPS), scanning transmission electron microscopy (STEM), energy-dispersive X-ray spectroscopy (EDS), electron energy loss spectroscopy (EELS), X-ray diffraction (XRD), and X-ray computed tomography (CT) characterization methods. It is found that the re-pyrolysis of M-N-C materials could result in the partial amorphization of the carbonaceous substrate. It causes the rearrangement and transformation of multitudinous N moieties, leading to optimization of their morphological display and association with atomically dispersed transition metal dopants. Ultimately, the re-pyrolysis results in an increase in uniformity of the active Fe-Nx sites distribution without the formation of nano-crystalline phases (metallic or carbide) and with overall preservation of the morphology of the carbonaceous framework achieved during the first formative pyrolysis step of the templated synthesis. Finally, these observations provide confirmation that empirically established re-pyrolysis is recommended to be used on all M-N-C materials despite the different synthesis routes to obtain a practical advanced catalytic material.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Open source implementation of glued sphere discrete element method and nonspherical biomass fast pyrolysis simulation

Here, in this research, a glued-sphere Discrete Element Method (DEM) method was implemented in the open-source, Computational Fluid Dynamics (CFD) software MFiX. The implementation was verified using a cylinder-wall collision and then validated by simulating the packing and fluidization of non-spherical particles. The validated code was applied to simulate fast pyrolysis of non-spherical biomass particles in a bubbling fluidized bed. The glued sphere occupancy ratio was proposed to quantify the quality of shape resolution using glued sphere. Shape resolution showed significant influence on the packing height in the simulation of particle packing and an occupancy ratio of 80 % was recommended. Its influence is minor in fully fluidized bed but can be 8 times higher in fixed fluidized bed. Three tested heat transfer models predicted similar yields of elongated biomass fast pyrolysis. The solver developed in this research can be used to simulate other multiphase reacting flows involving non-spherical particles.

42 ENGINEERING↗
NASA NTRSDate not supplied

Material Characterization and Modeling of Room Temperature Vulcanizing Silicone

Room Temperature Vulcanizing silicone (RTV) is a high-temperature adhesive that has successfully been used as a gap-filler between Thermal Protection System (TPS) tiles for heatshields on numerous missions. It is also used to bond instrumentation plugs such as temperature and pressure sensors into the heatshields. While RTV has been traditionally assumed to be a non-porous and non-ablating material, numerous experiments have shown that RTV pyrolyzes and becomes highly porous as it is heated. Heating RTV has also shown swelling, or intumescence, which can pose unique problems that lead to roughness induced boundary-layer transition, surface oxide formation and contamination of heat shield sensors. Therefore, it is crucial to understand and model the intumescence phenomenon of RTV. As data for RTV material properties is limited, the first step in modeling RTV is to collect material properties such as pyrolysis mass-loss, microstructure change, virgin and char porosity, etc. which was performed in our initial study. Additionally, thermomechanical properties such as Young’s modulus and Poisson ratio are required for modeling the intumescence of RTV, which were taken from literature and the coefficient of thermal expansion was collected using in-situ heating and Micro Computed Tomography (µ-CT) in previous studies. Finally, numerous other properties such as pyrolysis gas properties, virgin and char thermal conductivity and specific heat were compiled from previous experiments and literature into a material database that can be used for simulations. In Porous Material Analysis Toolbox based on OpenFOAM (PATO) [4], structural mechanics coupled with material response was used for simulating the intumescence of RTV as it is heated. However, since the permeability of the material is very low, the pyrolysis gas creates an internal pressure build-up as the material is being heated, significantly contributing to the deformation of the material. To correctly characterize this phenomenon, additional physics models were implemented into PATO's stress analysis solver, and results were compared with RTV dilatometry test data as a preliminary verification case. Future work will include experiments of RTV at the Plasmatron X facility and the in-situ heating cell with µ-CT, and improvement of simulation tools to more accurately model RTV intumescence.

TPS↗
NASA NTRSDate not supplied

Material Properties and Modeling of Room Temperature Vulcanizing Silicone

Room Temperature Vulcanizing silicone (RTV) is a high-temperature adhesive that has successfully been used as a gap-filler between Thermal Protection System (TPS) tiles for heatshields on numerous missions. It is also used to bond instrumentation plugs such as temperature and pressure sensors into the heatshields. While RTV has been traditionally assumed to be a non-porous and non-ablating material, numerous experiments have shown that RTV pyrolyzes and becomes highly porous as it is heated. Heating RTV has also shown swelling, or intumescence, which can pose unique problems that lead to roughness induced boundary-layer transition, surface oxide formation and contamination of heat shield sensors. Therefore, it is crucial to understand and model the intumescence phenomenon of RTV. As data for RTV material properties is limited, the first step in modeling RTV is to collect material properties such as pyrolysis mass-loss, microstructure change, virgin and char porosity, etc. which was performed in our initial study. Additionally, thermomechanical properties such as Young’s modulus and Poisson ratio are required for modeling the intumescence of RTV, which were taken from literature and the coefficient of thermal expansion was collected using in-situ heating and Micro Computed Tomography (µ-CT) in previous studies. Finally, numerous other properties such as pyrolysis gas properties, virgin and char thermal conductivity and specific heat were compiled from previous experiments and literature into a material database that can be used for simulations. In Porous Material Analysis Toolbox based on OpenFOAM (PATO) [4], structural mechanics coupled with material response was used for simulating the intumescence of RTV as it is heated. However, since the permeability of the material is very low, the pyrolysis gas creates an internal pressure build-up as the material is being heated, significantly contributing to the deformation of the material. To correctly characterize this phenomenon, additional physics models were implemented into PATO's stress analysis solver, and results were compared with RTV dilatometry test data as a preliminary verification case. Future work will include experiments of RTV at the Plasmatron X facility and the in-situ heating cell with µ-CT, and improvement of simulation tools to more accurately model RTV intumescence.

PATO↗
DOE OSTI2022

Mechanistic insights into metal, nitrogen doped carbon catalysts for oxygen reduction: progress in computational modeling

We report metal and nitrogen doped carbon materials (denoted as M–N–C) synthesized through high-temperature pyrolysis have been found to exhibit activity for oxygen reduction reaction (ORR) approaching that of Pt and electrochemical stability higher than previous MN 4 -containing macrocyclic molecular catalysts. Tremendous efforts have thus been devoted to the advancement of M–N–C catalysts as an economical alternative to Pt-based catalysts for proton exchange membrane fuel cell cathodes with a focus on simultaneously improving activity and stability. To this end, novel computational modeling techniques have been developed and applied to acquire knowledge crucial for accelerating the pace of M–N–C catalyst development. In this review, recent progress in computational method development, as well as the predictions of chemical structure of active sites, reaction pathways, ORR kinetics, and catalyst stability in electrochemical environments, are critically surveyed. Moreover, the crucial role of computational modeling to elucidate the functional mechanism of M–N–C catalysts for ORR in acid media and enable rational design of M–N–C catalysts is discussed with a visionary outlook for the field.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Experimental and Theoretical Study of Oxolan-3-one Thermal Decomposition

The thermal decomposition of oxolan-3-one, a common component of the bio-oil formed during biomass pyrolysis, has been studied using ab initio calculations and experiments employing pulsed gas-phase pyrolysis with matrix isolation FTIR product detection. Four pathways for unimolecular decomposition were predicted using computational methods. The dominant reaction channel led to carbon monoxide, formaldehyde, and ethylene, all of which were observed experimentally. Here, the other channels led to an assortment of products including ketene, water, propyne, and acetylene, which were all confirmed in the matrix-isolation FTIR spectra. There is also evidence for the production of substituted ketenes in pyrolysis, most likely hydroxyketene and methylketene.

09 BIOMASS FUELS↗
DOE OSTI2024

Mathematical Modeling of Multi-Phenomena Anisotropic Systems: Ejection of Primary Aerosols during the Fast Pyrolysis of Biomass

This study introduces a novel particle model for biomass fast pyrolysis, incorporating an anisotropic cylindrical particle to address mass and energy transport coupled with aerosol ejection, which previous models have overlooked. The main contribution lies in developing a model that considers aerosol generation in anisotropic cylindrical particles for the first time, addressing bubbling dynamics and bursting within the liquid phase. The population balance equation describes bubble dynamics and aerosol formation, capturing phenomena like nucleation, growth, coalescence, and bursting. The model employs the method of moments with bubble volume as an internal variable, substantially reducing computational costs by eliminating dependence on this variable. Results highlight the significant impact of anisotropy and particle size on aerosol ejection: smaller, less elongated particles experience faster heating, quicker conversion, and the increased accumulation of the liquid intermediate phase. Specifically, 1 mm diameter particles yield higher concentrations of metaplast and bio-oil aerosols, exceeding 15%, compared to concentrations below 11% for 3 mm particles. This model provides insights into aerosol structure (volume, surface area), aiding in understanding aerosol reactivity at the reactor scale.

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

Reinforcement induced microcracking during the conversion of polymer-derived ceramics

This study investigates the role of discontinuous reinforcement on the microcracking of a preceramic polymer matrix during polymer to ceramic conversion. The material system comprised a carbosilane-based, preceramic polymer reinforced with mullite particles. The preceramic resin slurry was formulated for photopolymerization on a digital light projection 3D printer. Micro X-ray computed tomography, using a synchrotron light source, monitored a printed composite part during pyrolysis. The conversion of the carbosilane matrix into Si(O)C generated microcracks in the matrix that radially extended from particles. The likelihood of microcrack formation positively correlated with particle size. The largest particles (>10 4 µm 3 volume or >60 µm side length) always abutted microcracks, while the matrix surrounding smaller particles (<5 µm) was free of microcracks. Numerical calculations showed the particles resist the shrinking matrix during its conversion. This created tensile stresses within the matrix. The driving force for microcrack growth was also analyzed with regard to the influence of matrix material parameters, particle form-factor and crack length. These results reveal the need to consider the form factor of discontinuous reinforcements, with regard to the material properties of the polymer-derived ceramic, in order to reduce or eliminate defects in the final part.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Bridging particle and reactor scales in the simulation of biomass fast pyrolysis by coupling particle resolved simulation and coarse grained CFD-DEM

Efficient utilization of biomass is challenging due to the complexity of the feedstock at molecular, particle, and reactor scales. Direct coupling of particle scale model and reactor scale simulation is computationally infeasible. In this work, the results from a particle resolved simulation were used to calibrate the heat transfer and reaction kinetics of particles. A hybrid drag model was used to consider the different fluidization properties of sands and biomass particles. The computation cost was reduced using a coarse-grained Discrete Element Method (DEM). The simulated conversion rate, char yield, tar yield, and pyrolysis gas yield compare well with experiment. The residence time and axial distributions of biomass were also analyzed. This novel multi-scale method provides an efficient and accurate tool for the modeling of biomass pyrolysis reactors.

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

Dense Ru single-atoms integrated with sulfoacids for cellulose valorization to isosorbide

Metal-acid bifunctional catalysts are the unity of two opposites (metal-acid repelling) for cellulosic biomass valorization to isosorbide. It is challenging to improve the selectivity of metal hydrogenation catalysts and their synergism with acids for catalytic hydrolysis and dehydration to achieve substantial isosorbide. Herein, dense Ru single-atoms (10.1 wt% of Ru SAs) are anchored on sulfoacid-functionalized hollow mesoporous carbon shells, designed by assembling silica and 8-hydroxyquinoline-modified chitosan (HQ-CTS) through in situ Stöber templating strategy before pyrolysis and acid treatment. Based on X-ray absorption fine structure and computational modeling results, the structure of Ru SAs is determined as RuN 4 , which is more selective for a transitional glucose hydrogenation to sorbitol than Ru001 of Ru clusters. A lower-energy barrier of 1.21 (0.72) eV is delivered over RuN 4 (Ru001) for glucose hydrogenation (isomerization). These Ru SAs are integrated with sulfoacids (SO 3 H) but resistent against acids, rendering enhanced isosorbide yield in water as compared to Ru clusters, via a one-pot cascade reaction under harsh conditions (220 °C, 6 MPa H 2 ). The elaborately fabricated dense Ru SAs and sulfoacids, achieved by varying the addition time of HQ-CTS during the in situ Stöber templating process, improve the synergism of glucose hydrogenation with cellulose hydrolysis and sorbitol dehydration. In conclusion, this study provides a new idea for rational design of high-performance metal-acid bifunctional catalysts toward one-pot conversion of cellulose to isosorbide.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Simulations of biomass pyrolysis using glued-sphere CFD-DEM with 3-D intra-particle models

The shape of biomass particles influences intra-particle heat transfer, particle-particle collisions, interphase momentum, heat transfer, and pyrolysis yield. In this research, the effects of these multi-scale phenomena are captured using a glued-sphere computational fluid dynamic discrete element method (CFD-DEM) with 3-D intra-particle models. The glued-spheres resolved the shapes of particles and intra-particle temperature and species distributions. The implementation was validated using the pyrolysis data of a spherical particle and a cylindrical particle. Then, the influences of biomass shapes and sizes on pyrolysis were investigated. The results revealed the significant influence of surface areas and mass distributions on interphase heat transfer and chemical reactions. The intra-particle heat transfer is neglectable for a Biot number smaller than 0.41. Finally, the pyrolysis of irregularly shaped biomass particles in a fixed bed, fluidized bed, and spouted bed was simulated. Analysis of conversion distributions showed the largest deviations in the fixed bed and the slowest conversion rate in the spouted bed. A good balance between uniform and fast conversion is achieved in the fluidized bed. The multi-scale method proposed in this research provides a tool for the evaluation of large irregular-shaped biomass particle pyrolysis in different types of reactors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1981

Experimental determination of ablation vapor species from carbon phenolic heat-shield materials

The relative concentrations of vapors produced from carbon phenolic composites under thermal loadings approximating those expected at peak heating during vehicle entry into the atmospheres of the outer planets have been determined. The technique of vaporizing the surface of bulk samples by laser irradiation while measuring in situ the vapor species by mass spectrometry is described. Results show that vapor composition varies with irradiance level and with depth of heating (or extent of pyrolysis). Attempts are made to compare these experimental results with the theoretical predictions from computer codes.

Lincoln, K. A.↗
NASA NTRS1975

A users manual for a revised version of the Langley charring ablator program

A computer program is described that will compute the transient response of a thermal protection material to a prescribed heat input at the surface. The program has the capability of analyzing pyrolysis gas chemical kinetics in detail and treating pyrolysis reactions-in-depth. Deposition of solid products produced by chemical reactions in the gas phase is included in the analysis. An outline is given for the theory. detailed operating instructions for the computer program are included.

Stroud, C. W.↗
DOE OSTI2023

A computational fluid dynamics model to estimate local quantities in firebrand char oxidation

Firebrand burning is a complex phenomenon that is influenced by several parameters which are difficult to fully explore experimentally. Computational fluid dynamics models capable of predicting local quantities are essential for accurate prediction of char oxidation in firebrands. This article presents a computational fluid dynamics model to estimate firebrand mass loss, diameter change, and surface temperature during char oxidation. The model was validated using previously conducted wind tunnel experiments. These experiments were conducted for firebrands of two different aspect ratios, which were arranged in three different configurations (single, horizontal array, and vertical array), and for four different wind speeds (0.5, 1, 1.5, and 2 m/s). The computational fluid dynamics results were compared with a previous 1 D model. In all the test cases, the computational fluid dynamics model predicted the physical phenomena with significantly improved accuracy compared to a 1 D model. The char oxidation model presented in this article can be coupled with other models to study firebrand generation and trajectory, biomass pyrolysis, fluidized bed reactors, and coal combustion.

Engineering↗