Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “slip”

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 1,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Time domain probabilistic seismic risk analysis using ground motion prediction equations of Fourier amplitude spectra

Modeling of Fourier amplitude spectra (FAS) of seismic motions has gained much attention in engineering seismology. In the past few years, several ground motion prediction equations (GMPEs) and inter-frequency correlation structure of FAS have been established. Due to many preferable characteristics of FAS, probabilistic seismic hazard/risk analysis is rapidly changing from ergodic, spectrum acceleration Sa(T 0 )-based approach to non-ergodic, site-specific, FAS-based approach. This paper presents time domain intrusive framework for probabilistic seismic risk analysis using GMPE of FAS. Herein, methodology for time domain stochastic ground motion modeling based on GMPEs of FAS is presented in some detail. The simulated uncertain motions are modeled as a random process and represented by polynomial chaos Karhunen-Loève expansion. The random process excitations are further propagated into the uncertain structural system using Galerkin stochastic finite element method (SFEM). Probabilistic evolution of structural response is solved, and such solution is used to develop seismic risk for any damage state. The presented framework is illustrated through seismic risk analysis of a four-story building subjected to possible earthquakes from two strike slip faults. The influences of the epistemic uncertainties in source stress drop Δσ and site attenuation κ0 on seismic risk are investigated. The need for non-ergodic seismic risk analysis with source-specific and site specific characterizations is emphasized.

58 GEOSCIENCES↗
DOE OSTI2020

Measurement of flow properties coupled to experimental and numerical analyses of dense, granular flows for solar thermal energy storage

Granular flows of sintered bauxite proppants were researched along an inclined plane for solar thermal energy storage applications. Granular flow properties needed to drive numeric granular models were measured for improved numerical model predictions for Carbobead CP 50/140, 40/100, and 30/60 particles. Particle shape and size distributions were determined by coupling optical microscopy to an in-house image processing algorithm. The impulse excitation technique was used to measure elastic and shear moduli, and compute Poisson’s ratio. The coefficient of static sliding friction was measured using the slip-stick method, and the static rolling friction was determined from measured shear on particles positioned between two hot-pressed plates. The coefficient of restitution was measured by dropping particles on a surface and determining the kinetic energy before and after impact with the surface using high resolution particle tracking velocimetry. Particle size did not significantly impact the coefficients of restitution and static rolling friction, however, particle shape distribution resulted in a large variation in measurements. An inclined flow experiment was performed to characterize granular flows of Carbobead CP 30/60 particles using particle image velocimetry. Numerical models of the experiment using discrete element method were generated with the measured mechanical properties as inputs for comparison with experimental results. A constant directional torque rolling friction model best predicted bulk granular flow behavior. Furthermore, good agreement between the model and experiment was achieved at ambient, steady state conditions, with average velocity differences <10%.

14 SOLAR ENERGY↗
DOE OSTI2023

The effect of current and lambda on white-etch-crack failures

White etching cracks (WECs) have been associated with premature failure of wind turbine roller bearings. Various drivers for the generation of WECs have been identified such as loading conditions, slip, steel quality, lubrication, hydrogen embrittlement, corrosion fatigue cracking, and stray electrical currents passing through the surface. Here, in this work, a benchtop test rig utilizing a three-ring-on-roller test configuration was used to investigate the effect of electrical current and operation in different lubricating regimes, defined by lambda (λ), on high-quality bearing steel samples tested in a commercially available power transmission EP gear lubricant. It was observed that there is an inverse correlation between the magnitude of electric current applied to the ring/ roller system and time-to-failure. Higher current magnitudes lead to shorter time-to-failure than lower current magnitudes, with macropitting as the main failure mode. Sub-surface investigation revealed the presence of WECs in all cases. For the same current magnitude, tests conducted in boundary and mixed lubrication regimes showed that time-to-failure increased as lambda increased, and the tests resulted in WEC related macropits, whereas tests conducted in near-hydrodynamic regime resulted in surface damage with no macropit. It was also noted that a shift toward near-hydrodynamic lubrication resulted in a distinct surface distress on the roller surface. Furthermore, there seems to be a transition in the mixed regime during which the surface distress occurred. The damage on the surface of the test samples resembled non-spatially, periodic, groove-like corrugations and, in some cases, crater-like depressions. Sub-surface imaging, performed by sequential sectioning, revealed the presence of WECs in all cases, and broad, branching cracks that were more prevalent under the more severe boundary conditions.

17 WIND ENERGY↗
DOE OSTI2021

Observations of modal coupling due to bolted joints in an experimental benchmark structure

The goal of this paper is to present a set of measurements from a benchmark structure containing two bolted joints to support future efforts to predict the damping due to the joints and to model nonlinear coupling between the first two elastic modes. Bolted joints introduce nonlinearities in structures, typically causing a softening in the natural frequency and an increase in damping because of frictional slip between the contact interfaces within the joint. These nonlinearities pose significant challenges when characterizing the response of the structure under a large range of load amplitudes, especially when the modal responses become coupled, causing the effective damping and natural frequency to not only depend on the excitation amplitude of the targeted mode, but also the relative amplitudes of other modes. In this work, two nominally identical benchmark structures, known in some prior works as the S4 beam, are tested to characterize their nonlinear properties for the first two elastic modes. Detailed surface measurements are presented and validated through finite element analysis and reveal distinct contact interactions between the two sets of beams. The free-free test structures are excited with an impact hammer and the transient response is analyzed to extract the damping and frequency backbone curves. A range of impact amplitudes and drive points are used to isolate a single mode or to excite both modes simultaneously. Differences in the nonlinear response correlate with the relative strength of the modes that are excited, allowing one to characterize mode coupling. Each of the beams shows different nonlinear properties for each mode, which is attributed to the different contact pressure distributions between the parts, although the mode coupling relationship is found to be consistent between the two. Lastly, the test data key finding are presented in this paper and the supporting data is available on a public repository for interested researchers.

42 ENGINEERING↗
DOE OSTI2021

A quantitative assessment of the model form error of friction models across different interface representations for jointed structures

Hysteretic models are widely used to model frictional interactions in joints to recreate experimental behavior. However, it is unclear which models are best suited for fitting or predicting the responses of structures. Here, the present study evaluates 26 friction model/interface representation combinations to quantify the model form error. A Quasi-Static Modal Analysis approach (termed Rayleigh Quotient Nonlinear Modal Analysis) is adopted to calculate the nonlinear system response, and a Multi-Objective Optimization is solved to fit experimental data of the first mode of the Brake-Reuß Beam. Optimized parameters from the first mode are applied to the second and third bending modes to quantify the predictive ability of the models. Formulations for both tracing full hysteresis loops and recreating hysteresis loops from a single loading curve (Masing assumptions) are considered. Smoothly varying models applied to a five patch representation showed the highest flexibility (for fitting mode 1) and good predictive potential (for modes 2 and 3). For a second formulation, which uses 152 frictional elements to represent the interface, the physically motivated spring in series with a Coulomb slip model (elastic dry friction) has high error for fitting mode 1 and performs near the middle for predicting higher modes. For both interface representation, the best fit models are not the most physical, but rather the ones with the most parameters (as expected); however, the more physical models perform somewhat better for predicting the higher modes.

42 ENGINEERING↗
DOE OSTI2023

Multi-mode quasi-static excitation for systems with nonlinear joints

Finite element models can be used to model and predict the hysteresis and energy dissipation exhibited by nonlinear joints in structures. As a result of the nonlinearity, the frequency and damping of a mode is dependent on excitation amplitude, and when the modes remain uncoupled, quasi-static modal analysis has been shown to efficiently predict this behavior. However, in some cases the modes have been observed to couple such that the frequency and damping of one mode is dependent on the amplitude of other modes. To model the interactions between modes, one must integrate the dynamic equations in time, which is several orders of magnitude more expensive than quasi-static analysis. This work explores an alternative where quasi-static forces are applied in the shapes of two or more modes of vibration simultaneously, and the resulting load–displacement curves are used to deduce the effect of other modes on the effective frequency and damping of the mode in question. This methodology is demonstrated on a simple 2D cantilever beam structure with a single bolted joint which exhibits micro-slip nonlinearity over a range of vibration amplitudes. The predicted frequency and damping are compared with those extracted from a few expensive dynamic simulations of the structure, showing that the quasi-static approach produces reasonable albeit highly conservative bounds on the observed dynamics. This framework is also demonstrated on a 3D structure where dynamic simulations are infeasible.

42 ENGINEERING↗
DOE OSTI2021

Thermoosmotic Transport in Nanochannels Grafted with pH-responsive Polyelectrolyte Brushes Modelled Using Augmented Strong Stretching Theory

In this paper, we develop a theory to establish that the thermoosmotic (TOS) effects, induced by the application of an axial temperature gradient, leads to a massive enhancement in liquid transport in nanochannels grafted with charged polyelectrolyte (PE) brushes. Here, we quantify the TOS transport by quantifying the induced electric field and the induced TOS flow field. The different components of the electric field, namely the ionic component, the thermal component, and the osmotic component, as well as the contributions of different ions to these components, are quantified. Furthermore, we express the TOS velocity as a combination of chemiosmotic (COS), thermal, and electroosmotic (EOS) components. The COS and the thermal components augment each other and the overall strength and direction of the TOS flow are dictated by the direction and the relative strength of the EOS component. Most importantly, we compare the cases of brush-grafted nanochannels with those of the brush-free nanochannels of identical surface charge densities: the TOS transport is massively augmented in the brush-grafted nanochannels attributed to the combination of the localization of the EDL (and hence any body force that depends on the EDL charge density) away from the nanochannel wall (i.e., the location of the maximum drag force) and the presence of a possible molecular slip (experienced by the liquid) along the brush surface.

42 ENGINEERING↗
DOE OSTI2022

Reciprocal swimming at intermediate Reynolds number

In Stokes flow, Purcell's scallop theorem forbids objects with time-reversible (reciprocal) swimming strokes from moving. In the presence of inertia, this restriction is eased and reciprocally deforming bodies can swim. A number of recent works have investigated dimer models that swim reciprocally at intermediate Reynolds numbers Re≈1–1000. These show interesting results (e.g. switches of the swim direction as a function of inertia) but the results vary and seem to be case specific. Here, in this work, we introduce a general model and investigate the behaviour of an asymmetric spherical dimer of oscillating length for small-amplitude motion at intermediate Re. In our analysis we make the important distinction between particle and fluid inertia, both of which need to be considered separately. We asymptotically expand the Navier–Stokes equations in the small-amplitude limit to obtain a system of linear partial differential equations. Using a combination of numerical (finite element) and analytical (reciprocal theorem, method of reflections) methods we solve the system to obtain the dimer's swim speed and show that there are two mechanisms that give rise to motion: boundary conditions (an effective slip velocity) and Reynolds stresses. Each mechanism is driven by two classes of sphere–sphere interactions, between one sphere's motion and (1) the oscillating background flow induced by the other's motion, and (2) a geometric asymmetry induced by the other's presence. We can thus unify and explain behaviours observed in other works. Our results show how sensitive, counterintuitive and rich motility is in the parameter space of finite inertia of particles and fluid.

42 ENGINEERING↗
DOE OSTI2022

Floquet stability analyses of stratified oscillating boundary layers on adiabatic slopes

The presence of a no-slip, impermeable, adiabatic, sloped boundary in an otherwise quiescent, stably stratified, Boussinesq flow generates baroclinic vorticity within a diffusive boundary layer. Such conditions are typical of the oscillating boundary layers on adiabatic abyssal slopes, sloped lake bathymetry and sloped coastal bathymetry in the absence of high-wavenumber internal waves, mean flows, far-field turbulence on larger scales, and resonant tidal–bathymetric interaction. We investigate the linear stability of the oscillating flow within non-dimensional parameter space typical of the $M_2$ tide and hydraulically smooth, middle-latitude abyssal slopes through Floquet linear stability analysis. The flow dynamics depends on three non-dimensional variables: the Reynolds number for Stokes’ second problem, the Prandtl number, and a frequency ratio that accounts for the resonance conditions ($C$, criticality) of the buoyant restoring force and the tidal forcing. The Floquet analysis results suggest that oscillating laminar boundary layers on adiabatic abyssal slopes are increasingly unstable as Reynolds number, criticality parameter and/or spanwise disturbance wavenumber are increased. We also show that the two-dimensional Floquet linear instability necessarily generates three-dimensional baroclinic vorticity, which suggests that the evolution of the gravitational instabilities may be nonlinear as $t→∞$.

42 ENGINEERING↗
DOE OSTI2023

Beyond the Stokes approximation: shallow visco-elastic ice-sheet models

Abstract The hypothesis that ice-sheet evolution is only controlled by the long-term non-Newtonian viscous behavior of ice has been challenged by observations indicating that effects like brittle failure, stick-slip sliding, tides and wave action may affect ice-sheet evolution on sub-daily timescales. Over these timescales, the quasi-static-creep approximation is no longer appropriate and elastic effects become important. Simulating elastic effects in ice-sheet models over relevant timescales, however, remains challenging. Here, we show that by including a visco-elastic rheology and reintroducing the oft neglected acceleration term back into the ice-sheet stress balance, we can create a visco-elastic system where the velocity is locally determined and information propagates at the elastic wave speed. Crucially, the elastic wave speed can be treated like an adjustable parameter and set to any value to reproduce a range of phenomena, provided the wave speed is large compared to the viscous velocity. We illustrate the system using three examples. The first two examples demonstrate that the system converges to the steady-state viscous and elastic limits. The third example examines ice-shelf rifting and iceberg calving. This final example hints at the utility of the visco-elastic formulation in treating both long-term evolution and short-term environmental effects.

58 GEOSCIENCES↗
DOE OSTI2023

Tuning the Interaction between Platinum Single Atoms and Ceria by Zirconia Doping for Efficient Catalytic Ammonia Oxidation

Aiming at the development of an efficient NH 3 oxidation catalyst to eliminate the harmful NH 3 slip from the stationary flue gas denitrification system and diesel exhaust aftertreatment system, a facile ZrO 2 doping strategy was proposed to construct Pt 1 /Ce x Zr 1–x O 2 catalysts with a tunable Pt-CeO 2 interaction strength and Pt–O–Ce coordination environment. According to the results of systematic characterizations, Pt species supported on Ce x Zr 1–x O 2 were mainly in the form of single atoms when x ≥ 0.7, and the strength of the Pt-CeO 2 interaction and the coordination number of Pt–O–Ce bond (CN Pt–O–Ce ) on Pt 1 /Ce x Zr 1–x O 2 showed a volcanic change as a function of the ZrO 2 doping amount. Here, it was proposed that the balance between the reasonable concentration of oxygen defects and limited surface Zr–O x species well accounted for the strongest Pt-CeO 2 interaction and the highest CN Pt–O–Ce on Pt/Ce 0.9 Zr 0.1 O 2 . It was observed that the Pt/Ce 0.9 Zr 0.1 O 2 catalyst exhibited much higher NH 3 oxidation activity than other Pt/Ce x Zr 1–x O 2 catalysts. The mechanism study revealed that the Pt 1 species with the stronger Pt-CeO 2 interaction and higher CN Pt–O–Ce within Pt/Ce 0.9 Zr 0.1 O 2 could better activate NH 3 adsorbed on Lewis acid sites to react with O 2 thus resulting in superior NH 3 oxidation activity. This work provides a new approach for designing highly efficient Pt/CeO 2 based catalysts for low-temperature NH 3 oxidation.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Development of a Filtered CFD-DEM Drag Model with Multiscale Markers Using an Artificial Neural Network and Nonlinear Regression

Here, the accuracy of coarse-grained Euler-Lagrangian simulations of fluidized beds heavily depends on the mesoscale drag models to account for the influences of the unresolved sub-grid structures. Traditional filtered drag models are regressed with mesoscale markers such as voidage and slip velocities. In this research, a filtered drag was regressed with both mesoscale and macro-scale markers using fine grid Computational Fluid Dynamics - Discrete Element Method (CFD-DEM) simulations. The traditional non-linear regression method was compared with machine learning regression using an Artificial Neural Network (ANN) implemented in PyTorch and coupled with MFiX. The new drag showed higher accuracy than the Wen-Yu drag and another filtered drag derived from the two-fluid model. The nonlinear regression shows slightly better results than ANN regression in cases with similar R 2 values. The utilization of the gas inlet velocity as an additional macro-scale marker reduced the errors by up to 55.3% in the tested cases.

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

Balancing Charge Transfer and Frenkel Exciton Coupling Leads to Excimer Formation in Molecular Dimers: Implications for Singlet Fission

Photoexcitation of molecular chromophore aggregates can form excimer states that play a significant role in photophysical processes such as charge and energy transfer as well as singlet fission. An excimer state is commonly defined as a superposition of Frenkel exciton and charge transfer states. In this work, we investigate the dynamics of excimer formation and decay in π-stacked 9,10-bis(phenylethynyl)anthracene (BPEA) covalent dimers appended to a xanthene spacer, where the electronic coupling between the two BPEA molecules is adjusted by changing their longitudinal molecular slip distances. Using exciton coupling calculations, we quantify the relative contributions of Frenkel excitons and charge transfer states and find that there is an upper and lower threshold of the charge transfer contribution for efficient excimer formation to occur. In conclusion, knowing these thresholds can aid the design of molecular aggregates that optimize singlet fission.

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

Boosting Photoelectric Conductivity in Porphyrin-Based MOFs Incorporating C 60

Electronic structure calculations show that guest C60 in the porphyrin-containing metal organic frameworks Zn 2 (TCPB)(DA-ZnP) (DA-MOF; H 4 TCPB = 1,2,4,5-tetrakis(4-carboxyphenyl)benzene, DA-ZnP = [5,15-bis[(4-pyridyl)ethynyl]-10,20-diphenylporphinato]zinc(II)) and Zn 2 (TCPB)(F-ZnP) (F-MOF; F-ZnP = [5,15-di(4-pyridyl)-10,20-bis(pentafluorophenyl)porphinato]zinc(II)) engenders high photoelectrical conductivity due to efficient donor–acceptor charge-transfer (CT) interactions. Structural modifications at the meso position of the porphyrin influence the preferred positions of C 60 within the frameworks, giving rise to host–guest interactions with different anisotropic structural, electronic, and optoelectronic properties. A preferred slipped-parallel π-stacked interaction of C 60 that is predicted for NH 2 -substituted DA-MOF and F-MOF fosters strong CT transitions and lowers band gaps by ~1.0 eV compared to the pristine DA-MOF and F-MOF. Hopping rates computed using Marcus theory are found to be anisotropic and accelerated by multiple orders of magnitude across π-stacked interfaces created by C 60 incorporation, a consequence of strong electronic coupling between initial and final diabatic states. Calculations indicate that photoinduced electron transfer, as well as direct CT from porphyrin to C 60 upon irradiation, triggers a charge separation process that leads to the formation of what should be long-lived electron-trapped states at the heterojunctions. Lastly, design principles revealed here for the control of photophysical and electron-transfer processes will be useful for constructing new MOF-derived visible- and infrared-based optoelectronics.

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

Understanding Strain and Failure of a Knot in Polyethylene Using Molecular Dynamics with Machine-Learned Potentials

A neural network potential (NNP) has been developed by fitting to ab initio electronic structure data on hydrocarbons and is used to study failure of linear and knotted polyethylene (PE) chains. A linear PE chain must be highly strained before breaking as the stress is equally distributed across the chain. In contrast, the stress in a PE chain with a 31 or overhand knot, accumulates at the knot’s entrance/exit. We find the strain energy is greatest when the bond length and angle are strained simultaneously, and that the knot weakens the chain by increasing the variance of the C–C–C angle, thereby allowing rupture at lower bond strains. Here, we extend our analysis to both 51 and 52 knots and find that both break at the entrance/exit of a loop. Notably, molecular scale PE knots exhibit many of the same characteristics as knots in a macroscopic rope, with stick–slip phenomena upon tightening and similar points of failure.

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

Photoreactive Capture and Conversion of Dilute Carbon Dioxide into Synthetic Natural Gas

This study introduces a photoreactive system that integrates the capture of dilute CO 2 streams with their catalytic conversion to synthetic natural gas (CH 4 ), utilizing a Ru nanoparticle (NP)-doped TiO 2 composite loaded with linear polyethylenimine (L-PEI) and enhanced with plasmonic titanium nitride (TiN). This light-driven approach mitigates challenges that have plagued traditional thermal reactive carbon capture (RCC) methods, such as CO 2 slip and amine degradation. We demonstrate that L-PEI enables stable CO 2 capture and conversion, achieving ~70% conversion of captured CO 2 to CH 4 across multiple reaction cycles using nonflammable forming gas (~5% H 2 ) as the reductant. In contrast, branched PEI (B-PEI)-loaded composites exhibited significant catalyst deactivation after several RCC cycles. Scanning transmission electron microscopy (STEM) imaging confirms that significant sintering of the Ru NPs occur in the B-PEI sample under RCC conditions, whereas their size remains stable in more rigid L-PEI composites. Technoeconomic analysis (TEA) estimates that CH 4 production using this system could cost less than $\$$5/kg based on current electrocatalytic H 2 prices. These results represent one of the most promising demonstrations of amine-based RCC employing dilute CO 2 sources to date.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Reversible Interlayer Sliding and Conductivity Changes in Adaptive Tetrathiafulvalene-Based Covalent Organic Frameworks

Ordered interlayer stacking is intrinsic in two-dimensional covalent organic frameworks (2D COFs) and has strong implications on COF's optoelectronic properties. Reversible interlayer sliding, corresponding to shearing of 2D layers along their basal plane, is an appealing dynamic control of both structures and properties, yet it remains unexplored in the 2D COF field. In this paper, we demonstrate that the reversible interlayer sliding can be realized in an imine-linked tetrathiafulvalene (TTF)-based COF TTF-DMTA. The solvent treatment induces crystalline phase changes between the proposed staircase-like sql net structure and a slightly slipped eclipsed sql net structure. The solvation-induced crystallinity changes correlate well with reversible spectroscopic and electrical conductivity changes as demonstrated in oriented COF thin films. In contrast, no reversible switching is observed in a related TTF-TA COF, which differs from TTF-DMTA in terms of the absence of methoxy groups on the phenylene linkers. This work represents the first 2D COF example of which eclipsed and staircase-like aggregated states are interchangeably accessed via interlayer sliding, an uncharted structural feature that may enable applications such as chemiresistive sensors.

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

Improving and Understanding the Hydrogen Evolving Activity of a Cobalt Dithiolene Metal–Organic Framework

Despite the promising previous reports on the development of electrocatalytic dithiolene-based metal–organic frameworks (MOFs) for the hydrogen evolution reaction (HER), these materials often display poor reproducibility of the HER performance because of their poor bulk properties upon integration with electrode materials. In this work, we demonstrate an in-depth investigation of the electrocatalytic HER activity of a cobalt 2,3,6,7,10,11-triphenylenehexathiolate (CoTHT) MOF. To enhance the durability and charge transport properties of the constructed CoTHT/electrode architecture, CoTHT is deposited as an ink composite (1) composed of Nafion and carbon black. We leverage here the well-established use of catalyst inks in the literature to increase adhesion of the catalyst to the electrode surface and to improve the overall electrical conductivity of the integrated catalyst/electrode. The utilization of the composite 1 leads to a significant improvement in the overpotential (η) to reach a current density of 10 mA/cm 2 (η = 143 mV) compared to prior reports, resulting in the most active MOF-based electrocatalyst for the HER that contains only earth-abundant elements. Extensive density functional theory (DFT) calculations were applied to understand the structure of CoTHT and the mechanistic pathways of the HER. The computational results suggest that an AB stacking geometry is energetically favorable, where one layer is slipped by 1.6 Å relative to the neighboring one along the a and b vectors. Additionally, the DFT calculations indicate that the catalytic cycle likely involves a Volmer discharge step to generate a cobalt hydride, followed by a Heyrovsky step to form a cobalt–H 2 intermediate, and finally the dissociation of H 2 .

DFT calculations↗