Search NASA⌕ Search

SEARCH · Search NASA

Results for “Kinetic methods”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Sulfur-functionalized solid-phase materials for the selective separation of arsenic and selenium

Radioactive arsenic (As) isotopes are of growing interest for applications in nuclear medicine, national security, and environmental research. Recent efforts at the Facility for Rare Isotope Beams (FRIB) have focused on aqueous harvesting of selenium-72,73 ( 72,73 Se) and their daughter isotopes, arsenic-72,73 ( 72,73 As), which are particularly valuable for medical applications and nuclear data studies, respectively. Both conventional isotope production and harvesting methods require chemical separations to purify radioactive As from parent and co-produced Se radioisotopes. While several solid-phase separation methods for As and Se exist, many depend on complex oxidation state control or highly acidic conditions. This study presents results for sulfur-based solid-phase materials selected to enable uptake at lower acidity and eliminate the need for intricate redox chemistry. Specifically, the performance of three covalently bound sulfur-based ligands were evaluated: (1) thiophenol-polystyrene, (2) propanethiol-silica, and (3) thiourea-silica. Uptake characteristics—including distribution coefficients (Dw), kinetics, and column separation behavior—were assessed using 75 Se and 73 As in hydrochloric (HCl) acid and nitric (HNO 3 ) solutions. The resins demonstrated high-yield (>95%) and high-purity As recovery across a range of HCl concentrations. Comparable results in HNO 3 were achieved when combined with anion exchange chromatography. Furthermore, the potential application of these materials for medical isotope generators was also investigated through ligand stability and repeated elution studies. Overall, sulfur leaching from the resins was negligible at the concentrations relevant for these separations but increased with higher acid concentrations.

Arsenic↗
DOE OSTI2025

Carbonation of Alkaline Earth Metal Hydroxides: Structure across Nano- to Mesoscales

The conversion of gaseous CO 2 into a solid constitution through mineralization is an active area of carbon capture, and alkaline earth metal hydroxides (M(OH) 2 , M = Ca 2+ , Mg 2+ ) are frontrunners in this area. As model systems, nanolime samples are excellent templates for the study of this reaction. Here, we have examined these under ambient pressure conditions with controlled humidity and CO 2 . Utilizing a broad range of analytical methods, we first established the purity and structures of the selected materials. We then examined the structural changes due to carbonation, using infrared spectroscopy, X-ray scattering, and neutron scattering. The resulting structural changes are resolved from nanoscale to mesoscale and from early-stage to late-stage carbonation. Ca(OH) 2 and Mg(OH) 2 are found to behave quite differently. As expected from prior work, the carbonation of Ca(OH) 2 is kinetically favored. Further, our structure studies suggest this is due to a facile reaction at the fractal interface of the particles. From early- to late-stage carbonation, there is a consistent increase in the fractal roughness. This is in contrast to Mg(OH) 2 where the same surface evolves into a smooth conformal coating. For this material the major reacting component is at the mesoscale, suggesting globular particle growth or evolving macro-porosity. Because neutron scattering is sensitive to hydrogen content, we expected a significant change as M(OH) 2 evolves to MCO 3 . Such a change is found for Ca(OH) 2 but not for Mg(OH) 2 , providing evidence for the formation of hydrated carbonates for the later material. The formation of a conformal layer along with water-rich carbonate formation is an impediment to the use of Mg(OH) 2 for carbon capture. For energy-efficient carbon capture, it would be desirable to enhance carbonation rates for Mg(OH) 2 , and one possible route would be the use of anhydrous fluids.

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

Blueprinting Electrified Transit System Implementation

To achieve a more affordable and reliable transportation system, we need to smartly upgrade our power systems and install a large number of charging stations, but conventional planning methods are not up to the task. By applying advanced simulation and optimization tools, we can design a smarter, more cost-effective electric transportation network. The initial focus was on public transit systems, demonstrating how this approach can deliver broader economic, reliability, and air quality benefits nationwide.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

A Local Macroscopic Conservative (LoMaC) Low Rank Tensor Method for the Vlasov Dynamics

Abstract In this paper, we propose a novel Local Macroscopic Conservative (LoMaC) low rank tensor method for simulating the Vlasov-Poisson (VP) system. The LoMaC property refers to the exact local conservation of macroscopic mass, momentum and energy at the discrete level. This is a follow-up work of our previous development of a conservative low rank tensor approach for Vlasov dynamics ( arXiv:2201.10397 ). In that work, we applied a low rank tensor method with a conservative singular value decomposition to the high dimensional VP system to mitigate the curse of dimensionality, while maintaining the local conservation of mass and momentum. However, energy conservation is not guaranteed, which is a critical property to avoid unphysical plasma self-heating or cooling. The new ingredient in the LoMaC low rank tensor algorithm is that we simultaneously evolve the macroscopic conservation laws of mass, momentum and energy using a flux-difference form with kinetic flux vector splitting; then the LoMaC property is realized by projecting the low rank kinetic solution onto a subspace that shares the same macroscopic observables by a conservative orthogonal projection. The algorithm is extended to the high dimensional problems by hierarchical Tuck decomposition of solution tensors and a corresponding conservative projection algorithm. Extensive numerical tests on the VP system are showcased for the algorithm’s efficacy.

Guo, Wei↗
DOE OSTI2024

Relativistic Douglas–Kroll–Hess calculations of hyperfine interactions within first-principles multireference methods

A relativistic magnetic hyperfine interaction Hamiltonian based on the Douglas–Kroll–Hess (DKH) theory up to the second order is implemented within the ab initio multireference methods, including spin–orbit coupling in the Molcas/OpenMolcas package. This implementation is applied to calculate relativistic hyperfine coupling (HFC) parameters for atomic systems and diatomic radicals with valence s or d orbitals by systematically varying active space size in the restricted active space self-consistent field formalism with restricted active space state interaction for spin–orbit coupling. The DKH relativistic treatment of the hyperfine interaction reduces the Fermi contact contribution to the HFC due to the presence of kinetic factors that regularize the singularity of the Dirac delta function in the nonrelativistic Fermi contact operator. This effect is more prominent for heavier nuclei. As the active space size increases, the relativistic correction of the Fermi contact contribution converges well to the experimental data for light and moderately heavy nuclei. The relativistic correction, however, does not significantly affect the spin-dipole contribution to the hyperfine interaction. In addition to the atomic and molecular systems, the implementation is applied to calculate the relativistic HFC parameters for large trivalent and divalent Tb-based single-molecule magnets (SMMs), such as Tb(III)Pc2 and Tb(II)(CpiPr5)2 without ligand truncation using well-converged basis sets. In particular, for the divalent SMM, which has an unpaired valence 6s/5d hybrid orbital, the relativistic treatment of HFC is crucial for a proper description of the Fermi contact contribution. Even with the relativistic hyperfine Hamiltonian, the divalent SMM is shown to exhibit strong tunability of HFC via an external electric field (i.e., strong hyperfine Stark effect).

Chemistry↗
DOE OSTI2025

Discovery of new isotopes in the fragmentation of 82 Se and insights into their production

The production cross sections of neutron-rich nuclei for elements just above 60 Ca were measured at the Facility for Rare Isotope Beams (FRIB). Four previously unobserved isotopes ( 63 Sc, 65,66 Ti, and 68 V) were produced, separated, and identified for the first time using the Advanced Rare Isotope Separator (ARIS). One event was found to be consistent with 61 Ca . The new isotopes were created through the interaction of an 82 Se beam with a carbon target at an energy of 228 MeV/𝑢 and a primary beam current of 1.07 p⁢µ⁢A. The event-by-event particle identification of the mass number (𝐴), atomic number (𝑍), and ionic charge state (𝑞) for the reaction products was achieved by combining measurements of energy loss, time of flight, magnetic rigidity, and total kinetic energy. This successful search for new isotopes, conducted at the beginning of FRIB's third year of operation, highlights the facility's discovery potential, which will continue to grow as the beam current increases. The production cross sections were analyzed in the context of the newly developed Δ⁢𝐵⁢𝐸 systematics that provides trends in neutron-rich regions and increased sensitivity to nuclear binding, an important nuclear structure observable. Furthermore, the new systematics addresses gaps in 𝑄 𝑔 -based methods commonly used to interpret the production cross sections of light neutron-rich elements, and provides an efficient empirical approach to describe cross-section trends across isotopic chains. Good agreement with the measured data was obtained when the Δ⁢𝐵⁢𝐸 systematics and abrasion-ablation calculations were performed using the Hartree-Fock-Bogoliubov (HFB-22) mass table.

59 ≤ A ≤ 89↗
DOE OSTI2024

Enhancing Lattice Kinetic Schemes for Fluid Dynamics with Lattice-Equivariant Neural Networks

A new class of equivariant neural networks is presented, hereby dubbed lattice-equivariant neural networks (LENNs), designed to satisfy local symmetries of a lattice structure. The approach develops within a recently introduced framework aimed at learning neural network-based surrogate models’ lattice Boltzmann collision operators. Whenever neural networks are employed to model physical systems, respecting symmetries and equivariance properties has been shown to be key for accuracy, numerical stability, and performance. Here, hinging on ideas from group representation theory, trainable layers are defined whose algebraic structure is equivariant with respect to the symmetries of the lattice cell. In this work, the presented method naturally allows for efficient implementations, in terms of both memory usage and computational costs, supporting scalable training/testing for lattices in two spatial dimensions and higher (in which the size of symmetry group grows). The approach is validated and tested considering 2D and 3D flowing dynamics, both in laminar and turbulent regimes. It is compared with group-averaged-based symmetric networks and with plain, nonsymmetric, networks, showing how the presented approach unlocks the (a posteriori) accuracy and training stability of the former models and the train/inference speed of the latter networks. (LENNs are about one order of magnitude faster than group-averaged networks in 3D.) The work in this paper opens toward practical use of machine learning-augmented lattice Boltzmann CFD in real-world simulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Initial Uncertainty Analysis of Carbon Tetrachloride Contamination and Remediation in the Ringold A and Lower Mud Units at the Central Plateau

The long-term effectiveness of groundwater cleanup at the Hanford Site Central Plateau depends on predictive models that can capture key uncertainties in contaminant fate and transport. Carbon tetrachloride (CCl 4 ), a persistent and toxic compound, presents particular challenges due to variability in degradation rates, uncertainty in initial plume distribution, and subsurface heterogeneity. These uncertainties directly influence plume persistence, migration pathways, and remedy performance, and thus must be systematically evaluated to support long-term remediation planning. To address these gaps, a large-scale Monte Carlo analysis was conducted using the Plateau to River (P2R) model framework. The modeling approach parameterized three primary uncertainty factors: (1) degradation rate, (2) initial plume distribution, and (3) hydraulic conductivity. Degradation was represented as a first-order process, with half-lives ranging from 70 to 700 years. Initial plume distributions were created using a geostatistical simulation method (sgsim), which generates many equally plausible versions of how contaminants might be distributed underground. From this, 100 different scenarios were mapped onto the P2R grid. Variability in hydraulic conductivity was represented in a similar way, with 100 scenarios each for the Ringold Lower Mud and Ringold A units (layers 6 and 7), based on fitted exponential variograms and conditioned to well data. In total, more than 1000 realizations were simulated to assess plume behavior under uncertainty. Results demonstrate that degradation kinetics exert the strongest control over plume persistence: Shorter half-lives produced rapid mass reduction, while longer half-lives yielded persistent plumes with limited attenuation. A nonlinear response was observed, with steep mass reductions at half-lives greater than 200 years and near-linear declines beyond this threshold, reflecting interactions between degradation and pumping. The initial plume distribution strongly influenced early transport patterns, with broader sources generating larger plume footprints, although pump-and-treat operations constrained plume migration to managed areas. By comparison, hydraulic conductivity variability in the Ringold units had only a secondary influence, modifying spreading behavior without altering the dominant migration pathways governed by source configuration and hydraulic controls. Overall, the analysis highlights that uncertainty in degradation rate and initial plume configuration are the primary drivers of variability in plume predictions, while conductivity heterogeneity plays a limited role. These findings underscore the need for improved site-specific data on degradation processes and source characterization to enhance the reliability of long-term performance assessments and to better inform remedial decision-making at the Central Plateau.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Supramolecular assembly of polycation/mRNA nanoparticles and in vivo monocyte programming

Size-dependent phagocytosis is a well-characterized phenomenon in monocytes and macrophages. However, this size effect for preferential gene delivery to these important cell targets has not been fully exploited because commonly adopted stabilization methods for electrostatically complexed nucleic acid nanoparticles, such as PEGylation and charge repulsion, typically arrest the vehicle size below 200 nm. Here, we bridge the technical gap in scalable synthesis of larger submicron gene delivery vehicles by electrostatic self-assembly of charged nanoparticles, facilitated by a polymer structurally designed to modulate internanoparticle Coulombic and van der Waals forces. Specifically, our strategy permits controlled assembly of small poly(β-amino ester)/messenger ribonucleic acid (mRNA) nanoparticles into particles with a size that is kinetically tunable between 200 and 1,000 nm with high colloidal stability in physiological media. We found that assembled particles with an average size of 400 nm safely and most efficiently transfect monocytes following intravenous administration and mediate their differentiation into macrophages in the periphery. When a CpG adjuvant is co-loaded into the particles with an antigen mRNA, the monocytes differentiate into inflammatory dendritic cells and prime adaptive anticancer immunity in the tumor-draining lymph node. This platform technology offers a unique ligand-independent, particle-size-mediated strategy for preferential mRNA delivery and enables therapeutic paradigmsviamonocyte programming.

Science & Technology - Other Topics↗
DOE OSTI2024

Unraveling the Pyrolytic Behavior and Kinetics of Single Polymers and Plastic-Rich Municipal Solid Waste Using Thermal Analysis

Pyrolysis is a highly promising thermochemical recycling technology for converting heterogenous plastic waste into sustainable fuels in a single step. Therefore, understanding the pyrolysis mechanism is essential for enabling rational reactor design and enhancing efficient recycling techniques. In this study, the thermal degradation behaviors and corresponding kinetics of pure polymers (PE, PP, and PET) and plastic-rich MSW were examined using simultaneous thermogravimetric analysis (TGA) and differential scanning calorimetry (DSC). Experiments were carried out in the temperature range of 30-800°C with variable heating rates from 5°C/min to 20°C/min in an ultra-high purity Argon atmosphere. Our results indicated that the plastic pyrolysis was an endothermic process, with varying decomposition temperature ranges depending on their structure and composition. Various iso-conversional model-free methods (Friedman, Flynn-Wall-Ozawa, Starink, and Kissinger-Akahira-Sunose) were utilized to determine the apparent activation energy of the plastic degradation, which increased in the following order: PET (214 kJ/mol), PP (218 kJ/mol), PE (245 kJ/mol), and MSW (249 kJ/mol). Finally, Criado’s master plots were employed to identify the best-fitting reaction model and the pre-exponential factor was subsequently determined.

Bashir, Muhammad Aamir↗
DOE OSTI2025

Turbulence in Core-collapse Supernovae

It is understood in a general sense that turbulent fluid motion below the shock front in a core-collapse supernova stiffens the effective equation of state of the fluid and aids in the revival of the explosion. However, when one wishes to be precise and quantify the amount of turbulence in a supernova simulation, one immediately encounters the problem that turbulence is difficult to define and measure. Using the 3D magnetohydrodynamic code ELEPHANT, we study how different definitions of turbulence change one’s conclusions about the amount of turbulence in a supernova and the extent to which it helps the explosion. We find that, while all the definitions of turbulence we use lead to a qualitatively similar growth pattern over time of the turbulent kinetic energy in the gain region, the total amount of turbulent kinetic energy, and especially the ratios of turbulent to total kinetic energy, distinguish them. Some of the definitions appear to indicate turbulence is a necessary contributor to the explosion, and others indicate it is not. The different definitions also produce turbulence maps with different correlations with maps of the enstrophy, a quantity widely regarded as also indicating the presence of turbulence. We also compute the turbulent adiabatic index and observe that, in regions of low enstrophy, this quantity is sensitive to the definition used. As a consequence, the effective adiabatic index depends upon the method used to measure the turbulence, and thus it alters one’s conclusions regarding the impact of turbulence within the supernova.

computational astronomy↗
DOE OSTI2025

Modeling MTS pyrolysis and SiC deposition kinetics using principal component analysis and neural networks

Accurate chemical kinetics modeling is crucial for improving the efficiency of chemical processing and synthesis of ceramic matrix composites. Detailed kinetic models are computationally expensive due to the large number of transported chemical species, while the simplified physics-based models, such as single-step global mechanisms, are efficient but often overlook key chemical intermediates and pathways. Recent deep learning approaches promise accurate and cost-effective models. Yet, they require additional closures for the transported nonlinear latent variables, complicating integration with existing solvers. In this work, we develop a hybrid linear—nonlinear reduced model for silicon carbide deposition from methyltrichlorosilane precursor by combining principal component analysis (PCA) and autoencoder (AE) neural network (NN) approaches. PCA is used to identify a smaller set of linear transport variables, enabling direct reuse of conventional transport solvers. NNs then reconstruct the full chemical state from these reduced variables. We demonstrate the method on a chemical vapor deposition reactor—comprising a gas-phase pyrolysis plug flow reactor and a heterogeneous surface reactor—over a wide range of temperatures, pressures, and residence times. Our PCA–AE model achieves high accuracy with only five transported scalars, achieving an eightfold cost reduction compared to detailed mechanisms, in both a priori (using data from the test set only) and a posteriori (coupled with a differential equation solver). In conclusion, notable errors arise primarily near training domain boundaries and for long residence times, indicating the need for domain shift indicators and better long-horizon predictions in future reduced chemistry model development.

autoencoder neural networks↗
DOE OSTI2026

Biopolymer-Templated Titania Film Formation for Nanostructured Coatings Revealed by Machine Learning-Supported Time-Resolved Analysis

This study presents a machine learning approach to derive the film formation of biopolymer-templated titania nanostructures during spray deposition, in combination with in situ grazing-incidence small-angle X-ray scattering (GISAXS). A neural network trained on synthetic GISAXS data directly predicts domain-size distributions from experimental two-dimensional scattering patterns, capturing the full kinetics of nanostructure evolution with high temporal resolution. The predictions reveal hierarchical size distributions and periodic growth features, consistent with layer-by-layer spray deposition and validated by complementary scanning electron microscopy (SEM) imaging. Quantitative comparison with conventional parametric GISAXS fits shows good qualitative agreement, with systematic differences explained by domain-shape assumptions and resolved by applying a geometric scaling factor. Simulated SEM-like surfaces derived from neural network outputs reproduce the porous, foam-like nanoscale morphology observed experimentally, reinforcing the method’s credibility. This integrated approach enables real-time, nondestructive, statistically averaged monitoring of bulk nanostructure development in functional coatings, offering a scalable methodology to accelerate the characterization and process control of sustainably manufactured nanostructured titania films for energy-related applications such as photocatalysis and photovoltaics.

Heger, JulianEliah↗
DOE OSTI2026

Aluminum Based Solvent-Free Organic–Inorganic Hybrid Materials

In emerging materials, molecular hybrids are especially promising, as they have molecular level mixing of the organic and inorganic components, producing homogeneous materials without interfaces that can deteriorate properties. However, the current methods of manufacturing molecular hybrids are based on solution processing, which is impractical for bulk materials such as may be used for optically clear radiation and electromagnetic shielding components or photonics. Here we examine molecular hybrids composed of aluminum isopropoxide (AIP) and epoxy resins aiming to understand the molecular scale chemistry and manufacturability of these hybrid materials. DSCmonitored cure revealed the ideal cure temperature for these materials is 160−170 °C and demonstrated that an AIP concentration of 16.7 wt % maximizes the extent of reaction. Kinetic analysis of the curing reaction showed the Sestak−Berggren autocatalytic model is effective at temperatures over 140 °C but the reaction has diffusion limitations at a temperature of 120 °C. Mechanical testing with custom resin molds revealed a decrease in properties of the bulk samples with increasing AIP content due to an increase in defects but further testing with nanoindentation demonstrated comparable or improved mechanical properties of AIP-epoxy hybrids compared to epoxy resin with a standard hardener. Ultimately, this work lays the foundation for hardener-free epoxy-aluminum inorganic/organic hybrids and presents opportunities to expand on properties for specific applications such as thermal conductivity, optical clarity, and dielectric constant.

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

Synergistic effects of Pd single atoms and nanoclusters boosting SnO 2 gas sensing performance

Tin(IV) oxide-supported Pd is a promising heterogenous catalyst for CO oxidation relevant for environmental cleanup reactions. In this study, an atomically dispersed Pd catalyst on SnO 2 (ADC Pd/SnO 2 ) hybrid material is successfully synthesized via a straightforward wet chemistry method and is found to exhibit superior performance toward CO sensing. Ex situ EXAFS analysis confirms the formation of single Pd atoms and small Pd nanoclusters stabilized on the SnO 2 (110) surface. Further, the material exhibits high efficiency in generating adsorbed O 2 - as well as high activity in catalyzing CO oxidation at low temperatures, resulting in exceptional sensitivity and selectivity toward CO in comparison to pure SnO 2 and Pd nanoparticles loaded on SnO 2 respectively. In situ FTIR measurements unravel CO adsorption kinetics on ADC Pd/SnO 2 under reaction conditions, and a possible sensing mechanism is put forth in which CO is transformed into CO 2 by reaction with active oxygen species; and concurrently, carbon-related species (bicarbonates and carbonates) are formed and decomposed into CO 2 .

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

Energy conservation in real-time nuclear–electronic orbital Ehrenfest dynamics

Real-time nuclear–electronic orbital Ehrenfest (RT-NEO-Ehrenfest) dynamics methods provide a first-principles approach for describing nonadiabatic molecular processes with nuclear quantum effects. For an efficient description of proton transfer within RT-NEO-Ehrenfest dynamics, the basis function center associated with the quantum proton can be allowed to move classically. Here, this traveling proton basis (TPB) approach effectively captures proton quantum dynamics, although its energy conservation behavior is not yet fully satisfactory. Two recently proposed TPB approaches, in principle, conserve the extended energy, which includes both the system energy and the kinetic energy associated with the proton basis function center. Herein, a thermostatted TPB approach is proposed to improve the conservation of the system energy, excluding the kinetic energy associated with the proton basis function center. In this approach, the quantum proton dynamics are modulated by dynamically rescaling the proton momentum operator to maintain the system energy conservation. With the excited-state intramolecular proton transfer of o-hydroxybenzaldehyde as an example, this approach is shown to significantly improve the system energy conservation while preserving the accuracy of the quantum proton dynamics as achieved in the original TPB approach.

Ehrenfest dynamics↗
DOE OSTI2026

Non-local transport in radiation-hydrodynamics codes for ICF by efficient coupling to an external Vlasov–Fokker–Planck code

Accurately incorporating non-local transport into radiation-hydrodynamics codes, and indeed any fluid system, has long been elusive. To date, a simplified and accurate theory that can be easily integrated has not been available. This limitation affects modeling in inertial confinement fusion (ICF) and magnetic confinement fusion systems, among others, where non-local transport is well-known to be present. Here, we present a coupling methodology between a full Vlasov–Fokker–Planck (VFP) electron kinetic code and radiation-hydrodynamics (rad-hydro) codes. The VFP code is used to adjust native electron transport in the rad-hydro code, thus enabling improved transport without the need to integrate a full electron VFP solver into the rad-hydro code. This approach necessitates only occasional invocation of the VFP code, reducing computational intensity compared to following the dynamic evolution entirely with the VFP code on fluid time scales. We illustrate that the methodology is more accurate than other simplified methods in thermal decay systems relevant to ICF and can replicate standard theoretical results with high accuracy.

Electronic transport↗
DOE OSTI2025

Estimation of extreme temperatures in direct solar methane pyrolysis within a porous medium

Porous media have wide application in renewable energy conversion processes, such as solar-thermal fuels production and decarbonization. Heat transport mechanisms within porous media can be highly complex, particularly under extreme conditions encountered in concentrated solar thermal reactors in which direct measurement of temperature is challenging. Here, we implement and report an inverse heat conduction model to estimate the temperature distribution throughout a porous substrate domain in a direct solar methane pyrolysis process. By solving a two-dimensional heat transfer problem and applying an inverse optimization algorithm, we estimate the quasi-steady state spatial temperature distribution in a fibrous porous carbon substrate. The results are validated indirectly by experimentally measured graphite deposition and a simplified reaction kinetic model.

finite difference method↗