Viscous Drag Research Articles

An Analytical Formulation for Correcting the Relative Permeability of Gas‐Water Flow in Propped Fractures Considering the Effect of Brinkman Flow

AbstractGas‐water flow in propped fractures can be commonly observed in various practical applications, including hydrocarbon development, geothermal exploitation, contaminant transport, and geological carbon storage. The fluid flow in a propped fracture can be regarded as Darcy type if only the resistance from the propping materials (e.g., cement in natural fractures and proppant‐pack in hydraulic fractures) accounts. However, if the fracture width is sufficiently small, the extra resistance from the viscous shear of fracture walls cannot be neglected, resulting in the appearance of Brinkman flow. In practice, during the development of a naturally fractured aquifer or a hydraulicly fractured hydrocarbon reservoir, the fracture width will be significantly reduced as the production proceeds. Therefore, the Brinkman flow can impose a strong effect on the fluid transportation within the fractures. However, the existing study about the two‐phase Brinkman flow in propped fractures is still far from adequate. In this work, on the basis of a modified two‐phase Brinkman equation, the authors derive a novel analytical formulation to correct the relative permeability of gas‐water two‐phase flow in propped fractures to account for the Brinkman effect. With the aid of the proposed formulation, the authors carry out a comprehensive investigation of the influence of Brinkman flow on the effective gas‐water relative permeability and well performance. The calculated results show that the effect of Brinkman flow on the water phase (or gas phase) transportation is more significant with a larger water (or gas) saturation. A narrower propped fracture is more likely to induce Brinkman flow, thus leading to a lower relative permeability for both water and gas phases. As the propping‐material permeability is increased, the fluid transportation bears more severe viscous drag from fracture walls, and the relative permeability will be consequently reduced. Only if the fracture width is significantly reduced during the production, the Brinkman flow demonstrates its influence on the well performance. Otherwise, Darcy's law can provide sufficiently accurate results in characterizing the two‐phase flow in propped fractures.

Three-dimensional stagnation point motion of bioconvection nanofluid via moving stretching sheet with convective and anisotropic slip condition

The current article deals with the SM of a NFs on a 3D moving surface containing bioconvection, microorganisms, and convective conditions over anisotropic slip. This work reports describes the characteristics of bioconvection aspects of nanofluid containing microorganisms with anisotropic slip condition. The study may be important from the application point of view in microfluidic devices, antibiotics, enhanced oil recovery and many other areas. Applying the similarity variables, the basic gov. Eqs are translated to ordinary ones, and such Eqs are calculate by R-K-F fourth order programming with shooting scheme by MATLAB software. The physical parameters on the moving surface are explained in detail via graphs. We found that the temperature increases with a higher revolution and larger numerical numbers (convection parameter), while the HTR declined for various values of the SFP. The HTR enhances with surface convection, and therefore, the convection condition is helpful in reducing the viscous drag on the moving SS.

Numerical simulation of particles distribution of environmental tobacco smoke (ETS) and its concentration in a closed environment

Humans are frequently exposed to various chemical hazards, with smoking being a significant source of indoor air pollution. This pollution is particularly pronounced in environments such as smoking rooms, conference halls, cafes, and coffee houses. This study focuses on analyzing the concentration, distribution, and behavior of emitted particles from both mainstream and sidestream smoke in a controlled smoking room environment. A critical aspect of the research is the consideration of the thermal plume generated by the smoker's body, particularly around the breathing zone, which plays a crucial role in the dispersion and inhalation of pollutants. To conduct this analysis, a computational fluid dynamics approach was employed, incorporating a displacement ventilation (DV) system integrated with a detailed surface mesh model of a standing manikin. A Lagrangian approach was utilized to track the particle paths, considering factors such as particle size, density, and initial velocity. This approach was complemented by an Eulerian method to simulate the complex airflow patterns within the smoking room, accounting for turbulent flow and temperature variations. The integration of these methods provides a comprehensive understanding of both the macroscopic and microscopic behaviors of smoke particles. The motion equations governing the particles incorporated several forces, including inertial drag force, viscous drag force, buoyancy force, and gravitational force. Additionally, the study examined the interaction between thermal plumes and smoke particles, particularly how temperature gradients around the body affect particle dispersion and concentration. The study meticulously evaluated the distribution and volume average of respirable suspended particles (RSPs) and nicotine by analyzing a plane opposite to the breathing zone of both heated and unheated manikins. It was observed that the thermal plume generated by the temperature gradient around the body could significantly alter the flow field of both sidestream and mainstream smoke between puff periods. The simulation results demonstrated that the thermal plume increased the concentration of RSPs in the breathing zone, highlighting a critical exposure pathway for smokers and nonsmokers alike. A significant increase of approximately 70% in the nicotine concentration is due to the thermal plume. The concentration of nicotine peaks at approximately 0.003 mg/m3 by the third second in smoker#1's (6 smokers' position) breathing zone. Specific configurations of the DV system showed varying efficiencies in reducing pollutant concentrations, emphasizing the need for targeted ventilation strategies. In a chamber with six smokers, researchers found that the time required for suspended particles to reach a cleaning ratio, considering the dynamics of particle movement due to smoking activity, was approximately 184 s per puff.

Research on the impact of hydrodynamic parameter design deviation on the motion performance of fully actuated AUV

Abstract This study, based on the Simulink simulation platform, investigates the influence of hydrodynamic parameters on the navigation speed of fully actuated Autonomous Underwater Vehicles (AUVs) and conducts a sensitivity analysis of AUV motion performance to different design deviations in hydrodynamic parameters. The results show that the velocity of underwater motion in each direction exponentially decreases as the corresponding viscous drag coefficient increases while being less influenced by viscous drag coefficients in other directions. Additionally, the sensitivity of the viscous drag coefficient to the underwater motion of AUVs is much greater than that of the inertial drag coefficient.

Performance assessment of a Triple Coaxial-Cylinder Wave-Energy Converter (TCWEC)

We propose a novel Triple Coaxial-cylinder Wave-Energy Converter (TCWEC) system, as an evolution of the commonly deployed dual-coaxial cylinder WEC (DCWEC). TCWEC consists of one inner cylinder and two concentric outer cylinders, characterized by two coupled resonant frequencies, improving wave power absorption. A semi-analytical model, based on potential flow theory and matched eigenfunction expansions, is developed to analyze the 3-Degrees of Freedom (DOF) system. The viscous drag coefficients are determined by Computational Fluid Dynamics (CFD) simulations and employed as equivalent linear damping. Analysis of impact of the viscous flow-separation factor fvis on the capture width indicates that the higher energy-absorption capability of TCWEC is not merely attributed to a smaller fvis value, but rather a planned split of the outer cylinder in DCWEC design into two cylinders in TCWEC, so as to produce an effectively broader resonance property. Under identical sea conditions, it is found that, compared to DCWEC, the optimal capture width of TCWEC is increased by as much as 77% in regular waves. In irregular waves, the optimal capture width is increased up to 40% across the entire frequency range. These findings suggest that the energy absorption efficiency of TCWEC is promising and its potential for real-world applications.

Quantifying CSF Dynamics disruption in idiopathic normal pressure hydrocephalus using phase lag between transmantle pressure and volumetric flow rate

Background and purposeIdiopathic normal pressure hydrocephalus (iNPH) is a cerebrospinal fluid (CSF) dynamics disorder as evidenced by the delayed ascent of radiotracers over the cerebral convexity on radionuclide cisternography. However, the exact mechanism causing this disruption remains unclear. Elucidating the pathophysiology of iNPH is crucial, as it is a treatable cause of dementia. Improving the diagnosis and treatment prognosis rely on the better understanding of this disease. In this study, we calculated the pulsatile transmantle pressure and investigated the phase lag between this pressure and the volumetric CSF flow rate as a novel biomarker of CSF dynamics disruption in iNPH. Methods44 iNPH patients and 44 age- and sex-matched cognitively unimpaired (CU) control participants underwent MRI scans on a 3T Siemens scanner. Pulsatile transmantle pressure was calculated analytically and computationally using volumetric CSF flow rate, cardiac frequency, and aqueduct dimensions as inputs. CSF flow rate through the aqueduct was acquired using phase-contrast MRI. The aqueduct length and radius were measured using 3D T1-weighted anatomical images. ResultsPeak pressure amplitudes and the pressure load (integrated pressure exerted over a cardiac cycle) were similar between the groups, but the non-dimensionalized pressure load (adjusted for anatomical factors) was significantly lower in the iNPH group (p<0.001, Welch's t-test). The phase lag between the pressure and the flow rate, arising due to viscous drag, was significantly higher in the iNPH group (p<0.001). ConclusionThe increased phase lag is a promising new biomarker for quantifying CSF dynamics dysfunction in iNPH. Statement of SignificanceThe exact mechanism causing the disruption of CSF circulation in idiopathic normal pressure hydrocephalus (iNPH) remains unclear. Elucidating the pathophysiology of iNPH is crucial, as it is a treatable cause of dementia. In this study, we provided an analytical and a computational method to calculate the pulsatile transmantle pressure and the phase lag between the pressure and the volumetric CSF flow rate across the cerebral aqueduct. The phase lag was significantly higher in iNPH patients than in controls and may serve as a novel biomarker of CSF dynamics disruption in iNPH.

Analysis of Surface Drag Reduction Characteristics of Non-Smooth Jet Coupled Structures

To enhance the service life of shipping equipment and minimize surface wear, this study employs biomimetic principles, integrating fitted structures with jet dynamics to model three configurations: non-smooth structures, single jet structures, and non-smooth jet-coupled structures. We utilized the SST k-ω turbulence model for numerical simulations to investigate the drag reduction characteristics of these structural models. By varying the jet angle and speed, we analyzed the changes in viscous resistance, pressure differential resistance, and drag reduction rates at the wall surface. Furthermore, the mechanisms of compressive stress, velocity fields, vortex structures, and shear stress on drag-reducing surfaces were elucidated, revealing how these factors contribute to drag reduction in non-smooth jet-coupled structures. The results indicate that the non-smooth jet-coupled structure exhibits superior drag reduction performance at a main flow field velocity of 20 m/s. As the jet velocity increases, the viscous drag on the surface of the non-smooth jet-coupled structure decreases, while the pressure differential drag increases. Conversely, variations in the jet angle have a minimal effect on viscous drag but lead to a reduction in pressure differential drag. Specifically, when the jet velocity is set at 1 m/s, and the jet angle is 60°, the drag reduction achieved by the non-smooth jet-coupled structure peaks at 7.48%. Additionally, the non-smooth jet-coupled structure features a larger area characterized by low shear stress, along with an increased boundary layer thickness at the bottom; this configuration effectively reduces surface velocity and consequent viscous drag.

Electronically actuated artificial hinged cilia for efficient bidirectional pumping.

Cilial pumping is a potent mechanism used to control and manipulate fluids on microscales. Recently, we introduced an electronically driven μ-cilial platform that can create arbitrary flow patterns in liquids near a surface with the potential for various engineering applications. This μ-cilial platform, however, utilized the coupling between elasticity and viscous drag to obtain pumping and had several limitations. For example, each cilium could only pump in one direction. Thus, to create bidirectional flows, it was necessary to fabricate and separately actuate two oppositely facing cilia. As another example, the generation of non-reciprocal cilial motions, a necessary condition for pumping at these scales, could only be achieved by matching the elastic stresses inherent in actuating the cilia with the viscous drag forces generated by the flows. This criterion severely restricted the frequency range over which the cilia could be operated and resulted in a small swept area, both of which restricted the volume of fluid being pumped in each cycle. These limitations contrast with the capabilities of natural cilia, which can achieve omnidirectional transport and operation over a broad range of frequencies. In natural cilia, these capabilities arise from their complex internal structure. Inspired by this strategy we designed hinged cilia and show they can achieve bidirectional pumping of larger fluid volumes over a broad range of frequencies. Finally, we demonstrate that even regular arrays of individually controlled hinged cilia can generate a variety of flow patterns using fewer cilia than in previous cilia metasurface designs.

A variable fidelity approach for predicting aerodynamic wall quantities of hypersonic vehicles using the ConvNeXt encoder-decoder framework

Computational fluid dynamics (CFD) simulations for obtaining three-dimensional hypersonic vehicle aerodynamic characteristics are resource-intensive. Deep learning offers a promising alternative for predicting aerodynamic wall quantities but typically requires a large dataset, which conflicts with the intention to reduce CFD costs. To address this, we propose a novel prediction model leveraging variable fidelity data to alleviate computational demands. The model utilizes an encoder-decoder architecture with ConvNeXt blocks as operators, processing variable fidelity data as inputs and outputs. We also developed a U-Net with residual blocks (ResUNet) for performance comparison. Transformation and topology patching techniques were applied to tackle the challenges posed by large grid volumes and complex topologies in three-dimensional vehicles. Results demonstrate that the ConvNeXt Encoder-Decoder predicts peak regions more accurately than ResUNet and aligns closely with CFD in low-value areas. The ConvNeXt Encoder-Decoder maintains maximum heat flux errors below 1 % and viscous drag errors below 3.81 % across varying angles of attack. It exhibits superior fitting performance with minimal deviation from CFD compared to ResUNet. Prediction accuracy decreases under multiple inflow changes compared to the angle of attack variations alone. Heat flux predictions show high consistency with CFD, with relative errors below 6.38 %, whereas friction predictions exhibit higher errors with 10.28 % maximum error and 0.28 % minimum error. The model accurately predicts friction variation trends in peak regions at the blunt edge's central plane but performs poorly in low-value areas. In summary, the ConvNeXt Encoder-Decoder accurately predicts the wall quantities under varying multiple inflow conditions.

The hydrodynamic noise suppression of the towed array by the combination of the flow and sustained release

Abstract This paper proposes a jet flow release method based on fish mucus to reduce the flow noise of the towed array. This method aims to alleviate the adverse effects of turbulent fluctuating pressure on the self-noise of a towed array at high speeds and address the problem of target signal masking caused by flow noise. A circular jet nozzle is installed at the front end of the towed array to spray mucus backward, reducing the turbulent fluctuating pressure on the surface of the towed array and improving the signal-to-noise ratio. Large eddy simulation using Fluent and finite element calculation of the acoustic field using Actran are employed to obtain flow and acoustic field results, respectively. The research findings demonstrate that after spraying mucus through the jet nozzle, the model’s viscous drag is reduced, and the excitation on the surface is significantly decreased, leading to a noticeable reduction in the self-noise sound pressure level in the frequency range of 150–2000 Hz.

Translations of multiple interacting bubbles in one-dimensional system

Abstract The pulsations and translations of N cavitation bubbles obey combined ordinary differential equations on one-dimensional chain by analyzing the pressures near the cavitation bubbles, and the translations of the cavitation bubbles are discussed. The results indicate that the cavitation bubbles with the same ambient radii at both ends move symmetrically towards or away from the middle cavitation bubble, and the secondary Bjerknes forces on the cavitation bubbles at both ends are equal in magnitude and opposite in direction and the secondary Bjerknes force on the middle cavitation bubble is zero. In addition, dynamic balance can be maintained between bubbles due to the viscous drag force, and the secondary Bjerknes forces on the cavitation bubbles at both ends are equal in magnitude and opposite in direction, and the secondary Bjerknes force on the middle cavitation bubble is almost zero. This is of immediate consequence for understanding cavitation structure formation processes in an acoustic field.

Mathematical derivation of a unified equations for adjoint lattice Boltzmann method in airfoil and wing aerodynamic shape optimization

Unified equations of the adjoint lattice Boltzmann method (ALBM) are derived for five applicable objective functions in 2D/3D aerodynamic shape optimization problems. The derived equations include the adjoint equation, boundary condition, terminal condition and gradient of the cost function. In this research, firstly, these relations are extracted for each objective in details and then the general form of ALBM equations are presented for all defined practical aerodynamic objective function. Five applicable cost functions which are the most important objectives in optimization of aerodynamic geometries include desired pressure and viscous shear stress (VSS) inverse design, drag and moment at fixed lift and finally lift to drag ratio at fixed angle of attack. The new extracted relations are based on the circular and spherical function scheme, and are valid for viscous/inviscid, compressible/incompressible and 2D/3D flows in all continuous flow regimes. Proof of new extracted general relations have been performed by authors.

Experimental study of wave diffraction loads on a vertical circular cylinder with heave plates at deep and shallow drafts

This study focuses on wave diffraction loads influenced by heave plates at deep and shallow submerged depths. An extensive experimental campaign is conducted on a fixed vertical and free-surface piercing circular cylinder with five different circular plates, including a perforated plate. The horizontal and vertical wave forces are studied in monochromatic and bichromatic waves. The results show that the heave plate significantly influences the forces varying in time and representative harmonics, including sum- and difference-frequency components. The study found that wave steepness and submerged depth of the plate contribute to the nonlinearity of the loads and that the porous plate reduces the loads compared to the solid plate. Corresponding results from a boundary element method (BEM) solver, HydroStar, provide a reasonable prediction of the harmonics. However, the first harmonic of the vertical force is underpredicted at low frequencies, specifically where the BEM results indicate the cancellation of the forces. Dedicated computational fluid dynamics (CFD) simulations are carried out to investigate the discrepancy, and it is observed that flow separation occurs around the edge of the plate. A simplified approach based on Morison’s equation is employed to model the flow separation effects using a quadratic drag force calculated with a constant drag coefficient. The approach shows that viscous drag is the dominant contributor to vertical wave forces on the heave plate in the low-frequency regime.

Pattern formation on hydrophobic PDMS substrates by evaporating droplets: combined effect of substrate stiffness, relative humidity and particle size

Desiccation patterns left by micro-droplets of water impregnated with particles on hydrophobic substrates have been analyzed with respect to variations in the elastic stiffness of the substrates, particle size and relative humidity. The complex and unique patterns obtained, have been analyzed and explained in terms of the time scales of moving Triple Phase Line (TPL) on substrate and substrate relaxation rate. The rate of TPL movement is found to depend on the relative humidity and substrate stiffness. In turn, this affects the contact angle hysteresis. Particle movement is a result of viscous drag and inertia apart from electrostatic interactions. We have successfully explained the myriad patterns obtained from drying droplets via systematic rheological measurements along with an understanding of the role of all the effective forces and their time scales of action.

Field-Directed Motion, Cargo Capture, and Closed-Loop Controlled Navigation of Microellipsoids.

Microrobots have the potential for diverse applications, including targeted drug delivery and minimally invasive surgery. Despite advancements in microrobot design and actuation strategies, achieving precise control over their motion remains challenging due to the dominance of viscous drag, system disturbances, physicochemical heterogeneities, and stochastic Brownian forces. Here, a precise control over the interfacial motion of model microellipsoids is demonstrated using time-varying rotating magnetic fields. The impacts of microellipsoid aspect ratio, field characteristics, and magnetic properties of the medium and the particle on the motion are investigated. The role of mobile micro-vortices generated is highlighted by rotating microellipsoids in capturing, transporting, and releasing cargo objects. Furthermore, an approach is presented for controlled navigation through mazes based on real-time particle and obstacle sensing, path planning, and magnetic field actuation without human intervention. The study introduces a mechanism of directing motion of microparticles using rotating magnetic fields, and a control scheme for precise navigation and delivery of micron-sized cargo using simple microellipsoids as microbots.

Nonmodal stability analysis of Poiseuille flow through a porous medium

We unravel the nonmodal stability of a three-dimensional nonstratified Poiseuille flow in a saturated hyperporous medium constrained by impermeable rigid parallel plates. The primary objective is to broaden the scope of previous studies that conducted modal stability analysis for two-dimensional disturbances. Here, we explore both temporal and spatial transient disturbance energy growths for three-dimensional disturbances when the Reynolds number and porosity of the material are high, based on evolution equations with respect to time and space, respectively. Modal stability analysis reveals that the critical Reynolds number for the onset of shear mode instability increases as porosity increases. Moreover, the Darcy viscous drag term stabilizes shear mode instability, resulting in a delay in the transition from laminar flow to turbulence. In addition, it demonstrates the suppression of three-dimensional shear mode instability as the spanwise wavenumber increases, thereby confirming the statement of Squire’s theorem. By contrast, nonmodal stability analysis discloses that both temporal and spatial transient disturbance energy growths curtail as the effect of the Darcy viscous drag force intensifies. But their maximum values behave like O(Re2) for a fixed porous material, where Re is the Reynolds number. However, for different porous materials, the scalings for both temporal and spatial transient disturbance energy growths are different. Furthermore, increasing porosity also suppresses both temporal and spatial disturbance energy growths. Finally, we observe that temporal transient disturbance energy growth becomes larger for a spanwise perturbation, while spatial transient disturbance energy growth becomes larger for a steady perturbation when angular frequency vanishes. The initial disturbance that excites the largest temporal energy amplification generates two sets of alternating high-speed and low-speed elongated streaks in the streamwise direction.

An improved delayed detached-eddy study on the aerodynamic braking technique based on blunting the streamlined section of the high-speed train

This paper utilizes the improved delayed detached-eddy simulation method to investigate an aerodynamic braking technique involving blunting the streamlined portion of a high-speed train (HST) at Re = 5.0 × 105. The accuracy of the numerical simulation method was validated through reduced-scale wind tunnel experiments at the same Reynolds number level. The study compares aerodynamic drag, pressure distribution, boundary layer, and wake flow characteristics between the original configuration and the braking configuration of the HST. Additionally, the impact of aerodynamic braking plates on the flow characteristics around the key components of the HST has also been discussed. The results indicate a significant increase in the pressure drag experienced by the HST with the application of aerodynamic braking plates to its streamlined sections, while there is a slight decrease in viscous drag. This leads to a remarkable 235.4% rise in the overall aerodynamic drag of the entire HST. The aerodynamic braking plates also have a substantial impact on the turbulent wake flow topology, significantly increasing turbulence levels in the near-wake region. Furthermore, the implementation of aerodynamic braking plates may affect pantograph current collection by significantly altering stream-wise and vertical velocity components, notably increasing velocity fluctuation around the contact position between the pantograph and power supply lines.

Motor domain of condensin and step formation in extruding loop of DNA.

During the asymmetric loop extrusion of DNA by a condensin complex, one domain of the complex stably anchors to the DNA molecule, and another domain reels in the DNA strand into a loop. The DNA strand in the loop is fully relaxed, or there is no tension in the loop. Just outside of the loop, there is a tension that resists the extrusion of DNA. To maintain the extrusion of the DNA loop, the condensin complex must have a domain capable of generating a force to overcome the tension outside of the loop. This study proposes that the groove-shaped HEAT repeat domain Ycg1 plays the role of a molecular motor. A DNA molecule may bind to the groove electrostatically, and the weak binding force facilitates the random thermal motion of DNA molecules. A mechanical model that random collisions between DNA and the nonparallel inner surfaces of the groove may generate a directional force which is required for the loop extrusion to sustain. The hinge domain binds to the DNA molecule and acts as an anchor during asymmetric DNA loop extrusion. When the effects of ATP hydrolysis and the viscous drag of the fluid environment are considered, the motor-anchor model for the condensin complex and the mechanical model might explain the asymmetric loop extrusion, the formation of steps, the step size distribution in the loop extrusion, the tension-dependent extrusion speed, the interaction between coexisting loops on the DNA strand, and untying the knots during extrusion. This model can also explain the observed formation of the Z-loop.

Low-frequency electrodynamics in the mixed state of superconducting NbN and a-MoGe films using two-coil mutual inductance technique

We investigate the low-frequency electrodynamics in the vortex state of two type-II superconducting films, namely, a moderate-to-strongly pinned Niobium Nitride (NbN) and a very weakly pinned amorphous Molybdenum Germanium (a-MoGe) film. We employ a two-coil mutual inductance technique to extract the complex penetration depth, λ~ . The sample response is studied through the temperature variation of λ~ in the mixed state, where we employ a model developed by Coffey and Clem (the CC model) to extract different vortex lattice (VL) parameters, such as the restoring pinning force constant (Labusch parameter), VL drag coefficient and pinning potential barrier. We observe that a consistent description of the inductive and dissipative part of the response is only possible when we take the viscous drag on the vortices to be several orders of magnitude larger than the viscous drag estimated from the Bardeen–Stephen model.

Microscale Diffractive Lenses Integrated into Microfluidic Devices for Size-Selective Optical Trapping of Particles.

Integration of optical components into microfluidic devices can enhance particle manipulations, separations, and analyses. We present a method to fabricate microscale diffractive lenses composed of aperiodically spaced concentric rings milled into a thin metal film to precisely position optical tweezers within microfluidic channels. Integrated thin-film microlenses perform the laser focusing required to generate sufficient optical forces to trap particles without significant off-device beam manipulation. Moreover, the ability to trap particles with unfocused laser light allows multiple optical traps to be powered simultaneously by a single input laser. We have optically trapped polystyrene particles with diameters of 0.5, 1, 2, and 4 μm over microlenses fabricated in chromium and gold films. Optical forces generated by these microlenses captured particles traveling at fluid velocities up to 64 μm/s. Quantitative trapping experiments with particles in microfluidic flow demonstrate size-based differential trapping of neutrally buoyant particles where larger particles required a stronger trapping force. The optical forces on these particles are identical to traditional optical traps, but the addition of a continuous viscous drag force from the microfluidic flow introduces tunable size selectivity across a range of laser powers and fluid velocities.