Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Shear Layers”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Plate motion in sheared granular fault system

Plate motion near the fault gouge layer, and the elastic interplay between the gouge layer and the plate under stick-slip conditions, is key to understanding the dynamics of sheared granular fault systems. Here, a two-dimensional implementation of the combined finite-discrete element method (FDEM), which merges the finite element method (FEM) and the discrete element method (DEM), is used to explicitly simulate a sheared granular gouge fault system. In this work we focus on investigating the influence of normal load, driving shear velocity and plate stiffness on the velocities and displacements in the direction parallel to the shear direction (x-direction) measured at locations on the upper and lower plates just adjacent to the gouge. The simulations show that during slip phases the magnitudes of the measured velocities on the upper and lower plates are proportional to the normal load and may be inversely proportional to the square root of the plate's shear modulus. Whereas, the driving shear velocity does not show distinct influence on the measured velocities. Additionally, large slip velocities are generally associated with large macroscopic friction coefficient drops. For the models subjected to smaller normal loads, larger shear velocities and with stiffer shear plates, the same magnitude of slip velocity could cause a larger drop of macroscopic friction coefficient. During stick phases, the velocities of the upper and lower plates are respectively slightly greater and slightly smaller than half of the driving shear velocity and are both in the same direction of shear. The shear strain rate of the gouge is calculated from this velocity difference between the upper and lower plate during stick phases and thus the gouge effective shear modulus can be calculated. The results show that the gouge effective shear modulus increases proportionally with normal load, while the influence of shear velocity and plate stiffness on gouge effective shear modulus is minor. The simulations address the dynamics of a laboratory-scale fault gouge system and may aid in revealing the complexities of earthquake frictional dynamics.

58 GEOSCIENCES↗
DOE OSTI2022

Shallow Cumulus Entrainment Dynamics in a Sheared Environment

Abstract Vertical wind shear has long been known to tilt convective towers and reduce thermal ascent rates. The purpose of this study is to better understand the physical mechanisms responsible for reduced ascent rates in shallow convection. In particular, the study focuses on cloud-edge mass flux to assess how shear impacts mass-flux profiles of both the ensemble and individual clouds of various depths. A compositing algorithm is used to distill large-eddy simulation (LES) output to focus on up- and down-shear cloud edges that are not influenced by complex cloud geometry or nearby clouds. A direct entrainment algorithm is used to estimate the mass flux through the cloud surface. We find that the dynamics on the up- and down-shear sides are fundamentally different, with the entrainment of environmental momentum and dilution of buoyancy being primarily responsible for the reduced down-shear ascent rates. Direct estimates of fluid flow through the cloud interface indicate a counter-shear organized flow pattern that entrains on the down-shear side and detrains on the up-shear side, resulting from the subcloud shear being lifted into the cloud layer by the updraft. In spite of organized regions of entrainment and detrainment, the overall net lateral mass flux remains unchanged with respect to the no shear run, with weak detrainment present throughout cloud depth.

Meteorology & Atmospheric Sciences↗
DOE OSTI2021

The diurnal cycle of entrainment and detrainment in LES of the Southern Ocean driven by observed surface fluxes and waves

Empirical rules for both entrainment and detrainment are developed from LES of the Southern Ocean boundary layer when the turbulence, stratification, and shear cannot be assumed to be in equilibrium with diurnal variability in surface flux and wave (Stokes drift) forcing. A major consequence is the failure of downgradient eddy viscosity, which becomes more serious with Stokes drift and is overcome by relating the angle between the stress and shear vectors to the orientations of Lagrangian shear to the surface and of local Eulerian shear over 5 m. Thus, the momentum flux can be parameterized as a stress magnitude and this empirical direction. In addition, the response of a deep boundary layer to sufficiently strong diurnal heating includes boundary layer collapse and the subsequent growth of a morning boundary layer, whose depth is empirically related to the time history of the forcing, as are both morning detrainment and afternoon entrainment into weak diurnal stratification. Below the boundary layer, detrainment rules give the maximum buoyancy flux and its depth, as well a specific stress direction. Another rule relates both afternoon and nighttime entrainment depth and buoyancy flux to surface layer turbulent kinetic energy production integrals. These empirical relationships are combined with rules for boundary layer transport to formulate two parameterizations; one based on eddy diffusivity and viscosity profiles and another on flux profiles of buoyancy and of stress magnitude. Evaluations against LES fluxes show the flux profiles to be more representative of the diurnal cycle, especially with Stokes drift.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Baroclinic Residual Circulation and Mass Transport Due to Internal Tides

Baroclinic (BC) tidal residual circulation due to internal tides is investigated around islands over a shallow ridge using a numerical ocean model. Internal tides enhance vertical mixing over shallow slopes, leading to horizontal density gradients that drive BC residual circulation along the main thermocline. For a strongly stratified summer case, the vertical diffusivity estimated by the Mellor and Yamada turbulence closure model exceeds 1 × 10 -2 m 2 s -1 , and the velocity of BC residual circulations reaches 0.2 m s -1 . The magnitude of BC residual circulation is larger than that of barotropic residual tidal circulation, implying that BC residual circulation due to internal tides plays an important role in coastal ocean circulation. Furthermore, BC residual flow accounts for an equal percentage (5%) of the total tidal kinetic energy as the barotropic residual flow under summer stratification conditions. Results from a coupled sediment resuspension and transport model show the growth of intermediate nepheloid layers formed by strong bottom shear stress and BC residual flow. The residual component contributes largely (23% of the total in summer) to the total suspended sediment flux. Seasonal variability is explored, with weaker winter stratification leading to reduced mixing, and thus weaker BC residual circulation and sediment flux. Lastly, the magnitude of the BC residual circulation is also shown to be proportional to the square of the tidal amplitude.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Wind Fields in Category 1–3 Tropical Cyclones Are Not Fully Represented in Wind Turbine Design Standards

Abstract Offshore wind energy deployment in the U.S. is expected to increase in the years to come, with proposed wind farm sites located in regions with high risk for tropical cyclones. Yet, the wind turbine design criteria outlined by the International Electrotechnical Commission for extreme events may not account for wind field characteristics unique to tropical cyclones. To evaluate if current design standards capture the extreme conditions of these storms, we perform idealized large‐eddy simulations of three tropical cyclones (one Category 1, one Category 2, and one Category 3 storms) using the Weather Research and Forecasting model. Wind conditions near the eyewall of Cat‐1, Cat‐2, and Cat‐3 storms can exceed current design standards for offshore wind turbines. Winds at 90 m can be faster than outlined in the design criteria for Class I and Class T turbines for 50‐year recurrence periods. Wind speed shear across the turbine rotor layer is also larger than assumed in design specifications. Moreover, vertical variations in wind direction across the turbine rotor layer are large for tropical cyclones of all intensity levels, suggesting design standards should include veer, which can amplify loads in wind turbines. The existence of hurricane wind field characteristics that are not represented in design standards does not imply that damage or failure will certainly occur. Engineering safety factors incorporated in the design of the turbine's blades and structural components may prevent damage from occurring.

17 WIND ENERGY↗
DOE OSTI2023

Near-wall model for compressible turbulent boundary layers based on an inverse velocity transformation

In this work, a near-wall model, which couples the inverse of a recently developed compressible velocity transformation (Griffin et al., Proc. Natl Acad. Sci., vol. 118, 2021, p. 34) and an algebraic temperature–velocity relation, is developed for high-speed turbulent boundary layers. As input, the model requires the mean flow state at one wall-normal height in the inner layer of the boundary layer and at the boundary-layer edge. As output, the model can predict mean temperature and velocity profiles across the entire inner layer, as well as the wall shear stress and heat flux. The model is tested in an a priori sense using a wide database of direct numerical simulation high-Mach-number turbulent channel flows, pipe flows and boundary layers (48 cases, with edge Mach numbers in the range 0.77–11, and semi-local friction Reynolds numbers in the range 170–5700). The present model is significantly more accurate than the classical ordinary differential equation (ODE) model for all cases tested. The model is deployed as a wall model for large-eddy simulations in channel flows with bulk Mach numbers in the range 0.7–4 and friction Reynolds numbers in the range 320–1800. When compared to the classical framework, in the a posteriori sense, the present method greatly improves the predicted heat flux, wall stress, and temperature and velocity profiles, especially in cases with strong heat transfer. In addition, the present model solves one ODE instead of two, and has a computational cost and implementation complexity similar to that of the commonly used ODE model.

42 ENGINEERING↗
DOE OSTI2020

Kink band and shear band localization in anisotropic perfectly plastic solids

Shear-driven strain localization has been observed in a wide variety of materials and may take the form of shear bands or kink bands. Based on observations of kink bands in plastically anisotropic metallic nanolaminates and single crystal metals, here we posit that, for the specific case of isochoric deformation, the kinematics of kink band formation are indistinguishable from those of plane strain shear band formation. The only distinction between shear bands and kink bands in these systems would then be that kink bands ‘lock up’ at a particular value of material rotation while shear bands may progress to arbitrarily high strains. In order to investigate whether strong material anisotropy is sufficient to arrest shear localization at a geometry that matches the classic kink band geometry, we model the development of a band of simple shear within an anisotropic perfectly plastic material. The resulting analytical model provides the stress state needed to maintain the kinematics of simple shear as a function of material anisotropy, deformation band orientation, and shear strain (or equivalently, material rotation). It is found that plastic anisotropy can promote either kink band or shear band formation depending on the loading orientation. When the deviatoric stress is positive parallel to the plane of anisotropy, shear localization may progress without bound and a shear band is produced. When the deviatoric stress is negative parallel to the plane of anisotropy, shear localization is arrested after a certain material rotation, resulting in a kink band. Examination of the requisite applied stress state during kink band formation provides an explanation for the experimentally-observed ‘lock up’ geometry. Solutions for the band boundary inclination angle are obtained and used to provide bounds on permissible band angles for both shear bands and kink bands.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Microstructure and mechanical properties of ultrasonic spot welding of AA7075-T6 and A380 casting aluminum alloy

Ultrasonic spot welding (USW) was employed to join AA7075-T6 and casting Al alloy A380. Under a constant welding power, the influence of the welding time duration (from 1 s to 3 s) on peak temperature, microstructure, hardness, and mechanical strength of the joints were studied. It was found that the peak temperature increased with welding time. A sound solid-state metallurgical bonding was obtained at the joint interface. The cracks occurred at the edge of the joints were eventually eliminated as the welding time increased. A layer of refined grains formed in A380 adjacent to the bonding line. Further, hardness reduction occurred in AA7075 at the joint interface and the degree of hardness reduction and the volume increased with the peak temperature. By increasing welding time, lap shear tensile tests show that the peak load increased from 5.3 kN to 6.55 kN. The corresponding interfacial shear strength ranged from 96.29 MPa to 99.72 MPa, lower than the shear strength of both base materials (331 MPa for AA7075 and 185 MPa for A380). All joint samples exhibited interfacial fracture after lap shear tensile testing, with a thin layer of A380 remnant attached to AA7075 at the center of the fractured surface.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A Study of Intermittent Turbulence in Stable Arctic Boundary Layers

Abstract Polar boundary layers are difficult to model due to the existence of intermittent turbulence within stable layers. Here we present a case study evaluation of coherent structures in a stable boundary layer observed during a series of flights with an uncrewed aircraft system (DataHawk2) on 19 October 2016 at Oliktok Point, Alaska as part of the ERASMUS (Evaluation of Routine Atmospheric Sounding Measurements using Unmanned Systems) field campaign. During a sequence of five flights over a nine-hour period, 57 profiles of atmospheric properties (0–400 m a.g.l) were collected. Turbulence was identified using derived Richardson Number, temperature structure function parameter, and turbulence kinetic energy dissipation. Throughout all flights on this strongly stable day, intermittent turbulence was observed. These turbulent layers showed well-mixed potential temperature profiles embedded within otherwise stable potential temperature profiles; often resulting in a characteristic staircase pattern. Turbulent layers ranged from 1 to 30 m deep, with most individual layers being 1–2 m deep. Vertical propagation velocities of layers in the lower atmosphere were on the order of a few cm s −1 , typical of non-convective environments. In different regions of the profile, turbulence was driven by a different balance of buoyancy and shear forces, with turbulence in the near surface environment driven by strong shear forces overcoming strong resistance to buoyancy, while turbulence in elevated layers characterized by weaker shear forces overcoming weaker resistance to buoyancy. We discuss the potential of such datasets for improving subgrid parameterizations of small-scale turbulence embedded within stable boundary layers.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Effects of Tethered Polymers on Dynamics of Nanoparticles in Unentangled Polymer Melts

Polymer-tethered nanoparticles (NPs) are commonly added to a polymer matrix to improve the material properties. Critical to the fabrication and processing of such composites is the mobility of the tethered NPs. Here, we study the motion of tethered NPs in unentangled polymer melts using molecular dynamics simulations, which offer a precise control of the grafted chain length N g and the number z of grafted chains per particle. As Ng increases, there is a crossover from particle-dominated to tethered-chain-dominated terminal diffusion of NPs with the same z. The mean squared displacement of loosely tethered NPs in the case of tethered-chain-dominated terminal diffusion exhibits two subdiffusive regimes at intermediate time scales for small z. The first one at shorter time scales arises from the dynamical coupling of the particle and matrix chains, while the one at longer time scales is due to the participation of the particle in the dynamics of the tethered chains. The friction of loosely grafted chains in unentangled melts scales linearly with the total number of monomers in the chains, as the friction of individual monomers is additive in the absence of hydrodynamic coupling. As more chains are grafted to a particle, hydrodynamic interactions between grafted chains emerge. As a result, there is a nondraining layer of hydrodynamically coupled chain segments surrounding the bare particle. Outside the nondraining layer is a free-draining layer of grafted chain segments with no hydrodynamic coupling. The boundary of the two layers is the stick surface where the shear stress due to the relative melt flow is balanced by the friction between the grafted and melt chains in the interpenetration layer. The stick surface is located further away from the bare surface of the particle with higher grafting density.

42 ENGINEERING↗
DOE OSTI2025

An Observational Evaluation of RKW Theory over the U.S. Southern Great Plains

The theory of Rotunno et al. (“RKW” theory) addresses the behavior of squall-line cold pools in vertically sheared flows. It predicts that, within a given thermodynamic environment, a balance between baroclinic vorticity generation by the cold pool and low-level environmental vertical wind shear induces an upright updraft along the gust front that maximizes the initiation of new convective cells. Although this theory has been evaluated numerically, its applicability to observed systems remains unclear and is limited by a lack of critical measurements, including high-frequency thermodynamic and wind profiles across the gust front. Herein, observations from the Atmospheric Radiation Measurement Southern Great Plains (ARM-SGP) observatory near Lamont, Oklahoma, are used to evaluate RKW theory for 10 well-observed squall lines over a 11-yr period. For this evaluation, RKW parameters including cold-pool intensity (c), low-level ambient, line-normal vertical shear (ΔV n ), subcloud and cloud-layer updraft tilts, and multiple measures of system intensity are estimated. Furthermore, the c estimates rely on thermodynamic retrievals from the Atmosphere Emitted Radiance Interferometer (AERI), which are uncertain but verify reasonably well against independent observations. As predicted by the theory, for c/ΔV n ≥ 1, c/ΔV n correlates positively with updraft tilt and negatively with system intensity, but these results are not always statistically significant and are also sensitive to the method by which ΔV n is evaluated. Specifically, ΔV n evaluations that extend above the cold-pool top yield greater consistency with RKW predictions. Also, some measures of intensity correlate more strongly with standard moist instability metrics than with RKW parameters.

Cold pools↗
DOE OSTI2022

The coexistence of the streaming instability and the vertical shear instability in protoplanetary disks

The streaming instability is a promising mechanism to induce the formation of planetesimals. Nonetheless, this process has been found in previous studies to require either a dust-to-gas surface density ratio or a dust size that is enhanced compared to observed values. Employing two-dimensional global simulations of protoplanetary disks, we show that the vertical shear instability and the streaming instability in concert can cause dust concentration that is sufficient for planetesimal formation for lower surface density ratios and smaller dust sizes than the streaming instability in isolation, and in particular under conditions that are consistent with observational constraints. This is because dust overdensities forming in pressure bumps induced by the vertical shear instability act as seeds for the streaming instability and are enhanced by it. While our two-dimensional model does not include self-gravity, we find that strong dust clumping and the formation (and dissolution) of gravitationally unstable overdensities can be robustly inferred from the evolution of the maximum or the mean dust-to-gas volume density ratio. The vertical shear instability puffs up the dust layer to an average mid-plane dust-to-gas density ratio that is significantly below unity. We therefore find that reaching a mid-plane density ratio of one is not necessary to trigger planetesimal formation via the streaming instability when it acts in unison with the vertical shear instability.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Mixing across stable density interfaces in forced stratified turbulence

Understanding how turbulence enhances irreversible scalar mixing in density-stratified fluids is a central problem in geophysical fluid dynamics. While isotropic overturning regions are commonly the focus of mixing analyses, we here investigate whether significant mixing may arise in anisotropic statically stable regions of the flow. Focusing on a single forced direct numerical simulation of stratified turbulence, we analyse spatial correlations between the vertical density gradient $\partial \rho /\partial z$ and the dissipation rates of kinetic energy $\epsilon$ and scalar variance $\chi$ , the latter quantifying scalar mixing. The domain is characterized by relatively well-mixed density layers separated by sharp stable interfaces that are correlated with high vertical shear. While static instability is most prevalent within the mixed layers, much of the scalar mixing is localized to the intervening interfaces, a phenomenon not apparent if considering local static instability or $\epsilon$ alone. While the majority of the domain is characterized by the canonical flux coefficient $\varGamma \equiv \chi /\epsilon =0.2$ , often assumed in ocean mixing parametrizations, extreme values of $\chi$ within the statically stable interfaces, associated with elevated $\varGamma$ , strongly skew the bulk statistics. Our findings suggest that current parametrizations of turbulent mixing may be biased by undersampling, such that the most common, but not necessarily the most significant, mixing events are overweighted. Having focused here on a single simulation of stratified turbulence, it is hoped that our results motivate a broader investigation into the role played by stable density interfaces in mixing, across a wider range of parameters and forcing schemes representative of ocean turbulence.

Couchman, Miles M. P. (ORCID:0000000246676829)↗
DOE OSTI2026

Alkali-Metal Interlocking of 2D V 4 O 10 Sheets Defines Discretized Interlayer Shear Relationships

Low-dimensional materials manifest structural anisotropy, quantum confinement, and tightly bound excitonic states, which make them attractive building blocks that can be assembled within three-dimensional laterally stitched heterostructures, stacked van der Waals solids, and complex moiré superlattices. Ion intercalation in the galleries between layered materials provides a means of modifying interlayer separation and coupling, but it is also known to drive the shearing of the layers. In this article, we explore the distinct ligand coordination environments afforded by vanadyl oxygens of singular [V 4 O 10 ] sheets and examine how the size, polarizability, and stoichiometry of Group I cations sandwiched between such layers determine the interlocking of the sheets in stacked structures. Based on the topochemical insertion of alkali-metal ions into the layered λ-V 2 O 5 , we identify seven types of guest ion coordination sites discretized into four distinct regimes of interlayer shear in units of half octahedral widths. The coordination preferences of intercalated cations govern how they interlock 2D [V 4 O 10 ] sheets and engender specific shear conformations. We present evidence that static and dynamic disorder in guest ion arrangement modulate the magnetic structure of the intercalated compounds based on electrostatic polarization, localization of charge and spin density, and lattice distortion. The results illustrate the use of topochemical ion insertion to modulate stacking relationships and magnetic transition characteristics.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Microstructure and mechanical properties of nanoscale Cu/(Ta50Nb25Mo25) multilayers

Nanoscale Cu/(Ta50Nb25Mo25) multilayers comprised of medium entropy alloy Ta50Nb25Mo25 and pure Cu are found to exhibit similar length-scale-dependent deformation mechanisms as metallic multilayers with elemental individual layers, i.e., confined layer slip at layer thickness of tens of nanometers and shear localization at several nanometers. The hardness of the nanoscale Cu/(Ta50Nb25Mo25) multilayers is higher as compared to the pure element FCC/BCC multilayers due to the extra solid solution strengthening contribution in the BCC layer. The synergy between the size effect and solid solution strengthening mechanism is discussed. In conclusion, this is an early report of the multi-component multilayers, which serves to provide directions for developing novel multilayer materials.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Large change of interlayer vibrational coupling with stacking in Mo 1-x W x Te 2

Stacking variations in quasi-two-dimensional materials can have an important influence on material properties, such as changing the topology of the band structure. Unfortunately, the weakness of van der Waals (vdW) interactions makes it difficult to compute the stacking dependence of properties, and even in a material as simple as graphite the stacking energetics remain unclear. Mo 1 - x W x Te 2 is a material in which three differently stacked phases are conveniently accessible by temperature changes: 1 T ' , T d * , and the reported Weyl semimetal phase T d . Additionally, the transitions proceed via layer sliding, and the corresponding interlayer shear mode (ISM) is relevant not just for the stacking energetics but also for understanding the relationship between Weyl physics and structural changes. However, the interlayer interactions of Mo 1 - x W x Te 2 are not well understood, with wide variation in computed properties. We report inelastic neutron scattering of the ISM in a Mo 0.91 W 0.09 Te 2 crystal. The ISM energies are generally consistent with the linear chain model, as expected given the weak interlayer interaction, though there are some discrepancies from predicted intensities. However, the interlayer force constants K x in the T d * and 1 T ' phases are substantially weaker than that of T d at 75(3) and 83(3)%, respectively. Considering that the relative positioning of atoms in neighboring layers is approximately the same regardless of overall stacking, our results suggest that longer-range influences, such as stacking-induced electronic band-structure changes, may be responsible for the substantial change in the interlayer vibrational coupling and thus the C 55 elastic constant. These findings should elucidate the stacking energetics of Mo 1 - x W x Te 2 and other vdW layered materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Self-similarity of scalar isosurface area density in a temporal mixing layer

Understanding the behaviour of scalar isosurfaces in a turbulent flow is of particular interest for many problems in turbulent mixing that contain sharp interfaces between regions of the flow. Common examples include combustion, where the chemical reactions occur in thin regions within the flow, and the turbulent/non-turbulent interface in shear flows, where a thin region separates the rotational, turbulent motions from the irrotational, non-turbulent background. Recent advances in computing technology allow for in-depth analysis of these interface problems that are difficult to quantify in a laboratory setting. In this paper, the results of a direct numerical simulation of a passive scalar $\varPhi$ evolving on a turbulent, temporally developing mixing layer are described. A novel approach has been taken to calculate the surface area of individual scalar isosurfaces, $A_{iso}$ , throughout the simulation, as well as the mean isosurface area density, $\varSigma$ , as a function of the cross-stream direction and time. A notable finding is that the profiles of $\varSigma$ develop in a self-similar manner when scaled by the Taylor scale of the scalar field, $\lambda _\phi$ . Remarkably, the scaling appears to hold for a wide range of isovalues. A rough scaling argument based on the formal definition of $\varSigma$ and properties of a temporal mixing layer is presented which also exposes a dependence on $\lambda _\phi$ . Based on these results, a possible scaling for the isosurface area is presented as $A_{iso}/A_0 \sim (Re \, Sc)^{1/2}$ , where $Re$ and $Sc$ are local Reynolds and Schmidt numbers, respectively.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Toward Optimizing Electrospun Nanofiber Fuel Cell Catalyst Layers: Polymer–Particle Interactions and Spinnability

We investigate the effect of the poly(acrylic acid) (PAA) carrier polymer concentration on the microstructure and rheological properties of catalyst inks for electrospun polymer–electrolyte membrane fuel-cell catalyst layers. Characterization of an ink microstructure using oscillatory shear rheology showed that the catalyst particles (platinum on carbon) are significantly agglomerated in the absence of PAA or an ionomer. Both the ionomer and PAA promoted the stability of the particles against agglomeration via electrosteric stabilization by adsorbing onto the particle surface. Increasing the PAA concentration increased the stability of the particles (or reduced the agglomerated structure) due to increasing PAA coverage onto the free surface area of the particles. However, beyond a certain increase in concentration, PAA was found to predominantly remain as an excess free polymer in the ink due to an insufficient free/available surface area on the particles for further PAA coverage. Extensional rheology measurements demonstrated that PAA enhances the extensional viscosities of the inks. Consequently, increasing the PAA concentration in the ink promoted the evolution of uniform nanofibers. However, beyond a certain concentration, a significant increase in the shear viscosities of the inks led to defective fiber morphologies because of the onset of flow instabilities. Electrochemical performance comparisons between catalyst layers with different PAA concentrations showed maximum performance at the PAA concentration that led to the least agglomerated structure of the catalyst, most uniform fiber morphologies, and low concentrations of free (non-adsorbing) PAA in the electrode. These results provide a rationale for optimization of electrospun catalyst nanofibers for both spinnability and electrochemical performance.

25 ENERGY STORAGE↗