Solids dispersion is an important part of hydraulic fracturing, both in helping to understand phenomena such as tip screen-out and spreading of the pad, and in new process variations such as cyclic pumping of proppant. Whereas many frac fluids have low viscosity, e.g. slickwater, others transport proppant through increased viscosity. In this context, one method for influencing both dispersion and solids-carrying capacity is to use a yield stress fluid as the frac fluid. We propose a model framework for this scenario and analyse one of the simplifications. A key effect of including a yield stress is to focus high shear rates near the fracture walls. In typical fracturing flows this results in a large variation in shear rates across the fracture. In using shear-thinning viscous frac fluids, flows may vary significantly on the particle scale, from Stokesian behaviour to inertial behaviour across the width of the fracture. Equally, according to the flow rates, Hele-Shaw style models give way at higher Reynolds number to those in which inertia must be considered. We develop a model framework able to include this range of flows, while still representing a significant simplification over fully three-dimensional computations. In relatively straight fractures and for fluids of moderate rheology, this simplifies into a one-dimensional model that predicts the solids concentration along a streamline within the fracture. We use this model to make estimates of the streamwise dispersion in various relevant scenarios. This model framework also predicts the transverse distributions of the solid volume fraction and velocity profiles as well as their evolutions along the flow part.
more »
« less
Fluid deformation in random steady three-dimensional flow
The deformation of elementary fluid volumes by velocity gradients is a key process for scalar mixing, chemical reactions and biological processes in flows. Whilst fluid deformation in unsteady, turbulent flow has gained much attention over the past half-century, deformation in steady random flows with complex structure – such as flow through heterogeneous porous media – has received significantly less attention. In contrast to turbulent flow, the steady nature of these flows constrains fluid deformation to be anisotropic with respect to the fluid velocity, with significant implications for e.g. longitudinal and transverse mixing and dispersion. In this study we derive an ab initio coupled continuous-time random walk (CTRW) model of fluid deformation in random steady three-dimensional flow that is based upon a streamline coordinate transform which renders the velocity gradient and fluid deformation tensors upper triangular. We apply this coupled CTRW model to several model flows and find that these exhibit a remarkably simple deformation structure in the streamline coordinate frame, facilitating solution of the stochastic deformation tensor components. These results show that the evolution of longitudinal and transverse fluid deformation for chaotic flows is governed by both the Lyapunov exponent and power-law exponent of the velocity probability distribution function at small velocities, whereas algebraic deformation in non-chaotic flows arises from the intermittency of shear events following similar dynamics as that for steady two-dimensional flow.
more »
« less
- Award ID(s):
- 1654009
- NSF-PAR ID:
- 10094696
- Date Published:
- Journal Name:
- Journal of Fluid Mechanics
- Volume:
- 855
- ISSN:
- 0022-1120
- Page Range / eLocation ID:
- 770 to 803
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
We study the mixing dynamics of solute blobs in the flow through saturated heterogeneous porous media. As the solute plume is advected through a heterogeneous porous medium it suffers a series of deformations that determine its mixing with the ambient fluid through diffusion. Key questions are the relation between the spatial disorder and the mixing dynamics and the effect of the initial solute distribution. To address these questions, we formulate the advection–diffusion problem in a coordinate system that moves and rotates along streamlines of the steady flow field. The impact of the medium heterogeneity is quantified systematically within a stochastic modelling approach. For a simple shear flow, the maximum concentration of a blob decays asymptotically as $t^{-2}$ . For heterogeneous porous media, the mixing of the solute blob is determined by the random sampling of flow and deformation heterogeneity along trajectories, a mechanism different from persistent shear. We derive explicit perturbation theory expressions for stretching-enhanced solute mixing that relate the medium structure and mixing behaviour. The solution is valid for moderate heterogeneity. The random sampling of shear along trajectories leads to a $t^{-3/2}$ decay of the maximum concentration as opposed to an equivalent homogeneous medium, for which it decays as $t^{-1}$ .more » « less
-
null (Ed.)more » « less
Abstract Flapping insect wings experience appreciable deformation due to aerodynamic and inertial forces. This deformation is believed to benefit the insect’s aerodynamic force production as well as energetic efficiency. However, the fluid-structure interaction (FSI) models used to estimate wing deformations are often computationally demanding and are therefore challenged by parametric studies. Here, we develop a simple FSI model of a flapping wing idealized as a two-dimensional pitching-plunging airfoil. Using the Lagrangian formulation, we derive the reduced-order structural framework governing wing’s elastic deformation. We consider two fluid models: quasi-steady Deformable Blade Element Theory (DBET) and Unsteady Vortex Lattice Method (UVLM). DBET is computationally economical but does not provide insight into the flow structure surrounding the wing, whereas UVLM approximates flows but requires more time to solve. For simple flapping kinematics, DBET and UVLM produce similar estimates of the aerodynamic force normal to the surface of a rigid wing. More importantly, when the wing is permitted to deform, DBET and UVLM agree well in predicting wingtip deflection and aerodynamic normal force. The most notable difference between the model predictions is a roughly 20° phase difference in normal force. DBET estimates wing deformation and force production approximately 15 times faster than UVLM for the parameters considered, and both models solve in under a minute when considering 15 flapping periods. Moving forward, we will benchmark both low-order models with respect to high fidelity computational fluid dynamics coupled to finite element analysis, and assess the agreement between DBET and UVLM over a broader range of flapping kinematics.
-
We develop a theory of fluid--structure interaction (FSI) between an oscillatory Newtonian fluid flow and a compliant conduit. We consider the canonical geometries of a 2D channel with a deformable top wall and an axisymmetric deformable tube. Focusing on the hydrodynamics, we employ a linear relationship between wall displacement and hydrodynamic pressure, which has been shown to be suitable for a leading-order-in-slenderness theory. The slenderness assumption also allows the use of lubrication theory, and the flow rate is related to the pressure gradient (and the tube/wall deformation) via the classical solutions for oscillatory flow in a channel and in a tube (attributed to Womersley). Then, by two-way coupling the oscillatory flow and the wall deformation via the continuity equation, a one-dimensional nonlinear partial differential equation (PDE) governing the instantaneous pressure distribution along the conduit is obtained, without \textit{a priori} assumptions on the magnitude of the oscillation frequency (\textit{i.e.}, at arbitrary Womersley number). We find that the cycle-averaged pressure (for harmonic pressure-controlled conditions) deviates from the expected steady pressure distribution, suggesting the presence of a streaming flow. An analytical perturbative solution for a weakly deformable conduit is obtained to rationalize how FSI induces such streaming. In the case of a compliant tube, the results obtained from the proposed reduced-order PDE and its perturbative solutions are validated against three-dimensional, two-way-coupled direct numerical simulations. We find good agreement between theory and simulations for a range of dimensionless parameters characterizing the oscillatory flow and the FSI, demonstrating the validity of the proposed theory of oscillatory flows in compliant conduits at arbitrary Womersley number.more » « less
-
null (Ed.)Abstract ‘Horizontal convection’ (HC) is the generic name for the flow resulting from a buoyancy variation imposed along a horizontal boundary of a fluid. We study the effects of rotation on three-dimensional HC numerically in two stages: first, when baroclinic instability is suppressed and, second, when it ensues and baroclinic eddies are formed. We concentrate on changes to the thickness of the near-surface boundary layer, the stratification at depth, the overturning circulation and the flow energetics during each of these stages. Our results show that, for moderate flux Rayleigh numbers ( $O(1{0}^{11} )$ ), rapid rotation greatly alters the steady-state solution of HC. When the flow is constrained to be uniform in the transverse direction, rapidly rotating solutions do not support a boundary layer, exhibit weaker overturning circulation and greater stratification at all depths. In this case, diffusion is the dominant mechanism for lateral buoyancy flux and the consequent buildup of available potential energy leads to baroclinically unstable solutions. When these rapidly rotating flows are perturbed, baroclinic instability develops and baroclinic eddies dominate both the lateral and vertical buoyancy fluxes. The resulting statistically steady solution supports a boundary layer, larger values of deep stratification and multiple overturning cells compared with non-rotating HC. A transformed Eulerian-mean approach shows that the residual circulation is dominated by the quasi-geostrophic eddy streamfunction and that the eddy buoyancy flux has a non-negligible interior diabatic component. The kinetic and available potential energies are greater than in the non-rotating case and the mixing efficiency drops from ${\sim }0. 7$ to ${\sim }0. 17$ . The eddies play an important role in the formation of the thermal boundary layer and, together with the negatively buoyant plume, help establish deep stratification. These baroclinically active solutions have characteristics of geostrophic turbulence.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1017/jfm.2018.654
