Computational Geodynamics: Advection equations and the art of numerical modeling
|
|
- Brianna Pierce
- 5 years ago
- Views:
Transcription
1 540 Geodynamcs Unversty of Southern Calforna, Los Angeles Computatonal Geodynamcs: Advecton equatons and the art of numercal modelng Sofar we manly focussed on dffuson equaton n a non-movng doman. Ths s maybe relevant for the case of a dke ntruson or for a lthosphere whch remans undeformed. However there are also cases where materal moves. An eample s a plume rsng through a convectng mantle. The plume s hot and hence ts densty s low compared to the colder mantle around t. The hot materal rses wth a gven velocty (when term-projects are fnshed we wll have codes that gve us those veloctes for a gven densty dstrbuton). If the numercal grd remans fed n the background, the hot temperatures should be moved to dfferent grdponts (see fg. 1 for an llustraton of ths effect). Fgure 1: Snapshots of a bottom heated thermal convecton model wth a Raylegh-number of and constant vscosty (no nternal heatng). Temperature s advected through a fed (euleran) grd (crcles) wth a velocty (arrows) that s computed wth a Stokes solver. Mathematcally, the temperature equaton gets an addtonal term and n 1-D the dffuson-advecton equaton becomes (see equaton cheat sheet) ( ) ρc p t + v = ( k ) (1) or n 2-D ( ) ρc p t + v + v z = ( k ) + ( k ) (2) z z z where v,v z are veloctes n -,respectvely z-drecton. If dffuson s absent (.e. k = 0), the advecton equatons are t + v = 0 (3) and t + v + v z z = 0 (4) These equatons should also be solved f you want to compute the settlng of a heavy partcle n a flud, or n weather predcton codes. In the current lecture we re gong to look at some optons on how to solve these equatons wth a fnte dfference scheme on a fed grd. Even though the equatons are farly smple, t s far from smple to solve them accurately. If not done carefully, you ll up wth strong numercal artfacts such as wggles and numercal dffuson.
2 % advecton_central_1d % clear all n = 201; W = 40; % wdth of doman Vel = -4; % velocty sgma = 1; Ampl = 2; nt = 500; % number of tmesteps dt = 1e-2; d = W/(n-1); = 0:d:W; % Intal gaussan T-profle c = 20; T = Ampl*ep(-(-c).ˆ2/sgmaˆ2); % Velocty V = ones(1,n)*vel; % Central fnte dfference dscretzaton for tme=1:nt % central fn. dff for =2:n-1 Tnew() =?? Tnew(1) = T(1); Tnew(n) = T(n); T = Tnew; tme = tme*dt; % Analytcal soluton for ths case T_anal = Ampl*ep(-(-c-tme*Vel).ˆ2/sgmaˆ2); fgure(1),clf, plot(,t,,t_anal), drawnow leg( Numercal, Analytcal ) Fgure 2: MATLAB scrpt to be used wth eercse 1. FTCS method In a 1-D case the smplest way to dscretze equaton 3 s by employng a central fnte dfference scheme (also called the forward-tme, central space or FTCS scheme): T n+1 T n = v, T n +1 T n 1 2 (5) Eercse 1 Program the FTCS method n the code of fgure 2 and watch what happens. Change the sgn of the velocty. Change the tmestep and grdspacng and compute the nondmensonal parameter α = v /. When do unstable results occur? As you saw from the eercse, the FTCS method does not work... In fact t s a nce eample of a scheme that looks nce (and logcal!) on paper, but sucks [you can actually mathematcally prove ths by a socalled von Neuman stablty analyss, f you re nterested you should read chapter 5 of Marc Spegelman s scrpt on the course webpage; t s related to the fact that ths scheme produces ant-dffuson, whch s numercally nstable]. La method Some day, long ago, there was a guy called La and he thought: well let s replace the T n n the tme-dervatve of equaton 5 wth (T+1 n + T 1 n )/2. The resultng equaton s T n+1 (T+1 n + T 1 n )/2 T+1 n = v T n, 2 1 (6) 2
3 Eercse 2 Program the La method by modfyng the scrpt of the last eercse. Play wth dfferent veloctes and compute the Courant number α, whch s gven by the followng equaton: α = v (7) Is the numercal scheme stable for all Courant numbers? As you saw from the eercse the La method doesn t blow up, but does have a lot of numercal dffuson. So t s an mprovement but t s stll not perfect, snce you may loose the plumes of fgure 1 around md-mantle purely due to numercal dffuson (and than you ll probably fnd some fancy geochemcal eplanaton for ths...). I hope that you also slowly start to see that numercal modelng s a bt of an art... Upwnd scheme The net scheme that people came up wth s the so-called upwnd scheme. Here, the spatal fnte dfference scheme deps on the sgn of the velocty: T n+1 T n = v, { T n T n 1, f v, > 0 T+1 n T n, f v, < 0 (8) Eercse 3 Program the upwnd scheme method. Play wth dfferent veloctes and compute the Courant number α Is the numercal scheme stable for all Courant numbers? Also the upwnd scheme suffers from numercal dffuson. Sofar we employed eplct dscretzatons. You re probably wonderng whether mplct dscretzatons wll save us agan ths tme. Bad news: doesn t work (try f you want). So we have to come up wth somethng else. Staggered Leapfrog The eplct dscretzatons dscussed sofar were second order accurate n tme, but only frst order n space. We can also come up wth a scheme that s second order n tme and space T n+1 T n 1 T+1 n = v T n, 2 2 the dffculty n ths scheme s that we ll have to store two tmestep levels: both T n 1 and T n (but m sure you ll manage). Eercse 4 Program the staggered leapfrog method (assume that at the frst tmestep T n 1 = T n ). Play around wth dfferent values of the Courant number α. Also make the wdth of the gaussan curve smaller. The staggered leapfrog method works qute well regardng the ampltude and wavespeed as long as α s close to one. If however α << 1 and the lengthscale of the to-be-transported quantty s small compared to the number of grdponts (e.g. we have a thn plume), numercal oscllatons agan occur. Ths s typcally the case n mantle convecton smulatons (see Fg. 1), so agan we have a numercal scheme that s not great for mantle convecton smulatons. 1 (9) 3
4 MPDATA Ths s a technque that s frequently appled n (older) mantle convecton codes. The basc dea s based on an dea of Smolarkewcz and s bascally a modfcaton of the upwnd scheme, where some teratve correctons are appled. The results are qute ok. However t s somewhat more complcated to mplement (the system that s solved s actually nonlnear and requres teratons, whch we dd not dscuss yet). Moreover we stll have a restrcton on the tmestep (gven by the Courant crterum). For ths reason we won t go nto detal here (see Spegelman s lecture notes, chapter 5 f you re nterested). Sem-Lagrangan So what we want s a scheme that doesn t blow up, has only small numercal dffuson and that s not lmted by the Courant crterum. It actually ests and s called the sem-lagrangan method. I m a bg fan of t and so s Marc Spegelman (see hs course notes), and people that do clmate modellng or numercal weather predcton. It s however a bt werd and has lttle to do wth the fnte dfference schemes we dscussed sofar. Snce ths scheme s however the one that s maybe most mportant for practcal purposes we wll go n more detal here. Basc dea The basc dea of the sem-lagrangan method s llustrated n fgure 3A and s gven by the followng For each pont 1. Assume that the future velocty v (n + 1,) s known. Under the assumpton that the velocty at the old tmestep s almost dentcal to the future velocty (v (n+1,) = v (n,)) and that velocty doesn t vary spatally (v (n, 1) = v (n,) = v (n,+1)), we can compute the locaton X where the partcle came from by X= () v (n+ 1,) 2. Interpolate your temperature from grdponts to the locaton X at tme n. Use cubc nterpolaton (n MATLAB use the command nterp1(,t, pont, cubc ) for nterpolaton). Now you smply say that T (n+1,) = T (n,x), and you re done. Note that ths scheme assumes that no heat-sources were actve durng the advecton of T from T (n,x) to T (n + 1,). If heat sources are present and are spatally varable, some etra stuff needs to be done (read Marc s lecture notes). Eercse 5 Program the sem-lagrangan advecton scheme llustrated n fgure 3A. Is there a Courant crterum for stablty? Some mprovements The algorthm descrbed n fgure 3A llustrates the basc dea of the sem-lagrangan scheme. However t has two problems. Frst t assumes that velocty s spatally constant (whch s clearly not the case n mantle convecton smulatons). Second, t assumes that velocty doesn t change between tme n and n + 1. We can overcome both problems by usng a more accurate tmesteppng algorthm. An eample s an teratve md-pont scheme whch works as follows (see fgure 3B): For each pont 1. Use the velocty v (n + 1,) to compute the locaton at tme n + 1/2 (n other words: take half a tmestep backwards n tme). 2. Fnd the velocty at the locaton at half tmestep n + 1/2. Assume that the velocty at the half-tmestep can be computed as v (n + 1/2,) = v (n+1,)+v (n,) 2. Use lnear nterpolaton for the spatal nterpolaton of velocty v (n + 1/2,). 4
5 A tme t n+1 n T(n,-2) v (n,-1) T(n,-1) v (n+1,) v (n,) T(n,) B tme t n+1 n+1/2 n v (n,-2) true path v (n,-1) v (n+1,) v (n,) space space Fgure 3: Bascs of the sem-lagrangan method. See tet for eplanaton. 3. Go back to pont 1, but use the velocty v (n + 1/2,) nstead of v (n + 1,) to move the pont (n + 1,) backwards n tme. Repeat ths process a number of tmes (for eample 5 tmes). Ths gves a farly accurate centered velocty. 4. Compute the locaton X wth the centered velocty v (n + 1/2,), wth X= () v (n + 1/2,). 5. Use cubc nterpolaton to fnd the temperature at pont X. Other methods are also possble. A popular one s the fourth order Runga-Kutta method, whch s mplemented smlarly (but takes a bt more work). Eercse 6 Program the sem-lagrangan advecton scheme wth the centered mdpont method as llustrated n fgure 3B. If you don t understand the recpe above check Marc s lecture notes, page 67 for a MATLAB pseudocode. Some care has to be taken f pont X s outsde of the computatonal doman, snce MATLAB wll return NaN for the velocty (or temperature) of ths pont. Use the velocty v (n + 1,) n ths case. A pseudocode s gven by f snan(velocty) Velocty = V() 2D advecton Snce I m convnced that the sem-lagrangan method s the only good advecton algorthm (ecept n the case of pseudospectral methods), ths s the one we re gong to mplement n 2D. Assume that velocty s gven by v (,z) = z (10) v z (,z) = (11) Moreover, assume that the ntal temperature dstrbuton s gaussan and gven by (( ( ) 2 + z 2) ) T (,z) = 2ep (12) wth [ 0.5,0.5], z [ 0.5,0.5]. Eercse 7 Program advecton n 2D usng the sem-lagrangan advecton scheme wth the centered mdpont method. Use the MATLAB routne nterp2 for nterpolaton and employ lnear nterpolaton for velocty and cubc nterpolaton for temperature. A MATLAB scrpt that ll get you started s shown on fgure 4. 5
6 % sem_lagrangan_2d: 2D sem-lagrangan wth center mdpont tme steppng method % clear all W = 40; % wdth of doman sgma =.1; Ampl = 2; nt = 500; % number of tmesteps dt = 5e-1; % Intal grd and velocty [,z] = meshgrd(-0.5:.025:0.5,-0.5:.025:0.5); nz = sze(,1); n = sze(,2); % Intal gaussan T-profle T = Ampl*ep(-((+0.25).ˆ2+z.ˆ2)/sgmaˆ2); % Velocty V = z; Vz = -; for tme=1:nt V_n = V; % Velocty at tme=n V_n1 = V; % Velocty at tme=n+1 % V_n1_2 =??; % Velocty at tme=n+1/2 % Vz_n =??; % Velocty at tme=n % Vz_n1 =??; % Velocty at tme=n+1 % Vz_n1_2 =??; % Velocty at tme=n+1/2 Tnew = zeros(sze(t)); for =2:n-1 for z=2:nz-1 V_cen = V(z,); Vz_cen = Vz(z,); % for?? % X =? % Z =? %lnear nterpolaton of velocty % V_cen = nterp2(,z,?,?,?, lnear ); % Vz_cen = nterp2(,z,?,?,?, lnear ); f snan(v_cen) V_cen = V(z,); f snan(vz_cen) Vz_cen = Vz(z,); % % X =?; % Z =?; % Interpolate temperature on X % T_X = nterp2(,z,?,?,?, cubc ); f snan(t_x) T_X = T(z,); Tnew(z,) = T_X; Tnew(1,:) = T(1,:); Tnew(n,:) = T(n,:); Tnew(:,1) = T(:,1); Tnew(:,n) = T(:,n); T = Tnew; tme = tme*dt; fgure(1),clf pcolor(,z,t), shadng nterp, hold on, colorbar contour(,z,t,[.1:.1:2], k ), hold on, quver(,z,v,vz, w ) as equal, as tght drawnow pause Fgure 4: MATLAB scrpt to be used wth eercse 7. 6
7 Advecton and dffuson: operator splttng Sofar we learned that the sem-lagrangan methods are probably the best methods for advecton domnated problems. We also learned how to solve the dffuson equaton n 1-D and 2-D. In geodynamcs, we often want to solve the coupled advecton-dffuson equaton, whch s gven by equaton 1 n 1-D and by equaton 2 n 2-D. We can solve ths pretty easly by takng the equaton apart and by computng the advecton part separately from the dffuson part. Ths s called operator-splttng, and what s done n 1-D s the followng. Frst solve the advecton equaton T n+1 T n + v = 0 (13) for eample by usng a sem-lagrangan advecton scheme. Then solve the dffuson equaton ρc p T t = ( k ) T + Q (14) and you re done (note that I assumed that Q s spatally constant; f not you should consder to slghtly mprove the advecton scheme by ntroducng source terms). In 2D the scheme s dentcal (but you ll have to fnd your old 2D codes...). Eercse 8 Program dffuson-advecton n 2D usng the sem-lagrangan advecton scheme coupled wth an mplct 2D dffuson code (from last week s eercse). Base your code on the scrpt of fgure 4. 7
2 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 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 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 informationCHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE
CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE Analytcal soluton s usually not possble when exctaton vares arbtrarly wth tme or f the system s nonlnear. Such problems can be solved by numercal tmesteppng
More informationNON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS
IJRRAS 8 (3 September 011 www.arpapress.com/volumes/vol8issue3/ijrras_8_3_08.pdf NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS H.O. Bakodah Dept. of Mathematc
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 informationχ x B E (c) Figure 2.1.1: (a) a material particle in a body, (b) a place in space, (c) a configuration of the body
Secton.. Moton.. The Materal Body and Moton hyscal materals n the real world are modeled usng an abstract mathematcal entty called a body. Ths body conssts of an nfnte number of materal partcles. Shown
More informationLecture Notes on Linear Regression
Lecture Notes on Lnear Regresson Feng L fl@sdueducn Shandong Unversty, Chna Lnear Regresson Problem In regresson problem, we am at predct a contnuous target value gven an nput feature vector We assume
More informationPhysics 5153 Classical Mechanics. Principle of Virtual Work-1
P. Guterrez 1 Introducton Physcs 5153 Classcal Mechancs Prncple of Vrtual Work The frst varatonal prncple we encounter n mechancs s the prncple of vrtual work. It establshes the equlbrum condton of a mechancal
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 informationLine Drawing and Clipping Week 1, Lecture 2
CS 43 Computer Graphcs I Lne Drawng and Clppng Week, Lecture 2 Davd Breen, Wllam Regl and Maxm Peysakhov Geometrc and Intellgent Computng Laboratory Department of Computer Scence Drexel Unversty http://gcl.mcs.drexel.edu
More informationFTCS Solution to the Heat Equation
FTCS Soluton to the Heat Equaton ME 448/548 Notes Gerald Recktenwald Portland State Unversty Department of Mechancal Engneerng gerry@pdx.edu ME 448/548: FTCS Soluton to the Heat Equaton Overvew 1. Use
More informationFinite Element Modelling of truss/cable structures
Pet Schreurs Endhoven Unversty of echnology Department of Mechancal Engneerng Materals echnology November 3, 214 Fnte Element Modellng of truss/cable structures 1 Fnte Element Analyss of prestressed structures
More informationE91: Dynamics. E91: Dynamics. Numerical Integration & State Space Representation
E91: Dnamcs Numercal Integraton & State Space Representaton The Algorthm In steps of Δ f ( new ) f ( old ) df ( d old ) Δ Numercal Integraton of ODEs d d f() h Δ Intal value problem: Gven the ntal state
More informationLecture 5.8 Flux Vector Splitting
Lecture 5.8 Flux Vector Splttng 1 Flux Vector Splttng The vector E n (5.7.) can be rewrtten as E = AU (5.8.1) (wth A as gven n (5.7.4) or (5.7.6) ) whenever, the equaton of state s of the separable form
More informationDefinition. Measures of Dispersion. Measures of Dispersion. Definition. The Range. Measures of Dispersion 3/24/2014
Measures of Dsperson Defenton Range Interquartle Range Varance and Standard Devaton Defnton Measures of dsperson are descrptve statstcs that descrbe how smlar a set of scores are to each other The more
More informationLinear Feature Engineering 11
Lnear Feature Engneerng 11 2 Least-Squares 2.1 Smple least-squares Consder the followng dataset. We have a bunch of nputs x and correspondng outputs y. The partcular values n ths dataset are x y 0.23 0.19
More informationNumerical Transient Heat Conduction Experiment
Numercal ransent Heat Conducton Experment OBJECIVE 1. o demonstrate the basc prncples of conducton heat transfer.. o show how the thermal conductvty of a sold can be measured. 3. o demonstrate the use
More informationChapter - 2. Distribution System Power Flow Analysis
Chapter - 2 Dstrbuton System Power Flow Analyss CHAPTER - 2 Radal Dstrbuton System Load Flow 2.1 Introducton Load flow s an mportant tool [66] for analyzng electrcal power system network performance. Load
More information36.1 Why is it important to be able to find roots to systems of equations? Up to this point, we have discussed how to find the solution to
ChE Lecture Notes - D. Keer, 5/9/98 Lecture 6,7,8 - Rootndng n systems o equatons (A) Theory (B) Problems (C) MATLAB Applcatons Tet: Supplementary notes rom Instructor 6. Why s t mportant to be able to
More informationNumerical Simulation of Wave Propagation Using the Shallow Water Equations
umercal Smulaton of Wave Propagaton Usng the Shallow Water Equatons Junbo Par Harve udd College 6th Aprl 007 Abstract The shallow water equatons SWE were used to model water wave propagaton n one dmenson
More informationConsideration of 2D Unsteady Boundary Layer Over Oscillating Flat Plate
Proceedngs of the th WSEAS Internatonal Conference on Flud Mechancs and Aerodynamcs, Elounda, Greece, August -, (pp-) Consderaton of D Unsteady Boundary Layer Over Oscllatng Flat Plate N.M. NOURI, H.R.
More informationImplicit Integration Henyey Method
Implct Integraton Henyey Method In realstc stellar evoluton codes nstead of a drect ntegraton usng for example the Runge-Kutta method one employs an teratve mplct technque. Ths s because the structure
More informationLab session: numerical simulations of sponateous polarization
Lab sesson: numercal smulatons of sponateous polarzaton Emerc Boun & Vncent Calvez CNRS, ENS Lyon, France CIMPA, Hammamet, March 2012 Spontaneous cell polarzaton: the 1D case The Hawkns-Voturez model for
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 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 informationMATH 5630: Discrete Time-Space Model Hung Phan, UMass Lowell March 1, 2018
MATH 5630: Dscrete Tme-Space Model Hung Phan, UMass Lowell March, 08 Newton s Law of Coolng Consder the coolng of a well strred coffee so that the temperature does not depend on space Newton s law of collng
More information6.3.7 Example with Runga Kutta 4 th order method
6.3.7 Example wth Runga Kutta 4 th order method Agan, as an example, 3 machne, 9 bus system shown n Fg. 6.4 s agan consdered. Intally, the dampng of the generators are neglected (.e. d = 0 for = 1, 2,
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 informationLecture 13 APPROXIMATION OF SECOMD ORDER DERIVATIVES
COMPUTATIONAL FLUID DYNAMICS: FDM: Appromaton of Second Order Dervatves Lecture APPROXIMATION OF SECOMD ORDER DERIVATIVES. APPROXIMATION OF SECOND ORDER DERIVATIVES Second order dervatves appear n dffusve
More informationMechanics Physics 151
Mechancs Physcs 5 Lecture 7 Specal Relatvty (Chapter 7) What We Dd Last Tme Worked on relatvstc knematcs Essental tool for epermental physcs Basc technques are easy: Defne all 4 vectors Calculate c-o-m
More informationDifference Equations
Dfference Equatons c Jan Vrbk 1 Bascs Suppose a sequence of numbers, say a 0,a 1,a,a 3,... s defned by a certan general relatonshp between, say, three consecutve values of the sequence, e.g. a + +3a +1
More informationChapter 3 Differentiation and Integration
MEE07 Computer Modelng Technques n Engneerng Chapter Derentaton and Integraton Reerence: An Introducton to Numercal Computatons, nd edton, S. yakowtz and F. zdarovsky, Mawell/Macmllan, 990. Derentaton
More informationErrors for Linear Systems
Errors for Lnear Systems When we solve a lnear system Ax b we often do not know A and b exactly, but have only approxmatons  and ˆb avalable. Then the best thng we can do s to solve ˆx ˆb exactly whch
More informationA Hybrid Variational Iteration Method for Blasius Equation
Avalable at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Vol. 10, Issue 1 (June 2015), pp. 223-229 Applcatons and Appled Mathematcs: An Internatonal Journal (AAM) A Hybrd Varatonal Iteraton Method
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 informationModule 9. Lecture 6. Duality in Assignment Problems
Module 9 1 Lecture 6 Dualty n Assgnment Problems In ths lecture we attempt to answer few other mportant questons posed n earler lecture for (AP) and see how some of them can be explaned through the concept
More informationModule 1 : The equation of continuity. Lecture 1: Equation of Continuity
1 Module 1 : The equaton of contnuty Lecture 1: Equaton of Contnuty 2 Advanced Heat and Mass Transfer: Modules 1. THE EQUATION OF CONTINUITY : Lectures 1-6 () () () (v) (v) Overall Mass Balance Momentum
More information2.29 Numerical Fluid Mechanics Fall 2011 Lecture 12
REVIEW Lecture 11: 2.29 Numercal Flud Mechancs Fall 2011 Lecture 12 End of (Lnear) Algebrac Systems Gradent Methods Krylov Subspace Methods Precondtonng of Ax=b FINITE DIFFERENCES Classfcaton of Partal
More information1-Dimensional Advection-Diffusion Finite Difference Model Due to a Flow under Propagating Solitary Wave
014 4th Internatonal Conference on Future nvronment and nergy IPCB vol.61 (014) (014) IACSIT Press, Sngapore I: 10.776/IPCB. 014. V61. 6 1-mensonal Advecton-ffuson Fnte fference Model ue to a Flow under
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 informationThe Geometry of Logit and Probit
The Geometry of Logt and Probt Ths short note s meant as a supplement to Chapters and 3 of Spatal Models of Parlamentary Votng and the notaton and reference to fgures n the text below s to those two chapters.
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 information12. The Hamilton-Jacobi Equation Michael Fowler
1. The Hamlton-Jacob Equaton Mchael Fowler Back to Confguraton Space We ve establshed that the acton, regarded as a functon of ts coordnate endponts and tme, satsfes ( ) ( ) S q, t / t+ H qpt,, = 0, and
More informationPhysics 181. Particle Systems
Physcs 181 Partcle Systems Overvew In these notes we dscuss the varables approprate to the descrpton of systems of partcles, ther defntons, ther relatons, and ther conservatons laws. We consder a system
More informationOn the Interval Zoro Symmetric Single-step Procedure for Simultaneous Finding of Polynomial Zeros
Appled Mathematcal Scences, Vol. 5, 2011, no. 75, 3693-3706 On the Interval Zoro Symmetrc Sngle-step Procedure for Smultaneous Fndng of Polynomal Zeros S. F. M. Rusl, M. Mons, M. A. Hassan and W. J. Leong
More informationAdditional Codes using Finite Difference Method. 1 HJB Equation for Consumption-Saving Problem Without Uncertainty
Addtonal Codes usng Fnte Dfference Method Benamn Moll 1 HJB Equaton for Consumpton-Savng Problem Wthout Uncertanty Before consderng the case wth stochastc ncome n http://www.prnceton.edu/~moll/ HACTproect/HACT_Numercal_Appendx.pdf,
More informationChapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems
Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons
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 2 Solution of Nonlinear Equations ( Root Finding Problems )
Lecture Soluton o Nonlnear Equatons Root Fndng Problems Dentons Classcaton o Methods Analytcal Solutons Graphcal Methods Numercal Methods Bracketng Methods Open Methods Convergence Notatons Root Fndng
More informationCIS526: Machine Learning Lecture 3 (Sept 16, 2003) Linear Regression. Preparation help: Xiaoying Huang. x 1 θ 1 output... θ M x M
CIS56: achne Learnng Lecture 3 (Sept 6, 003) Preparaton help: Xaoyng Huang Lnear Regresson Lnear regresson can be represented by a functonal form: f(; θ) = θ 0 0 +θ + + θ = θ = 0 ote: 0 s a dummy attrbute
More information(Online First)A Lattice Boltzmann Scheme for Diffusion Equation in Spherical Coordinate
Internatonal Journal of Mathematcs and Systems Scence (018) Volume 1 do:10.494/jmss.v1.815 (Onlne Frst)A Lattce Boltzmann Scheme for Dffuson Equaton n Sphercal Coordnate Debabrata Datta 1 *, T K Pal 1
More informationNote 10. Modeling and Simulation of Dynamic Systems
Lecture Notes of ME 475: Introducton to Mechatroncs Note 0 Modelng and Smulaton of Dynamc Systems Department of Mechancal Engneerng, Unversty Of Saskatchewan, 57 Campus Drve, Saskatoon, SK S7N 5A9, Canada
More informationLecture 2: Numerical Methods for Differentiations and Integrations
Numercal Smulaton of Space Plasmas (I [AP-4036] Lecture 2 by Lng-Hsao Lyu March, 2018 Lecture 2: Numercal Methods for Dfferentatons and Integratons As we have dscussed n Lecture 1 that numercal smulaton
More informationModified Mass Matrices and Positivity Preservation for Hyperbolic and Parabolic PDEs
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2000; 00:6 Prepared usng cnmauth.cls [Verson: 2000/03/22 v.0] Modfed Mass Matrces and Postvty Preservaton for Hyperbolc and
More informationUncertainty as the Overlap of Alternate Conditional Distributions
Uncertanty as the Overlap of Alternate Condtonal Dstrbutons Olena Babak and Clayton V. Deutsch Centre for Computatonal Geostatstcs Department of Cvl & Envronmental Engneerng Unversty of Alberta An mportant
More informationEcon107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4)
I. Classcal Assumptons Econ7 Appled Econometrcs Topc 3: Classcal Model (Studenmund, Chapter 4) We have defned OLS and studed some algebrac propertes of OLS. In ths topc we wll study statstcal propertes
More informationCHARACTERISTICS OF COMPLEX SEPARATION SCHEMES AND AN ERROR OF SEPARATION PRODUCTS OUTPUT DETERMINATION
Górnctwo Geonżynera Rok 0 Zeszyt / 006 Igor Konstantnovch Mladetskj * Petr Ivanovch Plov * Ekaterna Nkolaevna Kobets * Tasya Igorevna Markova * CHARACTERISTICS OF COMPLEX SEPARATION SCHEMES AND AN ERROR
More informationChapter 4 The Wave Equation
Chapter 4 The Wave Equaton Another classcal example of a hyperbolc PDE s a wave equaton. The wave equaton s a second-order lnear hyperbolc PDE that descrbes the propagaton of a varety of waves, such as
More informationHigher Order Wall Boundary Conditions for Incompressible Flow Simulations
THE 5 TH ASIAN COMPUTAITIONAL FLUID DYNAMICS BUSAN KOREA OCTOBER 7-30 003 Hgher Order Wall Boundary Condtons for Incompressble Flow Smulatons Hdetosh Nshda. Department of Mechancal and System Engneerng
More informationTensor Smooth Length for SPH Modelling of High Speed Impact
Tensor Smooth Length for SPH Modellng of Hgh Speed Impact Roman Cherepanov and Alexander Gerasmov Insttute of Appled mathematcs and mechancs, Tomsk State Unversty 634050, Lenna av. 36, Tomsk, Russa RCherepanov82@gmal.com,Ger@npmm.tsu.ru
More informationA new Approach for Solving Linear Ordinary Differential Equations
, ISSN 974-57X (Onlne), ISSN 974-5718 (Prnt), Vol. ; Issue No. 1; Year 14, Copyrght 13-14 by CESER PUBLICATIONS A new Approach for Solvng Lnear Ordnary Dfferental Equatons Fawz Abdelwahd Department of
More informationComputational Fluid Dynamics. Smoothed Particle Hydrodynamics. Simulations. Smoothing Kernels and Basis of SPH
Computatonal Flud Dynamcs If you want to learn a bt more of the math behnd flud dynamcs, read my prevous post about the Naver- Stokes equatons and Newtonan fluds. The equatons derved n the post are the
More informationDUE: WEDS FEB 21ST 2018
HOMEWORK # 1: FINITE DIFFERENCES IN ONE DIMENSION DUE: WEDS FEB 21ST 2018 1. Theory Beam bendng s a classcal engneerng analyss. The tradtonal soluton technque makes smplfyng assumptons such as a constant
More informationKernel Methods and SVMs Extension
Kernel Methods and SVMs Extenson The purpose of ths document s to revew materal covered n Machne Learnng 1 Supervsed Learnng regardng support vector machnes (SVMs). Ths document also provdes a general
More informationNice plotting of proteins II
Nce plottng of protens II Fnal remark regardng effcency: It s possble to wrte the Newton representaton n a way that can be computed effcently, usng smlar bracketng that we made for the frst representaton
More informationPhysics 207: Lecture 20. Today s Agenda Homework for Monday
Physcs 207: Lecture 20 Today s Agenda Homework for Monday Recap: Systems of Partcles Center of mass Velocty and acceleraton of the center of mass Dynamcs of the center of mass Lnear Momentum Example problems
More informationGlobal Sensitivity. Tuesday 20 th February, 2018
Global Senstvty Tuesday 2 th February, 28 ) Local Senstvty Most senstvty analyses [] are based on local estmates of senstvty, typcally by expandng the response n a Taylor seres about some specfc values
More informationOne-sided finite-difference approximations suitable for use with Richardson extrapolation
Journal of Computatonal Physcs 219 (2006) 13 20 Short note One-sded fnte-dfference approxmatons sutable for use wth Rchardson extrapolaton Kumar Rahul, S.N. Bhattacharyya * Department of Mechancal Engneerng,
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 informationKinematics of Fluids. Lecture 16. (Refer the text book CONTINUUM MECHANICS by GEORGE E. MASE, Schaum s Outlines) 17/02/2017
17/0/017 Lecture 16 (Refer the text boo CONTINUUM MECHANICS by GEORGE E. MASE, Schaum s Outlnes) Knematcs of Fluds Last class, we started dscussng about the nematcs of fluds. Recall the Lagrangan and Euleran
More informationFor now, let us focus on a specific model of neurons. These are simplified from reality but can achieve remarkable results.
Neural Networks : Dervaton compled by Alvn Wan from Professor Jtendra Malk s lecture Ths type of computaton s called deep learnng and s the most popular method for many problems, such as computer vson
More information: Numerical Analysis Topic 2: Solution of Nonlinear Equations Lectures 5-11:
764: Numercal Analyss Topc : Soluton o Nonlnear Equatons Lectures 5-: UIN Malang Read Chapters 5 and 6 o the tetbook 764_Topc Lecture 5 Soluton o Nonlnear Equatons Root Fndng Problems Dentons Classcaton
More informationCurve Fitting with the Least Square Method
WIKI Document Number 5 Interpolaton wth Least Squares Curve Fttng wth the Least Square Method Mattheu Bultelle Department of Bo-Engneerng Imperal College, London Context We wsh to model the postve feedback
More informationGaussian Mixture Models
Lab Gaussan Mxture Models Lab Objectve: Understand the formulaton of Gaussan Mxture Models (GMMs) and how to estmate GMM parameters. You ve already seen GMMs as the observaton dstrbuton n certan contnuous
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 informationStanford University CS359G: Graph Partitioning and Expanders Handout 4 Luca Trevisan January 13, 2011
Stanford Unversty CS359G: Graph Parttonng and Expanders Handout 4 Luca Trevsan January 3, 0 Lecture 4 In whch we prove the dffcult drecton of Cheeger s nequalty. As n the past lectures, consder an undrected
More informationMean Field / Variational Approximations
Mean Feld / Varatonal Appromatons resented by Jose Nuñez 0/24/05 Outlne Introducton Mean Feld Appromaton Structured Mean Feld Weghted Mean Feld Varatonal Methods Introducton roblem: We have dstrbuton but
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 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 informationHomework Assignment 3 Due in class, Thursday October 15
Homework Assgnment 3 Due n class, Thursday October 15 SDS 383C Statstcal Modelng I 1 Rdge regresson and Lasso 1. Get the Prostrate cancer data from http://statweb.stanford.edu/~tbs/elemstatlearn/ datasets/prostate.data.
More informationWeek 11: Chapter 11. The Vector Product. The Vector Product Defined. The Vector Product and Torque. More About the Vector Product
The Vector Product Week 11: Chapter 11 Angular Momentum There are nstances where the product of two vectors s another vector Earler we saw where the product of two vectors was a scalar Ths was called the
More informationCase A. P k = Ni ( 2L i k 1 ) + (# big cells) 10d 2 P k.
THE CELLULAR METHOD In ths lecture, we ntroduce the cellular method as an approach to ncdence geometry theorems lke the Szemeréd-Trotter theorem. The method was ntroduced n the paper Combnatoral complexty
More informationMA 323 Geometric Modelling Course Notes: Day 13 Bezier Curves & Bernstein Polynomials
MA 323 Geometrc Modellng Course Notes: Day 13 Bezer Curves & Bernsten Polynomals Davd L. Fnn Over the past few days, we have looked at de Casteljau s algorthm for generatng a polynomal curve, and we have
More informationThis column is a continuation of our previous column
Comparson of Goodness of Ft Statstcs for Lnear Regresson, Part II The authors contnue ther dscusson of the correlaton coeffcent n developng a calbraton for quanttatve analyss. Jerome Workman Jr. and Howard
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 informationCHAPTER 14 GENERAL PERTURBATION THEORY
CHAPTER 4 GENERAL PERTURBATION THEORY 4 Introducton A partcle n orbt around a pont mass or a sphercally symmetrc mass dstrbuton s movng n a gravtatonal potental of the form GM / r In ths potental t moves
More informationBasically, if you have a dummy dependent variable you will be estimating a probability.
ECON 497: Lecture Notes 13 Page 1 of 1 Metropoltan State Unversty ECON 497: Research and Forecastng Lecture Notes 13 Dummy Dependent Varable Technques Studenmund Chapter 13 Bascally, f you have a dummy
More informationFundamental loop-current method using virtual voltage sources technique for special cases
Fundamental loop-current method usng vrtual voltage sources technque for specal cases George E. Chatzaraks, 1 Marna D. Tortorel 1 and Anastasos D. Tzolas 1 Electrcal and Electroncs Engneerng Departments,
More informationA New Refinement of Jacobi Method for Solution of Linear System Equations AX=b
Int J Contemp Math Scences, Vol 3, 28, no 17, 819-827 A New Refnement of Jacob Method for Soluton of Lnear System Equatons AX=b F Naem Dafchah Department of Mathematcs, Faculty of Scences Unversty of Gulan,
More information18.1 Introduction and Recap
CS787: Advanced Algorthms Scrbe: Pryananda Shenoy and Shjn Kong Lecturer: Shuch Chawla Topc: Streamng Algorthmscontnued) Date: 0/26/2007 We contnue talng about streamng algorthms n ths lecture, ncludng
More informationTHE VIBRATIONS OF MOLECULES II THE CARBON DIOXIDE MOLECULE Student Instructions
THE VIBRATIONS OF MOLECULES II THE CARBON DIOXIDE MOLECULE Student Instructons by George Hardgrove Chemstry Department St. Olaf College Northfeld, MN 55057 hardgrov@lars.acc.stolaf.edu Copyrght George
More informationThe Finite Element Method
The Fnte Element Method GENERAL INTRODUCTION Read: Chapters 1 and 2 CONTENTS Engneerng and analyss Smulaton of a physcal process Examples mathematcal model development Approxmate solutons and methods of
More informationIntroduction to Vapor/Liquid Equilibrium, part 2. Raoult s Law:
CE304, Sprng 2004 Lecture 4 Introducton to Vapor/Lqud Equlbrum, part 2 Raoult s Law: The smplest model that allows us do VLE calculatons s obtaned when we assume that the vapor phase s an deal gas, and
More informationIntroduction to Computational Fluid Dynamics
Introducton to Computatonal Flud Dynamcs M. Zanub 1, T. Mahalakshm 2 1 (PG MATHS), Department of Mathematcs, St. Josephs College of Arts and Scence for Women-Hosur, Peryar Unversty 2 Assstance professor,
More informationExtension of Smoothed Particle Hydrodynamics (SPH), Mathematical Background of Vortex Blob Method (VBM) and Moving Particle Semi-Implicit (MPS)
Amercan Journal of Computatonal athematcs, 04, 5, 44-445 Publshed Onlne December 04 n ScRes. http://www.scrp.org/ournal/acm http://dx.do.org/0.436/acm.04.45036 Extenson of Smoothed Partcle Hydrodynamcs
More informationChapter Twelve. Integration. We now turn our attention to the idea of an integral in dimensions higher than one. Consider a real-valued function f : D
Chapter Twelve Integraton 12.1 Introducton We now turn our attenton to the dea of an ntegral n dmensons hgher than one. Consder a real-valued functon f : R, where the doman s a nce closed subset of Eucldean
More informationReview of Taylor Series. Read Section 1.2
Revew of Taylor Seres Read Secton 1.2 1 Power Seres A power seres about c s an nfnte seres of the form k = 0 k a ( x c) = a + a ( x c) + a ( x c) + a ( x c) k 2 3 0 1 2 3 + In many cases, c = 0, and the
More informationLecture 20: Noether s Theorem
Lecture 20: Noether s Theorem In our revew of Newtonan Mechancs, we were remnded that some quanttes (energy, lnear momentum, and angular momentum) are conserved That s, they are constant f no external
More informationBindel, Spring 2014 Applications of Parallel Computers (CS 5220)
Bndel, Sprng 2014 Applcatons of Parallel Computers (CS 5220) 1 Dervng SPH The Naver-Stokes equatons wth gravty are ρa = p + µ 2 v + ρg. The acceleraton s the materal dervatve of velocty, and we usually
More information