The Tangential Force Distribution on Inner Cylinder of Power Law Fluid Flowing in Eccentric Annuli with the Inner Cylinder Reciprocating Axially

Similar documents
(Online First)A Lattice Boltzmann Scheme for Diffusion Equation in Spherical Coordinate

Effects of Polymer Concentration and Molecular Weight on the Dynamics of Visco-Elasto- Capillary Breakup

A Solution of the Harry-Dym Equation Using Lattice-Boltzmannn and a Solitary Wave Methods

The Jacobsthal and Jacobsthal-Lucas Numbers via Square Roots of Matrices

Calculating the Quasi-static Pressures of Confined Explosions Considering Chemical Reactions under the Constant Entropy Assumption

The Exact Formulation of the Inverse of the Tridiagonal Matrix for Solving the 1D Poisson Equation with the Finite Difference Method

Irregular vibrations in multi-mass discrete-continuous systems torsionally deformed

A Solution of Porous Media Equation

The Order Relation and Trace Inequalities for. Hermitian Operators

THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS OF A TELESCOPIC HYDRAULIC CYLINDER SUBJECTED TO EULER S LOAD

Tensor Smooth Length for SPH Modelling of High Speed Impact

EVALUATION OF THE VISCO-ELASTIC PROPERTIES IN ASPHALT RUBBER AND CONVENTIONAL MIXES

A Numerical Study of Heat Transfer and Fluid Flow past Single Tube

Research Article Green s Theorem for Sign Data

General viscosity iterative method for a sequence of quasi-nonexpansive mappings

Air Age Equation Parameterized by Ventilation Grouped Time WU Wen-zhong

The Analysis of Convection Experiment

Similar Constructing Method for Solving the Boundary Value Problem of the Compound Kummer Equation

829. An adaptive method for inertia force identification in cantilever under moving mass

Implicit Integration Henyey Method

Kinematics of Fluids. Lecture 16. (Refer the text book CONTINUUM MECHANICS by GEORGE E. MASE, Schaum s Outlines) 17/02/2017

The Two-scale Finite Element Errors Analysis for One Class of Thermoelastic Problem in Periodic Composites

STUDY ON TWO PHASE FLOW IN MICRO CHANNEL BASED ON EXPERI- MENTS AND NUMERICAL EXAMINATIONS

One-sided finite-difference approximations suitable for use with Richardson extrapolation

Basic concept of reactive flows. Basic concept of reactive flows Combustion Mixing and reaction in high viscous fluid Application of Chaos

Asymptotics of the Solution of a Boundary Value. Problem for One-Characteristic Differential. Equation Degenerating into a Parabolic Equation

Supplemental Material: Causal Entropic Forces

Introduction to Computational Fluid Dynamics

On Graphs with Same Distance Distribution

Normally, in one phase reservoir simulation we would deal with one of the following fluid systems:

modeling of equilibrium and dynamic multi-component adsorption in a two-layered fixed bed for purification of hydrogen from methane reforming products

HEAT TRANSFER THROUGH ANNULAR COMPOSITE FINS

Uniqueness of Weak Solutions to the 3D Ginzburg- Landau Model for Superconductivity

COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD

An identification algorithm of model kinetic parameters of the interfacial layer growth in fiber composites

Visco-Rubber Elastic Model for Pressure Sensitive Adhesive

ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION

International Power, Electronics and Materials Engineering Conference (IPEMEC 2015)

CONTROLLED FLOW SIMULATION USING SPH METHOD

The Study of Teaching-learning-based Optimization Algorithm

2.29 Numerical Fluid Mechanics Fall 2011 Lecture 12

Sharp integral inequalities involving high-order partial derivatives. Journal Of Inequalities And Applications, 2008, v. 2008, article no.

A large scale tsunami run-up simulation and numerical evaluation of fluid force during tsunami by using a particle method

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

Note: Please use the actual date you accessed this material in your citation.

NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS

Grid Generation around a Cylinder by Complex Potential Functions

Thermal-Fluids I. Chapter 18 Transient heat conduction. Dr. Primal Fernando Ph: (850)

Mathematical Preparations

FORCED CONVECTION HEAT TRANSFER FROM A RECTANGULAR CYLINDER: EFFECT OF ASPECT RATIO

Numerical Heat and Mass Transfer

Back Analysis of Self-Compacting Concrete Rheological Parameters Based on H-B Model

DETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM

COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS

Comparison of the Population Variance Estimators. of 2-Parameter Exponential Distribution Based on. Multiple Criteria Decision Making Method

NONLINEAR NATURAL FREQUENCIES OF A TAPERED CANTILEVER BEAM

Homogeneous model: Horizontal pipe and horizontal well. Flow loops can't duplicate field conditions. Daniel D. Joseph. April 2001

Significance of Dirichlet Series Solution for a Boundary Value Problem

A Cartesian-grid integrated-rbf method for viscoelastic flows

Y. Guo. A. Liu, T. Liu, Q. Ma UDC

A new integrated-rbf-based domain-embedding scheme for solving fluid-flow problems

Aerodynamics. Finite Wings Lifting line theory Glauert s method

MATH 5630: Discrete Time-Space Model Hung Phan, UMass Lowell March 1, 2018

Extension of Smoothed Particle Hydrodynamics (SPH), Mathematical Background of Vortex Blob Method (VBM) and Moving Particle Semi-Implicit (MPS)

Research Article Cubic B-Spline Collocation Method for One-Dimensional Heat and Advection-Diffusion Equations

Study on Non-Linear Dynamic Characteristic of Vehicle. Suspension Rubber Component

Improved delay-dependent stability criteria for discrete-time stochastic neural networks with time-varying delays

Physics 5153 Classical Mechanics. Principle of Virtual Work-1

Solving Fractional Nonlinear Fredholm Integro-differential Equations via Hybrid of Rationalized Haar Functions

The unsteady flow characteristic research on the initial period flow of micro channel

2-π STRUCTURES ASSOCIATED TO THE LAGRANGIAN MECHANICAL SYSTEMS UDC 531.3: (045)=111. Victor Blãnuţã, Manuela Gîrţu

On Topological Indices of Circumcoronene Series of Benzenoid

Available online at ScienceDirect. Procedia Materials Science 10 (2015 )

Operating conditions of a mine fan under conditions of variable resistance

Principles of Food and Bioprocess Engineering (FS 231) Solutions to Example Problems on Heat Transfer

Effect of Structure of Internal Waves on Energy Transfer and Flow Control Between South and Middle Basins of Caspian Sea

Computing Nonequilibrium Conformational Dynamics of Structured Nucleic Acid Assemblies

The Finite Element Method

Case Study of Markov Chains Ray-Knight Compactification

Perfect Fluid Cosmological Model in the Frame Work Lyra s Manifold

MMA and GCMMA two methods for nonlinear optimization

Fast Simulation of Pyroshock Responses of a Conical Structure Using Rotation-Superposition Method

Convexity preserving interpolation by splines of arbitrary degree

Numerical Solution of Boussinesq Equations as a Model of Interfacial-wave Propagation

Numerical Solutions of a Generalized Nth Order Boundary Value Problems Using Power Series Approximation Method

Binomial transforms of the modified k-fibonacci-like sequence

Statistical inference for generalized Pareto distribution based on progressive Type-II censored data with random removals

Computational Modelling of the Unbalanced Magnetic Pull by Finite Element Method

Influence of the PCM Layer Location on the Multilayer Wall Thermal Performance

Study on thermal evolution of organic matter simulation of the four basins in the eastern of Songliao Basin

Discusion on Failure Mechanism of Piping

Lecture 13 APPROXIMATION OF SECOMD ORDER DERIVATIVES

ECONOMICS 351*-A Mid-Term Exam -- Fall Term 2000 Page 1 of 13 pages. QUEEN'S UNIVERSITY AT KINGSTON Department of Economics

POLYMER MELT FLOW IN SUDDEN EXPANSIONS: THE EFFECTS OF VISCOUS HEATING

Three views of mechanics

Consideration of 2D Unsteady Boundary Layer Over Oscillating Flat Plate

SIMULATION OF SOUND WAVE PROPAGATION IN TURBULENT FLOWS USING A LATTICE-BOLTZMANN SCHEME. Abstract

Solutions for Euler and Navier-Stokes Equations in Powers of Time

In this section is given an overview of the common elasticity models.

χ x B E (c) Figure 2.1.1: (a) a material particle in a body, (b) a place in space, (c) a configuration of the body

Transcription:

Open Journal of Flud Dynamcs, 2015, 5, 183-187 Publshed Onlne June 2015 n ScRes. http://www.scrp.org/journal/ojfd http://dx.do.org/10.4236/ojfd.2015.52020 The Tangental Force Dstrbuton on Inner Cylnder of Power Law Flud Flowng n Eccentrc Annul wth the Inner Cylnder Recprocatng Axally Junhe Ma, Yu We, Janan Zhang, Hongjan Lu, Xanjun Meng Enhanced Ol and Gas Recovery Key Laboratory of Mnstry of Educaton, Northeast Petroleum Unversty, Daqng, Chna Emal: 147974485@qq.com Receved 12 May 2015; accepted 12 June 2015; publshed 15 June 2015 Copyrght 2015 by authors and Scentfc Research Publshng Inc. Ths work s lcensed under the Creatve Commons Attrbuton Internatonal Lcense (CC BY). http://creatvecommons.org/lcenses/by/4.0/ Abstract Based on the governng equatons of the nner cynder of the unsteady flow of the power law flud n eccentrc annul wth the nner cylnder recprocatng axally n bpolar coordnate system, the calculaton formulae of tangental force were establshed, and the relevant numercal calculaton method was gven. Takng the aqueous soluton of partally hydrolyed polyacrylamdes (HPAM) for examples, the tangental forces were calculated by usng the formulae and numercal calculaton method mentoned above; the curves of the tangental force on the wall of the nner cylnder of HPAM aqueous soluton were plotted; and the effects on the tangental force of the flow behavor ndex of the power law flud, the stroke and the stroke frequency of the nner cylnder were analyed. Keywords Tangental Force, Power Law Flud, Eccentrc Annul, Recprocate Axally 1. Introducton At present n our country, there s few theoretcal research on the tangental force on the wall of the nner cylnder of the power law flud flowng n eccentrc annul wth the nner cylnder recprocatng axally. Ths force for sucker rod eccentrc wear research has an mportant nfluence. Therefore, research on the tangental force has a practcal sgnfcance. Based on the governng equatons of the nner cynder of the unsteady flow of the How to cte ths paper: Ma, J.Z., We, Y., Zhang, J.N., Lu, H.J. and Meng, X.J. (2015) The Tangental Force Dstrbuton on Inner Cylnder of Power Law Flud Flowng n Eccentrc Annul wth the Inner Cylnder Recprocatng Axally. Open Journal of Flud Dynamcs, 5, 183-187. http://dx.do.org/10.4236/ojfd.2015.52020

power law flud n eccentrc annul wth the nner cylnder recprocatng axally n the bpolar coordnate system [1], the calculaton formulae of tangental force are establshed; the relevant numercal calculaton method s gven; and takng the aqueous soluton of partally hydrolyed HPAM for examples, the tangental forces are calculated and analyed n ths artcle. 2. Mathematcal Models 2.1. Assumed Condton Power law flud wth the nner cylnder recprocatng axally, wt ( ) s axal velocty, T s the moton cycle, R o s the radus of the nner cylnder of annul, R s the radus of the nner cylnder, e s the eccentrcty of the eccentrc annulus. 2.2. Governng Equatons The governng equatons of the unsteady flow of the power law flud n eccentrc annul wth the nner cylnder recprocatng axally n the bpolar coordnate system by Haqng Cu, Tao Gao [1]-[3]. 2 η w w C w ( I2) η ( I2) P ξ ξ ζ ζ c s ρ + = + t (1a) 2.3. Intal Condtons and Boundary Condtons By assumpton flow condtons can be obtaned the boundary condtons of dynamc flow equaton. where ( ) ( ξζ,,0) ( ξζ, ) w = w p (1b) (,, ) sn ( 2π ) w ξ ζ t = W ft (1c) ( t ) w ξ, ζ, = 0 (1d) o η I 2 s the vscosty functon of the power law flud; I 2 s the second nvarant of the component of the frst order Rvln-Ercksen tensor; ξ, ζ are the dpolar coordnates, and ξ < 0, ζ 0 ; ρ s the densty of the power law flud; t s the tme; p s the pressure gradent; wp ( ξζ, ) s the velocty dstrbuton of the power law flud n eccentrc annul wth the motonless nner cylnder; W s the ampltude of the velocty of the nner cylnder recprocatng axally; f s the stroke frequency of the nner cylnder. 2.4. Tangental Force on the Wall of the Inner Cylnder From the formulae of shear stress on the wall of the nner cylnder ch ξ cos ζ τ ( I2) w ξ = η C ξ τ ch ξ cos ζ η( I ) w = C ζ ζ 2 (2) (3) τ = τ + τ (4) ξ ξ the calculaton formulae of tangental force are gven as follows whch can be expressed as where 0 ( ( ) ) π F = 2 R τ ξ, ζ, t cos θ d θ (5) π F = 2R f ( ζ ) dζ (6) 0 184

f ( ζ ) shξ chξ cosζ w chξ cosζ w = η( I2) + η( I2). ( cosζ chξ ) C ξ C ζ (7) 3. Calculaton Method From the densty of the power law flud ρ, the flow behavor ndex of the power law flud n, the consstency coeffcent of the power law flud k, the radus of the nner cylnder of annul R o, the radus of the nner cylnder R, the eccentrcty of the eccentrc annulus e, the stroke S, the stroke frequency of the nner cylnder f and the pressure gradent P, the dstrbuton of the tangental force F can be calculated: frstly, through the equatons (1a), (1b), (1c) and (1d), the velocty dstrbuton of the power law flud at anytme by usng the fnte dfference method s calculated; then, through the formula (6), the tangental force by usng the numercal ntegraton method s calculated. 4. Calculaton Examples Take HPAM aqueous soluton as an example, the radus of the nner cylnder of annulus R o = 2.960 10 2 m, the radus of the nner cylnder R = 0.885 10 2 m, the pressure gradent of the aqueous soluton P = 61.061 Pa/m and by usng the calculaton method mentoned above, the tangental forces on the wall of the nner cylnder of the HPAM aqueous soluton are calculated and the relevant dstrbuton curves are plotted as Fgure 1 to Fgure 4 (the mnus of the tangental force ndcates that the nner cylnder moves upwards). Under dfferent flow behavor ndex n, the curves of the tangental force to the cycle number N are shown as Fgure 1, where N= tt. It shows that, wth the eccentrcty, the stroke and the stroke frequency of the nner cylnder beng gven as constants, as the flow behavor ndex n ncreases, the peak values (or valley values) of the tangental force F change a lot. Under dfferent eccentrctes e, the curves of the tangental force to the cycle number are shown as Fgure 2. It shows that, wth the flow behavor ndex n, the stroke S and the stroke frequency of the nner cylnder f beng gven as constants, as the eccentrcty e ncreases, the peak values (or valley values) of the tangental force F change a lot. Under dfferent strokes S, the curves of the tangental force to the cycle number are shown as Fgure 3. It shows that, wth the flow behavor ndex n, the eccentrcty e and the stroke frequency of the nner cylnder f beng gven as constants, as the strokes S ncreases, the peak values (or valley values) of the tangental force F change a lot. Fgure 1. The curves of the dstrbuton of the tangental force of the HPAM aqueous soluton under dfferent flow behavor ndex (ρ = 998 Kg/m 3, k = 10.886 10 2 Pa s n, e = 0.965 10 2 m, S = 1.0 m, f = 0.167 H). 185

Fgure 2. The dstrbuton of the tangental force of the HPAM aqueous soluton of qualty concentraton 0.1% under dfferent eccentrctes (ρ = 998 Kg/m 3, n = 0.556, k = 10.886 10 2 Pa s n, S = 1.0 m, f = 0.167 H). Fgure 3. The dstrbuton of the tangental force of the HPAM aqueous soluton of qualty concentraton 0.1% under dfferent strokes (ρ = 998 Kg/m 3, n = 0.556, k = 10.886 10 2 Pa s n, e = 0.965 10 2 m, f = 0.167 H). Under dfferent frequences f, the curves of the tangental force dstrbuton to the cycle number are shown as Fgure 4. It shows that, wth the flow behavor ndex n, the eccentrcty e and the stroke S beng gven as constants, as the stroke frequency f ncreases, the peak values (or valley values) of the tangental force F change a lot. 5. Summary 1) Based on the governng equatons of the nner cynder of the unsteady flow of the power law flud n eccentrc annul wth the nner cylnder recprocatng axally n the bpolar coordnate system, the calculaton formulae of tangental force on the wall of the nner cylnder of the power law flud flowng n eccentrc annul wth the nner cylnder recprocatng axally n the bpolar coordnate system were establshed, and the relevant numercal calculaton method was gven. 186

Fgure 4. The dstrbuton curves of tangental force of the HPAM aqueous soluton of qualty concentraton 0.1% under dfferent frequency (ρ = 998 Kg/m 3, n = 0.556, k = 10.886 10 2 Pa s n, e = 0.965 10 2 m, S = 1.0 m). 2) Takng the HPAM aqueous soluton for examples, by usng the calculaton formulae and method mentoned above, the tangental forces of the HPAM aqueous soluton on the wall of the nner cylnder were calculated and the tangental force curves were plotted. 3) The effects on the tangental force of the flow behavor ndex of the power law flud, the eccentrcty of the eccentrc annulus, the stroke and the stroke frequency of the nner cylnder were analyed; the nfluence of the eccentrcty was more obvous. References [1] Cu, H.Q., Sun, Z. and Gao, T. (2003) Velocty Dstrbuton of Unsteady Flow of Non-Newtonan Flud n Eccentrc Annul wth the Inner Cylnder Recprocatng Axally. Journal of Hydrodynamcs, 18, 711-715. [2] Yang, Y.J., Cu, H.Q. and Gao, T. (2004) Flow Rate Dstrbuton of the Unsteady Flow of Power Law Flud n Eccentrc Annul wth Inner Cylnder Recprocatng Axally. Journal of Daqng Petroleum Insttute, 28, 17-19. [3] Meng, X.J., Lu, H.J., Zhang, H.Z., Gao, T. and Cu, H.Q. (2013) The Normal Force Dstrbuton on the Wall of the Inner Cylnder of the Power Law Flud Flowng n the Eccentrc Annul wth the Inner Cylnder Recprocatng Axally. Appled Mechancs and Materals, 487, 527-531. http://www.scentfc.net/amm.487.527 187