A Method for Setting the Artificial Boundary Conditions of Groundwater Model
|
|
- Hector Houston
- 5 years ago
- Views:
Transcription
1 Open Journal of Geology, 2013, 3, doi: /ojg b012 Published Online April 2013 (hp:// A Mehod for Seing he Arificial Boundary Condiions of Groundwaer Model Yipeng Zhou 1,2, Zhaoli Shen 1, Weijun Shi 2, Jinhui Liu 2, Yajie Liu 2 1 School of Waer Resources and Environmen, China Universiy of Geosciences, Beijing, China 2 Deparmen of Civil and Environmenal Engineering, Eas China Insiue of Technology, Nanchang, China zyp721@163.com, wjshieci@163.com Received 2013 ABSTRACT Numerical simulaion echnology is nowadays an imporan means for groundwaer issues because of is efficiency and economical advanages. Bu in case of naural hydrogeological boundaries are no wihin he ineres area, i may be a big rouble o se boundary condiions of he model arificially wihou enough field invesigaion informaion. This paper inroduced a mehod for solving such problem applying field pumping es and recovery es. The mehod was applied o build an in-siu leaching of uranium model. Resuls showed ha he model boundary condiions can be se saisfacorily, and also he calculaed heads mached he observed daa well in boh wo models. Keywords: Groundwaer; Numerical Model; Arificial Boundary Condiions 1. Inroducion Since he 1960s, wih he developmen of compuer echnology, he mehod of numerical simulaion had been widely used o solve groundwaer flow and solue ranspor problems because of is effeciveness, flexibiliy and relaively economical wih spend, and gradually become an imporan mehod for groundwaer issues [1-7]. However, alhough los of models have been buil in various applicaions, few people care abou he real effecs of hose models in pracices [8]. One of imporan facors influencing he reliabiliy of he groundwaer model is geology and hydrogeology invesigaion; and usually making reasonable undersandings on boundary condiions is a big challenge [9]. Once he boundary condiions are disored o he ruh, i is bound o lead o significan deviaion of he model calibraion parameers from acual values, and hen serious impac on he reliabiliy of he model would no be avoided. Model boundary condiions are usually se according o field invesigaions. When he ineres area is small ha he model boundaries are far away from naural hydrogeological boundaries, arificial model boundary condiions need o be se according o a long-erm observaion of groundwaer a hose boundaries. However, in many cases, he required observaion daa are ofen unavailable; in his dilemma, one alernaive way is o expand model exen so ha he groundwaer can be assumed no o be affeced by human aciviies (such as pumping es) aken placed in he ineres area; hen he boundary condiions of firs ype or of second ype can be se a model s boundaries [9]. Bu his kind of soluion also has shorcomings, one of which is ha o esablish of model hydrogeology configuraion beyond he ineres area wihou supplemenary geology invesigaion informaion may bring unexpeced serious error o he simulaion resuls [10,11]. In his paper, in order o build a flow model of groundwaer and leaching soluion during insiu leaching of uranium process, a mehod has been employed o se arificial model boundaries by combining heoreical calculaion according groundwaer unseady flow heory and he model ieraive calibraion using observaion daa of pumping es and recovery es independenly. 2. Mehods 2.1. Basic Principles The basic principle is using field pumping es and recovery es o calibrae he model parameers and boundary condiions simulaneously. Firs, ge he head drawdown funcion derived from Jacob formula of groundwaer unseady flow a he boundaries locaed wihin he cone depression; Then se iniial heads generalized from he head drawdown funcion of he model boundaries for model building and calibraion; finally, calibrae model parameers and boundary condiions ieraively using he observaion daa of pumping es and recovery es independenly and hen make he resuls fi he facs o he mos degree.
2 Y. P. ZHOU ET AL Funcions of he Model Boundary Heads During pumping es in he confined aquifer, draw down lead o he formaion of he cone of depression of pressure head. The range of cone expands coninuously wih pumping, and gradually achieves a relaive sable sae. When he whole model is locaed wihin he cone, heads a model boundaries would vary wih ime; so obviously, model boundary condiions need o be se according o head changes. Therefore, he funcion of head variaion mus be go a firs. According o he heory of confined waer s unseady flow owards o fully peneraing well, he variaion of he head drawdown wihin pumping influence scope can be approximaely described by Jacob Formula [12], as in Q 2.25T s ln (1) 2 * 4T r where s is he head drawdown; Q is he pumping rae; T is he coefficien of ransmissibiliy; is he pumping ime; r is he disance o pumping well; μ * is he coefficien of sorage. Thus, he live head can be calculaed via Equaion (2): H (,) r Hr (, ) sr (,) (2) S where H (,) r S is he head a he poin wih he disance of r o he pumping well and a he ime of ; H (, r 0 ) is he iniial head a he poin wih he disance of r o he pumping well; s(,) r is he head drawdown a he poin wih he disance of r o he pumping well and a he ime of. Since in Equaion (2) he head is a coninuiy funcion of ime, i can no o be applied o se model boundary condiions ye; i need o be emporally discreized o n periods, and in each period he head is a consan, hus, he variaion of he head drawdown can be described by Equaion (3): s1, 0 1 s2, 1 2 sr (,) (3)... sn, n 1 n So, he live head anywhere wihin he cone of depression during pumping es can be given by he piecewise consan funcion as Equaion (4): Hr (, 0) s1, 0 1 Hr (, 0) s2, 1 2 Hr (,) S (4)... Hr (, 0) sn, n 1 n Divide he simulaion ime o n sress periods, and in each period se he model boundary head according o he corresponding consan value of each definiion domain 0 of he funcion shown in Equaion (4) Calibraions of he Model Boundary Condiions Afer iniial boundary condiions being se hrough he heoreical calculaion menioned above, run he model o calibrae boundary condiions and parameers ieraively using observaion daa of he pumping es and recovery es independenly. The flow char of he calibraion process is shown in Figure 1. Firs, build groundwaer flow model of he pumping es and se he model iniial boundary condiions for ieraion according o Equaion (4). Then, run he model buil in he former sep o calibrae model parameers and boundary condiions by he observaion daa of he pumping es. If he sandard error of esimae (S.E.E) exceeds 5%, he hydrogeology parameers and boundary condiions would be adjused slighly and hen he calibraion repeas. When he S.E.E is below 5%, he calibraion process goes o he nex sep. Finally, build groundwaer model of he recovery es applying he calibraion resuls of he second sep as iniial sae; run i o calibrae model parameers and boundary condiions again using he observaion daa of he recovery es. The S.E.E of 5% also is applied as calibraion error crierion; if resuls mee he crierion, he calibraion process ends and he model parameers and boundary condiions are fixed; oherwise, reurn o he firs sep, modify he model parameers and boundary condiions and hen he whole process repeas again. 3. Applicaion Example 3.1. Backgrounds and Model Overview The sudy was conduced a he piedmon alluvial slope in he souhern region of he Turpan-Hami basin. The Figure 1. Flow char of he model calibraion process.
3 52 Y. P. ZHOU ET AL. ineres confined groundwaer sysem has sable impervious roof and base; he groundwaer flew from souhwes o norheas wih a hydraulic gradien of 0.02, and was mainly recharged indirecly by he Quaernary phreaic waer from he souhern mounainous bedrock fissure. The sudied issue was abou simulaion of groundwaer and leaching soluion flow in he ore-bearing aquifer a an in-siu leaching of uranium sie. There are five wells (Figure 2); he well CK1 was for exracion, and he ress were for leaching soluion injecion. The hydrogeology characerisic of he aquifer wihin he mining scope is as shown in Figure 3. The average hickness of he aquifer is abou 40m; he sable impervious roof and base are mainly of mudsone and sily mudsone (in gray); in he aquifer (in blue) here are four disconinuous inerlayer, one is of sily mudsone (shown in gray) wih he hickness of 1-3 meers, and he hree ohers are of calcareous sandsone wih he hickness of m (in whie). Field pumping es and recovery es were conduced employing well CK1 for pumping, well ZK1 and well ZK3 for observaion, he es resuls are shown in Table 1. Figure 2. The plan view of well disribuion. Figure 3. Cross secion of he aquifer. Table 1 field ess and he hydraulic parameers of he aquifer. Field es Tes ime (min) Pumping rae (m 3 /h) Pumping es Recovery es Coefficien of ransmissibiliy (m 2 /d) coefficien of sorage Based on hose field invesigaions, he concepual model of he groundwaer flow during he pumping es was buil as Equaion (5): H H H ( K xx ) ( Kyy ) ( Kzz ) W x x y y z z H Ss,( x, y, z) D, 0 (5) H( x, y, z, o) H0 ( x, y, z) ( x, y, z) D H ( xyz,,, ) f( xyz,,, ) ( xyz,, ) S S where K xx, K yy and K zz are respecively he conduciviies in x, y and z direcion in he hree-dimension space; H is he confined waer head; W is he flux rae per uni volume, which is used o describe he flow rae of wells; Ss is he specific soraiviy; is he ime; H 0 (x,y,z) is he iniial head a posiion wih he coordinae (x, y, z); H ( xyz,,, ) S is he head a he model boundary. The model area denoed D and boundary S. The iniial head of H 0 (x, y, z) was se according o saic waer level observed a he beginning of he pumping es. Since he modeling area was small, is edges are far away from he naural hydraulic boundaries known, hence arificial boundaries were necessary. A circle surrounding he well CK1 and wih radius of 100 meers was se as he boundary of he model Head Funcions of he Model Boundary In case of he coefficien of ransmissibiliy is 5.7 m 2 /d, he coefficien of sorage is and he pumping rae is 7.2 m 3 /h, he radius of influence of he pumping es is more han 800 meers; clearly, he heads a he pumping es model boundaries which were 100 meers away from he pumping well mus o be varying wih ime. Calculaion according o Jacob formula showed ha head drawdown sared o ake place a he model boundaries afer pumping for 3.65 hours. The drawdown funcion can be derived from Equaion (1), as in 0, s (6) 2.34 ln 2.95, Temporally discreized he coninuiy funcion of he drawdown o a piecewise consan one, as shown in Figure 4. A broken line of discree funcion s () was used
4 Y. P. ZHOU ET AL. 53 o approximaely replace he curve of he original funcion s(), making sure he difference beween he wo adjacen consans of he funcion s () was no greaer han 0.5 m. The funcion values are lised in Table 2. And hen, divided he whole simulaion ime of he pumping es ino 15 periods according o he funcion s (), in each period iniial boundary heads of he model were se as a corresponding consan value, which could be ge from he Equaion (4) by replacing he values of funcion s () lised in Table 2 for he drawdown. Then, he ieraive calibraion process sared Resuls Afer correcing he model parameers and boundary condiions repeaedly via he calibraion processes of he pumping es model and recovery es model, he model parameers and boundary condiions were fixed on. Resuls showed ha he calculaed heads mached he observed daa saisfacorily in boh wo models (Figure 5 and Figure 6); he mean absolue error beween he calculaed heads and observed daa of pumping es simulaion was m, and recovery es simulaion m; boh variances were less han 5%. Furhermore, he coefficien of ransmissibiliy was 5.3 m 2 /d; i was close o he resuls of field ess (5.7 m 2 /d). Figure 5. The resuls of calculaed heads mached o observed heads of he pumping es model. Figure 6. The resuls of calculaed heads mached o observed heads of he recovery es model. Figure 4. The curves of funcion s () and funcion s'(). Table 2. The values of funcion s'(). s' ( ) s' ( ) s' ( ) Conclusions The boundary condiion is one of key facors influencing he reliabiliy of groundwaer model. In case of arificial boundary condiions are needed, hey should be se reasonably using as much field invesigaion daa as we can ge, oherwise, i is prone o cause grea disorion o he ruh and make he model worhless. Sudy resuls show ha he mehod inroduced in his paper can be a feasible choice o se arificial boundary condiions of he groundwaer model. 5. Acknowledgemens Thanks for he fund of he Major Sae Basic Research Developmen Program of China (973 Program) (No. 2012CB723101). REFERENCES [1] H. R. Zhang, Unseady Groundwaer Flow Theory Developmen and Applicaion, Geology Publishing House,
5 54 Y. P. ZHOU ET AL. Beijing, 1975, pp [2] H. R. Zhang, Groundwaer Hydraulics Developmen, Geology Publishing House, Beijing, 1992, pp [3] E. M. Labolle, A. A. Ahmed and G. E. Fogg, Review of he Inegraed Groundwaer and Surface Waer Model (IGSM), Ground waer, Vol. 41, No. 2, 2003, pp doi: /j b02587.x [4] L. H. Wei, L. C. Shu and Z. C. Hao, The Presen Siuaion and Developmen Tendency of Groundwaer Flow Numerical Simulaion, Journal of Chongqing Universiy, Vol. 23, 2000, pp [5] J. H. Ding, D. L. Zhou and S. Z. Ma, The Sae-of -Ar and Trends of Developmen of Groundwaer Modeling Sofware Abroad, Sie Invesigaion Science and Technology, No. 1, 2002, pp [6] L. C. Xu, Inroducion o Common Sofware Producs Modeling Groundwaer, Uranium Mining and Meallurgy, Vol. 21, No. 1, 2002, pp [7] H. Q. Qu, S. Zeng and H. L. Liu, Research Saus and Developmen of Numerical Simulaion of In-siu Leaching of Uranium, Sci-ech Informaion Developmen & Economy, Vol. 21, No. 9, 2011, pp [8] Y. Q. Xue, Presen Siuaion and Prospec of Groundwaer Numerical Simulaion in China, Geological Journal of China Universiies, Vol. 16, No.1, 2010, pp [9] W.X. Lu, Approach on Boundary Condiion in Numerical Simulaion of Groundwaer Flows, Journal of Hydraulic engineering, No. 3, 2003, pp [10] Y. Y. Shen and Y. Z. Jiang, Research on Disposal Mehod of Arificial Boundary Condiion in Numerical Simulaion of Groundwaer Flow, Hydrogeology and Engineering Geology, No. 6, 2008, pp [11] De lange W. J. and A. Cauchy, Boundary Condiion for he Lumped Ineracion Beween an Arbirary Number of Surface Waers and a Regional Aquifer, Journal of Hydrology, 1999, pp [12] Y. Q. Xue, Groundwaer Dynamics, Geology Publishing House, Beijing, 1997, pp
Simulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010
Simulaion-Solving Dynamic Models ABE 5646 Week 2, Spring 2010 Week Descripion Reading Maerial 2 Compuer Simulaion of Dynamic Models Finie Difference, coninuous saes, discree ime Simple Mehods Euler Trapezoid
More informationVehicle 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 informationModule 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 informationPhysics 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 informationKriging Models Predicting Atrazine Concentrations in Surface Water Draining Agricultural Watersheds
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 Kriging Models Predicing Arazine Concenraions in Surface Waer Draining Agriculural Waersheds Paul L. Mosquin, Jeremy Aldworh, Wenlin Chen Supplemenal Maerial Number
More informationExponential Weighted Moving Average (EWMA) Chart Under The Assumption of Moderateness And Its 3 Control Limits
DOI: 0.545/mjis.07.5009 Exponenial Weighed Moving Average (EWMA) Char Under The Assumpion of Moderaeness And Is 3 Conrol Limis KALPESH S TAILOR Assisan Professor, Deparmen of Saisics, M. K. Bhavnagar Universiy,
More informationSome Basic Information about M-S-D Systems
Some Basic Informaion abou M-S-D Sysems 1 Inroducion We wan o give some summary of he facs concerning unforced (homogeneous) and forced (non-homogeneous) models for linear oscillaors governed by second-order,
More informationEcon107 Applied Econometrics Topic 7: Multicollinearity (Studenmund, Chapter 8)
I. Definiions and Problems A. Perfec Mulicollineariy Econ7 Applied Economerics Topic 7: Mulicollineariy (Sudenmund, Chaper 8) Definiion: Perfec mulicollineariy exiss in a following K-variable regression
More informationAir Traffic Forecast Empirical Research Based on the MCMC Method
Compuer and Informaion Science; Vol. 5, No. 5; 0 ISSN 93-8989 E-ISSN 93-8997 Published by Canadian Cener of Science and Educaion Air Traffic Forecas Empirical Research Based on he MCMC Mehod Jian-bo Wang,
More informationNavneet 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 information2.7. Some common engineering functions. Introduction. Prerequisites. Learning Outcomes
Some common engineering funcions 2.7 Inroducion This secion provides a caalogue of some common funcions ofen used in Science and Engineering. These include polynomials, raional funcions, he modulus funcion
More informationDevelopment of a new metrological model for measuring of the water surface evaporation Tovmach L. Tovmach Yr. Abstract Introduction
Developmen of a new merological model for measuring of he waer surface evaporaion Tovmach L. Tovmach Yr. Sae Hydrological Insiue 23 Second Line, 199053 S. Peersburg, Russian Federaion Telephone (812) 323
More informationA Dynamic Model of Economic Fluctuations
CHAPTER 15 A Dynamic Model of Economic Flucuaions Modified for ECON 2204 by Bob Murphy 2016 Worh Publishers, all righs reserved IN THIS CHAPTER, OU WILL LEARN: how o incorporae dynamics ino he AD-AS model
More information04. Kinetics of a second order reaction
4. Kineics of a second order reacion Imporan conceps Reacion rae, reacion exen, reacion rae equaion, order of a reacion, firs-order reacions, second-order reacions, differenial and inegraed rae laws, Arrhenius
More informationLecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still.
Lecure - Kinemaics in One Dimension Displacemen, Velociy and Acceleraion Everyhing in he world is moving. Nohing says sill. Moion occurs a all scales of he universe, saring from he moion of elecrons in
More informationSTATE-SPACE MODELLING. A mass balance across the tank gives:
B. Lennox and N.F. Thornhill, 9, Sae Space Modelling, IChemE Process Managemen and Conrol Subjec Group Newsleer STE-SPACE MODELLING Inroducion: Over he pas decade or so here has been an ever increasing
More informationKINEMATICS IN ONE DIMENSION
KINEMATICS IN ONE DIMENSION PREVIEW Kinemaics is he sudy of how hings move how far (disance and displacemen), how fas (speed and velociy), and how fas ha how fas changes (acceleraion). We say ha an objec
More informationdt = C exp (3 ln t 4 ). t 4 W = C exp ( ln(4 t) 3) = C(4 t) 3.
Mah Rahman Exam Review Soluions () Consider he IVP: ( 4)y 3y + 4y = ; y(3) = 0, y (3) =. (a) Please deermine he longes inerval for which he IVP is guaraneed o have a unique soluion. Soluion: The disconinuiies
More informationDiebold, 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 information3.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 informationSPH3U: Projectiles. Recorder: Manager: Speaker:
SPH3U: Projeciles Now i s ime o use our new skills o analyze he moion of a golf ball ha was ossed hrough he air. Le s find ou wha is special abou he moion of a projecile. Recorder: Manager: Speaker: 0
More informationOrdinary Differential Equations
Lecure 22 Ordinary Differenial Equaions Course Coordinaor: Dr. Suresh A. Karha, Associae Professor, Deparmen of Civil Engineering, IIT Guwahai. In naure, mos of he phenomena ha can be mahemaically described
More information15. Vector Valued Functions
1. Vecor Valued Funcions Up o his poin, we have presened vecors wih consan componens, for example, 1, and,,4. However, we can allow he componens of a vecor o be funcions of a common variable. For example,
More informationTime series Decomposition method
Time series Decomposiion mehod A ime series is described using a mulifacor model such as = f (rend, cyclical, seasonal, error) = f (T, C, S, e) Long- Iner-mediaed Seasonal Irregular erm erm effec, effec,
More informationClass Meeting # 10: Introduction to the Wave Equation
MATH 8.5 COURSE NOTES - CLASS MEETING # 0 8.5 Inroducion o PDEs, Fall 0 Professor: Jared Speck Class Meeing # 0: Inroducion o he Wave Equaion. Wha is he wave equaion? The sandard wave equaion for a funcion
More informationReliability of Technical Systems
eliabiliy of Technical Sysems Main Topics Inroducion, Key erms, framing he problem eliabiliy parameers: Failure ae, Failure Probabiliy, Availabiliy, ec. Some imporan reliabiliy disribuions Componen reliabiliy
More informationStability and Bifurcation in a Neural Network Model with Two Delays
Inernaional Mahemaical Forum, Vol. 6, 11, no. 35, 175-1731 Sabiliy and Bifurcaion in a Neural Nework Model wih Two Delays GuangPing Hu and XiaoLing Li School of Mahemaics and Physics, Nanjing Universiy
More informationFinal 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 informationPET467E-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 information1 Differential Equation Investigations using Customizable
Differenial Equaion Invesigaions using Cusomizable Mahles Rober Decker The Universiy of Harford Absrac. The auhor has developed some plaform independen, freely available, ineracive programs (mahles) for
More informationHamilton- 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 informationIn this chapter the model of free motion under gravity is extended to objects projected at an angle. When you have completed it, you should
Cambridge Universiy Press 978--36-60033-7 Cambridge Inernaional AS and A Level Mahemaics: Mechanics Coursebook Excerp More Informaion Chaper The moion of projeciles In his chaper he model of free moion
More informationd 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 information6.003 Homework 1. Problems. Due at the beginning of recitation on Wednesday, February 10, 2010.
6.003 Homework Due a he beginning of reciaion on Wednesday, February 0, 200. Problems. Independen and Dependen Variables Assume ha he heigh of a waer wave is given by g(x v) where x is disance, v is velociy,
More informationChapter 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 informationRandom Walk with Anti-Correlated Steps
Random Walk wih Ani-Correlaed Seps John Noga Dirk Wagner 2 Absrac We conjecure he expeced value of random walks wih ani-correlaed seps o be exacly. We suppor his conjecure wih 2 plausibiliy argumens and
More information20. 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 informationZürich. ETH Master Course: L Autonomous Mobile Robots Localization II
Roland Siegwar Margaria Chli Paul Furgale Marco Huer Marin Rufli Davide Scaramuzza ETH Maser Course: 151-0854-00L Auonomous Mobile Robos Localizaion II ACT and SEE For all do, (predicion updae / ACT),
More informationImproved Approximate Solutions for Nonlinear Evolutions Equations in Mathematical Physics Using the Reduced Differential Transform Method
Journal of Applied Mahemaics & Bioinformaics, vol., no., 01, 1-14 ISSN: 179-660 (prin), 179-699 (online) Scienpress Ld, 01 Improved Approimae Soluions for Nonlinear Evoluions Equaions in Mahemaical Physics
More informationWednesday, November 7 Handout: Heteroskedasticity
Amhers College Deparmen of Economics Economics 360 Fall 202 Wednesday, November 7 Handou: Heeroskedasiciy Preview Review o Regression Model o Sandard Ordinary Leas Squares (OLS) Premises o Esimaion Procedures
More informationExplaining Total Factor Productivity. Ulrich Kohli University of Geneva December 2015
Explaining Toal Facor Produciviy Ulrich Kohli Universiy of Geneva December 2015 Needed: A Theory of Toal Facor Produciviy Edward C. Presco (1998) 2 1. Inroducion Toal Facor Produciviy (TFP) has become
More information2. Nonlinear Conservation Law Equations
. Nonlinear Conservaion Law Equaions One of he clear lessons learned over recen years in sudying nonlinear parial differenial equaions is ha i is generally no wise o ry o aack a general class of nonlinear
More informationRobust estimation based on the first- and third-moment restrictions of the power transformation model
h Inernaional Congress on Modelling and Simulaion, Adelaide, Ausralia, 6 December 3 www.mssanz.org.au/modsim3 Robus esimaion based on he firs- and hird-momen resricions of he power ransformaion Nawaa,
More informationLecture 4 Kinetics of a particle Part 3: Impulse and Momentum
MEE Engineering Mechanics II Lecure 4 Lecure 4 Kineics of a paricle Par 3: Impulse and Momenum Linear impulse and momenum Saring from he equaion of moion for a paricle of mass m which is subjeced o an
More informationDEPARTMENT OF STATISTICS
A Tes for Mulivariae ARCH Effecs R. Sco Hacker and Abdulnasser Haemi-J 004: DEPARTMENT OF STATISTICS S-0 07 LUND SWEDEN A Tes for Mulivariae ARCH Effecs R. Sco Hacker Jönköping Inernaional Business School
More informationOn Measuring Pro-Poor Growth. 1. On Various Ways of Measuring Pro-Poor Growth: A Short Review of the Literature
On Measuring Pro-Poor Growh 1. On Various Ways of Measuring Pro-Poor Growh: A Shor eview of he Lieraure During he pas en years or so here have been various suggesions concerning he way one should check
More informationA First Course on Kinetics and Reaction Engineering. Class 19 on Unit 18
A Firs ourse on Kineics and Reacion Engineering lass 19 on Uni 18 Par I - hemical Reacions Par II - hemical Reacion Kineics Where We re Going Par III - hemical Reacion Engineering A. Ideal Reacors B. Perfecly
More informationA Specification Test for Linear Dynamic Stochastic General Equilibrium Models
Journal of Saisical and Economeric Mehods, vol.1, no.2, 2012, 65-70 ISSN: 2241-0384 (prin), 2241-0376 (online) Scienpress Ld, 2012 A Specificaion Tes for Linear Dynamic Sochasic General Equilibrium Models
More informationTraveling Waves. Chapter Introduction
Chaper 4 Traveling Waves 4.1 Inroducion To dae, we have considered oscillaions, i.e., periodic, ofen harmonic, variaions of a physical characerisic of a sysem. The sysem a one ime is indisinguishable from
More informationRC, RL and RLC circuits
Name Dae Time o Complee h m Parner Course/ Secion / Grade RC, RL and RLC circuis Inroducion In his experimen we will invesigae he behavior of circuis conaining combinaions of resisors, capaciors, and inducors.
More informationLet us start with a two dimensional case. We consider a vector ( x,
Roaion marices We consider now roaion marices in wo and hree dimensions. We sar wih wo dimensions since wo dimensions are easier han hree o undersand, and one dimension is a lile oo simple. However, our
More informationVariational Iteration Method for Solving System of Fractional Order Ordinary Differential Equations
IOSR Journal of Mahemaics (IOSR-JM) e-issn: 2278-5728, p-issn: 2319-765X. Volume 1, Issue 6 Ver. II (Nov - Dec. 214), PP 48-54 Variaional Ieraion Mehod for Solving Sysem of Fracional Order Ordinary Differenial
More informationMethodology. -ratios are biased and that the appropriate critical values have to be increased by an amount. that depends on the sample size.
Mehodology. Uni Roo Tess A ime series is inegraed when i has a mean revering propery and a finie variance. I is only emporarily ou of equilibrium and is called saionary in I(0). However a ime series ha
More informationLAPLACE 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 information2017 3rd International Conference on E-commerce and Contemporary Economic Development (ECED 2017) ISBN:
7 3rd Inernaional Conference on E-commerce and Conemporary Economic Developmen (ECED 7) ISBN: 978--6595-446- Fuures Arbirage of Differen Varieies and based on he Coinegraion Which is under he Framework
More informationWEEK-3 Recitation PHYS 131. of the projectile s velocity remains constant throughout the motion, since the acceleration a x
WEEK-3 Reciaion PHYS 131 Ch. 3: FOC 1, 3, 4, 6, 14. Problems 9, 37, 41 & 71 and Ch. 4: FOC 1, 3, 5, 8. Problems 3, 5 & 16. Feb 8, 018 Ch. 3: FOC 1, 3, 4, 6, 14. 1. (a) The horizonal componen of he projecile
More informationProblem Set 5. Graduate Macro II, Spring 2017 The University of Notre Dame Professor Sims
Problem Se 5 Graduae Macro II, Spring 2017 The Universiy of Nore Dame Professor Sims Insrucions: You may consul wih oher members of he class, bu please make sure o urn in your own work. Where applicable,
More informationACE 562 Fall Lecture 5: The Simple Linear Regression Model: Sampling Properties of the Least Squares Estimators. by Professor Scott H.
ACE 56 Fall 005 Lecure 5: he Simple Linear Regression Model: Sampling Properies of he Leas Squares Esimaors by Professor Sco H. Irwin Required Reading: Griffihs, Hill and Judge. "Inference in he Simple
More informationMathcad Lecture #8 In-class Worksheet Curve Fitting and Interpolation
Mahcad Lecure #8 In-class Workshee Curve Fiing and Inerpolaion A he end of his lecure, you will be able o: explain he difference beween curve fiing and inerpolaion decide wheher curve fiing or inerpolaion
More informationMath 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 informationChapter 15: Phenomena. Chapter 15 Chemical Kinetics. Reaction Rates. Reaction Rates R P. Reaction Rates. Rate Laws
Chaper 5: Phenomena Phenomena: The reacion (aq) + B(aq) C(aq) was sudied a wo differen emperaures (98 K and 35 K). For each emperaure he reacion was sared by puing differen concenraions of he 3 species
More informationHW6: MRI Imaging Pulse Sequences (7 Problems for 100 pts)
HW6: MRI Imaging Pulse Sequences (7 Problems for 100 ps) GOAL The overall goal of HW6 is o beer undersand pulse sequences for MRI image reconsrucion. OBJECTIVES 1) Design a spin echo pulse sequence o image
More informationKeywords: thermal stress; thermal fatigue; inverse analysis; heat conduction; regularization
Proceedings Inverse Analysis for Esimaing Temperaure and Residual Sress Disribuions in a Pipe from Ouer Surface Temperaure Measuremen and Is Regularizaion Shiro Kubo * and Shoki Taguwa Deparmen of Mechanical
More informationNumerical Dispersion
eview of Linear Numerical Sabiliy Numerical Dispersion n he previous lecure, we considered he linear numerical sabiliy of boh advecion and diffusion erms when approimaed wih several spaial and emporal
More informationLab #2: Kinematics in 1-Dimension
Reading Assignmen: Chaper 2, Secions 2-1 hrough 2-8 Lab #2: Kinemaics in 1-Dimension Inroducion: The sudy of moion is broken ino wo main areas of sudy kinemaics and dynamics. Kinemaics is he descripion
More informationSection 7.4 Modeling Changing Amplitude and Midline
488 Chaper 7 Secion 7.4 Modeling Changing Ampliude and Midline While sinusoidal funcions can model a variey of behaviors, i is ofen necessary o combine sinusoidal funcions wih linear and exponenial curves
More informationNature 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 informationSterilization 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 informationThe 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 information23.5. Half-Range Series. Introduction. Prerequisites. Learning Outcomes
Half-Range Series 2.5 Inroducion In his Secion we address he following problem: Can we find a Fourier series expansion of a funcion defined over a finie inerval? Of course we recognise ha such a funcion
More informationMath 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 informationA note on spurious regressions between stationary series
A noe on spurious regressions beween saionary series Auhor Su, Jen-Je Published 008 Journal Tile Applied Economics Leers DOI hps://doi.org/10.1080/13504850601018106 Copyrigh Saemen 008 Rouledge. This is
More informationt is a basis for the solution space to this system, then the matrix having these solutions as columns, t x 1 t, x 2 t,... x n t x 2 t...
Mah 228- Fri Mar 24 5.6 Marix exponenials and linear sysems: The analogy beween firs order sysems of linear differenial equaions (Chaper 5) and scalar linear differenial equaions (Chaper ) is much sronger
More informationLecture 4 Notes (Little s Theorem)
Lecure 4 Noes (Lile s Theorem) This lecure concerns one of he mos imporan (and simples) heorems in Queuing Theory, Lile s Theorem. More informaion can be found in he course book, Bersekas & Gallagher,
More information5.2. The Natural Logarithm. Solution
5.2 The Naural Logarihm The number e is an irraional number, similar in naure o π. Is non-erminaing, non-repeaing value is e 2.718 281 828 59. Like π, e also occurs frequenly in naural phenomena. In fac,
More informationCalculation of the Two High Voltage Transmission Line Conductors Minimum Distance
World Journal of Engineering and Technology, 15, 3, 89-96 Published Online Ocober 15 in SciRes. hp://www.scirp.org/journal/wje hp://dx.doi.org/1.436/wje.15.33c14 Calculaion of he Two High Volage Transmission
More informationSub Module 2.6. Measurement of transient temperature
Mechanical Measuremens Prof. S.P.Venkaeshan Sub Module 2.6 Measuremen of ransien emperaure Many processes of engineering relevance involve variaions wih respec o ime. The sysem properies like emperaure,
More informationSolutions 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 informationComparing Means: t-tests for One Sample & Two Related Samples
Comparing Means: -Tess for One Sample & Two Relaed Samples Using he z-tes: Assumpions -Tess for One Sample & Two Relaed Samples The z-es (of a sample mean agains a populaion mean) is based on he assumpion
More informationGMM - Generalized Method of Moments
GMM - Generalized Mehod of Momens Conens GMM esimaion, shor inroducion 2 GMM inuiion: Maching momens 2 3 General overview of GMM esimaion. 3 3. Weighing marix...........................................
More informationEE650R: Reliability Physics of Nanoelectronic Devices Lecture 9:
EE65R: Reliabiliy Physics of anoelecronic Devices Lecure 9: Feaures of Time-Dependen BTI Degradaion Dae: Sep. 9, 6 Classnoe Lufe Siddique Review Animesh Daa 9. Background/Review: BTI is observed when he
More informationINVERSE RESPONSE COMPENSATION BY ESTIMATING PARAMETERS OF A PROCESS COMPRISING OF TWO FIRST ORDER SYSTEMS
Inernaional Journal of Informaion Technology and nowledge Managemen July-December 0, Volume 5, No., pp. 433-438 INVERSE RESPONSE COMPENSATION BY ESTIMATING PARAMETERS OF A PROCESS COMPRISING OF TWO FIRST
More informationAn Invariance for (2+1)-Extension of Burgers Equation and Formulae to Obtain Solutions of KP Equation
Commun Theor Phys Beijing, China 43 2005 pp 591 596 c Inernaional Academic Publishers Vol 43, No 4, April 15, 2005 An Invariance for 2+1-Eension of Burgers Equaion Formulae o Obain Soluions of KP Equaion
More informationMath 333 Problem Set #2 Solution 14 February 2003
Mah 333 Problem Se #2 Soluion 14 February 2003 A1. Solve he iniial value problem dy dx = x2 + e 3x ; 2y 4 y(0) = 1. Soluion: This is separable; we wrie 2y 4 dy = x 2 + e x dx and inegrae o ge The iniial
More informationIntroduction D P. r = constant discount rate, g = Gordon Model (1962): constant dividend growth rate.
Inroducion Gordon Model (1962): D P = r g r = consan discoun rae, g = consan dividend growh rae. If raional expecaions of fuure discoun raes and dividend growh vary over ime, so should he D/P raio. Since
More informationIMPROVED HYPERBOLIC FUNCTION METHOD AND EXACT SOLUTIONS FOR VARIABLE COEFFICIENT BENJAMIN-BONA-MAHONY-BURGERS EQUATION
THERMAL SCIENCE, Year 015, Vol. 19, No. 4, pp. 1183-1187 1183 IMPROVED HYPERBOLIC FUNCTION METHOD AND EXACT SOLUTIONS FOR VARIABLE COEFFICIENT BENJAMIN-BONA-MAHONY-BURGERS EQUATION by Hong-Cai MA a,b*,
More informationSTRUCTURAL CHANGE IN TIME SERIES OF THE EXCHANGE RATES BETWEEN YEN-DOLLAR AND YEN-EURO IN
Inernaional Journal of Applied Economerics and Quaniaive Sudies. Vol.1-3(004) STRUCTURAL CHANGE IN TIME SERIES OF THE EXCHANGE RATES BETWEEN YEN-DOLLAR AND YEN-EURO IN 001-004 OBARA, Takashi * Absrac The
More informationMulti-scale 2D acoustic full waveform inversion with high frequency impulsive source
Muli-scale D acousic full waveform inversion wih high frequency impulsive source Vladimir N Zubov*, Universiy of Calgary, Calgary AB vzubov@ucalgaryca and Michael P Lamoureux, Universiy of Calgary, Calgary
More informationBias in Conditional and Unconditional Fixed Effects Logit Estimation: a Correction * Tom Coupé
Bias in Condiional and Uncondiional Fixed Effecs Logi Esimaion: a Correcion * Tom Coupé Economics Educaion and Research Consorium, Naional Universiy of Kyiv Mohyla Academy Address: Vul Voloska 10, 04070
More informationLicenciatura de ADE y Licenciatura conjunta Derecho y ADE. Hoja de ejercicios 2 PARTE A
Licenciaura de ADE y Licenciaura conjuna Derecho y ADE Hoja de ejercicios PARTE A 1. Consider he following models Δy = 0.8 + ε (1 + 0.8L) Δ 1 y = ε where ε and ε are independen whie noise processes. In
More information) were both constant and we brought them from under the integral.
YIELD-PER-RECRUIT (coninued The yield-per-recrui model applies o a cohor, bu we saw in he Age Disribuions lecure ha he properies of a cohor do no apply in general o a collecion of cohors, which is wha
More informationThe fundamental mass balance equation is ( 1 ) where: I = inputs P = production O = outputs L = losses A = accumulation
Hea (iffusion) Equaion erivaion of iffusion Equaion The fundamenal mass balance equaion is I P O L A ( 1 ) where: I inpus P producion O oupus L losses A accumulaion Assume ha no chemical is produced or
More information1. VELOCITY AND ACCELERATION
1. VELOCITY AND ACCELERATION 1.1 Kinemaics Equaions s = u + 1 a and s = v 1 a s = 1 (u + v) v = u + as 1. Displacemen-Time Graph Gradien = speed 1.3 Velociy-Time Graph Gradien = acceleraion Area under
More information18 Biological models with discrete time
8 Biological models wih discree ime The mos imporan applicaions, however, may be pedagogical. The elegan body of mahemaical heory peraining o linear sysems (Fourier analysis, orhogonal funcions, and so
More informationTesting for a Single Factor Model in the Multivariate State Space Framework
esing for a Single Facor Model in he Mulivariae Sae Space Framework Chen C.-Y. M. Chiba and M. Kobayashi Inernaional Graduae School of Social Sciences Yokohama Naional Universiy Japan Faculy of Economics
More informationNotes on Kalman Filtering
Noes on Kalman Filering Brian Borchers and Rick Aser November 7, Inroducion Daa Assimilaion is he problem of merging model predicions wih acual measuremens of a sysem o produce an opimal esimae of he curren
More informationSolutions 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 informationACE 564 Spring Lecture 7. Extensions of The Multiple Regression Model: Dummy Independent Variables. by Professor Scott H.
ACE 564 Spring 2006 Lecure 7 Exensions of The Muliple Regression Model: Dumm Independen Variables b Professor Sco H. Irwin Readings: Griffihs, Hill and Judge. "Dumm Variables and Varing Coefficien Models
More informationBifurcation Analysis of a Stage-Structured Prey-Predator System with Discrete and Continuous Delays
Applied Mahemaics 4 59-64 hp://dx.doi.org/.46/am..4744 Published Online July (hp://www.scirp.org/ournal/am) Bifurcaion Analysis of a Sage-Srucured Prey-Predaor Sysem wih Discree and Coninuous Delays Shunyi
More information6.2 Transforms of Derivatives and Integrals.
SEC. 6.2 Transforms of Derivaives and Inegrals. ODEs 2 3 33 39 23. Change of scale. If l( f ()) F(s) and c is any 33 45 APPLICATION OF s-shifting posiive consan, show ha l( f (c)) F(s>c)>c (Hin: In Probs.
More information