Search NASA⌕ Search

SEARCH · Search NASA

Results for “Surface Forces”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Overlooked cooling effects of albedo in terrestrial ecosystems

Radiative forcing (RF) resulting from changes in surface albedo is increasingly recognized as a significant driver of global climate change but has not been adequately estimated, including by Intergovernmental Panel on Climate Change (IPCC) assessment reports, compared with other warming agents. Here, we first present the physical foundation for modeling albedo-induced RF and the consequent global warming impact (GWI Δα ). We then highlight the shortcomings of available current databases and methodologies for calculating GWI Δα at multiple temporal scales. There is a clear lack of comprehensive in situ measurements of albedo due to sparse geographic coverage of ground-based stations, whereas estimates from satellites suffer from biases due to the limited frequency of image collection, and estimates from earth system models (ESMs) suffer from very coarse spatial resolution land cover maps and associated albedo values in pre-determined lookup tables. Field measurements of albedo show large differences by ecosystem type and large diurnal and seasonal changes. As indicated from our findings in southwest Michigan, GWI Δα is substantial, exceeding the RF Δα values of IPCC reports. Inclusion of GWI Δα to landowners and carbon credit markets for specific management practices are needed in future policies. We further identify four pressing research priorities: developing a comprehensive albedo database, pinpointing accurate reference sites within managed landscapes, refining algorithms for remote sensing of albedo by integrating geostationary and other orbital satellites, and integrating the GWI Δα component into future ESMs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Tuning Surface Adhesion Using Grayscale Electron-beam Lithography

Surface texturing of manufactured products tailors their properties, such as friction, adhesion, biocompatibility, or fluid interactions. However, advancements in this area are largely the result of trial-and-effort testing and generally lack a science-guided framework for determining the surface topography that will optimize performance. The present investigation explores grayscale electron-beam lithography as a means to create multiscale surface patterns to control surface performance. Here, we created and characterized a set of surface textures on a silicon wafer; the textures were superpositions of sine waves of varying wavelengths and amplitudes. First, the multiscale topography of the patterned surface was characterized, using profilometry and atomic force microscopy, to understand its fidelity to the designed-in pattern. The results of this analysis demonstrated how grayscale lithography accurately controlled the lateral size of features but was less precise on the vertical height of the surface, and also introduced inherent roughness below the scale of patterning. Second, a micromechanical tester was used to characterize the adhesion of the surfaces with large-scale polished silicon spheres. The results showed that adhesion could be tailored, with significant contribution from all of the designed-in length scales of topography. The strength of adhesion did not correlate with conventional roughness parameters but could be accurately modeled using simple numerical integration. Taken together, this investigation demonstrates the promise and challenges of grayscale e-beam lithography with multiscale patterns as a method for the tailoring of surface performance.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

On the additivity of climate responses to the volcanic and solar forcing in the early 19th century

The early 19th century was the coldest period over the past 500 years, when strong tropical volcanic events and a solar minimum coincided. The 1809 unidentified eruption and the 1815 Tambora eruption happened consecutively during the Dalton minimum of solar irradiance; however, the relative role of the two forcing (volcano and solar) agents is still unclear. In this study, we examine the responses from a set of early 19th century simulations with combined and separated volcanic and solar forcing agents, as suggested in the protocol for the past1000 experiment of the Paleoclimate Modelling Intercomparison Project – Phase 4 (PMIP4). From 20-member ensemble simulations with the Max Planck Institute Earth system model (MPI-ESM1.2-LR), we find that the volcano- and solar-induced surface cooling is additive in the global mean/large scale, regardless of combining or separating the forcing agents. The two solar reconstructions (SATIRE (Spectral and Total Irradiance REconstruction-Millennia model) and PMOD (Physikalisch-Meteorologisches Observatorium Davos)) contribute to a cooling before and after 1815 of ~0.05 and ~0.15 K monthly average near-surface air cooling, respectively, indicating a limited solar contribution to the early 19th century cold period. The volcanic events provide the main cooling contributions, inducing a surface cooling that peaks at ~0.82 K for the 1809 event and ~1.35 K for Tambora. After the Tambora eruption, the temperature in most regions increases toward climatology largely within 5 years, along with the reduction of volcanic forcing. In the northern extratropical oceans, the temperature increases slowly at a constant rate until 1830, which is related to the reduction of seasonality and the concurrent changes in Arctic sea-ice extent. The albedo feedback of Arctic sea ice is found to be the main contributor to the Arctic amplification of the cooling signal. Several non-additive responses to solar and volcanic forcing happen on regional scales. In the atmosphere, the stratospheric polar vortex tends to strengthen when combining both volcano and solar forcing, even though the two forcing agents separately induce opposite-sign changes in stratospheric temperatures and zonal winds. In the ocean, when combining the two forcings, additional surface cold water propagates to the northern extratropics from the additional solar cooling in the tropics, which results in regional cooling along the propagation. Overall, this study not only quantifies the surface responses from combinations of the volcano and solar forcing, but also highlights the components that cannot be simply added from the responses of the individual forcing agents, indicating that a relatively small forcing agent (such as solar in early 19th century) can impact the response from the large forcing (such as the 1815 Tambora eruption) when considering regional climates.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Validated Reactive Force Field Quantifies MXene Interfacial Properties, Mechanics, and Thermal Transport

MXenes combine rich surface chemistry, mechanical strength, and high conductivity for a multitude of emerging applications. Predictive modeling supports accelerated materials designs and has been limited by the absence of validated and transferable force fields. Here, we introduce an interpretable, reactive INTERFACE force field (IFF and IFF-R) for Ti 3 C 2 T x MXenes that is trained based on chemical knowledge and achieves quantitative agreement with experiments across lattice parameters (<0.5%), density (<0.2%), liquid contact angles, Raman spectra, and the in-plane elastic modulus (∼320 GPa). The models cover surface terminations from hydroxyl (−OH) to fluorine (−F) groups and are extensible to other chemistries. We introduce pH-resolved surface chemistry and identify dopamine adsorption mechanisms at MXene–aqueous interfaces supported by QCM-D and UV–Vis experiments. The data reveal coplanar and perpendicular binding modes and concentration-dependent multilayer assembly. We predict previously inaccessible properties, including termination-dependent cleavage energies, interlayer shear moduli and dynamic shear failure, nanoindentation and brittle fracture, anisotropic in-plane and out-of-plane thermal conductivities, including the role of defects. Agreement with available experimental data is consistently close and exceeds DFT accuracy across the benchmark properties examined. The IFF/IFF-R model is compatible with CHARMM, AMBER, OPLS, and CVFF force fields for simulations of MXenes with diverse surface terminations, electrolyte interfaces, biointerfaces, and polymer composites without additional parameters. Parameter sets, 3D models, and analysis scripts are provided for community use. The validated, reactive, and transferable IFF framework facilitates predictive design of MXene-based films, membranes, sensing interfaces, and composites.

MXene↗
DOE OSTI2024

Formation of field-induced breakdown precursors on metallic electrode surfaces

Understanding the underlying factors responsible for higher-than-anticipated local field enhancements required to trigger vacuum breakdown on pristine metal surfaces is crucial for the development of devices capable of withstanding intense operational fields. In this study, we investigate the behavior of nominally flat copper electrode surfaces exposed to electric fields of hundreds of MV/m. Our novel approach considers curvature-driven diffusion processes to elucidate the formation of sharp breakdown precursors. To do so, we develop a mesoscale finite element model that accounts for driving forces arising from both electrostatic and surface-tension-induced contributions to the free energy. Our findings reveal a dual influence: surface tension tends to mitigate local curvature, while the electric field drives mass transport toward regions of high local field density. This phenomenon can trigger the growth of sharper protrusions, ultimately leading to a rapid enhancement of local fields and, consequently, to a runaway growth instability. We delineate supercritical and subcritical regimes across a range of initial surface roughness. Our numerical results are in qualitative agreement with experimentally reported data, indicating the potential practical relevance of field-driven diffusion in the formation of breakdown precursors.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Interannual variability of spring and summer monsoon growing season carbon exchange at a semiarid savanna over nearly two decades

Eddy covariance measurements of land-atmosphere energy, carbon, and water exchange now span multiple decades at some sites, supporting an improved understanding of flux interannual variability (IAV) and its ecophysiological and physical controls. Most eddy covariance IAV studies have focused on temperate forest ecosystems, where carbon fluxes are large and flux records are longest – but also where IAV is much lower than in dryland regions, which have been identified as an essential driver of the trend and variability in the global terrestrial carbon sink. In this study, we leveraged 19 years of continuous micrometeorological measurements at the AmeriFlux US-SRM mesquite savanna site in southern Arizona, USA to quantify the IAV, trends, and drivers of carbon fluxes during the distinct spring and summer growing seasons. We also assessed the ability of modern satellite and land surface models to capture the IAV of seasonal water and carbon fluxes. Annual net ecosystem production (NEP) was small and highly variable (23 +/- 64 gC m –2 yr –1 ). Precipitation and associated measures of water availability determined most of the variability in NEP, largely through their influence on annual and seasonal gross ecosystem productivity (GEP) as opposed to ecosystem respiration (ER). Root-zone soil moisture captured between 73% (spring) and 85% (summer) of GEP variability and between 73% (spring) and 58% (summer) of ER variability. Throughout the study period, soil moisture and greenness increased with associated increases in GEP, ER and NEP. These trends were strongly influenced by very productive and wet summer growing seasons during the last two years, which were characterized by abundant understory grass cover. Typically, less than half of the variability in growing season GEP and evapotranspiration was captured by satellite-based estimates and land surface model simulations with local site forcing and calibration, highlighting the ongoing utility of long-term datasets to support careful model testing and improvement.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Wire-coil insert optimization for high-heat-load/flux synchrotron components

Many synchrotron components require high levels of internal-flow, forced-convection heat transfer to minimize surface temperatures, thermal gradients, and thermally induced stress on high-power beam-interacting surfaces. Wire-coil inserts, physically similar to a common spring, are mechanically fitted inside of component cooling passages to optimize heat transfer performance. They are routinely used in Advanced Photon Source (APS) front end and beamline high-heat-load/flux components to significantly enhance convection heat transfer-up to 400% compared to plain open passages. This has the additional benefit of greatly reducing coolant flow requirements for these components. Using several cooling passage sizes, five different wire sizes, and a range of pitch values, an experimental investigation conducted at the APS has determined the average heat transfer coefficient and resulting pressure loss as a function of water flow rate for 65 different wire-coil inserts. Data from this study have been non-dimensionalized and generalized to yield relationships that can be used to determine the heat transfer performance and resulting pressure loss for any given wire-coil insert that may be used at the APS. Through data reduction, the wire-coil insert characteristic dimensions have also been optimized to yield the highest heat transfer enhancement while minimizing the coolant flow requirements. These generalized expressions for wirecoil inserts will be presented, and they can be used by scientists and engineers during the component design process to evaluate achievable heat transfer performance and associated pressure loss, aiding in the establishment of optimized operating parameters and cooling passage flow distribution schemes.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Modeling Ni redistribution in the hydrogen electrode of solid oxide cells through Ni(OH) 2 diffusion and Ni-YSZ wettability change

Ni redistribution in the hydrogen electrodes of solid oxide cells is an important degradation mechanism. Its driving force is still under debate. This work focuses on the interplay between Ni(OH) 2 diffusion and Ni-YSZ wettability change. With a Ni(OH) 2 diffusion model developed in our previous work, we further employ three models for the Ni-YSZ contact angle to quantify the Ni-YSZ wettability change. The microstructure evolutions in a reconstructed Ni-YSZ electrode are simulated under selected experimental operating conditions. It is shown that the phenomenological model can capture the driving force of Ni spreading/detachment on YSZ surface and the Ni migration reported in experiments can be qualitatively reproduced through the competition between Ni(OH) 2 diffusion and Ni-YSZ wettability change. It is also found that the initial microstructure of the active layer and the microstructure characteristics in the support layer can strongly affect the distribution of steam partial pressure and overpotential, hence the Ni redistribution, in the active layer, which may explain the inconsistencies in experiments. Further, our results show that there are two missing pieces in current theory of Ni redistribution: a mechanism of fast Ni diffusion under humid atmosphere and the physical reason behind the Ni spreading/detachment.

25 ENERGY STORAGE↗
DOE OSTI2023

Towards a predictive, physics-based friction model for the dynamics of jointed structures

Bolted connections are ubiquitous in mechanical designs and pose a significant challenge to understanding and predicting the vibration response of assembled structures. The present paper develops a physics-based rough contact model of the frictional interactions within a joint. This model sums over the probable interactions of asperities – defined as locally maximum surface features – to determine the contact forces. Here, the tangential contact forces vary smoothly between sticking and slipping and allow the model to better capture the qualitative trends of experimental amplitude dependent frequency and damping than previous studies. Furthermore, the novel model is generalized to allow for arbitrarily varying normal pressure including potential separation to better represent the interfacial dynamics. This includes developing a new, computationally tractable approximation to the analytical Mindlin partial slip solution for tangential loading of contacting spheres. Further, the results highlight the importance of accurately characterizing the as-built topology of the interface, the plastic behavior of the contacting asperities, the relevant length scale of asperities, and the eccentricity of asperities. A predictive friction coefficient based on plasticity provides a poor match to experiments, so fitting the friction coefficient is also considered. Numerical results are compared to experiments on the Brake-Reuß Beam to assess the predictive potential of the models. While blind predictions over-predict the slip limit, the current model presents a significant improvement in physics-based modeling and highlights areas for ongoing research.

42 ENGINEERING↗
DOE OSTI2020

Influence of Terminal Carboxyl Groups on the Structure and Reactivity of Functionalized m-Carboranethiolate Self-Assembled Monolayers

The structure and function of self-assembled monolayers (SAMs) at the nanoscale are determined by the steric and electronic effects of their building blocks. Carboranethiol molecules form pristine monolayers that provide tunable two-dimensional systems to probe lateral and interfacial interactions. Additional ω-functionality, such as carboxyl groups, can be introduced to change the properties of the exposed surfaces. Here, two geometrically similar isomeric m-carborane analogues of m-mercaptobenzoic acid, 1-COOH-7-SH-1,7-C 2 B 10 H 10 and racem-1-COOH-9-SH-1,7-C 2 B 10 H 10 , are characterized and their SAMs on Au{111} are examined. The latter isomer belongs to the rare group of chiral cage molecules and becomes, to our knowledge, the first example assembled on Au{111}. Although different in symmetry, molecules of both isomers assemble into similar hexagonal surface patterns. The nearest-neighbor spacing of 8.4 ± 0.4 Å is larger than that of non-carboxylated isomers, consistent with the increased steric demands of the carboxyl groups. Computational modeling reproduced this spacing and suggests a tilt relative to the surface normal. However, tilt domains are not observed experimentally, suggesting the presence of strong lateral interactions. Analyses of the influence of the functional groups through the pseudo-aromatic m-carborane skeleton showed that the thiol group attached to either carbon or boron atoms increases the carboxyl group acidity in solution. In contrast, the acidity of the exposed carboxyl group in the SAMs decreases upon surface attachment; computational analyses suggest that the driving force of this shift is the dielectric of the environment in the monolayer as a result of confined intermolecular interactions, proximity to the Au surface, and partial desolvation.

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

Correlated Response of Protein Side-Chain Fluctuations and Conformational Entropy to Ligand Binding

The heterogeneous fast side-chain dynamics of proteins plays crucial roles in molecular recognition and binding. Site-specific NMR experiments quantify these motions by measuring the model-free order parameter (O axis 2 ) on a scale of 0 (most flexible) to 1 (least flexible) for each methyl-containing residue of proteins. In this work, we have examined ligand-induced variations in the fast side-chain dynamics and conformational entropy of calmodulin (CaM) using five different CaM–peptide complexes. O axis 2 of CaM in the ligand-free (O axis,U 2 ) and ligand-bound (O axis,B 2 ) states are calculated from molecular dynamics trajectories and conformational energy surfaces obtained using the adaptive biasing force (ABF) method. ΔO axis 2 = O axis,B 2 – O axis,U 2 follows a Gaussian-like unimodal distribution whose second moment is a potential indicator of the binding affinity of these complexes. The probability for the binding-induced O axis,U 2 → O axis,B 2 transition decreases with increasing magnitude of ΔO axis 2 , indicating that large flexibility changes are improbable for side chains of CaM after ligand binding. A linear correlation established between ΔOaxis2 and the conformational entropy change of the protein makes possible the determination of the conformational entropy of binding of protein–ligand complexes. The results not only underscore the functional importance of fast side-chain fluctuations but also highlight key motional and thermodynamic correlates of protein–ligand binding.

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

Observed and Projected Changes of Large‐Scale Environments Conducive to Spring MCS Initiation Over the US Great Plains

Abstract Mesoscale convective systems (MCSs) are frequent over the US Great Plains during spring. The link between large‐scale environments and spring MCS initiation were well established. Here, historical and future changes of spring large‐scale environments favorable for MCS initiation are investigated using an MCS tracking data set, ERA5 reanalysis, and 20 Coupled Model Intercomparison Project Phase 6 (CMIP6) models. The frequency of Great Plains low‐level jet (GPLLJ)‐related MCS environments is found to have increased by ∼41% from 1979 to 2019, consistent with the enhanced GPLLJ and more frequent MCSs. Comparing CMIP6 AMIP and historical experiments, we find that the observed GPLLJ strengthening and more frequent MCS environments are mainly due to the decadal sea‐surface temperature variations rather than external forcings. Under a high emission scenario, the frequency of GPLLJ‐related environments favorable for MCS initiation will increase by ∼65% during 2015–2100, along with a stronger GPLLJ, suggesting more frequent MCSs over the US Great Plains in a warming world.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Reduction of Pacific Double‐ITCZ Bias by Convection Parameterization in NCAR CESM2.2

Abstract The impact of convective closure on the double‐ITCZ bias in the NCAR CESM2.2 is investigated in this study. The standard CESM2.2 simulates a remarkable double‐ITCZ bias in the central and eastern Pacific, especially in boreal winter and spring. Modifications to the closure in convection parameterization scheme greatly reduce the double‐ITCZ bias in all seasons, demonstrating that convection parameterization can substantially influence the double‐ITCZ bias in CESM2.2. Further analyses suggest that convection parameterization can modulate the tropical atmosphere‐ocean feedback processes, through which it influences the SST in the southern ITCZ region and hence the double‐ITCZ bias. The changes in the upper ocean temperature advection induced by modified convective closure plays important roles in reducing the warm SST bias and double‐ITCZ precipitation bias in the southern ITCZ region. The modified convective closure improves the low‐level cloud and shortwave cloud radiative forcing in the southeastern Pacific. However, surface heat flux plays only a limited role in reducing warm SST bias and double ITCZ bias because the impacts of shortwave radiation changes are largely canceled by changes in longwave radiation and latent heat flux.

Song, Xiaoliang [Scripps Institution of Oceanograp↗
DOE OSTI2021

Enhanced resistance to decay of imprinted nanopatterns in thin films by bare nanoparticles compared to polymer-grafted nanoparticles

We extend a previous study on the influence of nanoparticles on the decay of nanoimprinted polymer film patterns to compare the effects of “bare” silica (SiO 2 ) nanoparticles and SiO 2 nanoparticles with grafted polymer layers having the same chemical composition as the polymer matrix. This method involves nanoimprinting substrate-supported polymer films using a pattern replicated from a digital versatile disc (DVD), and then annealing the patterned polymer nanocomposite films at elevated temperatures to follow the decay of the topographic surface pattern with time by atomic force microscopy imaging after quenching. We quantified the relaxation of the pattern height (“slumping”) and determined the relaxation time τ for this pattern decay process as a function of nanoparticle filler type and concentration to determine how nanoparticle additives influence relative film stability. Attractive interactions between the bare nanoparticles and the polymer matrix significantly enhance the thermal resilience of the nanopatterns to decay, compared to those of the particle brushes, wherein the particle core interactions are screened from the matrix via the brush layer. A novel aspect of this method is that it readily lends itself to in situ film relaxation measurements in a manufacturing context. We observe that the relaxation time of the pattern relaxation exhibits entropy–enthalpy compensation in the free energy parameters governing the pattern relaxation process as a function of temperature, irrespective of the NP system used, consistent with our previous experimental and computational studies.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Collision integrals within the Chapman–Enskog theory for a generalized Lennard-Jones potential

Here we report the values of the collision integrals, needed for the calculation of the macroscopic transport properties such as viscosity (η) and diffusion coefficient (D) of gases within the Chapman–Enskog kinetic gas theory, for a generalized Lennard-Jones potential (gLJ), a more general potential with an adjustable long range 1/r dependence that can describe a wide range of intermolecular interactions.

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

The myosin interacting-heads motif present in live tarantula muscle explains tetanic and posttetanic phosphorylation mechanisms

Striated muscle contraction involves sliding of actin thin filaments along myosin thick filaments, controlled by calcium through thin filament activation. In relaxed muscle, the two heads of myosin interact with each other on the filament surface to form the interacting-heads motif (IHM). A key question is how both heads are released from the surface to approach actin and produce force. Here we used time-resolved synchrotron X-ray diffraction to study tarantula muscle before and after tetani. The patterns showed that the IHM is present in live relaxed muscle. Tetanic contraction produced only a very small backbone elongation, implying that mechanosensing—proposed in vertebrate muscle—is not of primary importance in tarantula. Rather, thick filament activation results from increases in myosin phosphorylation that release a fraction of heads to produce force, with the remainder staying in the ordered IHM configuration. After the tetanus, the released heads slowly recover toward the resting, helically ordered state. During this time the released heads remain close to actin and can quickly rebind, enhancing the force produced by posttetanic twitches, structurally explaining posttetanic potentiation. Taken together, these results suggest that, in addition to stretch activation in insects, two other mechanisms for thick filament activation have evolved to disrupt the interactions that establish the relaxed helices of IHMs: one in invertebrates, by either regulatory light-chain phosphorylation (as in arthropods) or Ca 2+ -binding (in mollusks, lacking phosphorylation), and another in vertebrates, by mechanosensing.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

A computational analysis of effective R-values of buried ducts – the dynamic performance of buried ducts

Here, this paper evaluates the thermal performance of ducts partially or fully buried in loose-fill attic insulation. The overall thermal resistance between the ducts and the attic is referred to as an effective R-value. This paper shows a strong dependency of assumed attic temperature on the effective R-value. Based on the results, the effective R-value can be about twice as much with an attic temperature of 130 °F [54.4 °C], compared to when the attic temperature is 80 °F [26.7 °C]. Thus, this paper provides a polynomial regression equation based on a large set of simulations to determine the effective R-value of buried ducts depending on attic temperature and whether the HVAC system runs in cooling or heating mode. Further, the work presented in this paper investigated the potential impact of convective airflow within the attic insulation, particularly around the ducts. The analysis was based on computational fluid dynamics (CFD) and indicated that convectional forces are presented around the exterior surface of the ducts, but with negligible impact on the overall heat balance between the duct and the attic space.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Rheology of concentrated fiber suspensions with a load-dependent friction coefficient

Here we investigate the effects of fiber aspect ratio, roughness, flexibility, and volume fraction on the rheology of concentrated suspensions in a steady shear flow using direct numerical simulations. We model the fibers as inextensible continuous flexible slender bodies with the Euler-Bernoulli beam equation governing their dynamics suspended in an incompressible Newtonian fluid. The fiber dynamics and fluid flow coupling is achieved using the immersed boundary method. In addition, the fiber surface roughness might lead to interfiber contacts, resulting in normal and tangential forces between the fibers, which follow Coulomb’s law of friction. The surface roughness is modeled as hemispherical protrusions on the fiber surfaces. Their deformation results in a normal load-dependent friction coefficient. Our simulations accurately predict the experimentally observed shear thinning in fiber suspensions. Furthermore, we find that the suspension viscosity η increases with increasing the volume fraction, roughness, fiber rigidity, and aspect ratio. The increase in η is the macroscopic manifestation of a similar increase in the microscopic contact contribution to the total stress with these parameters. In addition, we observe positive and negative first N 1 and second N 2 normal stress differences, respectively, with |N 2 | < |N 1 |, in agreement with previous experiments. Last, we propose a modified Maron-Pierce law to quantify the reduction in the jamming volume fraction by increasing the fiber aspect ratio and roughness. Our results and analysis establish the use of fiber surface tribology to tune the suspension flow behavior.

42 ENGINEERING↗