SPE Improved Permeability Prediction Relations for Low Permeability Sands

Size: px
Start display at page:

Download "SPE Improved Permeability Prediction Relations for Low Permeability Sands"

Transcription

1 SPE Imroved Permebility Prediction Reltions for Low Permebility Snds Frncois-Andre Florence, Texs A&M University T.A. Blsingme, Texs A&M University Dertment of Petroleum Engineering Texs A&M University College Sttion, TX Slide /2

2 Concet: Gs Slige Gs Slige Phenomenon: "Gs slige" occurs when the men free th () of the gs molecules is "not negligible" comred to effective ore throt rdius: the gs molecules "sli" on the surfces of the orous medi (i.e., the velocity v of the gs molecules t the wll is non-zero) : Gs Flow v y v x ( y 0) v 0, with v 0 dv c dy y0 Gs Slige /2 v x From: Kundt, A. nd Wrburg, E.: "Über Reibung und Wärmeleitung verdünnter Gse, " Poggendorfs Annlen der Physi und Chemie (875), 55, 337. Ignoring the gs slige effect led to n overestimtion of the gs flowrte nd the core ermebility (Drcy's lw). Slide 2/2

3 Concet: Gs Slige Gs Slige Phenomenon: The "gs slige" effect for orous medi ws first documented by Klinenberg. The gs ermebility tends to limiting vlue t n infinite men ressure which is referred to s equivlent liquid ermebility or Klinenberg-corrected ermebility. b K m Gs Slige 2/2 Where is the Klinenbergcorrected ermebility (identified s the bsolute ermebility), nd b K is the gs slige fctor. From: Klinenberg, L.J.: "The Permebility of Porous Medi to Liquids nd Gses," API Drilling nd Production Prctice (94). Slide 3/2

4 Concet: Correltions Correltions between b K nd the etrohysicl dt: The objective of this wor is to ccurtely correlte the b K -term with the Klinenberg-corrected ermebility nd other etrohysicl dt with the form: b K f (,...) This form is then substituted in the Klinenberg eqution: For given ir of vlues of gs ermebility nd men ressure The Klinenberg-corrected ermebility is then comuted by solving: Two tye of correltions re reorted in the literture: b K vs. b K vs. / f (,...) m f (,...) m 0 (Heid et l, Jones nd Owens) (Smth nd Keighin, squre-root (this wor)) Slide 4/2

5 Correltions: Heid et l. b K vs. /2 Discussion: Heid et l. Correltion is "ccetble" for vlues of ermebility greter thn md. Dt sctter must not be neglected Slide 5/2

6 Correltions: Jones-Owens nd Heid et l b K vs. 2/2 Discussion: Jones-Owens Although the correltion ers ccurte, the correltion errors re most significnt for low ermebility dt. Slide 6/2

7 Correltions: b K vs. / Bsed on 0 tight snd core smles, the Smth nd Keighin correltion reltes b K to /: b K b K vs. /3 From: Smth, K., nd Keighin, C.W.: Fctors Affecting Gs Slige in Tight Sndstones, er SPE 9872 resented t the 98 SPE/DOE Low Permebility Symosium, Denver, CO., My (98) Discussion: Smth-Keighin The correltion is ccetble. Reltively smll dtbse. Slide 7/2

8 Correltions: b K vs. / The lter formultion is interesting s it is close to the theoreticl definition of cillry rdius r s function of the squre-root of ( /): r Klinenberg defines b K by: b K m 4 r m It is ossible to derive rigorously correltion between b K nd the squre-root of ( /) using the revious definitions nd the clssicl definition of the men free th of the gs molecules: This "squre-root" correltion hs the following form: / RT (, T ) / 2 [ (, T )] M 0.5 Where is rmeter (intercet) deending on the tye of gs used for the core b K Where m is the men free th of the gs molecules t the men ressure m flow exeriment usully nitrogen, for which = b K vs. 2/3 Slide 8/2

9 bk vs. 3/3 Correltions: Smth-Keighin, squre-root The squre-root model ers to give better results. The Smth-Keighin Model mtches minly their dt. Slide 9/2 Discussion: Smth-Keighin

10 Correltions: Comrison on Exmle Cse Errors re generlly less thn 30 ercent. Discussion: Lower Cotton Vlley Formtion No. The Jones-Owens nd Smth-Keighin overestimte the ermebility; the squre-root model underestimtes. Slide 0/2

11 Correltions: Comrison on Exmle Cse Errors re generlly less thn 25 ercent, lthough some outliers exhibit more thn 50 ercent error. Discussion: Lower Cotton Vlley Formtion No. 2 The Jones-Owens/Smth-Keighin correltions for > 0.00 md. The squre-root model gives good results for low ermebilities. Slide /2

12 Microflow model: Concet: The gs slige henomenon occurs s subset of lrger re of study nown s the theory of rrefied gses this theory exerienced substntil growth in the lst 50 yers (eronutics, erosce, MEMS ) The Klinenberg eqution is bsed on the results of this theory vilble in 94: b 4 K c m, with c m r The theory of rrefied gses describes the flow of gs t very low ressures, or in smll micro-chnnels. The flow regime for gs flowing in micro-chnnel is tyiclly determined by the vlue of the Knudsen number (Kn) : Kn l chr Where is the men free th of the gs molecules (i.e., the verge distnce (length) between 2 consecutive moleculr interctions), nd l chr is the chrcteristic length of the flow geometry (e.g., chnnel height, ie rdius). Slide 2/2

13 Microflow model: Concets The clssicl definition of the men free th from thermodynmics is: RT (, T ) / 2, where (, T ) M Limits of the different flow regimes, s function of the length chrcteristic of the geometry, l chr, nd the recirocl men free th normlized t tmosheric conditions nd 300 K. The lines defining the vrious Knudsen number regimes re bsed on ir t isotherml conditions (Modified from Krnidis nd Beso, 2002). Slide 3/2

14 Microflow model: Unified Flow Model: In 2002, Krnidis nd Beso develoed "unified flow model" for gs flows in micro-ies, vlid over the entire rnge of flow regimes: (for micro-ie of rdius r nd length L) r 4 4Kn 28 - q ( Kn) Kn, where( Kn) tn 4Kn 2 8 Kn L From: Krnidis, G.E. nd Beso, A.: Micro-flows, Fundmentls nd Simultion, Sringer-Verlg, New-Yor (2002). By nlogy with Klinenberg's derivtion using Poiseuille eqution (corrected for slige), we derive microflow model: Kn tn 4Kn Kn 2 5 Kn Where [+(Kn) Kn] is defined s the rrefction coefficient, ccounting for the chnge of gs density. Slide 4/2

15 Microflow model: Correltion of the Knudsen Number (Kn): The Knudsen number cnnot be mesured by direct lbortory mesurement we need to define "seudo"-knudsen number which is function of the mesurble rmeters ( m,, ) Assuming the isotherml flow of n idel gs, Kn is inversely roortionl to m nd l chr (tyiclly, the ore throt rdius, or function of nd ) For smle for which nd re nown, the "seudo"- Knudsen number is comuted by solving the following eqution: tn - 4Kn 0.4 4Kn Kn 0 Kn Slide 5/2

16 Microflow model: Aliction using Field Dt: To demonstrte the liction of the "seudo"-knudsen number (Kn ), we only use the Lower Cotton Vlley dt (No. 2): the seudo-knudsen numbers (Kn ) cn be relted to the vilble etrohysicl dt (,, nd ), yielding two reltions: Reltion : Kn 0.4 m Reltion 2: Kn 2.62 These two reltions cn be lied to the microflow model, yielding to n exlicit eqution (with Reltion ) or n imlicit eqution (with Reltion 2) tht cn be solved for : m tn - 4Kn 0.4 4Kn Kn 0 Kn Slide 6/2

17 Microflow model: Field Dt Results (2/2) Discussion: Imlicit nd Exlicit models The two models chieve firly good correltion the imlicit model seems to be more ccurte. Slide 7/2

18 "Fully imlicit" formultion : The Klinenberg-corrected ermebility my not be vilble in rctice. The rocess of solving simultneously for Kn nd requires n imlicit formultion with the form: The coefficients 0, nd 2 hve to be determined simultneously with the estimtion of. This cn be hndled (with sufficient quntity of dt) by using non-liner regression methods on the following eqution: The coefficients ( i ) comuted re not "universl" nd must be comuted for ech dtsets. Microflow model: Field Dt Results 2 0 m Kn / 4 4 tn m m m Slide 8/2

19 Microflow model: Field Dt Results (2/2) Discussion: "Fully imlicit" formultion The fully imlicit formultion chieves good mtch for this dt set. Slide 9/2

20 Summry, Issues, nd Chllenges: Summry: Primry objective: Jones-Owens, Smth-Keighin nd squre-root correltions cn be used Smth-Keighin nd squre-root correltions would be referred s they offer more consistency with the theory. Secondry objective: new theoreticl model (microflow model) for the estimtion of ermebility from gs flow mesurement hs been develoed nd needs to be further vlidted. Issues nd Chllenges: Continued vlidtion of the new microflow model requires lrge quntity of relible dt. Slide 20/2

21 SPE Imroved Permebility Prediction Reltions for Low Permebility Snds End of Presenttion Frncois-Andre Florence, Texs A&M University T.A. Blsingme, Texs A&M University Dertment of Petroleum Engineering Texs A&M University College Sttion, TX Slide 2/2

IPTC November 2005 IPTC (Doha, Qatar)

IPTC November 2005 IPTC (Doha, Qatar) November 005 IPTC (Doh, Qtr A Moifie Purcell/Burine Moel for Estimting Absolute Permebility from Mercury-Injection Cillry Pressure Dt CC Huet, Texs A&M University, JA Rushing, Anrko Petroleum Cor, KE Newshm,

More information

Part I: Basic Concepts of Thermodynamics

Part I: Basic Concepts of Thermodynamics Prt I: Bsic Concepts o Thermodynmics Lecture 4: Kinetic Theory o Gses Kinetic Theory or rel gses 4-1 Kinetic Theory or rel gses Recll tht or rel gses: (i The volume occupied by the molecules under ordinry

More information

A proof of the strong twin prime conjecture

A proof of the strong twin prime conjecture A roof of the strong twin rime conjecture Men-Jw Ho # (retired), Chou-Jung Hsu, *, Wi-Jne Ho b Dertment of Industril Mngement # Nn Ki University of Technology Nn-Tou 54, Tiwn b Dertment of Medicinl Botnicl

More information

Module 2: Rate Law & Stoichiomtery (Chapter 3, Fogler)

Module 2: Rate Law & Stoichiomtery (Chapter 3, Fogler) CHE 309: Chemicl Rection Engineering Lecture-8 Module 2: Rte Lw & Stoichiomtery (Chpter 3, Fogler) Topics to be covered in tody s lecture Thermodynmics nd Kinetics Rection rtes for reversible rections

More information

CHM Physical Chemistry I Chapter 1 - Supplementary Material

CHM Physical Chemistry I Chapter 1 - Supplementary Material CHM 3410 - Physicl Chemistry I Chpter 1 - Supplementry Mteril For review of some bsic concepts in mth, see Atkins "Mthemticl Bckground 1 (pp 59-6), nd "Mthemticl Bckground " (pp 109-111). 1. Derivtion

More information

Rel Gses 1. Gses (N, CO ) which don t obey gs lws or gs eqution P=RT t ll pressure nd tempertures re clled rel gses.. Rel gses obey gs lws t extremely low pressure nd high temperture. Rel gses devited

More information

Properties of Lorenz Curves for Transformed Income Distributions

Properties of Lorenz Curves for Transformed Income Distributions Theoreticl Economics etters 22 2 487-493 htt://ddoiorg/4236/tel22259 Published Online December 22 (htt://wwwscirporg/journl/tel) Proerties of orenz Curves for Trnsformed Income Distributions John Fellmn

More information

THERMAL EXPANSION COEFFICIENT OF WATER FOR VOLUMETRIC CALIBRATION

THERMAL EXPANSION COEFFICIENT OF WATER FOR VOLUMETRIC CALIBRATION XX IMEKO World Congress Metrology for Green Growth September 9,, Busn, Republic of Kore THERMAL EXPANSION COEFFICIENT OF WATER FOR OLUMETRIC CALIBRATION Nieves Medin Hed of Mss Division, CEM, Spin, mnmedin@mityc.es

More information

The Thermodynamics of Aqueous Electrolyte Solutions

The Thermodynamics of Aqueous Electrolyte Solutions 18 The Thermodynmics of Aqueous Electrolyte Solutions As discussed in Chpter 10, when slt is dissolved in wter or in other pproprite solvent, the molecules dissocite into ions. In queous solutions, strong

More information

1B40 Practical Skills

1B40 Practical Skills B40 Prcticl Skills Comining uncertinties from severl quntities error propgtion We usully encounter situtions where the result of n experiment is given in terms of two (or more) quntities. We then need

More information

Classical Mechanics. From Molecular to Con/nuum Physics I WS 11/12 Emiliano Ippoli/ October, 2011

Classical Mechanics. From Molecular to Con/nuum Physics I WS 11/12 Emiliano Ippoli/ October, 2011 Clssicl Mechnics From Moleculr to Con/nuum Physics I WS 11/12 Emilino Ippoli/ October, 2011 Wednesdy, October 12, 2011 Review Mthemtics... Physics Bsic thermodynmics Temperture, idel gs, kinetic gs theory,

More information

Ideal Gas behaviour: summary

Ideal Gas behaviour: summary Lecture 4 Rel Gses Idel Gs ehviour: sury We recll the conditions under which the idel gs eqution of stte Pn is vlid: olue of individul gs olecules is neglected No interctions (either ttrctive or repulsive)

More information

The Properties of Stars

The Properties of Stars 10/11/010 The Properties of Strs sses Using Newton s Lw of Grvity to Determine the ss of Celestil ody ny two prticles in the universe ttrct ech other with force tht is directly proportionl to the product

More information

Thermal Diffusivity. Paul Hughes. Department of Physics and Astronomy The University of Manchester Manchester M13 9PL. Second Year Laboratory Report

Thermal Diffusivity. Paul Hughes. Department of Physics and Astronomy The University of Manchester Manchester M13 9PL. Second Year Laboratory Report Therml iffusivity Pul Hughes eprtment of Physics nd Astronomy The University of nchester nchester 3 9PL Second Yer Lbortory Report Nov 4 Abstrct We investigted the therml diffusivity of cylindricl block

More information

Black oils Correlations Comparative Study

Black oils Correlations Comparative Study Reservoir Technologies Blck oils Correltions Comprtive Study Dr. Muhmmd Al-Mrhoun, Mnging Director Sturdy, 26 April, 2014 Copyright 2008, NExT, All rights reserved Blck oils Correltions Introduction to

More information

Math 131. Numerical Integration Larson Section 4.6

Math 131. Numerical Integration Larson Section 4.6 Mth. Numericl Integrtion Lrson Section. This section looks t couple of methods for pproimting definite integrls numericlly. The gol is to get good pproimtion of the definite integrl in problems where n

More information

Joule-Thomson effect TEP

Joule-Thomson effect TEP Joule-homson effect EP elted oics el gs; intrinsic energy; Gy-Lussc theory; throttling; n der Wls eqution; n der Wls force; inverse Joule- homson effect; inversion temerture. Princile A strem of gs is

More information

Tests for the Ratio of Two Poisson Rates

Tests for the Ratio of Two Poisson Rates Chpter 437 Tests for the Rtio of Two Poisson Rtes Introduction The Poisson probbility lw gives the probbility distribution of the number of events occurring in specified intervl of time or spce. The Poisson

More information

P 3 (x) = f(0) + f (0)x + f (0) 2. x 2 + f (0) . In the problem set, you are asked to show, in general, the n th order term is a n = f (n) (0)

P 3 (x) = f(0) + f (0)x + f (0) 2. x 2 + f (0) . In the problem set, you are asked to show, in general, the n th order term is a n = f (n) (0) 1 Tylor polynomils In Section 3.5, we discussed how to pproximte function f(x) round point in terms of its first derivtive f (x) evluted t, tht is using the liner pproximtion f() + f ()(x ). We clled this

More information

Physics 201 Lab 3: Measurement of Earth s local gravitational field I Data Acquisition and Preliminary Analysis Dr. Timothy C. Black Summer I, 2018

Physics 201 Lab 3: Measurement of Earth s local gravitational field I Data Acquisition and Preliminary Analysis Dr. Timothy C. Black Summer I, 2018 Physics 201 Lb 3: Mesurement of Erth s locl grvittionl field I Dt Acquisition nd Preliminry Anlysis Dr. Timothy C. Blck Summer I, 2018 Theoreticl Discussion Grvity is one of the four known fundmentl forces.

More information

Satellite Retrieval Data Assimilation

Satellite Retrieval Data Assimilation tellite etrievl Dt Assimiltion odgers C. D. Inverse Methods for Atmospheric ounding: Theor nd Prctice World cientific Pu. Co. Hckensck N.J. 2000 Chpter 3 nd Chpter 8 Dve uhl Artist depiction of NAA terr

More information

Intro to Nuclear and Particle Physics (5110)

Intro to Nuclear and Particle Physics (5110) Intro to Nucler nd Prticle Physics (5110) Feb, 009 The Nucler Mss Spectrum The Liquid Drop Model //009 1 E(MeV) n n(n-1)/ E/[ n(n-1)/] (MeV/pir) 1 C 16 O 0 Ne 4 Mg 7.7 14.44 19.17 8.48 4 5 6 6 10 15.4.41

More information

63. Representation of functions as power series Consider a power series. ( 1) n x 2n for all 1 < x < 1

63. Representation of functions as power series Consider a power series. ( 1) n x 2n for all 1 < x < 1 3 9. SEQUENCES AND SERIES 63. Representtion of functions s power series Consider power series x 2 + x 4 x 6 + x 8 + = ( ) n x 2n It is geometric series with q = x 2 nd therefore it converges for ll q =

More information

LECTURE 14. Dr. Teresa D. Golden University of North Texas Department of Chemistry

LECTURE 14. Dr. Teresa D. Golden University of North Texas Department of Chemistry LECTURE 14 Dr. Teres D. Golden University of North Texs Deprtment of Chemistry Quntittive Methods A. Quntittive Phse Anlysis Qulittive D phses by comprison with stndrd ptterns. Estimte of proportions of

More information

Properties of Integrals, Indefinite Integrals. Goals: Definition of the Definite Integral Integral Calculations using Antiderivatives

Properties of Integrals, Indefinite Integrals. Goals: Definition of the Definite Integral Integral Calculations using Antiderivatives Block #6: Properties of Integrls, Indefinite Integrls Gols: Definition of the Definite Integrl Integrl Clcultions using Antiderivtives Properties of Integrls The Indefinite Integrl 1 Riemnn Sums - 1 Riemnn

More information

Estimation of the particle concentration in hydraulic liquid by the in-line automatic particle counter based on the CMOS image sensor

Estimation of the particle concentration in hydraulic liquid by the in-line automatic particle counter based on the CMOS image sensor Glyndŵr University Reserch Online Conference Presenttion Estimtion of the prticle concentrtion in hydrulic liquid by the in-line utomtic prticle counter bsed on the CMOS imge sensor Kornilin, D.V., Kudryvtsev,

More information

Maximum Likelihood Estimation for Allele Frequencies. Biostatistics 666

Maximum Likelihood Estimation for Allele Frequencies. Biostatistics 666 Mximum Likelihood Estimtion for Allele Frequencies Biosttistics 666 Previous Series of Lectures: Introduction to Colescent Models Comuttionlly efficient frmework Alterntive to forwrd simultions Amenble

More information

Math Calculus with Analytic Geometry II

Math Calculus with Analytic Geometry II orem of definite Mth 5.0 with Anlytic Geometry II Jnury 4, 0 orem of definite If < b then b f (x) dx = ( under f bove x-xis) ( bove f under x-xis) Exmple 8 0 3 9 x dx = π 3 4 = 9π 4 orem of definite Problem

More information

Lecture 13 - Linking E, ϕ, and ρ

Lecture 13 - Linking E, ϕ, and ρ Lecture 13 - Linking E, ϕ, nd ρ A Puzzle... Inner-Surfce Chrge Density A positive point chrge q is locted off-center inside neutrl conducting sphericl shell. We know from Guss s lw tht the totl chrge on

More information

A027 Uncertainties in Local Anisotropy Estimation from Multi-offset VSP Data

A027 Uncertainties in Local Anisotropy Estimation from Multi-offset VSP Data A07 Uncertinties in Locl Anisotropy Estimtion from Multi-offset VSP Dt M. Asghrzdeh* (Curtin University), A. Bon (Curtin University), R. Pevzner (Curtin University), M. Urosevic (Curtin University) & B.

More information

Synoptic Meteorology I: Finite Differences September Partial Derivatives (or, Why Do We Care About Finite Differences?

Synoptic Meteorology I: Finite Differences September Partial Derivatives (or, Why Do We Care About Finite Differences? Synoptic Meteorology I: Finite Differences 16-18 September 2014 Prtil Derivtives (or, Why Do We Cre About Finite Differences?) With the exception of the idel gs lw, the equtions tht govern the evolution

More information

Probabilistic Investigation of Sensitivities of Advanced Test- Analysis Model Correlation Methods

Probabilistic Investigation of Sensitivities of Advanced Test- Analysis Model Correlation Methods Probbilistic Investigtion of Sensitivities of Advnced Test- Anlysis Model Correltion Methods Liz Bergmn, Mtthew S. Allen, nd Dniel C. Kmmer Dept. of Engineering Physics University of Wisconsin-Mdison Rndll

More information

2 b. , a. area is S= 2π xds. Again, understand where these formulas came from (pages ).

2 b. , a. area is S= 2π xds. Again, understand where these formulas came from (pages ). AP Clculus BC Review Chpter 8 Prt nd Chpter 9 Things to Know nd Be Ale to Do Know everything from the first prt of Chpter 8 Given n integrnd figure out how to ntidifferentite it using ny of the following

More information

The Velocity Factor of an Insulated Two-Wire Transmission Line

The Velocity Factor of an Insulated Two-Wire Transmission Line The Velocity Fctor of n Insulted Two-Wire Trnsmission Line Problem Kirk T. McDonld Joseph Henry Lbortories, Princeton University, Princeton, NJ 08544 Mrch 7, 008 Estimte the velocity fctor F = v/c nd the

More information

Families of Solutions to Bernoulli ODEs

Families of Solutions to Bernoulli ODEs In the fmily of solutions to the differentil eqution y ry dx + = it is shown tht vrition of the initil condition y( 0 = cuses horizontl shift in the solution curve y = f ( x, rther thn the verticl shift

More information

Deriving hydraulic conductivity function from soil column tests

Deriving hydraulic conductivity function from soil column tests Deriving hydrulic conductivity function from soil column tests Yvonne Lins, Mri Dtchev & Tom Schnz DFG Fo 444 - TP 4 Experimentelle und theoretische Untersuchungen teilgesättigter Reibungsmterilien Content

More information

Electron Correlation Methods

Electron Correlation Methods Electron Correltion Methods HF method: electron-electron interction is replced by n verge interction E HF c E 0 E HF E 0 exct ground stte energy E HF HF energy for given bsis set HF Ec 0 - represents mesure

More information

The Regulated and Riemann Integrals

The Regulated and Riemann Integrals Chpter 1 The Regulted nd Riemnn Integrls 1.1 Introduction We will consider severl different pproches to defining the definite integrl f(x) dx of function f(x). These definitions will ll ssign the sme vlue

More information

+ x 2 dω 2 = c 2 dt 2 +a(t) [ 2 dr 2 + S 1 κx 2 /R0

+ x 2 dω 2 = c 2 dt 2 +a(t) [ 2 dr 2 + S 1 κx 2 /R0 Notes for Cosmology course, fll 2005 Cosmic Dynmics Prelude [ ds 2 = c 2 dt 2 +(t) 2 dx 2 ] + x 2 dω 2 = c 2 dt 2 +(t) [ 2 dr 2 + S 1 κx 2 /R0 2 κ (r) 2 dω 2] nd x = S κ (r) = r, R 0 sin(r/r 0 ), R 0 sinh(r/r

More information

Influence of Morphology and Packing Properties on the Permeability of Fine Particle Agglomerates

Influence of Morphology and Packing Properties on the Permeability of Fine Particle Agglomerates Influence of Morhology nd Pcking Proerties on the Permebility of Fine Prticle Agglomertes Alberto Scurti Ic Mns-Zloczower * Det. Mcromoleculr Science nd Engineering Donld L. Feke Det. Mcromoleculr Science

More information

13: Diffusion in 2 Energy Groups

13: Diffusion in 2 Energy Groups 3: Diffusion in Energy Groups B. Rouben McMster University Course EP 4D3/6D3 Nucler Rector Anlysis (Rector Physics) 5 Sept.-Dec. 5 September Contents We study the diffusion eqution in two energy groups

More information

Math 1B, lecture 4: Error bounds for numerical methods

Math 1B, lecture 4: Error bounds for numerical methods Mth B, lecture 4: Error bounds for numericl methods Nthn Pflueger 4 September 0 Introduction The five numericl methods descried in the previous lecture ll operte by the sme principle: they pproximte the

More information

A Case Study of Fluid Transport in Shale Crushed Samples: Experiment and Interpretation

A Case Study of Fluid Transport in Shale Crushed Samples: Experiment and Interpretation SCA2014-071 1/12 A Cse Study of Fluid Trnsport in Shle Crushed Smples: Experiment nd Interprettion Tinho Wu, Zheng Jing, Dongxio Zhng Energy nd Resources Engineering, College of Engineering Peking University,

More information

Definition of Continuity: The function f(x) is continuous at x = a if f(a) exists and lim

Definition of Continuity: The function f(x) is continuous at x = a if f(a) exists and lim Mth 9 Course Summry/Study Guide Fll, 2005 [1] Limits Definition of Limit: We sy tht L is the limit of f(x) s x pproches if f(x) gets closer nd closer to L s x gets closer nd closer to. We write lim f(x)

More information

W. We shall do so one by one, starting with I 1, and we shall do it greedily, trying

W. We shall do so one by one, starting with I 1, and we shall do it greedily, trying Vitli covers 1 Definition. A Vitli cover of set E R is set V of closed intervls with positive length so tht, for every δ > 0 nd every x E, there is some I V with λ(i ) < δ nd x I. 2 Lemm (Vitli covering)

More information

AQA Further Pure 1. Complex Numbers. Section 1: Introduction to Complex Numbers. The number system

AQA Further Pure 1. Complex Numbers. Section 1: Introduction to Complex Numbers. The number system Complex Numbers Section 1: Introduction to Complex Numbers Notes nd Exmples These notes contin subsections on The number system Adding nd subtrcting complex numbers Multiplying complex numbers Complex

More information

ADVANCEMENT OF THE CLOSELY COUPLED PROBES POTENTIAL DROP TECHNIQUE FOR NDE OF SURFACE CRACKS

ADVANCEMENT OF THE CLOSELY COUPLED PROBES POTENTIAL DROP TECHNIQUE FOR NDE OF SURFACE CRACKS ADVANCEMENT OF THE CLOSELY COUPLED PROBES POTENTIAL DROP TECHNIQUE FOR NDE OF SURFACE CRACKS F. Tkeo 1 nd M. Sk 1 Hchinohe Ntionl College of Technology, Hchinohe, Jpn; Tohoku University, Sendi, Jpn Abstrct:

More information

Section 5.1 #7, 10, 16, 21, 25; Section 5.2 #8, 9, 15, 20, 27, 30; Section 5.3 #4, 6, 9, 13, 16, 28, 31; Section 5.4 #7, 18, 21, 23, 25, 29, 40

Section 5.1 #7, 10, 16, 21, 25; Section 5.2 #8, 9, 15, 20, 27, 30; Section 5.3 #4, 6, 9, 13, 16, 28, 31; Section 5.4 #7, 18, 21, 23, 25, 29, 40 Mth B Prof. Audrey Terrs HW # Solutions by Alex Eustis Due Tuesdy, Oct. 9 Section 5. #7,, 6,, 5; Section 5. #8, 9, 5,, 7, 3; Section 5.3 #4, 6, 9, 3, 6, 8, 3; Section 5.4 #7, 8,, 3, 5, 9, 4 5..7 Since

More information

Chapter 4 Contravariance, Covariance, and Spacetime Diagrams

Chapter 4 Contravariance, Covariance, and Spacetime Diagrams Chpter 4 Contrvrince, Covrince, nd Spcetime Digrms 4. The Components of Vector in Skewed Coordintes We hve seen in Chpter 3; figure 3.9, tht in order to show inertil motion tht is consistent with the Lorentz

More information

Tom BLASINGAME Texas A&M U. Slide 1

Tom BLASINGAME Texas A&M U. Slide 1 Petroleum Engineering 620 Fluid Flow in Petroleum Reservoirs Petrophysics Lecture 1 Introduction to Porosity and Permeability Concepts Tom BLASINGAME t-blasingame@tamu.edu Texas A&M U. Slide 1 From: Jorden,

More information

( ) as a fraction. Determine location of the highest

( ) as a fraction. Determine location of the highest AB Clculus Exm Review Sheet - Solutions A. Preclculus Type prolems A1 A2 A3 A4 A5 A6 A7 This is wht you think of doing Find the zeros of f ( x). Set function equl to 0. Fctor or use qudrtic eqution if

More information

CHEMICAL KINETICS

CHEMICAL KINETICS CHEMICAL KINETICS Long Answer Questions: 1. Explin the following terms with suitble exmples ) Averge rte of Rection b) Slow nd Fst Rections c) Order of Rection d) Moleculrity of Rection e) Activtion Energy

More information

AB Calculus Review Sheet

AB Calculus Review Sheet AB Clculus Review Sheet Legend: A Preclculus, B Limits, C Differentil Clculus, D Applictions of Differentil Clculus, E Integrl Clculus, F Applictions of Integrl Clculus, G Prticle Motion nd Rtes This is

More information

5.7 Improper Integrals

5.7 Improper Integrals 458 pplictions of definite integrls 5.7 Improper Integrls In Section 5.4, we computed the work required to lift pylod of mss m from the surfce of moon of mss nd rdius R to height H bove the surfce of the

More information

Follow this and additional works at: Part of the Materials Science and Engineering Commons

Follow this and additional works at:  Part of the Materials Science and Engineering Commons Engineering Conferences Interntionl ECI Digitl Archives 5th Interntionl Conference on Porous Medi nd Their Applictions in Science, Engineering nd Industry Refereed Proceedings Summer 6-4-4 Numericl simultion

More information

MODELLING OF MAGNESIUM OXIDE APPLICATION FOR STRUVITE PRECIPITATION

MODELLING OF MAGNESIUM OXIDE APPLICATION FOR STRUVITE PRECIPITATION MDELLING F MAGNESIUM XIDE APPLIATIN FR STRUVITE PREIPITATIN J. Mzierski * nd S. Połwski ** * Polish Acdemy of Science, Institute of Environmentl Engineering, ul. M. Skłodowskiej-urie, -89 Zbrze, Polnd

More information

( ) where f ( x ) is a. AB Calculus Exam Review Sheet. A. Precalculus Type problems. Find the zeros of f ( x).

( ) where f ( x ) is a. AB Calculus Exam Review Sheet. A. Precalculus Type problems. Find the zeros of f ( x). AB Clculus Exm Review Sheet A. Preclculus Type prolems A1 Find the zeros of f ( x). This is wht you think of doing A2 A3 Find the intersection of f ( x) nd g( x). Show tht f ( x) is even. A4 Show tht f

More information

Math 8 Winter 2015 Applications of Integration

Math 8 Winter 2015 Applications of Integration Mth 8 Winter 205 Applictions of Integrtion Here re few importnt pplictions of integrtion. The pplictions you my see on n exm in this course include only the Net Chnge Theorem (which is relly just the Fundmentl

More information

Math 113 Exam 1-Review

Math 113 Exam 1-Review Mth 113 Exm 1-Review September 26, 2016 Exm 1 covers 6.1-7.3 in the textbook. It is dvisble to lso review the mteril from 5.3 nd 5.5 s this will be helpful in solving some of the problems. 6.1 Are Between

More information

Machine Design II Prof. K.Gopinath & Prof. M.M.Mayuram. Drum Brakes. Among the various types of devices to be studied, based on their practical use,

Machine Design II Prof. K.Gopinath & Prof. M.M.Mayuram. Drum Brakes. Among the various types of devices to be studied, based on their practical use, chine Design II Prof. K.Gointh & Prof...yurm Drum Brkes Among the vrious tyes of devices to be studied, bsed on their rcticl use, the discussion will be limited to Drum brkes of the following tyes which

More information

Vadose Zone Hydrology

Vadose Zone Hydrology Objectives Vdose Zone Hydrology 1. Review bsic concepts nd terminology of soil physics. 2. Understnd the role of wter-tble dynmics in GW-SW interction. Drcy s lw is useful in region A. Some knowledge of

More information

Summary Information and Formulae MTH109 College Algebra

Summary Information and Formulae MTH109 College Algebra Generl Formuls Summry Informtion nd Formule MTH109 College Algebr Temperture: F = 9 5 C + 32 nd C = 5 ( 9 F 32 ) F = degrees Fhrenheit C = degrees Celsius Simple Interest: I = Pr t I = Interest erned (chrged)

More information

INTRODUCTION. The three general approaches to the solution of kinetics problems are:

INTRODUCTION. The three general approaches to the solution of kinetics problems are: INTRODUCTION According to Newton s lw, prticle will ccelerte when it is subjected to unblnced forces. Kinetics is the study of the reltions between unblnced forces nd the resulting chnges in motion. The

More information

New data structures to reduce data size and search time

New data structures to reduce data size and search time New dt structures to reduce dt size nd serch time Tsuneo Kuwbr Deprtment of Informtion Sciences, Fculty of Science, Kngw University, Hirtsuk-shi, Jpn FIT2018 1D-1, No2, pp1-4 Copyright (c)2018 by The Institute

More information

Haplotype Frequencies and Linkage Disequilibrium. Biostatistics 666

Haplotype Frequencies and Linkage Disequilibrium. Biostatistics 666 Hlotye Frequencies nd Linkge isequilirium iosttistics 666 Lst Lecture Genotye Frequencies llele Frequencies Phenotyes nd Penetrnces Hrdy-Weinerg Equilirium Simle demonstrtion Exercise: NO2 nd owel isese

More information

New Expansion and Infinite Series

New Expansion and Infinite Series Interntionl Mthemticl Forum, Vol. 9, 204, no. 22, 06-073 HIKARI Ltd, www.m-hikri.com http://dx.doi.org/0.2988/imf.204.4502 New Expnsion nd Infinite Series Diyun Zhng College of Computer Nnjing University

More information

The Atwood Machine OBJECTIVE INTRODUCTION APPARATUS THEORY

The Atwood Machine OBJECTIVE INTRODUCTION APPARATUS THEORY The Atwood Mchine OBJECTIVE To derive the ening of Newton's second lw of otion s it pplies to the Atwood chine. To explin how ss iblnce cn led to the ccelertion of the syste. To deterine the ccelertion

More information

Measuring Electron Work Function in Metal

Measuring Electron Work Function in Metal n experiment of the Electron topic Mesuring Electron Work Function in Metl Instructor: 梁生 Office: 7-318 Emil: shling@bjtu.edu.cn Purposes 1. To understnd the concept of electron work function in metl nd

More information

International ejournals

International ejournals Avilble online t www.interntionlejournls.com ISSN 0976-4 Interntionl ejournls Interntionl ejournl of Mthemtics nd Engineering 93 (0) 85 855 RADIAL VIBRATIONS IN MICROPOLAR THIN SPHERICAL SHELL R.Srinivs*,

More information

Designing Information Devices and Systems I Discussion 8B

Designing Information Devices and Systems I Discussion 8B Lst Updted: 2018-10-17 19:40 1 EECS 16A Fll 2018 Designing Informtion Devices nd Systems I Discussion 8B 1. Why Bother With Thévenin Anywy? () Find Thévenin eqiuvlent for the circuit shown elow. 2kΩ 5V

More information

Introduction to Determinants. Remarks. Remarks. The determinant applies in the case of square matrices

Introduction to Determinants. Remarks. Remarks. The determinant applies in the case of square matrices Introduction to Determinnts Remrks The determinnt pplies in the cse of squre mtrices squre mtrix is nonsingulr if nd only if its determinnt not zero, hence the term determinnt Nonsingulr mtrices re sometimes

More information

SUMMER KNOWHOW STUDY AND LEARNING CENTRE

SUMMER KNOWHOW STUDY AND LEARNING CENTRE SUMMER KNOWHOW STUDY AND LEARNING CENTRE Indices & Logrithms 2 Contents Indices.2 Frctionl Indices.4 Logrithms 6 Exponentil equtions. Simplifying Surds 13 Opertions on Surds..16 Scientific Nottion..18

More information

Applications of Bernoulli s theorem. Lecture - 7

Applications of Bernoulli s theorem. Lecture - 7 Applictions of Bernoulli s theorem Lecture - 7 Prcticl Applictions of Bernoulli s Theorem The Bernoulli eqution cn be pplied to gret mny situtions not just the pipe flow we hve been considering up to now.

More information

Lecture 14: Quadrature

Lecture 14: Quadrature Lecture 14: Qudrture This lecture is concerned with the evlution of integrls fx)dx 1) over finite intervl [, b] The integrnd fx) is ssumed to be rel-vlues nd smooth The pproximtion of n integrl by numericl

More information

The International Association for the Properties of Water and Steam. Release on the Ionization Constant of H 2 O

The International Association for the Properties of Water and Steam. Release on the Ionization Constant of H 2 O IAPWS R-7 The Interntionl Assocition for the Properties of Wter nd Stem Lucerne, Sitzerlnd August 7 Relese on the Ioniztion Constnt of H O 7 The Interntionl Assocition for the Properties of Wter nd Stem

More information

Kepler's Three LAWS. Universal Gravitation Chapter 12. Heliocentric Model. Geocentric Model. Other Models. Johannes Kepler

Kepler's Three LAWS. Universal Gravitation Chapter 12. Heliocentric Model. Geocentric Model. Other Models. Johannes Kepler Universl Grvittion Chpter 1 Johnnes Kepler Johnnes Kepler ws Germn mthemticin, stronomer nd strologer, nd key figure in the 17th century Scientific revolution. He is best known for his lws of plnetry motion,

More information

Consequently, the temperature must be the same at each point in the cross section at x. Let:

Consequently, the temperature must be the same at each point in the cross section at x. Let: HW 2 Comments: L1-3. Derive the het eqution for n inhomogeneous rod where the therml coefficients used in the derivtion of the het eqution for homogeneous rod now become functions of position x in the

More information

Lecture XVII. Vector functions, vector and scalar fields Definition 1 A vector-valued function is a map associating vectors to real numbers, that is

Lecture XVII. Vector functions, vector and scalar fields Definition 1 A vector-valued function is a map associating vectors to real numbers, that is Lecture XVII Abstrct We introduce the concepts of vector functions, sclr nd vector fields nd stress their relevnce in pplied sciences. We study curves in three-dimensionl Eucliden spce nd introduce the

More information

A-Level Mathematics Transition Task (compulsory for all maths students and all further maths student)

A-Level Mathematics Transition Task (compulsory for all maths students and all further maths student) A-Level Mthemtics Trnsition Tsk (compulsory for ll mths students nd ll further mths student) Due: st Lesson of the yer. Length: - hours work (depending on prior knowledge) This trnsition tsk provides revision

More information

Math 100 Review Sheet

Math 100 Review Sheet Mth 100 Review Sheet Joseph H. Silvermn December 2010 This outline of Mth 100 is summry of the mteril covered in the course. It is designed to be study id, but it is only n outline nd should be used s

More information

Math& 152 Section Integration by Parts

Math& 152 Section Integration by Parts Mth& 5 Section 7. - Integrtion by Prts Integrtion by prts is rule tht trnsforms the integrl of the product of two functions into other (idelly simpler) integrls. Recll from Clculus I tht given two differentible

More information

The p-folded Cumulative Distribution Function and the Mean Absolute Deviation from the p-quantile

The p-folded Cumulative Distribution Function and the Mean Absolute Deviation from the p-quantile The -folded Cumultive Distribution Function nd the Men Absolute Devition from the -quntile Jing-Ho Xue,, D. Michel Titterington b Dertment of Sttisticl Science, University College London, London WC1E 6BT,

More information

5.1 How do we Measure Distance Traveled given Velocity? Student Notes

5.1 How do we Measure Distance Traveled given Velocity? Student Notes . How do we Mesure Distnce Trveled given Velocity? Student Notes EX ) The tle contins velocities of moving cr in ft/sec for time t in seconds: time (sec) 3 velocity (ft/sec) 3 A) Lel the x-xis & y-xis

More information

List all of the possible rational roots of each equation. Then find all solutions (both real and imaginary) of the equation. 1.

List all of the possible rational roots of each equation. Then find all solutions (both real and imaginary) of the equation. 1. Mth Anlysis CP WS 4.X- Section 4.-4.4 Review Complete ech question without the use of grphing clcultor.. Compre the mening of the words: roots, zeros nd fctors.. Determine whether - is root of 0. Show

More information

CHAPTER 20: Second Law of Thermodynamics

CHAPTER 20: Second Law of Thermodynamics CHAER 0: Second Lw of hermodynmics Responses to Questions 3. kg of liquid iron will hve greter entropy, since it is less ordered thn solid iron nd its molecules hve more therml motion. In ddition, het

More information

Lecture 6: Diffusion and Reaction kinetics

Lecture 6: Diffusion and Reaction kinetics Lecture 6: Diffusion nd Rection kinetics 1-1-1 Lecture pln: diffusion thermodynmic forces evolution of concentrtion distribution rection rtes nd methods to determine them rection mechnism in terms of the

More information

4 The dynamical FRW universe

4 The dynamical FRW universe 4 The dynmicl FRW universe 4.1 The Einstein equtions Einstein s equtions G µν = T µν (7) relte the expnsion rte (t) to energy distribution in the universe. On the left hnd side is the Einstein tensor which

More information

A little harder example. A block sits at rest on a flat surface. The block is held down by its weight. What is the interaction pair for the weight?

A little harder example. A block sits at rest on a flat surface. The block is held down by its weight. What is the interaction pair for the weight? Neton s Ls of Motion (ges 9-99) 1. An object s velocit vector v remins constnt if nd onl if the net force cting on the object is zero.. hen nonzero net force cts on n object, the object s velocit chnges.

More information

Rates of chemical reactions

Rates of chemical reactions Rtes of chemicl rections Mesuring rtes of chemicl rections Experimentl mesuring progress of the rection Monitoring pressure in the rection involving gses 2 NO( g) 4 NO ( g) + O ( g) 2 5 2 2 n(1 α) 2αn

More information

Math 113 Exam 2 Practice

Math 113 Exam 2 Practice Mth Em Prctice Februry, 8 Em will cover sections 6.5, 7.-7.5 nd 7.8. This sheet hs three sections. The first section will remind you bout techniques nd formuls tht you should know. The second gives number

More information

Simulation of Eclipsing Binary Star Systems. Abstract

Simulation of Eclipsing Binary Star Systems. Abstract Simultion of Eclipsing Binry Str Systems Boris Yim 1, Kenny Chn 1, Rphel Hui 1 Wh Yn College Kowloon Diocesn Boys School Abstrct This report briefly introduces the informtion on eclipsing binry str systems.

More information

The Algebra (al-jabr) of Matrices

The Algebra (al-jabr) of Matrices Section : Mtri lgebr nd Clculus Wshkewicz College of Engineering he lgebr (l-jbr) of Mtrices lgebr s brnch of mthemtics is much broder thn elementry lgebr ll of us studied in our high school dys. In sense

More information

Multivariate problems and matrix algebra

Multivariate problems and matrix algebra University of Ferrr Stefno Bonnini Multivrite problems nd mtrix lgebr Multivrite problems Multivrite sttisticl nlysis dels with dt contining observtions on two or more chrcteristics (vribles) ech mesured

More information

FEM ANALYSIS OF ROGOWSKI COILS COUPLED WITH BAR CONDUCTORS

FEM ANALYSIS OF ROGOWSKI COILS COUPLED WITH BAR CONDUCTORS XIX IMEKO orld Congress Fundmentl nd Applied Metrology September 6 11, 2009, Lisbon, Portugl FEM ANALYSIS OF ROGOSKI COILS COUPLED ITH BAR CONDUCTORS Mirko Mrrcci, Bernrdo Tellini, Crmine Zppcost University

More information

First Semester Review Calculus BC

First Semester Review Calculus BC First Semester Review lculus. Wht is the coordinte of the point of inflection on the grph of Multiple hoice: No lcultor y 3 3 5 4? 5 0 0 3 5 0. The grph of piecewise-liner function f, for 4, is shown below.

More information

Purpose of the experiment

Purpose of the experiment Newton s Lws II PES 6 Advnced Physics Lb I Purpose of the experiment Exmine two cses using Newton s Lws. Sttic ( = 0) Dynmic ( 0) fyi fyi Did you know tht the longest recorded flight of chicken is thirteen

More information

APPROXIMATE INTEGRATION

APPROXIMATE INTEGRATION APPROXIMATE INTEGRATION. Introduction We hve seen tht there re functions whose nti-derivtives cnnot be expressed in closed form. For these resons ny definite integrl involving these integrnds cnnot be

More information

USA Mathematical Talent Search Round 1 Solutions Year 25 Academic Year

USA Mathematical Talent Search Round 1 Solutions Year 25 Academic Year 1/1/5. Alex is trying to oen lock whose code is sequence tht is three letters long, with ech of the letters being one of A, B or C, ossibly reeted. The lock hs three buttons, lbeled A, B nd C. When the

More information

13.3 CLASSICAL STRAIGHTEDGE AND COMPASS CONSTRUCTIONS

13.3 CLASSICAL STRAIGHTEDGE AND COMPASS CONSTRUCTIONS 33 CLASSICAL STRAIGHTEDGE AND COMPASS CONSTRUCTIONS As simple ppliction of the results we hve obtined on lgebric extensions, nd in prticulr on the multiplictivity of extension degrees, we cn nswer (in

More information

Presentation Problems 5

Presentation Problems 5 Presenttion Problems 5 21-355 A For these problems, ssume ll sets re subsets of R unless otherwise specified. 1. Let P nd Q be prtitions of [, b] such tht P Q. Then U(f, P ) U(f, Q) nd L(f, P ) L(f, Q).

More information