Diffusivity Equations for Flow in Porous Media

Size: px
Start display at page:

Download "Diffusivity Equations for Flow in Porous Media"

Transcription

1 Peroleum Engineering 324 Reservoir Performance Texas A&M Universiy T.A. Blasingame, Texas A&M U. Dearmen of Peroleum Engineering Texas A&M Universiy College Saion, TX Slide 1

2 Lecure: : "Black Oil" (> b ) "Soluion-Gas Drive" (valid for all, referenced for < b ) "Dry Gas" (> d ) Mulihase Flow Slide 2

3 Diffusiviy Equaion: Black Oil (> b ) for a Black Oil: Slighly Comressible Liquid: (General Form) 2 2 c c ( ) k 2 c k Slighly Comressible Liquid: (Small and c form) Slide 3

4 Diffusiviy Equaion: Black Oil o and B o vs. Behavior of he o and B o variables as funcions of ressure for an examle black oil case. Noe behavior for > b boh variables should be considered o be "aroximaely consan" for he sake of develoing flow relaions. Such an assumion (i.e., o and B o consan) is no an absolue requiremen, bu his assumion is fundamenal for he develomen of "liquid" flow soluions in reservoir engineering. Slide 4

5 Diffusiviy Equaion: Black Oil c o vs. Behavior of he c o variable as a funcion of ressure examle black oil case. Noe he "jum" a = b, his behavior is due o he gas exansion a he bubbleoin. Slide 5

6 Diffusiviy Equaion: Soluion Gas Drive (< b ) for Soluion-Gas Drive: (< b ) Oil Pseudoressure Form: (Accouns for o and B o ) 2 o c k o Oil Pseudoressure Definiion: ( n is any reference ressure) 1 B B d o o o n base o o Slide 6

7 Diffusiviy Equaion: Soln Gas Drive 1/( o B o ) vs. 1/( o B o ) vs. ( b =5000 sia, Various T) "Soluion-Gas Drive" PVT Proeries: (1/( o B o ), < b, b =5000 sia) Aem o illusrae ha 1/( o B o ) consan for < b. This would allow us o aroximae behavior using "liquid" equaions. Slide 7

8 Diffusiviy Equaion: Soln Gas Drive 1/( o B o ) vs. 1/( o B o ) vs. ( b =5000 sia, T=175 Deg F) 1 o o Bo d n B base o o "Soluion-Gas Drive" Pseudoressure Condiion: (1/( o B o ) vs. ) Conce: IF 1/( o B o ) consan, THEN oil seudoressure NOT required. 1/( o B o ) is NEVER "consan" bu does no vary significanly wih. Oil seudoressure calculaion sraighforward, bu robably no necessary. Slide 8

9 Diffusiviy Equaion: Soln Gas Drive ( oco) vs. ( oco) vs. (b=5000 sia, Various T) "Soluion-Gas Drive" PVT Proeries: ( oco, <b, b=5000 sia) Aem o illusrae ha oco is NEVER consan. CAN NOT aroximae behavior using "liquid" equaions (or so i seems). Slide 9

10 Diffusiviy Equaion: Soln Gas Drive ( oco) vs. ( oco) vs. (b=5000 sia, T=175 Deg F) 1 a,o o c d n 0 o ( )c ( ) "Soluion-Gas Drive" Pseudoressure Condiion: (( oco) vs. ) Conce: IF ( oco) consan, THEN oil seudoime NOT required. ( oco) is NEVER "consan" BUT, oil seudoime would be very difficul. Oher evidence suggess ha ignoring ( oco) variance is acceable. Slide 10

11 Diffusiviy Equaion: Soln Gas Drive (c/ ) vs. D Dominaed Flow in Soluion-Gas-Drive Reservoirs," SPERE (November 1989) Camacho-V., R.G. and Raghavan, R.: "Boundary- Soluion-Gas Drive Mobiliy/Comressibiliy (Camacho) "Soluion-Gas Drive" Behavior: ((c/ ) vs. ime) Observaion: (c/ ) consan for >b and laer, for <b. wf = consan bu robably valid for any roducion/ressure scenario. Slide 11

12 Diffusiviy Equaion: Soln Gas Drive (2-hase ) Hisorical Noe Evinger-Muska Conce (1942) Why no use liquid seudoressure? Evinger and Muska (1942) noe ha: The indefinie inegral may be evaluaed, as was done for he wo-hase sysem, and he ressure disribuion may be deermined. However, i will be sufficien for he calculaion of he roduciviy facor o consider only he limiing form... (i.e., he consan roery liquid relaion). Slide 12

13 Diffusiviy Equaion: Dry Gas Relaions for a "Dry Gas:" General Form for Gas: c [ ] g z k z Diffusiviy Relaions: Pseudoressure/Time: g c k 2 Pseudoressure/Pseudoime: ( g c ) n k a 2 Definiions: Pseudoressure: g z g Pseudoime: d base g z n 1 a g c d n 0 g ( )c ( ) Slide 13

14 Diffusiviy Equaion: Pseudoime ( gcg vs. ) Dry Gas Pseudoime Condiion ( gcg vs., Various T) "Dry Gas" PVT Proeries: ( gcg vs. ) Conce: If gcg consan, seudoime NOT required. Readily observe ha gcg is NEVER consan, seudoime required. Slide 14

15 Diffusiviy Equaion: Pseudoime ( gcg vs. ) Dry Gas Pseudoime Condiion ( gcg vs., T=200 Deg F) 1 a g c d n 0 g ( )c ( ) "Dry Gas" Pseudoime Condiion: ( gcg vs. ) Conce: IF gcg consan, THEN seudoime NOT required. gcg is NEVER consan seudoime is always required (for liquid eq.). However, can generae numerical soluion for gas cases (no seudoime). Slide 15

16 Diffusiviy Equaion: 2 Relaions Dry Gas 2 Relaions for a "Dry Gas:" 2 Relaions 2 Form Full Formulaion: g c 2 ( ) 2 [ln( g z )] ( ) ( ) k Form Aroximaion: g c 2 ( ) ( ) k 2 2 Slide 16

17 Diffusiviy Equaion: Gas 2 ( gz vs., Various T) Dry Gas 2 Condiion ( gz vs., Various T) "Dry Gas" PVT Proeries: ( gz vs. ) Basis for he "ressure-squared" aroximaion (i.e., use of 2 variable). Conce: ( gz) = consan, valid only for <2000 sia. Slide 17

18 Diffusiviy Equaion: Gas 2 ( gz vs., T=200 Deg F) Dry Gas 2 Condiion ( gz vs., T=200 Deg F) g z g d base g z n "Dry Gas" PVT Proeries: ( gz vs. ) Conce: IF ( gz) = consan, THEN 2-variable valid. ( gz) consan for <2000 sia. Even wih numerical soluions, 2 formulaion would no be aroriae. Slide 18

19 Diffusiviy Equaion: Relaions Dry Gas Relaions for a "Dry Gas:" Relaions Form Full Formulaion: g z 2 g c ln ( ) k 2 Form Aroximaion: g c k 2 Slide 19

20 Diffusiviy Equaion: Gas (/( gz) vs., Various T) Dry Gas Condiion (/( gz) vs., Various T) "Dry Gas" PVT Proeries: (/( gz) vs. ) Basis for he "ressure" aroximaion (i.e., use of variable). Conce: (/ gz) = consan (never valid). Slide 20

21 Diffusiviy Equaion: Gas (/( gz) vs., T=200 Deg F) Dry Gas Condiion (/( gz) vs., T=200 Deg F) g z g d base g z n "Dry Gas" PVT Proeries: (/( gz) vs. ) Conce: IF /( gz) = consan, THEN -variable is valid. /( gz) is NEVER consan seudoressure required (for liquid eq.). formulaion is never aroriae (even if generaed numerically). Slide 21

22 Diffusiviy Equaion: Mulihase Relaions Mulihase Case -Form Relaions (Perrine-Marin) Gas Equaion: kg ko k w S g So Sw Rso Rsw Rso Rsw B B B B B B o o w w o w g g g Oil Equaion: k S o o B o o Bo Waer Equaion: k S w w w Bw Bw Mulihase Equaion: c ko k g k w o g w 2 Comressibiliy Terms: 1 dbo B g drso co Bo d Bo d 1 dbw B g drsw cw Bw d Bw d 1 db g c g B g d c co So cw S w c g S g c f Slide 22

23 Peroleum Engineering 324 Reservoir Performance Texas A&M Universiy (End of Lecure) T.A. Blasingame, Texas A&M U. Dearmen of Peroleum Engineering Texas A&M Universiy College Saion, TX Slide 23

Tom BLASINGAME Texas A&M U. Slide 1

Tom BLASINGAME Texas A&M U. Slide 1 Fundamenal Fl Lecure 7 The Diffusiviy Equain fr Mulihase Fl Slide 1 Fundamenal Fl Lecure 7 The Diffusiviy Equain fr Mulihase Fl Slide Fundamenal Fl Lecure 7 The Diffusiviy Equain fr Mulihase Fl Slide 3

More information

Diffusivity Equations (Governing Flow Relations)

Diffusivity Equations (Governing Flow Relations) Diffusiviy Equains (Gvernin Fl Relains) Thmas A. lasiname, Ph.D., P.E. Dearmen f Perleum Enineerin Texas A&M Universiy Cllee Sain, TX 77843-3116 (USA) +1.979.845.9 -blasiname@amu.edu Orienain Diffusiviy

More information

Diffusivity Equation

Diffusivity Equation Perleum Enineerin 34 Reservir Perfrmane Diffusiviy Equain 13 February 008 Thmas A. lasiname, Ph.D., P.E. Dilhan Il Dearmen f Perleum Enineerin Dearmen f Perleum Enineerin Texas A&M Universiy Texas A&M

More information

Module 2 F c i k c s la l w a s o s f dif di fusi s o i n

Module 2 F c i k c s la l w a s o s f dif di fusi s o i n Module Fick s laws of diffusion Fick s laws of diffusion and hin film soluion Adolf Fick (1855) proposed: d J α d d d J (mole/m s) flu (m /s) diffusion coefficien and (mole/m 3 ) concenraion of ions, aoms

More information

- The whole joint distribution is independent of the date at which it is measured and depends only on the lag.

- The whole joint distribution is independent of the date at which it is measured and depends only on the lag. Saionary Processes Sricly saionary - The whole join disribuion is indeenden of he dae a which i is measured and deends only on he lag. - E y ) is a finie consan. ( - V y ) is a finie consan. ( ( y, y s

More information

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle Chaper 2 Newonian Mechanics Single Paricle In his Chaper we will review wha Newon s laws of mechanics ell us abou he moion of a single paricle. Newon s laws are only valid in suiable reference frames,

More information

Production Analysis Reserves Estimation in Unconventional Reservoirs New Rate-Time Relations

Production Analysis Reserves Estimation in Unconventional Reservoirs New Rate-Time Relations Naural Gas Engineering Producion Analysis Reserves Esimaion in Unconvenional Reservoirs New Rae-Time Relaions T.A. Blasingame, Texas A&M U. Deparmen of Peroleum Engineering Texas A&M Universiy College

More information

20. Applications of the Genetic-Drift Model

20. Applications of the Genetic-Drift Model 0. Applicaions of he Geneic-Drif Model 1) Deermining he probabiliy of forming any paricular combinaion of genoypes in he nex generaion: Example: If he parenal allele frequencies are p 0 = 0.35 and q 0

More information

Challenge Problems. DIS 203 and 210. March 6, (e 2) k. k(k + 2). k=1. f(x) = k(k + 2) = 1 x k

Challenge Problems. DIS 203 and 210. March 6, (e 2) k. k(k + 2). k=1. f(x) = k(k + 2) = 1 x k Challenge Problems DIS 03 and 0 March 6, 05 Choose one of he following problems, and work on i in your group. Your goal is o convince me ha your answer is correc. Even if your answer isn compleely correc,

More information

ME 391 Mechanical Engineering Analysis

ME 391 Mechanical Engineering Analysis Fall 04 ME 39 Mechanical Engineering Analsis Eam # Soluions Direcions: Open noes (including course web posings). No books, compuers, or phones. An calculaor is fair game. Problem Deermine he posiion of

More information

4.6 One Dimensional Kinematics and Integration

4.6 One Dimensional Kinematics and Integration 4.6 One Dimensional Kinemaics and Inegraion When he acceleraion a( of an objec is a non-consan funcion of ime, we would like o deermine he ime dependence of he posiion funcion x( and he x -componen of

More information

Vehicle Arrival Models : Headway

Vehicle Arrival Models : Headway Chaper 12 Vehicle Arrival Models : Headway 12.1 Inroducion Modelling arrival of vehicle a secion of road is an imporan sep in raffic flow modelling. I has imporan applicaion in raffic flow simulaion where

More information

ON DETERMINATION OF SOME CHARACTERISTICS OF SEMI-MARKOV PROCESS FOR DIFFERENT DISTRIBUTIONS OF TRANSIENT PROBABILITIES ABSTRACT

ON DETERMINATION OF SOME CHARACTERISTICS OF SEMI-MARKOV PROCESS FOR DIFFERENT DISTRIBUTIONS OF TRANSIENT PROBABILITIES ABSTRACT Zajac, Budny ON DETERMINATION O SOME CHARACTERISTICS O SEMI MARKOV PROCESS OR DIERENT DISTRIBUTIONS O R&RATA # 2(3 ar 2 (Vol. 2 29, June ON DETERMINATION O SOME CHARACTERISTICS O SEMI-MARKOV PROCESS OR

More information

Predator - Prey Model Trajectories and the nonlinear conservation law

Predator - Prey Model Trajectories and the nonlinear conservation law Predaor - Prey Model Trajecories and he nonlinear conservaion law James K. Peerson Deparmen of Biological Sciences and Deparmen of Mahemaical Sciences Clemson Universiy Ocober 28, 213 Ouline Drawing Trajecories

More information

Well Performance Prediction of a Well Experiencing Changes in Completion T. Marhaendrajana, and J. Desroches, Schlumberger Oilfield Services

Well Performance Prediction of a Well Experiencing Changes in Completion T. Marhaendrajana, and J. Desroches, Schlumberger Oilfield Services SPE 75667 Well Performance Predicion of a Well Exeriencing Changes in Comleion T Marhaendrajana, and J esroches, Schlumberger Oilfield Services Coyrigh 2002, Sociey of Peroleum Engineers Inc This aer was

More information

Petroleum Engineering 324 Reservoir Performance. Material Balance 16 February 2007

Petroleum Engineering 324 Reservoir Performance. Material Balance 16 February 2007 Petroleum Engineering 324 Reservoir Performance Material Balance 16 February 2007 Thomas A. Blasingame, Ph.D., P.E. Deartment of Petroleum Engineering Texas A&M University College Station, TX 77843-3116

More information

ODEs II, Lecture 1: Homogeneous Linear Systems - I. Mike Raugh 1. March 8, 2004

ODEs II, Lecture 1: Homogeneous Linear Systems - I. Mike Raugh 1. March 8, 2004 ODEs II, Lecure : Homogeneous Linear Sysems - I Mike Raugh March 8, 4 Inroducion. In he firs lecure we discussed a sysem of linear ODEs for modeling he excreion of lead from he human body, saw how o ransform

More information

Comparison between the Discrete and Continuous Time Models

Comparison between the Discrete and Continuous Time Models Comparison beween e Discree and Coninuous Time Models D. Sulsky June 21, 2012 1 Discree o Coninuous Recall e discree ime model Î = AIS Ŝ = S Î. Tese equaions ell us ow e populaion canges from one day o

More information

The Asymptotic Behavior of Nonoscillatory Solutions of Some Nonlinear Dynamic Equations on Time Scales

The Asymptotic Behavior of Nonoscillatory Solutions of Some Nonlinear Dynamic Equations on Time Scales Advances in Dynamical Sysems and Applicaions. ISSN 0973-5321 Volume 1 Number 1 (2006, pp. 103 112 c Research India Publicaions hp://www.ripublicaion.com/adsa.hm The Asympoic Behavior of Nonoscillaory Soluions

More information

Linear Response Theory: The connection between QFT and experiments

Linear Response Theory: The connection between QFT and experiments Phys540.nb 39 3 Linear Response Theory: The connecion beween QFT and experimens 3.1. Basic conceps and ideas Q: How do we measure he conduciviy of a meal? A: we firs inroduce a weak elecric field E, and

More information

MA 214 Calculus IV (Spring 2016) Section 2. Homework Assignment 1 Solutions

MA 214 Calculus IV (Spring 2016) Section 2. Homework Assignment 1 Solutions MA 14 Calculus IV (Spring 016) Secion Homework Assignmen 1 Soluions 1 Boyce and DiPrima, p 40, Problem 10 (c) Soluion: In sandard form he given firs-order linear ODE is: An inegraing facor is given by

More information

Lecture 20: Riccati Equations and Least Squares Feedback Control

Lecture 20: Riccati Equations and Least Squares Feedback Control 34-5 LINEAR SYSTEMS Lecure : Riccai Equaions and Leas Squares Feedback Conrol 5.6.4 Sae Feedback via Riccai Equaions A recursive approach in generaing he marix-valued funcion W ( ) equaion for i for he

More information

Section 4.4 Logarithmic Properties

Section 4.4 Logarithmic Properties Secion. Logarihmic Properies 59 Secion. Logarihmic Properies In he previous secion, we derived wo imporan properies of arihms, which allowed us o solve some asic eponenial and arihmic equaions. Properies

More information

23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes

23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes Represening Periodic Funcions by Fourier Series 3. Inroducion In his Secion we show how a periodic funcion can be expressed as a series of sines and cosines. We begin by obaining some sandard inegrals

More information

MAE143A Signals & Systems - Homework 2, Winter 2014 due by the end of class Thursday January 23, 2014.

MAE143A Signals & Systems - Homework 2, Winter 2014 due by the end of class Thursday January 23, 2014. MAE43A Signals & Sysems - Homework, Winer 4 due by he end of class Thursday January 3, 4. Quesion Zener diode malab [Chaparro Quesion.] A zener diode circui is such ha an oupu corresponding o an inpu v

More information

1 1 + x 2 dx. tan 1 (2) = ] ] x 3. Solution: Recall that the given integral is improper because. x 3. 1 x 3. dx = lim dx.

1 1 + x 2 dx. tan 1 (2) = ] ] x 3. Solution: Recall that the given integral is improper because. x 3. 1 x 3. dx = lim dx. . Use Simpson s rule wih n 4 o esimae an () +. Soluion: Since we are using 4 seps, 4 Thus we have [ ( ) f() + 4f + f() + 4f 3 [ + 4 4 6 5 + + 4 4 3 + ] 5 [ + 6 6 5 + + 6 3 + ]. 5. Our funcion is f() +.

More information

ln 2 1 ln y x c y C x

ln 2 1 ln y x c y C x Lecure 14 Appendi B: Some sample problems from Boas Here are some soluions o he sample problems assigned for Chaper 8 8: 6 Soluion: We wan o find he soluion o he following firs order equaion using separaion

More information

4.5 Constant Acceleration

4.5 Constant Acceleration 4.5 Consan Acceleraion v() v() = v 0 + a a() a a() = a v 0 Area = a (a) (b) Figure 4.8 Consan acceleraion: (a) velociy, (b) acceleraion When he x -componen of he velociy is a linear funcion (Figure 4.8(a)),

More information

Unsteady Mass- Transfer Models

Unsteady Mass- Transfer Models See T&K Chaper 9 Unseady Mass- Transfer Models ChEn 6603 Wednesday, April 4, Ouline Conex for he discussion Soluion for ransien binary diffusion wih consan c, N. Soluion for mulicomponen diffusion wih

More information

Chapter 2. First Order Scalar Equations

Chapter 2. First Order Scalar Equations Chaper. Firs Order Scalar Equaions We sar our sudy of differenial equaions in he same way he pioneers in his field did. We show paricular echniques o solve paricular ypes of firs order differenial equaions.

More information

14 Autoregressive Moving Average Models

14 Autoregressive Moving Average Models 14 Auoregressive Moving Average Models In his chaper an imporan parameric family of saionary ime series is inroduced, he family of he auoregressive moving average, or ARMA, processes. For a large class

More information

Final Spring 2007

Final Spring 2007 .615 Final Spring 7 Overview The purpose of he final exam is o calculae he MHD β limi in a high-bea oroidal okamak agains he dangerous n = 1 exernal ballooning-kink mode. Effecively, his corresponds o

More information

1 Review of Zero-Sum Games

1 Review of Zero-Sum Games COS 5: heoreical Machine Learning Lecurer: Rob Schapire Lecure #23 Scribe: Eugene Brevdo April 30, 2008 Review of Zero-Sum Games Las ime we inroduced a mahemaical model for wo player zero-sum games. Any

More information

MATH 31B: MIDTERM 2 REVIEW. x 2 e x2 2x dx = 1. ue u du 2. x 2 e x2 e x2] + C 2. dx = x ln(x) 2 2. ln x dx = x ln x x + C. 2, or dx = 2u du.

MATH 31B: MIDTERM 2 REVIEW. x 2 e x2 2x dx = 1. ue u du 2. x 2 e x2 e x2] + C 2. dx = x ln(x) 2 2. ln x dx = x ln x x + C. 2, or dx = 2u du. MATH 3B: MIDTERM REVIEW JOE HUGHES. Inegraion by Pars. Evaluae 3 e. Soluion: Firs make he subsiuion u =. Then =, hence 3 e = e = ue u Now inegrae by pars o ge ue u = ue u e u + C and subsiue he definiion

More information

Sterilization D Values

Sterilization D Values Seriliaion D Values Seriliaion by seam consis of he simple observaion ha baceria die over ime during exposure o hea. They do no all live for a finie period of hea exposure and hen suddenly die a once,

More information

Lecture #31, 32: The Ornstein-Uhlenbeck Process as a Model of Volatility

Lecture #31, 32: The Ornstein-Uhlenbeck Process as a Model of Volatility Saisics 441 (Fall 214) November 19, 21, 214 Prof Michael Kozdron Lecure #31, 32: The Ornsein-Uhlenbeck Process as a Model of Volailiy The Ornsein-Uhlenbeck process is a di usion process ha was inroduced

More information

Math 2142 Exam 1 Review Problems. x 2 + f (0) 3! for the 3rd Taylor polynomial at x = 0. To calculate the various quantities:

Math 2142 Exam 1 Review Problems. x 2 + f (0) 3! for the 3rd Taylor polynomial at x = 0. To calculate the various quantities: Mah 4 Eam Review Problems Problem. Calculae he 3rd Taylor polynomial for arcsin a =. Soluion. Le f() = arcsin. For his problem, we use he formula f() + f () + f ()! + f () 3! for he 3rd Taylor polynomial

More information

Module 3: The Damped Oscillator-II Lecture 3: The Damped Oscillator-II

Module 3: The Damped Oscillator-II Lecture 3: The Damped Oscillator-II Module 3: The Damped Oscillaor-II Lecure 3: The Damped Oscillaor-II 3. Over-damped Oscillaions. This refers o he siuaion where β > ω (3.) The wo roos are and α = β + α 2 = β β 2 ω 2 = (3.2) β 2 ω 2 = 2

More information

EEEB113 CIRCUIT ANALYSIS I

EEEB113 CIRCUIT ANALYSIS I 9/14/29 1 EEEB113 CICUIT ANALYSIS I Chaper 7 Firs-Order Circuis Maerials from Fundamenals of Elecric Circuis 4e, Alexander Sadiku, McGraw-Hill Companies, Inc. 2 Firs-Order Circuis -Chaper 7 7.2 The Source-Free

More information

2.160 System Identification, Estimation, and Learning. Lecture Notes No. 8. March 6, 2006

2.160 System Identification, Estimation, and Learning. Lecture Notes No. 8. March 6, 2006 2.160 Sysem Idenificaion, Esimaion, and Learning Lecure Noes No. 8 March 6, 2006 4.9 Eended Kalman Filer In many pracical problems, he process dynamics are nonlinear. w Process Dynamics v y u Model (Linearized)

More information

Section 4.4 Logarithmic Properties

Section 4.4 Logarithmic Properties Secion. Logarihmic Properies 5 Secion. Logarihmic Properies In he previous secion, we derived wo imporan properies of arihms, which allowed us o solve some asic eponenial and arihmic equaions. Properies

More information

Math 36. Rumbos Spring Solutions to Assignment #6. 1. Suppose the growth of a population is governed by the differential equation.

Math 36. Rumbos Spring Solutions to Assignment #6. 1. Suppose the growth of a population is governed by the differential equation. Mah 36. Rumbos Spring 1 1 Soluions o Assignmen #6 1. Suppose he growh of a populaion is governed by he differenial equaion where k is a posiive consan. d d = k (a Explain why his model predics ha he populaion

More information

Cash Flow Valuation Mode Lin Discrete Time

Cash Flow Valuation Mode Lin Discrete Time IOSR Journal of Mahemaics (IOSR-JM) e-issn: 2278-5728,p-ISSN: 2319-765X, 6, Issue 6 (May. - Jun. 2013), PP 35-41 Cash Flow Valuaion Mode Lin Discree Time Olayiwola. M. A. and Oni, N. O. Deparmen of Mahemaics

More information

Nature Neuroscience: doi: /nn Supplementary Figure 1. Spike-count autocorrelations in time.

Nature Neuroscience: doi: /nn Supplementary Figure 1. Spike-count autocorrelations in time. Supplemenary Figure 1 Spike-coun auocorrelaions in ime. Normalized auocorrelaion marices are shown for each area in a daase. The marix shows he mean correlaion of he spike coun in each ime bin wih he spike

More information

SOLUTIONS TO ECE 3084

SOLUTIONS TO ECE 3084 SOLUTIONS TO ECE 384 PROBLEM 2.. For each sysem below, specify wheher or no i is: (i) memoryless; (ii) causal; (iii) inverible; (iv) linear; (v) ime invarian; Explain your reasoning. If he propery is no

More information

Sample Autocorrelations for Financial Time Series Models. Richard A. Davis Colorado State University Thomas Mikosch University of Copenhagen

Sample Autocorrelations for Financial Time Series Models. Richard A. Davis Colorado State University Thomas Mikosch University of Copenhagen Sample Auocorrelaions for Financial Time Series Models Richard A. Davis Colorado Sae Universiy Thomas Mikosch Universiy of Copenhagen Ouline Characerisics of some financial ime series IBM reurns NZ-USA

More information

PET467E-Analysis of Well Pressure Tests/2008 Spring Semester/İTÜ Midterm Examination (Duration 3:00 hours) Solutions

PET467E-Analysis of Well Pressure Tests/2008 Spring Semester/İTÜ Midterm Examination (Duration 3:00 hours) Solutions M. Onur 03.04.008 PET467E-Analysis of Well Pressure Tess/008 Spring Semeser/İTÜ Miderm Examinaion (Duraion 3:00 hours) Soluions Name of he Suden: Insrucions: Before saring he exam, wrie your name clearly

More information

Linear Surface Gravity Waves 3., Dispersion, Group Velocity, and Energy Propagation

Linear Surface Gravity Waves 3., Dispersion, Group Velocity, and Energy Propagation Chaper 4 Linear Surface Graviy Waves 3., Dispersion, Group Velociy, and Energy Propagaion 4. Descripion In many aspecs of wave evoluion, he concep of group velociy plays a cenral role. Mos people now i

More information

Excel-Based Solution Method For The Optimal Policy Of The Hadley And Whittin s Exact Model With Arma Demand

Excel-Based Solution Method For The Optimal Policy Of The Hadley And Whittin s Exact Model With Arma Demand Excel-Based Soluion Mehod For The Opimal Policy Of The Hadley And Whiin s Exac Model Wih Arma Demand Kal Nami School of Business and Economics Winson Salem Sae Universiy Winson Salem, NC 27110 Phone: (336)750-2338

More information

d 1 = c 1 b 2 - b 1 c 2 d 2 = c 1 b 3 - b 1 c 3

d 1 = c 1 b 2 - b 1 c 2 d 2 = c 1 b 3 - b 1 c 3 and d = c b - b c c d = c b - b c c This process is coninued unil he nh row has been compleed. The complee array of coefficiens is riangular. Noe ha in developing he array an enire row may be divided or

More information

Chapter 5: Control Volume Approach and Continuity Principle Dr Ali Jawarneh

Chapter 5: Control Volume Approach and Continuity Principle Dr Ali Jawarneh Chaper 5: Conrol Volume Approach and Coninuiy Principle By Dr Ali Jawarneh Deparmen of Mechanical Engineering Hashemie Universiy 1 Ouline Rae of Flow Conrol volume approach. Conservaion of mass he coninuiy

More information

F This leads to an unstable mode which is not observable at the output thus cannot be controlled by feeding back.

F This leads to an unstable mode which is not observable at the output thus cannot be controlled by feeding back. Lecure 8 Las ime: Semi-free configuraion design This is equivalen o: Noe ns, ener he sysem a he same place. is fixed. We design C (and perhaps B. We mus sabilize if i is given as unsable. Cs ( H( s = +

More information

PVT Concepts (Reservoir Fluids)

PVT Concepts (Reservoir Fluids) PVT Concepts (Reservoir Fluids) Thomas A. Blasingame, Ph.D., P.E. Department of Petroleum Engineering Texas A&M University College Station, TX 77843-3116 (USA) +1.979.845.2292 t-blasingame@tamu.edu Orientation

More information

Object tracking: Using HMMs to estimate the geographical location of fish

Object tracking: Using HMMs to estimate the geographical location of fish Objec racking: Using HMMs o esimae he geographical locaion of fish 02433 - Hidden Markov Models Marin Wæver Pedersen, Henrik Madsen Course week 13 MWP, compiled June 8, 2011 Objecive: Locae fish from agging

More information

LAPLACE TRANSFORM AND TRANSFER FUNCTION

LAPLACE TRANSFORM AND TRANSFER FUNCTION CHBE320 LECTURE V LAPLACE TRANSFORM AND TRANSFER FUNCTION Professor Dae Ryook Yang Spring 2018 Dep. of Chemical and Biological Engineering 5-1 Road Map of he Lecure V Laplace Transform and Transfer funcions

More information

Fault Tolerant Computing CS 530 Reliability Analysis. Yashwant K. Malaiya Colorado State University

Fault Tolerant Computing CS 530 Reliability Analysis. Yashwant K. Malaiya Colorado State University CS 53 eliabiliy Analysis Yashwan K. Malaiya Colorado Sae Universiy February 6, 27 eliabiliy Analysis: Ouline eliabiliy measures: eliabiliy, availabiliy, Transacion eliabiliy, MTTF and (, MTBF Basic Cases

More information

The Quantum Theory of Atoms and Molecules: The Schrodinger equation. Hilary Term 2008 Dr Grant Ritchie

The Quantum Theory of Atoms and Molecules: The Schrodinger equation. Hilary Term 2008 Dr Grant Ritchie e Quanum eory of Aoms and Molecules: e Scrodinger equaion Hilary erm 008 Dr Gran Ricie An equaion for maer waves? De Broglie posulaed a every paricles as an associaed wave of waveleng: / p Wave naure of

More information

EECE 301 Signals & Systems Prof. Mark Fowler

EECE 301 Signals & Systems Prof. Mark Fowler EECE 3 Signals & Sysems Prof. Mark Fowler Noe Se # Wha are Coninuous-Time Signals??? /6 Coninuous-Time Signal Coninuous Time (C-T) Signal: A C-T signal is defined on he coninuum of ime values. Tha is:

More information

Pseudosteady-State Flow Relations for a Radial System from Department of Petroleum Engineering Course Notes (1997)

Pseudosteady-State Flow Relations for a Radial System from Department of Petroleum Engineering Course Notes (1997) Pseudoseady-Sae Flow Relaions fo a Radial Sysem fom Deamen of Peoleum Engineeing Couse Noes (1997) (Deivaion of he Pseudoseady-Sae Flow Relaions fo a Radial Sysem) (Deivaion of he Pseudoseady-Sae Flow

More information

Solution of Integro-Differential Equations by Using ELzaki Transform

Solution of Integro-Differential Equations by Using ELzaki Transform Global Journal of Mahemaical Sciences: Theory and Pracical. Volume, Number (), pp. - Inernaional Research Publicaion House hp://www.irphouse.com Soluion of Inegro-Differenial Equaions by Using ELzaki Transform

More information

t + t sin t t cos t sin t. t cos t sin t dt t 2 = exp 2 log t log(t cos t sin t) = Multiplying by this factor and then integrating, we conclude that

t + t sin t t cos t sin t. t cos t sin t dt t 2 = exp 2 log t log(t cos t sin t) = Multiplying by this factor and then integrating, we conclude that ODEs, Homework #4 Soluions. Check ha y ( = is a soluion of he second-order ODE ( cos sin y + y sin y sin = 0 and hen use his fac o find all soluions of he ODE. When y =, we have y = and also y = 0, so

More information

MATHEMATICAL DESCRIPTION OF THEORETICAL METHODS OF RESERVE ECONOMY OF CONSIGNMENT STORES

MATHEMATICAL DESCRIPTION OF THEORETICAL METHODS OF RESERVE ECONOMY OF CONSIGNMENT STORES MAHEMAICAL DESCIPION OF HEOEICAL MEHODS OF ESEVE ECONOMY OF CONSIGNMEN SOES Péer elek, József Cselényi, György Demeer Universiy of Miskolc, Deparmen of Maerials Handling and Logisics Absrac: Opimizaion

More information

8. Basic RL and RC Circuits

8. Basic RL and RC Circuits 8. Basic L and C Circuis This chaper deals wih he soluions of he responses of L and C circuis The analysis of C and L circuis leads o a linear differenial equaion This chaper covers he following opics

More information

3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon

3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon 3..3 INRODUCION O DYNAMIC OPIMIZAION: DISCREE IME PROBLEMS A. he Hamilonian and Firs-Order Condiions in a Finie ime Horizon Define a new funcion, he Hamilonian funcion, H. H he change in he oal value of

More information

System of Linear Differential Equations

System of Linear Differential Equations Sysem of Linear Differenial Equaions In "Ordinary Differenial Equaions" we've learned how o solve a differenial equaion for a variable, such as: y'k5$e K2$x =0 solve DE yx = K 5 2 ek2 x C_C1 2$y''C7$y

More information

Navneet Saini, Mayank Goyal, Vishal Bansal (2013); Term Project AML310; Indian Institute of Technology Delhi

Navneet Saini, Mayank Goyal, Vishal Bansal (2013); Term Project AML310; Indian Institute of Technology Delhi Creep in Viscoelasic Subsances Numerical mehods o calculae he coefficiens of he Prony equaion using creep es daa and Herediary Inegrals Mehod Navnee Saini, Mayank Goyal, Vishal Bansal (23); Term Projec

More information

Richard A. Davis Colorado State University Bojan Basrak Eurandom Thomas Mikosch University of Groningen

Richard A. Davis Colorado State University Bojan Basrak Eurandom Thomas Mikosch University of Groningen Mulivariae Regular Variaion wih Applicaion o Financial Time Series Models Richard A. Davis Colorado Sae Universiy Bojan Basrak Eurandom Thomas Mikosch Universiy of Groningen Ouline + Characerisics of some

More information

Solutions: Wednesday, November 14

Solutions: Wednesday, November 14 Amhers College Deparmen of Economics Economics 360 Fall 2012 Soluions: Wednesday, November 14 Judicial Daa: Cross secion daa of judicial and economic saisics for he fify saes in 2000. JudExp CrimesAll

More information

The motions of the celt on a horizontal plane with viscous friction

The motions of the celt on a horizontal plane with viscous friction The h Join Inernaional Conference on Mulibody Sysem Dynamics June 8, 18, Lisboa, Porugal The moions of he cel on a horizonal plane wih viscous fricion Maria A. Munisyna 1 1 Moscow Insiue of Physics and

More information

FORECASTS GENERATING FOR ARCH-GARCH PROCESSES USING THE MATLAB PROCEDURES

FORECASTS GENERATING FOR ARCH-GARCH PROCESSES USING THE MATLAB PROCEDURES FORECASS GENERAING FOR ARCH-GARCH PROCESSES USING HE MALAB PROCEDURES Dušan Marček, Insiue of Comuer Science, Faculy of Philosohy and Science, he Silesian Universiy Oava he Faculy of Managemen Science

More information

Math 2214 Solution Test 1B Fall 2017

Math 2214 Solution Test 1B Fall 2017 Mah 14 Soluion Tes 1B Fall 017 Problem 1: A ank has a capaci for 500 gallons and conains 0 gallons of waer wih lbs of sal iniiall. A soluion conaining of 8 lbsgal of sal is pumped ino he ank a 10 galsmin.

More information

IB Physics Kinematics Worksheet

IB Physics Kinematics Worksheet IB Physics Kinemaics Workshee Wrie full soluions and noes for muliple choice answers. Do no use a calculaor for muliple choice answers. 1. Which of he following is a correc definiion of average acceleraion?

More information

fakultät für informatik informatik 12 technische universität dortmund Petri Nets Peter Marwedel TU Dortmund, Informatik /10/10

fakultät für informatik informatik 12 technische universität dortmund Petri Nets Peter Marwedel TU Dortmund, Informatik /10/10 2 Peri Nes Peer Marwedel TU Dormund, Informaik 2 2008/0/0 Grahics: Alexandra Nole, Gesine Marwedel, 2003 Generalizaion of daa flow: Comuaional grahs Examle: Peri nes Inroduced in 962 by Carl Adam Peri

More information

Solutions from Chapter 9.1 and 9.2

Solutions from Chapter 9.1 and 9.2 Soluions from Chaper 9 and 92 Secion 9 Problem # This basically boils down o an exercise in he chain rule from calculus We are looking for soluions of he form: u( x) = f( k x c) where k x R 3 and k is

More information

Guest Lectures for Dr. MacFarlane s EE3350 Part Deux

Guest Lectures for Dr. MacFarlane s EE3350 Part Deux Gues Lecures for Dr. MacFarlane s EE3350 Par Deux Michael Plane Mon., 08-30-2010 Wrie name in corner. Poin ou his is a review, so I will go faser. Remind hem o go lisen o online lecure abou geing an A

More information

Solutions of Sample Problems for Third In-Class Exam Math 246, Spring 2011, Professor David Levermore

Solutions of Sample Problems for Third In-Class Exam Math 246, Spring 2011, Professor David Levermore Soluions of Sample Problems for Third In-Class Exam Mah 6, Spring, Professor David Levermore Compue he Laplace ransform of f e from is definiion Soluion The definiion of he Laplace ransform gives L[f]s

More information

The Contradiction within Equations of Motion with Constant Acceleration

The Contradiction within Equations of Motion with Constant Acceleration The Conradicion wihin Equaions of Moion wih Consan Acceleraion Louai Hassan Elzein Basheir (Daed: July 7, 0 This paper is prepared o demonsrae he violaion of rules of mahemaics in he algebraic derivaion

More information

Chapter 3 Common Families of Distributions

Chapter 3 Common Families of Distributions Chaer 3 Common Families of Disribuions Secion 31 - Inroducion Purose of his Chaer: Caalog many of common saisical disribuions (families of disribuions ha are indeed by one or more arameers) Wha we should

More information

Homework 10 (Stats 620, Winter 2017) Due Tuesday April 18, in class Questions are derived from problems in Stochastic Processes by S. Ross.

Homework 10 (Stats 620, Winter 2017) Due Tuesday April 18, in class Questions are derived from problems in Stochastic Processes by S. Ross. Homework (Sas 6, Winer 7 Due Tuesday April 8, in class Quesions are derived from problems in Sochasic Processes by S. Ross.. A sochasic process {X(, } is said o be saionary if X(,..., X( n has he same

More information

Kinematics Vocabulary. Kinematics and One Dimensional Motion. Position. Coordinate System in One Dimension. Kinema means movement 8.

Kinematics Vocabulary. Kinematics and One Dimensional Motion. Position. Coordinate System in One Dimension. Kinema means movement 8. Kinemaics Vocabulary Kinemaics and One Dimensional Moion 8.1 WD1 Kinema means movemen Mahemaical descripion of moion Posiion Time Inerval Displacemen Velociy; absolue value: speed Acceleraion Averages

More information

Math Final Exam Solutions

Math Final Exam Solutions Mah 246 - Final Exam Soluions Friday, July h, 204 () Find explici soluions and give he inerval of definiion o he following iniial value problems (a) ( + 2 )y + 2y = e, y(0) = 0 Soluion: In normal form,

More information

KEY. Math 334 Midterm I Fall 2008 sections 001 and 003 Instructor: Scott Glasgow

KEY. Math 334 Midterm I Fall 2008 sections 001 and 003 Instructor: Scott Glasgow 1 KEY Mah 4 Miderm I Fall 8 secions 1 and Insrucor: Sco Glasgow Please do NOT wrie on his eam. No credi will be given for such work. Raher wrie in a blue book, or on our own paper, preferabl engineering

More information

Mathcad Lecture #7 In-class Worksheet "Smart" Solve Block Techniques Handout

Mathcad Lecture #7 In-class Worksheet Smart Solve Block Techniques Handout Mahcad Lecure #7 In-class Workshee "Smar" Solve Block echniques Handou A he end of his lecure, you will be able o: use funcions in solve block equaions o improve convergence consruc solve blocks wih minimal

More information

Diebold, Chapter 7. Francis X. Diebold, Elements of Forecasting, 4th Edition (Mason, Ohio: Cengage Learning, 2006). Chapter 7. Characterizing Cycles

Diebold, Chapter 7. Francis X. Diebold, Elements of Forecasting, 4th Edition (Mason, Ohio: Cengage Learning, 2006). Chapter 7. Characterizing Cycles Diebold, Chaper 7 Francis X. Diebold, Elemens of Forecasing, 4h Ediion (Mason, Ohio: Cengage Learning, 006). Chaper 7. Characerizing Cycles Afer compleing his reading you should be able o: Define covariance

More information

Hamilton- J acobi Equation: Weak S olution We continue the study of the Hamilton-Jacobi equation:

Hamilton- J acobi Equation: Weak S olution We continue the study of the Hamilton-Jacobi equation: M ah 5 7 Fall 9 L ecure O c. 4, 9 ) Hamilon- J acobi Equaion: Weak S oluion We coninue he sudy of he Hamilon-Jacobi equaion: We have shown ha u + H D u) = R n, ) ; u = g R n { = }. ). In general we canno

More information

An Analysis of the Impact of Valve Closure Time on the Course of Water Hammer

An Analysis of the Impact of Valve Closure Time on the Course of Water Hammer Archives of Hydro-Engineering and Environmenal Mechanics Vol. 63 (216), No. 1,. 35 45 DOI: 1.1515/heem-216-3 IBW PAN, ISSN 1231 3726 An Analysis of he Imac of Valve Closure Time on he Course of Waer Hammer

More information

Most Probable Phase Portraits of Stochastic Differential Equations and Its Numerical Simulation

Most Probable Phase Portraits of Stochastic Differential Equations and Its Numerical Simulation Mos Probable Phase Porrais of Sochasic Differenial Equaions and Is Numerical Simulaion Bing Yang, Zhu Zeng and Ling Wang 3 School of Mahemaics and Saisics, Huazhong Universiy of Science and Technology,

More information

( ) = Q 0. ( ) R = R dq. ( t) = I t

( ) = Q 0. ( ) R = R dq. ( t) = I t ircuis onceps The addiion of a simple capacior o a circui of resisors allows wo relaed phenomena o occur The observaion ha he ime-dependence of a complex waveform is alered by he circui is referred o as

More information

THE 2-BODY PROBLEM. FIGURE 1. A pair of ellipses sharing a common focus. (c,b) c+a ROBERT J. VANDERBEI

THE 2-BODY PROBLEM. FIGURE 1. A pair of ellipses sharing a common focus. (c,b) c+a ROBERT J. VANDERBEI THE 2-BODY PROBLEM ROBERT J. VANDERBEI ABSTRACT. In his shor noe, we show ha a pair of ellipses wih a common focus is a soluion o he 2-body problem. INTRODUCTION. Solving he 2-body problem from scrach

More information

The general Solow model

The general Solow model The general Solow model Back o a closed economy In he basic Solow model: no growh in GDP per worker in seady sae This conradics he empirics for he Wesern world (sylized fac #5) In he general Solow model:

More information

Math 2214 Solution Test 1A Spring 2016

Math 2214 Solution Test 1A Spring 2016 Mah 14 Soluion Tes 1A Spring 016 sec Problem 1: Wha is he larges -inerval for which ( 4) = has a guaraneed + unique soluion for iniial value (-1) = 3 according o he Exisence Uniqueness Theorem? Soluion

More information

Solutions to Assignment 1

Solutions to Assignment 1 MA 2326 Differenial Equaions Insrucor: Peronela Radu Friday, February 8, 203 Soluions o Assignmen. Find he general soluions of he following ODEs: (a) 2 x = an x Soluion: I is a separable equaion as we

More information

BOX-JENKINS MODEL NOTATION. The Box-Jenkins ARMA(p,q) model is denoted by the equation. pwhile the moving average (MA) part of the model is θ1at

BOX-JENKINS MODEL NOTATION. The Box-Jenkins ARMA(p,q) model is denoted by the equation. pwhile the moving average (MA) part of the model is θ1at BOX-JENKINS MODEL NOAION he Box-Jenkins ARMA(,q) model is denoed b he equaion + + L+ + a θ a L θ a 0 q q. () he auoregressive (AR) ar of he model is + L+ while he moving average (MA) ar of he model is

More information

Math 315: Linear Algebra Solutions to Assignment 6

Math 315: Linear Algebra Solutions to Assignment 6 Mah 35: Linear Algebra s o Assignmen 6 # Which of he following ses of vecors are bases for R 2? {2,, 3, }, {4,, 7, 8}, {,,, 3}, {3, 9, 4, 2}. Explain your answer. To generae he whole R 2, wo linearly independen

More information

Sections 2.2 & 2.3 Limit of a Function and Limit Laws

Sections 2.2 & 2.3 Limit of a Function and Limit Laws Mah 80 www.imeodare.com Secions. &. Limi of a Funcion and Limi Laws In secion. we saw how is arise when we wan o find he angen o a curve or he velociy of an objec. Now we urn our aenion o is in general

More information

f(s)dw Solution 1. Approximate f by piece-wise constant left-continuous non-random functions f n such that (f(s) f n (s)) 2 ds 0.

f(s)dw Solution 1. Approximate f by piece-wise constant left-continuous non-random functions f n such that (f(s) f n (s)) 2 ds 0. Advanced Financial Models Example shee 3 - Michaelmas 217 Michael Tehranchi Problem 1. Le f : [, R be a coninuous (non-random funcion and W a Brownian moion, and le σ 2 = f(s 2 ds and assume σ 2

More information

L07. KALMAN FILTERING FOR NON-LINEAR SYSTEMS. NA568 Mobile Robotics: Methods & Algorithms

L07. KALMAN FILTERING FOR NON-LINEAR SYSTEMS. NA568 Mobile Robotics: Methods & Algorithms L07. KALMAN FILTERING FOR NON-LINEAR SYSTEMS NA568 Mobile Roboics: Mehods & Algorihms Today s Topic Quick review on (Linear) Kalman Filer Kalman Filering for Non-Linear Sysems Exended Kalman Filer (EKF)

More information

Stochastic Reservoir Systems with Different Assumptions for Storage Losses

Stochastic Reservoir Systems with Different Assumptions for Storage Losses American Journal of Operaions Research, 26, 6, 44-423 hp://www.scirp.org/journal/ajor ISSN Online: 26-8849 ISSN Prin: 26-883 Sochasic Reservoir Ssems wih Differen Assumpions for Sorage Losses Carer Browning,

More information

2 int T. is the Fourier transform of f(t) which is the inverse Fourier transform of f. i t e

2 int T. is the Fourier transform of f(t) which is the inverse Fourier transform of f. i t e PHYS67 Class 3 ourier Transforms In he limi T, he ourier series becomes an inegral ( nt f in T ce f n f f e d, has been replaced by ) where i f e d is he ourier ransform of f() which is he inverse ourier

More information

ψ(t) = V x (0)V x (t)

ψ(t) = V x (0)V x (t) .93 Home Work Se No. (Professor Sow-Hsin Chen Spring Term 5. Due March 7, 5. This problem concerns calculaions of analyical expressions for he self-inermediae scaering funcion (ISF of he es paricle in

More information