Simulation of Rough Electrodes Coupled with Structural Dynamics
Explore the source record for details and available documents.
SEARCH · Engineering Papers
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.
Explore the source record for details and available documents.
Abstract Turbine blades generally present surface roughness introduced in the manufacturing process or caused by in-service degradation, which can have a significant impact on aero-thermal performance. A better understanding of the fundamental physical mechanisms arising from the interaction between the roughness and the turbine flow at engine-relevant conditions can provide insights for the design of blades with improved efficiency and longer operational life. To this end, a high-fidelity numerical framework combining a well-validated solver for direct numerical simulation and a second-order accurate immersed boundary method is employed to predict roughness-induced aero-thermal effects on an LS89 high-pressure turbine (HPT) blade at engine-relevant conditions. Different amplitudes and distributions of surface roughness are investigated and a reference smooth-blade simulation under the same flow conditions is conducted for comparison. Roughness of increasing amplitude progressively shifts the blade suction side boundary layer transition upstream, producing larger values of the turbulent kinetic energy and higher total wake losses. The on-surface data-capturing capabilities of the numerical framework provide direct measurements of the heat flux and the skin friction coefficient, hence offering quantitative information between the surface topology and engineering-relevant performance parameters. This work may provide a benchmark for future numerical studies of turbomachinery flows with roughness.
In this paper, an experimental and numerical investigation of internal cooling channels with rib turbulators is presented with sCO 2 as the working fluid at process conditions (pressure-20.7 MPa and temperature up to 150 °C). The effect of channel aspect ratio up to 2:1 on thermal-hydraulic performance is explored in additively manufactured rectangular channels and square channels, both with and without 60 deg ribs on the top and bottom sides. The Wilson-plot method is employed to experimentally measure channel-averaged Nusselt number over a Reynolds number range up to 370,000. The friction factor is calculated from pressure drop and mass flow rate and additionally, the overall thermal performance factor (TPF) is reported. A companion computational fluid dynamics (CFD) simulation is performed for the rib turbulated cooling configurations reported in the experiments using the Reynolds average Navier–Stokes-based turbulence model. The objective of the numerical study is to gain insight into the local heat transfer augmentation in the ribbed channels as a result of varying the aspect ratio, channel configuration (square versus rectangular), operating conditions (Reynolds number) and the surface roughness, an inherent outcome of the additive manufacturing process. Surface roughness is simulated using sand grain roughness height (K S ) calculated from the experimental data, and a comparison is presented with the corresponding channel configuration with varying surface roughness heights starting from smooth surfaces (K S = 0). Experimental results indicate that the heat transfer augmentation is negligible in the rectangular channels with ribs on the long side compared to the square channel. However, it is enhanced by 60% in comparison to placing ribs on the shorter side. The TPF remains constant at around 1 for the entire range of Reynolds numbers consistent with prior work at the National Energy Technology Laboratory (NETL). Here, the simulation results highlight that increased surface roughness can have a favorable considerable influence on Nusselt number and overall thermal performance enhancement.
In a concentrated suspension, particles come into contact due to the presence of asperities on their surfaces. As a result, the contact forces and interparticle friction become one of the important factors governing the rheology of rigid particle suspensions at high concentrations. We show that a load-dependent friction model can be used to reproduce the experimentally observed shear thickening [ST—continuous and discontinuous (DST)] behavior with an increase in the asperity size. Increasing the particle surface roughness size leads to a decrease in the critical shear rate above which shear-thickening takes place, a reduction in the critical volume fraction for DST and an increase in the viscosity jump across non-shear-thickened to shear-thickened regime. In this paper, we propose a constitutive model to quantify the effect of increasing the roughness size on the rheology of dense ST suspensions as well as on the critical shear rate for ST and the critical volume fraction for DST. We fit this model to our simulation data for stress controlled shear flow of dense rough particle suspensions. Once the fitting is complete, these equations are used to predict exact volume fractions and shear stress values for transitions between three regimes on the shear stress-shear rate flow state diagram for different roughness values. Furthermore, the results of this study can be used to tune the particle surface roughness for manipulating the dense suspension rheology according to different applications.
Rough surfaces in contact with a flow of fluid exhibit alternating no-slip and free shear boundary conditions at the solid–liquid and air–liquid interfaces, respectively, thereby potentially offering drag reduction benefits. The balance between the dynamic pressure in the flow and the restoring capillary pressure in the interasperity spaces determines the stability of the Cassie state of wettability and is a function of the relative extent of no-slip and free shear regions per unit surface area. In the present study, using a fractal representation of rough surface topography, an analytical model is developed to quantify the stability of the Cassie state of wettability as well as drag reduction and the friction factor for laminar flow in a rectangular channel between nonwetting multiscale rough surfaces. A systematic study is conducted to quantify the effects of fractal parameters of the surfaces and the flow Reynolds number on drag reduction and the friction factor. The studies are used to develop friction factor curves extending the classical Moody diagram to hydrophobic and superhydrophobic surfaces. On the basis of the studies, regime maps are derived for estimating the extent of drag reduction offered by hydrophobic and superhydrophobic surfaces, revealing that superhydrophobic surfaces do not always offer the best drag reduction performance. The application of the fractal model to practical topographies of nonwetting surfaces of copper, aluminum, and zinc oxide fabricated via electrodeposition and etching is also discussed.
In this work, coaxial conductor–ceramic direct ink writing enables the printing of sensitive or encapsulated materials onto heterogeneous and rough substrates. While encasing the core fluid within a stiff ceramic shell, continuity may be maintained, even while printing onto conventionally challenging substrates. Here, we report the development of a coaxial ceramic direct ink writing suite and explore coflow interrelationships based on microfluidic principles. A coaxial nozzle is designed to facilitate the coextrusion of an alumina shell, whereas indium–tin-oxide inks constitute the core. In this manner, a core–shell ceramic element may be printed onto rough substrates for future high-temperature applications. Colloidal inks are engineered to provide the required rheological and sintering performance. Moreover, flow simulations in conjunction with microfluidic coflow principles are used to explore the coaxial printing processing space, thus controlling the core–shell architectures. Physical modeling is further used to analyze core deformations and eccentricity. As a result, simulations are validated experimentally, and the analyses are used to deposit coaxial ceramic features onto heterogeneous, high-temperature ceramic substrates.
Nanoscale device fabrication requires control over film growth at the atomic scale. Growth conditions must be tuned in consideration of interface parameters like chemical bonding, surface free energy, and lattice matching. In metals, electronic properties may also be utilized for control of physical parameters. Quantum size effects can induce metals to spontaneously form specific shapes and sizes according to their electronic structure. Unfortunately, such electronic growth is generally known only for a few systems and is typically only stable under cryogenic conditions. In this work, we explore a recently discovered class of electronic growth systems in which metal films are grown upon the relatively inert surfaces of van der Waals crystals. In this class of materials, the electronic growth is highly stable at room temperature and actually requires higher temperature annealing to achieve proper equilibrium. We work with the Au/MoS 2 system, which shows excellent stability and can readily form discrete and atomically flat nanostructures. Here, we show how the electronic growth modes facilitate the formation of atomically flat films with nanometer scale thickness. Here, the surface roughness of these films was found to be less than a single atom over several square microns, creating nearly perfect surfaces for studies of self-assembled monolayers or other applications.
Metal additive manufacturing has enabled geometrically complex internal cooling channels for turbine and heat exchanger applications, but the process gives rise to large-scale roughness whose size is comparable to the channel height (which is 500 $\mathrm {\mu }$ m). These super-rough channels pose previously unseen challenges for experimental measurements, data interpretation and roughness modelling. First, it is not clear if measurements at a particular streamwise and spanwise location still provide accurate representation of the mean (time- and plane-averaged) flow. Second, we do not know if the logarithmic layer survives. Third, it is unknown how well previously developed rough-wall models work for these large-scale roughnesses. To answer the above practical questions, we conduct direct numerical simulations of flow in additively manufactured super-rough channels. Three rough surfaces are considered, all of which are obtained from computed tomography scans of additively manufactured surfaces. The roughness’ trough to peak sizes are 0.1 $h$ , 0.3 $h$ and 0.8 $h$ , respectively, where $h$ is the intended half-channel height. Each rough surface is placed opposite a smooth wall and the other two rough surfaces, leading to six rough-wall channel configurations. Two Reynolds numbers are considered, namely $Re_\tau =180$ and $Re_\tau =395$ . We show first that measurements at one streamwise and spanwise location are insufficient due to strong mean flow inhomogeneity across the entire channel, second that the logarithmic law of the wall survives despite the mean flow inhomogeneity and third that the established roughness sheltering model remains accurate.
Describing flow resistance from the properties of an underlying surface is a challenge in surface hydrology. Runoff models must specify a resistance formulation or “roughness scheme”—describing the functional relationship between flow resistance and flow depth/velocity—and its parameters. Uncertainty in runoff predictions derives from both the selected roughness scheme (e.g., Darcy Weisbach, Manning's, or laminar flow equations), and its parameterization with a roughness coefficient (e.g., Manning's n). Both choices are informed by model calibration to data, usually discharge, and, if available, velocity. In this study, a Saint Venant Equation-based runoff model is calibrated to discharge and velocity data from 112 rainfall simulator experiments. The results are used to identify the optimal roughness scheme among four widely-used options for each experiment, and to explore whether surface properties can be used to select the optimal roughness scheme and its coefficient. Among the tested roughness schemes, a transitional flow equation provided the best fit to the plurality of experiments. The most suitable roughness scheme for a given experiment was not related to measured surface properties. Regression models predicted the calibrated roughness coefficients with adjusted values between 0.48 and 0.54, depending on the roughness scheme used. Litter cover was the best predictor of the roughness coefficient, followed by soil cover and average canopy gap size. The results suggest that selection of an optimal roughness scheme based on surface properties alone remains difficult, but that once a scheme is selected, roughness coefficients can be estimated from surface properties.
Abstract Design for cooling effectiveness in turbine blades relies on accurate models for dynamic losses and heat transfer of internal cooling passages. Metal additive manufacturing (AM) has expanded the design space for these configurations, but can give rise to large-scale roughness features. The range of roughness length scales in these systems makes morphology resolved computational fluid dynamics (CFD) impractical. However, volumetric roughness models can be leveraged, as they have computational costs orders of magnitude lower. In this work, a discrete element roughness model (DERM), based on the double-averaged Navier–Stokes equations, is presented and applied to additively manufactured rough channels, representative of gas turbine blade cooling passages. Unique to this formulation of DERM is a generalized sheltering-based treatment of drag, a two-layer model for spatially averaged Reynolds stresses, and explicit treatment of dispersion. Six different AM rough surface channel configurations are studied, with roughness trough to peak sizes ranging from 15% to 60% nominal channel passage half-width, and the roughness Reynolds number ranges from Rek = 60 to 300. DERM predictions for spatially and temporally averaged mean flow quantities are compared to previously reported direct numerical simulation results. Good agreement in the mean velocity profiles, stress balances, and drag partitions are observed. While DERM models are typically calibrated to specific deterministic roughness morphologies at comparatively small roughness Reynolds numbers, the present more generalized DERM formulation has wider applicability. Here, it is demonstrated that the model can accommodate random roughness of large scale, typical of AM.
Heat exchangers are used in many industries to transfer thermal energy between two or more objects with an associated thermal and hydraulic performance indicated by thermal effectiveness and pressure drop respectively. Many investigations have aimed to enhance the thermal effectiveness of heat exchangers by introducing turbulent flow inducers that increase the heat transfer rate at the expense of increasing pressure drop. Flow inducers like surface roughness can be controlled with advanced manufacturing methods like additive manufacturing for heat exchanger production. Selecting and manufacturing an optimal surface roughness for heat exchangers requires a more comprehensive approach that accounts for the surface roughness effect on the thermal and hydraulic performance of a sample. Since fabricating samples with additive manufacturing to a desired surface roughness can be an extensive experimental process, this study will fundamentally investigate the effect of surface roughness on the thermal and velocity boundary layers for various surface roughness profiles with physics-based computational modeling. To accurately model the surface roughness effect on thermohydraulic performance, correlations like roughness ratio, friction factor, and Nusselt number are investigated. This study expects to indicate correlations for different flow regimes which can successfully predict the effect of surface roughness on thermohydraulic performance. Results generated from modeling surface roughness effect with unique correlations are then compared to experimental data obtained from literature for validation. The results of this study will provide insight to understand the thermal fluid physics associated with surface roughness and design guidelines for heat exchangers produced with additive manufacturing.
Abstract Additive manufacturing processes, such as direct metal laser sintering (DMLS), enable the creation of novel turbine cooling internal passages and systems. However, the DMLS method produces a significant and unique surface roughness. Previous work in scaled passages analyzed pressure losses and friction factors associated with the rough surfaces, as well as investigated the velocity profiles and turbulent flow characteristics within the passage. In this study, the heat transfer characteristics of scaled additively manufactured surfaces were measured using infrared (IR) thermography. Roughness panels were CNC machined from plates of aluminum 6061 to create near isothermal roughness elements when heated. Fluid resistance differences between the aluminum roughness panels and roughness panels constructed from ABS plastic using the same roughness patterns from McClain et al. (2020, “Flow in a Simulated Turbine Blade Cooling Channel With Spatially Varying Roughness Caused by Additive Manufacturing Orientation,” ASME Turbo Expo 2020, Turbomachinery Technical Conference and Exposition, Virtual Conference, Sept. 21–25, GT2020-16069) were investigated. Finally, the overall thermal performance enhancements and friction losses were assessed through the calculation of surface averaged “global thermal performance” ratios. The global thermal performance characterizations indicate results in-line with those found for traditional commercial roughness and slightly below traditional internal passage convection enhancement methods such as swirl chambers, dimples, and ribs. The passages investigated in this study do not include compressibility effects or the long-wavelength artifacts and channel geometric deviations observed by Wildgoose et al. (2020, “Impact of Additive Manufacturing on Internal Cooling Channels with Varying Diameters and Build Directions,” ASME Turbo Expo 2020, Turbomachinery Technical Conference and Exposition, Virtual Conference, Sept. 21–25, GT2020-15049). However, the results of this study indicate that, based on the roughness augmentation alone, artificial convective cooling enhancers such as turbulators or dimples may still be required for additively manufactured turbine component cooling.
At the Savannah River Site (SRS), multiple facilities are routinely producing and working with radiological constituents which require safety basis modeling by the National Nuclear Safety Administration. The models used to determine the safety basis employ roughness length as an input parameter which describes the roughness of the ground surface which influences momentum exchange near the surface and the subsequent generation of mechanical turbulence. Because roughness length is dependent on surface characteristics, including vegetation, any changes to the surface would have a subsequent impact on the roughness length and require a new safety basis analysis to be performed. An analysis of measurements prior to and after the dear-cuts and thinning was performed and roughness lengths were calculated for each of sixteen wind sectors. Not all sectors contained dear-cuts, so this provides comparison of natural variations in roughness. Based on the results of the estimated changes, consideration will be given to how these changes may affect atmospheric dispersion modeling. It was found that clear-cutting activities generally led to a reduction in roughness. Reductions in roughness length also occurred in some of the wind sectors containing only undisturbed forest, but those in sectors with dear-cuts were found to be slightly greater in magnitude than natural variation. Only in one instance was the change in roughness length found to exceed a decrease of 10%, indicating that there is no significant safety basis to the site as a whole based on previous clear-cutting activities.
Metal electrodes with rough surfaces are frequently found to convert CO or CO 2 to hydrocarbons and oxygenates with high selectivity and at high reaction rates in comparison with their smooth counterparts. The atomic-level morphology of a rough electrode is likely one key factor responsible for its comparatively high catalytic selectivity and activity. However, few methods are capable of probing the atomic-level structure of rough metal electrodes under electrocatalytic conditions. As a result, the nuances in the atomic-level surface morphology that control the catalytic characteristics of these electrodes have remained largely unexplored. Because the C≡O stretching frequency of atop-bound CO (CO atop ) depends on the coordination of the underlying metal atom, the IR spectrum of this reaction intermediate on the copper electrode could, in principle, provide structural information about the catalytic surface during electrolysis. However, other effects, such as dynamic dipole coupling, easily obscure the dependence of the frequency on the surface morphology. Further, in the limit of low CO atop coverage, where coupling effects are small, the C≡O stretching frequencies of CO atop on Cu(111) and Cu(100) facets are virtually identical. Therefore, on the basis of the C≡O stretching frequency, it is not straightforward to distinguish between these two ubiquitous surface facets, which exhibit vastly different CO reduction activities. Herein, we show that key features of the atomic-level surface morphology of rough copper electrodes can be inferred from the potential dependence of the line shape of the C≡O stretching band of CO atop . Specifically, we compared two types of rough copper thin-film electrodes that are routinely employed in the context of surface-enhanced infrared absorption spectroscopy (SEIRAS). We found that copper films that are electrochemically deposited on Si-supported Au films (CuAu–Si) are poor catalysts for the reduction of CO to ethylene in comparison to copper films (Cu–Si) that are electrolessly deposited onto Si crystals. As quantified by differential electrochemical mass spectrometry (DEMS), the onset potential for ethylene is ~200 ± 65 mV more cathodic for CuAu–Si than that for Cu–Si. To reveal the origin of the disparate catalytic properties of Cu–Si and CuAu–Si, we probed the surfaces of the electrodes with cyclic voltammetry (CV) and SEIRAS. The CV characterization suggests that the (111) surface facet predominates on CuAu–Si, whereas the (100) facet is more common on Cu–Si. SEIRAS reveals that the line shape of the C≡O stretching of CO atop is composed of two bands that are attributable to CO atop on terrace and defect sites. The different surface structures manifest themselves in the form of starkly different potential dependences of the line shape of the C≡O stretching mode of CO atop on the two types of electrodes. With a simple Boltzmann model that considers the different adsorption energies of CO atop on terrace and defect sites, and the resulting CO atop populations on terrace and defect sites, we deduced that the observed electrode-specific potential dependence of the line shape is consistent with the presence of different predominant terrace sites on the two types of films. This strategy for assessing the atomic-level morphology is not restricted to SEIRAS but could also be applied to the C≡O stretching bands recorded with surface-enhanced Raman spectroscopy (SERS), which is suitable for probing a wide range of rough copper electrodes. Therefore, with this work, we establish the potential dependence of the C≡O stretching band of CO atop as a probe of the atomic-level surface structure of rough metal electrodes under electrochemical conditions. When it is coupled with complementary techniques, this methodology provides essential structural information for further improvement in the reaction selectivity of rough metal electrodes.
Abstract Surface roughness is a well-known consequence of additive manufacturing methods, particularly powder bed fusion processes. To properly design parts for additive manufacturing, a comprehensive understanding of the inherent roughness is necessary. While many researchers have measured different surface roughness resultant from a variety of parameters in the laser powder bed fusion process, few have succeeded in determining causal relationships due to the large number of variables at play. To assist the community in understanding the roughness in laser powder bed fusion processes, this study explored several studies from the literature to identify common trends and discrepancies amongst roughness data. Then, an experimental study was carried out to explore the influence of certain process parameters on surface roughness. Through these comparisons, certain local and global roughness trends have been identified and discussed, as well as a new framework for considering the effect of process parameters on surface roughness.
Block copolymer (BCP) directed self-assembly (DSA) is a promising route to enhance lithography resolution by multiplying nanopattern density and reducing feature roughness. Eliminating kinetically trapped self-assembly defects requires fast self-assembly. However, acceleration strategies like solvent vapor annealing or homopolymer blending broaden domain interfaces, implying a trade-off in increased feature roughness. In this work, we experimentally investigate this apparent dilemma between self-assembly kinetics and line roughness for solvent vapor-annealed thin films of a lamellar poly(styrene-block-2-vinylpyridine) (PS-b-P2VP) BCP blended with PS and P2VP homopolymers. Binary blends with PS or P2VP homopolymers and ternary blends incorporating both in equal weight fractions were solvent vapor annealed using acetone, a near-neutral solvent for PS and P2VP, followed by P2VP-selective vapor-phase infiltration with alumina (AlOx) and polymer etching. Binary blends with P2VP exhibit a modest kinetic enhancement but also higher line-edge and -width roughness due to the increased frequency of P2VP protrusions and bridge defects in the alumina line patterns. In contrast, binary blends with PS self-assemble noticeably faster, while domain asymmetry from the added PS homopolymer reduces roughness by curbing the number of alumina protrusions and bridge defects. Ternary blends maintain DSA line patterns across a wider composition window and, at higher homopolymer loadings, reduce roughness at length scales near the lamellar period, consistent with a reduced impact of intradomain compositional fluctuations. These findings provide important insights for codesigning blend compositions and process flows to achieve high-resolution, defect free patterns with minimal roughness through BCP DSA.
Turbulent flows over a large surface area (S) covered by n obstacles experience an overall drag due to the presence of the ground and the protruding obstacles into the flow. The drag partition between the roughness obstacles and the ground is analyzed using an analytical model proposed by Raupach (Boundary-Layer Meteorol 60:375-395, 1992) and is hereafter referred to as R92. The R92 is based on the premise that the wake behind an isolated roughness element can be described by a shelter area A and a shelter volume V. The individual sizes of A and V without any interference from other obstacles can be determined from scaling analysis for the spread of wakes. To upscale from an individual roughness element to n/S elements where wakes may interact, R92 adopted a background stress re-normalizing instead of reducing A or V with each element addition. This work demonstrates that R92’s approach results in a linear background stress reduction in A and V only when the ratio of n/S is small, due to a low probability of wake interactions. This probabilistic nature suggests that up-scaling from individual to multiple roughness elements can be re-formulated using stochastic averaging methods proposed here. The two approaches are shown to recover R92 under plausible conditions. An alternative scaling for the shelter volume is also proposed here using thermodynamic arguments of work and dissipation though the final outcome remains similar to R92. Comparisons between R92 and available data spanning more than two decades after R92 on blocks and vegetation-like roughness elements confirm the practical utility of R92. The agreement between R92 and this updated databases of experiments and simulations confirm the potential use of R92 in large-scale models provided that the relevant parameters accommodate certain features of the roughness element type (cube versus vegetation-like) and, to a lesser extent, their configuration throughout S. Last, a comparison between R92 and models based on first-order closure principles with constant mixing length suggests that R92 can outperform such models when evaluated across a wide range of roughness densities.