Modelling rheological cone-plate test conditions

Similar documents
On composite conformal mapping of an annulus to a plane with two holes

Lecture 24: Observability and Constructibility

Auchmuty High School Mathematics Department Sequences & Series Notes Teacher Version

Technical Report: Bessel Filter Analysis

MATH Midterm Solutions

CHAPTER 5 : SERIES. 5.2 The Sum of a Series Sum of Power of n Positive Integers Sum of Series of Partial Fraction Difference Method

MATH /19: problems for supervision in week 08 SOLUTIONS

Mapping Radius of Regular Function and Center of Convex Region. Duan Wenxi

a) The average (mean) of the two fractions is halfway between them: b) The answer is yes. Assume without loss of generality that p < r.

( ) 1 Comparison Functions. α is strictly increasing since ( r) ( r ) α = for any positive real number c. = 0. It is said to belong to

Chapter 2 Sampling distribution

KEY. Math 334 Midterm II Fall 2007 section 004 Instructor: Scott Glasgow

= 5! 3! 2! = 5! 3! (5 3)!. In general, the number of different groups of r items out of n items (when the order is ignored) is given by n!

Effect of Material Gradient on Stresses of Thick FGM Spherical Pressure Vessels with Exponentially-Varying Properties

Multivector Functions

Supplementary materials. Suzuki reaction: mechanistic multiplicity versus exclusive homogeneous or exclusive heterogeneous catalysis

BINOMIAL THEOREM An expression consisting of two terms, connected by + or sign is called a

BINOMIAL THEOREM NCERT An expression consisting of two terms, connected by + or sign is called a

ME 354, MECHANICS OF MATERIALS LABORATORY MECHANICAL PROPERTIES AND PERFORMANCE OF MATERIALS: TORSION TESTING*

DANIEL YAQUBI, MADJID MIRZAVAZIRI AND YASIN SAEEDNEZHAD

LESSON 15: COMPOUND INTEREST

Using Counting Techniques to Determine Probabilities

SOME ARITHMETIC PROPERTIES OF OVERPARTITION K -TUPLES

Some Properties of the K-Jacobsthal Lucas Sequence

Using Difference Equations to Generalize Results for Periodic Nested Radicals

The Application of a Maximum Likelihood Approach to an Accelerated Life Testing with an Underlying Three- Parameter Weibull Model

FAR FIELD SOLUTION OF SH-WAVE BY CIRCULAR INCLUSION AND LINEAR CRACK

Counting Functions and Subsets

Progression. CATsyllabus.com. CATsyllabus.com. Sequence & Series. Arithmetic Progression (A.P.) n th term of an A.P.

Rotational symmetry applied to boundary element computation for nuclear fusion plasma

Volcanic activity and tidal heating of Saturn s moon Enceladus

Introduction to the Theory of Inference

By the end of this section you will be able to prove the Chinese Remainder Theorem apply this theorem to solve simultaneous linear congruences

Greatest term (numerically) in the expansion of (1 + x) Method 1 Let T

Lecture 2: Stress. 1. Forces Surface Forces and Body Forces

ECEN 5014, Spring 2013 Special Topics: Active Microwave Circuits and MMICs Zoya Popovic, University of Colorado, Boulder

L8b - Laplacians in a circle

Range Symmetric Matrices in Minkowski Space

THE GRAVITATIONAL POTENTIAL OF A MULTIDIMENSIONAL SHELL

Homologous Stars: Simple Scaling Relations

THE ANALYSIS OF SOME MODELS FOR CLAIM PROCESSING IN INSURANCE COMPANIES

EVALUATION OF SUMS INVOLVING GAUSSIAN q-binomial COEFFICIENTS WITH RATIONAL WEIGHT FUNCTIONS

Generalized Fibonacci-Lucas Sequence

2012 GCE A Level H2 Maths Solution Paper Let x,

SVD ( ) Linear Algebra for. A bit of repetition. Lecture: 8. Let s try the factorization. Is there a generalization? = Q2Λ2Q (spectral theorem!

The Pigeonhole Principle 3.4 Binomial Coefficients

Strong Result for Level Crossings of Random Polynomials

Ground Rules. PC1221 Fundamentals of Physics I. Uniform Circular Motion, cont. Uniform Circular Motion (on Horizon Plane) Lectures 11 and 12

A NOTE ON DOMINATION PARAMETERS IN RANDOM GRAPHS

Generalized Near Rough Probability. in Topological Spaces

Chapter 8 Complex Numbers

EDEXCEL NATIONAL CERTIFICATE UNIT 28 FURTHER MATHEMATICS FOR TECHNICIANS OUTCOME 2- ALGEBRAIC TECHNIQUES TUTORIAL 1 - PROGRESSIONS

Department of Mathematics, IST Probability and Statistics Unit

AN ANALYTICAL SOLUTION OF THE REYNOLDS EQUATION FOR THE FINITE JOURNAL BEARING AND EVALUATION OF THE LUBRICANT PRESSURE

FINITE ELEMENT ANALYSIS OF A BWR FEED WATER DISTRIBUTOR UNDER EXTREME TRANSIENT PRESSURE LOAD

LOCUS OF THE CENTERS OF MEUSNIER SPHERES IN EUCLIDEAN 3-SPACE. 1. Introduction

Advanced Physical Geodesy

Minimization of the quadratic test function

PROGRESSION AND SERIES

At the end of this topic, students should be able to understand the meaning of finite and infinite sequences and series, and use the notation u

ANSWERS, HINTS & SOLUTIONS HALF COURSE TEST VII (Main)

physicsandmathstutor.com

Lecture 6: October 16, 2017

Strong Result for Level Crossings of Random Polynomials. Dipty Rani Dhal, Dr. P. K. Mishra. Department of Mathematics, CET, BPUT, BBSR, ODISHA, INDIA

Minimal order perfect functional observers for singular linear systems

Ch 3.4 Binomial Coefficients. Pascal's Identit y and Triangle. Chapter 3.2 & 3.4. South China University of Technology

On ARMA(1,q) models with bounded and periodically correlated solutions

FIXED POINT AND HYERS-ULAM-RASSIAS STABILITY OF A QUADRATIC FUNCTIONAL EQUATION IN BANACH SPACES

Conditional Convergence of Infinite Products

Complementary Dual Subfield Linear Codes Over Finite Fields

( ) ( ) ( ) ( ) Solved Examples. JEE Main/Boards = The total number of terms in the expansion are 8.

Apply change-of-basis formula to rewrite x as a linear combination of eigenvectors v j.

Inertia Effects in an Externally Pressurized Thrust Bearing Using Herschel - Bulkley Lubricants

Journal of Non-Newtonian Fluid Mechanics

N-BODY SOLUTIONS AND COMPUTING GALACTIC MASSES

ON EUCLID S AND EULER S PROOF THAT THE NUMBER OF PRIMES IS INFINITE AND SOME APPLICATIONS

r, this equation is graphed in figure 1.

On randomly generated non-trivially intersecting hypergraphs

Math 7409 Homework 2 Fall from which we can calculate the cycle index of the action of S 5 on pairs of vertices as

MATHS FOR ENGINEERS ALGEBRA TUTORIAL 8 MATHEMATICAL PROGRESSIONS AND SERIES

14th European Signal Processing Conference (EUSIPCO 2006), Florence, Italy, September 4-8, 2006, copyright by EURASIP

2. Characteristics of Synchrotron Radiation

Relating to, connected or concerned with, quality or qualities. Now usually in implied or expressed opposition to quantitative.

Chaoyang University of Technology -- Radial Consolidation -- APPLICATION OF TERZAGHI'S THEORY OF CONSOLIDATION TO PROBLEMS INVOLVING RADIAL FLOW

THE ANALYTIC LARGE SIEVE

12.6 Sequential LMMSE Estimation

CRACK DETECTION IN EULER-BERNOULLI BEAMS ON ELASTIC FOUNDATION USING GENETIC ALGORITHM BASED ON DISCRETE ELEMENT TECHNIQUE

The Multivariate-t distribution and the Simes Inequality. Abstract. Sarkar (1998) showed that certain positively dependent (MTP 2 ) random variables

The number of r element subsets of a set with n r elements

JJMIE Jordan Journal of Mechanical and Industrial Engineering

ICS141: Discrete Mathematics for Computer Science I

Masses and orbits of minor planets with the GAIA mission

Einstein Classes, Unit No. 102, 103, Vardhman Ring Road Plaza, Vikas Puri Extn., Outer Ring Road New Delhi , Ph. : ,

Calculation of Matrix Elements in the Foldy-Wouthuysen Representation

ON CERTAIN CLASS OF ANALYTIC FUNCTIONS

Chem 453/544 Fall /08/03. Exam #1 Solutions

VIII.3 Method of Separation of Variables. Transient Initial-Boundary Value Problems

Journal of Molecular Graphics and Modelling

Core Variation in the Entrance Region Flow of Herschel- Bulkley Fluid in an Annuli

ELEMENTARY AND COMPOUND EVENTS PROBABILITY

Transcription:

ANNUAL TRANSACTIONS OF THE NORDIC RHEOLOGY SOCIETY, VOL. 16, 28 Modellig heological coe-plate test coditios Reida Bafod Schülle 1 ad Calos Salas-Bigas 2 1 Depatmet of Chemisty, Biotechology ad Food Sciece, Nowegia Uivesity of Life Scieces, N-1432 Ås, Noway. 2 Depatmet of Mathematical Scieces ad Techology, Nowegia Uivesity of Life Scieces, N-1432 Ås, Noway. ABSTRACT Cotol with flow coditios is of paamout impotace i elatio to heological tests. Special focus is made o the shea ate distibutio i the flow volume of coe/plate test geometies ad o the shea stess depedece o the cetal gap. Small vaiatios i the cetal gap due to tempeatue vaiatios may cause sigificat measuemet eos especially fo dilatat fluids. INTRODUCTION It is ofte desied to test o-newtoia fluids at coditios of costat shea ate sice the viscosity fo such fluids may be highly depedet o the shea ate. Coe/plate test geometies ae ofte chose fo these types of tests sice this test geomety appoaches the costat shea ate coditio i the whole flow domai. The flow coditios i coe/plate geometies have peviously bee descibed by Waltes 1. Complex flow coditios ca be peset if the coe agle is lage the 4º. Howeve, as will be show i this study, the shea ate is ot costat, ad the effects of this ca be sigificat. The cetal gap cleaace is ofte set at the begiig of a test whe the tempeatue equals oom tempeatue. If the test tempeatue is much highe o iceases duig the test it has bee show 2 that the coe shaft tempeatue also iceases ad this may cause the cetal gap distace to chage. The peset study assumes that the iitial gap is set to 5 μm, ad the effects of gap cleaace vaiatios ae quatified. BASIC THEORY Rheological behaviou fo fluids with o yield stess is ofte descibed by the well kow powe law equatio: τ = Kγ& (1) whee τ is the shea stess ad γ& is the shea ate. K is called the cosistecy coefficiet ad is the powe law idex 3. Diffeetiatio of the equatio with espect to shea ate yields a expessio fo the dyamic viscosity: dτ η = = K & γ d & γ 1 (2) Fo o-newtoia fluids, whee 1, it is clealy see that the viscosity is a fuctio of shea ate. MATERIALS AND METHODS The geomety studied i this wok is a well kow coe/plate system (Paa Physica MK22) whee the coe agle, α, is 1º ad the adius is 25 mm. The default cete distace fom the bottom plate duig measuemets is.5 mm. The zeo gap

settig is omally doe duig istumet iitiatio at oom tempeatue. It is ofte stated the viscosity measuemets i coe/plate systems ae pefomed with a costat shea ate eveywhee i the test volume. This is, howeve, ot quite tue. Cosideig the geomety i Fig. 1 we ca wite the followig equatio fo the agula plate velocity, u, ad fo the gap, δ : u = ω (3) δ = δ + taα (4) I = τ 2πd = 2π K & γ d (6) = 2πK ω d δ + taα Aalytical solutio to this itegal is ot kow, but it is possible to use fiite elemet modellig to calculate the flow domai. Howeve, if we let the gap distace ted to zeo, the ideal case with close to costat shea ate, we obtai the simplified appoximatio: I = 2πK = 2πK ω d taα ω d taα ω = πk taα 2 (7) Figue 1: Sketch (ot to scale) of coe plate test geomety. Radius is 25 mm, coe agle is 1 ad the cetal gap is 5 µm. The shea ate ca the be expessed by: u ω & γ = = (5) δ δ + taα The shea ate calculated fom this equatio fo diffeet values of δ is show i Fig. 2 with a otatioal speed of 1 ad/s. The toque o the shaft ca be expessed as the itegal of the shea stess multiplied by the adius ove the coe suface, thus: This equatio ca be solved eadily givig the equied combiatios of K ad fo simila values of toque o the divig shaft, I. The adius,, is 25 mm. The solutio of Eq. 7 is show i Fig. 3. K K 2 1 ω 1 2 = ta α (8) Eq. 8 has bee used to detemie the values of K to give the same toque fo a pseudo plastic ad a dilatat fluid; see Table 1. Table 1: Values fo compaative calculatios givig a toque of.56 Nm with zeo gap cleaace. K Newtoia case 5. 1. Pseudo plastic case 37.85.5 Dilatat case.87 2.

The values fo the Newtoia case have bee abitaily chose. The shea stess distibutio o the coe suface ca be detemied by combiatio of Eq. 1 ad Eq. 5, thus: ω τ = K δ + taα Shea Rate (1/s) 6 5 4 3 2 1 Coe/Plate Geomety, 1 coe agle Rotatioal speed = 1 ad/s Nomal opeatio (9) Cete cleaace (mm).1.5.1.15.2.25 5 1 15 2 25 Radius (mm) Figue 2: Calculated shea ate distibutio i a MK22 coe/plate geomety fo diffeet cetal gap cleaaces. FINITE ELEMENT MODELLING Fiite elemet modellig has peviously bee used to quatify the lack of tempeatue cotol i plate/plate heomete tests 2. The calculatios wee i this study pefomed with the compute pogam COMSOL Multiphysics 3.4. The basic geomety used i the calculatio was the MK22 set up with a.5 mm cetal gap cleaace. All the calculatios wee made with oe otatioal speed, ω = 1 ad/s. Both the coe ad the bottom plate bouday coditios wee specified as o-slip while the fluid to ai bouday coditio was specified as slip. The coe was specified as a movig wall with tagetial velocity, ω.the mai focus was ot effects of vaiatio i Reyolds umbe, Re, beig i the lamia egio with low values of Re. The umbe of degees of feedom solved fo was typically 5 ad the calculatio time was appoximately 1 miute. A figue of the axis-symmetic flow domai is show i Fig. 4. Toque (Nm) 5 1 1 1.1.1.1 Pseudoplastic Case Newtoia Case K =.87 K = 37.85 K = 5 Dilatat Case.2.2.4.6.8 1 1.2 1.4 1.6 1.8 2 Figue 3: The solutio to Eq. 7 ad the cases descibed i fo a commo toque of.56 Nm. The aim of the fiite elemet modellig was to see if vaiatios i K ad affected the flow ad to ivestigate how the shea stess distibutio o the coe suface was affected by itoducig a gap cleaace fo diffeet types of fluids as specified i Table 1. RESULTS The esults show that thee is o lage effect of vaiatios i the cosistecy coefficiet ad the powe law idex o the velocity distibutio ad theefoe oly a distibutio plot fo the Newtoia case is show i Fig. 5 fo the egio ea the peiphey. Calculatios wee made fo a gap cleaace of.5 mm ad fo a gap cleaace equal to zeo. The esults fo zeo gap show that the shea stess was the same at all adii ad the same fo pseudoplastic, Newtoia ad dilatat flow (Fig. 6).

Figue 4: Nomal filled geomety (to scale) of axis-symmetic 2D fluid domai with cetelie o the left. The coe adius is 25 mm ad thee is a.5 mm cetal gap cleaace. Howeve, whe a gap cleaace was itoduced, the shea stess distibutios became those show i Fig. 6. The sesitivity to vaiatios i the cetal gap cleaace fo the dilatat test fluid is show i Fig.7. Figue 5: Newtoia iso-velocity lies ea the peiphey of the coe. 3 27 24 DISCUSSION The legth of the otatig shaft is appoximately 1 mm, ad the pat of this that ca be affected by chages i tempeatue is 5 mm, say. The coefficiet of liea themal expasio of stailess steel at 2 ºC is 17.3 1-6 K -1. Calculatios show that the chages i the cetal gap due to a tempeatue chage of 5 K is 43 μm. The cetal gap cleaace must theefoe ot be set too small to pevet mechaical cotact, ad the stadad value is 5 μm. We see that quite small chages i tempeatue ca give appeciable chage i gap size. The mai focus of this aticle is to make a evaluatio of the shea ate distibutios foud i omal coe/plate heological measuemets. It is omally accepted o stated that the shea ate is costat i the measuig volume if a coe/plate geomety is chose. Calculatios show that this is ot the case (Fig. 2) ad the fiite elemet calculatios pefomed with COMSOL Multiphysics give the same esult. 21 Shea stess (Pa) 18 15 12 9 6 3 Newtoia case Pseudoplastic case Dilatat case Pseudoplastic case Newtoia case Dilatat case Zeo gap coditio.5.1.15.2.25 Radius (m) Figue 6: Wall shea stess distibutio o the coe as a fuctio of adius fo the thee basic cases whe a.5 mm gap is itoduced. Shea stess (Pa) 3 27 24 21 18 15 12 9 6 3 1 μm 2 μm 1 μm 25 μm 5 μm.5.1.15.2.25 Radius (m) Figue 7: Sesitivity to vaiatio i cetal gap cleaace fo the dilatat fluid of Table 1.

Reductio i wall shea stess (%) -2-4 -6-8 -1-12 -14-16 -18 Dilatat case Newtoia case Pseudoplastic case -2.1.2.3.4.5 Cetal gap (mm) Figue 8: Reductio i shea stess magitude at = 25 mm as a fuctio of cetal gap cleaace fo the thee base cases i Table 1. We see that the shea ate geeally iceases with adius fo all the cases whee the gap is lage tha zeo. The egio of maximum shea ate is located at the coe peiphey ad the shea ate is zeo at the cete. The esults fom this study show clealy that coe/plate test geometies do ot esue costat shea ate i the test volume. This would oly be tue if the cetal gap cleaace was zeo, but i pactise this gap is sigificat to make sue that themal effects o shaft elogatio do ot cause mechaical cotact betwee the coe ad the plate. Whe a cetal gap cleaace is itoduced, the shea stesses o the coe ae geeally educed. This, of couse, is a esult of the educed shea ate esultig fom iceasig the distace betwee the coe ad the bottom plate at all adii. Howeve, it is obseved lage diffeeces betwee the pseudoplastic, Newtoia ad dilatat fluids. The lagest eductio i shea stess is see with the dilatat fluid. The sesitivity to chages i the cetal gap is illustated i Fig. 7 fo a dilatat fluid beig most sesitive to chages i the cetal gap. The chose efeece coditio is fo zeo gap size sice the shea stess is the same at all adii i this special case. The effects would, howeve, be simila if the 5 μm gap was used as the efeece positio. Eve small chages i the gap cleaace cause sigificat chage i coe wall shea stess. Themal expasio effects may cause sigificat vaiatios i the cetal gap, so it seems vital that this is cotolled i ode to miimize measuemet eos. Fig. 8 shows the pecetage eductio i shea stess level as fuctio of cetal gap cleaace fom Eq. 9. The shea stesses closest to the peiphey give the lagest cotibutio to the toque, but it is see that the shea stesses ae also sigificatly affected i this egio whe the gap cleaace is vaied. CONCLUSIONS The coclusios fom this wok ca be summaized as follows: The shea ate is ot costat i omal coe/plate test geometies whe a cetal gap cleaace exists. It is zeo at the cete ad the maximum value is at the peiphey. The cetal gap cleaace ca be affected by a chage i tempeatue of the lowe coe shaft causig measuemet eos. The calculatios show that cotol with the cetal gap cleaace is vey impotat cotollig the shea stess distibutio o the coe. Dilatat fluids ae moe sesitive to vaiatios i the cetal gap cleaace tha Newtoia ad pseudoplastic fluids.

REFERENCES 1. Waltes, K., (1975) "Rheomety". Lodo: Chapma & Hall ISBN: 412 129 9 2. Schülle, R.B., R. O, ad C. Salas- Bigas, (27), "Fiite elemet modellig of the fluid tempeatue i a plate plate otatioal heomete i oscillatoy tests", Aual Tasactios of the Nodic Rheology Society, 15: p. 261-265. 3. Steffe, J.F., (1996) "Rheological methods i food pocess egieeig (Secod editio)". East Lacig: Feema Pess ISBN: -963236-1-4