FUNDAMENTALS OF FINITE DIFFERENCE METHODS
|
|
- Paul Briggs
- 6 years ago
- Views:
Transcription
1 FUNDAMENTALS OF FINITE DIFFERENCE METHODS By, Varn Khatan 3 rd year Undergradate IIT Kanpr Spervsed by, Professor Gatam Bswas, Mechancal Engneerng IIT Kanpr
2 We wll dscss. Classfcaton of Partal Dfferental eqatons nvolved n fld mechancs Fndamentals of Dscretzaton. Fnte Dfference Qotents, startng from Taylor s seres. Fnte Dfference eqatons from the pont of vew of ther consstency and convergence. Explct and Implct methods of dscretzaton. Errors and stablty analyss. Fndamentals of Fld Flow Modelng: Transportve and Conservatve propertes. Artfcal Vscosty. Upwnd Scheme of Dscretzaton (Frst and Second).
3 Partal Dfferental Eqatons Consder a second order dfferental eqn φ φ φ φ φ φ =, ( ) A B C D E F G x y x x y y x y A,B,C,D,E and F are constants or fnctons of only x and y, then lnear eqaton Qaslnear eqaton Homogenos eqaton Classfcaton: 4 = 0, B AC parabolc B 4AC < 0, ellptc B 4AC > 0, hyperbolc
4 Partal Dfferental Eqatons (examples) Laplace eqaton: Posson eqaton: o φ x φ x φ y + = 0 φ y + + S = 0 Both are ellptc eqatons Heat condcton eqaton s parabolc. φ φ = B t x
5 Contned Fld flow problems have nonlnear terms de to nerta and energy terms known as advecton and convecton respectvely. φ φ φ φ φ + + v = B ( + ) + S t x y x y φ denotes some transported property lke velocty, temperatre etc. and v denote velocty components. B s dffsvty for momentm or heat. S s sorce term. Parabolc n tme and ellptc n space. For very hgh speeds t becomes hyperbolc.
6 Bondary and Intal condtons No. of B.C.s=order of hghest dervatves n dfferental eqatons. Unsteady problems reqre ntal condtons for tme ntegraton. Dffson terms need two BC s. Spatal Bondary Condtons: Drchlet: Φ=Φ (r) A Nemann: Φ n=φ (r) A Mxed:a(r)Φ+b(r) Φ n=φ 3 (r) A 3
7 Fnte Dfference Methods Dscretzaton s a process by whch a partal dfferental eqaton s approxmated by analogos expressons whch prescrbe vales at only a fnte no. of dscrete ponts or volmes n doman. Basc phlosophy of fnte dfference methods s to replace dervatves wth algebrac dfference qotents, resltng n a system of algebrac eqatons whch can then be solved.
8 Grd Generaton x and y need not be same. When they are dfferent n the physcal plane, the plane s transformed. In transformed comptatonal plane the spacng s nform n all drectons.
9 Fnte dfference qotent Taylor seres expanson Frst order Forward dfference
10 Fnte dfference qotent (Contd )
11 Fnte dfference qotent (Contd ) Second order central dfference If mxed dervatves are nvolved the second order dfference looks lke xy φ + +, j+, j +, j, j+, j φ φ φ φ ( ) ( ) = +Ο x, y 4 x y Δ Δ ΔΔ
12 Fnte Dfference Eqatons The basc concept n solvng the PDE sng ths approach s to replace each dervatve term wth ts dfference element. Consder the nsteady -D heat condcton: t = α x Usng the concept of FTCS (Forward Tme Central Space) We transform the above eqaton to n+ Δt n ( n n n ) + + α + TE Δx ( )
13 Consstency and Convergence A fnte dfference representaton s sad to be consstent f ( PDE FDE ) lm = lm ( TE ) = 0 mesh 0 mesh 0 o A necessary condton s that Δ t Δ x 0 As x ->0 and t->0 o A FDE s sad to be convergent f the approxmate solton approaches the exact one as the sze of the mesh tends to zero. n = x, t asδx, Δt 0 ( ) n
14 Explct and mplct methods Explct method: Consder -D nsteady state heat condcton eqaton. o n+ Δt t n = α x The correspondng Fnte Dfference eqaton s: = α n n n ( + + ) ( Δx ) n + The only nknown n the above eqaton s Its vale can be obtaned from the vales of the varables at tme n. Ths the eqaton takes the form of a marchng procedre n steps of tme.
15 Explct and mplct methods (contd ) Implct method: Consderng the earler eqaton wth a dfferent dscretzaton algorthm we get the form n+ Δt n = α ( n+ n n+ n n+ n ) ( Δx) Here each term has been replaced by the tme average of ts neghborng terms. Ths s known as the Crank-Ncholson Implct Scheme.
16 Explct and mplct methods (contd ) Implct Method (cont.): The eqaton can be wrtten as, n+ n r ( n+ n n+ n n+ n ) = α where r= ( Δt) ( Δx) or, + r r + = + + r r r n+ n+ n+ n n n + + These eqatons when appled at all grd pts. from = to =k+ reslt n a trdagonal matrx whch can be solved to obtan the vales at tme n+. n+ n B() 0 0 K 0 T ( C() + A) n n B() K T3 C() n+ n 0 B(3) K 0 T4 = C(3) M M M M M M M n+ n K Bk ( ) T k ( Ck ( ) + D) Here A and D come from the bondary condtons and C(k) s the vale of the RHS at tme k+.
17 Explct and mplct methods (contd ) o o o Applyng the Crank Ncholson dscretzaton scheme to the -D heat eqaton gven by = α + t x y where we get δ δ Ths system has fve nknowns n+ n, j, j n n Δt α = δ + δ + + ( )( + ) x y, j, j n n n n +, j, j, j x T, j = ( Δ x ) + n n n n, j +, j, j y T, j = ( Δ y ),,, & n+ n+ n+ n+ n+, j +, j, j, j+, j On smplfcaton ths gves a matrx eqaton whch s mch more complcated to solve and consmes more tme.
18 Explct and mplct methods (contd ) Usng alternate drecton mplct method (e) splttng the eqaton nto two steps we obtan the followng new set of eqatons: Frst step: Second step: + / + / + / ( ) ( ) n+ / n, j + +, j = α + Δt / ( Δx) ( Δy) n n n n n n +, j, j, j, j+, j, j ( ) ( ) n+ n+ /, j + +, j = α + Δt / ( Δx) ( Δy) n+ / n+ / n+ / n+ n+ n+ +, j, j, j, j+, j, j Ths system gves two trdagonal matrx eqatons, one for each step. These can be solved to get the solton.
19 Explct and mplct methods (contd ) Implct Method Advantage: Takes lesser tme for calclatons over a tme nterval. Dsadvantage: - Complcated program. -Large comptaton tme n each step -For larger Δt, trncaton error s hgh. Explct method Advantage: The algorthm s smpler Dsadvantage: Reqres many tme steps to carry ot the calclatons over a gven tme nterval, de to restrctons on Δt mposed by stablty constrants.
20 Errors and stablty analyss Errors: A=analytcal sol of P.D.E. D=exact sol. of F.D.E. N=nmercal sol. from a real compter wth fnte accracy. Then Dscretzaton Error =A-D DE=Trncaton Error+ error ntrodced de to treatment of bondary condton. Rond-off Error = ε =N-D Nmercal error ntrodced by the compter as t ronds off the nmbers to some fxed decmal ponts. Ths, N=D+ε
21 Errors and stablty analyss (contd..) Consder the FDE of the one dmensonal heat eqaton: n + n n n n ( ) α + + = Δ t Δ x The Nmercal solton (N) and the exact solton (D) both satsfy the above eqaton. Owng to the lnearty of the eqns ε=n-d wll also satsfy t. Hence we get ε n+ ε Δt n = ( n n n ) + + α ε ε ε ( Δx ) We say the solton s stable f the error shrnks or remans same as we progress n tme. If the error magnfes wth tme then the solton s nstable. Ths the condton for stablty s: ( ) ε ε n + n
22 Errors and stablty analyss (contd..) Expandng the error expresson as a Forer seres, we get Snce the FDE s lnear t wold sffce to consder the term Sbstttng ths n the FDE we get: ε n + n ε 4αΔt e = sn k Δx/ Δx aδt aδt = e where ( ) m
23 Errors and stablty analyss (contd..) Ths we reqre: 4αΔt Δx ( k x ) sn / mδ Ths gves αδt Δx αδt Δx sn ( k / ) mδx / / Ths s called Von Nemann stablty analyss.
24 Stablty of hyperbolc eqatons Applyng Von Nemann stablty analyss to the frst order eqaton = t t + c = x We replace tme dervatve wth frst order dfference and spatal dervatve by central dfference to. Also, we replace (t) by the average of neghborng terms. Ths + Δt ( ) + n+ n n 0 ( n n ) n n n+ + + Δt + = c Δx e Δ = cos k Δx Csn k Δx The amplfcaton factor s at ( ) ( ) m m
25 Stablty of hyperbolc eqatons Ths the reqred condton for stablty s: Δt C = c Δx Here C s called the Corant nmber and the condton s CFL (Corant - Fredrchs - Lewy) condton. Applyng the same procedre to the second order eqaton also gves the same condton for stablty. In actal problems stablty s checked locally. For varable coeffcents Von Nemann analyss s necessary bt not sffcent. To assre stablty we need to carry ot expermentaton and nmercal checks. Rotne stablty analyss only gves gdance.
26 Fld Flow modellng Fld flow eqatons are complex as they form a nonlnear system. To model fld flow we consder the Brger s eqaton gven by: ζ ζ + = υ ζ t x x Here s the velocty, v s the vscosty and T s any property whch can be transported and dffsed. If vscos term on the rght sde s neglected the eqaton redces to: ζ ζ + = 0 t x An FDE s sad to be conservatve f t preserves the ntegral conservatve relatons. ω t = V. ω + υ ω Consder the Vortcty transport eq. ( ) Integratng over fxed regon, fnally get,
27 Fld Flow modellng (contd..) ( ) ( ) ω d R= V ω. nda + υ ω. nda t R A0 A0 Ths mples that the tme rate of accmlaton of ω n R eqals net advectve flx rate of ω across A0 nto Rand net dffsve flx rate of ω across A0 nto. Ths relaton needs to be preserved for conservatness. R Consder the conservatve form of the Brger eqaton: ω = ( ω ) t x Usng FTCS the correspondng FDE s ω ω ω ω = Δt Δx n+ n n n n n + + Evalatng the ntegral over a regon between =I R and =I as = I ω Δ x Δ t = I I I I we get n+ n n n ω Δx ω Δ x = ( ω) ( ω) + Δt = I = I = I
28 Fld Flow modellng (contd..) Smmaton of R.H.S. fnally gves, Integral ( ω) ( ω) = I I + Ths the accmalaton of w n s eqal to the net advectve flx rate across the bondary of R. Ths the FDE preserves the ntegral relaton and s conservatve n natre. If we consder the non-conservatve form of the Brger s eqaton ω ω =, t x then followng the earler procedre the ntegral comes ot to be I n n n n = I ω + Here t fals to preserve the Gass-Dvergence property and hence t s not conservatve. R ω
29 Upwnd scheme and Transportve property Von Nemann s stablty analyss appled to normal dscretzaton of the Brger s eqaton ndcates ncondtonal nstablty. Ths we apply an alternate scheme to the convectve term n the conservatve form: n+ n n n ζ ζ ζ ζ = + Δt Δx n+ n n n ζ ζ ζ + ζ = + Δt Δx vscos terms, >0 vscos terms, <0 Ths s known as the pwnd scheme. Transportve property: An FDE of a flow eqaton s sad to possess ths property f the effect of a pertrbaton s convected or advected only n the drecton of the velocty. Reglar formlaton of the FDE volates ths property
30 Transportve property (contd.) Consder a pertrbaton ε m =δ n ζ at the locaton m at tme=n. For >0, at the downstream locaton m+ we get o o n+ n ζ m+ ζ m+ 0 δ δ = =+ Δt Δx Δx Ths s acceptable. At the pont of pertrbaton,(e) m we get ζ n + n m ζ m δ 0 = = δ Δt Δx Δx Ths ndcates that the pertrbaton s beng transported ot of the affected regon. At the pstream staton, m- we see that ζ n + m n ζ m 0 0 = = Δ t Δ x Ths flow s ndrectonal and property s satsfed. 0
31 Upwnd scheme and artfcal vscosty Usng Taylor seres expanson of the terms n the Upwnd Scheme eqaton, we get or ( Δt) Ths eqaton may be wrtten as ( Δx) ζ ζ 3 ζ ζ 3 Δ t + +Ο ( Δ t) = Δx +Ο ( Δx) +[dffsve terms] Δt t t Δx x x ζ ζ Δt ζ ζ = + Δx + ν +Ο Δx t x x Δ x x ζ = ζ + ν ζ e + ν ζ + hgher-order terms t x x x ( ) Where ν ( ) s the artfcal vscosty and C e = Δx C 3 s the Corant nmber. In dervng ths we have taken ζ t ζ x =
32 Upwnd scheme and artfcal vscosty In steady state,.e., ζ/ t=0,νe= x/. For a two dmensonal convectve-dffsve eqaton, on followng smlar steps the eqaton obtaned s ζ ζ v ζ ( ) ζ = + ν ex + ν + ( ν ) ey + ν ζ t x y x y Where ν ex = ( ) wth C x = t/ x, Δx Cx ν ( ) ey = Δy C y C y =v t/ x
33 Hgher order pwnd dfferencng In order to mantan the transportve property the comptatons are naccrate de to the presence of false dffson. To get rd of ths we se hgher order dfferencng. Accordng to second order dfferencng ( ζ ) Rζ R Lζ L = x Δx, j where, R=(,j++,j)/ ; L=(,j+-,j)/ ; and, ζr=ζ,j for R>0; ζr=ζ+,j for R<0; ζl=ζ-,j for L>0; ζl=ζ,j for L<0; For R>0 and L>0 we get ( ζ ), +,,, j+ j j+ j = ζ, j ζ, j x Δx
34 x Hybrd scheme We can also represent t as a weghted combnaton of central and pwnd dfferencng, (, j +, j) η (, j +, j) = 4Δ + + ( ) ( ), j+, j η, j, j, j x Where 0<η<. The accracy of the scheme can always be ncreased by a stable adjstment of η vale.
35 Acknowledgements I take ths opportnty to thank my spervsor Prof. Gatam Bswas (Mechancal Engneerng,IIT Kanpr) for hs gdance and enthsastc spport, whch has played a pvotal role n the preparaton of ths talk.
AE/ME 339. K. M. Isaac. 8/31/2004 topic4: Implicit method, Stability, ADI method. Computational Fluid Dynamics (AE/ME 339) MAEEM Dept.
AE/ME 339 Comptatonal Fld Dynamcs (CFD) Comptatonal Fld Dynamcs (AE/ME 339) Implct form of dfference eqaton In the prevos explct method, the solton at tme level n,,n, depended only on the known vales of,
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 information3/31/ = 0. φi φi. Use 10 Linear elements to solve the equation. dx x dx THE PETROV-GALERKIN METHOD
THE PETROV-GAERKIN METHO Consder the Galern solton sng near elements of the modfed convecton-dffson eqaton α h d φ d φ + + = α s a parameter between and. If α =, we wll have the dscrete Galern form of
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 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 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 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 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 informationSolution of the Navier-Stokes Equations
Numercal Flud Mechancs Fall 2011 Lecture 25 REVIEW Lecture 24: Soluton of the Naver-Stokes Equatons Dscretzaton of the convectve and vscous terms Dscretzaton of the pressure term Conservaton prncples Momentum
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 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 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 information2.29 Numerical Fluid Mechanics Spring 2015 Lecture 10
REVIEW Lectre 9: Nmercal Fld Mechancs Srng 015 Lectre 10 End of (Lnear Algebrac Systems Gradent Methods Krylo Sbsace Methods Precondtonng of A=b FINITE DIFFERENCES Classfcaton of Partal Dfferental Eqatons
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 informationMEMBRANE ELEMENT WITH NORMAL ROTATIONS
9. MEMBRANE ELEMENT WITH NORMAL ROTATIONS Rotatons Mst Be Compatble Between Beam, Membrane and Shell Elements 9. INTRODUCTION { XE "Membrane Element" }The comple natre of most bldngs and other cvl engneerng
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 informationBAR & TRUSS FINITE ELEMENT. Direct Stiffness Method
BAR & TRUSS FINITE ELEMENT Drect Stness Method FINITE ELEMENT ANALYSIS AND APPLICATIONS INTRODUCTION TO FINITE ELEMENT METHOD What s the nte element method (FEM)? A technqe or obtanng approxmate soltons
More informationBezier curves. Michael S. Floater. August 25, These notes provide an introduction to Bezier curves. i=0
Bezer curves Mchael S. Floater August 25, 211 These notes provde an ntroducton to Bezer curves. 1 Bernsten polynomals Recall that a real polynomal of a real varable x R, wth degree n, s a functon of the
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 informationIntroduction to Turbulence Modelling
Introdcton to Trblence Modellng 1 Nmercal methods 0 1 t Mathematcal descrpton p F Reslts For eample speed, pressre, temperatre Geometry Models for trblence, combston etc. Mathematcal descrpton of physcal
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 informationPrinciples of Food and Bioprocess Engineering (FS 231) Solutions to Example Problems on Heat Transfer
Prncples of Food and Boprocess Engneerng (FS 31) Solutons to Example Problems on Heat Transfer 1. We start wth Fourer s law of heat conducton: Q = k A ( T/ x) Rearrangng, we get: Q/A = k ( T/ x) Here,
More information2.29 Numerical Fluid Mechanics
REVIEW Lecture 10: Sprng 2015 Lecture 11 Classfcaton of Partal Dfferental Equatons PDEs) and eamples wth fnte dfference dscretzatons Parabolc PDEs Ellptc PDEs Hyperbolc PDEs Error Types and Dscretzaton
More informationSalmon: Lectures on partial differential equations. Consider the general linear, second-order PDE in the form. ,x 2
Salmon: Lectures on partal dfferental equatons 5. Classfcaton of second-order equatons There are general methods for classfyng hgher-order partal dfferental equatons. One s very general (applyng even to
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 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 informationMMA and GCMMA two methods for nonlinear optimization
MMA and GCMMA two methods for nonlnear optmzaton Krster Svanberg Optmzaton and Systems Theory, KTH, Stockholm, Sweden. krlle@math.kth.se Ths note descrbes the algorthms used n the author s 2007 mplementatons
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 informationHigh resolution entropy stable scheme for shallow water equations
Internatonal Symposum on Computers & Informatcs (ISCI 05) Hgh resoluton entropy stable scheme for shallow water equatons Xaohan Cheng,a, Yufeng Ne,b, Department of Appled Mathematcs, Northwestern Polytechncal
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 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 informationResearch Article Cubic B-Spline Collocation Method for One-Dimensional Heat and Advection-Diffusion Equations
Appled Mathematcs Volume 22, Artcle ID 4587, 8 pages do:.55/22/4587 Research Artcle Cubc B-Splne Collocaton Method for One-Dmensonal Heat and Advecton-Dffuson Equatons Joan Goh, Ahmad Abd. Majd, and Ahmad
More informationThe equation of motion of a dynamical system is given by a set of differential equations. That is (1)
Dynamcal Systems Many engneerng and natural systems are dynamcal systems. For example a pendulum s a dynamcal system. State l The state of the dynamcal system specfes t condtons. For a pendulum n the absence
More informationPHYS 705: Classical Mechanics. Canonical Transformation II
1 PHYS 705: Classcal Mechancs Canoncal Transformaton II Example: Harmonc Oscllator f ( x) x m 0 x U( x) x mx x LT U m Defne or L p p mx x x m mx x H px L px p m p x m m H p 1 x m p m 1 m H x p m x m m
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 informationNew Method for Solving Poisson Equation. on Irregular Domains
Appled Mathematcal Scences Vol. 6 01 no. 8 369 380 New Method for Solvng Posson Equaton on Irregular Domans J. Izadan and N. Karamooz Department of Mathematcs Facult of Scences Mashhad BranchIslamc Azad
More informationBézier curves. Michael S. Floater. September 10, These notes provide an introduction to Bézier curves. i=0
Bézer curves Mchael S. Floater September 1, 215 These notes provde an ntroducton to Bézer curves. 1 Bernsten polynomals Recall that a real polynomal of a real varable x R, wth degree n, s a functon of
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 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 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 informationON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION
Advanced Mathematcal Models & Applcatons Vol.3, No.3, 2018, pp.215-222 ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EUATION
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 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 informationFUNDAMENTALS OF FINITE DIFFERENCE METHODS
FUNDAMENTALS OF FINITE DIFFERENCE METHODS By Deep Gupta 3 rd Year undergraduate, Mechanical Engg. Deptt., IIT Bombay Supervised by: Prof. Gautam Biswas, IIT Kanpur Acknowledgements It has been a pleasure
More informationExact Solutions for Nonlinear D-S Equation by Two Known Sub-ODE Methods
Internatonal Conference on Compter Technology and Scence (ICCTS ) IPCSIT vol. 47 () () IACSIT Press, Sngapore DOI:.7763/IPCSIT..V47.64 Exact Soltons for Nonlnear D-S Eqaton by Two Known Sb-ODE Methods
More informationSome modelling aspects for the Matlab implementation of MMA
Some modellng aspects for the Matlab mplementaton of MMA Krster Svanberg krlle@math.kth.se Optmzaton and Systems Theory Department of Mathematcs KTH, SE 10044 Stockholm September 2004 1. Consdered optmzaton
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 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 informationTHE STURM-LIOUVILLE EIGENVALUE PROBLEM - A NUMERICAL SOLUTION USING THE CONTROL VOLUME METHOD
Journal of Appled Mathematcs and Computatonal Mechancs 06, 5(), 7-36 www.amcm.pcz.pl p-iss 99-9965 DOI: 0.75/jamcm.06..4 e-iss 353-0588 THE STURM-LIOUVILLE EIGEVALUE PROBLEM - A UMERICAL SOLUTIO USIG THE
More informationGeorgia Tech PHYS 6124 Mathematical Methods of Physics I
Georga Tech PHYS 624 Mathematcal Methods of Physcs I Instructor: Predrag Cvtanovć Fall semester 202 Homework Set #7 due October 30 202 == show all your work for maxmum credt == put labels ttle legends
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 informationPhysics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1
P. Guterrez Physcs 5153 Classcal Mechancs D Alembert s Prncple and The Lagrangan 1 Introducton The prncple of vrtual work provdes a method of solvng problems of statc equlbrum wthout havng to consder the
More informationHandout # 6 (MEEN 617) Numerical Integration to Find Time Response of SDOF mechanical system. and write EOM (1) as two first-order Eqs.
Handout # 6 (MEEN 67) Numercal Integraton to Fnd Tme Response of SDOF mechancal system State Space Method The EOM for a lnear system s M X + DX + K X = F() t () t = X = X X = X = V wth ntal condtons, at
More informationALGORITHM FOR THE CALCULATION OF THE TWO VARIABLES CUBIC SPLINE FUNCTION
ANALELE ŞTIINŢIFICE ALE UNIVERSITĂŢII AL.I. CUZA DIN IAŞI (S.N.) MATEMATICĂ, Tomul LIX, 013, f.1 DOI: 10.478/v10157-01-00-y ALGORITHM FOR THE CALCULATION OF THE TWO VARIABLES CUBIC SPLINE FUNCTION BY ION
More informationModelli Clamfim Equazione del Calore Lezione ottobre 2014
CLAMFIM Bologna Modell 1 @ Clamfm Equazone del Calore Lezone 17 15 ottobre 2014 professor Danele Rtell danele.rtell@unbo.t 1/24? Convoluton The convoluton of two functons g(t) and f(t) s the functon (g
More informationChapter 4: Root Finding
Chapter 4: Root Fndng Startng values Closed nterval methods (roots are search wthn an nterval o Bsecton Open methods (no nterval o Fxed Pont o Newton-Raphson o Secant Method Repeated roots Zeros of Hgher-Dmensonal
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 informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 74 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS The Fnte Dfference Method These sldes are based on the recommended textbook: Culbert B. Laney. Computatonal
More informationFormal solvers of the RT equation
Formal solvers of the RT equaton Formal RT solvers Runge- Kutta (reference solver) Pskunov N.: 979, Master Thess Long characterstcs (Feautrer scheme) Cannon C.J.: 970, ApJ 6, 55 Short characterstcs (Hermtan
More informationPHYS 705: Classical Mechanics. Calculus of Variations II
1 PHYS 705: Classcal Mechancs Calculus of Varatons II 2 Calculus of Varatons: Generalzaton (no constrant yet) Suppose now that F depends on several dependent varables : We need to fnd such that has a statonary
More informationNUMERICAL MODELLING OF TSUNAMI WAVE EQUATIONS MASTER OF SCIENCE IN MATHEMATICS
NUMERICAL MODELLING OF TSUNAMI WAVE EQUATIONS A thess submtted n partal fulfllment of the requrements for the award of the degree of MASTER OF SCIENCE IN MATHEMATICS submtted by HARI SHANKAR SHAW Roll
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 informationDiscretization. Consistency. Exact Solution Convergence x, t --> 0. Figure 5.1: Relation between consistency, stability, and convergence.
Chapter 5 Theory The numercal smulaton of PDE s requres careful consderatons of propertes of the approxmate soluton. A necessary condton for the scheme used to model a physcal problem s the consstency
More informationInner Product. Euclidean Space. Orthonormal Basis. Orthogonal
Inner Product Defnton 1 () A Eucldean space s a fnte-dmensonal vector space over the reals R, wth an nner product,. Defnton 2 (Inner Product) An nner product, on a real vector space X s a symmetrc, blnear,
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 informationNumerical Simulation of One-Dimensional Wave Equation by Non-Polynomial Quintic Spline
IOSR Journal of Matematcs (IOSR-JM) e-issn: 78-578, p-issn: 319-765X. Volume 14, Issue 6 Ver. I (Nov - Dec 018), PP 6-30 www.osrournals.org Numercal Smulaton of One-Dmensonal Wave Equaton by Non-Polynomal
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 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 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 informationModule 3: Element Properties Lecture 1: Natural Coordinates
Module 3: Element Propertes Lecture : Natural Coordnates Natural coordnate system s bascally a local coordnate system whch allows the specfcaton of a pont wthn the element by a set of dmensonless numbers
More informationAdvanced Quantum Mechanics
Advanced Quantum Mechancs Rajdeep Sensarma! sensarma@theory.tfr.res.n ecture #9 QM of Relatvstc Partcles Recap of ast Class Scalar Felds and orentz nvarant actons Complex Scalar Feld and Charge conjugaton
More informationcoordinates. Then, the position vectors are described by
Revewng, what we have dscussed so far: Generalzed coordnates Any number of varables (say, n) suffcent to specfy the confguraton of the system at each nstant to tme (need not be the mnmum number). In general,
More informationNumerical Solution of two dimensional coupled viscous Burgers Equation using the Modified Cubic B-Spline Differential Quadrature Method
umercal Soluton of two dmensonal coupled vscous Burgers Equaton usng the odfed Cubc B-Splne Dfferental Quadrature ethod H. S. Shukla 1, ohammad Tamsr 1*, Vneet K. Srvastava, Ja Kumar 3 1 Department of
More informationPHYS 705: Classical Mechanics. Newtonian Mechanics
1 PHYS 705: Classcal Mechancs Newtonan Mechancs Quck Revew of Newtonan Mechancs Basc Descrpton: -An dealzed pont partcle or a system of pont partcles n an nertal reference frame [Rgd bodes (ch. 5 later)]
More informationIn this section is given an overview of the common elasticity models.
Secton 4.1 4.1 Elastc Solds In ths secton s gven an overvew of the common elastcty models. 4.1.1 The Lnear Elastc Sold The classcal Lnear Elastc model, or Hooean model, has the followng lnear relatonshp
More informationThe Exact Formulation of the Inverse of the Tridiagonal Matrix for Solving the 1D Poisson Equation with the Finite Difference Method
Journal of Electromagnetc Analyss and Applcatons, 04, 6, 0-08 Publshed Onlne September 04 n ScRes. http://www.scrp.org/journal/jemaa http://dx.do.org/0.46/jemaa.04.6000 The Exact Formulaton of the Inverse
More informationA new integrated-rbf-based domain-embedding scheme for solving fluid-flow problems
Home Search Collectons Journals About Contact us My IOPscence A new ntegrated-rbf-based doman-embeddng scheme for solvng flud-flow problems Ths artcle has been downloaded from IOPscence. Please scroll
More informationn α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0
MODULE 2 Topcs: Lnear ndependence, bass and dmenson We have seen that f n a set of vectors one vector s a lnear combnaton of the remanng vectors n the set then the span of the set s unchanged f that vector
More informationAERODYNAMICS I LECTURE 6 AERODYNAMICS OF A WING FUNDAMENTALS OF THE LIFTING-LINE THEORY
LECTURE 6 AERODYNAMICS OF A WING FUNDAMENTALS OF THE LIFTING-LINE THEORY The Bot-Savart Law The velocty nduced by the sngular vortex lne wth the crculaton can be determned by means of the Bot- Savart formula
More informationMore metrics on cartesian products
More metrcs on cartesan products If (X, d ) are metrc spaces for 1 n, then n Secton II4 of the lecture notes we defned three metrcs on X whose underlyng topologes are the product topology The purpose of
More information3 Basic boundary value problems for analytic function in the upper half plane
3 Basc boundary value problems for analytc functon n the upper half plane 3. Posson representaton formulas for the half plane Let f be an analytc functon of z throughout the half plane Imz > 0, contnuous
More informationArbitrary Lagrangian Eulerian Electromechanics in 3D
Progress In Electromagnetcs Research Symposum 2006, Cambrdge, USA, March 26-29 265 Arbtrary Lagrangan Euleran Electromechancs n 3D R. Reben, B. Walln, and D. Whte Lawrence Lvermore Natonal Laboratory,
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 informationEPR Paradox and the Physical Meaning of an Experiment in Quantum Mechanics. Vesselin C. Noninski
EPR Paradox and the Physcal Meanng of an Experment n Quantum Mechancs Vesseln C Nonnsk vesselnnonnsk@verzonnet Abstract It s shown that there s one purely determnstc outcome when measurement s made on
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 informationConsistency & Convergence
/9/007 CHE 374 Computatonal Methods n Engneerng Ordnary Dfferental Equatons Consstency, Convergence, Stablty, Stffness and Adaptve and Implct Methods ODE s n MATLAB, etc Consstency & Convergence Consstency
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 informationResearch Note NONLINEAR ANALYSIS OF SEMI-RIGID FRAMES WITH RIGID END SECTIONS * S. SEREN AKAVCI **
Iranan Jornal of Scence & Technology, Transacton, Engneerng, Vol., No., pp 7-7 rnted n The Islamc Repblc of Iran, 7 Shraz Unversty Research Note NONLINER NLYSIS OF SEMI-RIGID FRMES WIT RIGID END SECTIONS
More informationCHAPTER 6. LAGRANGE S EQUATIONS (Analytical Mechanics)
CHAPTER 6 LAGRANGE S EQUATIONS (Analytcal Mechancs) 1 Ex. 1: Consder a partcle movng on a fxed horzontal surface. r P Let, be the poston and F be the total force on the partcle. The FBD s: -mgk F 1 x O
More informationpage 2 2 dscretzaton mantans ths stablty under a sutable restrcton on the tme step. SSP tme dscretzaton methods were frst developed by Shu n [20] and
page 1 A Survey of Strong Stablty Preservng Hgh Order Tme Dscretzatons Ch-Wang Shu Λ 1 Introducton Numercal soluton for ordnary dfferental equatons (ODEs) s an establshed research area. There are many
More informationResearch Article A Multilevel Finite Difference Scheme for One-Dimensional Burgers Equation Derived from the Lattice Boltzmann Method
Appled Mathematcs Volume 01, Artcle ID 9590, 13 pages do:10.1155/01/9590 Research Artcle A Multlevel Fnte Dfference Scheme for One-Dmensonal Burgers Equaton Derved from the Lattce Boltzmann Method Qaoe
More informationCalculus of Variations Basics
Chapter 1 Calculus of Varatons Bascs 1.1 Varaton of a General Functonal In ths chapter, we derve the general formula for the varaton of a functonal of the form J [y 1,y 2,,y n ] F x,y 1,y 2,,y n,y 1,y
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 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 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 informationOptimal Control of Temperature in Fluid Flow
Kawahara Lab. 5 March. 27 Optmal Control of Temperature n Flud Flow Dasuke YAMAZAKI Department of Cvl Engneerng, Chuo Unversty Kasuga -3-27, Bunkyou-ku, Tokyo 2-855, Japan E-mal : d33422@educ.kc.chuo-u.ac.jp
More information