Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonlinear correlation”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Heteroanionic Control of Exemplary Second-Harmonic Generation and Phase Matchability in 1D LiAsS 2– x Se x

The isostructural heteroanionic compounds β-LiAsS 2–x Se x (x = 0, 0.25, 1, 1.75, 2) show a positive correlation between selenium content and second-harmonic response and greatly outperform the industry standard AgGaSe 2 . These materials crystallize in the noncentrosymmetric space group Cc as one-dimensional $\frac{1}{∞}$ [AsQ 2 ] - (Q = S, Se, S/Se) chains consisting of corner-sharing AsQ3 trigonal pyramids with charge-balancing Li + atoms interspersed between the chains. LiAsS 2–x Se x melts congruently for 0 ≤ x ≤ 1.75, but when the Se content exceeds x = 1.75, crystallization is complicated by a phase transition. This behavior is attributed to the β- to α-phase transition present in LiAsSe 2 , which is observed in the Se-rich compositions. The band gap decreases with increasing Se content, starting at 1.63 eV (LiAsS 2 ) and reaching 1.06 eV (β-LiAsSe 2 ). Second-harmonic generation measurements as a function of wavelength on powder samples of β-LiAsS 2–x Se x show that these materials exhibit significantly higher nonlinearity than AgGaSe 2 (d 36 = 33 pm/V), reaching a maximum of 61.2 pm/V for LiAsS 2 . In comparison, single-crystal measurements for LiAsSSe yielded a d eff = 410 pm/V. LiAsSSe, LiAsS 0.25 Se 1.75 , and β-LiAsSe 2 show phase-matching behavior for incident wavelengths exceeding 3 μm. The laser-induced damage thresholds from two-photon absorption processes are on the same order of magnitude as AgGaSe 2 , with S-rich materials slightly outperforming AgGaSe 2 and Se-rich materials slightly underperforming AgGaSe 2 .

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

Necklace-structured high-harmonic generation for low-divergence, soft x-ray harmonic combs with tunable line spacing

The extreme nonlinear optical process of high-harmonic generation (HHG) makes it possible to map the properties of a laser beam onto a radiating electron wave function and, in turn, onto the emitted x-ray light. Bright HHG beams typically emerge from a longitudinal phased distribution of atomic-scale quantum antennae. Here, we form a transverse necklace-shaped phased array of linearly polarized HHG emitters, where orbital angular momentum conservation allows us to tune the line spacing and divergence properties of extreme ultraviolet and soft x-ray high-harmonic combs. The on-axis HHG emission has extremely low divergence, well below that obtained when using Gaussian driving beams, which further decreases with harmonic order. This work provides a new degree of freedom for the design of harmonic combs—particularly in the soft x-ray regime, where very limited options are available. Such harmonic beams can enable more sensitive probes of the fastest correlated charge and spin dynamics in molecules, nanoparticles, and materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

On the mechanism of ionization oscillations in Hall thrusters

Low-frequency ionization oscillations involving plasma and neutral density (breathing modes) are the most violent perturbations in Hall thrusters for electric propulsion. Because of its simplicity, the zero-dimensional (0D) predator–prey model of two nonlinearly coupled ordinary differential equations for plasma and neutral density has often been used for the characterization of such oscillations and scaling estimates. In this work, we investigate the properties of its continuum analog, the one-dimensional (1D) system of two nonlinearly coupled equations in partial derivatives (PDEs) for plasma and neutral density. This is a more general model, of which the standard 0D predator–prey model is a special limit case. We show that the 1D model is stable and does not show any oscillations for the boundary conditions relevant to Hall thrusters and the uniform ion velocity. We then propose a reduced 1D model based on two coupled PDEs for plasma and neutral densities that is unstable and exhibit oscillations if the ion velocity profile with the near-the-anode back-flow (toward the anode) region is used. Comparisons of the reduced model with the predictions of the full model that takes into account the self-consistent plasma response show that the main properties of the breathing mode are well captured. In particular, it is shown that the frequency of the breathing mode oscillations is weakly dependent on the final ion velocity but shows a strong correlation with the width of the ion back-flow region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

A Better Method to Calculate Fuel Burnup in Pebble Bed Reactors Using Machine Learning

Burnup measurement is an important step in material control and accountancy (MC&A) at nuclear reactors, and may be done by examining gamma spectra of fuel samples. Traditional approaches rely on known correlations to specific photopeaks (e.g. 137 Cs) and operate via a standard linear regression method. However, the quality of these regression methods is limited even in the best case, and is significantly poorer at short fuel cool-down times, due to the elevated radiation background by short life-time isotopes, and self-shielding effect of the fuel. For practical operation of pebble bed reactors (PBRs), quick measurements (in minutes) and short cooling times (in hours) are required from a safety and security perspective. We investigated the efficacy and performance of machine learning (ML) methods to predict the burnup of the pebble fuel from full gamma spectra (rather than specific discrete photopeaks) and found a full-spectrum ML approach to far outperform baseline regression predictions in all measurement and cooling conditions - including in operational-like measurement conditions. We also performed model and data ablation experiments to determine the relative performance impact of our ML methods' capacity to model data nonlinearities and the inherent additional information in full spectra. Applying our ML methods, we found a number of surprising results, including improved accuracy at shorter fuel cooling times (the opposite of the norm), remarkable robustness to spectrum compression (via rebinning), and competitive burnup predictions even when using background signal only (i.e. explicitly omitting known isotope photopeaks).

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2023

Dependence of the kinetic energy absorption capacity of bistable mechanical metamaterials on impactor mass and velocity

Using an alternative mechanism to dissipation or scattering, bistable structures and mechanical metamaterials have shown promise for mitigating the detrimental effects of impact by reversibly locking energy into strained material. Herein, we extend prior works on impact absorption via bistable metamaterials to computationally explore the dependence of kinetic energy transmission on the velocity and mass of the impactor, with strain rates exceeding 10 2 s –1 . We observe a large dependence on both impactor parameters, ranging from significantly better to worse performance than a comparative linear material. We then correlate the variability in performance to solitary wave formation in the system and give analytical estimates of idealized energy absorption capacity under dynamic loading. In addition, we find a significant dependence on damping accompanied by a qualitative difference in solitary wave propagation within the system. The complex dynamics revealed in this study offer potential future guidance for the application of bistable metamaterials to applications including human and engineered system shock and impact protection devices.

42 ENGINEERING↗
DOE OSTI2025

Improving Trustworthiness of Data-Driven Power Grid Contingency Analysis With Bayesian Residual Graph Neural Networks

The evolving energy landscape requires novel tools to efficiently perform contingency analysis and reliability assessment of power grids, potentially in real-time. The high computational cost of traditional power flow solvers limits their applicability in practice. Machine learning (ML) surrogates such as deep neural networks (NNs) accelerate power flow solvers computations, enabling high-order contingency analysis and real-time decision-making by learning highly nonlinear functions and integrating grid topology via graph architectures. However, (graph) NNs lack predictive power away from training data and do not provide predictive confidence estimates. Here, we present a Bayesian residual graph NN that integrates knowledge from low-fidelity data via residual training and embeds granular quantification of uncertainties, improving trustworthiness critical for high-consequence decision-making. Applying Bayesian concepts to NNs is challenging due to the high-dimensionality of both the parameter space, complicating derivation of a meaningful prior, and the output space in large grid systems, requiring enhanced techniques to assess the predicted high-dimensional uncertainties. Our contributions include: (1) Deriving a prior for fully connected and graph NNs that leverages low-fidelity data to guide mean predictions and appropriately control prior predictive uncertainty. (2) Integrating this prior within an ensembling with anchoring scheme for efficient approximate posterior inference. (3) Deriving enhanced metrics to assess accuracy of both the mean and uncertainty predictions in high dimensions, appropriately accounting for correlations propagated through graph layers. The resulting Bayesian residual graph NN is tested on a contingency analysis task for 14-bus and 118-bus grids.

24 - POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Opto-twistronic Hall effect in a three-dimensional spiral lattice

Studies of moire systems have elucidated the exquisite effect of quantum geometry on the electronic bands and their properties, leading to the discovery of new correlated phases. However, most experimental studies have been confined to a few layers in the 2D limit. The extension of twistronics to its 3D limit, where the twist is extended into the third dimension between adjacent layers, remains underexplored due to the challenges in precisely stacking layers. Here, in this study, we focus on 3D twistronics on a platform of self-assembled spiral superlattice of multilayered WS 2 . Our findings reveal an opto-twistronic Hall effect in the spiral superlattice. This mesoscopic response is an experimental manifestation of the noncommutative geometry that arises when translational symmetry is replaced by a non-symmorphic screw operation. We also discover signatures of altered laws of optical excitation, manifested as an unconventional photon momentum-lattice interaction owing to moire of moire modulations in the 3D twistronic system. Crucially, our findings mark the initial identification of higher-order quantum geometrical tensors in light-matter interactions. This breakthrough opens new avenues for designing quantum materials-based optical lattices with large nonlinearities, paving the way for the development of advanced quantum nanophotonic devices.

Ji, Zhurun↗
DOE OSTI2024

Nonlinear Lifshitz photon theory in condensed matter systems

We present an interacting theory of a U(1) gauge boson with a quadratic dispersion relation, which we call the "nonlinear Lifshitz photon theory.'' The Lifshitz photon is a three-dimensional generalization of the Tkachenko mode in rotating superfluids. Starting from the Wigner crystal of charged particles coupled to a dynamical U(1) gauge field, after integrating out gapped degrees of freedom, we arrive at the Lagrangian for the nonlinear Lifshitz photon. The symmetries of the theory include a global U(1) 1-form symmetry and nonlinearly realized "magnetic" translation and rotation symmetries. The interaction terms in the theory lead to the decay of the Lifshitz photon, the rate of which we estimate. We show that the Wilson loop, which plays the role of the order parameter of the spontaneous breaking of the 1-form global symmetry, deviates from the perimeter law by an additional logarithmic factor. Here, we explore potential connections to other condensed matter systems, with a particular focus on quantum spin ice and ferromagnets. Finally, we generalize our theory to higher dimensions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Robust cycloid crossover driven by anisotropy in the skyrmion host GaV 4 S 8

We report on the anomalous magnetization dynamics of the cycloidally modulated spin textures under the influence of uniaxial anisotropy in multiferroic GaV 4 S 8 . The temperature and field dependence of the linear ac susceptibility [$χ'_{1ω}$ (T, H)], AC magnetic loss [$χ''_{1ω}$ (T, H)], and nonlinear AC magnetic response [$M_{3ω}$ (T, H)] are examined across the magnetic phase diagram in the frequency range f = 10 - 10 000 Hz. According to recent theory, skyrmion vortices under axial crystal symmetry are confined along specific orientations, resulting in enhanced robustness against oblique magnetic fields and altered spin dynamics. We characterize the magnetic response of each spin texture and find that the dynamic rigidity of the Néel skyrmion lattice appears enhanced compared to Bloch-type skyrmions in cubic systems, even in the multidomain state. Anomalous M 3 ω and strong dissipation emerge over the same phase regime where strong variations in the cycloid pitch were observed on lowering temperature in recent small-angle neutron-scattering experiments [White et al., Phys. Rev. B 97, 020401(R) (2018)]. Here, we show that strong anisotropy also drives an extended crossover of the zero-field cycloid texture in GaV 4 S 8 . The frequency dependence of these dynamic signatures is consistent with that of a robust anharmonic spin texture exhibiting a correlated domain arrangement. The results underpin the essential role of magnetic anisotropy in enhancing the rigidity of topological spin textures for diverse applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Vibrational sum frequency generation spectroscopy reveals the inertness of chromium oxide (001) surfaces

Nanoengineered metal oxides such as Cr(III)-oxide (chromia) films have diverse potential applications in corrosion inhibition, remediation, energy generation, catalysis, data storage, and biological and environmental systems. Concerns about material degradation or oxidation to toxic chromate necessitate an understanding of chromia/aqueous interfaces, beginning with their hydroxylation and hydration behavior. Vibrational sum-frequency generation spectroscopy (vSFG) provides specific molecular-level information about water at the oxide/aqueous junction with high surface selectivity. To overcome the strong absorber problem typical of certain metal oxides in the UV–visible range during nonlinear optical studies, we employed molecular beam epitaxy to deposit transparent, epitaxial nanofilms of Cr 2 O 3 with (001) crystalline orientation on sapphire (Al 2 O 3 (001)) substrates, as confirmed by atomic force microscopy and X-ray diffraction. vSFG spectra of the air and water interfaces of the Cr 2 O 3 (001) films reveal hydroxyl features corresponding to both dissociated and molecular water on the surface. In contrast to the dangling hydroxyls found on the bare Al 2 O 3 (001) substrate, the hydroxyl groups on the deposited Cr 2 O 3 (001) nanofilm do not readily undergo isotopic H/D exchange when exposed to varying forms of D 2 O under ambient conditions. When considering the chemistries of the corresponding trivalent cations in aqueous solution, the finding is at variance with their similar acidities and proton exchange dynamics but consistent with markedly slower inner-sphere water ligand exchange of hexaquo Cr 3+ . This finding challenges the idea that proton and water exchange at the water/Cr₂O₃(001) interface are solely correlated, driven by the strength of metal-bridging oxygen bonds and surface hydroxyl distribution, without direct causation.

Boamah-Agyemang, Mavis D. [Pacific Northwest Natio↗
DOE OSTI2022

Optics-Free Chip-Scale Intraoperative Imaging Using NIR-Excited Upconverting Nanoparticles

Here we present an optics-free CMOS image sensor that incorporates a novel time-gated dual-photodiode pixel design to allow filter- and lens-less image acquisition of near-infrared-excited (NIR-excited) upconverting nanoparticles. Recent biomedical advances have highlighted the benefits of NIR excitation, but NIR interaction with silicon has remained a challenge, even with high-performance optical blocking filters. Using a secondary diode and a dual-photodiode design, this sensor is able to remove the 100s of mV of NIR background on pixels and bring it down to single-digit mV level, nearing its noise floor of 2.2 mV rms, not achievable with any optical filter. Non-linear effects of background cancellation using the diode pair has been mitigated using an initial one-time pixel-level curve fitting and calibration in a post-processing setting. This imager comprises a highly linear 11 fF metal-oxide-metal (MOM) capacitor and includes integrated angle-selective gratings to reject oblique light and enhance sharpness. Each pixel also includes two distinct correlated double sampling schemes, to remove low frequency flicker noise and systematic offset in the datapath. We demonstrate the performance of this imager using pulsed NIR-excited upconverting nanoparticles on standard United-States-Air-Force (USAF) resolution targets and achieve an SNR of 15 dB, while keeping NIR background below 6 mV. This 36-by-80-pixel array measures only 2.3 mm by 4.8 mm and can be thinned down to 25 µm, allowing it to become surgically compatible with intraoperative instruments and equipment, while remaining optics-free.

42 ENGINEERING↗
DOE OSTI2023

BaCu 2 SiS 4 : A New Member of the A II B I 2 M IV Q 4 Chalcogenide Family with a Chiral Crystal Structure

Abstract Noncentrosymmetric ternary and quaternary chalcogenides are studied as promising nonlinear optical (NLO) materials in the mid‐infrared region. Here, we report the synthesis of a new material BaCu 2 SiS 4 in the A II B I 2 M IV Q 4 family ( A =divalent metal; B =monovalent metal; M =tetrel, Q =chalcogen), and discuss its crystal structure, thermal stability, optical behavior, and electronic structure. BaCu 2 SiS 4 crystallizes in the noncentrosymmetric chiral space group P 3 2 21 with lattice parameters a =6.1440(3) Å, c =15.3542(8) Å, V =501.95(6) Å 3 , Z =3. The structure features helical channels formed by corner‐sharing [CuS 4 ] and [SiS 4 ] tetrahedral units. Synthesis was carried out in a molten salt flux, as opposed to a traditional solid‐state route from elements, to minimize the formation of a competing ternary phase, Ba 2 SiS 4 . BaCu 2 SiS 4 is a semiconductor with an experimentally‐determined direct bandgap of ~2.2 eV. The material exhibits second harmonic generation (SHG) activity, confirming the noncentrosymmetric nature of the structure. Analysis of reported A II B I 2 M IV Q 4 crystal structures pointed out a correlation among potential structure types and the radii of the constituent elements. Total energy calculations were carried out to explore the relative stability of several reported crystal structures in this family of compounds.

Sarkar, Arka↗
DOE OSTI2025

Understanding latent timescales in neural ordinary differential equation models of advection-dominated dynamical systems

The neural ordinary differential equation (ODE) framework has shown considerable promise in recent years in developing highly accelerated surrogate models for complex physical systems characterized by partial differential equations (PDEs). For PDE-based systems, state-of-the-art neural ODE strategies leverage a two-step procedure to achieve this acceleration: a nonlinear dimensionality reduction step provided by an autoencoder, and a time integration step provided by a neural-network based model for the resultant latent space dynamics (the neural ODE). This work explores the applicability of such autoencoder-based neural ODE strategies for PDEs in which advection terms play a critical role. More specifically, alongside predictive demonstrations, physical insight into the sources of model acceleration (i.e., how the neural ODE achieves its acceleration) is the scope of the current study. Such investigations are performed by quantifying the effects of both autoencoder and neural ODE components on latent system time-scales using eigenvalue analysis of dynamical system Jacobians. To this end, the sensitivity of various critical training parameters – de-coupled versus end-to-end training, latent space dimensionality, and the role of training trajectory length, for example – to both model accuracy and the discovered latent system timescales is quantified. Furthermore, this work specifically uncovers the key role played by the training trajectory length (the number of rollout steps in the loss function during training) on the latent system timescales: larger trajectory lengths correlate with an increase in limiting neural ODE time-scales, and optimal neural ODEs are found to recover the largest time-scales of the full-order (ground-truth) system. Demonstrations are performed across fundamentally different unsteady fluid dynamics configurations influenced by advection: (1) the Kuramoto–Sivashinsky equations (2) Hydrogen-Air channel detonations (the compressible reacting Navier–Stokes equations with detailed chemistry), and (3) 2D Atmospheric flow.

Advection-dominated dynamical systems↗
DOE OSTI2022

Effective drift velocity from turbulent transport by vorticity

Here, we highlight the differing roles of vorticity and strain in the transport of coarse-grained scalars at length scales larger than ℓ by smaller scale (subscale) turbulence. We use the first term in a multiscale gradient expansion due to which exhibits excellent correlation with the exact subscale physics when the partitioning length ℓ is any scale smaller than that of the spectral peak. We show that unlike subscale strain, which acts as an anisotropic diffusion/antidiffusion tensor, subscale vorticity's contribution is solely a conservative advection of coarse-grained quantities by an eddy-induced nondivergent velocity, v * , that is proportional to the curl of vorticity. Therefore, material (Lagrangian) advection of coarse-grained quantities is accomplished not by the coarse-grained flow velocity, u¯ ℓ , but by the effective velocity, u¯ ℓ + v * , the physics of which may improve commonly used LES models.

42 ENGINEERING↗
DOE OSTI2025

Exploring the transition from continuous turbulence fluctuations to bursting ELMs in high SOL density regimes

BOUT++ turbulence simulations of the DIII-D reveal that the density profile between the separatrix and pedestal plays a crucial role in the dynamics of edge localized modes (ELMs) and edge plasma turbulent transport. Nonlinear simulations demonstrate that small ELMs in the DIII-D hybrid scenario under high SOL density conditions are predominantly driven by local ballooning modes near the separatrix, stabilizing global instabilities while enhancing localized pressure fluctuations. A key control parameters for ELM dynamics is the separatrix-to-pedestal density ratio, n e,sep /n e,ped . A high ratio indicates a shallow gradient, favoring small ELMs, while a lower ratio signals a steep gradient, which increases the likelihood of large ELMs. Comprehensive parameter scans, including n e,sep /n e,ped , density gradient profiles near the separatrix, and resistivity, reveal the critical role of these parameters in shaping transitions between turbulence-driven transport and ELM bursting. The scans demonstrate that in high SOL density regimes, small ELMs can result from either global resistive MHD instabilities or local ballooning modes near the separatrix, depending on the steepness of the separatrix density gradient. These findings also highlight the transition from continuous turbulence to small ELMs. The post-crash peak in pressure fluctuations, δP rms serves as a critical metric for identifying transition from continuous turbulence fluctuations to ELM bursting. Larger δP rms values correlate with ELM bursts driven by local or global instabilities, whereas smaller values indicate turbulence-dominated transport. Drift-Alfvén and resistive ballooning turbulence enhance the entrainment of fluctuations from the pedestal to the SOL, contributing to the complex interplay of dynamics in this regime. These findings emphasize the importance of separatrix density shaping and pedestal gradient control for optimizing ELM behavior in ITER and future fusion devices.

Li, Nami [Lawrence Livermore National Laboratory (↗
DOE OSTI2020

Quantum scars of bosons with correlated hopping

Abstract Recent experiments on Rydberg atom arrays have found evidence of anomalously slow thermalization and persistent density oscillations, which have been interpreted as a many-body analog of the phenomenon of quantum scars. Periodic dynamics and atypical scarred eigenstates originate from a “hard” kinetic constraint: the neighboring Rydberg atoms cannot be simultaneously excited. Here we propose a realization of quantum many-body scars in a 1D bosonic lattice model with a “soft” constraint in the form of density-assisted hopping. We discuss the relation of this model to the standard Bose-Hubbard model and possible experimental realizations using ultracold atoms. We find that this model exhibits similar phenomenology to the Rydberg atom chain, including weakly entangled eigenstates at high energy densities and the presence of a large number of exact zero energy states, with distinct algebraic structure.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Local charge-screening and polarization switching in a ferroelectric ionic-superconductor

For ferroelectric ionic-conductors, polarization switching is complicated by the interplay between ion mobility and charge screening effects. When the ionic charge carriers also play a key role in the domain reversal, such as in Rb-doped KTiOPO 4 (RKTP), a higher level of complexity is introduced. RKTP provides an ideal platform for investigating the relationship between ionic conductivity and polarization reversal because its highly anisotropic crystal properties allow selective modification of material characteristics through diffusive cation doping. Here, we use indiffused Ba/K doping to create a significant increase in the ionic conductivity. Time-of-flight secondary ion mass spectrometry is employed to map Ba/K doping within the RKTP crystal and correlate it to changes in ionic mobility and polarization switching characteristics under an external field applied to the nonpolar face. Using band-excitation piezoresponse force microscopy, we demonstrate a selective switching-inhibition mechanism driven by the enhanced charge screening.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Quantifying transport and electrocatalytic reaction processes in a gastight rotating cylinder electrode reactor via integration of Computational Fluid Dynamics modeling and experiments

Understanding the complexity of the multiple processes of mass, momentum, charge, and heat transport, and how these affect reaction kinetics at the electrode/electrolyte interface is one of the major challenges in the field of energy and catalysis. The rapid and rational scale-up of electrocatalytic systems to industrial scales require a detailed understanding of nonlinear transport-reaction processes, accessible only through the building of multi-physics models that capture with high fidelity the complexity of real-world devices. The gastight rotating cylinder electrode (RCE) reactor is a promising lab-scale tool that can decouple transport from intrinsic kinetics to generate data for first-principle models useful in the design of industrial, electrochemical reactors. Computational Fluid Dynamics (CFD) studies have previously been used to investigate the bulk flow in RCE reactors for simple corrosion and electroplating processes. However, the quantification of changes in local concentration within the viscous layer where catalysis takes place requires capturing the correct flow conditions inside the hydrodynamic boundary layer near the surface of the electrode. Further, this requires simulations with spatial resolution in the nm and μm scale and temporal resolutions between ms and s scales that are similar to the timescales for reactions on the electrode surface. In this study, experimental electrocatalysis is combined with CFD modeling to elucidate and parameterize the hydrodynamics in a gastight RCE reactor. CFD simulations of the electrochemical ferricyanide reduction reaction under mass transport limited conditions are used to evaluate the validity of the CFD model parameters by comparing calculated dimensionless mass transport descriptors to dimensionless correlations obtained experimentally. Justifications for assumptions and details of the simulation methods used in this study are presented to provide a detailed understanding of the effect that each model parameter has on the ability to accurately simulate electrocatalysis in RCE systems. The simulation methodology reported here is a first step towards the development of multi-scale models for the study of transport dependent electrocatalytic processes, such as the electrochemical transformation of CO 2 to fuels and chemicals.

42 ENGINEERING↗