scholarly journals Response of Viscoelastic Turbulent Pipeflow Past Square Bar Roughness: The Effect on Mean Flow

Computation ◽  
2021 ◽  
Vol 9 (8) ◽  
pp. 85
Author(s):  
Shubham Goswami ◽  
Arman Hemmati

The influence of viscoelastic polymer additives on response and recovery of turbulent pipeflow over square bar roughness elements was examined using Direct Numerical Simulations at a Reynolds number of 5×103. Two different bar heights for the square bar roughness elements were examined, h/D=0.05 and 0.1. A Finitely Extensible Non-linear Elastic-Peterlin (FENE-P) rheological model was employed for modeling viscoelastic fluid features. The rheological parameters for the simulation corresponded to a high concentration polymer of 160 ppm. Recirculation regions formed behind the bar elements by the viscoelastic fluid were shorter than those associated with Newtonian fluid, which was attributed to mixed effects of viscous and elastic forces due to the added polymers. The recovery of the mean viscoelastic flow was faster. The pressure losses on the surface of the roughness were larger compared to the Newtonian fluid, and the overall contribution to local drag was reduced due to viscoelastic effects.

Author(s):  
Y. Kagawa ◽  
B. Yu ◽  
Y. Kawaguchi

For the purpose of elucidating the mechanism of drag reduction by additives and finding a way to judge optimum drag-reducing additives through a simple rheological test, we performed DNS analysis of viscoelastic fluid turbulent flow in a two-dimensional channel. In this calculation, we employed the Giesekus constitutive equation to model the interaction between water-soluble polymer, or the elastic micellar network structure, and solvent. We calculated the fluid flow by varying the rheological parameters of the model. We examined the turbulent kinetic energy budget and studied the “viscoelastic contribution” term in the budget equation for turbulent intensity, which is not apparent in normal Newtonian fluid turbulence. Viscoelastic contribution has a characteristic effect on viscoelastic fluid turbulence. We concluded that the viscoelastic contribution plays a major role in turbulent frictional drag reduction. Dissipation and viscoelastic contribution serve as a key factor of turbulent kinetic energy loss in most areas of the channel. From the visualization of local and instantaneous eddy behavior, we discussed the relationship between viscoelastic contribution, elastic energy and turbulent production. We found that viscoelastic contribution serves as a direct local source of turbulent production, and that energy is stored in the elasticity.


2021 ◽  
Vol 36 (3) ◽  
pp. 165-176
Author(s):  
Kirill Nikitin ◽  
Yuri Vassilevski ◽  
Ruslan Yanbarisov

Abstract This work presents a new approach to modelling of free surface non-Newtonian (viscoplastic or viscoelastic) fluid flows on dynamically adapted octree grids. The numerical model is based on the implicit formulation and the staggered location of governing variables. We verify our model by comparing simulations with experimental and numerical results known from the literature.


2006 ◽  
Author(s):  
Shriram Pillapakkam ◽  
Pushpendra Singh ◽  
Denis L. Blackmore ◽  
Nadine Aubry

A finite element code based on the level set method is developed for performing two and three dimensional direct numerical simulations (DNS) of viscoelastic two-phase flow problems. The Oldroyd-B constitutive equation is used to model the viscoelastic liquid and both transient and steady state shapes of bubbles in viscoelastic buoyancy driven flows are studied. The influence of the governing dimensionless parameters, namely the Capillary number (Ca), the Deborah Number (De) and the polymer concentration parameter c, on the deformation of the bubble is also analyzed. Our simulations demonstrate that the rise velocity oscillates before reaching a steady value. The shape of the bubble, the magnitude of velocity overshoot and the amount of damping depend mainly on the parameter c and the bubble radius. Simulations also show that there is a critical bubble volume at which there is a sharp increase in the bubble terminal velocity as the increasing bubble volume increases, similar to the behavior observed in experiments. The structure of the wake of a bubble rising in a Newtonian fluid is strikingly different from that of a bubble rising in a viscoelastic fluid. In addition to the two recirculation zones at the equator of the bubble rising in a Newtonian fluid, two more recirculation zones exist in the wake of a bubble rising in viscoelastic fluids which influence the shape of a rising bubble. Interestingly, the direction of motion of the fluid a short distance below the trailing edge of a bubble rising in a viscoelastic fluid is in the opposite direction to the direction of the motion of the bubble, thus creating a “negative wake”. In this paper, the velocity field in the wake of the bubble, the effect of the parameters on the velocity field and their influence on the shape of the bubble are also investigated.


Water ◽  
2020 ◽  
Vol 12 (5) ◽  
pp. 1341
Author(s):  
Irene Daprà ◽  
Giambattista Scarpi

The purpose of this work is to study the motion of a non-Newtonian fluid in a rock fracture, generated by a constant pressure gradient to which a pulsating component is superposed. The momentum equation is faced analytically by adopting a logarithmic constitutive law; the velocity is expressed as a power series of the amplitude of the pulsating component, up to the second order, easily usable for numerical calculations. The results obtained are compared with those provided in the past by the authors, using a three-parameter Williamson model. The comparison highlights that the value of the mean flow rate in a period differs by less than 10% even if the velocity profiles look quite different.


Author(s):  
ROSANE DA SILVA RODRIGUES ◽  
ÂNGELA MARIA GOZZO ◽  
ROBERTO HERMÍNIO MORETTI

Estudou-se o comportamento reológico de extratos elaborados com grãos, farinha integral e isolado protéico de soja. Extratos com 3% de proteína foram obtidos de grãos de soja (em equipamento conhecido como “vaca mecânica”), de farinha integral e de isolado protéico (por dissolução em água) e pasteurizados a 74±2°C por 15 seg. Os parâmetros reológicos estudados foram o coeficiente de consistência, o índice de comportamento do fluxo e a viscosidade aparente. Foram determinados, também, sólidos totais, proteínas, lipídios, cinzas, fibra bruta, carboidratos e índice de sedimentação (AOAC, 1995). Verificou-se que a viscosidade aparente, a 4 e 25°C, do extrato obtido de grãos de soja (fluido nãonewtoniano com comportamento pseudoplástico) foi maior que a do extrato da farinha (fluido não-newtoniano com comportamento pseudoplástico a 4°C e dilatante a 25°C) e do isolado (fluido newtoniano). O maior índice de sedimentação constatado no extrato de farinha integral e o baixo teor de sólidos solúveis no isolado protéico foram determinantes, entre outros fatores, no comportamento reológico verificado. AbstractRHEOLOGICAL BEHAVIOR OF SOYMILK, WHOLE SOY FLOUR AND SOY ISOLATED PROTEIN The rheological properties of extracts elaborated with soybeans, whole soy flour and isolated protein were studied. Extracts with 3% of protein were obtained from soybeans (in a equipment called “mechanic cow”), whole soy flour and isolated soy protein (by dissolution in water) and pasteurized at 74±2°C for 15 seg. The rheological parameters studied were the flow behavior index, consistency coefficient and apparent viscosity. Total solid, proteins, lipids, ashes, brute fiber, carbohydrates and sedimentation index (AOAC, 1995) was also determined. The apparent viscosity, at 4 and 25°C, of the extract obtained from soybeans (non-newtonian fluid with pseudoplastic behavior), was greater than the one of whole soy flour extract (nonnewtonian fluid with pseudoplastic behavior at 4°C and dilatant at 25°C) and of the isolated soy protein extract (Newtonian fluid). The highest sedimentation index verified in the whole soy flour extract and the low content of soluble solids in the extract of isolated soy protein defined, among others factors, the rheological behavior verified.


1989 ◽  
Vol 111 (2) ◽  
pp. 149-153 ◽  
Author(s):  
E. Logan ◽  
P. Phataraphruk

The response of a fully developed pipe flow to wall mounted roughness elements of rectangular cross section was investigated experimentally using a probe with a single hot-wire. Four heights of rectangular, ring-type elements were installed rigidly in a 63.5-mm diameter, smooth-walled, circular pipe in which air was flowing at a Reynolds number of 50,000. After passing over the roughness element, the flow recovery occurred in three stages. The three flow regions are delineated, and the velocity profiles for each are correlated.


2011 ◽  
Vol 689 ◽  
pp. 279-316 ◽  
Author(s):  
Xuesong Wu

AbstractIn two previous papers (Wu, J. Fluid Mech., vol. 453, 2002, p. 289, and Wu & Hogg, J. Fluid Mech., vol. 550, 2006, p. 307), a formal asymptotic procedure was developed to calculate the sound radiated by unsteady boundary-layer flows that are described by the triple-deck theory. That approach requires lengthy calculations, and so is now improved to construct a simpler composite theory, which retains the capacity of systematically identifying and approximating the relevant sources, but also naturally includes the effect of mean-flow refraction and more importantly the back action of the emitted sound on the source itself. The combined effect of refraction and back action is represented by an ‘impedance coefficient’, and the present analysis yields an analytical expression for this parameter, which was usually introduced on a semi-empirical basis. The expression indicates that for Mach number $M= O(1)$, the mean-flow refraction and back action of the sound have a leading-order effect on the acoustic field within the shallow angles to the streamwise directions. A parametric study suggests that the back effect of sound is actually appreciable in a sizeable portion of the acoustic field for $M\gt 0. 5$, becomes more pronounced, and eventually influences the entire acoustic field in the transonic limit. In the supersonic regime, the acoustic field is characterized by distinctive Mach-wave beams, which exert a leading-order influence on the source. The analysis also indicates that acoustic radiation in the subsonic and supersonic regimes is fundamentally different. In the subsonic regime, the sound is produced by small-wavenumber components of the hydrodynamic motion, and can be characterized by acoustic multipoles, whereas in the supersonic regime, broadband finite-wavenumber components of the hydrodynamic motion contribute and the concept of a multipolar source becomes untenable. The global acoustic feedback loop is investigated using a model consisting of two well-separated roughness elements, in which the sound wave emitted due to the scattering of a Tollmien–Schlichting (T–S) wave by the downstream roughness propagates upstream and impinges on the upstream roughness to regenerate the T–S wave. Numerical calculations suggest that at high Reynolds numbers and for moderate roughness heights, the long-range acoustic coupling may lead to global instability, which is characterized by self-sustained oscillations at discrete frequencies. The dominant peak frequency may jump from one value to another as the Reynolds number or the distance between the roughness elements is varied gradually.


2003 ◽  
Author(s):  
Ayten S. Bakhtiyarova

In this paper we present the results of the 3-D simulations of the power law model fluid flow in rectangular cavity. FLOW-3D commercial, software (Flow-Simulations, Inc.) has been used to predict flow velocity distributions and pressure losses when non-Newtonian fluid flows through the square cavity.


SPE Journal ◽  
2007 ◽  
Vol 12 (01) ◽  
pp. 100-107 ◽  
Author(s):  
Zhu Chen ◽  
Ramadan Mohammed Ahmed ◽  
Stefan Z. Miska ◽  
Nicholas E. Takach ◽  
Mengjiao Yu ◽  
...  

Summary An experimental investigation on polymer-based drilling foams was carried out. Rheology tests were performed with foams that have different concentrations of hydroxylethylcellulose (HEC) and 1% commercial surfactant. Experiments were conducted in a large-scale flow loop that permits foam flow through 2-, 3-, and 4-in. pipe sections, and a 6×3.5-in. annular section. During the experiments, frictional pressure losses across the pipe and annular sections were measured for different gas/liquid flow rates, polymer concentrations (0, 0.25, and 0.5%), and foam qualities (70, 80, and 90%). Significant rheological variations were observed between aqueous foams containing no polymers and polymer-thickened foams. Experimental data show three distinct flow curves for the 2-, 3-, and 4-in. pipe sections, which indicates the presence of wall slip. The Oldroyd-Jastrzebski approach was used to calculate the wall slip velocity and determine the true shear rate. It has been found that wall slip decreases as the foam quality or polymer concentration increases. Two foam hydraulic models, which use slip-corrected and slip-uncorrected rheological parameters, have been proposed. These models are applicable for predicting pressure loss in pipes and annuli. Model predictions for the annular test section are compared with the measured data. A satisfactory agreement between the model predictions and measured data is obtained. This paper will help to better design foam drilling and cleanup operations. Introduction The use of drilling foams is increasing because foams exhibit properties that are desirable in many drilling operations. In practice, aqueous and polymer-based foams have been used with commercial success. However, drilling-foam rheology and hydraulics are still not sufficiently understood to minimize the risk and costs associated with foam drilling. It is generally accepted that the addition of polymers to the liquid phase affects the viscosity and stability of foams. However, the degree to which the bulk properties of drilling foams are enhanced by polymers has not been well understood and is difficult to predict. For safe and economical foam drilling, accurate knowledge of bottomhole pressure is essential. However, foam rheology and pressure drop predictions are not accurate enough to provide adequate hydraulic design information such as equivalent circulation density. This problem is more pronounced when polymers are added, because the apparent foam viscosity of polymer-thickened foams can be significantly higher than aqueous foams. It becomes apparent that there is a need for polymer foam rheological characterization in order to improve the knowledge of foam rheology and hydraulics. Foam rheological characterization was carried out using large-scale, single-pass pipe viscometers (composed of 2-, 3-, and 4-in. pipe sections). Foam qualities were varied from 70 to 90%. Test pressure and temperature were 100 psig and 80°F. Two foam hydraulic models were considered, assuming both no-slip condition at the wall and slip condition at the wall. The first model assumes no-slip boundary conditions in both pipes and annulus. By assuming no slip condition at the wall, slip-uncorrected foam rheological parameters were obtained from the pipe viscometer measurements. It has been found that if we plot friction factors vs. Reynolds numbers for all test data, regardless of pipe diameters, foam qualities, and flow rates, a single curve is obtained. This curve is similar to that obtained for incompressible fluid flow. Pressure drop in the annulus is calculated with the proposed model, and satisfactory predictions are obtained. The second model is based on the assumption that there is wall slip in both pipes and annulus. Rheological parameters and wall-slip coefficient corrections were first obtained using Oldroyd-Jastrzebski approach. The annular pressure losses are predicted based on slip-corrected rheological parameters and wall-slip coefficient correlations.


Sign in / Sign up

Export Citation Format

Share Document