SEARCH · Engineering Papers

Results for “convergence rate”

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.

SourceResearch connectionFrom yearThrough year

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

At least 199 records · Page 11

All sources NASA NTRS DOE OSTI

DOE OSTI2026

A lumped particle direct simulation Monte-Carlo method combined with the collisional-radiative model for simulations of non-equilibrium laser-induced plasma plumes

Collisional plasma plumes induced by laser irradiation of material targets exhibit large variations in local density as well as ionization and excitation states, making purely hydrodynamic or kinetic simulations inaccurate or infeasible. To address this challenge and capture non-equilibrium effects in laser-induced plasma plumes at arbitrary degrees of ionization, we develop a hybrid computational approach that combines the kinetic direct simulation Monte Carlo (DSMC) method with a collisional-radiative model (CRM). This ℓDSMC-CRM approach utilizes a lumped particle method to represent minor fractions of excited ions in particle-based simulations and a special coarse-graining technique for atomic spectra and photoionization rates, ensuring numerical convergence at reduced computational cost. The hybrid approach is applied to simulate spatially homogeneous relaxation as well as one- and two-dimensional expansions of plasma plumes induced by irradiation of a copper target by a nanosecond laser pulse in a vacuum or background gas. The comparison with an equilibrium model, where local Saha-Boltzmann equilibrium is enforced, shows that the non-equilibrium effects play a dominant role. The equilibrium model can fail to predict the flow structure and strongly underestimate the degree of absorption of laser radiation by the plume. The ℓDSMC-CRM approach is validated against experimental data demonstrating reasonable agreement with the experimental electron density and temperature, while the equilibrium model is found to dramatically underestimate electron density and temperature. The flexibility of the ℓDSMC-CRM approach allows for its seamless integration into existing DSMC frameworks, making it a valuable tool for high-fidelity plasma modeling in laser-material interactions, laser-based manufacturing, and beyond.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Self-regulating behavior of hybrid membrane systems as demonstrated in an element-scale forward osmosis-reverse osmosis hybrid system

Hybrid membrane systems can be difficult to design due to the requisite flow rate matching between up- and downstream unit operations. In this work, we use a forward osmosis-reverse osmosis (FO-RO) hybrid system to demonstrate how some membrane systems can exhibit self-regulating behavior due to osmotic coupling. This can reduce the need for complex control systems for flow balancing. We show this behavior using a module-scale test bed that can mimic the behavior of larger scale operations. The system shows permeate flow rate near-convergence between the FO and RO modules after startup or when perturbed by a change in RO module pressure. The behavior of this hybrid system demonstrates that some membrane operations can exploit osmotic interdependence, rather than expensive control systems, to achieve steady state operation.

Debottlenecking↗

DOE OSTI2024

Direct, simple and efficient computation of all components of the virtual-casing magnetic field in axisymmetric geometries with Kapur–Rokhlin quadrature

In a recent publication (Toler et al., J. Plasma Phys., vol. 89, issue 2, 2023, p. 905890210), we demonstrated that for axisymmetric geometries, the Kapur–Rokhlin quadrature rule provided an efficient and high-order accurate method for computing the normal component, on the plasma surface, of the magnetic field due to the toroidal current flowing in the plasma, via the virtual-casing principle. The calculation was indirect, as it required the prior computation of the magnetic vector potential from the virtual-casing principle, followed by the computation of its tangential derivative by Fourier differentiation, to obtain the normal component of the magnetic field. Our approach did not provide the other components of the virtual-casing magnetic field. In this letter, we show that a more direct and more general approach is available for the computation of the virtual-casing magnetic field. The Kapur–Rokhlin quadrature rule accurately calculates the principal value integrals in the expression for all the components of the magnetic field on the plasma boundary, and the numerical error converges at a rate nearly as high as the indirect method we presented previously.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Long term performance of vacuum insulation panels integrated into building components

Vacuum insulation panels provide much greater resistance to heat transfer than conventional insulation materials, making them ideal for building applications, especially where space is a constraint. Of particular interest is energy retrofits where high insulation values are required with very little or no modifications to existing building envelope elements such as fenestrations and service penetrations. To facilitate the use of vacuum insulation panels (VIPs) in residential and commercial construction, this study investigates the performance of VIPs integrated into polyisocyanurate foam board and vinyl siding after several years of service. Vacuum insulation panels integrated with polyisocyanurate boards were studied in mixed humid and cold climates. The integration of VIPs with vinyl siding was investigated in a mixed humid climate. Periods of service were between two and five years. In all cases, the integrated building components outperformed conventional materials by almost 40 percent during heating and cooling periods. However, the percent reduction in heat gains and losses of the integrated building components decreased with time at a rate between approximately 4% and 10% per year during the service period. In some cases, the rate seems to converge to a steady -state value. Measurements are ongoing to determine if a steady-state value is attained or the panel's insulation value continues to decrease with time. This study is critical in addressing the service life and durability of the building envelope with VIP integrated building components. The effect of these changes will be reviewed in the context of an exterior retrofit.

Aldykiewicz Jr, Antonio↗

DOE OSTI2020

Simulating combustion of a seven-component surrogate for a gasoline/ethanol blend including soot formation and comparison with experiments

Combustion of a seven-component surrogate for a research grade 87 octane gasoline mixed with 10% ethanol is investigated experimentally and numerically from the perspective of an isolated droplet burning under conditions that promote one-dimensional gas transport. The numerical analysis included a kinetic mechanism comprised of 398 species and 24,814 reactions and a soot model that accounted for nucleation, surface growth, coalescence/aggregation of soot particles, and luminous flame radiation. Measurements of droplet and flame diameters were made for an initial droplet diameter (D o ) of approximately 0.63 mm. The simulations agreed well with the measurements including the location of the soot shell. Preferential vaporization was revealed by simulations of the liquid concentrations in the droplet. Predicted peak soot volume fractions coincided with temperatures between 1300 K and 1400 K as a soot inception temperature. Simulations were also carried out for D o between 0.25 mm and 5 mm to explore the effect of radiation and D o on burning. Below 0.25 mm radiation was negligible and burning rates and flame temperatures converged to a single value. Increasing D o up to 1.8 mm lowered the burning rate with luminous radiation having a strong effect. When radiation was entirely removed from the model the burning rate was nearly constant. Above D o = 2 mm droplets extinguished almost immediately after ignition. The flame temperature decreased with increasing D o while it increased when radiation was omitted. As a result, the simulations show that soot precursors including polyaromatic hydrocarbons were concentrated around the soot shell.

42 ENGINEERING↗

DOE OSTI2026

Convergence Criteria for Multiphysics Simulations

The behavior of engineered systems is often influenced by multiple physical phenomena, such as mechanical deformation, heat transfer, and chemical species transport and reactions. There are often strong interactions between these phenomena, and there is increasing interest in applying coupled-physics models to improve understanding of physical behavior under complex environmental conditions. Multiple simulation frameworks that facilitate coupled-physics simulations are in widespread use, and these employ a variety of techniques to account for interactions between those physics. Many frameworks solve the physics models independently and transfer results between them. Alternatively, a single monolithic system of equations for every physics model can be formed and solved. Each of these approaches has its benefits and drawbacks, and the optimal approach varies depending on the nature of the problem. The open-source MOOSE framework was developed targeting solution of large-scale multiphysics problems. Although it provides options for all these coupling approaches, its standard approach for multiphysics solutions is to form and solve a single monolithic system of equations containing the unknowns for all physics models. MOOSE provides a streamlined approach for users to define the solution variables, the terms in the partial differential equations pertaining to each variable, and interactions between solution variables. One aspect of the monolithic solution approach that can be problematic, however, is defining appropriate convergence criteria for the nonlinear system. A standard approach is to determine convergence is to simply take a norm of the residual vector corresponding to the full vector of unknowns. However, if the residual vector contains variables for multiple physics models, the magnitudes of those variables can differ significantly, and the variables can converge at significantly different rates from each other. It is important to ensure that the variables for each of the physics are converged, and also ensure that the convergence criteria are not excessively stringent in cases when there is little change in the solution. This talk presents representative multiphysics problems to highlight these issues, and shows strategies for convergence criteria in MOOSE that are robust for multiphysics models under a variety of conditions.

97 - MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Linear convergence of accelerated conditional gradient algorithms in spaces of measures

A class of generalized conditional gradient algorithms for the solution of optimization problem in spaces of Radon measures is presented. The method iteratively inserts additional Dirac-delta functions and optimizes the corresponding coefficients. Under general assumptions, a sub-linear [see formula in PDF] rate in the objective functional is obtained, which is sharp in most cases. To improve efficiency, one can fully resolve the finite-dimensional subproblems occurring in each iteration of the method. We provide an analysis for the resulting procedure: under a structural assumption on the optimal solution, a linear [see formula in PDF] convergence rate is obtained locally.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

High-throughput quantification of quasistatic, dynamic and spall strength of materials across 10 orders of strain rates

Abstract The response of metals and their microstructures under extreme dynamic conditions can be markedly different from that under quasistatic conditions. Traditionally, high strain rates and shock stresses are achieved using cumbersome and expensive methods such as the Kolsky bar or large spall experiments. These methods are low throughput and do not facilitate high-fidelity microstructure–property linkages. In this work, we combine two powerful small-scale testing methods, custom nanoindentation, and laser-driven microflyer (LDMF) shock, to measure the dynamic and spall strength of metals. The nanoindentation system is configured to test samples from quasistatic to dynamic strain-rate regimes. The LDMF shock system can test samples through impact loading, triggering spall failure. The model material used for testing is magnesium alloys, which are lightweight, possess high-specific strengths, and have historically been challenging to design and strengthen due to their mechanical anisotropy. We adopt two distinct microstructures, solutionized (no precipitates) and peak-aged (with precipitates) to demonstrate interesting upticks in strain-rate sensitivity and evolution of dynamic strength. At high shock-loading rates, we unravel an interesting paradigm where the spall strength vs. strain rate of these materials converges, but the failure mechanisms are markedly different. Peak aging, considered to be a standard method to strengthen metallic alloys, causes catastrophic failure, faring much worse than solutionized alloys. Our high-throughput testing framework not only quantifies strength but also teases out unexplored failure mechanisms at extreme strain rates, providing valuable insights for the rapid design and improvement of materials for extreme environments.

Eswarappa Prameela, Suhas (ORCID:0000000334530184)↗

DOE OSTI2020

Ab initio electron-two-phonon scattering in GaAs from next-to-leading order perturbation theory

Electron-phonon (e–ph) interactions are usually treated in the lowest order of perturbation theory. Here we derive next-to-leading order e–ph interactions, and compute from first principles the associated electron-two-phonon (2ph) scattering rates. The derivations involve Matsubara sums of two-loop Feynman diagrams, and the numerical calculations are challenging as they involve Brillouin zone integrals over two crystal momenta and depend critically on the intermediate state lifetimes. Using Monte Carlo integration together with a self-consistent update of the intermediate state lifetimes, we compute and converge the 2ph scattering rates, and analyze their energy and temperature dependence. We apply our method to GaAs, a weakly polar semiconductor with dominant optical-mode long-range e–ph interactions. We find that the 2ph scattering rates are as large as nearly half the value of the one-phonon rates, and that including the 2ph processes is necessary to accurately predict the electron mobility in GaAs from first principles.

36 MATERIALS SCIENCE↗

DOE OSTI2024

A Particle-in-Cell Method for Plasmas with a Generalized Momentum Formulation, Part II: Enforcing the Lorenz Gauge Condition

In a previous paper Christlieb et al. (A particle-in-cell method for plasmas with a generalized momentum formulation, part I: Model formulation, 2024), we developed a new particle-in-cell (PIC) method for the relativistic Vlasov–Maxwell system in which the electromagnetic fields and the equations of motion for the particles were cast in terms of scalar and vector potentials through a Hamiltonian formulation. This new method evolved the potentials under the Lorenz gauge using integral equation methods. New methods to construct spatial derivatives of the potentials that converge at the same rates as the fields were also presented. The new particle method was compared against standard explicit discretizations, including the well-known FDTD-PIC method, for a range of applications involving sheaths and particle beams. Here, this paper extends this new class of methods by focusing on the enforcement the Lorenz gauge condition in both exact and approximate forms using co-located meshes. A time-consistency property of the proposed field solver for the vector potential form of Maxwell’s equations is established, which is shown to preserve the equivalence between the semi-discrete Lorenz gauge condition and the analogous semi-discrete continuity equation. Using this property, we present three methods to enforce a semi-discrete gauge condition. The first method introduces an update for the continuity equation that is consistent with the discretization of the Lorenz gauge condition. Both the finite difference and spectral implementations satisfy this discrete gauge condition to machine precision. The second approach we propose enforces a semi-discrete continuity equation using the boundary integral solution to the field equations. The potential benefit of this approach is that it eliminates spatial derivatives that appear on the particle data, namely the current density, which is often calculated by linear combinations of low-order spline basis functions. This method is ideally suited to boundary integral equation methods that invert multi-dimensional operators without dimensional splitting techniques and will be the subject of future work. The third approach introduces a gauge correcting method that makes direct use of the gauge condition to modify the scalar potential and uses local maps for both the charge and current densities. This results in a gauge error, as the maps do not enforce the continuity equation. The vector potential coming from the current density is taken to be exact, and using the Lorenz gauge, we compute a correction to the scalar potential that makes the two potentials satisfy the gauge condition. This method also enforces the gauge condition to machine precision. We demonstrate two of the proposed methods in the context of periodic domains. Problems defined on bounded domains, including those with complex geometric features remain an ongoing effort. However, this work shows that it is possible to design computationally efficient methods that can effectively enforce the Lorenz gauge condition in a non-staggered PIC formulation.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

A Particle-in-Cell Method for Plasmas with a Generalized Momentum Formulation, Part I: Model Formulation

Here, this paper formulates a new particle-in-cell method for the Vlasov–Maxwell system. Under the Lorenz gauge condition, Maxwell’s equations for the electromagnetic fields can be written as a collection of scalar and vector wave equations. The use of potentials for the fields motivates the adoption of a Hamiltonian formulation for particles that employs the generalized (conjugate) momentum. A notable advantage offered by the Hamiltonian formulation is the elimination of time derivatives in the Lorenz gauge formulation that are required by the standard Newton–Lorentz treatment of the particles. This allows the fields to retain the full time-accuracy guaranteed by the field solver. The resulting updates for particles require only knowledge of the fields and their spatial derivatives. An analytical method for constructing these spatial derivatives is presented that exploits the underlying integral solution used in the field solver for the wave equations. Moreover, these derivatives are demonstrated to converge at the same rate as the fields in both time and space. The Method of Lines Transpose field solver we consider in this work is globally first-order accurate in time and high-order accurate in space (e.g., fourth- and fifth-order) and belongs to a larger class of methods which are unconditionally stable, can address geometry, and leverage $\mathcal {O}(N)$ fast summation methods for efficiency. We demonstrate the method on several well-established benchmark problems on bounded domains, including a plasma sheath as well as a relativistic particle beam. The efficacy of the proposed formulation is established by comparing with a second-order accurate finite-difference time-domain method that employs a leapfrog time advance for particles and a charge conserving map suitable for bounded domains. The new method shows mesh-independent numerical heating properties even in cases where the plasma Debye length is smaller than the grid spacing. This is an important feature of the new method for problems defined on bounded domains, because it permits the use of coarser grids in space in the representation of the fields. Such a capability has significant implications for the simulation of plasmas in bounded domains with complex geometry, where the ratio between the largest and smallest cells can vary significantly. The use of high-order spatial approximations in the new method also means that fewer grid points are required in order to achieve a fixed accuracy. Our results also suggest that the new method can be used with fewer simulation particles per cell compared to the benchmark explicit method, which permits further computational savings.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Development and Implementation of a CTF Code Verification Suite.

CTF is a thermal hydraulic subchannel code developed to predict light water reactor (LWR) core behavior. It is a version of Coolant Boiling in Rod Arrays (COBRA) developed by Oak Ridge National Laboratory (ORNL) and North Carolina State University (NCSU) and used in the Consortium for the Advanced Simulation of LWRs (CASL). Here, the existing CTF code verification matrix is expanded, which ensures that the code is a faithful representation of the underlying mathematical model. The suite of code verification tests are mapped to the underlying conservation equations of CTF and significant gaps are addressed. As such, five new problems are incorporated: isokinetic advection, conduction, pressure drop, convection, and pipe boiling. Convergence behavior and numerical errors are quantified for each of the tests and all tests converge at the correct rate to their corresponding analytic solution. A new verification utility that generalizes the code verification process is used to incorporate these problems into the CTF automated test suite.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2020

Improvements to CTF Code Verification and Unit Testing (FY2020)

In 2010, the U.S. Department of Energy created its first Energy Innovation Hub, which focuses on improving Light Water Reactors (LWRs) through Modeling and Simulation. This hub, named the Consortium for the Advanced Simulation of LWRs (CASL), attempts to characterize and understand LWR behavior under normal operating conditions and use any gained insights to improve their efficiency. In collaboration with North Carolina State University (NCSU), CASL has worked extensively on the thermal-hydraulic subchannel code Coolant Boiling in Rod Arrays—Three Field (COBRA-TF). The NCSU/CASL version of COBRA-TF has been rebranded as CTF. This document focuses on code verification test problems that ensure CTF converges to the correct answer for the intended application. The suite of code verification tests are mapped to the underlying conservation equations of CTF, and significant gaps are addressed. Convergence behavior and numerical errors are quantified for each of the tests. Tests that converge at the correct rate to the corresponding analytic solution are incorporated into the CTF automated regression suite. A new verification utility is created for this purpose, which enables code verification by generalizing the process. For problems that do not behave correctly, the results are reported but the problem is not included in the regression suite. In addition to verification studies, this document also quantifies the existing tests of constitutive models. A few existing gaps are addressed by adding new unit tests.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2020

FY20 Improvements to CTF Code Verification and Unit Testing

In 2010, the U.S. Department of Energy created its first Energy Innovation Hub, which focuses on improving Light Water Reactors (LWRs) through Modeling and Simulation. This hub, named the Consortium for the Advanced Simulation of LWRs (CASL), attempts to characterize and understand LWR behavior under normal operating conditions and use any gained insights to improve their efficiency. In collaboration with North Carolina State University (NCSU), CASL has worked extensively on the thermal-hydraulic subchannel code Coolant Boiling in Rod Arrays–Three Field (COBRA-TF). The NCSU/CASL version of COBRA-TF has been rebranded as CTF. This document focuses on code verification test problems that ensure CTF converges to the correct answer for the intended application. The suite of code verification tests are mapped to the underlying conservation equations of CTF, and significant gaps are addressed. Convergence behavior and numerical errors are quantified for each of the tests. Tests that converge at the correct rate to the corresponding analytic solution are incorporated into the CTF automated regression suite. A new verification utility is created for this purpose, which enables code verification by generalizing the process. For problems that do not behave correctly, the results are reported but the problem is not included in the regression suite. In addition to verification studies, this document also quantifies the existing tests of constitutive models. A few existing gaps are addressed by adding new unit tests.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2023

Smoothing tool design and performance during subaperture glass polishing

During subaperture tool grinding and polishing, overlaps of the tool influence function can result in undesirable mid-spatial frequency (MSF) errors in the form of surface ripples, which are often corrected using a smoothing polishing step. Here, in this study, flat multi-layer smoothing polishing tools are designed and tested to simultaneously (1) reduce or remove MSF errors, (2) minimize surface figure degradation, and (3) maximize the material removal rate. A time-dependent convergence model in which spatial material removal varies with a workpiece-tool height mismatch, combined with a finite element mechanical analysis to determine the interface contact pressure distribution, was developed to evaluate various smoothing tool designs as a function of tool material properties, thicknesses, pad textures, and displacements. An improvement in smoothing tool performance is achieved when the gap pressure constant, $\overline{h}$ (which describes the inverse rate at which the pressure drops with a workpiece-tool height mismatch), is minimized for smaller spatial scale length surface features (namely, MSF errors) and maximized for large spatial scale length features (i.e., surface figure). Five specific smoothing tool designs were experimentally evaluated. A two-layer smoothing tool using a thin, grooved IC1000 polyurethane pad (with a high elastic modulus, E pad =360 MPa ), thicker blue foam (with an intermediate modulus, E foam =5.3 MPa ) underlayer, and an optimized displacement (d t =1 mm ) provided the best overall performance (namely, high MSF error convergence, minimal surface figure degradation, and high material removal rate).

47 OTHER INSTRUMENTATION↗

DOE OSTI2022

Black Hole Flares: Ejection of Accreted Magnetic Flux through 3D Plasmoid-mediated Reconnection

Magnetic reconnection can power bright, rapid flares originating from the inner magnetosphere of accreting black holes. We conduct extremely high-resolution (5376 × 2304 × 2304 cells) general-relativistic magnetohydrodynamics simulations, capturing plasmoid-mediated reconnection in a 3D magnetically arrested disk for the first time. We show that an equatorial, plasmoid-unstable current sheet forms in a transient, nonaxisymmetric, low-density magnetosphere within the inner few Schwarzschild radii. Magnetic flux bundles escape from the event horizon through reconnection at the universal plasmoid-mediated rate in this current sheet. The reconnection feeds on the highly magnetized plasma in the jets and heats the plasma that ends up trapped in flux bundles to temperatures proportional to the jet's magnetization. The escaped flux bundles can complete a full orbit as low-density hot spots, consistent with Sgr A* observations by the GRAVITY interferometer. Reconnection near the horizon produces sufficiently energetic plasma to explain flares from accreting black holes, such as the TeV emission observed from M87. The drop in the mass accretion rate during the flare and the resulting low-density magnetosphere make it easier for very-high-energy photons produced by reconnection-accelerated particles to escape. The extreme-resolution results in a converged plasmoid-mediated reconnection rate that directly determines the timescales and properties of the flare.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Degrees of Rate Control in Interconnected Reaction Networks

Overall reactions in interconnected networks exhibit net, forward, and reverse rates that are governed by both constitutive elementary steps in the pathway of interest and branching elementary steps that lead to alternative products. Accordingly, steps in branching pathways exhibit negative net, forward, and reverse degrees of rate control, as they reduce reaction flux to the desired product. We here contextualize the forward and reverse degrees of rate control in terms of kinetic resistances (inverse of rates) and leverage the additive nature of kinetic resistance to decouple kinetic driving forces contributed by constitutive elementary steps and branching points (nodal species) in interconnected networks. Regardless of the network connectivity, forward and reverse degrees of rate control are shown to converge at equilibrium. Away from equilibrium, we identify two critical features of interconnected networks: stoichiometric regularity─condition where all stoichiometric numbers are unity─and pathway symmetry around nodal species─condition where branching pathways share the same rate constants, stoichiometry, and species concentrations/activities─that result in (i) equal forward, reverse, and consequently net degrees of rate control and (ii) forward and reverse degrees of rate control that exhibit constant offsets, respectively, across all extents of reaction. Furthermore, our discourse further provides a mathematical description for the influence of stoichiometric irregularity and pathway asymmetry on forward and reverse degrees of rate control. Altogether, the presented work details the effects of network (inter)connectivity and stoichiometry on reaction kinetics and, in doing so, establishes general protocols for capturing these effects as additive terms in the formulation of forward and reverse degrees of rate control.

10 SYNTHETIC FUELS↗

DOE OSTI2026

Quantifying the Impact of Metal Population Distribution in MFI-Supported Mo Catalysts for Methane Dehydroaromatization

Precise evaluation of intrinsic kinetic behavior in Mo/MFI catalysts for methane dehydroaromatization (MDA) is confounded by variations in Mo dispersion and speciation, which are influenced by metal loading and zeolite acidity. This work advances the utility of H 2 -temperature programmed reduction (H 2 -TPR) for characterizing Mo/MFI catalysts by enabling quantitative comparison of MoO x populations distinguished by reduction behavior and linked to initial catalytic performance. A comprehensive H 2 reduction pathway is established through systematic H 2 - TPR studies varying catalyst composition (1–10 wt % Mo, Si/Al = 15, 40, ∞), supplemented by UV-Raman spectroscopy, X-ray powder diffraction, N 2 physisorption, NH 3 -TPD, and advanced spectroscopic analysis (in situ XAS with principal component analysis/multivariate curve resolution—alternating least squares). Two low-temperature H 2 -TPR regions capture distinct Mo populations undergoing initial Mo(VI)→Mo(IV) reduction: a lower-temperature population (Mo-RI) associated primarily with highly dispersed, anchored MoO x species expected to predominantly reside within MFI channels, and a higher-temperature population (Mo-RII) corresponding to a broader set of MoO x species that becomes increasingly bulk-like/extrazeolitic at higher Mo loading. Quantification of these populations provides practical, kinetically relevant descriptors for comparing initial MDA rates across catalysts with varying Mo loading, Si/Al ratio, and MoO x heterogeneity. Application of lower-temperature Mo-RI estimates to kinetic measurements reveals a minimum threshold of ∼0.12 × 10 –3 mol Mo-RI/g cat , above which initial forward benzene rates normalized to this population converge despite differences in metal loading and zeolite Brønsted acidity. This threshold coincides with a transition toward a common C 2 -mediated benzene-forming regime, as indicated by approach-to-equilibrium analysis of methane-to-ethane, ethane-to-ethylene, and ethylene-to-benzene reaction steps. Above this threshold, normalized initial benzene rates are nearly invariant with increasing Mo-RII/Mo-RI population ratio, indicating that excess Mo-RII populations, including bulk-like/extrazeolitic MoO x domains present at higher loading, do not measurably suppress benzene formation associated with anchored and mostly channel-confined Mo population under the initial-rate conditions examined.

Mo/MFI↗