Proppant Settling Correlations for Non-Newtonian Fluids Under Static and Dynamic Conditions

1982 ◽  
Vol 22 (02) ◽  
pp. 164-170 ◽  
Author(s):  
Subhash N. Shah

Abstract This paper presents a new approach for analysis of proppant sealing data in non-Newtonian pseudoplastic fracturing fluids and develops drag coefficient correlations as a function of fluid model parameter n', Results of experiments with these fluids under static as well as dynamic conditions are discussed. A wide range of n' and particle Reynolds number is investigated. It is shown that at low-particle Reynolds numbers the fluid model parameter n' has a significant effect on proppant settling velocity. This effect diminishes at higher particle Reynolds numbers. The dynamic settling velocity data agree reasonably well with the correlations developed from static velocity data. Further, experimental results of earlier investigations agree well with the correlations of this study. Introduction Many examples of flow around submerged objects appear in engineering work. In the oilfield industry, one example is particle or proppant transport in a fracture during hydraulic fracturing. Hydraulic fracturing has been in use commercially for 30 years. The formations are fractured hydraulically by pumping a slurry of viscous fracturing fluid and proppant at high pressures. The purpose of the proppant is to hold the fracture open at the end of the treatment. The production increase resulting from the treatment depends on fracture conductivity and final proppant distribution. Knowledge of settling velocity of a single particle in the fracture and streamlines around the particle is of great value in understanding the complex transport process leading to a particular proppant distribution in the fracture. A thorough understanding of proppant transpose would help design better fracturing treatments. The purpose of this investigation is two fold:to gather experimental data of proppant settling velocity in several non-Newtonian fracturing fluids, andto develop correlations between drag coefficient and particle Reynolds number. These correlations later can be used to predict the settling velocity of a particular particle in a non-Newtonian fracturing fluid of known rheological properties. Experience has shown that most fracturing fluids exhibit highly non-Newtonian fluid characteristics. They may possess completely different properties under shear than when at rest. For simplicity, most of the studies reported in the literature are undertaken by dropping proppant in a stagnant fluid. In the actual fracturing process, however, the proppant is settling while fluid is moving in the fracture. Thus, to validate the correlations derived from the settling velocity data in a stagnant fluid, some dynamic experiments also are conducted. Literature Review The subject of particle motion in Newtonian fluids has been studied extensively by many investigators. In contrast, very little has been accomplished in the case of particle motion in non-Newtonian fluids, particularly on proppant transport in the fracture. Most recently, a preliminary work on dynamic proppant transport in fracturing gels has been reported by Hannah and Harrington. Experiments were conducted with a concentric cylinder tester to measure the fall rate of proppant in non-Newtonian fracturing gels. The test device used was similar to one reported previously by Novotny. Their experimental data did not agree with the theoretical predictions, and no explanation has been given for the discrepancy. Later, using a similar tester, Harrington et al. measured the settling rate of proppant in cross-linked fracturing gels. SPEJ P. 164^

2021 ◽  
Author(s):  
Silin Jing ◽  
Xianzhi Song ◽  
Zhaopeng Zhu ◽  
Buwen Yu ◽  
Shiming Duan

Abstract Accurate description of cuttings slippage in the gas-liquid phase is of great significance for wellbore cleaning and the control accuracy of bottom hole pressure during MPD. In this study, the wellbore bubble flow environment was simulated by a constant pressure air pump and the transparent wellbore, and the settling characteristics of spherical particles under different gas volume concentrations were recorded and analyzed by highspeed photography. A total of 225 tests were conducted to analyze the influence of particle diameter (1–12mm), particle density (2700–7860kg/m^3), liquid viscosity and bubble volume concentration on particle settling velocity. Gas drag force is defined to quantitatively evaluate the bubble’s resistance to particle slippage. The relationship between bubble drag coefficient and particle Reynolds number is obtained by fitting the experimental results. An explicit settling velocity equation is established by introducing Archimedes number. This explicit equation with an average relative error of only 8.09% can directly predict the terminal settling velocity of the sphere in bubble containing Newtonian fluids. The models for predicting bubble drag coefficient and the terminal settling velocity are valid with particle Reynolds number ranging from 0.05 to 167 and bubble volume concentration ranging from 3.0% to 20.0%. Besides, a trial-and-error procedure and an illustrative example are presented to show how to calculate bubble drag coefficient and settling velocity in bubble containing fluids. The results of this study will provide the theoretical basis for wellbore cleaning and accurate downhole pressure to further improve the performance of MPD in treating gas influx.


Author(s):  
Patricia Streufert ◽  
Terry X. Yan ◽  
Mahdi G. Baygloo

Local turbulent convective heat transfer from a flat plate to a circular impinging air jet is numerically investigated. The jet-to-plate distance (L/D) effect on local heat transfer is the main focus of this study. The eddy viscosity V2F turbulence model is used with a nonuniform structured mesh. Reynolds-Averaged Navier-Stokes equations (RANS) and the energy equation are solved for axisymmetric, three-dimensional flow. The numerical solutions obtained are compared with published experimental data. Four jet-to-plate distances, (L/D = 2, 4, 6 and 10) and seven Reynolds numbers (Re = 7,000, 15,000, 23,000, 50,000, 70,000, 100,000 and 120,000) were parametrically studied. Local and average heat transfer results are analyzed and correlated with Reynolds number and the jet-to-plate distance. Results show that the numerical solutions matched experimental data best at low jet-to-plate distances and lower Reynolds numbers, decreasing in ability to accurately predict the heat transfer as jet-to-plate distance and Reynolds number was increased.


SPE Journal ◽  
2020 ◽  
Vol 25 (06) ◽  
pp. 2962-2983 ◽  
Author(s):  
Zhengming Xu ◽  
Xianzhi Song ◽  
Zhaopeng Zhu

Summary Viscoelastic fluids are frequently used as drilling or fracturing fluids to enhance cuttings or proppant transport efficiency. The solid transport performance of these fluids largely depends on the settling behaviors of suspended particles. Different from viscoinelastic fluids, the elastic and viscous characteristics of viscoelastic fluids both affect particle settling behaviors. In this study, to separately quantify the contribution degrees of the shear viscosity and fluid elasticity on the terminal settling velocity, we decompose the total drag force into a viscous drag force and an elastic drag force. Based on the experimental data from the available literature, it is concluded that the elastic drag force is a function of the fluid elasticity, particle diameter, particle terminal settling velocity, and density difference between the fluid and particle. The formula for the elastic drag force is determined on the basis of the force analysis, and a relationship between the elastic drag coefficient and particle Reynolds number (Re) is developed. An explicit equation that directly predicts the terminal settling velocity in viscoelastic fluids is determined by correlating the dimensionless particle diameter and Re. To validate the proposed model, a total of 108 settling experiments in viscoelastic fluids are conducted. The absolute percentage error (APE) between the predicted and measured terminal settling velocities is 15.26%, which indicates that the proposed explicit terminal settling velocity equation can provide satisfactory prediction accuracy of the terminal settling velocity for particles in viscoelastic fluids. Furthermore, an illustrative example is provided to show that the proposed model can be used to calculate the required fluid elasticity to obtain the desired terminal settling velocity when the fluid shear viscosity is fixed. The proposed models are valid with a consistency index range of approximately 0.16 to 1.2 Pa⋅sn, flow behavior index range of approximately 0.282 to 0.579, an Re range of approximately 0.005 to 30, and a fluid relaxation time range of approximately 0.183 to 110 seconds. This study can help operators choose proper drilling/fracturing fluids to enhance the cuttings/proppant transport and maximize drilling/fracturing performance.


1991 ◽  
Vol 226 ◽  
pp. 475-496 ◽  
Author(s):  
F. T. Pinho ◽  
J. H. Whitelaw

Measurements of wall pressure, and mean and r.m.s. velocities of the confined flow about a disk of 50 % area blockage have been carried out for two Newtonian fluids and four concentrations of a shear-thinning weakly elastic polymer in aqueous solution encompassing a Reynolds-number range from 220 to 138000. The flows of Newtonian and non-Newtonian fluids were found to be increasingly dependent on Reynolds numbers below 50000, with a decrease in the length of the recirculation region and dampening of the normal Reynolds stresses. At Reynolds numbers less than 25000, the recirculation bubble lengthened and all turbulence components were suppressed with increased polymer concentration so that, at a Reynolds number of 8000, the maximum values of turbulent kinetic energy were 35 and 45% lower than that for water, with 0.2% and 0.4% solutions of the polymer. Non-Newtonian effects were found to be important in regions of low local strain rates in low-Reynolds-number flows, especially inside the recirculation bubble and close to the shear layer, and are represented by both an increase in viscous diffusion and a decrease in turbulent diffusion to, respectively, 6% and 18% of the largest term of the momentum balance with a 0.4 % polymer solution at a Reynolds number of 7700. The asymmetry and unsteadiness of the flow at Reynolds numbers between 400 and 6000 is shown to be an aerodynamic effect which increases in range and amplitude with the more concentrated polymer solutions.


1994 ◽  
Vol 116 (1) ◽  
pp. 10-15 ◽  
Author(s):  
L. Jin ◽  
M. E. Chenevert

A drag force measurement method is presented which makes it possible to study the settling of particles in transparent and opaque fluids. A dimensionless treatment that takes into account the shear thinning effects of fluids was applied to normalize the measured drag force data. A wide range of particle Reynolds numbers can be covered by this method and a profile of friction factor versus Reynolds number can be established by the proposed dimensionless treatment. An algorithm for the prediction of settling of particles in non-Newtonian fluids was introduced. It can be executed by a computer program. With a good set of experimental data, the settling velocities predicted by the computer model are very close to the measured ones in the fluids tested. This method can be used to study the suspension properties of drilling and fracturing fluids, transparent or opaque. The wide coverage of Reynolds number range simplifies the experiment.


Author(s):  
Matthias Ku¨rner ◽  
Carsten Schneider ◽  
Martin G. Rose ◽  
Stephan Staudacher ◽  
Jochen Gier

The new LP turbine test rig “ATRD” at the Institute of Aircraft Propulsion Systems (ILA) at Stuttgart University has been used to study the detailed effects of Reynolds number variation. The two-stage LP turbine has been developed in a cooperation of ILA and MTU Aero Engines GmbH. Changes in the turbine characteristics are discussed. Five hole probe area traverse data has been acquired at exit from each row of aerofoils across a broad range of Reynolds numbers, over 88,000 down to 35,000. The experimental data is supported by multi-row steady CFD predictions. The behaviour of wakes, loss cores and secondary deviations is identified across the Reynolds number range. The present study is focusing on the effects of Reynolds number variation on the vane of the second stage.


2016 ◽  
Vol 139 (2) ◽  
Author(s):  
Sergio L. D. Kfuri ◽  
Edson J. Soares ◽  
Roney L. Thompson ◽  
Renato N. Siqueira

Industrial processes with non-Newtonian fluids are common in many segments such as petroleum, cosmetic, and food industries. Slurries, emulsions, and gas–liquid dispersions are some examples with industrial relevance. When a fluid flows in a pipe system, pressure losses are always present. For Newtonian fluids, a quite reasonable understanding of this phenomenon was already achieved and is available in the literature. The same cannot be stated for non-Newtonian fluids owing to their complex characteristics, such as pseudoplasticity, viscoplasticity, elasticity, and thixotropy. The understanding of the influence of these characteristics on flow behavior is very important in order to design efficient pipeline systems. The design of such systems requires the estimation of the pressure drop due to friction effects. However, there are few works regarding friction losses for non-Newtonian fluids in pipeline systems, making this task a difficult one. In this study, two classes of fluids are investigated and compared with the Newtonian results. The first category of fluids are the ones that exhibits pseudoplastic behavior and can be modeled as a power-law fluid, and the second category are the ones that possesses a yield stress and can be modeled as a Bingham fluid. Polyflow was used to compute the friction losses in both abrupt contractions and expansions laminar flow conditions. It shows that for the expansion cases, the aspect ratio affects more the local friction coefficients than for the contraction cases. The influence of the power index n on local friction losses is similar for both cases, abrupt contractions and abrupt expansions. At low Reynolds numbers, dilatant fluids present the lowest values of the friction coefficient, K, independent of geometry. At high Reynolds numbers, a reversal of the curves occurs, and the dilatant fluid presents larger values of K coefficient. For the cases investigated, there is also a Reynolds number in which all the curves exhibit the same value of K for any value of the power-law index. The effect of τy′ shows a different behavior between contractions and expansions. In the case of contractions, the material with the highest dimensionless yield stress has the highest K value. In the case of the expansions, the behavior is the opposite, i.e., the higher the yield stress, the lower is the values of the K coefficient. Equations for each accessory as a function of the rheological parameters of the fluid and the Reynolds number of the flow are also proposed. The data were adjusted according to two main equations: the two Ks method proposed by Hooper (1981, “The Two-K Method Predicts Head Losses in Pipe Fittings,” Chem. Eng., 81, pp. 96–100.) is used for all the contractions cases, and the equation proposed by Oliveira et al. (1997, “A General Correlation for the Local Coefficient in Newtonian Axisymmetric Sudden Expansions,” Int. J. Heat Fluid Flow, 19(6), pp. 655–660.) is used for all the expansions cases. The equations found were compared with the numerical results and showed satisfactory precision and thus can be used for engineering applications.


2021 ◽  
Vol 69 (3) ◽  
pp. 255-262
Author(s):  
Václav Matoušek ◽  
Andrew Chryss ◽  
Lionel Pullum

Abstract Vertical concentration distributions of solids conveyed in Newtonian fluids can be modelled using Rouse-Schmidt type distributions. Observations of solids conveyed in turbulent low Reynolds number visco-plastic carriers, suggest that solids are more readily suspended than their Newtonian counterparts, producing higher concentrations in the centre of the pipe. A Newtonian concentration profile model was adapted to include typical turbulent viscosity distributions within the pipe and particle motion calculated using non-Newtonian sheared settling. Predictions from this and the unmodified model, using the same wall viscosity, are compared with the chord averaged profile extracted from tomographic data obtained using a 50 mm horizontal pipe.


2004 ◽  
Vol 127 (2) ◽  
pp. 388-392 ◽  
Author(s):  
Ho-Yun Nam ◽  
Jong-Man Kim ◽  
Kyung-Won Seo ◽  
Seok-Ki Choi

An experimental study has been carried out to measure the pressure loss at the side orifice of a liquid metal reactor fuel assembly. The characteristics of the pressure loss at the side orifice are investigated using the experimental data measured from 17 different types of side orifices that have different geometric shapes, dimensions, and arrangements of nozzles, and a correlation that covers the whole flow range by one equation is developed. The error range of the correlation is within ±10%, and most of the errors occurred in a region where the Reynolds number is small. The range of Reynolds numbers based on the hydraulic diameter of the orifice is 2000–350,000. It is found that the geometric factor is the most important parameter for the pressure loss when the Reynolds number is >30,000. As the Reynolds number becomes smaller, its effect becomes larger, and when the Reynolds number is small, it is the most important parameter for the pressure loss at the side orifices. The measured data shows a trend that the pressure loss coefficient increases as the number of orifices increases, and the effect of the longitudinal arrangement is small.


1937 ◽  
Vol 4 (2) ◽  
pp. A53-A54
Author(s):  
W. E. Howland

Abstract The author presents a figure in which the coefficient of discharge Cd, velocity Cv, and contraction Cc determined by several investigators are plotted logarithmically as points against Reynolds’ numbers. Curves for the coefficients drawn by the author, based on theoretical considerations, show good agreement with the experimental data, thus throwing some light upon the basic phenomena of the discharge of sharp-edged orifices. The variation of the coefficient of discharge of a circular orifice as a function of the Reynolds number is explained as a purely viscous phenomenon for low Reynolds numbers, and by means of a momentum analysis for higher speeds. The analysis presented by the author leads to the development of several formulas for the discharge coefficient, which formulas are in fair agreement with experimental results.


Sign in / Sign up

Export Citation Format

Share Document