RP 2.4. SEG/Houston 2005 Annual Meeting 1513

Size: px
Start display at page:

Download "RP 2.4. SEG/Houston 2005 Annual Meeting 1513"

Transcription

1 P 2.4 Measurement of sear ave velocity of eavy oil De-ua Han, Jiajin Liu, University of Houston Micael Batzle, Colorado Scool of Mines Introduction It is ell knon tat te fluids ave no sear modulus and terefore no sear ave can propagate troug fluids. But eavy oils ave properties tat are muc complex tan ligter oils. t lo temperatures, eavy oils are extremely viscous and egin to act like solids. n effective sear modulus appears alloing propagation of sear aves. Te ulk moduli of eavy oils increase since te ulk and sear moduli are related it eac oter. Terefore, it is oserved tat te compressional ave velocity of eavy oils increases faster tan tat of ligt oils at lo temperatures. noise if using te S-ave transducers instead of P- ave transducers. It is true even for very eavy oil at lo temperatures at ic te oil acts like a solid. It means it may not e feasile for direct measurement of sear velocity in eavy oil. Water S Transducer Buffer oil Heavy oil On te oter and, e kno from te viscoelastic model of Maxell tat, te sear impedance of te viscous fluid is related it its viscosity ic is one of te important control factors for ave attenuation. (a) Fig. 1 () * 2 2 (1 ) G 1/ 2 1/ 2 V S ere, * is te magnitude of te complex sear impedance, is viscosity, is density, is angle frequency and G is effective ig frequency sear modulus and is called relaxation time. It is apparent tat te sear ave velocity information is useful along it P-ave velocity and density to predict te viscosity and attenuation of te eavy oil. Tis paper descries te metod, facilities and accuracy of S-ave velocity measurements for eavy oil. Measuring Metod Te traditional metod for measurement of fluid velocity is to measure te travel time of te transmission ave and ten te velocity can e derived it te knon travel distance. specially designed camer it to transducers located at its ot ends, te distance of ic ere calirated at different temperature and pressure y distilled ater, as een used and is good for P-ave measurement for a lot of fluid samples. But te transmission signal is very ard to e recognized due to presence of strong metod is designed for indirect measurement of sear ave velocity. Figure 1 sos te principle of tis metod. Te sear ave transducer is coupled it a uffer made of some kind of plastic and te transducer is in te mode of self emission and receiving. Te measurement is done in to steps. Te first is to record te sear reflection from te ottom of te uffer tat is in contact it te reference sample (ater) and te second is to record te same reflection ave en te uffer is in contact it te real sample (eavy oil). Te difference eteen and is caused y te different reflection coefficient (eteen uffer and ater) and (eteen uffer and eavy oil) if all measurement conditions can e controlled. It is ovious tat te ratio of amplitude of reflected ave equals to te ratio of te reflection coefficient. ere, (f) and (f) are te magnitude spectra (function of frequency) of aves and. Te reflection coefficient is -1 ecause te sear impedance of ater is SEG/Houston 2005 nnual Meeting

2 P 2.4 ere, and are sear impedances of ater and uffer respectively. If te impedance of te uffer and te density of te eavy oil are knon, ten te sear velocity V s of te eavy oil can e derived. Pipe P transducer Pressure vessel Vs 1 1 Water or eavy oil Buffer S transducer ere, is te impedance of te eavy oil sample. Water is used as te reference material, since ater can e easily kept in te desired temperature and pressure. In practice, te peak-to-peak amplitude can e taken to replace te spectrum analysis in order to simplify te calculation. Experimental Facility Te scematic of experimental facility is son in figure 2. pressure vessel is used for measuring ater or eavy oil sample. Te sample pressure can e controlled y a pump troug te pipe. Te pressure vessel is put into a at tank filled it anti-freezing liquid ic can e adjusted to te desired temperature. Te S-ave transducer is used to record te reflection aves from te uffer. Meanile, a P- ave transducer is used to receive te transmitted P- ave tat is converted from sear ave. Terefore, te compressional velocity can e measured simultaneously during sear velocity measurement. It is necessary to calirate te facility first. Distilled ater is filled in te vessel and te sear reflection ave and compressional transmission ave are recorded separately at different pressures and temperatures. Te reflection aves are used to uild te relationsip of reflection aves it pressure and temperature. Te travel time of te compressional ave is used to calirate te propagation distance of compressional ave in te sample at various pressure and temperature since te velocity of te distilled ater and te uffer along it te tickness of te uffer are knon. Te calirated distance is needed to calculate P-ave velocity en te eavy oil is filled in te vessel. Virtual Instrument (softare) ritten in LBVIEW for signal spectral analysis as developed in our la. Te S-ave velocity it frequency can e derived y tis softare. Examples Fig. 2 Several eavy oil samples it different PI gravity ave een measured in te la. Figure 3 sos te P- ave and S-ave velocities versus temperature at pressure of aout 0.7 MPa for a sample of eavy oil of PI gravity of 8.0. It is clear ot velocities decreases it increasing temperature ut its a near linear relation for S-ave and non-linear for P-ave. Oter samples tat e measured support same rule. Fig. 3 Te spectral analysis can give canges in te sear velocity it frequency, ile te peak-to-peak value of te signal can only so te sear velocity near te main frequency of te transducer. Figure 4 sos te sear velocity versus frequency tat indicates sear velocity increases sligtly it increasing frequency SEG/Houston 2005 nnual Meeting

3 P 2.4 for a sample it PI of For comparison, te velocity calculated from peak-to-peak is pointed in te figure. We can find from figure 4 tat te sear velocity increases gradually it increasing frequency at various temperatures and pressures. Meanile, te value of sear velocity aout 900 khz from spectral analysis is consistent it tat using peak-to-peak ratio. Te dominant frequency of te reflection ave from spectral analysis for various temperatures is in te range of 800 to 1000 khz it 1MHz transducer. Fig. 4 Te modulus of te sample can e derived if its P-and S-ave velocities and density are knon. Figure 5 sos te ulk and sear moduli of te first sample (PI=8.0) at pressure of aout 0.7MPa. It can e seen tat te cange of te ulk modulus it temperature is very similar to tat of sear modulus. It means ot moduli are dependent on eac oter. We tink it is possile to estimate S-ave velocity from P-ave velocity. Suc relationsip ill e developed in future. Fig. 5 Velocities calculated itout spectral analysis n interesting sample is ax-ric oil it PI of It seems solid at room temperature since it contains muc ax even toug te oil is not very eavy. Te sear velocity is rater lo compared it te real eavy oil. We ave also measured sear velocity of glycerol tat is relatively lo altoug te glycerol as very ig density of gm/cc. Factors Tat ffect ccuracy of Measurement It is ell knon tat te accuracy of P-ave velocity depends on te accuracies of te travel time and calirated lengt of te camer. Te relative error of P ave velocity is less tan ±1% as per our experience. Te sear velocity measurement of eavy oil is an indirect measurement and many more factors may affect te accuracy of te measurement. Te first prolem e met in measurements is te staility of te coupling eteen te transducer and te uffer. Te state of coupling may strongly affect te amplitude of te reflection aves. Because e ave to measure ater and oil sample separately at different temperatures and pressures, e ope te canges due coupling may reverse ile canging te temperature and pressure. We tried several different materials for coupling and finally a fast epoxy is selected and te coupling staility is acceptale. Te main error is proaly from te amplitudes for aves and caused y te stailities of electronic equipment and poer. verage difference of te amplitudes of te reflection aves from statistics is in te range of 3% as per our experience. Terefore, te maximum error may occur en as 3% error and as +3% error or vice versa. s an example, e use te first sample to ceck te maximum error in te measurement of S-ave velocity. Figure 6 sos te maximum asolute and relative errors in S-ave velocity, respectively. It can e seen tat eiter one of te asolute or relative error is iger if te S-ave velocity is loer. In tis figure, te relative error of S-ave velocity may e iger tan 10% (asolute error may e larger ten 36m/s) if te S-ave velocity is less tan 0.36km/s. It means te velocity value is reliale en S-ave velocity is iger tan 0.4km/s, ut may not e accurate enoug if te S-ave velocity is relatively small. On te oter and, te material of uffer may affect te sensitivity of te measurement. Te reflection coefficient ill e close to 1 if te impedance of te uffer is very different from ater or oil; no matter te uffer is in contact eiter ater or oil. Te small difference eteen and may reduce te accuracy of te velocity calculation. Te est uffer sould ave similar impedance as ater and oil. Te density and S-ave velocity of te uffer e used are 1.277g/cm 3 and 1.05km/s respectively. In order to reduce te ave attenuation of te uffer, e select a special plastic material as te uffer it small tickness. Measurement of P-ave Velocity Te P-ave velocity can e otained simultaneously y our facility (figure 2). Te accuracy of te P-ave velocity may e questionale since te transmitted conpressional ave is te convert ave. Our measured data so te differences from te velocities y traditional metod are less tan 2% for various temperatures and pressures. Te accuracy may raise if SEG/Houston 2005 nnual Meeting

4 P 2.4 (a) te correction of travel time for delay of converted ave is taken. It is interesting tat te amplitude of te conpressional ave canges oviously it te temperature (see figure 7 as an example). Te ave amplitude can indicate te attenuation of te P-ave troug te eavy oil and te maximum attenuation in tis sample is at te temperature aout 20C since te amplitude reaces te minimum at te aove temperature. It is knon tat te maximum attenuation may occur en te ave frequency is at te resonance frequency of te material. Te resonance frequency of te eavy oil depends mainly on its viscosity tat may vary it temperature. Fig. 6 Te measurements of S-ave velocity of eavy oil demonstrate te indirect metod for measuring S-ave velocity. In order to otain certain accuracy, special care sould e taken to keep all conditions stale during measurements. Te measurement error may e rater large if te value of S-ave velocity is small and terefore, measurements of S-ave velocity of eavy oil it ig PI gravity at ig temperature may not ave ig accuracy. eferences () Batzle, M and Wang, W-P, 2004, Heavy Oils Glass Point: Fluids/DHI annual report. Beggs, H.D., and oinson, J.., 1975, Estimating te viscosity of crude oil systems: J. Petr. Tec., 27, Blitz, J., 1967, Foundamentals of ultrasonics, Plenum Press, Ne York. Conclusion Fig. 7 Greenood, M. S. and Bamerger, J.., 2002, Measurement of viscosity and sear velocity of a liquid or slurry for on0line process control, Ultrasonics, 39, Liu, J. and Han, D-H., 2004, Measurements of P-ave and S-ave velocity for eavy oil samples: Fluids/DHI annual report. SEG/Houston 2005 nnual Meeting

5 EDITED EFEENCES Note: Tis reference list is a copy-edited version of te reference list sumitted y te autor. eference lists for te 2005 SEG Tecnical Program Expanded stracts ave een copy edited so tat references provided it te online metadata for eac paper ill acieve a ig degree of linking to cited sources tat appear on te We. Measurement of sear ave velocity of eavy oil EFEENCES Batzle, M. and W.-P. Wang, 2004, Heavy Oils Glass Point: Fluids/DHI annual report. Beggs, H. D., and J.. oinson, 1975, Estimating te viscosity of crude oil systems: Journal of Petroleum Tecnology, 27, Blitz, J., 1967, Foundamentals of ultrasonics: Plenum Press. Greenood, M. S., and J.. Bamerger, 2002, Measurement of viscosity and sear velocity of a liquid or slurry for on-line process control: Ultrasonics, 39, Liu, J., and D-H. Han, 2004, Measurements of P-ave and S-ave velocity for eavy oil samples: Fluids/DHI annual report.

LIMITATIONS OF EULER S METHOD FOR NUMERICAL INTEGRATION

LIMITATIONS OF EULER S METHOD FOR NUMERICAL INTEGRATION LIMITATIONS OF EULER S METHOD FOR NUMERICAL INTEGRATION LAURA EVANS.. Introduction Not all differential equations can be explicitly solved for y. Tis can be problematic if we need to know te value of y

More information

Journal of Chemical and Pharmaceutical Research, 2013, 5(12): Research Article

Journal of Chemical and Pharmaceutical Research, 2013, 5(12): Research Article Available online.jocpr.com Journal of emical and Parmaceutical Researc, 013, 5(1):55-531 Researc Article ISSN : 0975-7384 ODEN(USA) : JPR5 Performance and empirical models of a eat pump ater eater system

More information

The Open Petroleum Engineering Journal

The Open Petroleum Engineering Journal Send Orders for Reprints to reprints@bentamscience.ae Te Open Petroleum Engineering Journal, 16, 9, 169-177 169 Te Open Petroleum Engineering Journal Content list available at:.bentamopen.com/topej/ DOI:

More information

CFD calculation of convective heat transfer coefficients and validation Part I: Laminar flow. Annex 41 Kyoto, April 3 rd to 5 th, 2006

CFD calculation of convective heat transfer coefficients and validation Part I: Laminar flow. Annex 41 Kyoto, April 3 rd to 5 th, 2006 CFD calculation of convective eat transfer coefficients and validation Part I: Laminar flow Annex 41 Kyoto, April 3 rd to 5 t, 2006 Adam Neale 1, Dominique Derome 1, Bert Blocken 2 and Jan Carmeliet 2,3

More information

Tribology in Industry

Tribology in Industry Vol. 4, No. (9) 34-4, DOI:.4874/ti.9.4..5 Triology in Industry www.triology.rs RESEARCH A Metod for Predicting Contact Strengt and Life of Arcimedes and Involute Worm Gears, Considering te Effect of Wear

More information

Taylor Series and the Mean Value Theorem of Derivatives

Taylor Series and the Mean Value Theorem of Derivatives 1 - Taylor Series and te Mean Value Teorem o Derivatives Te numerical solution o engineering and scientiic problems described by matematical models oten requires solving dierential equations. Dierential

More information

CALCULATION OF COLLAPSE PRESSURE IN SHALE GAS FORMATION AND THE INFLUENCE OF FORMATION ANISOTROPY

CALCULATION OF COLLAPSE PRESSURE IN SHALE GAS FORMATION AND THE INFLUENCE OF FORMATION ANISOTROPY CALCULATION OF COLLAPSE PRESSURE IN SHALE GAS FORMATION AND THE INFLUENCE OF FORMATION ANISOTROPY L.Hu, J.Deng, F.Deng, H.Lin, C.Yan, Y.Li, H.Liu, W.Cao (Cina University of Petroleum) Sale gas formations

More information

The Basics of Vacuum Technology

The Basics of Vacuum Technology Te Basics of Vacuum Tecnology Grolik Benno, Kopp Joacim January 2, 2003 Basics Many scientific and industrial processes are so sensitive tat is is necessary to omit te disturbing influence of air. For

More information

Precalculus Test 2 Practice Questions Page 1. Note: You can expect other types of questions on the test than the ones presented here!

Precalculus Test 2 Practice Questions Page 1. Note: You can expect other types of questions on the test than the ones presented here! Precalculus Test 2 Practice Questions Page Note: You can expect oter types of questions on te test tan te ones presented ere! Questions Example. Find te vertex of te quadratic f(x) = 4x 2 x. Example 2.

More information

Laboratory experiments and numerical simulation on Bitumen Saturated Carbonates: A Rock Physics Study for 4D Seismology

Laboratory experiments and numerical simulation on Bitumen Saturated Carbonates: A Rock Physics Study for 4D Seismology Laboratory experiments and numerical simulation on Bitumen Saturated Carbonates: A Rock Physics Study for 4D Seismology Arif Rabbani Jason Nycz* Zizhen Wong Doug Schmitt Ken Gray Department of Physics

More information

2.3 Product and Quotient Rules

2.3 Product and Quotient Rules .3. PRODUCT AND QUOTIENT RULES 75.3 Product and Quotient Rules.3.1 Product rule Suppose tat f and g are two di erentiable functions. Ten ( g (x)) 0 = f 0 (x) g (x) + g 0 (x) See.3.5 on page 77 for a proof.

More information

EOQ and EPQ-Partial Backordering-Approximations

EOQ and EPQ-Partial Backordering-Approximations USING A ONSTANT RATE TO APPROXIMATE A LINEARLY HANGING RATE FOR THE EOQ AND EPQ WITH PARTIAL BAKORDERING David W. Pentico, Palumo-Donaue Scool of Business, Duquesne University, Pittsurg, PA 158-18, pentico@duq.edu,

More information

Homework 1. Problem 1 Browse the 331 website to answer: When you should use data symbols on a graph. (Hint check out lab reports...

Homework 1. Problem 1 Browse the 331 website to answer: When you should use data symbols on a graph. (Hint check out lab reports... Homework 1 Problem 1 Browse te 331 website to answer: Wen you sould use data symbols on a grap. (Hint ceck out lab reports...) Solution 1 Use data symbols to sow data points unless tere is so muc data

More information

Forces, centre of gravity, reactions and stability

Forces, centre of gravity, reactions and stability Forces, centre of gravity, reactions and stability Topic areas Mecanical engineering: Centre of gravity Forces Moments Reactions Resolving forces on an inclined plane. Matematics: Angles Trigonometric

More information

Explicit Interleavers for a Repeat Accumulate Accumulate (RAA) code construction

Explicit Interleavers for a Repeat Accumulate Accumulate (RAA) code construction Eplicit Interleavers for a Repeat Accumulate Accumulate RAA code construction Venkatesan Gurusami Computer Science and Engineering University of Wasington Seattle, WA 98195, USA Email: venkat@csasingtonedu

More information

The derivative function

The derivative function Roberto s Notes on Differential Calculus Capter : Definition of derivative Section Te derivative function Wat you need to know already: f is at a point on its grap and ow to compute it. Wat te derivative

More information

Combining functions: algebraic methods

Combining functions: algebraic methods Combining functions: algebraic metods Functions can be added, subtracted, multiplied, divided, and raised to a power, just like numbers or algebra expressions. If f(x) = x 2 and g(x) = x + 2, clearly f(x)

More information

DEBONDING FAILURES OF RC BEAMS STRENGTHENED WITH EXTERNALLY BONDED FRP REINFORCEMENT: BEHAVIOUR AND MODELLING

DEBONDING FAILURES OF RC BEAMS STRENGTHENED WITH EXTERNALLY BONDED FRP REINFORCEMENT: BEHAVIOUR AND MODELLING Asia-Pacific Conference on FRP in Structures (APFIS 2007) S.T. Smit (ed) 2007 International Institute for FRP in Construction DEBONDING FAILURES OF RC BEAMS STRENGTHENED WITH EXTERNALLY BONDED FRP REINFORCEMENT:

More information

SIMG Solution Set #5

SIMG Solution Set #5 SIMG-303-0033 Solution Set #5. Describe completely te state of polarization of eac of te following waves: (a) E [z,t] =ˆxE 0 cos [k 0 z ω 0 t] ŷe 0 cos [k 0 z ω 0 t] Bot components are traveling down te

More information

Microstrip Antennas- Rectangular Patch

Microstrip Antennas- Rectangular Patch April 4, 7 rect_patc_tl.doc Page of 6 Microstrip Antennas- Rectangular Patc (Capter 4 in Antenna Teory, Analysis and Design (nd Edition) by Balanis) Sown in Figures 4. - 4.3 Easy to analyze using transmission

More information

UNIVERSITY OF SASKATCHEWAN Department of Physics and Engineering Physics

UNIVERSITY OF SASKATCHEWAN Department of Physics and Engineering Physics UNIVERSITY O SASKATCHEWAN Department of Pysics and Engineering Pysics Pysics 117.3 MIDTERM EXAM Regular Sitting NAME: (Last) Please Print (Given) Time: 90 minutes STUDENT NO.: LECTURE SECTION (please ceck):

More information

Theoretical Analysis of Flow Characteristics and Bearing Load for Mass-produced External Gear Pump

Theoretical Analysis of Flow Characteristics and Bearing Load for Mass-produced External Gear Pump TECHNICAL PAPE Teoretical Analysis of Flow Caracteristics and Bearing Load for Mass-produced External Gear Pump N. YOSHIDA Tis paper presents teoretical equations for calculating pump flow rate and bearing

More information

Consider a function f we ll specify which assumptions we need to make about it in a minute. Let us reformulate the integral. 1 f(x) dx.

Consider a function f we ll specify which assumptions we need to make about it in a minute. Let us reformulate the integral. 1 f(x) dx. Capter 2 Integrals as sums and derivatives as differences We now switc to te simplest metods for integrating or differentiating a function from its function samples. A careful study of Taylor expansions

More information

How to Find the Derivative of a Function: Calculus 1

How to Find the Derivative of a Function: Calculus 1 Introduction How to Find te Derivative of a Function: Calculus 1 Calculus is not an easy matematics course Te fact tat you ave enrolled in suc a difficult subject indicates tat you are interested in te

More information

WINKLER PLATES BY THE BOUNDARY KNOT METHOD

WINKLER PLATES BY THE BOUNDARY KNOT METHOD WINKLER PLATES BY THE BOUNARY KNOT ETHO Sofía Roble, sroble@fing.edu.uy Berardi Sensale, sensale@fing.edu.uy Facultad de Ingeniería, Julio Herrera y Reissig 565, ontevideo Abstract. Tis paper describes

More information

Lecture XVII. Abstract We introduce the concept of directional derivative of a scalar function and discuss its relation with the gradient operator.

Lecture XVII. Abstract We introduce the concept of directional derivative of a scalar function and discuss its relation with the gradient operator. Lecture XVII Abstract We introduce te concept of directional derivative of a scalar function and discuss its relation wit te gradient operator. Directional derivative and gradient Te directional derivative

More information

Sedimentation. Treatment Processes. Sedimentation. Sedimentation. Sedimentation. Sedimentation. CIVL 1112 Water Treatment - Sedimentation 1/7

Sedimentation. Treatment Processes. Sedimentation. Sedimentation. Sedimentation. Sedimentation. CIVL 1112 Water Treatment - Sedimentation 1/7 CIVL 111 Water Treatment - 1/7 Treatment Processes is te donards movement of an object relative to its surrounding medium, due to te force of gravity. Screening Aeration Preclorination locculation Coagulation

More information

The total error in numerical differentiation

The total error in numerical differentiation AMS 147 Computational Metods and Applications Lecture 08 Copyrigt by Hongyun Wang, UCSC Recap: Loss of accuracy due to numerical cancellation A B 3, 3 ~10 16 In calculating te difference between A and

More information

Derivation Of The Schwarzschild Radius Without General Relativity

Derivation Of The Schwarzschild Radius Without General Relativity Derivation Of Te Scwarzscild Radius Witout General Relativity In tis paper I present an alternative metod of deriving te Scwarzscild radius of a black ole. Te metod uses tree of te Planck units formulas:

More information

Investigating Euler s Method and Differential Equations to Approximate π. Lindsay Crowl August 2, 2001

Investigating Euler s Method and Differential Equations to Approximate π. Lindsay Crowl August 2, 2001 Investigating Euler s Metod and Differential Equations to Approximate π Lindsa Crowl August 2, 2001 Tis researc paper focuses on finding a more efficient and accurate wa to approximate π. Suppose tat x

More information

Practice Problem Solutions: Exam 1

Practice Problem Solutions: Exam 1 Practice Problem Solutions: Exam 1 1. (a) Algebraic Solution: Te largest term in te numerator is 3x 2, wile te largest term in te denominator is 5x 2 3x 2 + 5. Tus lim x 5x 2 2x 3x 2 x 5x 2 = 3 5 Numerical

More information

The effects of shear stress on the lubrication performances of oil film of large-scale mill bearing

The effects of shear stress on the lubrication performances of oil film of large-scale mill bearing Universit of Wollongong Researc Online Facult of Engineering - Papers (Arcive) Facult of Engineering and Information Sciences 9 Te effects of sear stress on te lubrication performances of oil film of large-scale

More information

Notes on wavefunctions II: momentum wavefunctions

Notes on wavefunctions II: momentum wavefunctions Notes on wavefunctions II: momentum wavefunctions and uncertainty Te state of a particle at any time is described by a wavefunction ψ(x). Tese wavefunction must cange wit time, since we know tat particles

More information

1watt=1W=1kg m 2 /s 3

1watt=1W=1kg m 2 /s 3 Appendix A Matematics Appendix A.1 Units To measure a pysical quantity, you need a standard. Eac pysical quantity as certain units. A unit is just a standard we use to compare, e.g. a ruler. In tis laboratory

More information

Quaternion Dynamics, Part 1 Functions, Derivatives, and Integrals. Gary D. Simpson. rev 01 Aug 08, 2016.

Quaternion Dynamics, Part 1 Functions, Derivatives, and Integrals. Gary D. Simpson. rev 01 Aug 08, 2016. Quaternion Dynamics, Part 1 Functions, Derivatives, and Integrals Gary D. Simpson gsim1887@aol.com rev 1 Aug 8, 216 Summary Definitions are presented for "quaternion functions" of a quaternion. Polynomial

More information

Dynamics and Relativity

Dynamics and Relativity Dynamics and Relativity Stepen Siklos Lent term 2011 Hand-outs and examples seets, wic I will give out in lectures, are available from my web site www.damtp.cam.ac.uk/user/stcs/dynamics.tml Lecture notes,

More information

Empirical models for estimating liquefaction-induced lateral spread displacement

Empirical models for estimating liquefaction-induced lateral spread displacement Empirical models for estimating liquefaction-induced lateral spread displacement J.J. Zang and J.X. Zao Institute of Geological & Nuclear Sciences Ltd, Lower Hutt, New Zealand. 2004 NZSEE Conference ABSTRACT:

More information

Exam 1 Review Solutions

Exam 1 Review Solutions Exam Review Solutions Please also review te old quizzes, and be sure tat you understand te omework problems. General notes: () Always give an algebraic reason for your answer (graps are not sufficient),

More information

Average Rate of Change

Average Rate of Change Te Derivative Tis can be tougt of as an attempt to draw a parallel (pysically and metaporically) between a line and a curve, applying te concept of slope to someting tat isn't actually straigt. Te slope

More information

Problem Solving. Problem Solving Process

Problem Solving. Problem Solving Process Problem Solving One of te primary tasks for engineers is often solving problems. It is wat tey are, or sould be, good at. Solving engineering problems requires more tan just learning new terms, ideas and

More information

THE IDEA OF DIFFERENTIABILITY FOR FUNCTIONS OF SEVERAL VARIABLES Math 225

THE IDEA OF DIFFERENTIABILITY FOR FUNCTIONS OF SEVERAL VARIABLES Math 225 THE IDEA OF DIFFERENTIABILITY FOR FUNCTIONS OF SEVERAL VARIABLES Mat 225 As we ave seen, te definition of derivative for a Mat 111 function g : R R and for acurveγ : R E n are te same, except for interpretation:

More information

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

Nondestructive Monitoring of Setting and Hardening of Portland Cement Mortar with Sonic Methods Nondestructive Monitoring of Setting and Hardening of Portland Cement Mortar ith Sonic Methods Thomas Voigt, Northestern University, Evanston, USA Surendra P. Shah, Northestern University, Evanston, USA

More information

Mathematics 5 Worksheet 11 Geometry, Tangency, and the Derivative

Mathematics 5 Worksheet 11 Geometry, Tangency, and the Derivative Matematics 5 Workseet 11 Geometry, Tangency, and te Derivative Problem 1. Find te equation of a line wit slope m tat intersects te point (3, 9). Solution. Te equation for a line passing troug a point (x

More information

Analysis of Static and Dynamic Load on Hydrostatic Bearing with Variable Viscosity and Pressure

Analysis of Static and Dynamic Load on Hydrostatic Bearing with Variable Viscosity and Pressure Indian Journal of Science and Tecnology Supplementary Article Analysis of Static and Dynamic Load on Hydrostatic Bearing wit Variable Viscosity and Pressure V. Srinivasan* Professor, Scool of Mecanical

More information

3 Minority carrier profiles (the hyperbolic functions) Consider a

3 Minority carrier profiles (the hyperbolic functions) Consider a Microelectronic Devices and Circuits October 9, 013 - Homework #3 Due Nov 9, 013 1 Te pn junction Consider an abrupt Si pn + junction tat as 10 15 acceptors cm -3 on te p-side and 10 19 donors on te n-side.

More information

Solve exponential equations in one variable using a variety of strategies. LEARN ABOUT the Math. What is the half-life of radon?

Solve exponential equations in one variable using a variety of strategies. LEARN ABOUT the Math. What is the half-life of radon? 8.5 Solving Exponential Equations GOAL Solve exponential equations in one variable using a variety of strategies. LEARN ABOUT te Mat All radioactive substances decrease in mass over time. Jamie works in

More information

NCCI: Simple methods for second order effects in portal frames

NCCI: Simple methods for second order effects in portal frames NCC: Simple metods for second order effects in portal frames NCC: Simple metods for second order effects in portal frames NCC: Simple metods for second order effects in portal frames Tis NCC presents information

More information

Notes: Most of the material in this chapter is taken from Young and Freedman, Chap. 12.

Notes: Most of the material in this chapter is taken from Young and Freedman, Chap. 12. Capter 6. Fluid Mecanics Notes: Most of te material in tis capter is taken from Young and Freedman, Cap. 12. 6.1 Fluid Statics Fluids, i.e., substances tat can flow, are te subjects of tis capter. But

More information

Physically Based Modeling: Principles and Practice Implicit Methods for Differential Equations

Physically Based Modeling: Principles and Practice Implicit Methods for Differential Equations Pysically Based Modeling: Principles and Practice Implicit Metods for Differential Equations David Baraff Robotics Institute Carnegie Mellon University Please note: Tis document is 997 by David Baraff

More information

A STUDY ON THE GROUND MOTION CHARACTERISTICS OF TAIPEI BASIN, TAIWAN, BASED ON OBSERVED STRONG MOTIONS AND MEASURED MICROTREMORS

A STUDY ON THE GROUND MOTION CHARACTERISTICS OF TAIPEI BASIN, TAIWAN, BASED ON OBSERVED STRONG MOTIONS AND MEASURED MICROTREMORS A STUDY ON THE GROUND MOTION CHARACTERISTICS OF TAIPEI BASIN, TAIWAN, BASED ON OBSERVED STRONG MOTIONS AND MEASURED MICROTREMORS Ying Liu 1, Kentaro Motoki 2 and Kazuo Seo 2 1 Eartquake Engineer Group,

More information

Differential Calculus (The basics) Prepared by Mr. C. Hull

Differential Calculus (The basics) Prepared by Mr. C. Hull Differential Calculus Te basics) A : Limits In tis work on limits, we will deal only wit functions i.e. tose relationsips in wic an input variable ) defines a unique output variable y). Wen we work wit

More information

Summary. Simple model for kerogen maturity (Carcione, 2000)

Summary. Simple model for kerogen maturity (Carcione, 2000) Malleswar Yenugu* and De-hua Han, University of Houston, USA Summary The conversion of kerogen to oil/gas will build up overpressure. Overpressure is caused by conversion of solid kerogen to fluid hydrocarbons

More information

FINITE DIFFERENCE APPROXIMATION TO WAVE EQUATION TO SIMULATE SLOSHING OF GROUND SUPPORTED TANKS FOR EARTHQUAKE LOADINGS

FINITE DIFFERENCE APPROXIMATION TO WAVE EQUATION TO SIMULATE SLOSHING OF GROUND SUPPORTED TANKS FOR EARTHQUAKE LOADINGS e 7 t International Conference on Sustainale uilt Environment, Earl s egency Hotel, Kandy, Sri anka from 16 t to 18 t Decemer 16 ICSE16-81 FINIE DIFFEENCE APPOXIAION O WAVE EQUAION O SIUAE SOSHING OF GOUND

More information

LIMITS AND DERIVATIVES CONDITIONS FOR THE EXISTENCE OF A LIMIT

LIMITS AND DERIVATIVES CONDITIONS FOR THE EXISTENCE OF A LIMIT LIMITS AND DERIVATIVES Te limit of a function is defined as te value of y tat te curve approaces, as x approaces a particular value. Te limit of f (x) as x approaces a is written as f (x) approaces, as

More information

= 0 and states ''hence there is a stationary point'' All aspects of the proof dx must be correct (c)

= 0 and states ''hence there is a stationary point'' All aspects of the proof dx must be correct (c) Paper 1: Pure Matematics 1 Mark Sceme 1(a) (i) (ii) d d y 3 1x 4x x M1 A1 d y dx 1.1b 1.1b 36x 48x A1ft 1.1b Substitutes x = into teir dx (3) 3 1 4 Sows d y 0 and states ''ence tere is a stationary point''

More information

Cubic Functions: Local Analysis

Cubic Functions: Local Analysis Cubic function cubing coefficient Capter 13 Cubic Functions: Local Analysis Input-Output Pairs, 378 Normalized Input-Output Rule, 380 Local I-O Rule Near, 382 Local Grap Near, 384 Types of Local Graps

More information

A New Model for the Prediction of the Dog-bone Shape in Steel Mills

A New Model for the Prediction of the Dog-bone Shape in Steel Mills ISIJ International, Vol. (), No. 6, pp. 9 7 New Model for te Prediction of te Dog-one Sape in Steel Mills Duckjoong YUN, ) Dongun L, ) Jaeoo KIM ) and Sangmoo HNG ) ) Department of Mecanical ngineering,

More information

A Reconsideration of Matter Waves

A Reconsideration of Matter Waves A Reconsideration of Matter Waves by Roger Ellman Abstract Matter waves were discovered in te early 20t century from teir wavelengt, predicted by DeBroglie, Planck's constant divided by te particle's momentum,

More information

The Hangingstone steam-assisted gravity drainage

The Hangingstone steam-assisted gravity drainage SPECIAL Heavy SECTION: oil H e a v y o i l Elastic property changes in a bitumen reservoir during steam injection AYATO KATO, University of Houston, USA SHIGENOBU ONOZUKA, JOGMEC, Chiba, Japan TORU NAKAYAMA,

More information

Simulation and verification of a plate heat exchanger with a built-in tap water accumulator

Simulation and verification of a plate heat exchanger with a built-in tap water accumulator Simulation and verification of a plate eat excanger wit a built-in tap water accumulator Anders Eriksson Abstract In order to test and verify a compact brazed eat excanger (CBE wit a built-in accumulation

More information

Announcements. Exam 4 - Review of important concepts

Announcements. Exam 4 - Review of important concepts Announcements 1. Exam 4 starts Friday! a. Available in esting Center from Friday, Dec 7 (opening time), up to Monday, Dec 10 at 4:00 pm. i. Late fee if you start your exam after 4 pm b. Covers C. 9-1 (up

More information

A general articulation angle stability model for non-slewing articulated mobile cranes on slopes *

A general articulation angle stability model for non-slewing articulated mobile cranes on slopes * tecnical note 3 general articulation angle stability model for non-slewing articulated mobile cranes on slopes * J Wu, L uzzomi and M Hodkiewicz Scool of Mecanical and Cemical Engineering, University of

More information

SECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY

SECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY (Section 3.2: Derivative Functions and Differentiability) 3.2.1 SECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY LEARNING OBJECTIVES Know, understand, and apply te Limit Definition of te Derivative

More information

Velocity Measurements of Pore Fluids at Pressure and Temperature: Application to bitumen

Velocity Measurements of Pore Fluids at Pressure and Temperature: Application to bitumen Velocity Measurements of Pore Fluids at Pressure and Temperature: Application to bitumen Arif Rabbani 1*, Douglas R Schmitt 1, Jason Nycz 2, and Ken Gray 3 1 Institute for Geophysical Research, Department

More information

De-Coupler Design for an Interacting Tanks System

De-Coupler Design for an Interacting Tanks System IOSR Journal of Electrical and Electronics Engineering (IOSR-JEEE) e-issn: 2278-1676,p-ISSN: 2320-3331, Volume 7, Issue 3 (Sep. - Oct. 2013), PP 77-81 De-Coupler Design for an Interacting Tanks System

More information

Evaluation and Accurate Estimation from Petrophysical Parameters of a Reservoir

Evaluation and Accurate Estimation from Petrophysical Parameters of a Reservoir American Journal of Environmental Engineering and Science 2016; 3(2): 68-74 ttp://www.aascit.org/journal/ajees ISSN: 2381-1153 (Print); ISSN: 2381-1161 (Online) Evaluation and Accurate Estimation from

More information

Measurement of cement s particle size distribution by the buoyancy weighing-bar method

Measurement of cement s particle size distribution by the buoyancy weighing-bar method Measurement of cement s particle size distribution by te buoyancy weiging-bar metod *Rondang Tambun, Nofriko ratama, Ely, and Farida Hanum Department of Cemical Engineering, Faculty of Engineering, University

More information

Math 34A Practice Final Solutions Fall 2007

Math 34A Practice Final Solutions Fall 2007 Mat 34A Practice Final Solutions Fall 007 Problem Find te derivatives of te following functions:. f(x) = 3x + e 3x. f(x) = x + x 3. f(x) = (x + a) 4. Is te function 3t 4t t 3 increasing or decreasing wen

More information

232 Calculus and Structures

232 Calculus and Structures 3 Calculus and Structures CHAPTER 17 JUSTIFICATION OF THE AREA AND SLOPE METHODS FOR EVALUATING BEAMS Calculus and Structures 33 Copyrigt Capter 17 JUSTIFICATION OF THE AREA AND SLOPE METHODS 17.1 THE

More information

Analysis: The speed of the proton is much less than light speed, so we can use the

Analysis: The speed of the proton is much less than light speed, so we can use the Section 1.3: Wave Proerties of Classical Particles Tutorial 1 Practice, age 634 1. Given: 1.8! 10 "5 kg # m/s; 6.63! 10 "34 J #s Analysis: Use te de Broglie relation, λ. Solution:! 6.63 " 10#34 kg $ m

More information

Logic Effort Revisited

Logic Effort Revisited Logic Effort Revisited Mark This note ill take another look at logical effort, first revieing the asic idea ehind logical effort, and then looking at some of the more sutle issues in sizing transistors..0

More information

The entransy dissipation minimization principle under given heat duty and heat transfer area conditions

The entransy dissipation minimization principle under given heat duty and heat transfer area conditions Article Engineering Termopysics July 2011 Vol.56 No.19: 2071 2076 doi: 10.1007/s11434-010-4189-x SPECIAL TOPICS: Te entransy dissipation minimization principle under given eat duty and eat transfer area

More information

HOW TO DEAL WITH FFT SAMPLING INFLUENCES ON ADEV CALCULATIONS

HOW TO DEAL WITH FFT SAMPLING INFLUENCES ON ADEV CALCULATIONS HOW TO DEAL WITH FFT SAMPLING INFLUENCES ON ADEV CALCULATIONS Po-Ceng Cang National Standard Time & Frequency Lab., TL, Taiwan 1, Lane 551, Min-Tsu Road, Sec. 5, Yang-Mei, Taoyuan, Taiwan 36 Tel: 886 3

More information

Chapters 19 & 20 Heat and the First Law of Thermodynamics

Chapters 19 & 20 Heat and the First Law of Thermodynamics Capters 19 & 20 Heat and te First Law of Termodynamics Te Zerot Law of Termodynamics Te First Law of Termodynamics Termal Processes Te Second Law of Termodynamics Heat Engines and te Carnot Cycle Refrigerators,

More information

Pumping Heat with Quantum Ratchets

Pumping Heat with Quantum Ratchets Pumping Heat wit Quantum Ratcets T. E. Humprey a H. Linke ab R. Newbury a a Scool of Pysics University of New Sout Wales UNSW Sydney 5 Australia b Pysics Department University of Oregon Eugene OR 9743-74

More information

Numerical Analysis MTH603. dy dt = = (0) , y n+1. We obtain yn. Therefore. and. Copyright Virtual University of Pakistan 1

Numerical Analysis MTH603. dy dt = = (0) , y n+1. We obtain yn. Therefore. and. Copyright Virtual University of Pakistan 1 Numerical Analysis MTH60 PREDICTOR CORRECTOR METHOD Te metods presented so far are called single-step metods, were we ave seen tat te computation of y at t n+ tat is y n+ requires te knowledge of y n only.

More information

Order of Accuracy. ũ h u Ch p, (1)

Order of Accuracy. ũ h u Ch p, (1) Order of Accuracy 1 Terminology We consider a numerical approximation of an exact value u. Te approximation depends on a small parameter, wic can be for instance te grid size or time step in a numerical

More information

Math 124. Section 2.6: Limits at infinity & Horizontal Asymptotes. 1 x. lim

Math 124. Section 2.6: Limits at infinity & Horizontal Asymptotes. 1 x. lim Mat 4 Section.6: Limits at infinity & Horizontal Asymptotes Tolstoy, Count Lev Nikolgevic (88-90) A man is like a fraction wose numerator is wat e is and wose denominator is wat e tinks of imself. Te larger

More information

KEY CONCEPT: THE DERIVATIVE

KEY CONCEPT: THE DERIVATIVE Capter Two KEY CONCEPT: THE DERIVATIVE We begin tis capter by investigating te problem of speed: How can we measure te speed of a moving object at a given instant in time? Or, more fundamentally, wat do

More information

DEFORMABLE GROUND INFLUENCE ON THE FRICTION DRIVE TRANSMISSION BETWEEN DRIVE PULLEY AND ELASTOMETRIC TRACK

DEFORMABLE GROUND INFLUENCE ON THE FRICTION DRIVE TRANSMISSION BETWEEN DRIVE PULLEY AND ELASTOMETRIC TRACK Journal of KONE Powertrain and Transport, Vol. 8, No. 0 EFORMLE GROUN INFLUENCE ON THE FRICTION RIVE TRNMIION ETWEEN RIVE PULLEY N ELTOMETRIC TRCK Roert Czaanowski Wrocaw University of Tecnology, epartment

More information

Numerical Differentiation

Numerical Differentiation Numerical Differentiation Finite Difference Formulas for te first derivative (Using Taylor Expansion tecnique) (section 8.3.) Suppose tat f() = g() is a function of te variable, and tat as 0 te function

More information

Vidmantas Jokūbaitis a, Linas Juknevičius b, *, Remigijus Šalna c

Vidmantas Jokūbaitis a, Linas Juknevičius b, *, Remigijus Šalna c Availale online at www.sciencedirect.com Procedia Engineering 57 ( 203 ) 466 472 t International Conference on Modern Building Materials, Structures and Tecniques, MBMST 203 Conditions for Failure of Normal

More information

Preface. Here are a couple of warnings to my students who may be here to get a copy of what happened on a day that you missed.

Preface. Here are a couple of warnings to my students who may be here to get a copy of what happened on a day that you missed. Preface Here are my online notes for my course tat I teac ere at Lamar University. Despite te fact tat tese are my class notes, tey sould be accessible to anyone wanting to learn or needing a refreser

More information

Exercises for numerical differentiation. Øyvind Ryan

Exercises for numerical differentiation. Øyvind Ryan Exercises for numerical differentiation Øyvind Ryan February 25, 2013 1. Mark eac of te following statements as true or false. a. Wen we use te approximation f (a) (f (a +) f (a))/ on a computer, we can

More information

Time (hours) Morphine sulfate (mg)

Time (hours) Morphine sulfate (mg) Mat Xa Fall 2002 Review Notes Limits and Definition of Derivative Important Information: 1 According to te most recent information from te Registrar, te Xa final exam will be eld from 9:15 am to 12:15

More information

Introduction. Learning Objectives. On completion of this chapter you will be able to:

Introduction. Learning Objectives. On completion of this chapter you will be able to: Introduction Learning Objectives On completion of tis capter you will be able to: 1. Define Compton Effect. 2. Derive te sift in incident ligt wavelengt and Compton wavelengt. 3. Explain ow te Compton

More information

Comment on Experimental observations of saltwater up-coning

Comment on Experimental observations of saltwater up-coning 1 Comment on Experimental observations of saltwater up-coning H. Zang 1,, D.A. Barry 2 and G.C. Hocking 3 1 Griffit Scool of Engineering, Griffit University, Gold Coast Campus, QLD 4222, Australia. Tel.:

More information

Optimal Shape Design of a Two-dimensional Asymmetric Diffsuer in Turbulent Flow

Optimal Shape Design of a Two-dimensional Asymmetric Diffsuer in Turbulent Flow THE 5 TH ASIAN COMPUTAITIONAL FLUID DYNAMICS BUSAN, KOREA, OCTOBER 7 ~ OCTOBER 30, 003 Optimal Sape Design of a Two-dimensional Asymmetric Diffsuer in Turbulent Flow Seokyun Lim and Haeceon Coi. Center

More information

Firooz Arash (a,b) , which have caused confusion in the normalization of ZEUS data is discussed and resolved. PCAC I. INTRODUCTION

Firooz Arash (a,b) , which have caused confusion in the normalization of ZEUS data is discussed and resolved. PCAC I. INTRODUCTION PION STRUCTURE FUNCTION F π 2 In te Valon Model Firooz Aras (a,b) (a) Center for Teoretical Pysics and Matematics, AEOI P.O.Box 365-8486, Teran, Iran (b) Pysics Department, AmirKabir University (Tafres

More information

A method of Lagrange Galerkin of second order in time. Une méthode de Lagrange Galerkin d ordre deux en temps

A method of Lagrange Galerkin of second order in time. Une méthode de Lagrange Galerkin d ordre deux en temps A metod of Lagrange Galerkin of second order in time Une métode de Lagrange Galerkin d ordre deux en temps Jocelyn Étienne a a DAMTP, University of Cambridge, Wilberforce Road, Cambridge CB3 0WA, Great-Britain.

More information

Excursions in Computing Science: Week v Milli-micro-nano-..math Part II

Excursions in Computing Science: Week v Milli-micro-nano-..math Part II Excursions in Computing Science: Week v Milli-micro-nano-..mat Part II T. H. Merrett McGill University, Montreal, Canada June, 5 I. Prefatory Notes. Cube root of 8. Almost every calculator as a square-root

More information

Click here to see an animation of the derivative

Click here to see an animation of the derivative Differentiation Massoud Malek Derivative Te concept of derivative is at te core of Calculus; It is a very powerful tool for understanding te beavior of matematical functions. It allows us to optimize functions,

More information

. If lim. x 2 x 1. f(x+h) f(x)

. If lim. x 2 x 1. f(x+h) f(x) Review of Differential Calculus Wen te value of one variable y is uniquely determined by te value of anoter variable x, ten te relationsip between x and y is described by a function f tat assigns a value

More information

Applied Linear Statistical Models. Simultaneous Inference Topics. Simultaneous Estimation of β 0 and β 1 Issues. Simultaneous Inference. Dr.

Applied Linear Statistical Models. Simultaneous Inference Topics. Simultaneous Estimation of β 0 and β 1 Issues. Simultaneous Inference. Dr. Applied Linear Statistical Models Simultaneous Inference Dr. DH Jones Simultaneous Inference Topics Simultaneous estimation of β 0 and β 1 Bonferroni Metod Simultaneous estimation of several mean responses

More information

The Complexity of Computing the MCD-Estimator

The Complexity of Computing the MCD-Estimator Te Complexity of Computing te MCD-Estimator Torsten Bernolt Lerstul Informatik 2 Universität Dortmund, Germany torstenbernolt@uni-dortmundde Paul Fiscer IMM, Danisc Tecnical University Kongens Lyngby,

More information

The structure of the atoms

The structure of the atoms Te structure of te atoms Atomos = indivisible University of Pécs, Medical Scool, Dept. Biopysics All tat exists are atoms and empty space; everyting else is merely tougt to exist. Democritus, 415 B.C.

More information

Lecture 15. Interpolation II. 2 Piecewise polynomial interpolation Hermite splines

Lecture 15. Interpolation II. 2 Piecewise polynomial interpolation Hermite splines Lecture 5 Interpolation II Introduction In te previous lecture we focused primarily on polynomial interpolation of a set of n points. A difficulty we observed is tat wen n is large, our polynomial as to

More information

Quantization of electrical conductance

Quantization of electrical conductance 1 Introduction Quantization of electrical conductance Te resistance of a wire in te classical Drude model of metal conduction is given by RR = ρρρρ AA, were ρρ, AA and ll are te conductivity of te material,

More information

Development of new and validation of existing convection correlations for rooms with displacement ventilation systems

Development of new and validation of existing convection correlations for rooms with displacement ventilation systems Energy and Buildings 38 (2006) 163 173 www.elsevier.com/locate/enbuild Development of new and validation of existing convection correlations for rooms wit displacement ventilation systems Atila Novoselac

More information

Fast Computation of Capacitance Matrix and Potential Distribution for Multiconductor in Non-Homogenous Multilayered Dielectric Media

Fast Computation of Capacitance Matrix and Potential Distribution for Multiconductor in Non-Homogenous Multilayered Dielectric Media Excerpt from te Proceedings of te OMSOL onference 2009 Boston Fast omputation of apacitance Matrix and Potential Distribution for Multiconductor in Non-Homogenous Multilayered Dielectric Media S. M. Musa

More information