Introduction to groundwater flow modeling: finite difference methods. Tyson Strand
|
|
- Emory Martin
- 6 years ago
- Views:
Transcription
1 Introducton to groundwater flow odelng: fnte dfference etods Tson Strand
2 Darc s law contnut and te groundwater flow equaton Fundaentals of fnte dfference etods 3 FD soluton of Laplace s equaton 4 FD soluton of Posson s equaton 5 Transent flow
3 Darc s law contnut and te groundwater flow equaton q d q q d q
4 d q q d q q d d 0 Evertng tat goes n ust coe out q d d q d q d q 0 dd 0 After dvdng b te volue area let d and d go to zero q q 0 Contnut
5 Darc s law: sotropc edu K q K q z K q z Scalar for K k z K ˆ ˆ ˆ Vector for
6 Darc s law contnut n D K q K q 0 q q 0 K K K 0 Laplace s equaton
7 Te groundwater flow equaton t S W z K z K K s zz Sources and snks Transent flow ter Darc s law contnut
8 Splfng assuptons Isotropc edu: K K K K zz K W t W ft Stead-state flow Posson s equaton W f Stead-state flow no sources/snks Laplace s equaton S W K s 0
9 Fundaentals of fnte dfference etods Dscretzaton of space Dscretzaton of contnuous quanttes Dscretzaton of te Te frst spatal dervatve Te second spatal dervatve Boundar condtons and ntal condtons Solvng te proble
10 Dsctretzaton of space L N d - d - d
11 Dscretzaton of quanttes Eac spatal locaton as assocated wt t a ead value and a vector darc flu - q q -
12 Dscretzaton of te t n nt nt n3t n4t n5t t t t t 3 t 4 t 5 Increasng te
13 Te frst spatal dervatve - - Upstrea dfference Central dfference
14 Te second spatal dervatve: step - NOTE: Ts s essentall a central dfference approaton of te frst dervatve evaluated at te d-postons -
15 Te second spatal dervatve: step Te second dervatve s ust te dervatve of te dervatve
16 Boundar condtons Wat are boundar condtons and w are te necessar? Consder te sple proble: 0. We can solve ts drectl b separatng varables and ntegratng: C
17 0. C Tere are an nfnte nuber of solutons to ts equaton. Te are referred to as a fal of curves We ust specf a value of at a known poston n order to solve for C.
18 For a second spatal dervatve two boundar condtons ust be specfed. Tere are 3 tpes of bc s tat we can appl Head s specfed at a boundar - Called Drclet condtons Flow frst dervatve of ead s specfed at a boundar - Called Neuann condtons 3 Soe cobnaton of and - Called ed condtons
19 Intal condtons Equvalent to boundar condtons ecept tat te boundar s now teporal nstead of spatal Tpcall te state of te sste.e. ead values are specfed at te t 0 K K z K zz z W S s t Queston: How an boundar and ntal condtons Would be needed to solve te above governng equaton?
20 Solvng te proble Set of dfferental equatons Mateatcal odel Fnte Dfference Fnte Eleent Set of algebrac equatons dscrete odel Calculus tecnques Iteratve or drect etods Analtcal soluton Not alwas possble Copare If possble Approate soluton Feld observatons After Wang & Anderson 98
21 Recall ow we arrved at Laplace s equaton K q K q 0 q q 0 K K K 0 Laplace s equaton We are focusng on teratve etods
22 Solvng Laplace s equaton D 0 0
23 Let [ ] 4
24 [ ] 4 Set equal to zero 0 4 Te above equaton s te basc fnte dfference soluton to Laplace s equaton
25 Now we rearrange te prevous equaton so tat we can pleent t nto our regular grd Solve for 4
26 Now wat do we do? Specf boundar condtons Guess at ntal values for ead at all locatons Begn teratons
27 Wat does t ean to terate? Gven: Do: Wle: We know te boundar condtons and ave ntal guess values for eadposton Update ead values at ever locaton based on prevous equaton ncludng boundares f applcable Convergence crtera are not et
28 Convergence Copare obtaned ead values at a gven teraton to tose at te prevous teraton Contnue to terate as long as ead values contnue to cange wtn soe preset lt 3 Test te soluton to be descrbed
29 Testng a coputatonal odel Coparson to oter analtcal results/feld data valdaton? Cangng te grd spacng Cangng te convergence crteron Mass balance ceck
30 Valdaton/verfcaton Analtcal soluton s known uncoon Coparson to feld data 3 Test on sple stuatons for wc or are known 4 Ceck odel predctons
31 Possble odel probles Stablt Robustness
32 Saple proble: Laplace s equaton
33 Lecture Posson s equaton Dgresson: Inflow outflow and sgn conventons Fnte dfference for for Posson s equaton Eaple progras solvng Posson s equaton Transent flow Dgresson: Storage paraeters Fnte dfference for for transent gw flow equaton eplct etods & stablt Eaple transent flow progra Iplct teratve etods Eaple transent flow progra full plct
34 Solvng Posson s equaton T R K W Wat s W? K W b R W Were b s te tckness of te aqufer z denson
35 Inflow outflow and sgn conventons Tnk about te agntudes of q out and q n Wat does t ean n ters of sources/snks? 0 > 0
36 Sgn conventons < 0 < 0 Tere ust be recarge or soe oter source at te REV stead-state assupton > 0 Tere ust be dscarge or soe oter snk at te REV stead-state assupton
37 W K R T You need to be careful as to ow tngs are defned W MUST be defned dscarge snk postve It s ore ntutve to defne a recarge source functon R tat s postve for recarge sources and negatve for dscarge snks Ts wll also elp gve us an ntutve understandng of te transent gw flow equaton
38 Back to Posson s equaton T R Proceedng as we dd for Laplace s equaton [ ] T R 4 Solvng for T R 4
39 Solvng Posson s equaton nuercall Bascall we can proceed eactl as we dd for Laplace s equaton usng te prevous fnte dfference approaton for Defne boundar condtons Set ntal guess values Iterate Ceck results
40 Saple proble: Posson s equaton wt unfor recarge/dscarge
41 Saple proble: Posson s equaton wt pont source recarge/dscarge
42 Solvng transent flow probles b t R t S K s We ve learned ow to take care of evertng ecept te te dervatve t t t Forward dfference approaton t t t Backward dfference approaton
43 Central dfference approaton t t t BAD NEWS Te central dfference approaton s uncondtonall unstable for te dervatve Ts etod sould NOT be used to estate te frst te dervatve
44 Storage paraeters S storage coeffcent or storatvt Densonless Measures te volue of water epelled absorbed per unt surface area per unt ead cange S s specfc storage S/b Unts [/L] also called elastc storage coeff. Measures te water volue per unt aqufer volue tat s epelled stored due to copressblt atr & water per unt ead cange Used for confned unts or te saturated parts of an unconfned unt b aqufer tckness for confned unts use b s tckness of te saturated regon for unconfned sstes
45 S specfc eld Densonless Specfc eld Used to descrbe unconfned sstes Rato of water volue tat drans fro a saturated regon of aqufer due to gravt forces per unt aqufer volue Generall several orders of agntude larger tan b s S s ecept n ver fne graned sstes Relatonsp between S S s and S S S b s S s Fetter 994
46 Storage paraeters concl. In unconfned sstes te specfc storage can generall be neglected Ecept n ver fne graned eda In confned sstes onl te specfc storage s pertnent so long as ead values rean above te upper confnng unt In general for coputatonal odels t s best to use te storage coeffcent storatvt Defntons of storage paraeters are not alwas consstent n lterature be careful
47 Transent flow So let s consder te transent flow equaton n te for S R t T t T Assue te left and sze s negatve out > n: ten b te conventon we ve set eter recarge s postve or te te rate of cange of storage s negatve rgt?
48 T t R t T S Let s wrte te fnte dfference approaton for ts equaton T R t T S Ts s called te eplct representaton snce all spatal ters are evaluated at te old te
49 Let and solve for te ead at at te nde [ ] S t T S R t S t T 4 Te eplct equaton s eas to solve: Gven a set of ntal condtons & bc s Te ead at all locatons at te nde can be calculated drectl fro te ead values at te prevous te nde Notce tere s no teraton once ou ave an ntal condton t s a sple atter to step forward n te b an aount t
50 Stablt of te eplct etod One ust be careful to keep te te steps sall enoug suc tat te eplct etod reans stable If te te step s too large fluctuatons wll develop n te ead as a functon of te wc wll be aplfed as te odel contnues to step troug te Te stablt crtera are: Tt < 0.5 S Tt < 0.5 S -D -D
51 Eaple transent proble: Eplct soluton D Intal Fnal We are nterested n wat appens n between te ntal and fnal states Te ead at te rgt boundar drops fro blue to red at te t 0 Te red potentoetrc surface s te fnal state
52 Iplct etods Notce n te eplct forulaton tat we ave evaluated all spatal dervatves at te old te t Is tat te best? Wat f we evaluated spatal dervatves at te new te t? Wat f we evaluated spatal dervatves alfwa between?
53 α α Were α can var between 0 full eplct and full plct A slar epresson can be wrtten for te second dervatve wt respect to We splf te algebra b defnng te followng and lettng 4 ~
54 [ ] [ ] ~ 4 ~ 4 α α Put ts approaton nto te gw flow equaton [ ] [ ] T R t T S ~ 4 ~ 4 α α Rearrange [ ] T R t T S t T S 4 ~ 4 ~ 4 α α α
55 Solvng te prevous equaton for te new ead value at [ ] T R t T S t T S 4 ~ 4 ~ 4 α α α Ts s te fundaental plct fnte dfference approaton α 0 Full eplct α ½ Crank-Ncolson α Full plct Te above epresson can be solved b teratve or drect etods We wll contnue to focus on teratve etods Te advantage of plct etods s STABILITY
56 To solve teratvel: - Sweep troug te lattce updatng te ead value at te new te b te above forula - Contnue to sweep troug te lattce updatng as long as te new ead values keep cangng - Stop wen new ead values do not cange wtn preset lts - update as ou go along: called Gauss-Sedel teraton - Repeat for te net te step
57 Full plct D saple proble t T S t T S ~ ~ Were s now defned b ~
58 References H. F. Wang and M. P. Anderson 98. Introducton to Groundwater Modelng Acadec Press Inc. C. W. Fetter 994. Appled Hdrogeolog Macllan College Publsng Copan. F. W. Scwartz and H. Zang 003. Fundaentals of Ground Water Wle publsng.
On Pfaff s solution of the Pfaff problem
Zur Pfaff scen Lösung des Pfaff scen Probles Mat. Ann. 7 (880) 53-530. On Pfaff s soluton of te Pfaff proble By A. MAYER n Lepzg Translated by D. H. Delpenc Te way tat Pfaff adopted for te ntegraton of
More informationOutline. Review Numerical Approach. Schedule for April and May. Review Simple Methods. Review Notation and Order
Sstes of Ordnar Dfferental Equatons Aprl, Solvng Sstes of Ordnar Dfferental Equatons Larr Caretto Mecancal Engneerng 9 Nuercal Analss of Engneerng Sstes Aprl, Outlne Revew bascs of nuercal solutons of
More informationRectilinear motion. Lecture 2: Kinematics of Particles. External motion is known, find force. External forces are known, find motion
Lecture : Kneatcs of Partcles Rectlnear oton Straght-Lne oton [.1] Analtcal solutons for poston/veloct [.1] Solvng equatons of oton Analtcal solutons (1 D revew) [.1] Nuercal solutons [.1] Nuercal ntegraton
More informationWhat is LP? LP is an optimization technique that allocates limited resources among competing activities in the best possible manner.
(C) 998 Gerald B Sheblé, all rghts reserved Lnear Prograng Introducton Contents I. What s LP? II. LP Theor III. The Splex Method IV. Refneents to the Splex Method What s LP? LP s an optzaton technque that
More informationExcess Error, Approximation Error, and Estimation Error
E0 370 Statstcal Learnng Theory Lecture 10 Sep 15, 011 Excess Error, Approxaton Error, and Estaton Error Lecturer: Shvan Agarwal Scrbe: Shvan Agarwal 1 Introducton So far, we have consdered the fnte saple
More informationChapter 12. Ordinary Differential Equation Boundary Value (BV) Problems
Chapter. Ordnar Dfferental Equaton Boundar Value (BV) Problems In ths chapter we wll learn how to solve ODE boundar value problem. BV ODE s usuall gven wth x beng the ndependent space varable. p( x) q(
More informationITERATIVE ESTIMATION PROCEDURE FOR GEOSTATISTICAL REGRESSION AND GEOSTATISTICAL KRIGING
ESE 5 ITERATIVE ESTIMATION PROCEDURE FOR GEOSTATISTICAL REGRESSION AND GEOSTATISTICAL KRIGING Gven a geostatstcal regresson odel: k Y () s x () s () s x () s () s, s R wth () unknown () E[ ( s)], s R ()
More informationBoundary Value Problems. Lecture Objectives. Ch. 27
Boundar Vaue Probes Ch. 7 Lecture Obectves o understand the dfference between an nta vaue and boundar vaue ODE o be abe to understand when and how to app the shootng ethod and FD ethod. o understand what
More informationLeast Squares Fitting of Data
Least Squares Fttng of Data Davd Eberly Geoetrc Tools, LLC http://www.geoetrctools.co/ Copyrght c 1998-2014. All Rghts Reserved. Created: July 15, 1999 Last Modfed: February 9, 2008 Contents 1 Lnear Fttng
More informationShuai Dong. Isaac Newton. Gottfried Leibniz
Computatonal pyscs Sua Dong Isaac Newton Gottred Lebnz Numercal calculus poston dervatve ntegral v velocty dervatve ntegral a acceleraton Numercal calculus Numercal derentaton Numercal ntegraton Roots
More information2 Finite difference basics
Numersche Methoden 1, WS 11/12 B.J.P. Kaus 2 Fnte dfference bascs Consder the one- The bascs of the fnte dfference method are best understood wth an example. dmensonal transent heat conducton equaton T
More informationLinear discriminants. Nuno Vasconcelos ECE Department, UCSD
Lnear dscrmnants Nuno Vasconcelos ECE Department UCSD Classfcaton a classfcaton problem as to tpes of varables e.g. X - vector of observatons features n te orld Y - state class of te orld X R 2 fever blood
More informationOur focus will be on linear systems. A system is linear if it obeys the principle of superposition and homogenity, i.e.
SSTEM MODELLIN In order to solve a control syste proble, the descrptons of the syste and ts coponents ust be put nto a for sutable for analyss and evaluaton. The followng ethods can be used to odel physcal
More informationBAYESIAN CURVE FITTING USING PIECEWISE POLYNOMIALS. Dariusz Biskup
BAYESIAN CURVE FITTING USING PIECEWISE POLYNOMIALS Darusz Bskup 1. Introducton The paper presents a nonparaetrc procedure for estaton of an unknown functon f n the regresson odel y = f x + ε = N. (1) (
More informationELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM
ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM An elastc wave s a deformaton of the body that travels throughout the body n all drectons. We can examne the deformaton over a perod of tme by fxng our look
More informationSolving Fuzzy Linear Programming Problem With Fuzzy Relational Equation Constraint
Intern. J. Fuzz Maeatcal Archve Vol., 0, -0 ISSN: 0 (P, 0 0 (onlne Publshed on 0 Septeber 0 www.researchasc.org Internatonal Journal of Solvng Fuzz Lnear Prograng Proble W Fuzz Relatonal Equaton Constrant
More informationMultigrid Methods and Applications in CFD
Multgrd Metods and Applcatons n CFD Mcael Wurst 0 May 009 Contents Introducton Typcal desgn of CFD solvers 3 Basc metods and ter propertes for solvng lnear systems of equatons 4 Geometrc Multgrd 3 5 Algebrac
More informationSolutions to GW Flow Equations Part 2- Finite Difference Approximations
مرکز تحقیقات آببهای زیرزمینی متا ب Groundwater Researc Center GRC Solutons to GW Flow Equatons Part - Fnte Dfference ppromatons حاکم بر جریان آبهای زیرزمینی معادلات نیومریکال نیومریکال عددی حل Readng:
More informationCENTROID (AĞIRLIK MERKEZİ )
CENTOD (ĞLK MEKEZİ ) centrod s a geometrcal concept arsng from parallel forces. Tus, onl parallel forces possess a centrod. Centrod s tougt of as te pont were te wole wegt of a pscal od or sstem of partcles
More informationDenote the function derivatives f(x) in given points. x a b. Using relationships (1.2), polynomials (1.1) are written in the form
SET OF METHODS FO SOUTION THE AUHY POBEM FO STIFF SYSTEMS OF ODINAY DIFFEENTIA EUATIONS AF atypov and YuV Nulchev Insttute of Theoretcal and Appled Mechancs SB AS 639 Novosbrs ussa Introducton A constructon
More informationLeast Squares Fitting of Data
Least Squares Fttng of Data Davd Eberly Geoetrc Tools, LLC http://www.geoetrctools.co/ Copyrght c 1998-2015. All Rghts Reserved. Created: July 15, 1999 Last Modfed: January 5, 2015 Contents 1 Lnear Fttng
More informationThe Finite Element Method: A Short Introduction
Te Fnte Element Metod: A Sort ntroducton Wat s FEM? Te Fnte Element Metod (FEM) ntroduced by engneers n late 50 s and 60 s s a numercal tecnque for solvng problems wc are descrbed by Ordnary Dfferental
More informationCENTROID (AĞIRLIK MERKEZİ )
CENTOD (ĞLK MEKEZİ ) centrod s a geometrcal concept arsng from parallel forces. Tus, onl parallel forces possess a centrod. Centrod s tougt of as te pont were te wole wegt of a pscal od or sstem of partcles
More informationLecture 21: Numerical methods for pricing American type derivatives
Lecture 21: Numercal methods for prcng Amercan type dervatves Xaoguang Wang STAT 598W Aprl 10th, 2014 (STAT 598W) Lecture 21 1 / 26 Outlne 1 Fnte Dfference Method Explct Method Penalty Method (STAT 598W)
More informationMultistage Homotopy Analysis Method for Solving Nonlinear Integral Equations
Avalable at ttp://pvau.edu/aa Appl. Appl. Mat. ISSN: 9-966 Specal Issue No. (August ) pp. 5 Applcatons and Appled Mateatcs: An Internatonal Journal (AAM) Multstage Hootop Analss Metod for Solvng Nonlnear
More informationMachine Learning. What is a good Decision Boundary? Support Vector Machines
Machne Learnng 0-70/5 70/5-78 78 Sprng 200 Support Vector Machnes Erc Xng Lecture 7 March 5 200 Readng: Chap. 6&7 C.B book and lsted papers Erc Xng @ CMU 2006-200 What s a good Decson Boundar? Consder
More informationMain components of the above cycle are: 1) Boiler (steam generator) heat exchanger 2) Turbine generates work 3) Condenser heat exchanger 4) Pump
Introducton to Terodynacs, Lecture -5 Pro. G. Cccarell (0 Applcaton o Control olue Energy Analyss Most terodynac devces consst o a seres o coponents operatng n a cycle, e.g., stea power plant Man coponents
More informationLecture 12: Discrete Laplacian
Lecture 12: Dscrete Laplacan Scrbe: Tanye Lu Our goal s to come up wth a dscrete verson of Laplacan operator for trangulated surfaces, so that we can use t n practce to solve related problems We are mostly
More informationLECTURE :FACTOR ANALYSIS
LCUR :FACOR ANALYSIS Rta Osadchy Based on Lecture Notes by A. Ng Motvaton Dstrbuton coes fro MoG Have suffcent aount of data: >>n denson Use M to ft Mture of Gaussans nu. of tranng ponts If
More informationMachine Learning. Support Vector Machines. Eric Xing. Lecture 4, August 12, Reading: Eric CMU,
Machne Learnng Support Vector Machnes Erc Xng Lecture 4 August 2 200 Readng: Erc Xng @ CMU 2006-200 Erc Xng @ CMU 2006-200 2 What s a good Decson Boundar? Wh e a have such boundares? Irregular dstrbuton
More informationApplied Mathematics Letters
Appled Matheatcs Letters 2 (2) 46 5 Contents lsts avalable at ScenceDrect Appled Matheatcs Letters journal hoepage: wwwelseverco/locate/al Calculaton of coeffcents of a cardnal B-splne Gradr V Mlovanovć
More information1 Definition of Rademacher Complexity
COS 511: Theoretcal Machne Learnng Lecturer: Rob Schapre Lecture #9 Scrbe: Josh Chen March 5, 2013 We ve spent the past few classes provng bounds on the generalzaton error of PAClearnng algorths for the
More informationComputational and Statistical Learning theory Assignment 4
Coputatonal and Statstcal Learnng theory Assgnent 4 Due: March 2nd Eal solutons to : karthk at ttc dot edu Notatons/Defntons Recall the defnton of saple based Radeacher coplexty : [ ] R S F) := E ɛ {±}
More informationME 501A Seminar in Engineering Analysis Page 1
umercal Solutons of oundary-value Problems n Os ovember 7, 7 umercal Solutons of oundary- Value Problems n Os Larry aretto Mechancal ngneerng 5 Semnar n ngneerng nalyss ovember 7, 7 Outlne Revew stff equaton
More informationSolution of Linear System of Equations and Matrix Inversion Gauss Seidel Iteration Method
Soluton of Lnear System of Equatons and Matr Inverson Gauss Sedel Iteraton Method It s another well-known teratve method for solvng a system of lnear equatons of the form a + a22 + + ann = b a2 + a222
More informationLimit Cycle Bifurcations in a Class of Cubic System near a Nilpotent Center *
Appled Mateatcs 77-777 ttp://dxdoorg/6/a75 Publsed Onlne July (ttp://wwwscrporg/journal/a) Lt Cycle Bfurcatons n a Class of Cubc Syste near a Nlpotent Center * Jao Jang Departent of Mateatcs Sanga Marte
More informationRelaxation Methods for Iterative Solution to Linear Systems of Equations
Relaxaton Methods for Iteratve Soluton to Lnear Systems of Equatons Gerald Recktenwald Portland State Unversty Mechancal Engneerng Department gerry@pdx.edu Overvew Techncal topcs Basc Concepts Statonary
More informationCOS 511: Theoretical Machine Learning
COS 5: Theoretcal Machne Learnng Lecturer: Rob Schapre Lecture #0 Scrbe: José Sões Ferrera March 06, 203 In the last lecture the concept of Radeacher coplexty was ntroduced, wth the goal of showng that
More informationSlobodan Lakić. Communicated by R. Van Keer
Serdca Math. J. 21 (1995), 335-344 AN ITERATIVE METHOD FOR THE MATRIX PRINCIPAL n-th ROOT Slobodan Lakć Councated by R. Van Keer In ths paper we gve an teratve ethod to copute the prncpal n-th root and
More information1 GSW Iterative Techniques for y = Ax
1 for y = A I m gong to cheat here. here are a lot of teratve technques that can be used to solve the general case of a set of smultaneous equatons (wrtten n the matr form as y = A), but ths chapter sn
More informationME 412 Heat Transfer Laboratory
ME Heat Transfer Laboratory ME Heat Transfer Laboratory andle owered Desalnator nalytcal Modelng roperty Inforaton oposed of paraffn wax wc s a petroleu product consstng anly of noral alkanes. Molecular
More informationPART 8. Partial Differential Equations PDEs
he Islamc Unverst of Gaza Facult of Engneerng Cvl Engneerng Department Numercal Analss ECIV 3306 PAR 8 Partal Dfferental Equatons PDEs Chapter 9; Fnte Dfference: Ellptc Equatons Assocate Prof. Mazen Abualtaef
More informationCOMP4630: λ-calculus
COMP4630: λ-calculus 4. Standardsaton Mcael Norrs Mcael.Norrs@ncta.com.au Canberra Researc Lab., NICTA Semester 2, 2015 Last Tme Confluence Te property tat dvergent evaluatons can rejon one anoter Proof
More informationPhysics 3A: Linear Momentum. Physics 3A: Linear Momentum. Physics 3A: Linear Momentum. Physics 3A: Linear Momentum
Recall that there was ore to oton than just spee A ore coplete escrpton of oton s the concept of lnear oentu: p v (8.) Beng a prouct of a scalar () an a vector (v), oentu s a vector: p v p y v y p z v
More informationElastic Collisions. Definition: two point masses on which no external forces act collide without losing any energy.
Elastc Collsons Defnton: to pont asses on hch no external forces act collde thout losng any energy v Prerequstes: θ θ collsons n one denson conservaton of oentu and energy occurs frequently n everyday
More informationFinite Difference Method
7/0/07 Instructor r. Ramond Rump (9) 747 698 rcrump@utep.edu EE 337 Computatonal Electromagnetcs (CEM) Lecture #0 Fnte erence Method Lecture 0 These notes ma contan coprghted materal obtaned under ar use
More informationSeveral generation methods of multinomial distributed random number Tian Lei 1, a,linxihe 1,b,Zhigang Zhang 1,c
Internatonal Conference on Appled Scence and Engneerng Innovaton (ASEI 205) Several generaton ethods of ultnoal dstrbuted rando nuber Tan Le, a,lnhe,b,zhgang Zhang,c School of Matheatcs and Physcs, USTB,
More informationStanford University CS254: Computational Complexity Notes 7 Luca Trevisan January 29, Notes for Lecture 7
Stanford Unversty CS54: Computatonal Complexty Notes 7 Luca Trevsan January 9, 014 Notes for Lecture 7 1 Approxmate Countng wt an N oracle We complete te proof of te followng result: Teorem 1 For every
More informationOn the number of regions in an m-dimensional space cut by n hyperplanes
6 On the nuber of regons n an -densonal space cut by n hyperplanes Chungwu Ho and Seth Zeran Abstract In ths note we provde a unfor approach for the nuber of bounded regons cut by n hyperplanes n general
More informationSpecial Relativity and Riemannian Geometry. Department of Mathematical Sciences
Tutoral Letter 06//018 Specal Relatvty and Reannan Geoetry APM3713 Seester Departent of Matheatcal Scences IMPORTANT INFORMATION: Ths tutoral letter contans the solutons to Assgnent 06. BAR CODE Learn
More informationModule 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur
Module 3 LOSSY IMAGE COMPRESSION SYSTEMS Verson ECE IIT, Kharagpur Lesson 6 Theory of Quantzaton Verson ECE IIT, Kharagpur Instructonal Objectves At the end of ths lesson, the students should be able to:
More informationGrid Generation around a Cylinder by Complex Potential Functions
Research Journal of Appled Scences, Engneerng and Technolog 4(): 53-535, 0 ISSN: 040-7467 Mawell Scentfc Organzaton, 0 Submtted: December 0, 0 Accepted: Januar, 0 Publshed: June 0, 0 Grd Generaton around
More information1. Statement of the problem
Volue 14, 010 15 ON THE ITERATIVE SOUTION OF A SYSTEM OF DISCRETE TIMOSHENKO EQUATIONS Peradze J. and Tsklaur Z. I. Javakhshvl Tbls State Uversty,, Uversty St., Tbls 0186, Georga Georgan Techcal Uversty,
More informationThe Fundamental Theorem of Algebra. Objective To use the Fundamental Theorem of Algebra to solve polynomial equations with complex solutions
5-6 The Fundamental Theorem of Algebra Content Standards N.CN.7 Solve quadratc equatons wth real coeffcents that have comple solutons. N.CN.8 Etend polnomal denttes to the comple numbers. Also N.CN.9,
More informationFinal Exam Solutions, 1998
58.439 Fnal Exa Solutons, 1998 roble 1 art a: Equlbru eans that the therodynac potental of a consttuent s the sae everywhere n a syste. An exaple s the Nernst potental. If the potental across a ebrane
More informationSome Comments on Accelerating Convergence of Iterative Sequences Using Direct Inversion of the Iterative Subspace (DIIS)
Some Comments on Acceleratng Convergence of Iteratve Sequences Usng Drect Inverson of the Iteratve Subspace (DIIS) C. Davd Sherrll School of Chemstry and Bochemstry Georga Insttute of Technology May 1998
More informationResource Allocation and Decision Analysis (ECON 8010) Spring 2014 Foundations of Regression Analysis
Resource Allocaton and Decson Analss (ECON 800) Sprng 04 Foundatons of Regresson Analss Readng: Regresson Analss (ECON 800 Coursepak, Page 3) Defntons and Concepts: Regresson Analss statstcal technques
More informationCME 302: NUMERICAL LINEAR ALGEBRA FALL 2005/06 LECTURE 13
CME 30: NUMERICAL LINEAR ALGEBRA FALL 005/06 LECTURE 13 GENE H GOLUB 1 Iteratve Methods Very large problems (naturally sparse, from applcatons): teratve methods Structured matrces (even sometmes dense,
More informationAppendix B. The Finite Difference Scheme
140 APPENDIXES Appendx B. The Fnte Dfference Scheme In ths appendx we present numercal technques whch are used to approxmate solutons of system 3.1 3.3. A comprehensve treatment of theoretcal and mplementaton
More informationNumerical Solution of Ordinary Differential Equations
College of Engneerng and Computer Scence Mecancal Engneerng Department otes on Engneerng Analss Larr Caretto ovember 9, 7 Goal of tese notes umercal Soluton of Ordnar Dfferental Equatons ese notes were
More informatione i is a random error
Chapter - The Smple Lnear Regresson Model The lnear regresson equaton s: where + β + β e for,..., and are observable varables e s a random error How can an estmaton rule be constructed for the unknown
More informationNumerical Heat and Mass Transfer
Master degree n Mechancal Engneerng Numercal Heat and Mass Transfer 06-Fnte-Dfference Method (One-dmensonal, steady state heat conducton) Fausto Arpno f.arpno@uncas.t Introducton Why we use models and
More informationLecture 26 Finite Differences and Boundary Value Problems
4//3 Leture 6 Fnte erenes and Boundar Value Problems Numeral derentaton A nte derene s an appromaton o a dervatve - eample erved rom Talor seres 3 O! Negletng all terms ger tan rst order O O Tat s te orward
More information6.3.4 Modified Euler s method of integration
6.3.4 Modfed Euler s method of ntegraton Before dscussng the applcaton of Euler s method for solvng the swng equatons, let us frst revew the basc Euler s method of numercal ntegraton. Let the general from
More informationLinear Approximation with Regularization and Moving Least Squares
Lnear Approxmaton wth Regularzaton and Movng Least Squares Igor Grešovn May 007 Revson 4.6 (Revson : March 004). 5 4 3 0.5 3 3.5 4 Contents: Lnear Fttng...4. Weghted Least Squares n Functon Approxmaton...
More informationPreference and Demand Examples
Dvson of the Huantes and Socal Scences Preference and Deand Exaples KC Border October, 2002 Revsed Noveber 206 These notes show how to use the Lagrange Karush Kuhn Tucker ultpler theores to solve the proble
More information8.6 The Complex Number System
8.6 The Complex Number System Earler n the chapter, we mentoned that we cannot have a negatve under a square root, snce the square of any postve or negatve number s always postve. In ths secton we want
More information1.3 Hence, calculate a formula for the force required to break the bond (i.e. the maximum value of F)
EN40: Dynacs and Vbratons Hoework 4: Work, Energy and Lnear Moentu Due Frday March 6 th School of Engneerng Brown Unversty 1. The Rydberg potental s a sple odel of atoc nteractons. It specfes the potental
More informationAN ANALYSIS OF A FRACTAL KINETICS CURVE OF SAVAGEAU
AN ANALYI OF A FRACTAL KINETIC CURE OF AAGEAU by John Maloney and Jack Hedel Departent of Matheatcs Unversty of Nebraska at Oaha Oaha, Nebraska 688 Eal addresses: aloney@unoaha.edu, jhedel@unoaha.edu Runnng
More informationCorrelation and Regression. Correlation 9.1. Correlation. Chapter 9
Chapter 9 Correlaton and Regresson 9. Correlaton Correlaton A correlaton s a relatonshp between two varables. The data can be represented b the ordered pars (, ) where s the ndependent (or eplanator) varable,
More informationChapter Newton s Method
Chapter 9. Newton s Method After readng ths chapter, you should be able to:. Understand how Newton s method s dfferent from the Golden Secton Search method. Understand how Newton s method works 3. Solve
More informationDigital Signal Processing
Dgtal Sgnal Processng Dscrete-tme System Analyss Manar Mohasen Offce: F8 Emal: manar.subh@ut.ac.r School of IT Engneerng Revew of Precedent Class Contnuous Sgnal The value of the sgnal s avalable over
More informationMultivariate Ratio Estimator of the Population Total under Stratified Random Sampling
Open Journal of Statstcs, 0,, 300-304 ttp://dx.do.org/0.436/ojs.0.3036 Publsed Onlne July 0 (ttp://www.scrp.org/journal/ojs) Multvarate Rato Estmator of te Populaton Total under Stratfed Random Samplng
More informationand Daniel B. Szyld Report October 2017, Revised April 2018
Syncronous and Asyncronous Optzed Scwarz Metods for Posson s Equaton n Rectangular Doans José Garay, Frédérc Magoulès, and Danel B. Szyld Report 7--8 October 7, Revsed Aprl 8 Departent of Mateatcs Teple
More informationNUMERICAL DIFFERENTIATION
NUMERICAL DIFFERENTIATION 1 Introducton Dfferentaton s a method to compute the rate at whch a dependent output y changes wth respect to the change n the ndependent nput x. Ths rate of change s called the
More informationNumerical Solution of Ordinary Differential Equations
Numercal Methods (CENG 00) CHAPTER-VI Numercal Soluton of Ordnar Dfferental Equatons 6 Introducton Dfferental equatons are equatons composed of an unknown functon and ts dervatves The followng are examples
More informationWeek3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity
Week3, Chapter 4 Moton n Two Dmensons Lecture Quz A partcle confned to moton along the x axs moves wth constant acceleraton from x =.0 m to x = 8.0 m durng a 1-s tme nterval. The velocty of the partcle
More information1 Review From Last Time
COS 5: Foundatons of Machne Learnng Rob Schapre Lecture #8 Scrbe: Monrul I Sharf Aprl 0, 2003 Revew Fro Last Te Last te, we were talkng about how to odel dstrbutons, and we had ths setup: Gven - exaples
More informationA PROCEDURE FOR SIMULATING THE NONLINEAR CONDUCTION HEAT TRANSFER IN A BODY WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY.
Proceedngs of the th Brazlan Congress of Thermal Scences and Engneerng -- ENCIT 006 Braz. Soc. of Mechancal Scences and Engneerng -- ABCM, Curtba, Brazl,- Dec. 5-8, 006 A PROCEDURE FOR SIMULATING THE NONLINEAR
More informationFermi-Dirac statistics
UCC/Physcs/MK/EM/October 8, 205 Fer-Drac statstcs Fer-Drac dstrbuton Matter partcles that are eleentary ostly have a type of angular oentu called spn. hese partcles are known to have a agnetc oent whch
More informationXII.3 The EM (Expectation-Maximization) Algorithm
XII.3 The EM (Expectaton-Maxzaton) Algorth Toshnor Munaata 3/7/06 The EM algorth s a technque to deal wth varous types of ncoplete data or hdden varables. It can be appled to a wde range of learnng probles
More informationPHYS 1443 Section 002 Lecture #20
PHYS 1443 Secton 002 Lecture #20 Dr. Jae Condtons for Equlbru & Mechancal Equlbru How to Solve Equlbru Probles? A ew Exaples of Mechancal Equlbru Elastc Propertes of Solds Densty and Specfc Gravty lud
More informationPre-Calculus Summer Assignment
Pre-Calculus Summer Assgnment Dear Future Pre-Calculus Student, Congratulatons on our successful completon of Algebra! Below ou wll fnd the summer assgnment questons. It s assumed that these concepts,
More informationNormally, in one phase reservoir simulation we would deal with one of the following fluid systems:
TPG4160 Reservor Smulaton 2017 page 1 of 9 ONE-DIMENSIONAL, ONE-PHASE RESERVOIR SIMULATION Flud systems The term sngle phase apples to any system wth only one phase present n the reservor In some cases
More information5 The Laplace Equation in a convex polygon
5 Te Laplace Equaton n a convex polygon Te most mportant ellptc PDEs are te Laplace, te modfed Helmoltz and te Helmoltz equatons. Te Laplace equaton s u xx + u yy =. (5.) Te real and magnary parts of an
More informationLecture Note 3. Eshelby s Inclusion II
ME340B Elastcty of Mcroscopc Structures Stanford Unversty Wnter 004 Lecture Note 3. Eshelby s Incluson II Chrs Wenberger and We Ca c All rghts reserved January 6, 004 Contents 1 Incluson energy n an nfnte
More informationParallel VLSI CAD Algorithms. Lecture 7 Iterative methods for IC thermal analysis Zhuo Feng
Parallel VLSI CAD Algortms Lecture 7 Iteratve metods for IC termal analyss Zuo Feng 7. We ve taled about usng LU to solve a sparse lnear matr t problem Anoter way of dealng wt sparse matr problems s to
More informationUsing T.O.M to Estimate Parameter of distributions that have not Single Exponential Family
IOSR Journal of Mathematcs IOSR-JM) ISSN: 2278-5728. Volume 3, Issue 3 Sep-Oct. 202), PP 44-48 www.osrjournals.org Usng T.O.M to Estmate Parameter of dstrbutons that have not Sngle Exponental Famly Jubran
More informationUNIVERSITY OF BOLTON RAK ACADEMIC CENTRE BENG(HONS) MECHANICAL ENGINEERING SEMESTER TWO EXAMINATION 2017/2018 FINITE ELEMENT AND DIFFERENCE SOLUTIONS
OCD0 UNIVERSITY OF BOLTON RAK ACADEMIC CENTRE BENG(HONS) MECHANICAL ENGINEERING SEMESTER TWO EXAMINATION 07/08 FINITE ELEMENT AND DIFFERENCE SOLUTIONS MODULE NO. AME6006 Date: Wednesda 0 Ma 08 Tme: 0:00
More informationTransfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system
Transfer Functons Convenent representaton of a lnear, dynamc model. A transfer functon (TF) relates one nput and one output: x t X s y t system Y s The followng termnology s used: x y nput output forcng
More informationPhysics 2A Chapter 3 HW Solutions
Phscs A Chapter 3 HW Solutons Chapter 3 Conceptual Queston: 4, 6, 8, Problems: 5,, 8, 7, 3, 44, 46, 69, 70, 73 Q3.4. Reason: (a) C = A+ B onl A and B are n the same drecton. Sze does not matter. (b) C
More information4DVAR, according to the name, is a four-dimensional variational method.
4D-Varatonal Data Assmlaton (4D-Var) 4DVAR, accordng to the name, s a four-dmensonal varatonal method. 4D-Var s actually a drect generalzaton of 3D-Var to handle observatons that are dstrbuted n tme. The
More informationImage classification. Given the bag-of-features representations of images from different classes, how do we learn a model for distinguishing i them?
Image classfcaton Gven te bag-of-features representatons of mages from dfferent classes ow do we learn a model for dstngusng tem? Classfers Learn a decson rule assgnng bag-offeatures representatons of
More informationCHAPTER 4. Vector Spaces
man 2007/2/16 page 234 CHAPTER 4 Vector Spaces To crtcze mathematcs for ts abstracton s to mss the pont entrel. Abstracton s what makes mathematcs work. Ian Stewart The man am of ths tet s to stud lnear
More informationThe Karush-Kuhn-Tucker. Nuno Vasconcelos ECE Department, UCSD
e Karus-Kun-ucker condtons and dualt Nuno Vasconcelos ECE Department, UCSD Optmzaton goal: nd mamum or mnmum o a uncton Denton: gven unctons, g, 1,...,k and, 1,...m dened on some doman Ω R n mn w, w Ω
More informationInductance Calculation for Conductors of Arbitrary Shape
CRYO/02/028 Aprl 5, 2002 Inductance Calculaton for Conductors of Arbtrary Shape L. Bottura Dstrbuton: Internal Summary In ths note we descrbe a method for the numercal calculaton of nductances among conductors
More information3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X
Statstcs 1: Probablty Theory II 37 3 EPECTATION OF SEVERAL RANDOM VARIABLES As n Probablty Theory I, the nterest n most stuatons les not on the actual dstrbuton of a random vector, but rather on a number
More information1 Proving the Fundamental Theorem of Statistical Learning
THEORETICAL MACHINE LEARNING COS 5 LECTURE #7 APRIL 5, 6 LECTURER: ELAD HAZAN NAME: FERMI MA ANDDANIEL SUO oving te Fundaental Teore of Statistical Learning In tis section, we prove te following: Teore.
More informationDirect Methods for Solving Macromolecular Structures Ed. by S. Fortier Kluwer Academic Publishes, The Netherlands, 1998, pp
Drect Metods for Solvng Macromolecular Structures Ed. by S. Forter Kluwer Academc Publses, Te Neterlands, 998, pp. 79-85. SAYRE EQUATION, TANGENT FORMULA AND SAYTAN FAN HAI-FU Insttute of Pyscs, Cnese
More informationLectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix
Lectures - Week 4 Matrx norms, Condtonng, Vector Spaces, Lnear Independence, Spannng sets and Bass, Null space and Range of a Matrx Matrx Norms Now we turn to assocatng a number to each matrx. We could
More information