skip to main content


Title: Boundary Conditions and Polymeric Drag Reduction for the Navier–Stokes Equations
Reducing wall drag in turbulent pipe and channel flows is an issue of great practical importance. In engineering applications, end-functionalized polymer chains are often employed as agents to reduce drag. These are polymers which are floating in the solvent and attach (either by adsorption or through irreversible chemical binding) at one of their chain ends to the substrate (wall). We propose a PDE model to study this setup in the simple setting where the solvent is a viscous incompressible Navier–Stokes fluid occupying the bulk of a smooth domain Ω⊂ℝ𝑑, and the wall-grafted polymer is in the so-called mushroom regime (inter-polymer spacing on the order of the typical polymer length). The microscopic description of the polymer enters into the macroscopic description of the fluid motion through a dynamical boundary condition on the wall-tangential stress of the fluid, something akin to (but distinct from) a history-dependent slip-length. We establish the global well-posedness of strong solutions in two-spatial dimensions and prove that the inviscid limit to the strong Euler solution holds with a rate. Moreover, the wall-friction factor ⟨𝑓⟩ and the global energy dissipation ⟨𝜀⟩ vanish inversely proportional to the Reynolds number 𝐑𝐞. This scaling corresponds to Poiseuille’s law for the friction factor ⟨𝑓⟩∼1/𝐑𝐞 for laminar flow and thereby quantifies drag reduction in our setting. These results are in stark contrast to those available for physical boundaries without polymer additives modeled by, for example, no-slip conditions, where no such results are generally known even in two-dimensions.  more » « less
Award ID(s):
1703997
NSF-PAR ID:
10271659
Author(s) / Creator(s):
;
Date Published:
Journal Name:
Archive for Rational Mechanics and Analysis
ISSN:
0003-9527
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. The effect of slip surfaces on the laminar–turbulent separatrix of plane Poiseuille flow is studied by direct numerical simulation. In laminar flows, the inclusion of the slip surfaces results in a drag reduction of over 10 %, which is in good agreement with previous studies and the theory of laminar slip flows. Turbulence lifetimes, the likelihood that turbulence is sustained, is investigated for transitional flows with various slip lengths. We show that slip surfaces decrease the likelihood of sustained turbulence compared to the no-slip case, and the likelihood is further decreased as slip length is increased. A more deterministic analysis of the effects of slip surfaces on a transition to turbulence is performed by using nonlinear travelling-wave solutions to the Navier–Stokes equations, also known as exact coherent solutions. Two solution families, dubbed P3 and P4, are used since their lower-branch solutions are embedded on the boundary of the basin of attraction of laminar and turbulent flows (Park & Graham, J. Fluid Mech. , vol. 782, 2015, pp. 430–454). Additionally, they exhibit distinct flow structures – the P3 and P4 are denoted as core mode and critical-layer mode, respectively. Distinct effects of slip surfaces on the solutions are observed by the skin-friction evolution, linear growth rate and phase-space projection of transitional trajectories. The slip surface appears to modify the transition dynamics very little for the core mode, but quite considerably for the critical-layer mode. Most importantly, the slip surface promotes different transition dynamics – an early and bypass-like transition for the core mode and a delayed and H- or K-type-like transition for the critical-layer mode. We explain these distinct transition dynamics based on spatio-temporal and quadrant analyses. It is found that slip surfaces promote the prevalence of strong wall-toward motions (sweep-like events) near vortex cores close to the channel centre, inducing an early transition, while long sustained ejection events are present in the region of the $\unicode[STIX]{x1D6EC}$ -shaped vortex cores close to the critical layer, resulting in a delayed transition. This should motivate flow control strategies to fully exploit these distinct transition dynamics for transition to turbulence. 
    more » « less
  2. We present direct numerical simulation results of a wave-current boundary layer in a current-dominated flow regime (wave driven to steady current ratio of 0.34) over bumpy walls for hydraulically smooth flow conditions (wave orbital excursion to roughness ratio of 10). The turbulent, wave-current channel flow has a friction Reynolds number of $350$ and a wave Reynolds number of $351$ . At the lower boundary, a bumpy wall is introduced with a direct forcing immersed boundary method, while the top wall has a free-slip boundary condition. Despite the hydraulically smooth nature of the wave-driven flow, the phase variations of the turbulent statistics for the bumpy wall case were found to vary substantially when compared with the flat wall case. Results show that the addition of weak waves to a steady current over flat walls has a negligible effect on the turbulence or bottom drag. However, the addition of weak waves to a steady current over bumpy walls has a significant effect through enhancement of the Reynolds stress (RS) accompanied by a drag coefficient increase of $11\,\%$ relative to the steady current case. This enhancement occurs just below the top of the roughness elements during the acceleration portion of the wave cycle: Turbulent kinetic energy (TKE) is subsequently transported above the roughness elements to a maximum height of roughly twice the turbulent Stokes length. We analyse the TKE and RS budgets to understand the mechanisms behind the alterations in the turbulence properties due to the bumpy wall. The results provide a mechanistic picture of the differences between bumpy and flat walls in wave-current turbulent boundary layers and illustrate the importance of bumpy features even in weakly energetic wave conditions. 
    more » « less
  3. Abstract

    The engineering problems involving clay under non‐isothermal conditions (e.g., geothermal energy harvest, landfill cover system, and nuclear waste disposal) are multiscale and multiphysics by nature. The nanoscale hydrodynamics of clay at elevated temperature is essential in developing a physics‐based multiscale model for clay under non‐isothermal conditions. The nonequilibrium molecular dynamics (NEMD) is a useful tool to study the nanoscale hydrodyndamics of clay. This article presents an NEMD modeling of hydrodynamics of clay nanopores at elevated temperatures. Water flow confined in pyrophyllite and montmorillonite clay nanopores is investigated. The nonequilibrium state is maintained by uniformly exerting an external force on each water molecule. The NEMD simulations have provided a molecular‐scale perspective of temperature effect on clay‐water density, water flow velocity, shear viscosity, clay‐water slip length, hydraulic conductivity, and clay‐water friction coefficient. The numerical results have shown a strong temperature dependence of fluid flow velocity, shear viscosity, clay‐water slip length, and hydraulic conductivity at the nanoscale. We have validated the applicability of cubic law in determining hydraulic conductivity at the nanopore scale at elevated temperatures. It is found from our numerical results that slip clay‐water boundary condition is an essential factor in properly determining nanoscale fluid flow velocity. By numerical examples, we also study the impact of nanopore size and clay layer thickness on the hydrodynamics of the clay‐water system.

     
    more » « less
  4. The topic of friction reduction has been studied through the decades for numerous engineering applications that involve internal and external flows. Inspired by the natural surface structure of different plants and animals, engineered microtexturing of surfaces is one of the effective ways of reducing the drag. By introducing different texture geometries, the flow behavior close to the solid boundary can be altered and thus manipulated towards achieving a reduced net drag force on the surface. Despite considerable research on the subject, most works have concentrated on optimization of the surface texturing for maximizing the friction reduction and minimizing the pumping power requirements, and less attention has been paid to characterization of the flow and boundary layer in the vicinity of the wall, especially in laminar regime. In this work we investigate the role that microtexturing has on friction reduction under low to moderate Reynolds numbers (Re). We perform a parametric study on the shape and dimensions of the surface textures and investigate the boundary layer and streamline behavior as well as the local shear stress and pressure distribution along the solid-fluid interface under different flow conditions. The outcomes of this work will provide a guideline for optimal design of artificial textures with major implications for many engineering applications such as microfluidic systems used in thermal management and biochemical diagnostics. 
    more » « less
  5. Many of the cell membrane's vital functions are achieved by the self-organization of the proteins and biopolymers embedded in it. The protein dynamics is in part determined by its drag. A large number of these proteins can polymerize to form filaments.In vitrostudies of protein–membrane interactions often involve using rigid beads coated with lipid bilayers, as a model for the cell membrane. Motivated by this, we use slender-body theory to compute the translational and rotational resistance of a single filamentous protein embedded in the outer layer of a supported bilayer membrane and surrounded on the exterior by a Newtonian fluid. We first consider the regime where the two layers are strongly coupled through their inter-leaflet friction. We find that the drag along the parallel direction grows linearly with the filament's length and quadratically with the length for the perpendicular and rotational drag coefficients. These findings are explained using scaling arguments and by analysing the velocity fields around the moving filament. We then present and discuss the qualitative differences between the drag of a filament moving in a freely suspended bilayer and a supported membrane as a function of the membrane's inter-leaflet friction. Finally, we briefly discuss how these findings can be used in experiments to determine membrane rheology. In summary, we present a formulation that allows computation of the effects of membrane properties (its curvature, viscosity and inter-leaflet friction), and the exterior and interior three-dimensional fluids’ depth and viscosity on the drag of a rod-like/filamentous protein, all in a unified theoretical framework.

     
    more » « less