Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density distribution”

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,189 records · Page 66

Reproducible benchmark for the SNAP 8 experimental reactor at operating conditions

This work presents fully reproducible multiphysics benchmark models of the Systems for Nuclear Auxiliary Power (SNAP) 8 Experimental Reactor at operating conditions with coolant flow. Wet experiment (with coolant, at power) validation benchmarks are presented using both deterministic (Serpent-Griffin) and Monte-Carlo (OpenMC-Cardinal) multiphysics frameworks coupled with thermal-hydraulic solvers in MOOSE. Reactivity coefficient measurements including fuel temperature, isothermal temperature, and power coefficients show good agreement with experiments, with discrepancies within experimental uncertainty. Reactivity worth experiments for coolant, samarium, and xenon poisoning are reproduced with differences under 200 pcm. Comparison between Serpent-Griffin and OpenMC-Cardinal frameworks reveal multiphysics coupling introduces positive reactivity effects (100-200 pcm) compared to uniform temperature and density fields at nominal operating conditions. Comparison between Serpent-Griffin and reference Serpent solution shows that power distributions maintain consistent radial and axial peaking behavior. All models, assumptions, thermophysical and thermomechanical properties, and material definitions are thoroughly documented with cited references; model inputs and model generating scripts are stored in the snapReactors GitHub repository.

SNAP↗

Modeling the microplastic distribution along the Delaware River Estuary: Accumulation patterns and hydrodynamic influences

Microplastic pollution is an escalating environmental concern, particularly in densely populated estuary regions, where it poses significant threats to aquatic life and human health. The dispersion patterns of microplastic particles along estuaries are influenced and complicated by multiple environmental factors such as river flow, tidal mixing, salt intrusion, and estuarine circulation. This study examines the accumulation and dispersion patterns by modeling three typical classes of microplastics in the Delaware River Estuary: synthetic fibers, sinking plastic films, and rising plastic pellets. Our findings reveal specific areas with high microplastic accumulation for each type. Notably, the upper estuary regions exhibit significant retention of rising microplastics, associated with a region with reduced along-thalweg velocities downstream of Trenton, NJ and upstream of Philadelphia, PA. Conversely, synthetic fibers and sinking plastic films accumulate in the flow convergence zone near the bottom salinity front, typically downstream of Philadelphia. All of the microplastic accumulation hot spot locations are controlled by the balance of river discharge and salinity intrusions. During the dry season, microplastic accumulation hot spots shift upstream in the estuary, whereas in the wet season, the strong river discharge pushes them downstream. Furthermore, on the other hand, tidal mixing, settling, and resuspension processes strongly impact the spreading of microplastics along the river.

Delaware River Estuary↗

Direct numerical simulations of three-component Rayleigh–Taylor mixing and an improved model for multicomponent reacting mixtures

We present direct numerical simulations of a three-layer Rayleigh–Taylor instability (RTI) problem with a configuration based on the experiments of Suchandra & Ranjan ( J. Fluid Mech. , vol. 974, 2023, A35) and Jacobs & Dalziel ( J. Fluid Mech. , vol. 542, 2005, pp. 251–279). The problem consists of a layer of light fluid between two layers of heavy fluid with an Atwood number of 0.3. These simulations are first validated through comparison with available experimental data. The validated simulations are then utilized to analyse statistics in this three-component flow. First, length scales are examined utilizing spectra and two-point spatial correlations of velocity and species concentration fluctuations. Next, joint probability density functions (p.d.f.s) of species concentration are compared against several model p.d.f.s representing generalizations of the bivariate beta distribution. Notably, the joint p.d.f.s do not appear to be accurately described by a Dirichlet distribution, indicating the marginal distributions do not conform to a beta distribution. Finally, similarity of the present configuration to three-component mixing found in inertial confinement fusion (ICF) applications is exploited to develop and validate an improved model for the impact of multicomponent mixing on thermonuclear (TN) reaction rates. A single time instant from the present simulations is chosen for a TN burn calculation under the hypothetical assumption of ICF materials and temperatures. Total TN output from this second calculation is then compared against the prediction of the improved model. The new model is found to accurately predict TN reaction rates in both premixed and non-premixed configurations.

42 ENGINEERING↗

Hydrogen Dynamics in Supercritical Water Probed by Neutron Scattering and Computer Simulations

In this work, an investigation of supercritical water is presented combining inelastic and deep inelastic neutron scattering experiments and molecular dynamics simulations based on a machine-learned potential of ab initio quality. The local hydrogen dynamics is investigated at 250 bar and in the temperature range of 553–823 K, covering the evolution from subcritical liquid to supercritical gas-like water. The evolution of libration, bending, and stretching motions in the vibrational density of states is studied, analyzing the spectral features by a mode decomposition. Moreover, the hydrogen nuclear momentum distribution is measured, and its anisotropy is probed experimentally. It is shown that hydrogen bonds survive up to the higher temperatures investigated, and we discuss our results in the framework of the coupling between intramolecular modes and intermolecular librations. Results show that the local potential affecting hydrogen becomes less anisotropic within the molecular plane in the supercritical phase, and we attribute this result to the presence of more distorted hydrogen bonds.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Task Sharing of Proton Incorporation in Vertically Aligned Nanocomposite Triple Conductors: Growth, Structure, and Surface Exchange Kinetics

As protonic ceramic electrolysis cells emerge for efficient H 2 production, there is a need to develop air electrode materials enabling fast, durable steam splitting and proton incorporation. Single-phase triple conductors may fail to satisfy the myriad performance/stability requirements, and their critical charge-carriers (holes, oxygen vacancies, and protons) are in competition, limiting their concentrations. Instead, we propose task-sharing, vertically aligned nanocomposites (VANs), comprising a proton conductor (BaZr 0.9 Y 0.1 O 3-δ ) and a redox-active mixed ionic electronic conductor (Ce 0.9 Pr 0.1 O 2-δ ), that may enable rapid proton surface exchange at the solid–gas interface and transport along the solid–solid heterointerfaces. We grew VANs by pulsed laser deposition and investigated the interplay between their processing conditions, structure, and proton and oxygen surface exchange kinetics. We varied the substrate temperature, laser repetition rate, laser fluence, and processing oxygen pressure. The crystallinity and phases were characterized by grazing-incidence X-ray diffraction, and the strain and structural order as a function of depth were evaluated by angle-dependent synchrotron X-ray pair distribution function analysis. To evaluate the potential for interdiffusion, the formation energies of substitutional defects were simulated with density functional theory. Corresponding structural analysis and elemental mapping were performed by scanning/transmission electron microscopy, energy-dispersive X-ray spectroscopy, and electron energy-loss spectroscopy, indicating distinct nanoscale compositional regions with a hierarchical structure embedded in individual VANs columns and minimal interdiffusion across a bilayer film. Proton and oxygen surface exchange coefficients (k H , k O ) and polarization resistances were evaluated by electrical and optical relaxations and impedance spectroscopy of VAN-incorporated protonic ceramic electrochemical cells, respectively, at 400–500 °C, demonstrating values comparable to some of the best-known triple and mixed conductors.

36 MATERIALS SCIENCE↗

Effect of Formation Heterogeneity on CO 2 Dissolution in Subsurface Porous Media

Dissolution trapping is one of the most dominant mechanisms for the secure storage of CO 2 injected in porous subsurface formations saturated with brine. This trapping mechanism is enhanced by convective mixing, which occurs due to the gravitational instability between the different fluid layers in the aquifer. The reservoir permeability also plays a crucial role in the dissolution rate and overall fluid flow dynamics during the density-driven convection in porous media. Herein this study investigates the role of complex heterogeneity, i.e., irregular permeability distribution in CO 2 dissolution, using a novel experimental approach to create medium permeability heterogeneity in Hele–Shaw cells. Complex subsurface transport phenomena such as a preferential dissolution path, CO 2 sweep efficiency, changes in finger morphology, and CO 2 concentration distribution are visualized by creating heterogeneous media. Experimental results showed that reservoir permeability heterogeneity causes significant channeling effects and poor sweep efficiency. A scaling relationship between average finger growth rate (Gr) and permeability (k) was obtained as Gr [m s –1 ] = 266.8k [m 2 ] + 1.20 × 10 –6 . Furthermore, the mass of CO 2 dissolved is calculated using the spectrophotometric method to characterize the convective instability. The convective flux was analyzed by comparing the experimental dissolution flux with the theoretical diffusion flux, calculating a maximum Sherwood number of 6.8. The study’s findings improve the current understanding of the CO 2 convection morphology in heterogeneous media, allowing better assessment of long-term CO 2 storage.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Quantitative Determination of Metal Ion Adsorption on Cellulose Nanocrystals Surfaces

Nanocellulose is a bio-based material that holds significant potential in the field of water purification. Of particular interest is their potential use as a key sorbent material for the removal of metal ions from solution. However, the structure of metal ions adsorbed onto cellulose surfaces is not well understood. The focus of this work is to determine quantitatively the three-dimensional distribution of metal ions of different valencies surrounding negatively charged carboxylate functionalized cellulose nanocrystals (CNCs) using anomalous small-angle X-ray scattering (ASAXS). These distributions can affect the water and ionic permeability in these materials. The data show that increasing the carboxylate density on the surface of the CNCs from 740 to 1100 mmol/kg changed the nature of the structure of the adsorbed ions from a monolayer into a multilayer structure. The monolayer was modeled as a Stern layer around the CNC nanoparticles, whereas the multilayer structure was modeled as a diffuse layer on top of the Stern layer around the nanoparticles. Within the Stern layer, the maximum ion density increases from 1680 to 4350 mmol of Rb + /(kg of CNC) with the increase in the carboxylate density on the surface of the nanoparticles. Additionally, the data show that CNCs can leverage multiple mechanisms, such as electrostatic attraction and the chaotropic effect, to adsorb ions of different valencies. By understanding the spatial organization of the adsorbed metal ions, the design of cellulose-based sorbents can be further optimized to improve the uptake capacity and selectivity in separation applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Competing electronic states emerging on polar surfaces

Excess charge on polar surfaces of ionic compounds is commonly described by the two-dimensional electron gas (2DEG) model, a homogeneous distribution of charge, spatially-confined in a few atomic layers. Here, by combining scanning probe microscopy with density functional theory calculations, we show that excess charge on the polar TaO 2 termination of KTaO 3 (001) forms more complex electronic states with different degrees of spatial and electronic localization: charge density waves (CDW) coexist with strongly-localized electron polarons and bipolarons. These surface electronic reconstructions, originating from the combined action of electron-lattice interaction and electronic correlation, are energetically more favorable than the 2DEG solution. They exhibit distinct spectroscopy signals and impact on the surface properties, as manifested by a local suppression of ferroelectric distortions.

36 MATERIALS SCIENCE↗

Synergistic effects of mixing and strain in high entropy spinel oxides for oxygen evolution reaction

Developing stable and efficient electrocatalysts is vital for boosting oxygen evolution reaction (OER) rates in sustainable hydrogen production. High-entropy oxides (HEOs) consist of five or more metal cations, providing opportunities to tune their catalytic properties toward high OER efficiency. This work combines theoretical and experimental studies to scrutinize the OER activity and stability for spinel-type HEOs. Density functional theory confirms that randomly mixed metal sites show thermodynamic stability, with intermediate adsorption energies displaying wider distributions due to mixing-induced equatorial strain in active metal-oxygen bonds. The rapid sol-flame method is employed to synthesize HEO, comprising five 3d-transition metal cations, which exhibits superior OER activity and durability under alkaline conditions, outperforming lower-entropy oxides, even with partial surface oxidations. The study highlights that the enhanced activity of HEO is primarily attributed to the mixing of multiple elements, leading to strain effects near the active site, as well as surface composition and coverage.

08 HYDROGEN↗

SIT 45: An interacting, compact, and star-forming isolated galaxy triplet

The underlying scenario of the formation and evolution of galaxy triplets is still uncertain. Mergers of galaxies in isolated triplets give us the opportunity to study the already complex merging process, with minimal contamination of other environmental effects that potentially allow and accelerate galaxy transitions from active star-forming to passive galaxies. The merging system SIT 45 (UGC 12589) is one of 315 systems in the SDSS-based catalogue of Isolated Triplets (SIT); it is an unusual isolated galaxy triplet, consisting of three merging late-type galaxies. The main aims of this work are to study the dynamical evolution and star formation history (SFH) of SIT 45, as well as its dependence on its local and large-scale environment. Methods. To study its dynamics, parameters such as the velocity dispersion (σ v ), the harmonic radius (R H ), the crossing time (H 0 t c ), and the virial mass (M vir ), along with the compactness of the triplet (S) were considered. To investigate the possible dependence of these dynamical parameters on the environment, the tidal force Q parameters (both local and large-scale) and the projected local density (η k ) were used. To constrain the SFH, we used CIGALE to fit its observed spectral energy distribution using multiwavelength data from the ultraviolet to the infrared. SIT 45 is one of the most compact triplets in the SIT, and it is also more compact than triplets in other samples. According to its SFH, SIT 45 presents star formation, where the galaxies also present recent (~200 Myr) star formation increase, indicating that this activity may have been triggered by the interaction. Its dynamical configuration suggests that the system is highly evolved in comparison to the SIT. However, this is not expected for systems composed of star-forming late-type galaxies, based on observations in compact groups. We conclude that SIT 45 is a system of three interacting galaxies that are evolving within the same dark matter halo; its compact configuration is a consequence of the ongoing interaction rather than being due to a long-term evolution (as suggested from its H 0 t c value). We consider two scenarios for the present configuration of the triplet, one where one of the members is a tidal galaxy, and another where this galaxy arrives to the system after the interaction. Both scenarios need further exploration. The isolated triplet SIT 45 is therefore an ideal system to study short timescale mechanisms (~10 8 yr), such as starbursts triggered by interactions which are more frequent at higher redshift.

79 ASTRONOMY AND ASTROPHYSICS↗

Deterministic modeling of hybrid nonlinear effects in epsilon-near-zero thin films

In nonlinear optics, significant effort is concentrated on improving the strength and efficiency of interactions; however, experimentally investigating nonlinear materials is a complex, time-consuming, and costly investment. Moreover, it is often challenging to isolate, study, and optimize material parameters in an experiment due to complexities in the growth process. Recently, epsilon-near-zero materials have received a great deal of attention as promising nonlinear optical materials, but like many up-and-coming materials, the ability to explore and optimize their properties has been challenging. Here, we establish a framework to rapidly evaluate the performance of nonlinear epsilon-near-zero materials for both inter- and intraband effects in silico, requiring only an energy-momentum (E-k) diagram, linear optical properties, and experimental conditions. Measured nonlinear reflection and transmission in gallium-doped zinc oxide films are compared to the numerical framework for both intra- and interband excitation to verify accuracy across wavelength and irradiance while two figures of merit (FoMs) are introduced to quickly evaluate the performance of films without a full numerical framework. This capability is used to predict the performance of highly doped gallium nitride, cadmium oxide, zinc oxide, and indium tin oxide films, and efficient intra- and interband operation conditions are identified. Through this numerical framework and the FoMs, the exploration of unstudied epsilon-near-zero materials is enabled without the need for a nonlinear experiment, thereby accelerating the search for more efficient nonlinear materials and excitation conditions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Axisymmetric gyrokinetic simulation of ASDEX-Upgrade scrape-off layer using a conservative implicit BGK collision operator

Collisions play an important role in turbulence and transport of fusion plasmas. For kinetic simulations, as the collisionality increases in the domain of interest, the size of the time step to resolve the collisional physics can become overly restrictive in an explicit time integration scheme, leading to high computational cost. With the aim of overcoming such restriction, we have implemented an implicit Bhatnagar–Gross–Krook (BGK) collision operator for use in the discontinuous Galerkin full-f gyrokinetic solver within the Gkeyll framework, which, when combined with Gkeyll's traditional explicit time integrator for collisionless advection, can significantly increase the time step in gyrokinetic simulations of highly collisional regimes. To ensure conservation of density, momentum, and energy, we utilize an iterative scheme to correct the discretized approximation to the equilibrium Maxwellian distribution to which the BGK collision operator relaxes. We have further generalized the BGK infrastructure, both the implicit scheme and the correction routine, to handle cross-species collisions. This improved implicit and conservative BGK operator is benchmarked against the more accurate but more computationally expensive Lenard–Bernstein–Dougherty (LBD) operator, which has been utilized in prior studies with Gkeyll. The implicit BGK operator enables 2D axisymmetric simulations of the ASDEX-Upgrade scrape-off layer to run 56 times faster to completion than the simulations with the LBD operator, because the BGK operator is more robust and converges at a lower resolution than is required by the LBD operator. Additionally, in this more collisional limit, we demonstrate that the results of our simulations utilizing the implicit BGK operator agreed well with simulations utilizing the more computationally expensive LBD operator.

Gyrokinetic simulations↗

Computational Thermal Hydraulics of a High-Performance Low-Enriched-Uranium Annular Target for HFIR Irradiation

Molybdenum-99 has historically been generated via isolation from fissioned highly enriched uranium (HEU) targets. Here, this isotope is in high demand due to its daily use across the world in radiopharmaceutical medical procedures. The primary objective of this work was to design and analyze an experimental target assembly containing one low-enriched-uranium (LEU) annular target for irradiation at the High Flux Isotope Reactor (HFIR). Efforts included incorporating spatially dependent energy sources from neutron and gamma interactions, quantifying thermal contact conductance at material interfaces, performing grid-independent studies, comparing turbulence models, and simulating various steady-state and transient scenarios relevant for irradiation qualification and eventual insertion. These models provide velocity, pressure, and temperature distributions in both space and time. Such results enable the selection of an appropriate irradiation location, fission rate density, and flow-limiting orifice size and demonstrate compliance with HFIR safety requirements such that insertion into the reactor can be approved. This analysis shows that across all scenarios, wetted surface temperatures remain below the coolant saturation temperature with no net vapor formation in the coolant. In every scenario, all components stay below 30% of the aluminum 6061 melting temperature. Computational fluid dynamics and system-level models predict peak target temperatures that agree within 4%, though the predicted axial location of the peak differs by about 10% of the heated length due to differences in flow development length. These results de-risk the irradiation of LEU (annular targets) and strengthen a domestic, HEU-independent 99 Mo supply by providing important fuel performance data to form the foundation for a robust licensing basis.

Molybdenum-99↗

Development of HELIOS/BIPR/PARCS/MCNP6 Computation Route for WWER RPV Neutron Fluence Analysis and Validation Against Ex-Vessel Detector Measurement Data

Here, a computational route was developed for precise calculation of fast neutron fluence on a WWER-type reactor pressure vessel (RPV). The method is based on the transfer of neutronics data from HELIOS-2 lattice calculations and nodal diffusion neutronics data (power, density, and temperature) from BIPR7.1 and PARCS 3.36/PATHS core calculations into a three-dimensional (pinwise axially distributed) fixed neutron source for modeling of transport of fast neutrons from the reactor core to the outer surface of the RPV using MCNP6.2. Validation of the proposed computational method was carried out based on comparative analysis of MCNP6.2-predicted and neutron dosimetry–measured reaction rates [ 54 Fe(n,p) 54 Мn, 93 Nb(n,n') 93 mNb, and 58 Ni(n,р) 58 Со] on the outer surface of the Armenian Nuclear Power Plant (ANPP) Unit 2 RPV. Validation revealed that the MCNP6.2-predicted fast neutron fluence results are very sensitive to the ENDF-B neutron data. Particularly, MCNP6.2 with ENDF/B-VIII.0 significantly underpredicts (20% to 30%) fast neutron fluence while using ENDF/B-VII.1 data overpredicts it. Adding revised beta-released evaluations of 54 Fe, 56 Fe, 57 Fe, and 16 O from the International Nuclear Data Evaluation Network (INDEN) to ENDF/B-VIII.0 allows one to obtain reasonable agreement with measurement results for all types of measured reaction rates.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Velocity reconstruction in the era of DESI and Rubin/LSST. I. Exploring spectroscopic, photometric, and hybrid samples

Peculiar velocities of galaxies and halos can be reconstructed from their spatial distribution alone. This technique is analogous to the baryon acoustic oscillations reconstruction, using the continuity equation to connect density and velocity fields. The resulting reconstructed velocities can be used to measure imprints of galaxy velocities on the cosmic microwave background like the kinematic Sunyaev-Zel’dovich effect or the moving lens effect. As the precision of these measurements increases, characterizing the performance of the velocity reconstruction becomes crucial to allow unbiased and statistically optimal inference. In this paper, we quantify the relevant performance metrics: the variance of the reconstructed velocities and their correlation coefficient with the true velocities. We show that the relevant velocities to reconstruct for kSZ and moving lens are actually the halo—rather than galaxy—velocities. We quantify the impact of redshift-space distortions, photometric redshift errors, satellite galaxy fraction, incorrect cosmological parameter assumptions and smoothing scale on the reconstruction performance. Here, we also investigate hybrid reconstruction methods, where velocities inferred from spectroscopic samples are evaluated at the positions of denser photometric samples. We find that using exclusively the photometric sample is better than performing a hybrid analysis. The 2 Gpc/ℎ length simulations from abacussummit with realistic galaxy samples for DESI and Rubin LSST allow us to perform this analysis in a controlled setting. In the companion paper [B. Hadzhiyska, S. Ferraro, B. Ried Guachalla, and E. Schaan, companion paper, Phys. Rev. D 109, 103534 (2024).], we further include the effects of evolution along the light cone and give realistic performance estimates for DESI luminous red galaxies, emission line galaxies, and Rubin LSST-like samples.

79 ASTRONOMY AND ASTROPHYSICS↗

Computing nuclear response functions with time-dependent coupled-cluster theory

We compute nuclear response functions by solving the time-dependent 𝐴-body Schrödinger equation, recording the time-dependent transition moment and extracting spectral information via Fourier transforms. The solution of the time-dependent many-body problem accounts for correlations on top of the mean field by taking advantage of a time-dependent formulation of coupled-cluster theory. As a validation, we focus on electric dipole transitions in 4 He and 16 O and compare moments of the response function distribution to the results of an equivalent static framework, finding negligible discrepancies. We investigate how proton and neutron densities evolve in time, and we see the traditional picture of soft and giant dipole resonances as collective oscillations of protons and neutrons emerging from our calculations in 16 O and 24 O. Furthermore, this method also allows us to investigate the behavior of the nucleus in the presence of a strong electric field. In that regime, the behavior of the system becomes chaotic. Qualitatively, the spectral information obtained in this limit is in line with previous time-dependent mean-field results.

Ab initio calculations↗

A Multi-Fidelity Gaussian Process Regression Method for Probabilistic Wind Farm Power Curve Estimation

Accurate estimation of the power curve for wind turbines or wind farms is crucial to ensure their efficient operation and management. However, conventional methods for power curve estimation rely either on expensive and infrequent measurements or on low-quality numerical simulations. Moreover, the majority of previous studies on power curve estimation for wind turbines or wind farms focused on deterministic estimation, which provides a point estimate of the relationship between wind speed and power generation. Nevertheless, the deterministic approach fails to consider the inherent uncertainty associated with wind energy production resulting from varying turbine characteristics. This can lead to inaccurate power generation estimation and suboptimal decisions regarding energy management. In this paper, a kernel density estimation (KDE) based Multi-Fidelity Gaussian Process Regression (MFGPR) model is proposed to fuse theoretical power curve data and the ground true measurements to create a mapping of wind speed and wind power. By conducting a case study on an actual wind farm in China, the efficacy of the proposed MFGPR model was demonstrated in characterizing the variability of wind power. The probabilistic MFGPR model was also able to generate confidence intervals that encompassed the measured power, thereby improving the accuracy and confidence in wind power estimation or wind resource assessment. Overall, the proposed MFGPR model offers a reliable approach to integrate high-fidelity ground measurements and theoretical power curve data, resulting in precise wind resource assessment and power estimation.

Gaussian process regression↗

Predictive Modeling of Electrodeposition in a Single Pore Flow-Through Electrode: From Electronucleation to Coating Thickness Uniformity

The emergence of advanced manufacturing methods capable of producing porous three-dimensional structures has expanded the design space for next-generation functional components. The ability to fabricate ordered 3D foams for use in electrocatalysis reactors has increased the need for controlled deposition of catalytic metals onto porous support materials, such as carbon. However, there is a lack of clear design guidelines for electrodeposition onto 3D substrates, due to the geometric complexity and multi-scale nature of the problem. Furthermore, electro-nucleation phenomena are often overlooked in macro-scale models of current distribution during deposition. Here, a graphite flow-through electrode (FTE) is used as a model system for copper deposition within a single pore. Potential distributions and electro-nucleation phenomena are coupled in a continuum level model by incorporating nucleation size and density as a function of overpotential, determined experimentally using in situ atomic force microscopy (AFM). The model predictions are validated by measuring the coating uniformity in the pore using micro-computed X-ray tomography (μCT). A scaling analysis comprising dimensionless parameters such as the Wagner number is presented. The simplified scaling relationship framework can guide the electrodeposition process and electrode design to optimize plating of porous substrates under fluid flow conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗