Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Fluid interfaces”

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 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Advanced Computational Modeling of High-Level Waste Vitrification at the Hanford Site

The U.S. Department of Energy (DOE) has selected vitrification for stabilizing legacy tank waste at the Hanford site, where radioactive waste from plutonium production was historically stored in underground tanks. This waste will be separated into low-activity waste (LAW) and high-level waste (HLW) fractions and processed at the Waste Treatment and Immobilization Plant (WTP). At WTP, glass melters are used for the vitrification of radioactive tank waste, transforming it into a stable borosilicate glass form for safe long-term storage. The melter vessel is constructed from highly durable and heat-resistant materials, where the vitrification process occurs. The main regions that are modeled are the melt pool, plenum, cold cap, riser/discharge chamber, and surrounding structure with insulation layers. Forced convection induced by air bubblers at the base of the melter ensure uniform temperature distribution and provide heat to the cold cap layer. The cold cap is a region of reacting batch feed that floats on top of the molten glass and is where the batch-to-glass reactions occur. Joule heating provided by electrodes mounted along the vertical walls of the melter and immersed directly in the glass, generates the necessary heat for the net endothermic conversion processes that occur in the cold cap. The high temperatures, radioactivity, and opaque nature of the glass prevent direct observation inside the melters. Therefore, computational models are essential for providing insight into factors that affect melter throughput. Thermocouples in the plenum provide operators with plenum temperature measurements. Operational adjustments include bubbling rate, voltage supplied to the electrodes, feed adjustments, and glass removal rate. Different computational fluid dynamics (CFD) models have been developed, each serving a specific purpose. There are CFD models of different scale melters, as well as models that capture the two-phase flow interfaces of rising bubbles in the molten glass or models with a simplified molten glass region so that the surrounding structure and plenum can be feasibly incorporated. Pilot-scale melter models have been developed to serve as validation of the methods employed in the simulation of the full-scale WTP melters. Models incorporating resolved bubbling are used to develop momentum source terms to implement into a single phase, multi-region, steady-state flow model that is being validated by measured process parameters such as glass production rate, voltage, input power, plenum temperatures, etc. The resolved bubbling model uses the multiphase volume of fluid approach to model the system with a high-resolution interface capturing scheme to maintain sharp interfaces between the molten glass and the air phase. The suite of CFD models is continually being improved to incorporate more realistic physics and achieve faster turnaround time. For example, an incremental controller is implemented to automatically adjust electrode voltage within the simulation to a molten glass set point temperature of 1150°C. Newer models feature improved meshes to ensure conformal meshes between regions and eliminate unnecessary mesh refinement in areas that are not of interest (such as boundary layers in offgas ports). Instead of explicitly modeling the structural, refractory, and insulation layers of the melter, a thermal resistance approach is used with published correlations used for boundary conditions. The development of robust and efficient CFD models will be instrumental in enabling the WTP to successfully fulfill its mission of safely stabilizing legacy nuclear waste.

12 - MGMT OF RADIOACTIVE AND NON-RADIOACTIVE WASTE↗
NASA NTRS1992

CONE - An STS-based cryogenic fluid management experiment

An overview of the CONE program is presented which includes a definition of the technology addressed by CONE and a baseline experiment set, a description of the experimental and support subsystems, interface requirements between the STS and the experiment carrier (Hitchhiker M), and the reusability and expansion capacity for additional experiment flights. CONE evaluates three primary technologies: the active thermodynamic vent system, the passive thermodynamic vent system, and liquid acquisition device performance. The cryogenic fluid management technology database that the system offers will allow for efficient subcritical cryogenic system designs for operation in a low-gravity environment. This system maximizes the balance between existing component technology and the need for the development of a cryogenic-fluid-management (CFM) test bed to investigate and demonstrate methods of storage and handling arenas.

Bell, R. S.↗
NASA NTRS1994

Dynamic Characteristics of the Partially Filled Rotating Dewar of the Gravity Probe-B Spacecraft

The dynamical behavior of fluids, in particular the effect of surface tension on partially-filled rotating fluids (cryogenic liquid helium and helium vapor) in a full scale Gravity Probe-B Spacecraft propellant dewar tank imposed by various frequencies of gravity jitters have been investigated. The study of the liquid-vapor interface oscillations due to various frequencies of gravity jitter under different dewar rotating speeds and different levels of background gravity have been numerically simulated. Results disclose the time sequence evolution for the excitation of large amplitude oscillation waves at the liquid-vapor interface. The study shows that slosh waves excited inside the spacecraft propellant tank are characterized by the lowest frequency of the waves initiated, frequencies of the gravity jitters imposed on the propellant system, the levels of background gravity environment, and dewar rotating speeds. Conditions for suppression and amplification of the slosh waves are discussed. It also shows that fluid stress distribution exerted on the walls of the rotating dewar are closely related to the characteristics of slosh waves excited on the liquid-vapor interface in the rotating dewar tank. This can provide a set of data leading toward the control of spacecraft imbalance caused by the uneven fluid stress distribution from slosh waves.

Hung, R. J.↗
DOE OSTI2023

Utah FORGE: Fault Shear Reactivation Experimental Data for Fluid Injection-Rate Controls on Seismic Moment

Included are experimental data recorded from shear experiments that specifically explore the link between fluid-injection rate and seismic moment resulting from shear reactivation of laboratory faults. Raw mechanical data from three experiments are included alongside corresponding MATLAB scripts that import and plot the data, as well as use it to calculate shear and normal stress. Experiments are performed on 2.5-3 inch long granitoid cores from the Utah FORGE EGS demonstration site, containing a single inclined fracture with small-scale roughness added to the fracture surface. The raw data included here were recorded from an aluminum triaxial pressure vessel (TEMCO) configured with three independent servo-controlled pumps, with distilled water used as the working fluid. The pumps control confining pressure, upstream pore pressure, and axial pressure, with each pump connected to a LabView interface to record applied pressures, cumulative injected water volumes, and pump flow rates. The downstream outlet from the fracture is closed to allow pressurization, which is measured by an external pressure transducer. A linear variable differential transformer (LVDT) attached to the axial piston measures axial displacement, from which we calculate shear displacement along the fracture. Additionally, P-wave transducers are used to record acoustic signatures, where acoustic emission events and maximum amplitudes are compared against seismic moment and shear slip velocity. Fluid injection rates range between 0.05 mL/min, 0.25 mL/min, and 0.75 mL/min for each experiment. Along-fault pressure distributions are progressively less uniform as injection rates increase, representing a switch from steady-state to transient conditions. Triggered shear displacement is used as a proxy for seismic moment and is indexed against cumulative injection volume and rate. Each experiment is performed under constant shear stress conditions, and the sample is fully saturated with DI water. Axial and confining stresses are applied to 3 MPa through pressure-stepping in 500 kPa increments. The pore pressure is held constant at 200 kPa prior to initiating the experiment, and initial axial displacement is recorded. The axial stress is then increased to initiate shear mobilization during the loading phase (run-in) until a peak steady state is achieved. The initial shear stress is reduced to approximately 80% of the peak shear stress by decreasing the axial stress, then held constant for the duration of each experiment.

15 GEOTHERMAL ENERGY↗
DOE OSTI2026

Shear-free, inhomogeneous turbulence in a stably stratified fluid

High-resolution large eddy simulations are conducted of locally forced, shear-free turbulence in the presence of an initially sharp density interface. The simulations are reminiscent of oscillating grid turbulence experiments used to isolate the effect of turbulent diffusion and entrainment from background shear. By simulating such a flow we avoid common challenges of the experiments such as secondary-flow contamination due to sidewall effects and the inevitable interaction of the stratifying agent and forcing region. To address the latter concern, we add a heating term (potential energy sink) to the governing equations in the forcing layer, thereby preventing a heat flux through the source region. This modification sets up a continuous stratification in the mixed layer that is often assumed to be negligible in experiments. Despite this difference, we are able to make meaningful comparisons in terms of the overall entrainment rate, which varies as a power law with a turbulent Richardson number. Two exponents, −2 and −1, are measured depending on the definition of the Richardson number and entrainment rate used. The definition leading to −1 is consistent with most experiments, and we argue it is the superior choice if one is able to measure the relevant quantities. We also verify the self-similar scaling of turbulence velocity and length scales in the homogeneous fluid and propose ‘inner’ and ‘outer’ scalings for the stratified cases based on a local Froude number. The detailed scaling results are useful for turbulence model validation.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1981

Study of Thermal Control Systems for orbiting power systems

Thermal control system designs were evaluated for the 25 kW power system. Factors considered include long operating life, high reliability, and meteoroid hazards to the space radiator. Based on a cost advantage, the bumpered pumped fluid radiator is recommended for the initial 25 kW power system and intermediate versions up to 50 kW. For advanced power systems with heat rejection rates above 50 kW the lower weight of the advanced heat pipe radiator offsets the higher cost and this design is recommended. The power system payloads heat rejection allocations studies show that a centralized heat rejection system is the most weight and cost effective approach. The thermal interface between the power system and the payloads was addressed and a concept for a contact heat exchanger that eliminates fluid transfer between the power system and the payloads was developed. Finally, a preliminary design of the thermal control system, with emphasis on the radiator and radiator deployment mechanism, is presented.

Howell, H. R.↗
DOE OSTI2024

Heat transfer and transport property contrast effects on the compressible Rayleigh-Taylor instability

In extreme environments, the Rayleigh-Taylor instability (RTI) may occur under large variations in density and temperature and with fluid transport properties strongly dependent on temperature. Direct numerical simulations of the 3D fully compressible RTI are conducted, examining the idealized configuration of a hotter, less dense fluid pushing against a colder, denser fluid. Various temperature ratios and transport property configurations are explored to examine how heat conduction, large variations in transport properties, and sudden changes in transport properties can affect the evolution of the mixing layer. Nonuniform fluid expansion and contraction induced by heat transfer can significantly affect local density differences and overall instability growth, causing profile asymmetries about the initial interface for flow and mixing statistics. The departures from classical self-similar development of the instability along with misalignment between regions of mixing and regions of most intense turbulent activity caused by both heat transfer and transport property contrasts are examined. After sudden changes in fluid transport properties, which may occur as a result of rapid heating (e.g., in inertial confinement fusion), the flow quickly responds and begins to relax towards quasi-self-similar late-time evolution. For many dynamical quantities such as vorticity and dissipation, this late-time evolution resembles that of the configuration that already started with the final transport property magnitudes, suggesting that these quantities depend only on the transport properties and not on past flow history, provided that the density field distributions for the flows remain similar. On the other hand, the mixing evolution after the transport property change is unique, implying that both property magnitudes and previous history are impactful on the mixing. In conclusion, these simulations demonstrate how various temperature-related effects are extremely important to consider in compressible RTI flows with large temperature variations.

42 ENGINEERING↗
DOE OSTI2025

Incorporating the molecular-scale into a hydrodynamic description of confined aqueous systems

Hydrodynamics provides a continuum-level description of fluid motion, but its applicability at the nanoscale becomes uncertain due to the emerging importance of molecular-level effects such as spatial heterogeneity. Hydrodynamic boundary conditions that incorporate molecular details allow us to partition the system into a near-wall region and a bulk fluid region. We identify a hydrodynamic wall located inside the fluid that determines where slip begins. By extending the hydrodynamic wall with the slip length, the position of the extrapolated wall is established. This offers a unified description of both slip and stagnant flow behaviors, with wall hydrophobicity characterized by the relative location of the extrapolated wall with respect to the physical wall. Employing this concept in analyses of equilibrium molecular dynamics (MD) and non-equilibrium MD simulations of Couette and Poiseuille flows, our results demonstrate consistency between equilibrium and non-equilibrium approaches across different flow types and confinement levels. This demonstrates the robust nature of linear response theory. We then explore the effects of fluid-wall and bulk fluid interactions on the hydrodynamic properties. These findings enhance the effectiveness of molecular-based simulations for investigating complex confined systems in nanofluidics, biology, and colloidal science, offering a complementary molecular-scale perspective to traditional continuum approaches.

flow simulations↗
DOE OSTI2024

SMART Deliverable 6.1.2a: Application of the ORION tool to the IBDP Carbon Storage Site

Forecasting and managing potential induced seismic activity is one of the challenges facing commercialscale geologic carbon sequestration (GCS), as well as other geologic energy extraction and byproduct disposal technologies. Historically, the process to develop robust, science-based forecasts of induced seismicity has required an integrated effort from experts in seismology, geomechanics, and reservoir engineering to manage data, develop and evaluate models of subsurface processes, and to calibrate and interpret the results from a range of models to understand site behavior relative to prescribed standards and in the context of uncertainty in geologic characterization data, forecasting models, and operational scenario uncertainty. The Operational Forecasting of Induced Seismicity (ORION) toolkit is an open-source, observation-based forecasting toolkit that is being co-developed by two U.S. DOE-funded initiatives: the National Risk Assessment Partnership (NRAP) and the Scienceinformed Machine Learning for Accelerating Real Time Decisions in Subsurface Applications (SMART) Initiative. ORION is designed to provide functionality to support decision making about seismic hazard analysis and risk management for GCS stakeholders ranging from the public to site operators to expert seismologists. The tool, which is written as open-source code in the Python programming language, is composed of a desktop graphical user interface (GUI) and an underlying forecasting engine. The forecasting engine uses available reservoir properties, well and fluid injection scenario details, and observed seismic catalog data as inputs to produce a set of temporal and spatio-temporal seismic forecasts.

58 GEOSCIENCES↗
NASA NTRS1981

Numerical modeling of diffusive-convective physical vapor transport in cylindrical vertical ampoules

Diffusive-convective physical vapor transport (PVT) in cylindrical, vertical ampoules of aspect ratio (length/radius) between 0.5 and 10 was modeled numerically. The transport of a crystal forming component through an inert component that undergoes zero net transport was considered. Systems were treated in which: (1) with unequal molecular weight of the components and with temperature gradients typically employed in PVT, convective flow arises dominantly from solutal density gradients; and (2) with equal molecular weight of the components, convective flow can arise only from thermal expansion. It was found that, due to the diffusion-induced horizontal density gradients, buoyancy-driven convective flows are superimposed on the diffusive-advective fluxes without threshold. Net recirculation sets in adjacent to the growing interface, in contrast to the corresponding monocomponent situation where marginally stable convective modes fill the whole fluid space. Depending on the orientation of the main transport direction with respect to gravity, convection can either reduce or enhance the diffusion-induced radial concentration gradients. Significant enhancement of the net transport rate was found to occur only when the whole vapor space between source and growing crystal is filled by a convective recirculation roll. Solutal and thermal convection results are similar; yet for quantitative discussions, thermal and solutal Rayleigh numbers are not interchangeable in contrast to convective situations that lack net mass transport across the fluid space.

Markham, B. L.↗
NASA NTRS1988

Space erectable radiator system development

The NASA Space Station's Space Erectable Radiator System features modular radiator panels with high-capacity tapered artery heat pipes bonded within their aluminum honeycomb structures. Simple, dry aluminum-to-aluminum thermal contact surfaces are used for the connections, thereby requiring no fluid joints; a uniformly distributed clamping force at the radiator panel-to-thermal transport loop interface heat exchange surface yields high thermal contact conductance as well as minimum area and weight for the requisite performance. The design has been optimized for weight and cost.

Oren, J. A.↗
NASA NTRS1986

Probes and monitors for the study of solidification of molten semiconductors

The purpose is to examine solidification in the LiCl-KCl system to determine if phenomena such as solute rejection can be obseved by laser schlieren imaging. Molten salts have attributes that make them attractive as physical models in solidification studies. With optical techniques of investigation such as schlieren imaging, it is possible to study fluid flow phenomena in molten salts and to watch the trajectory of the solid-liquid interface.

Sadoway, D. R.↗
NASA NTRS2017

Sample Processor for Life on Icy Worlds (SPLIce): Design and Test Results

We report the design, development, and testing of the Sample Processor for Life on Icy Worlds (SPLIce) system, a microfluidic sample processor to enable autonomous detection of signatures of life and measurements of habitability parameters in Ocean Worlds. This monolithic fluid processing-and-handling system (Figure 1; mass 0.5 kg) retrieves a 50-L-volume sample and prepares it to supply a suite of detection instruments, each with unique preparation needs. SPLIce has potential applications in orbiter missions that sample ocean plumes, such as found in Saturns icy moon Enceladus, or landed missions on the surface of icy satellites, such as Jupiters moon Europa. Answering the question Are we alone in the universe? is captivating and exceptionally challenging. Even general criteria that define life very broadly include a significant role for water [1,2]. Searches for extinct or extant life therefore prioritize locations of abundant water whether in ancient (Mars), or present (Europa and Enceladus) times. Only two previous planetary missions had onboard fluid processing: the Viking Biology Experiments [3] and Phoenixs Wet Chemistry Laboratory (WCL) [4]. SPLIce differs crucially from those systems, including its capability to process and distribute L-volume samples and the integration autonomous control of a wide range of fluidic functions, including: 1) retrieval of fluid samples from an evacuated sample chamber; 2) onboard multi-year storage of dehydrated reagents; 3) integrated pressure, pH, and conductivity measurement; 4) filtration and retention of insoluble particles for microscopy; 5) dilution or vacuum-driven concentration of samples to accommodate instrument working ranges; 6) removal of gas bubbles from sample aliquots; 7) unidirectional flow (check valves); 8) active flow-path selection (solenoid-actuated valves); 9) metered pumping in 100 nL volume increments. The SPLIce manifold, made of three thermally fused layers of precision-machined cyclo-olefin polymer, supports all fluidic components (Figure 1) and integrated microchannels (125 x 250 m). Fluid is pumped by a stepper-motor-driven pump (Lee Co.). The functionality of the integrated MEMS pressure sensor (Honeywell) and passive check valves (Figure 2) were tested in conjunction with our newly designed integral bubble traps (Figure 3) and hydrophobic membrane-based concentrator (Figure 4). The concentrator (initially tested as a standalone component) demonstrated 5-fold vacuum-evaporative concentration. Polyethylene fused bead beds (PEFBBs; 50 porosity) store drylyophilized buffers, calibrants, and fluorescent dyes, and also promote mixing of sample with calibrant, dye, or H2O. Software-controlled automated tests demonstrated successful 1) fluid delivery to each component 2) valve and pump synchronization 3) sample aliquot delivery to instrument interface ports, and 4) rehydration of vacuum-dried fluorescent dye. In Figure 5, fluorescein on PEFBBs was rehydrated for 15 min using a pump-delivered water aliquot; it is displaced as H2O enters the bottom of the channel and pushes the dye into a check valve. Ultimately, SPLIce will fluorescently label amino acids in the sample for microchip-based electrophoretic (MCE) chiral separation and detection to seek and quantify key organic bio-signatures [5]; it will also deliver sample to a microfluidic version of WCL (mWCL) to measure soluble ions and redox-active species.

Life detection↗
NASA NTRS2020

Mineralogy of Vera Rubin Ridge in Gale Crater from the Mars Science Laboratory CheMin instrument

Gale crater was selected as the landing site for the Mars Science Laboratory Curiosity rover because of orbital evidence for a variety of secondary minerals in the lower slopes of Aeolis Mons (aka Mount Sharp) that indicate changes in aqueous conditions over time. Distinct units demonstrate orbital spectral signatures of hematite, phyllosilicate (smectite), and sulfate minerals, which suggest that ancient aqueous environments in Gale crater varied in oxidation potential, pH, and water activity. Vera Rubin ridge (VRR) is the first of these units identified from orbit to have been studied by Curiosity. Orbital near-infrared data from VRR show a strong band at 860 nm indicative of hematite. Before Curiosity arrived at VRR, the hypotheses to explain the formation of hematite included (1) precipitation at a redox interface where aqueous Fe2+ was oxidized to Fe3+, and (2) acidic alteration of olivine in oxic fluids. Studying the composition and sedimentology of the rocks on VRR allow us to test and refine these hypotheses and flesh out the depositional and diagenetic history of the ridge. Here, we focus on the mineralogical results of four rock powders drilled from and immediately below VRR as determined by CheMin.

Rampe, E. B.↗
NASA NTRS2019

In Situ Observation of Growth Dynamics in DECLIC Directional Solidification Insert Onboard ISS: DSI-R Flight Campaign

The study of solidification microstructure formation is of utmost importance for materials design and processing, as solid-liquid interface patterns largely govern mechanical and physical properties. Pattern selection occurs under dynamic conditions of growth in which the initial morphological instability evolves nonlinearly and undergoes a reorganization process. The dynamic and nonlinear nature of this instability renders in situ observation of the interface an invaluable tool to gain knowledge on the time-evolution of the interface pattern. Transparent organic analogs, which solidify like metallic alloys, allow direct visualization of interface dynamics. Extensive ground-based studies of both metallic and organic bulk samples have established the presence of significant convection during solidification processes that alters the formation of interfacial microstructures. A reduced-gravity environment is therefore mandatory for fluid flow elimination in bulk samples. In the framework of the CNES project MISOL3D (MIcrostrutures de SOLidification 3D) and the NASA projects DSIP (Dynamical Selection of 3D Interface Patterns), SPADES (SPAtiotemporal Evolution of three-dimensional DEndritic array Structures) and CAMUS (ComputAtional Studies of MicrostrUcture Formation During Alloy Solidification in Microgravity), we participated in the development of the Directional Solidification Insert (DSI) of the DEvice for the study of Critical Liquids and Crystallization (DECLIC). The DECLIC-DSI is dedicated to in situ and real-time characterization of solid-liquid interface patterns during directional solidification of transparent alloys in diffusive transport regime. Between April 2010 and March 2011, the first ISS campaign (DSI) explored the entire range of microstructures resulting in unprecedented observations. A second campaign (DSI-R), performed between October 2017 and December 2018, in which the insert contained an alloy of higher solute concentration, allowed to complete the benchmark database. The increase of solute concentration resulted in well-developed dendritic patterns at lower velocities (lower interface curvature and larger tip radius). The microstructure resulting from dendritic growth is dominant in metallurgy so that it is fundamental to understand the mechanisms of its formation. The main aims of this experimental campaign are to understand: the mechanisms of the cell to dendrite transition, the fundamental mechanisms of sidebranching formation, the dependence of dendrite tip shapes on growth conditions, the interaction of primary array and secondary sidebranches, and the influence of subgrain boundaries on the spatiotemporal organization of the array structure. Preparation, analysis and interpretation of the experiments performed onboard ISS are considerably enhanced by experiments performed on ground using thin-samples (Pr. Trivedi’s group) and phase-field simulations of microstructure formation in a diffuse growth regime (Pr. Karma’s group). In this summary, we will present an initial assessment of the results obtained during the DSI-R campaign. Nomenclature Vp: pulling velocity G: thermal gradient L: solidified length rtip: tip curvature radius Acronyms/Abbreviations CADMOS: Centre d'Aide au Développement des activités en Micro-pesanteur et des Opérations Spatiales CNES: Centre National d’Études Spatiales DECLIC: Device for the study of Critical LIquids and Crystallization DSI: Directional Solidification Insert DSI-R: Directional Solidification Insert - Refurbish ISS: International Space Station NASA: National Aeronautics and Space Administration PF: phase-field SCN: Succinonitrile 3D: three-dimensional

Mota, F.L.↗
DOE OSTI2021

Amphilic Water-Lean Carbon Capture Solvent Wetting Behavior through Decomposition by Stainless-Steel Interfaces

In this work, we report on a combined experimental and theoretical study of the wetting and reactivity of water-lean carbon capture solvents on the surface of common column packing materials. Paradoxically, we find that these solvents are equally able to wet hydrophobic and hydrophilic surfaces. The solvents are amphiphilic and are able to adapt to any interfacial environment due to their inherent heterogeneous molecular structure. Ab initio molecular dynamics indicates that these structures enable the formation of a strong adlayer on the surface of hydrophilic surfaces like oxidized steel which promotes solvent decomposition akin to hydrolysis from surface oxides and hydroxides. This decomposition passivates the surface, making it effectively hydrophobic, and the decomposed solvent promotes leaching of the iron into the bulk fluid. This study links the wetting behavior to the observed corrosion of the steels via decomposition of solvent at steel interfaces. The overall affect is strongly dependent on the chemical composition of the solvent, in that amines are stable, whereas imines and alcohols are not. We also observed that plastic packing shows little to no solvent degradation, but an equal degree of wetting.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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↗