Data analysis from capillary rheometry can be enhanced by a method that is an alternative to the Rabinowitsch correction

Similar documents
A calibration method for a new type of rheometer

CONVERGING FLOW ON-LINE RHEOMETRY FOR AN ENGINEERING EXTENSIONAL VISCOSITY OF UPVC.

CENG 5210 Advanced Separation Processes. Reverse osmosis

Modelling of dispersed, multicomponent, multiphase flows in resource industries Section 4: Non-Newtonian fluids and rheometry (PART 1)

The Flow Characteristics of Fly Ash Slurry for Plugging Abandoned Wells Using Coiled Tubing

Fluid temperature control in rotational rheometers with plate-plate measuring systems

Abstract The large extensional viscosity of dilute polymer solutions has been shown to

Computation of turbulent natural convection at vertical walls using new wall functions

Effect of Molecular Weight on Viscosity of Polypropylene Melt by Capillary Rheometry

Why do Golf Balls have Dimples on Their Surfaces?

Rheology in Feed Production

PET467E-Analysis of Well Pressure Tests 2008 Spring/İTÜ HW No. 5 Solutions

3-D FINITE ELEMENT MODELING OF THE REMOTE FIELD EDDY CURRENT EFFECT

A New Method for Calculating Oil-Water Relative Permeabilities with Consideration of Capillary Pressure

How to measure the shear viscosity properly?

Rheology of strongly sedimenting magnetite suspensions

Corrections to flow data in polymer melts

CHARACTERIZATION OF ULTRASONIC IMMERSION TRANSDUCERS

MODIFICATION OF THE DYKSTRA-PARSONS METHOD TO INCORPORATE BUCKLEY-LEVERETT DISPLACEMENT THEORY FOR WATERFLOODS. A Thesis RUSTAM RAUF GASIMOV

CALCULATION OF STEAM AND WATER RELATIVE PERMEABILITIES USING FIELD PRODUCTION DATA, WITH LABORATORY VERIFICATION

Hydroplaning and submarine debris flows

arxiv: v3 [physics.bio-ph] 2 May 2018

Numerical Computation of Mixed Convection Past a Heated Vertical Plate within a Saturated Porous Medium with Variable Permeability

CHAPTER 3. CONVENTIONAL RHEOMETRY: STATE-OF-THE-ART. briefly introduces conventional rheometers. In sections 3.2 and 3.

Nondestructive Monitoring of Setting and Hardening of Portland Cement Mortar with Sonic Methods

A turbulence closure based on the maximum entropy method

Critical Conditions for Water-based Suppression of Plastic Pool Fires. H. Li 1, A. S. Rangwala 1 and J.L. Torero 2

A ROBUST BEAMFORMER BASED ON WEIGHTED SPARSE CONSTRAINT

On Relationship between PVT and Rheological Measurements of Polymer Melts

Shape tilting in viscoelastic fluids

Enhancing Generalization Capability of SVM Classifiers with Feature Weight Adjustment

Rheological Measurements of Cementitious Suspensions Using a Grooved Measuring Device

Finite difference solution of the mixed convection flow of MHD micropolar fluid past a moving surface with radiation effect

Novel method for on-line rheology measurement in manufacturing of non- Newtonian liquids

MASSACHUSETTS INSTITUTE OF TECHNOLOGY SPRING 2007

Introduction To Resonant. Circuits. Resonance in series & parallel RLC circuits

Group-invariant solutions of nonlinear elastodynamic problems of plates and shells *

y(x) = x w + ε(x), (1)

EFFECT OF TYPICAL MELT TEMPERATURE NON-UNIFORMITY ON FLOW DISTRIBUTION IN FLAT DIES

Science Degree in Statistics at Iowa State College, Ames, Iowa in INTRODUCTION

Heat Transfer to Sub- and Supercritical Water at Low Mass Fluxes: Numerical Analysis and Experimental Validation

Contraction flow measurements of extensional properties

Vapor Pressure Prediction for Stacked-Chip Packages in Reflow by Convection-Diffusion Model

Ithamar R. Rangel. Roney L. Thompson

In-line ultrasound based rheology of concentrated suspensions

Available online at ScienceDirect. Procedia Engineering 90 (2014 )

5. MODELING OF NON-STRATIFIED MIXTURE FLOWS (Pseudo-homogeneous flows)

Variable Viscosity Effect on Heat Transfer over a. Continuous Moving Surface with Variable Internal. Heat Generation in Micropolar Fluids

Rheological Properties of ABS at Low Shear Rates: Effects of Phase Heterogeneity

HT FACTOR ANALYSIS FOR FORCED AND MIXED CONVECTION LAMINAR HEAT TRANSFER IN A HORIZONTAL TUBE USING ARTIFICIAL NEURAL NETWORK

, Sathyamangalam, 2 Department of Mathematics, Institute of Road and Transport, , Erode

RH7-D & RH10-D CAPILLARY RHEOMETERS

A Generalization of a result of Catlin: 2-factors in line graphs

MEASUREMENTS OF TIME-SPACE DISTRIBUTION OF CONVECTIVE HEAT TRANSFER TO AIR USING A THIN CONDUCTIVE-FILM

CM4655 Polymer Rheology Lab. Torsional Shear Flow: Parallel-plate and Cone-and-plate

Effect of Viscous dissipation on Heat Transfer of Magneto- Williamson Nanofluid

File Name: Supplementary Information Description: Supplementary Figures and Supplementary Note. File Name: Peer Review File Description:

A RIGOROUS COMPRESSIBLE STREAMLINE FORMULATION FOR BLACK OIL AND COMPOSITIONAL SIMULATION

Concept of an Integrated Workflow for Geothermal Exploration in Hot Sedimentary Aquifers

Chapter 3 Non-Newtonian fluid

Further validation to the variational method to obtain flow relations for generalized Newtonian fluids

FLUID MECHANICS PROF. DR. METİN GÜNER COMPILER

Figure 3: Problem 7. (a) 0.9 m (b) 1.8 m (c) 2.7 m (d) 3.6 m

Wall-Slip of Highly Filled Powder Injection Molding Compounds: Effect of Flow Channel Geometry and Roughness

Numerical Treatment of Two-Phase Flow in Porous Media Including Specific Interfacial Area

A Study on Contact Force Between Wraps of Scroll Compressor for CO2 Refrigerant

2D TIME AVERAGED FLOW MAPPING OF DIE ENTRY IN FLOW OF HIGHLY CONCENTRATED SHEAR-THINNING AND SHEAR-THICKENING SUSPENSIONS

EARTH PRESSURES AGAINST RIGID RETAINING WALLS IN THE MODE OF ROTATION ABOUT BASE BY APPLYING ZERO-EXTENSION LINE THEORY

HIGH PRESSURE CAPILLARY RHEOMETRY ON WOOD PLASTIC COMPOSITES WITH VARIATION OF WOOD CONTENT AND MATRIX POLYMER ABSTRACT

UNCERTAINTY SCOPE OF THE FORCE CALIBRATION MACHINES. A. Sawla Physikalisch-Technische Bundesanstalt Bundesallee 100, D Braunschweig, Germany

Analytical and Experimental Validation of Electromagnetic Simulations Using COMSOL, re Inductance, Induction Heating and Magnetic Fields

TWO-DIMENSIONAL SIMULATIONS OF THE EFFECT OF THE RESERVOIR REGION ON THE PRESSURE OSCILLATIONS OBSERVED IN THE STICK-SLIP INSTABILITY REGIME

MATRIX-FRACTURE TRANSFER FUNCTIONS FOR PARTIALLY AND COMPLETELY IMMERSED FRACTURES

Analysis of fluid induced vibration of cryogenic pipes in consideration of the cooling effect

Modeling the Time-Dependent Rheological Properties of Pistachio Butter

H. W. Müllner (Sp), J. Eberhardsteiner, Technische Universität Wien (A); W. Fidi, Semperit Technische Produkte Ges.m.b.H. & Co. KG, Wimpassing (A)

Prompt Alpha Calculation With Monte Carlo Code JMCT

Polymer Engineering (MM3POE)

Polynomial and Rational Functions

Characteristic Variables and Entrainment in 3-D Density Currents

Hydraulics. B.E. (Civil), Year/Part: II/II. Tutorial solutions: Pipe flow. Tutorial 1

Early & Quick COSMIC-FFP Analysis using Analytic Hierarchy Process

Polynomial and Rational Functions

Phase of acoustic impedance and performance of standing wave thermoacoustic coolers

Fluid Flow. Fundamentals of Rheology. Rheology is the science of deformation and flow. Food rheology is the material science of food

Theoretical Design and Analysis of Gravity Assisted Heat Pipes

Rheometry. II.1 Introduction

SPRING 2011 Phys 450 Solution set 2

Design and Modeling of Fluid Power Systems ME 597/ABE Lecture 7

Hydrate Inhibition with Methanol A Review and New Concerns over Experimental Data Presentation

Performance evaluation of different model mixers by numerical simulation

Effects of Asymmetric Distributions on Roadway Luminance

Lecture Notes 8

Stat 8112 Lecture Notes Weak Convergence in Metric Spaces Charles J. Geyer January 23, Metric Spaces

Measurement and Prediction of Fluid Viscosities at High Shear Rates

USING AN ALL-SPEED ALGORITHM TO SIMULATE UNSTEADY FLOWS THROUGH MIXED HYDRO-PNEUMATIC SYSTEMS

Advances in Water Resources

Lecture #05. Methods for local Flow Velocity Measurements. AerE 545 class notes

Energy and Power Engineering, 2009, doi: /epe Published Online August 2009 (

Minimize Cost of Materials

Transcription:

ANNUAL TRANSACTIONS OF THE NORDIC RHEOLOGY SOCIETY, VOL. 18, 2010 Data analysis from capillary rheometry can be enhanced by a method that is an alternative to the Rabinoitsch correction Carlos Salas-Bringas1, Odd-Ivar Lekang1, and Reidar Barfod Schüller2 1Dep. Of Mathematical Sciences and Technology, Noregian University of Life Sciences, P.O. Box 5003, N-1432 Ås, Noray. 2Dep. Of Chemistry, Biotechnology and Food Science, Noregian University of Life Sciences, P.O. Box 5003, N-1432 Ås, Noray. ABSTRACT The folloing article compares a method (presented by Salas-Bringas et al 1 ) using a flo rate equation ith the Rabinoitsch method to estimate the consistency and flo behavior index of a Non-Netonian fluid in capillary flo. The differences is that the method presented by Salas-Bringas et al 1 gives the uncertainty of the measurements and requires less number of experiments hen different flo rates are measured. INTRODUCTION The simplicity of capillary rheometers 2, 3 has resulted in its ide usage ithin the process industries for applications as diverse as iron or slurries, foams, distilled ater, glycerol and polymer solutions 3. Due to their inherent similarity to many process flos, hich typically involve pipes, capillary rheometers are idely employed in process engineering applications 2. Measurements of pressure drop in the capillary are required to estimate absolute values of viscosity. Pressure estimations are currently performed by measuring the pressure in the reservoir or barrel above the capillary. Pressure losses occur at the capillary entry, and thus corrections must be taken to avoid rong estimations of shear stresses 1, 4, 5. Salas-Bringas et al 1 is developing a capillary rheometer that uses direct pressure measurements in the capillary. The increased developments in sensor technology ill soon allo commercial capillary rheometers to account ith pressure measurements in the capillary, and thus pressure corrections ill no longer be required. Today, the most important method to estimate the shear rate at the capillary alls. ( γ ) for Non-Netonian fluid behavior is the often termed Weissenberg - Rabinoitsch or Rabinoitsch - Mooney correction 2, 4. Salas-Bringas et al 1 presented an alternative method based in the flo rate information that does not uses the Rabinoitsch correction to estimate the consistency (K) and the flo behavior index (n) of Non-Netonian fluids. Ho the method presented by Salas- Bringas et al 1 performs hen comparing to the Rabinoitsch correction is a question that has not been analysed yet and thus, there is no reference from here to make the decision of selecting one method over the other. This article presents a study that analyses and compares the performance of both methods using the same capillary data.

Table 1. Capillary data for dough (34.7%) at room temperature made from defatted soy flour treated to cause protein denaturation. Using a capillary of D = 0.00318 m (data from Morgan 6 and used by Steffe 5, the first and the last to columns here calculated for this study). * Measured pressure drop ** Entrance loss pressure correction *** Corrected pressure drop MATERIALS AND METHODS To compare ho both methods performs to determine K and n values in a fluid ithout yield stress, K and n ill be estimated using the same set of capillary data hich is given in literature 5, 6 (Table 1), here: Q is the volumetric flo rate (m 3 /s), R the capillary radius (m), L the capillary length (m), D the capillary diameter (m), P the pressure drop (MPa) and τ the shear stress at the capillary alls (kpa). Each method ill be presented describing step by step the solving procedure. Consistency and flo behavior index by Rabinoitsch-Mooney equation The Rabinoitsch-Mooney equation is used to determine the shear rate at the capillary alls as follos:. 3Q γ = = f ( τ ) 3 π R + d ( Q/( π R σ dτ 3 )) (1) The experimental data shon in Table 1 is plotted in Fig. 1. Figure 1. Capillary data for soy dough corrected for pressure loss. The derivative of the equation created from the poer la curve fitting (Fig. 1) is expressed as follos:

d ( Q /( π R dτ 3 )) 5 2.44 = 1.06 10 ( τ ) (2)) This equation can be incorporated into the Rabinoitsch-Mooney equation (Eq. (1)) to find out the shear rate at the capillary alls:. γ 5 2.44 [ 1.06 10 ( τ ) ] 3Q = τ 3 π R + (3) Once these calculations are performed, a shear rate versus shear stress plot can be made. Consistency and flo behavior index by a flo rate equation To calculate the flo rate (Q) in capillaries, for fluids ithout yield stress that satisfy the poer la model, the folloing equation can be used 1, 2, 5 : performed if more values (different Q) are taken. The number of crosses or intersections (J N ) can be calculated by: I = 1 N i= 1 J i (5) here I is the number of curves at different Q. Averages of K and n can be estimated ith their errors around a mean. Since three repetitions ere done by Morgan 6 for a single Q, the crosses cannot be plotted in one figure as the plot demands to be built using different Q. This ill result in a large number of plots (27) that cannot be shon in this article. Hoever one of the plots is included (Fig. 2) as an example. 1 n P n (3n+ 1) n Q = π R 2l K (4) 3n + 1 Since Q, R, l P and are knon from Table 1, it is possible to estimate K and n values by solving Eq. (4). To obtain the to unknon variables (K, n), at least to equations are needed. The data from Table 1 includes four different Q, each of them repeated three times. K and n can be estimated by re-arranging Eq. (4) as: Q 1 -f(r, P, l, K, n) = 0 and iteratively find the K values for a number of n values. These values can be plotted positioning n in the independent axis and K on the dependent axis, the result ill be a curve. By doing the same procedure for a ne Q 2 and plotting the ne curve together ith the one generated using Q 1, it ill produce an intercept in a point (n, K) indicating the to unknon values. Statistical analysis (e.g. standard deviations, standard errors, etc) can be Figure 2. Six intersections for four different flo rates. The intersections are indicated ith an arro. RESULTS Consistency and flo behavior index by Rabinoitsch-Mooney equation From the plot presented in Fig. 3, is possible to obtain the consistency K, and flo behavior index n.

Figure 3. Rheogram for soy dough shoing the K (23.3 Pa s n ) and n values (0.29). As shon in Fig. 3, the Rabinoitsch- Mooney equation resulted in a K value of 23.3 Pa s n and a n value of 0.29. Consistency and flo behavior index by a flo rate equation Fig. 4 indicates the distribution of intersections for the different K and n values for all data (four flo rates repeated three times). The average of these points ill represent the K and n values of the soy dough. All intersections together (Fig. 4) presented the same trend as the one shon in Fig. 3 for the six intersections. Figure 5. Frequency plot for the distribution of n values. Figure 6. Frequency plot for the distribution of K values. The averages of K and n values ith their standard deviations are: K = 24.28 ± 9.56 and n = 0.29 ± 0.09. Figure 4. Distribution of intersections (K and n) for all plots. Each point represents one intersection. The data points or crosses displayed in Fig. 4 here re-plotted in a frequency plot to observe ho the data ere distributed for the K and n values (see Fig.5 and Fig. 6). DISCUSSIONS Similar assumptions have to be made hen using the method using the flo rate information 1 compared to the Rabinoitsch- Mooney correction 1, 2, 5 : flo is laminar and steady, fluid is incompressible, properties are time and pressure independent, temperature is constant, no slip occur at the all of the capillary, radial and tangential velocity components are zero and the pressure drop should be linear in the capillary. Therefore both methods demand the same experimental conditions. The method from Rabinoitsch demanded less plots using Morgan s data 6,

because each experiment as repeated three times. The Rabinoitsch method delivers only one n and K value ithout informing the absolute uncertainty of the measurements. Instead the method using the flo rate equation demanded more plots for these specific data, but is able to give absolute values of the uncertainty of the measurements and thus comparisons ith more confidence is possible hen using different fluids. If uncertainty is to be estimated from the Rabinoitsch method, additional statistical analyses have to be made like the Root Mean Square Error of Prediction (RMSEP). This ill result in a prolonged calculation. With Morgan s data and using RMSEP, the uncertainty can be estimated based on 12 data points (see Fig. 3). Instead, hen using the method based on the flo rate equation, the uncertainty (e.g. standard deviation, standard error, etc) can be estimated based on 42 data points (see Fig. 4). The large number of plots using the flo rate equation as due to Morgan 6 repeated three times his experiments at each flo rate. Instead, a minimum of to nonrepeated flo rates are required to obtain a single n and K value using the method based on the flo rate equation, but the uncertainty of the measurement can be done using three or more flo rates. Therefore, to obtain uncertainty is preferable to do three or more experiments at different flo rates ithout been repeated. Both methods gave the same n values (0.29), but not the same K values. Rabinoitsch method resulted in a K of 23.3 Pa s n, and the flo rate method in a K of 24.28 Pa s n. According to the estimated uncertainty of measurements, the consistency of the fluid can be found in an interval 14.72 K 33.84 hich cover the value obtained by the Rabinoitsch method. The flo behavior index of the fluid can be found in an interval 0.2 n 0.38. CONCLUSIONS Using different flo rates, the method based on the flo rate equation requires a minimum of to flo rates to obtain the n and K values hich is less than the minimum necessary to build a curve hen using the Rabinoitsch method. For the set of data given by Morgan 6 that uses a flo rate that is repeated three times, the Rabinoitsch method is faster to perform because it requires only one plot. Rabinoitsch does not give information about the uncertainty of the measurements. If uncertainty is anted, other statistical methods have to be made in additionally (e.g. RMSEP). Using the same number of experiments at different flo rates, the method based on the flo rate equation gives the uncertainty based on more data points compared to the Rabinoitsch method. The method based on the flo rate equation gives both the K and n values together ith the uncertainty in the same calculative step. The Rabinoitsch method and the method using the flo rate equation resulted in the same n value and a similar K value. REFERENCES 1. Salas-Bringas, C., O.-I. Lekang, E.O. Rukke, and R.B. Schüller, (2009), "Development of a ne capillary rheometer that uses direct pressure measurements in the capillary", Annual Transactions of the Nordic Rheology Society, 17: p. 39-47. 2. Chhabra, R.P. and J.F. Richardson, (1999) "Non-Netonian flo in the process industries". Oxford: Butterorth-Heinemann. xiii, 436 s., 3. Roberts, I., In-line and on-line rheology measurement, in Instrumentation and sensors for the food industry, E. Kress-Rogers and C.J.B. Brimelo, Editors. 2001, CRC Press: Boca Raton, Fla. p. 403-422. 1-59124-341-6

4. Rauendaal, C., (2001) "Polymer extrusion". Krugzell, Germany: Hanser Gardner Publications, Inc. XIV, 777., 1-56990-321-2 5. Steffe, J.F., (1996) "Rheological methods in food process engineering". East Lansing, Mich.: Freeman Press. XIII, 418., 0-9632036- 1-4 6. Morgan, R.G., Modeling the effects of temperature-time history, temperature, shear rate and moisture on the viscosity of defatted soy flour. Thesis dissertation. 1979, Texas A&M University, College Station: TX.