Shear Stress Model Research Articles

Wall-Modeled Large-Eddy Simulation Method for Unstructured-Grid Navier–Stokes Solvers

This paper reports on the implementation and assessment of a wall-modeled large-eddy simulation (WMLES) methodology in an unstructured-grid, node-centered flow solver, FUN3D, that is developed and supported at the NASA Langley Research Center. Finite-volume (FV) and finite-element (FE) discretization schemes considered in the study provide formal second-order spatial accuracy. Large-eddy simulations (LES) resolve large-scale turbulent-flow features, and small-scale effects are modeled using the Vreman subgrid-scale model. At solid-wall boundaries, a shear-stress model is employed to provide a proper boundary-flux closure. The nonlinear equations are integrated in time using either an optimized backward difference formula or an implicit multistage Runge–Kutta temporal scheme. The implicit equations at each time step are solved by strong nonlinear iteration schemes. WMLES demonstrations are shown for two high-lift configurations, namely, the McDonnell Douglas 30P30N multi-element airfoil and a NASA High-Lift Common Research Model. Results show that the WMLES approaches implemented in the FV and FE discretization methods produce consistent solutions and are capable of capturing key aerodynamic characteristics and flow structures for high-lift configurations at a wide range of angles of attack, including maximum-lift conditions. In the 30P30N example, correct trends in the variations of integrated aerodynamic forces and moments, surface pressure distributions, and boundary-layer profiles are captured as the Reynolds number is increased.

A 44-year balanced fertilizer application affected rill erosion resistance by changing humus, aggregates, and polyvalent cation

Fertilizer application can affect the physicochemical properties of soil, such as the contents of large aggregates, humus, and exchangeable cations, thereby influencing soil erosion resistance. However, the rill erosion resistance and its key influencing factors of soil following long-term balanced fertilizer application remains unclear. This study aimed to analyze the change in rill erosion resistance following 44 years of balanced fertilizer application (No changes in the type of fertilizer) to non-calcareous soils and establish a Partial Least Squares Regression (PLSR) model of rill erodibility (kd) and soil critical shear stress (τc) to changes in the physical and chemical properties of the soil. Five treatments were designed: (1) CK (no fertilizer applied), (2) N (nitrogen), (3) NP (nitrogen plus phosphorus), (4) NK (nitrogen plus potassium), and (5) NPK (nitrogen, phosphorus, and potassium). Compared to CK, the N, NP, NK, and NPK application significantly increased τc by 49.6, 96.7, 73.6, and 36.2 %, respectively. Whereas, kd increased significantly only in the NPK treatment group. The optimal partial least squares regression model showed that mean weight diameter (MWD) had the greatest positive influence on soil critical shear stress, followed by fulvic acid (FA) content, whereas water-dissolved substances had a negative influence. Long-term balanced fertilizer application can increase MWD and τc by combining micro-aggregates and humus into large aggregates. Ca2+ content had the greatest positive effect on kd. Compared with that of CK, exchangeable Ca2+ content increased significantly with NP and NPK application (7.9 and 34.3 %, respectively). Ca2+ can increase kd by binding to the polar functional groups in FA to promote the shedding of hydration shells in large aggregates. Among all treatments, the NP treatment showed the best performance for reducing kd and increasing τc. This study could contribute to the understanding of the rill erosion process and modeling in non-calcareous soils and offer a reference for agricultural erosion control treatments.

Computational fluid dynamics simulating of the FDA benchmark blood pump with different coefficient sets and scaler shear stress models used in the power-law hemolysis model.

Hemolysis is the most important issue to consider in the design and optimization of blood-contacting devices. Although the use of Computational Fluid Dynamics (CFD) in hemolysis prediction studies provides convenience and has promising potential, it is an extremely challenging process. Hemolysis predictions with CFD depend on the mesh, implementation method, coefficient set, and scalar-shear-stress model. To this end, an attempt was made to find the combination that would provide the most accurate result in hemolysis prediction with the commonly cited power-law based hemolysis model. In the hemolysis predictions conducted using CFD on the Food and Drug Administration (FDA) benchmark blood pump, 3 different scalar-shear-stress models, and 5 different coefficient sets with the power-law based hemolysis model were used. Also, a mesh independence test based on hemolysis and pressure head was performed. The pressure head results of CFD simulations were compared with published pressure head of the FDA benchmark blood pump and a good agreement was observed. In addition, results of CFD-hemolysis predictions which are conducted with scalar-shear-stress model and coefficient set combinations were compared with experimental hemolysis data at three operating conditions such as 6-7L/min flow rates at 3500rpm rotational speeds and 6L/min at 2500rpm. One of the combinations of the scalar-shear-stress model and the coefficient set was found to be within the error limits of the experimental measurements, while all other combinations overestimated hemolysis.

Abstract Mo029: Angiotensin Receptor Neprilysin Inhibitor improves cardiac hemodynamics and arrhythmogenesis in mitral regurgitation-induced heart failure—from bench to bedside approach

Background: An effective pharmacological therapy yet exists to mitigate mitral regurgitation (MR) progression before surgical interventions. Angiotensin receptor-neprilysin inhibitor (ARNI) has shown benefits in patients with heart failure (HF) with reduced ejection fraction but whether it can reduce degenerative MR induced HF remains unknown. Aim: To test the efficacy of ARNI in improving cardiac function in degenerative MR -induced HF. Methods and Results: Using a novel rat model of degenerative MR, we found that in degenerative MR rats treated with ARNI, left ventricular function and atrial fibrillation (AF) inducibility significantly improved compared with those treated with valsartan only (fractional shortening 61.78±7.26% v.s. 55.07±6.63%). Also, the AF inducibility dropped from 83% to 25% in MR rats treated with ARNI than to 58.3% in those treated with valsartan. ARNI treatment demonstrated significant attenuation of fibrosis, apoptosis and expressions of endoplasmic reticulum (ER) stress and inflammatory-associated proteins in left ventricular tissues compared to valsartan. In a shear stress model, the treatment of ARNI on cardiomyocytes mitigated the flow triggered activities of both ER released calcium signaling and store-operated calcium entry. Parallelly, we prospectively included 65 patients with degenerative MR receiving valsartan or ARNI and followed Holter electrocardiography and echocardiography. Notably, compared with degenerative MR patients receiving valsartan, those prescribed with ARNI demonstrated significant improvement of AF burden and left ventricular function. Conclusions: Degenerative MR-induced cardiovascular dysfunction including fibrosis and ER stress, were more effectively attenuated by ARNI treatment than by valsartan alone.

Numerical modeling on wave–current flows and bed shear stresses over an algal reef

AbstractCurrents in coastal zones under multiple mechanisms in terms of tides, waves, wind, and high roughness are difficult to model; bed shear stresses under wave–current flows are particularly challenging yet not being well studied. Few studies reported the modeling and validation of the bed shear stress in reef environments. In this paper, we present the first direct assessment of numerical modeling on depth-averaged currents and bed shear stresses over an algal reef using a coupled wave–current model (Delft-3D). The modeled results were validated and compared to the field observed data. The model considers hydrodynamic forcing in terms of tides, waves, wind stresses, and bed friction. Results show that the model generally reproduces the depth-averaged currents and bed shear stresses when considering all the mechanisms. Two numerical cases with and without wind forcing were tested to examine the effects of the winds. We found that the tide is mostly the primary factor driving the current, even in shallow waters within a depth of 3 m; however, the currents are also significantly affected by wind speeds and wind directions during high-wind events. When the wind direction is in the same direction as the tidal current, the current speed increases, suggesting the importance of the wind stress on the coastal currents. In addition, two models were chosen to study the nonlinear enhancement of bed shear stress by waves. We found a significant difference between the two models in predicting the bed shear stresses compared to the observed data. Nonlinear contribution from wave enhances the magnitude of bed shear stresses, which reduces the model error. The results highlight the nonlinear interaction between waves and currents is meaningful in predicting the bed shear stresses during high-wave-orbital motions; improvement of the present wave-current nonlinear interaction model for predicting the bed shear stresses may be needed.

Microscopic modeling and experimental investigation of inner surface magnetorheological polishing based on particle micromechanics

Traditional surface polishing methods are no longer able to meet the ultra-precision requirements of the high-tech industry for the inner surface of the pipe, but magnetorheological polishing technology is very suitable due to its advantages of high precision, fast controllability and good deposition stability. However, there is even less investigation on the microforce analysis, chaining mechanism and micromodeling methods of magnetorheological polishing fluid (MRPF), and the polishing mechanism of MRPF has not been explored yet. As a step to completely develop the magnetorheological polishing (MRP) technique, this paper proposed the simulation method of MRPF based on particle dynamics, and the shear stress model of magnetorheological fluid (MRF) is optimized under the action of the magnetic field after performing the chain simulation. On the basis of the optimized shear stress model and the hexagonal close-packed structure of MRF, the holding mechanism of polishing abrasive particles is explored for the MRPF and the corresponding holding models are proposed. Then, the shear yield stress models and material removal model are also established for the inner surface polishing, respectively. Eventually, the above theoretical analysis and related models have been verified though the polishing experiment of the titanium alloy pipe.

Enhancing Wind Turbine Blade Preventive Maintenance Procedure through Computational Fluid Dynamics-Based Prediction of Wall Shear Stress

Wind turbine blades are essential parts of wind energy systems and are frequently exposed to harsh environmental elements, such as strong winds, turbulence, and corrosive atmospheric elements. Over time, these circumstances may result in serious harm to blades, such as delamination and erosion, which may negatively affect the wind turbine’s functionality and durability. Accurate prediction of various types of damage is crucial to improve the toughness and lifespan of wind turbine blades and to maximize the overall effectiveness of wind energy systems. This article presents a novel computational fluid dynamics (CFDs)-based method for analyzing the distribution of wall shear stress on turbine blades, aimed at publicizing the yearly maintenance procedure. The investigation results from the CFDs, when compared with the current situation in a wind turbine farm in Thailand, confirmed that our wall shear stress modeling accurately predicted wind turbine damage. A maximum wall shear stress level higher than 5.00 Pa in the case of PA 90°, incoming air velocity 10.00 m/s, and 15 rpm was the main contribution to presenting the erosion and delamination from current drone inspection in wind turbine farms. In conclusion, these findings demonstrated the potential of using CFDs to predict wind turbine blade delamination and erosion, thereby significantly contributing to the development of specific and accurate yearly preventive maintenance. The proposed CFDs-based approach should serve as a sustainability tool for local human development, benefiting wind turbine engineers and operating technicians by providing them with a deeper understanding of the local flow conditions and wall shear stress distribution along wind turbine blades. This enables them to make informed decisions regarding blade design and maintenance.

A simple prediction of time-mean and wave orbital velocities in submerged canopy

Flow within submerged canopies influences the transport of nutrients, sediment, pollutants, plant seeds and the settlement of larvae. To improve our understanding of mass transport within canopies, a simple model is proposed to predict the total time-varying velocity within submerged rigid canopies (representing coral reefs) and flexible canopies (representing seagrasses and saltmarshes). The model divides the momentum equations into a canopy layer and free-stream layer. The difference in the time derivative of the velocity between the two layers is balanced by the sum of the shear stress and canopy drag, both of which depend on the in-canopy total velocity. The present model extended the shear stress model developed for steady current to combined current and wave conditions without additional calibrating coefficients. The model agreed well with the in-canopy velocity measured in the present and several previous studies. Importantly, the proposed model significantly improved the accuracy of canopy time-mean velocity prediction, which reduced the root mean square error by more than 50 %, compared with previous models. The model revealed that the addition of waves can significantly decrease the in-canopy time-mean velocity.

Low shear stress induces macrophage infiltration and aggravates aneurysm wall inflammation via CCL7/CCR1/TAK1/ NF-κB axis

BackgroundThis study aimed to elucidate the mechanism by which wall shear stress (WSS) influences vascular walls, accounting for the susceptibility of intracranial aneurysms (IAs) to rupture. MethodWe collected blood samples from the sacs of 24 ruptured and 28 unruptured IAs and analyzed the expression of chemokine CCL7 using enzyme-linked immunosorbent assay (ELISA). Univariate and multivariate logistic regression analyses were employed to assess clinical data, aneurysm morphology, and hemodynamics in both groups. Pearson correlation analysis investigated the relationship between CCL7 expression in aneurysm sac blood and WSS. Additionally, we established a bionic cell parallel plate co-culture shear stress model and a mouse low shear stress (LSS) model. The model was modulated using CCL7 recombinant protein, CCR1 inhibitor, and TAK1 inhibitor. We further evaluated CCL7 expression in endothelial cells and the levels of TAK1, NF-κB, IL-1β, and TNF-α in macrophages. Subsequently, the intergroup differences in expression were calculated. ResultsCCL7 expression was significantly higher in the ruptured group compared to the unruptured group. Hemodynamic analysis indicated that WSS was an independent predictor of the risk of aneurysm rupture. A negative linear correlation was observed between CCL7 expression and WSS. Upon addition of CCL7 recombinant protein, upregulation of CCR1 expression and increased levels of p-TAK1 and p-p65 were observed. Treatment with CCR1 and TAK1 inhibitors reduced inflammatory cytokine expression in macrophages under LSS conditions. Overexpression of TAK1 significantly alleviated the inhibitory effects of CCR1 inhibitors on p-p65 and inflammatory cytokines. ConclusionLSS prompts endothelial cells to secrete CCL7, which, upon binding to the macrophage surface receptor CCR1, stimulates the release of macrophage inflammatory factors via the TAK1/NF-κB signaling pathway. This process exacerbates aneurysm wall inflammation and increases the risk of aneurysm rupture.

Estimating the contribution of thermal processes to generation of seismicity in the junction zone of the Chuya Depression and the Kyrgyz Ridge of the Northern Tien Shan

This paper presents the attempt to estimate the contribution of thermal processes in the Earth’s crust to the generation of seismicity by the example of the seismically active region of the Chuya depression and the Kyrgyz ridge of the Northern Tien Shan. For this purpose, we use the models of temperature, petrophysical characteristics, and elastic moduli constructed in the previous works. The silica content model based on seismic tomography data is used to build a lithotype model. The constructed thermal conductivity model is utilized, along with the temperature model, to create a depth model of heat flow density. Based on the density, elastic moduli, and temperature models, the shear and thermoelastic stress models in the rocks are constructed. Their comparison with the distribution of earthquake hypocenters suggests that on the scale of the seismically active zone of the Chuya depression and the Kyrgyz ridge of the Northern Tien Shan, seismicity is mainly determined by the thermomechanical effect resulting from the hot ascending flow of acid magma from the upper mantle beneath the Muyunkum–Narat median mass.

Histone demethylase PHF8 protected against chondrocyte injury and alleviated posttraumatic osteoarthritis by epigenetically enhancing WWP2 expression.

Aberrant mechanical forces were considered as an important factor for osteoarthritis (OA) pathogenesis. Plant homeodomain finger-containing protein 8 (PHF8) participated in osteogenic differentiation and inflammatory progression. However, the role of PHF8 in aberrant force-related OA remains to be elucidated. In this study, a fluid shear stress (FSS) model in ATDC5 cells and an anterior cruciate ligament transection (ACLT) animal model were constructed. The results revealed the decrease of PHF8 in aberrant force-induced cartilage damage in vitro and in vivo. PHF8 overexpression alleviated the aberrant force-induced cell apoptosis, extracellular matrix degradation, and inflammation. Chromatin immunoprecipitation (ChIP) assays demonstrated that PHF8 epigenetically regulated WWP2 expression through demethylating H3K9me2 at WWP2 promoter, which was influenced by FSS treatment. C-X-C chemokine receptor type 4 (CXCR4) was identified as a potential substrate of WWP2. Co-immunoprecipitation (Co-IP) and ubiquitination experiments further demonstrated WWP2 decreased the stability of CXCR4 via the ubiquitination pathway. Subsequently, rescue experiments validated reintroduction of WWP2 significantly attenuated the effects of PHF8 deletion on FSS-induced chondrocyte injury, and CXCR4 overexpression reversed the protective effects of WWP2 overexpression on chondrocyte injury in FSS-treated ATDC5 cells. Moreover, delivery of a PHF8 adeno-associated virus (AAV) into articular cartilage remarkably ameliorated the breakdown of cartilage matrix by ACLT in mice. In conclusion, our findings highlighted the importance of PHF8/WWP2/CXCR4 signaling pathway in aberrant force-induced cartilage injury, which might provide a novel insight on future epigenetic-based treatment of posttraumatic OA.

EGCG, a Green Tea Compound, Increases NO Production and Has Antioxidant Action in a Static and Shear Stress In Vitro Model of Preeclampsia.

Preeclampsia (PE) is a gestational hypertensive disease characterized by endothelial dysfunction. Epigallocatechin-3-gallate (EGCG), the main compound in green tea, is a promising therapeutic target for the disease. By activating eNOS, EGCG increased NO production and exerted an important antioxidant action, but its specific impact in the context of PE remains understudied. The aim of this study is to evaluate the effects of EGCG on endothelial function in static and shear stress in in vitro models of PE. Endothelial cells were incubated with healthy (HP) and preeclamptic (PE) pregnant women's plasma, and the latter group was treated with EGCG. Additionally, NOS (L-NAME) and PI3K protein (LY249002) inhibitors were also used. The levels of NO, ROS, and O2•- were evaluated, as well as the antioxidant potential. These investigations were also carried out in a shear stress model. We found that EGCG increases the NO levels, which were reduced in the PE group. This effect was attenuated with the use of L-NAME and LY249002. Furthermore, EGCG increased the antioxidant capacity of PE, but its action decreased with LY294002. In cells subjected to shear stress, EGCG increased nitrite levels in the PE group and maintained its action on the antioxidant capacity. This is the first study of the effects of EGCG in this experimental model, as well as the investigation of its effects along with shear stress. Our findings suggest that EGCG improves parameters of endothelial dysfunction in vitro, making it a promising target in the search for treatments for the disease.

A novel approach for wall-boundary immersed flow simulation (part 2: modeling of wall shear stress)

A novel approach for wall-boundary immersed flow simulation (part 2: modeling of wall shear stress)

Errata: A novel approach for wall-boundary immersed flow simulation (part 2: modeling of wall shear stress) [Journal of Fluid Science and Technology Vol. 19, No. 3 (2024), JSFT0026

Errata: A novel approach for wall-boundary immersed flow simulation (part 2: modeling of wall shear stress) [Journal of Fluid Science and Technology Vol. 19, No. 3 (2024), JSFT0026

An improved formula for bed-load rate in open channel flows with emergent vegetation

There is an urgent need to predict the bed-load transport rate in vegetated river ecosystems to support restoration efforts. In response, we have developed a novel model for estimating the effective shear stress acting on the riverbed. This model is based on the energy equation and considers the intrinsic relationship between energy loss in the mean flow and turbulence generated by vegetation in open channel flows with emergent vegetation. Using this bed shear stress model, we assessed the performance of the Meyer-Peter–Müller (MPM) formula in predicting the bed-load transport rate in vegetated flows by comparing it with collected literature experimental data. The results revealed that the MPM formula does not provide accurate predictions. It tends to overestimate the bed-load transport rate when the dimensionless effective shear stress is approximately less than one and underestimate them when the dimensionless effective shear stress is approximately greater than one. This suggests that vegetation enhances and decreases the sediment transport rate when the dimensionless effective shear stress is approximately larger or lower than one, respectively. Consequently, we modified the coefficients of the MPM formula using extensive experimental data, leading to the development of a novel predictive formula for the bed-load transport rate in vegetated flows. This new formula outperforms existing literature equations and is effective for predicting the bedload transport rate, even for umbrella-like vegetation.

Numerical Investigation of Flow and Scour around Complex Bridge Piers in Wind–Wave–Current Conditions

A sea-crossing bridge is typically constructed in a marine environment with complex piers, and is susceptible to severe scour at the foundation. This study presents a numerical investigation on flow and scour around a complex pier, specifically focusing on a real-world sea-crossing bridge in China. A comprehensive CFD model incorporating hydrodynamic, free surface, sediment transport, and morphological models is employed for numerical modeling. Additionally, a wind shear stress model is considered to accurately simulate wind generation. The validation of the CFD model is achieved through comparison with experimental data of scour around a cylinder, demonstrating its capability to accurately replicate scour morphology and the temporal evolution of scour depth. Subsequently, the validated model is utilized for full-scale simulation of scour around the complex bridge pier under different wind, wave, and current conditions. The results indicate that compared to single piers with uniform cross-sectional shapes, flow patterns around complex piers are much more complicated. Scour predominantly occurs around the first row of group piles, while downstream piles experience less scour due to the sheltering effect from upstream piles. Furthermore, it becomes evident that the current exerts greater influence on pier scour than waves and wind, while the latter two factors primarily influence the superstructure of the bridge.

A comparative study of machine learning approaches for predicting viscosity in Sacran/CNF solutions

This study explores the predictive capabilities of machine learning algorithms for Sacran/CNF solution viscosity. Leveraging data on initial CNF amount, shear rate, shear stress, speed, torque, and temperature, decision tree (DT) and random forest (RF) models were employed. Both models demonstrated strong alignment with experimental results, achieving a coefficient of determination (R2) of 0.99. The main features of those models established that 60 % of the importance is attributed to speed or shear rate, while 31 % is attributed to the ratio of CNF content in the solutions. This research underscores the potential of machine learning in effectively analyzing viscosity datasets.

Evaluation of theoretical models describing the scratch adhesion test based on experimental results

ABSTRACT The article presents theoretical analyzes and results of research work on the verification of a commonly used method for measuring the adhesion of thin hard and superhard coatings applied to the blades of cutting tools made of various materials. The paper contains descriptions of the scratch adhesion test, taking into account, inter alia, interfacial shear stress model, energy model and friction model. This method of measuring the adhesion of the coating to the substrate has been presented in theoretical analyzes as a means of transmitting the appropriate amount of mechanical force or energy necessary to detach the coating from the substrate. The article discusses three approaches to describe the scratch adhesion test. The presented approaches concern: 1) “delivering” the critical shear stresses in the transition layer between the coating and the substrate according to the Benjamin and Weaver model, 2) “providing” additional elastic energy according to the Lauger model (according to the energy balance equation), 3) taking into account the criterion of deformation and matching of the coating to the substrate in the place of the indentation in the Burnett and Rickerby model by using the interfacial restraint parameter . The above-mentioned theoretical models have been experimentally verified for coatings deposited on substrates (cutting blades) made of high-speed steel, sintered carbides, white (oxide) tool ceramics, black tool ceramics (mixed, oxide-carbide), and cubic boron nitride. The tests of the adhesion force were performed with the use of proprietary equipment. This device is equipped with two channels for acoustic wave measurement by using a sound level meter and vibration measurement by a piezoelectric transducer. The adhesion force at which detaches the coating from the substrate was estimated by measuring the actual value of the amplitude of the vibration acceleration signal, its absolute value and the vibration signal intensity. The obtained values were verified using images from the profilometer and scanning electron microscope (SEM). In the final part of the article, the described models were verified on the basis of own experimental research and conclusions were formulated resulting from the analyzes and comparisons made. So far, mentioned models have not been assessed in terms of which of them allows the most accurate determination of the actual value of the critical load for the selected coating-substrate system. This was done in this article on the basis of own research. This innovative approach is a basic achievement.

Blood flow patterns switch VEGFR2 activity through differential S-nitrosylation and S-oxidation

Vascular endothelial growth factor receptor-2 (VEGFR2) plays a key role in maintaining vascular endothelial homeostasis. Here, we show that blood flows determine activation and inactivation of VEGFR2 through selective cysteine modifications. VEGFR2 activation is regulated by reversible oxidation at Cys1206 residue. H2O2-mediated VEGFR2 oxidation is induced by oscillatory flow in vascular endothelial cells through the induction of NADPH oxidase-4 expression. In contrast, laminar flow induces the expression of endothelial nitric oxide synthase and results in the S-nitrosylation of VEGFR2 at Cys1206, which counteracts the oxidative inactivation. The shear stress model study reveals that disturbed blood flow operated by partial ligation in the carotid arteries induces endothelial damage and intimal hyperplasia in control mice but not in knock-in mice harboring the oxidation-resistant mutant (C1206S) of VEGFR2. Thus, our findings reveal that flow-dependent redox regulation of the VEGFR2 kinase is critical for the structural and functional integrity of the arterial endothelium.

Form Shear Stress Correction in Bulk Flow Analysis of Grooved Seals Based on Effective Film Thickness

Abstract Bulk flow models for grooved annular seals provide computationally efficient static and dynamic response predictions, though heavy reliance on empirical relationships often leads to undesirable levels of uncertainty. The flow complexity caused by the grooves adds difficulty to shear stress modeling for these seals. This study seeks to improve shear stress modeling for grooved seals through the identification and quantification of the additional bulk flow shear stress contributions within the groove region. Through single groove computational fluid dynamics (CFD) simulations and an effective film thickness analysis framework, the additional groove shear stress component is identified as a form shear stress (FSS) due to its clear relationship to the effective film thickness behavior. The FSS is quantified as a correction to traditional shear stress definitions. Predictive models for the FSS are developed as functions of the ratio of circumferential to axial Reynolds number and the total resultant Reynolds number. Implementation of the FSS models into a simplified bulk flow method delivers leakage predictions for three seal cases within 10% of the experimental results and qualitative agreement in predicted circumferential velocity profiles while eliminating the need for an assumed groove loss coefficient. This is the first paper to utilize an effective film thickness-based procedure to quantify and model the FSS component in grooved seal bulk flow analysis. The demonstrated predictive capability and widespread applicability of the models and approach presented in this paper provide an avenue for significant improvements in grooved seal bulk flow prediction accuracy through improved shear stress modeling.