Shear Flow Research Articles

The minimal seed for transition to convective turbulence in heated pipe flow

It is well known that buoyancy suppresses, and can even laminarise, turbulence in upward heated pipe flow. Heat transfer seriously deteriorates in this case. A new direct numerical simulation model is established to simulate flow-dependent heat transfer in an upward heated pipe. The model shows good agreement with experimental results. Three flow states are simulated for different values of the buoyancy parameter $C$ : shear turbulence, laminarisation and convective turbulence. The latter two regimes correspond to the heat transfer deterioration regime and the heat transfer recovery regime, respectively (Jackson & Li 2002; Bae et al., Phys. Fluids, vol. 17, issue 10, 2005; Zhang et al., Appl. Energy, vol. 269, 2020, 114962). We confirm that convective turbulence is driven by a linear instability (Su & Chung, J. Fluid Mech., vol. 422, 2000, pp. 141–166) and that the deteriorated heat transfer within convective turbulence is related to a lack of rolls near the wall, which leads to weak mixing between the flow near the wall and the centre of the pipe. Having surveyed the fundamental properties of the system, we perform a nonlinear non-modal stability analysis, which seeks the minimal perturbation that triggers a transition from the laminar state. Given the differences between shear and convective turbulence, we aim to determine how the nonlinear optimal (NLOP) changes as the buoyancy parameter $C$ increases. We find that at first, the NLOP becomes thinner and closer to the wall. Most importantly, the critical initial energy $E_0$ required to trigger turbulence keeps increasing, implying that attempts to trigger it artificially may not be an efficient means to improve heat transfer at larger $C$ . At $C=6$ , a new type of NLOP is discovered, capable of triggering convective turbulence from lower energy, but over a longer time. It is active only in the centre of the pipe. We next compare the transition processes, from linear instability and by the nonlinear non-modal excitation. At $C=4$ , linear instability leads to a state that approaches a travelling wave solution or periodic solutions, while the minimal seed triggers shear turbulence before decaying to convective turbulence. Deeper into the parameter space for convective turbulence, at $C=6$ , the new nonlinear optimal triggers convective turbulence directly. Detailed analysis of the periodic solution at $C=4$ reveals three stages: growth of the unstable eigenfunction, the formation of streaks, and the decay of the streaks. The stages of the cycle correspond to changes in the linear instability of the turbulent mean velocity profile. Unlike the self-sustaining process for classical shear flows, where the streak is disrupted via instability, here, decay of the streak is more closely linked to suppression of the linear instability of the mean flow, and hence suppression of the rolls. Flow visualisations at $C$ up to $10$ also show similar processes, suggesting that the convective turbulence in the heat transfer recovery regime is sustained by these three typical processes.

Infiltration Characteristics and Hydrodynamic Parameters in Response to Topographic Factors in Bare Soil Surfaces, Laboratory Experiments Based on Cropland Fields of Purple Soil in Southwest China

Topography is an important factor that impacts the hydrological processes on sloping farmlands. Yet, few studies have reported the combined influences of slope gradient and slope position on infiltration characteristics and hydrodynamic parameters on sloping croplands in purple soil regions, an important area for agricultural productivity in Southwest China. Here, laboratory-simulated rainfall experiments were conducted in a steel trough (5 m long, 2 m wide, and 0.45 m deep), and rainfall lasted for 1 h at a rate of 90 mm h−1 to examine the variations in the infiltration rates and hydrodynamic parameters under varying slope gradients (i.e., 3°, 6°, 10°, 15°, 21°, and 27°) and slope positions (i.e., upper, middle, and lower), and explore the relationships between the infiltration rate and the soil detachment rate. The results showed that the infiltration rate decreased gradually with duration rainfall and ultimately approached a steady state in the six slope treatments. Cumulative infiltration ranged from 15.54 to 39.32 mm during rainfall, and gradually reduced with the increase of slope gradient. The Horton’s model outperforms other models for predicting the infiltration rate with an R2 value of 0.86. Factors such as Darcy–Weisbach friction, flow shear force, Manning friction coefficient, unit energy, and runoff depth varied in the following order: upper slope > middle slope > lower slope, whilst the Reynolds number and Froude number gradually increased along the slope transect from the upper to lower slope positions. A significant linear function was fitted between the soil detachment rate and the infiltration rate at the gentle slopes (3°, 6°, 10°), whereas an exponential relationship was observed at the steep slopes (15°, 21°, and 27°). Observation also suggested that 15° was the critical slope gradient of sediment detachment, infiltration characteristics, and hydrodynamic parameters. Our results provide theoretical insight for developing models that predict the impacts of topographic factors on hydrological characteristic and soil erosion in hilly agricultural landscapes of purple soil fields.

Shearmetry of Fluids with Tunable Rheology by Polarized Luminescence of Rare Earth-Doped Nanorods.

Shear stress plays a critical role in regulating physiological processes within microcirculatory systems. While particle imaging velocimetry is a standard technique for quantifying shear flow, uncertainty near boundaries and low resolution remain severe restrictions. Additionally, shear stress determination is particularly challenging in biofluids due to their significant non-Newtonian behaviors. The present study develops a shearmetry technique in physiological settings using a biomimetic fluid containing rare earth-doped luminescent nanorods acting in two roles. First, they are used as colloidal additives adjusting rheological properties in physiological media. Their anisotropic morphology and interparticle interaction synergistically induce a non-Newtonian shear-thinning effect emulating real biofluids. Second, they can probe shear stress due to the shear-induced alignment. The polarized luminescence of the nanorods allows for quantifying their orientational order parameter and thus correlated shear stress. Using scanning confocal microscopy, we demonstrate the tomographic mapping of the shear stress distribution in microfluidics. High shear stress is evident near the constriction and the cellular periphery, in which non-Newtonian effects can have a significant impact. This emerging shearmetry technique is promising for implementation in physiological and rheological environments of biofluids.

Stokes drift and particle trajectories induced by surface waves atop a shear flow

Surface waves and currents are crucial to the mass transfer in the air-sea interaction as they can drive a variety of dynamical processes. How mass can be transported by surface waves and current coupling is addressed through a study of their induced motions of fluid parcels. To this end, a weakly nonlinear wavetrain is imposed on the background flow whose direction and magnitude are permitted to vary with water depth and second-order features of this configuration are investigated. A leading-order approximation to the Stokes drift is derived, correct to the second order in wave steepness, and applicable to an arbitrarily depth-dependent background flow. The reduced forms of the approximate Stokes drift are provided in a few limiting cases such as a current with an exponential profile or propagating in an orthogonal direction to the wave propagation. Novel features related to the Stokes drift and particle trajectories have been reported for the first time as a result of the rotation induced by the wave and current coupling. A non-vanishing component of the Stokes drift velocity and net-mean displacement of fluid parcels in the span-wise direction to the wave propagation are observed in the cases where a shear current propagates obliquely to the waves direction. A non-monotonic dependence on water depth of the stream-wise component of the Stokes drift is shown, and thereby the largest mass transport induced no longer occurs on the still water surface but some depth beneath. The non-monotonic behavior occurs beyond the regime of the near-irrotational assumption of wave-induced motions. It can also lead to the change of the signs for the stream-wise Stokes drift throughout the water column, and thus an overall cancellation of the integrated mass transport by waves over the water column, indicating that the depth-integrated models can likely lead to underestimated effects of the mass transport which is non-trivial at a local depth. The results from this study have far-reaching impact. The Stokes drift profile is a direct input to the parametrization of the surface waves forcing in ocean circulations and the obliquely propagating Stokes drift can be plausibly responsible for the formation of oblique Langmuir rolls to wave propagation in the open ocean.

Aerodynamic study of a leading-edge rotating cylinder in a cambered aerofoil (NACA2412)

Purpose Flow separation over an aircraft’s wing beyond a specific angle of attack is challenging. Flow boundary layer manipulation has been investigated to improve aerofoil lift and mitigate flow separation difficulties including stall and drag. This is solved via active or passive flow control. Active flow control method moving surface boundary (MSB) enhances shear flow momentum, making it effective. MSB is easier than suction and blowing. Asymmetrical airfoils, which generate lift in aircraft wings, have received less MSB research than symmetrical ones. The purpose of this study is to asses the design efficacy of MSB’s NACA 2412 aerofoil. Design/methodology/approach To observe the performance of MSB in NACA 2412, a computational model has been created, and aerodynamical performance has been analyzed. Findings The study results show that the NACA 2412 with MSB has better aerodynamic efficiency than the NACA 2412 base design. It works best when it reaches its optimal speed and the delay in flow separation works well. Research limitations/implications Limitations may include specific aerofoil applicability, external factors and simulation constraints. Implications guide future research for broader insights and applications. Practical implications Improving asymmetrical aerofoil performance, mitigating stall effects, reducing drag and optimizing designs with moving surface boundary. Insights gained can enhance overall aircraft efficiency and flow control techniques. Originality/value The MSB flow control in a chambered aerofoil is less explored and not explored enough in wing-based aerofoils, and the optimal cylinder speed ratio trend has been discussed at each angle of attack studied.

I-mode plasma confinement improvement by real-time lithium injection and its classification on EAST tokamak

Abstract I-mode is a promising regime for future fusion reactors due to the high energy confinement and the moderate particle confinement. However, the effect of lithium, which has been widely applied for particle recycling and impurity control, on I-mode plasma is still unclear. Recently, experiments of real-time lithium powder injection on I-mode plasma have been carried out in EAST Tokamak. It was found that the confinement performance of the I-mode can be improved by the lithium powder injection, which can strongly reduce electron turbulence (ET) and then trigger ion turbulence (IT). And it was observed that the ET intensity is inversely proportional to the velocity shear, which suggests that the injection of lithium powder leads to a gradual enhancement of the shear flow, whereby the turbulence is reduced and consequently the confinement is improved. Four different regimes of I-mode have been identified in EAST. The Type I I-mode plasma is characterized by the weakly coherent mode (WCM) and the geodesic-acoustic mode (GAM). The Type II I-mode is featured as the WCM and the edge temperature ring oscillation (ETRO). The Type III I-mode corresponds to the plasma with the co-existence of ETRO, GAM, and WCM. The Type IV I-mode denotes the plasma with only WCM but without ETRO and GAM. It was observed that the WCM intensity is increased with lithium powder injection by the confinement improvement/pedestal temperature increase. EAST experiments demonstrate that lithium powder injection is an effective tool for real-time control and confinement improvement of I-mode plasma.

First Principles Variational Formulation and Analytical Solutions for Third Order Shear Deformable Beams with Simple Supports using Fourier Series Method

The Fourier series method for solving bending problems of third-order shear deformable beams (TSBD) is presented in this paper. The theory accounts for transverse shear deformation and is suitable for moderately thick beams. Transverse shear stress-free conditions are valid at the beam surfaces for TSDB. The field equations are two coupled ordinary differential equations in terms of two unknown displacement functions – transverse deflection w(x) and warping function For simply supported ends considered, the loading and unknown functions w(x) and are represented in Fourier series that satisfies the boundary conditions. The problem is reduced to a system of algebraic equations in terms of the Fourier coefficients wn and which is solved to obtain wn and Axial and transverse displacements fields, axial bending stress and transverse shear stress are then determined for the cases of point load at midspan, uniformly and linearly distributed loads over the span. The results are identical with previous results obtained by other scholars who used Ritz variational methods and other mathematical tools. For moderately thick beams under uniform load with l/h = 4, the results obtained for is 0.516% greater than the exact result by Timoshenko and Goodier while for thick beams under uniform load with l/h = 2, the result for is 1.96% greater than the exact result by Timoshenko and Goodier. Similar acceptable variation was obtained for for beams under linearly distributed load with the variations being less than 2% for l/h = 2. However, unacceptable variations of 68.37% were found for in thick beams under point load, for l/h = 2.0. The variation for however reduced to 12.513% for moderately thick beams under point load for l/h =4.

Development of Asymmetric Convection in a Tropical Cyclone under Environmental Uniform Flow and Vertical Wind Shear

AbstractIdealized numerical simulations have been carried out to reveal the complexity in the development of asymmetric convection in a tropical cyclone (TC) under the influence of an environment with either uniform flow, vertical wind shear (VWS), or both. Results show that rainwater is enhanced to the right of the motion in the outer rainband, but such enhancement occurs in the upshear-left area of the inner-core region. Additionally, due to the asymmetries introduced by environmental flow, wavenumber-1 temperature and height anomalies develop at a radius of ~1000 km in the upper levels. A sub-vortex aside from the TC center encompassing the wavenumber-1 warm center appears, and asymmetric horizontal winds emerge, which, in turn, changes the storm-scale (within 400 km) VWS. Deep convection in the inner core closely follows the changing storm-scale VWS when its magnitude is larger than 2 m s–1 and is located downshear of the storm-scale VWS in all the experiments with environmental flow. In the outer rainbands, the maximum boundary layer convergence is mainly controlled by the direction of motion and is located in the rear-right quadrant. These results extend upon the findings of previous studies in three aspects: (1) The discovery of the roughly linear combination effect from the uniform flow and large-scale VWS; (2) The development of upper-level asymmetric winds on a 1000-km scale through the interaction between the TC vortex and environmental flow, resulting in changes in the storm-scale VWS pattern within the TC area; (3) The revelation that TC asymmetric convection closely aligns with the direction-varying storm-scale VWS instead of the initially designated VWS.

Shear-driven magnetic buoyancy in the solar tachocline: Dependence of the mean electromotive force on diffusivity and latitude

Abstract The details of the dynamo process that is responsible for driving the solar magnetic activity cycle are still not fully understood. In particular, whilst differential rotation provides a plausible mechanism for the regeneration of the toroidal (azimuthal) component of the large-scale magnetic field, there is ongoing debate regarding the process that is responsible for regenerating the Sun’s large-scale poloidal field. Our aim is to demonstrate that magnetic buoyancy, in the presence of rotation, is capable of producing the necessary regenerative effect. Building upon our previous work, we carry out numerical simulations of a local Cartesian model of the tachocline, consisting of a rotating, fully compressible, electrically conducting fluid with a forced shear flow. An initially weak, vertical magnetic field is sheared into a strong, horizontal magnetic layer that becomes subject to magnetic buoyancy instability. By increasing the Prandtl number we lessen the back reaction of the Lorentz force onto the shear flow, maintaining stronger shear and a more intense magnetic layer. This in turn leads to a more vigorous instability and a much stronger mean electromotive force, which has the potential to significantly influence the evolution of the mean magnetic field. These results are only weakly dependent upon the inclination of the rotation vector, i.e. the latitude of the local Cartesian model. Although further work is needed to confirm this, these results suggest that magnetic buoyancy in the tachocline is a viable poloidal field regeneration mechanism for the solar dynamo.

Scaling Relationships of the Structural and Rheological Behavior of Tadpole Polymer Chains in Dilute Solution Systems Using Brownian Dynamics Simulations.

Tadpole polymers, also known as lasso polymers, feature molecular structures that combine a single ring with a single linear side branch, leading to distinct conformational, dynamical, and rheological characteristics compared to their corresponding counterparts, particularly pure linear and pure ring polymers. To elucidate the mechanisms underlying these distinctive behaviors, comprehensive mesoscopic Brownian dynamics (BD) simulations of dilute solution systems of tadpole polymers were conducted using a bead-rod chain model under both equilibrium and flow conditions. Three types of tadpole polymer chains were prepared by varying the ring-to-linear ratio within the tadpole chain and comparing them with the corresponding linear and ring chains. Depending on this ratio, tadpole polymer chains exhibit entirely different structural properties and rotational dynamics, both in equilibrium and under shear flow. As the linear proportion within the tadpole chain increased, the structural, dynamic, and rheological properties of the tadpole polymer chains became more similar to those of pure linear polymers. Conversely, with an increasing ring proportion, these properties began to resemble those of pure ring polymers. Based on these observed tendencies, a simple general scaling expression is proposed for tadpole polymer properties that integrates scaling expressions for both pure linear and pure ring polymers. Our results indicate that the conformational, dynamic, and rheological properties of tadpole polymers, as predicted by these simple scaling expressions, are in good agreement with the simulated values, a result we consider statistically significant.

Comprehensive transverse shear stress analysis for FG-GRC sandwich plates based on a refined higher-order zigzag theory

Comprehensive transverse shear stress analysis for FG-GRC sandwich plates based on a refined higher-order zigzag theory

Heat transfer of MHD Eyring–Prandtl fluid flow past a moving permeable plate with shear flow

Heat transfer of MHD Eyring–Prandtl fluid flow past a moving permeable plate with shear flow

Effects of Silica Nanoparticles on the Piezoelectro-Elastic Response of PZT-7A-Polyimide Nanocomposites: Micromechanics Modeling Technique.

By using piezoelectric materials, it is possible to convert clean and renewable energy sources into electrical energy. In this paper, the effect on the piezoelectro-elastic response of piezoelectric-fiber-reinforced nanocomposites by adding silica nanoparticles into the polyimide matrix is investigated by a micromechanical method. First, the Ji and Mori-Tanaka models are used to calculate the properties of the nanoscale silica-filled polymer. The nanoparticle agglomeration and silica-polymer interphase are considered in the micromechanical modeling. Then, considering the filled polymer as the matrix and the piezoelectric fiber as the reinforcement, the Mori-Tanaka model is used to estimate the elastic and piezoelectric constants of the piezoelectric fibrous nanocomposites. It was found that adding silica nanoparticles into the polymer improves the elastic and piezoelectric properties of the piezoelectric fibrous nanocomposites. When the fiber volume fraction is 60%, the nanocomposite with the 3% silica-filled polyimide exhibits 39%, 31.8%, and 37% improvements in the transverse Young's modulus ET, transverse shear modulus GTL, and piezoelectric coefficient e31 in comparison with the composite without nanoparticles. Furthermore, the piezoelectro-elastic properties such as ET, GTL, and e31 can be improved as the nanoparticle diameter decreases. However, the elastic and piezoelectric constants of the piezoelectric fibrous nanocomposites decrease once the nanoparticles are agglomerated in the polymer matrix. A thick interphase with a high stiffness enhances the nanocomposite's piezoelectro-elastic performance. Also, the influence of volume fractions of the silica nanoparticles and piezoelectric fibers on the nanocomposite properties is studied.

Shear flow of nanofluid past a vertical plate in presence of Thompson–Troian slip and Stefan blowing

Shear flow of nanofluid past a vertical plate in presence of Thompson–Troian slip and Stefan blowing

Atomistic insights into scratch-induced structural evolution of silica glass

Understanding the mechanism of scratch damage is vital to developing better scratch-resistant glasses. To this extent, employing molecular dynamics simulations and experiments, we investigate the scratch damage of silica glass against a rigid diamond indenter. The glass surface is pre-indented to a constant depth and then dragged to simulate a linear scratch, and the structural impact in the indent-to-scratch transitioning phase is examined. We observe that despite the differences in length and timescales, the simulated values of indentation hardness and coefficient of friction exhibit excellent agreement with experimental values from nanoindentation and atomic force microscopy experiments. Interestingly, analysis of the subsurface deformation in the scratched region reveals densification and shear flow, in contrast to pure densification, as in the case of indentation. Furthermore, similar percentages of recovery from experiments and simulation reveal that the reversible component of plastic deformation owing to densification is comparable in both cases. Finally, in contrast to the common hypothesis, we demonstrate that while the bond angles and lengths recover significantly, the ring structure does not recover upon annealing, although they exhibit some relaxation. Thus, the present study sheds new light on the crucial role of the medium-range structure of glasses subjected to scratch deformation.

Free-swimming bacteria transcriptionally respond to shear flow

Surface-attached cells can sense and respond to shear flow, but planktonic (free-swimming) cells are typically assumed to be oblivious to any flow that carries them. Here, we find that planktonic bacteria can transcriptionally respond to flow, inducing expression changes that are beneficial in flow. Specifically, we use microfluidic experiments and quantitative modeling to show that in the presence of flow, planktonic Pseudomonas aeruginosa induce shear rate-dependent genes that promote growth in low-oxygen environments. Untangling this mechanism revealed that in flow, motile P. aeruginosa spatially redistribute, leading to cell density changes that activate quorum sensing, which in turn enhances the oxygen uptake rate. In diffusion-limited environments, including those commonly encountered by bacteria, flow-induced cell density gradients also independently generate oxygen gradients that alter gene expression. Mutants deficient in this flow-responsive mechanism exhibit decreased fitness in flow, suggesting that this dynamic coupling of biological and mechanical processes can be physiologically significant.

Rotational Taylor dispersion in linear flows

The coupling between advection and diffusion in position space can often lead to enhanced mass transport compared with diffusion without flow. An important framework used to characterize the long-time diffusive transport in position space is the generalized Taylor dispersion theory. In contrast, the dynamics and transport in orientation space remains less developed. In this work we develop a rotational Taylor dispersion theory that characterizes the long-time orientational transport of a spheroidal particle in linear flows that is constrained to rotate in the velocity-gradient plane. Similar to Taylor dispersion in position space, the orientational distribution of axisymmetric particles in linear flows at long times satisfies an effective advection–diffusion equation in orientation space. Using this framework, we then calculate the long-time average angular velocity and dispersion coefficient for both simple shear and extensional flows. Analytic expressions for the transport coefficients are derived in several asymptotic limits including nearly spherical particles, weak flow and strong flow. Our analysis shows that at long times the effective rotational dispersion is enhanced in simple shear and suppressed in extensional flow. The asymptotic solutions agree with full numerical solutions of the derived macrotransport equations and results from Brownian dynamics simulations. Our results show that the interplay between flow-induced rotations and Brownian diffusion can fundamentally change the long-time transport dynamics.

Can Microcavitated Slippery Surfaces Outperform Micropillared and Untextured?

Surface features' morphology is crucial in designing lubricant-infused slippery surfaces (LIS). Microcavities were hypothesized to provide lower physical pinning, reduced droplet normal adhesion, and superior lubricant retention as compared to micropillars and untextured surfaces. Micropillars and microcavities (h = 10 ± 3 μm, d = 8 ± 1 μm, p = 17 ± 3 μm, rw = 1.4 ± 0.2) were replicated on polydimethylsiloxane from Lotus leaf and were coated with 1000 cSt silicone oil films (530 nm-27 μm thick). Water wetting, water-oil thermodynamic stability, droplet's normal adhesion and oil shear drainage properties were investigated to evaluate the relative performance of microcavitated, micropillared and untextured LIS. For ≤7 μm thick oil films, cavitated and untextured LIS displayed superior slippery properties than micropillared LIS (16 ± 1°, 7 ± 1°, 30 ± 4° slide-off angles respectively). Also, normal adhesion is of the order: cavities < untextured < pillars, and smaller than their dry counterparts. Furthermore, the oil retention efficiency under the action of centrifugal forces and continuous shear flow of water is of the order: pillars > cavities > untextured. Thus, it can be concluded that microcavitated LIS can outperform micropillared and untextured LIS.

Two-species model for nonlinear flow of wormlike micelle solutions. Part I: Model

We develop a rheological model to approximate the nonlinear rheology of wormlike micelles using two constitutive models to represent a structural transition at high shear rates. The model is intended to describe the behavior of semidilute wormlike micellar solutions over a wide range of shear rates whose parameters can be determined mainly from small-amplitude equilibrium measurements. Length evolution equations are incorporated into reactive Rolie-Poly entangled-polymer rheology and dilute reactive-rod rheology, with a kinetic exchange between the two models. Although the micelle length is remarkably reduced during flow, surprisingly, we propose that they are not shortened by stress-enhanced breakage, which remains thermally driven. Instead, we hypothesize that stretching energy introduces a linear potential that decreases the rate of recombination and reduces the mean micelle length. This stress-hindered recombination approach accurately describes transient stress-growth upon start-up shear flow, and it predicts a transition of shear viscosity and alignment response observed at high shear rates. The proposed mechanism applies only when self-recombination occurs frequently. The effect of varying the relative rate of self-recombination on the rheology of wormlike micelles at high shear rates is yet to be explored.

Low Intraocular Pressure Induces Fibrotic Changes in the Trabecular Meshwork and Schlemm's Canal of Sprague Dawley Rats.

Continuous artificial aqueous humor drainage in the eyes of patients with glaucoma undergoing trabeculectomy likely exerts abnormal shear stress. However, it remains unknown how changes in intraocular pressure (IOP) can affect aqueous humor outflow (AHO). Here, we induced and maintained low intraocular pressure (L-IOP) in healthy Sprague Dawley (SD) rats by puncturing their eyes using a tube (200-µm diameter) for 2 weeks. After the rats were euthanized, their eyes were removed, fixed, embedded, stained, and scanned to analyze the physiological and pathological changes in the trabecular meshwork (TM) and Schlemm's canal (SC). We measured SC parameters using ImageJ software and assessed the expression of various markers related to flow shear stress (KLF4), fibrosis (TGF-β1, TGF-β2, α-SMA, pSmad1/5, pSmad2/3, and fibronectin), cytoskeleton (integrin β1 and F-actin), diastolic function (nitric oxide synthase and endothelial nitric oxide synthase [eNOS]), apoptosis (cleaved caspase-3), and proliferation (Ki-67) using immunofluorescence or immunohistochemistry. L-IOP eyes showed a larger SC area, higher eNOS expression, and lower KLF4 and F-actin expression in the TM and SC (both P < 0.05) than control eyes. The aqueous humor of L-IOP eyes had a higher abundance of fibrotic proteins and apoptotic cells than that of control eyes, with significantly higher TGF-β1, α-SMA, fibronectin, and cleaved caspase-3 expression (all P < 0.05). In conclusion, a persistence of L-IOP for 2 weeks may contribute to fibrosis in the TM and SC and might be detrimental to conventional AHO in SD rat eyes. Clinicians should consider that aberrant shear force induced by aqueous humor fluctuation may damage AHO outflow channel when treating patients.