Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Closure”

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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Uncertainty quantification for Multiphase-CFD simulations of bubbly flows: a machine learning-based Bayesian approach supported by high-resolution experiments

In this paper, we developed a machine learning-based Bayesian approach to inversely quantify and reduce the uncertainties of multiphase computational fluid dynamics (MCFD) simulations for bubbly flows. The proposed approach is supported by high-resolution two-phase flow measurements, including those by double-sensor conductivity probes, high-speed imaging, and particle image velocimetry. Local distributions of key physical quantities of interest (QoIs), including the void fraction and phasic velocities, are obtained to support the Bayesian inference. In the process, the epistemic uncertainties of the closure relations are inversely quantified while the aleatory uncertainties from stochastic fluctuations of the system are evaluated based on experimental uncertainty analysis. The combined uncertainties are then propagated through the MCFD solver to obtain uncertainties of the QoIs, based on which probability-boxes are constructed for validation. The proposed approach relies on three machine learning methods: feedforward neural networks and principal component analysis for surrogate modeling, and Gaussian processes for model form uncertainty modeling. The whole process is implemented within the framework of an open-source deep learning library PyTorch with graphics processing unit (GPU) acceleration, thus ensuring the efficiency of the computation. The results demonstrate that with the support of high-resolution data, the uncertainties of MCFD simulations can be significantly reduced. The proposed approach has the potential for other applications that involve numerical models with empirical parameters.

42 ENGINEERING↗
DOE OSTI2026

Computational simulation of asymmetric phase transformation in cracked Li₇La₃Zr₂O₁₂: Variant selection and chemo-mechanical implications

Coupling between microcracks and phase transformation in ion-conducting ceramics can jointly affect mechanical responses and ion transport. In this work we investigate the cubic-to-tetragonal phase transformation in Li₇La₃Zr₂O₁₂ in the presence of a microcrack under hydrostatic loading and quantify its implications for crack-tip stress concentration and Li-ion transport using phase-field and molecular dynamics simulations. The phase transformation exhibits a strong asymmetry between hydrostatic tension and compression. Under tension, the crack edge nucleates one tetragonal variant that amplifies the crack-tip stress intensity and promotes crack opening. Under compression, the crack tip nucleates a different tetragonal variant that enhances the stress-induced crack-closure tendency. Effective Li diffusivity analysis shows faster transport degradation under compression due to accelerated transformation kinetics, exposing a trade-off between mechanical stability and ionic conductivity. These results highlight the intertwined nature of cracking, phase transformation, and ionic transport in ion-conducting oxides and provide mechanistic insights into chemo-mechanical degradation of solid electrolytes.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Developing a model for the impact of non-conformal lithium contact on electro-chemo-mechanics and dendrite growth

Lithium dendrite growth hinders the use of lithium metal anodes in commercial batteries. We present a 3D model to study the mechanical and electrochemical mechanisms that drive microscale plating. With this model, we investigate electrochemical response across a lithium protrusion characteristic of rough anode surfaces, representing the separator as a porous polymer in non-conformal contact with a lithium anode. The impact of pressure on separator morphology and electrochemical response is of particular interest, as external pressure can improve cell performance. We explore the relationships between plating propensity, stack pressure, and material properties. External pressure suppresses lithium plating due to interfacial stress and separator pore closure, leading to inhomogeneous plating rates. For moderate pressures, dendrite growth is completely suppressed, as plating will occur in the electrolyte-filled gaps between anode and separator. In fast-charging conditions and systems with low electrolyte diffusivities, the benefits of pressure are overridden by ion transport limitations.

25 ENERGY STORAGE↗
DOE OSTI2021

Investigation of high-pressure turbulent jets using direct numerical simulation

Direct numerical simulations of free round jets at a Reynolds number (Re D ) of 5000, based on jet diameter (D) and jet-exit bulk velocity (U e ), are performed to study jet turbulence characteristics at supercritical pressures. The jet consists of nitrogen (N 2 ) that is injected into N 2 at the same temperature. To understand turbulent mixing, a passive scalar is transported with the flow at unity Schmidt number. Two sets of inflow conditions that model jets issuing from either a smooth contraction nozzle (laminar inflow) or a long pipe nozzle (turbulent inflow) are considered. By changing one parameter at a time, the simulations examine the jet-flow sensitivity to the thermodynamic condition (characterized in terms of the compressibility factor ($Z$) and the normalized isothermal compressibility), inflow condition and ambient pressure (p ∞ ) spanning perfect- to real-gas conditions. The inflow affects flow statistics in the near field (containing the potential core closure and the transition region) as well as further downstream (containing fully developed flow with self-similar statistics) at both atmospheric and supercritical p ∞ . The sensitivity to inflow is larger in the transition region, where the laminar-inflow jets exhibit dominant coherent structures that produce higher mean strain rates and higher turbulent kinetic energy than in turbulent-inflow jets. Decreasing $Z$ at a fixed supercritical p ∞ enhances pressure and density fluctuations (non-dimensionalized by local mean pressure and density, respectively), but the effect on velocity fluctuations depends also on the local flow dynamics. When $Z$ is reduced, large mean strain rates in the transition region of laminar-inflow jets significantly enhance velocity fluctuations (non-dimensionalized by local mean velocity) and scalar mixing, whereas the effects are minimal in jets from turbulent inflow.

42 ENGINEERING↗
DOE OSTI2023

Transition to turbulence in viscoelastic channel flow of dilute polymer solutions

The transition to turbulence in a plane Poiseuille flow of dilute polymer solutions is studied by direct numerical simulations of a finitely extensible nonlinear elastic fluid with the Peterlin closure. The range of Reynolds number ($Re$)$2000 \le Re \le 5000$is studied but with the same level of elasticity in viscoelastic flows. The evolution of a finite-amplitude perturbation and its effects on the transition dynamics are investigated. A viscoelastic flow begins transition at an earlier time than its Newtonian counterparts, but the transition time appears to be insensitive to polymer concentration in the dilute and semi-dilute regimes studied. Increasing polymer concentration, however, decreases the maximum attainable energy growth during the transition process. The critical or minimum perturbation amplitude required to trigger transition is computed. Interestingly, both Newtonian and viscoelastic flows follow almost the same power-law scaling of$Re^\gamma$with the critical exponent$\gamma \approx -1.25$, which is in close agreement with previous studies. However, a shift downward is observed for viscoelastic flow, suggesting that smaller perturbation amplitudes are required for the transition. A mechanism of the early transition is investigated by the evolution of wall-normal and spanwise velocity fluctuations and flow structure. The early growth of these fluctuations and the formation of quasi-streamwise vortices around low-speed streaks are promoted by polymers, hence causing an early transition. These vortical structures are found to support the critical exponent$\gamma \approx -1.25$. Once the transition process is completed, polymers play a role in dampening the wall-normal and spanwise velocity fluctuations and vortices to attain a drag-reduced state in viscoelastic turbulent flows.

Mechanics↗
DOE OSTI2025

Accurate models of the added mass force of a uniform random distribution of spherical particles or bubbles

The added mass force resulting from the acceleration of a body in a fluid is of fundamental and practical interest in dispersed multiphase flows. Euler–Lagrange (EL) and Euler–Euler (EE) simulations require closure terms for the added mass force in order to accurately couple the conserved variables between phases. Presently, a more thorough understanding of the added mass force in a multi-particle system is developed based on potential flow resulting in a resistance matrix formulation analogous to Stokesian dynamics. This formulation is then used to generate a dataset of added mass resistance matrices for large systems of randomly generated particles. This methodology is used to create a volume fraction corrected binary model for predicting the added mass force in large systems as well as generate statistics of the added mass force in such systems. This work provides clarification to the theory of the added mass force for particle clouds, and modelling options that may be implemented in existing EL and EE codes.

42 ENGINEERING↗
DOE OSTI2024

A streamline coordinate analysis of a turbulent boundary layer subject to pressure gradients and curvature on the windward side of a bump

Direct numerical simulation (DNS) of a turbulent boundary layer over the Gaussian (Boeing) bump is performed. This boundary layer exhibits a series of adverse and favourable pressure gradients and convex and concave curvature effects before separating. These effects on turbulent boundary layers are characterised and compared with a lower-Reynolds-number flow over the same geometry. The momentum budgets are analysed in the streamline-aligned coordinate system upstream of the separation region. These momentum budgets allow the simplification of equations to facilitate an integral analysis. Integral-analysis-based approximations for Reynolds stresses in the inner and outer regions of the boundary layer are also formulated. The shear and wall-normal Reynolds stress profiles normalised by these approximations exhibit a better collapse compared with friction velocity and Zagarola–Smits normalisations in the strong favourable pressure gradient region and in the mild adverse pressure region that precedes it in this flow. Simplification of these Reynolds stress approximations along with results from the DNS are used to obtain semi-empirical approximations that are able to provide stress closure in terms of wall solution fields for the turbulent boundary layer under consideration.

42 ENGINEERING↗
DOE OSTI2024

The vertical-velocity skewness in the inertial sublayer of turbulent wall flows

Empirical evidence is provided that within the inertial sublayer (i.e. logarithmic region) of adiabatic turbulent flows over smooth walls, the skewness of the vertical-velocity component 𝑆 𝑤 displays universal behaviour, being a positive constant and constrained within the range 𝑆 𝑤 ≈ 0.1–0.16, regardless of flow configuration and Reynolds number. A theoretical model is then proposed to explain this behaviour, including the observed range of variations of 𝑆 𝑤 . The proposed model clarifies why 𝑆 𝑤 cannot be predicted from down-gradient closure approximations routinely employed in large-scale meteorological and climate models. The proposed model also offers an alternative and implementable approach for such large-scale models.

turbulent boundary layers↗
DOE OSTI2025

Molecular fluctuations inhibit intermittency in compressible turbulence

In the standard picture of fully developed turbulence, highly intermittent hydrodynamic fields are nonlinearly coupled across scales, where local energy cascades from large scales into dissipative vortices and large density gradients. Microscopically, however, constituent fluid molecules are in constant thermal (Brownian) motion, but the role of molecular fluctuations in large-scale turbulence is largely unknown, and with rare exceptions, it has historically been considered irrelevant at scales larger than the molecular mean free path. Recent theoretical and computational investigations have shown that molecular fluctuations can impact energy cascade at Kolmogorov length scales. Here, we show that molecular fluctuations not only modify energy spectrum at wavelengths larger than the Kolmogorov length in compressible turbulence, but also significantly inhibit spatio-temporal intermittency across the entire dissipation range. Using large-scale direct numerical simulations of computational fluctuating hydrodynamics, we demonstrate that the extreme intermittency characteristic of turbulence models is replaced by nearly Gaussian statistics in the dissipation range. These results demonstrate that the compressible Navier–Stokes equations should be augmented with molecular fluctuations to accurately predict turbulence statistics across the dissipation range. Our findings have significant consequences for turbulence modelling in applications such as astrophysics, reactive flows and hypersonic aerodynamics, where dissipation-range turbulence is approximated by closure models.

compressible turbulence↗
DOE OSTI2024

Self-organization in collisionless, high- β turbulence

The magnetohydrodynamic (MHD) equations, as a collisional fluid model that remains in local thermodynamic equilibrium (LTE), have long been used to describe turbulence in myriad space and astrophysical plasmas. Yet, the vast majority of these plasmas, from the solar wind to the intracluster medium (ICM) of galaxy clusters, are only weakly collisional at best, meaning that significant deviations from LTE are not only possible but common. Recent studies have demonstrated that the kinetic physics inherent to this weakly collisional regime can fundamentally transform the evolution of such plasmas across a wide range of scales. Here, we explore the consequences of pressure anisotropy and Larmor-scale instabilities for collisionless, β $\gg$ 1, turbulence, focusing on the role of a self-organizational effect known as ‘magneto-immutability’. We describe this self-organization analytically through a high-β, reduced ordering of the Chew–Goldberger–Low-MHD (CGL-MHD) equations, finding that it is a robust inertial-range effect that dynamically suppresses magnetic-field-strength fluctuations, anisotropic-pressure stresses and dissipation due to heat fluxes. As a result, the turbulent cascade of Alfvénic fluctuations continues below the putative viscous scale to form a robust, nearly conservative, MHD-like inertial range. These findings are confirmed numerically via Landau-fluid CGL-MHD turbulence simulations that employ a collisional closure to mimic the effects of microinstabilities. We find that microinstabilities occupy a small (~5%) volume-filling fraction of the plasma, even when the pressure anisotropy is driven strongly towards its instability thresholds. We discuss these results in the context of recent predictions for ion-vs-electron heating in low-luminosity accretion flows and observations implying suppressed viscosity in ICM turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Biomimetic Metal-Free Hydride Donor Catalysts for CO 2 Reduction

The catalytic reduction of carbon dioxide to fuels and value-added chemicals is of significance for the development of carbon recycling technologies. One of the main challenges associated with catalytic CO 2 reduction is product selectivity: the formation of carbon monoxide, molecular hydrogen, formate, methanol, and other products occurs with similar thermodynamic driving forces, making it difficult to selectively reduce CO 2 to the target product. Significant scientific effort has been aimed at the development of catalysts that can suppress the undesired hydrogen evolution reaction and direct the reaction toward the selective formation of the desired products, which are easy to handle and store. Inspired by natural photosynthesis, where the CO 2 reduction is achieved using NADPH cofactors in the Calvin cycle, in this study we explore biomimetic metal-free hydride donors as catalysts for the selective reduction of CO 2 to formate. Here, we outline our recent findings on the thermodynamic and kinetic parameters that control the hydride transfer from metal-free hydrides to CO 2 . By experimentally measuring and theoretically calculating the thermodynamic hydricities of a range of metal-free hydride donors, we derive structural and electronic factors that affect their hydride-donating abilities. Two dominant factors that contribute to the stronger hydride donors are identified to be (i) the stabilization of the positive charge formed upon HT via aromatization or by the presence of electron-donating groups and (ii) the destabilization of hydride donors through the anomeric effect or in the presence of significant structural constrains in the hydride molecule. Hydride donors with appropriate thermodynamic hydricities were reacted with CO 2 , and the formation of the formate ion (the first reduction step in CO 2 reduction to methanol) was confirmed experimentally, providing an important proof of principle that organocatalytic CO 2 reduction is feasible. The kinetics of hydride transfer to CO 2 were found to be slow, and the sluggish kinetics were assigned in part to the large self-exchange reorganization energy associated with the organic hydrides in the DMSO solvent. Finally, we outline our approaches to the closure of the catalytic cycle via the electrochemical and photochemical regeneration of the hydride (R–H) from the conjugate hydride acceptors (R + ). We illustrate how proton-coupled electron transfer can be efficiently utilized not only to lower the electrochemical potential at which the hydride regeneration takes place but also to suppress the unwanted dimerization that neutral radical intermediates tend to undergo. Overall, this account provides a summary of important milestones achieved in organocatalytic CO 2 reduction and provides insights into the future research directions needed for the discovery of inexpensive catalysts for carbon recycling.

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

Predictive Assessment of the Chemical Composition of Coal Ash in Reserve at U.S. Disposal Sites

In the United States, more than 2 Gt of coal combustion residuals (i.e., coal ash) are stored in hundreds of disposal units. Recent federal regulations mandate the closure or retrofitting of most coal ash impoundments, presenting significant challenges for waste management. These regulatory pressures also present opportunities to reuse coal ash. However, the quality and quantity of discarded coal ash across the U.S. are not well known, even though this information is crucial for spurring its reuse for conventional and new material applications. This study describes a predictive model for the major element composition of coal ash in reserve at disposal sites of major U.S. coal-fired power plants. This model was constructed from coal purchase records of 705 power stations from 1973 to 2022 and was trained on coal ash composition data, showing that coal ash elemental composition is strongly associated with the source of feedstock coal. The model showed regional shifts in the major element contents of ash produced by power plants in the last 50 years, particularly for calcium and iron (expressed as %CaO and %Fe2O3), as power stations changed their source of coal over this time frame. Our approach enables an estimation of chemical composition for ash stored in waste impoundments at individual power stations. Such information can help to delineate the regional market resource potential of supplementary cements for concrete and other material innovations that would utilize coal ash harvested from disposal sites across the U.S.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2020

What Factors Drive the Changes in Water Withdrawals in the U.S. Agriculture and Food Manufacturing Industries between 1995 and 2010?

Climate change and increasing world population will directly impact the global food supply chain linkages. In the United States, agricultural production requires less irrigated water than before but it still accounts for a third of total water withdrawals. To better understand the evolution of its water use, we perform a structural decomposition analysis of water withdrawals across eight different crops and six livestock categories and differentiate the trends over 1995–2005 vs 2005–2010 to account for the role of the economic crisis in the second period. Based on USGS data, the results show that both periods experienced an overall decline in water withdrawals in the production of all crops except oilseeds. This trend is driven by a decrease in water intensity, reflecting greater efficiency of irrigation systems, and by reduced local per capita income in the second period. However, increased foreign demand for water-intensive sectors like oilseeds from NAFTA and Asian partners mitigated the decline. Results indicate also a decreasing water use in livestock production partially due to a shift from red to white meat consumption in the country. Arguably, recent tariff wars and border closures have greatly reduced the virtual water embodied in American exports.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2020

Experimental Evidence of Dioxole Unimolecular Decay Pathway for Isoprene-Derived Criegee Intermediates

Ozonolysis of isoprene, one of the most abundant volatile organic compounds emitted into the Earth’s atmosphere, generates two four-carbon unsaturated Criegee intermediates, methyl vinyl ketone oxide (MVK-oxide) and methacrolein oxide (MACR-oxide). The extended conjugation between the vinyl substituent and carbonyl oxide groups of these Criegee intermediates facilitates rapid electrocyclic ring closures that form 5-membered cyclic peroxides, known as dioxoles. This research discusses the first experimental evidence of this novel decay pathway, which is predicted to be the dominant atmospheric sink for specific conformational forms of MVK-oxide (anti) and MACR-oxide (syn) with the vinyl substituent adjacent to the terminal O atom. The resulting dioxoles are predicted to undergo rapid unimolecular decay to oxygenated hydrocarbon radical products, including acetyl, vinoxy, formyl, and 2-methyl-vinoxy radicals. In the presence of O 2 , these radicals rapidly react to form peroxy radicals (ROO), which quickly decay via carbon-centered radical intermediates (QOOH) to stable carbonyl products that are identified in this work. The carbonyl products are detected under thermal conditions (298 K, 10 torr He) using multiplexed photoionization mass spectrometry (MPIMS). The main products (and associated relative abundances) originating from unimolecular decay of anti-MVK-oxide and subsequent reaction with O 2 are formaldehyde (88 ± 5%), ketene (9 ± 1%) and glyoxal (3 ± 1%). Those identified from the unimolecular decay of syn-MACR-oxide and subsequent reaction of O 2 are acetaldehyde (37 ± 7%), vinyl alcohol (9 ± 1%), methylketene (2 ± 1%), and acrolein (52 ± 5%). In addition to the stable carbonyl products, the secondary peroxy chemistry also generates OH or HO 2 radical co-products.

09 BIOMASS FUELS↗
DOE OSTI2025

A Combined Crossed Molecular Beam and Theoretical Investigation of the Elementary Reaction of Tricarbon (C 3 (X 1 Σ g + )) with Diacetylene (C 4 H 2 (X 1 Σ g + )): Gas Phase Formation of the Heptatriynylidyne Radical ( l -C 7 H(X 2 Π))

An elucidation of the underlying formation pathways to acyclic hydrocarbons such as polyynes (C n H 2 ), cumulenes (C n H 2 ), and linear resonantly stabilized linear radicals (l-C n H) is indispensable to understand the hydrocarbon chemistry in extreme low and high temperature environments. In this study, we exploited the crossed molecular beam technique to investigate the reaction of tricarbon C 3 (X 1 Σ g + ) with diacetylene (butadiyne; HCCCCH; X 1 Σ g + ) at a collision energy of 47 ± 1 kJ mol ⁻1 . The experimental data were merged with ab initio calculations of the singlet C 7 H 2 potential energy surface (PES) revealing that the reaction is initiated via the formation of an initial van der Waals reactant complex in the entrance channel. Subsequent rearrangements lead to various carbene-type and cyclic intermediates via ring-opening, ring-closure, and hydrogen migration processes eventually forming acyclic C 7 H 2 isomers prior to their barrierless unimolecular decomposition to the most stable linear isomer, heptatriynylidyne (C 7 H, X 2 Π) in an overall endoergic reaction (+57 kJ mol ⁻1 ). The reaction exhibits strong similarities to the tricarbon – acetylene (C 3 – C 2 H 2 ). Furthermore, the significant energy threshold suggests that the tricarbon reaction with (poly)acetylenes forming resonantly stabilized linear radicals are open in high-temperature environments such as combustion flames and circumstellar envelopes of carbon stars and planetary nebulae as their descendants; however, these reactions are closed in low-temperature environments as in cold molecular clouds and hydrocarbon-rich atmospheres of planets and their moons such as in Titan.

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

Nonadiabatic Dynamics of Photoexcited cis-Stilbene Using Ab Initio Multiple Spawning

The photochemistry of cis-stilbene proceeds through two pathways: cis-trans isomerization and ring closure to 4a,4b-dihydrophenanthrene (DHP). Despite serving for many decades as a model system for photoisomerization, the photodynamics of cis-stilbene is still not fully understood. In this work, we use ab initio multiple spawning on a SA-2-CASSCF(2,2) potential energy surface to simulate the nonadiabatic dynamics of isolated cis-stilbene. We find the cyclization (to DHP and cis-stilbene) and isomerization (to trans- and cis-stilbene) reaction coordinates to be orthogonal; branching between the two pathways is determined on the S1 excited state within 150 fs of photoexcitation. Furthermore ,trajectory basis functions (TBFs) undergoing cyclization decay rapidly to the ground state in 250 fs, while TBFs moving along the isomerization coordinate remain on the excited state longer, with the majority decaying between 300 and 500 fs. We observe three avoided crossing regions in the dynamics: two along the isomerization coordinate (displaying pyramidalization and migration of an ethylenic hydrogen or phenyl group), and one DHP-like conical intersection along the cyclization coordinate. The isomeric form of the vibrationally hot photoproducts (as determined by measurement 2 ps after photoexcitation) is determined within less than 50 fs of decay to the ground state mediated by passage through a conical intersection. Excess vibrational energy of ground state cis- and trans-stilbene is channelled into phenyl torsions (with mostly opposing directionality). Finally, our simulations are validated by direct comparison to experiment for the absorption spectrum, branching ratio of the three photoproducts (44:52:4 cis-stilbene:trans-stilbene:DHP), and excited state lifetime (520 ± 40 fs).

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

Comment on “Localized and Delocalized States of a Diamine Cation: Resolution of a Controversy”

Since its appearance in and recent re-investigation in the dimethylpiperazine cation (DMP + ) has generated considerable discussion and controversy in the scientific literature over the existence of stable, local energy minima in this molecular system. Specifically, prior assumptions that the Rydberg state and radical cation of DMP are similar have led to significant confusion and debate regarding the accuracy of various quantum chemistry methods and the existence of stable configurations of DMP+ itself. The purpose of this Viewpoint is to highlight recent studies that call into question the main findings in the previously mentioned works as well as present new CCSDT (Coupled-Cluster with Single, Double, and Triple excitations) calculations to finally bring closure to this controversy.

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

The Hidden Path to the Resonance-Stabilized Fulvenallenyl Radical (C 7 H 5 ) via the Bimolecular Reaction of Tricarbon (C 3 , X 1 Σ g + ) with 1,3-Butadiene (C 4 H 6 ; X 1 A g )

Here, we report a novel bimolecular reaction pathway forming the resonance-stabilized fulvenallenyl radical (C 7 H 5 ) via the gas-phase reaction of tricarbon (C 3 , X 1 Σ g + ) with 1,3-butadiene (C 4 H 6 , X 1 A g ). Crossed molecular beam experiments combined with high-level electronic structures and statistical calculations reveal a rich potential energy surface. The reaction proceeds via tricarbon addition to a double bond of 1,3-butadiene, overcoming a 30 kJ mol –1 barrier, followed by ring closure, isomerizations (ring opening/closures, hydrogen shifts), and eventual hydrogen atom loss, yielding the fulvenallenyl radical (p1) almost exclusively via two dominant pathways. Considering the entrance barrier, this reaction may occur in high-temperature environments, like circumstellar envelopes of carbon stars, but is suppressed in colder regions such as molecular clouds and Titan’s atmosphere. These findings highlight the unexpected reactivity of small carbon clusters in shaping the molecular complexity of our universe, from the flicker of a flame to the death of a star.

Aromatic compounds↗