Finite element discretization of Laplace and Poisson equations
|
|
- Sarah Jenkins
- 6 years ago
- Views:
Transcription
1 Finit lmnt discrtization of Laplac and Poisson quations Yashwanth Tummala Tutor: Prof S.Mittal 1
2 Outlin Finit Elmnt Mthod for 1D Introduction to Poisson s and Laplac s Equations Finit Elmnt Mthod for 2D-Discrtization tc Wak formulation Intrpolation functions Boundary Conditions Assmbly of lmnt quations Exampls Summary 2 Finit Elmnt Discrtization of Laplac and Poisson Equations
3 Finit Elmnt Mthod For 1D Domain rprsntd by simpl domain/finit lmnts Approximation functions constructd ovr lmnts 1D quation looks lik: Construct Finit Elmnt Msh Driv th lmnt quations Assmbl th quations Impos th Boundary conditions Find th solutions d du a + cu f = 0, 0<x<L dx dx 3 Finit Elmnt Discrtization of Laplac and Poisson Equations
4 Introduction 4 What ar Poisson s and Laplac s Equations? ( k u) = f( x, y) in domain Ω, whr is th gradint oprator In Cartsian coordinat systm r r = i + j x y and th Poisson's Equation taks th form: u u k k = f( x, y) x x y y Whn f ( x, y) is zro it bcoms a Laplac's Equation! Finit Elmnt Discrtization of Laplac and Poisson Equations
5 Introduction Som xampls: Hat Transfr-Conduction, Convction (Tmp. distribution) Elctrostatics- Scalar potntial Magntostatics- Magntic potntial Fluid mchanics, Strss Distribution tc. Problm Statmnt: -Find th solution u(x,y) of th scond ordr partial diffrntial quation namly th Poisson s Equation -Givn k, f(x,y) and spcifid boundary conditions in a domain 5 Finit Elmnt Discrtization of Laplac and Poisson Equations
6 Finit Elmnt Discrtization Discrtization/Msh Gnration: 2D gomtric shaps- triangl, rctangl, quadrilatral tc. Gnral ruls: Elmnts should charactriz govrning quations Numbr, shap, typ, ar as rquird (accuracy) Larg gradints ar adquatly modld ( dnsity of lmnts) Msh rfinmnts should vary gradually Ω Ω $n Γ 6 Finit Elmnt Discrtization of Laplac and Poisson Equations
7 Wak Form 7 For a typical lmnt Ω Finit lmnt modl ovr Ω Ω,, u u w k k f dxdy = 0 x x y x W also hav th following, u w u u w k = k wk x x x x x x u w u u w k = k wk y y y y y y From gradint (or divrgnc thorm) in componnt form: Ω Ω u u wk dxdy = w k nxds x x x u u wk dxdy = w k nyds y y y Γ Γ Finit Elmnt Discrtization of Laplac and Poisson Equations
8 Wak Form 8 Whr n$ = n $ i + n $ j and n$ = unit normal vctor Thus w gt th wak form of Poisson's Equation: Ω x w u w u u u k + k wf dxdy w n k + n k ds = 0 x y x x y y Γ x y u u W say qn = nx k + x ny k y y Finit Elmnt Modl: u is approximatd n u( x, y) uh( x, y) = u jψ j( x, y) j = 1 u x y u nod x y th h(, ) = valu of h at j ( j, j) n = no. of nods of th lmnt. Ψ j = Lagrang intrpolation function/shap function. Finit Elmnt Discrtization of Laplac and Poisson Equations
9 Finit Elmnt Modl 9 Chos n j=1 { Ψ1 Ψ2 Ψn} linarly indpndnt functions w=,,... Th wak form bcoms: Or Ψ Ψ k + k dxdy u fψ dxdy Ψ q ds = 0 Ψi j Ψi j j i i n x x y y Ω Ω Γ n j=1 quation th this i ca ku = f + Q ( i= 1,2,..., n) ij j i i Unknowns : u,u,...,u In 1 2 n { K } { uj } = { f } + { Q } n b rwrittn as: th matrix form th st of Equations is: nxn nx1 nx1 nx1 Finit Elmnt Discrtization of Laplac and Poisson Equations
10 Intrpolation Functions 10 For a triangular lmnt: u( x, y) = c + c x + c y = u Ψ + u Ψ + u Ψ h u ( x, y ) = u = c + c x + c y, h h h i { } { }{ } u ( x, y ) = u = c + c x + c y, or u = A c 1 {} { } {} u ( x, y ) = u = c + c x + c y, or c = A u Aftr a bit of math w gt : α = xy x j k y, β = y y, γ = ( x x ) k j i j k i j k 1 c1 = ( α1u1+ α2u2 + α3u3),... A=ara of th triangl 2A From th first Equation abov w gt: 1 Ψ i = ( αi + βi x+ γi y),( i = 1,2,3) 2A Finit Elmnt Discrtization of Laplac and Poisson Equations
11 Intrpolation Functions 11 { } W calculat K and f for linar triangular lmnt. Th K ij associatd lmnt cofficint matrix is: k = + 4A 1 fi = fa 3 An xampl : ( βi βj γi γ j ) a + b b a 1 k fab K = 0, f 1 2ab = a 0 a [ ] b b { } Finit Elmnt Discrtization of Laplac and Poisson Equations
12 Intrpolation Functions 12 For a rctangular lmnt: u ( x, y) = c + c x + c y + c xy h u ( x, y) = c + c x + c y + c xy h u = u (0,0) = c, u = u ( a,0) = c + c a 1 h 1 2 h 1 2 u = u ( a, b) = c + c a + c b + c ab, u = u (0, b) = c + c b 3 h h 1 3 Thus w gt th intrpolation functions as: x y x y Ψ 1 = 1 1, Ψ 2 = 1, a b a b x y x y Ψ 3 =, Ψ 4 = 1 or ab a b i + 1 x i + x y i + y Ψ 1 ( x, y) = ( 1) 1 1 a b Finit Elmnt Discrtization of Laplac and Poisson Equations
13 Intrpolation Functions Th lmnt cofficint matrix for linar rctangular lmnt: ( a + b ) a b ( a + b ) b a ( + ) 2 ( + ) ( ) ( ) ( + ) 2 2( + ) a b a b b a a b k K = 6ab a + b b 2a 2 a + b a 2b b a a b a b a b Whn a=b, w furthr simplify to gt : k K = Finit Elmnt Discrtization of Laplac and Poisson Equations
14 Assmbly of Elmnt Equations Basd on two principls: Continuity of primary variabls Equilibrium or balanc of scondary variabls (lik flux) Infrncs: Continuity of u at nods continuity along intrlmnt boundary. At th intrfac, flux should b qual in magnitud and opposit in sign. Rlation/Cofficints of th stiffnss matrix ar thus drivd, taking global and local nod numbring into considration. 14 Finit Elmnt Discrtization of Laplac and Poisson Equations
15 Assmbly of Elmnt Equations For bttr undrstanding: With two lmnts w hav u = U, u = u = U, u = u = U, u = U, u = U continuity of u ( n) ( n) q = q Equilibrium of scondary variabls So, th assmbld quations ar givn by: K11 K12 K U1 F K21 K22 + K11 K23 + K14 K12 K U F F K3 1 K32+ K41 K33 + K44 K42 K U 43 3 = F3 + F K21 K24 K22 K23 U4 F K 2 31 K34 K32 K 33 U 5 F 3 15 Finit Elmnt Discrtization of Laplac and Poisson Equations
16 Boundary Conditions Th boundary intgral is of th form: Q i = q n Ψ i ( s ) d s Γ Whn is it not ncssary to comput??? Equilibrium of intrnal flux In gnral th intgral ovr th boundary is: Q = Ψ ( sq ) ( sds ) + Ψ ( sq ) ( sds ) + Ψ ( sq ) ( sds ) i i n i n i n i = i1 + i2 + i3 Q Q Q Q Q 1 What is th contribution of sid 2-3 to What about th rctangular lmnts? How many non-zro trms will thy hav? 16 Finit Elmnt Discrtization of Laplac and Poisson Equations
17 Exampl Considr th problm u = f0 i = f0 in a x y squar rgion. Using FEM find u(x,y), u=0 on th ntir boundary, is uniformly distributd sourc. f 0 u u Along th diagonal : u n u = = 0 x u 1 u u = = 0 n 2 y x u n u = = 0 y 17 On th lins of symmtry, normal drivativ of u(x,y)=0 Finit Elmnt Discrtization of Laplac and Poisson Equations
18 Exampl 18 Solution by Linar Triangular Elmnts: Symmtry along diagonal-modl as shown Msh as 4 linar triangls - local and global nods Elmnt 1 = Typical Elmnt W gt th lmnt coff matrix and sourc vctor: 2 2 b b fab 0 K = b a b a { f } 1 2ab + = a a 1 whr a = b = A/2 = 0.5 Thus w hav : { } { } { } { } K = K = K = K, f = f = f = f Finit Elmnt Discrtization of Laplac and Poisson Equations
19 Exampl 19 So, finally : f0 K = 1 2 1, { f } 1 2 = Th assmbld systm of quations look lik: U U 2 3 Q2 + Q3 + Q U 3 f 3 Q Q2 + Q1 = + 3 U 4 Q U Q1 + Q3 + Q U Q From th boundary conditions w gt; U = U = U = Finit Elmnt Discrtization of Laplac and Poisson Equations Q
20 Exampl 20 By boundary condition and balancing of intrnal Flux : Q = Q = 0, Q = Q + Q + Q = Q3 = Q3 + Q2 + Q1 = 0 Th abov st of quations can b writtn as: K11 K12 K13 U1 f K21 K22 + K33 + K11 K23 + K32 U2 = f2 + f3 + f K K32 K23 K33 K22 K 11 U f3 + f2 + f1 Us in g numrical valus for K and f ( with f = 1), W gt : 1 ij U U 2 = U 3 3 Finit Elmnt Discrtization of Laplac and Poisson Equations
21 Exampl 21 Solution by Linar Rctangular Elmnts: fa 1 0 K =, { f } = Obsrv that U = U = U = U = U = Boundary conditions on scondary variabls ar: Q = 0, Q + Q = 0, Q + Q = Balancing at global nod 5 givs: Q + Q + Q + Q = By t aking f = 1, a=0.5, w can solv for U 0 i u = 0 x Finit Elmnt Discrtization of Laplac and Poisson Equations u=0 u=0 u = 0 y
22 Conduction k Conduction in 1D is givn by: Conduction in 2D is givn by: u(x,y)= T = Tmpratur x q n, k f y T k x = f ( x, y ) x x T T k x k y = f( x, y ) x x y y = ngativ of hat flux = thrmal conductivitis along x and y = intrnal hat gnration Considr hat transfr in two-dimnsional plan. Find Tmpratur distribution in th plan. 22 Finit Elmnt Discrtization of Laplac and Poisson Equations
23 Convction 23 For a convction boundary: Equation for nrgy tranfsfr: T T k n + k n + β ( T T ) = q x y x x y y n β = convctiv conductanc, T = Ambint tmpratur q = spcifid hat flux Th W ak Form: Ω n w T w T kx + ky wf dxdy w qn β T T x x y y Γ Finit Elmnt Modl: n j=1 j ( x y) T= T Ψ, j Finit Elmnt Discrtization of Laplac and Poisson Equations ( ) ds = 0
24 Conduction/Convction So, w finally gt : n j = 1 ( ) K + H T = F + P ij ij j i i, whr, Ψ Ψ i j Ψ Ψ i j K ij = k x + k y dxdy x x y y Ω F = f Ψ dxdy + q Ψ ds = f + Q i i n i i i Ω Γ H = β Ψ Ψ ds, P = Ψ T ds ij i j i i Γ Γ 24 Finit Elmnt Discrtization of Laplac and Poisson Equations
25 Summary Finit Elmnt Modl for 1D problms Poisson s and Laplac s Equations (Modl Equation) Finit Elmnt Discrtization Wak Form Finit Elmnt Modl and Algbraic quations Intrpolation Functions for Triangular lmnts Intrpolation Functions for Rctangular lmnts Exampls For th Intrpolation/Shap functions 25 Finit Elmnt Discrtization of Laplac and Poisson Equations
26 Summary Associatd Elmnt Cofficint Matrix Boundary Conditions and Intgrals Assmbly of lmnt Equations and th principls Exampl with a squar domain and triangular lmnts Exampl with a squar domain and rctangular lmnts Application of FEM to Conduction, Conduction & Convction tc. 26 Finit Elmnt Discrtization of Laplac and Poisson Equations
27 Rfrncs An Introduction to th Finit Elmnt Mthod- 3 rd Edition by J.N.REDDY Th Finit Elmnt Mthod- 4 th Edition by Zinkiwicz - THANK YOU 27 Finit Elmnt Discrtization of Laplac and Poisson Equations
AS 5850 Finite Element Analysis
AS 5850 Finit Elmnt Analysis Two-Dimnsional Linar Elasticity Instructor Prof. IIT Madras Equations of Plan Elasticity - 1 displacmnt fild strain- displacmnt rlations (infinitsimal strain) in matrix form
More informationFEM FOR HEAT TRANSFER PROBLEMS دانشگاه صنعتي اصفهان- دانشكده مكانيك
FEM FOR HE RNSFER PROBLEMS 1 Fild problms Gnral orm o systm quations o D linar stady stat ild problms: For 1D problms: D D g Q y y (Hlmholtz quation) d D g Q d Fild problms Hat transr in D in h h ( D D
More informationMiddle East Technical University Department of Mechanical Engineering ME 413 Introduction to Finite Element Analysis
Middl East Tchnical Univrsity Dpartmnt of Mchanical Enginring ME Introduction to Finit Elmnt Analysis Chaptr 5 Two-Dimnsional Formulation Ths nots ar prpard by Dr. Cünyt Srt http://www.m.mtu.du.tr/popl/cunyt
More informationMiddle East Technical University Department of Mechanical Engineering ME 413 Introduction to Finite Element Analysis
Middl East Tchnical Univrsity Dpartmnt of Mchanical Enginring ME 43 Introduction to Finit Elmnt Analysis Chaptr 3 Computr Implmntation of D FEM Ths nots ar prpard by Dr. Cünyt Srt http://www.m.mtu.du.tr/popl/cunyt
More information16. Electromagnetics and vector elements (draft, under construction)
16. Elctromagntics (draft)... 1 16.1 Introduction... 1 16.2 Paramtric coordinats... 2 16.3 Edg Basd (Vctor) Finit Elmnts... 4 16.4 Whitny vctor lmnts... 5 16.5 Wak Form... 8 16.6 Vctor lmnt matrics...
More informationME469A Numerical Methods for Fluid Mechanics
ME469A Numrical Mthods for Fluid Mchanics Handout #5 Gianluca Iaccarino Finit Volum Mthods Last tim w introducd th FV mthod as a discrtization tchniqu applid to th intgral form of th govrning quations
More information1 Isoparametric Concept
UNIVERSITY OF CALIFORNIA BERKELEY Dpartmnt of Civil Enginring Spring 06 Structural Enginring, Mchanics and Matrials Profssor: S. Govindj Nots on D isoparamtric lmnts Isoparamtric Concpt Th isoparamtric
More informationCOMPUTATIONAL NUCLEAR THERMAL HYDRAULICS
COMPUTTIONL NUCLER THERML HYDRULICS Cho, Hyoung Kyu Dpartmnt of Nuclar Enginring Soul National Univrsity CHPTER4. THE FINITE VOLUME METHOD FOR DIFFUSION PROBLEMS 2 Tabl of Contnts Chaptr 1 Chaptr 2 Chaptr
More informationOTHER TPOICS OF INTEREST (Convection BC, Axisymmetric problems, 3D FEM)
OTHER TPOICS OF INTEREST (Convction BC, Axisymmtric problms, 3D FEM) CONTENTS 2-D Problms with convction BC Typs of Axisymmtric Problms Axisymmtric Problms (2-D) 3-D Hat Transfr 3-D Elasticity Typical
More informationMiddle East Technical University Department of Mechanical Engineering ME 413 Introduction to Finite Element Analysis
Middl East Tchnical Univrsity Dpartmnt of Mchanical Enginring ME 4 Introduction to Finit Elmnt Analysis Chaptr 4 Trusss, Bams and Frams Ths nots ar prpard by Dr. Cünyt Srt http://www.m.mtu.du.tr/popl/cunyt
More informationChapter 5. Introduction. Introduction. Introduction. Finite Element Modelling. Finite Element Modelling
Chaptr 5 wo-dimnsional problms using Constant Strain riangls (CS) Lctur Nots Dr Mohd Andi Univrsiti Malasia Prlis EN7 Finit Elmnt Analsis Introction wo-dimnsional init lmnt ormulation ollows th stps usd
More informationNumerical methods for PDEs FEM implementation: element stiffness matrix, isoparametric mapping, assembling global stiffness matrix
Platzhaltr für Bild, Bild auf Titlfoli hintr das Logo instzn Numrical mthods for PDEs FEM implmntation: lmnt stiffnss matrix, isoparamtric mapping, assmbling global stiffnss matrix Dr. Nomi Fridman Contnts
More informationDerivation of Eigenvalue Matrix Equations
Drivation of Eignvalu Matrix Equations h scalar wav quations ar φ φ η + ( k + 0ξ η β ) φ 0 x y x pq ε r r whr for E mod E, 1, y pq φ φ x 1 1 ε r nr (4 36) for E mod H,, 1 x η η ξ ξ n [ N ] { } i i i 1
More informationHydrogen Atom and One Electron Ions
Hydrogn Atom and On Elctron Ions Th Schrödingr quation for this two-body problm starts out th sam as th gnral two-body Schrödingr quation. First w sparat out th motion of th cntr of mass. Th intrnal potntial
More information3 Finite Element Parametric Geometry
3 Finit Elmnt Paramtric Gomtry 3. Introduction Th intgral of a matrix is th matrix containing th intgral of ach and vry on of its original componnts. Practical finit lmnt analysis rquirs intgrating matrics,
More informationHigher-Order Discrete Calculus Methods
Highr-Ordr Discrt Calculus Mthods J. Blair Prot V. Subramanian Ralistic Practical, Cost-ctiv, Physically Accurat Paralll, Moving Msh, Complx Gomtry, Slid 1 Contxt Discrt Calculus Mthods Finit Dirnc Mimtic
More informationFinite Element Models for Steady Flows of Viscous Incompressible Fluids
Finit Elmnt Modls for Stad Flows of Viscous Incomprssibl Fluids Rad: Chaptr 10 JN Rdd CONTENTS Govrning Equations of Flows of Incomprssibl Fluids Mid (Vlocit-Prssur) Finit Elmnt Modl Pnalt Function Mthod
More informationDifference -Analytical Method of The One-Dimensional Convection-Diffusion Equation
Diffrnc -Analytical Mthod of Th On-Dimnsional Convction-Diffusion Equation Dalabav Umurdin Dpartmnt mathmatic modlling, Univrsity of orld Economy and Diplomacy, Uzbistan Abstract. An analytical diffrncing
More informationIntroduction to Computational Fluid Dynamics: Governing Equations, Turbulence Modeling Introduction and Finite Volume Discretization Basics.
Introduction to Computational Fluid Dynamics: Govrning Equations, Turbulnc Modling Introduction and Finit Volum Discrtization Basics. Jol Gurrro Fbruary 13, 2014 Contnts 1 Notation and Mathmatical rliminaris
More informationFinite Element Analysis
Finit Elmnt Analysis L4 D Shap Functions, an Gauss Quaratur FEA Formulation Dr. Wiong Wu EGR 54 Finit Elmnt Analysis Roamap for Dvlopmnt of FE Strong form: govrning DE an BCs EGR 54 Finit Elmnt Analysis
More information2.3 Matrix Formulation
23 Matrix Formulation 43 A mor complicatd xampl ariss for a nonlinar systm of diffrntial quations Considr th following xampl Exampl 23 x y + x( x 2 y 2 y x + y( x 2 y 2 (233 Transforming to polar coordinats,
More informationElements of Statistical Thermodynamics
24 Elmnts of Statistical Thrmodynamics Statistical thrmodynamics is a branch of knowldg that has its own postulats and tchniqus. W do not attmpt to giv hr vn an introduction to th fild. In this chaptr,
More informationMAE4700/5700 Finite Element Analysis for Mechanical and Aerospace Design
MAE4700/5700 Finit Elmnt Analysis for Mchanical and Arospac Dsign Cornll Univrsity, Fall 2009 Nicholas Zabaras Matrials Procss Dsign and Control Laboratory Sibly School of Mchanical and Arospac Enginring
More informationDirect Approach for Discrete Systems One-Dimensional Elements
CONTINUUM & FINITE ELEMENT METHOD Dirct Approach or Discrt Systms On-Dimnsional Elmnts Pro. Song Jin Par Mchanical Enginring, POSTECH Dirct Approach or Discrt Systms Dirct approach has th ollowing aturs:
More informationNEW APPLICATIONS OF THE ABEL-LIOUVILLE FORMULA
NE APPLICATIONS OF THE ABEL-LIOUVILLE FORMULA Mirca I CÎRNU Ph Dp o Mathmatics III Faculty o Applid Scincs Univrsity Polithnica o Bucharst Cirnumirca @yahoocom Abstract In a rcnt papr [] 5 th indinit intgrals
More informationMath 34A. Final Review
Math A Final Rviw 1) Us th graph of y10 to find approimat valus: a) 50 0. b) y (0.65) solution for part a) first writ an quation: 50 0. now tak th logarithm of both sids: log() log(50 0. ) pand th right
More informationEinstein Equations for Tetrad Fields
Apiron, Vol 13, No, Octobr 006 6 Einstin Equations for Ttrad Filds Ali Rıza ŞAHİN, R T L Istanbul (Turky) Evry mtric tnsor can b xprssd by th innr product of ttrad filds W prov that Einstin quations for
More informationA Propagating Wave Packet Group Velocity Dispersion
Lctur 8 Phys 375 A Propagating Wav Packt Group Vlocity Disprsion Ovrviw and Motivation: In th last lctur w lookd at a localizd solution t) to th 1D fr-particl Schrödingr quation (SE) that corrsponds to
More informationNonlinear Bending of Strait Beams
Nonlinar Bnding of Strait Bams CONTENTS Th Eulr-Brnoulli bam thory Th Timoshnko bam thory Govrning Equations Wak Forms Finit lmnt modls Computr Implmntation: calculation of lmnt matrics Numrical ampls
More informationEngineering 323 Beautiful HW #13 Page 1 of 6 Brown Problem 5-12
Enginring Bautiful HW #1 Pag 1 of 6 5.1 Two componnts of a minicomputr hav th following joint pdf for thir usful liftims X and Y: = x(1+ x and y othrwis a. What is th probability that th liftim X of th
More informationAddition of angular momentum
Addition of angular momntum April, 0 Oftn w nd to combin diffrnt sourcs of angular momntum to charactriz th total angular momntum of a systm, or to divid th total angular momntum into parts to valuat th
More informationExercise 1. Sketch the graph of the following function. (x 2
Writtn tst: Fbruary 9th, 06 Exrcis. Sktch th graph of th following function fx = x + x, spcifying: domain, possibl asymptots, monotonicity, continuity, local and global maxima or minima, and non-drivability
More informationVSMN30 FINITA ELEMENTMETODEN - DUGGA
VSMN3 FINITA ELEMENTMETODEN - DUGGA 1-11-6 kl. 8.-1. Maximum points: 4, Rquird points to pass: Assistanc: CALFEM manual and calculator Problm 1 ( 8p ) 8 7 6 5 y 4 1. m x 1 3 1. m Th isotropic two-dimnsional
More informationBasic Polyhedral theory
Basic Polyhdral thory Th st P = { A b} is calld a polyhdron. Lmma 1. Eithr th systm A = b, b 0, 0 has a solution or thr is a vctorπ such that π A 0, πb < 0 Thr cass, if solution in top row dos not ist
More informationCHAPTER 2 LAGRANGIAN AND EULERIAN FINITE ELEMENTS IN ONE DIMENSION
CHAPTER 2 LAGRANGIAN AND EULERIAN FINITE ELEMENTS IN ONE DIMENSION by Td Blytschko Northwstrn Univrsity @ Copyright 1997 2.1 Introduction In this chaptr, th quations for on-dimnsional modls of nonlinar
More informationConstruction of Mimetic Numerical Methods
Construction of Mimtic Numrical Mthods Blair Prot Thortical and Computational Fluid Dynamics Laboratory Dltars July 17, 013 Numrical Mthods Th Foundation on which CFD rsts. Rvolution Math: Accuracy Stability
More informationAnalysis of potential flow around two-dimensional body by finite element method
Vol. 7(2), pp. 9-22, May, 2015 DOI: 10.5897/JMER2014.0342 rticl Numbr: 20E80053033 ISSN 2141 2383 Copyright 2015 uthor(s) rtain th copyright of this articl http://www.acadmicjournals.org/jmer Journal of
More informationCoupled Pendulums. Two normal modes.
Tim Dpndnt Two Stat Problm Coupld Pndulums Wak spring Two normal mods. No friction. No air rsistanc. Prfct Spring Start Swinging Som tim latr - swings with full amplitud. stationary M +n L M +m Elctron
More informationWhere k is either given or determined from the data and c is an arbitrary constant.
Exponntial growth and dcay applications W wish to solv an quation that has a drivativ. dy ky k > dx This quation says that th rat of chang of th function is proportional to th function. Th solution is
More informationData Assimilation 1. Alan O Neill National Centre for Earth Observation UK
Data Assimilation 1 Alan O Nill National Cntr for Earth Obsrvation UK Plan Motivation & basic idas Univariat (scalar) data assimilation Multivariat (vctor) data assimilation 3d-Variational Mthod (& optimal
More informationThe Matrix Exponential
Th Matrix Exponntial (with xrciss) by D. Klain Vrsion 207.0.05 Corrctions and commnts ar wlcom. Th Matrix Exponntial For ach n n complx matrix A, dfin th xponntial of A to b th matrix A A k I + A + k!
More informationThe pn junction: 2 Current vs Voltage (IV) characteristics
Th pn junction: Currnt vs Voltag (V) charactristics Considr a pn junction in quilibrium with no applid xtrnal voltag: o th V E F E F V p-typ Dpltion rgion n-typ Elctron movmnt across th junction: 1. n
More information(1) Then we could wave our hands over this and it would become:
MAT* K285 Spring 28 Anthony Bnoit 4/17/28 Wk 12: Laplac Tranform Rading: Kohlr & Johnon, Chaptr 5 to p. 35 HW: 5.1: 3, 7, 1*, 19 5.2: 1, 5*, 13*, 19, 45* 5.3: 1, 11*, 19 * Pla writ-up th problm natly and
More informationu 3 = u 3 (x 1, x 2, x 3 )
Lctur 23: Curvilinar Coordinats (RHB 8.0 It is oftn convnint to work with variabls othr than th Cartsian coordinats x i ( = x, y, z. For xampl in Lctur 5 w mt sphrical polar and cylindrical polar coordinats.
More informationThe Matrix Exponential
Th Matrix Exponntial (with xrciss) by Dan Klain Vrsion 28928 Corrctions and commnts ar wlcom Th Matrix Exponntial For ach n n complx matrix A, dfin th xponntial of A to b th matrix () A A k I + A + k!
More informationCHAPTER 1. Introductory Concepts Elements of Vector Analysis Newton s Laws Units The basis of Newtonian Mechanics D Alembert s Principle
CHPTER 1 Introductory Concpts Elmnts of Vctor nalysis Nwton s Laws Units Th basis of Nwtonian Mchanics D lmbrt s Principl 1 Scinc of Mchanics: It is concrnd with th motion of matrial bodis. odis hav diffrnt
More informationMEEN 617 Handout #12 The FEM in Vibrations A brief introduction to the finite element method for modeling of mechanical structures
MEEN 67 Handout # T FEM in Vibrations A brif introduction to t finit lmnt mtod for modling of mcanical structurs T finit lmnt mtod (FEM) is a picwis application of a variational mtod. Hr I provid you wit
More informationBSc Engineering Sciences A. Y. 2017/18 Written exam of the course Mathematical Analysis 2 August 30, x n, ) n 2
BSc Enginring Scincs A. Y. 27/8 Writtn xam of th cours Mathmatical Analysis 2 August, 28. Givn th powr sris + n + n 2 x n, n n dtrmin its radius of convrgnc r, and study th convrgnc for x ±r. By th root
More informationnd the particular orthogonal trajectory from the family of orthogonal trajectories passing through point (0; 1).
Eamn EDO. Givn th family of curvs y + C nd th particular orthogonal trajctory from th family of orthogonal trajctoris passing through point (0; ). Solution: In th rst plac, lt us calculat th di rntial
More information10. The Discrete-Time Fourier Transform (DTFT)
Th Discrt-Tim Fourir Transform (DTFT Dfinition of th discrt-tim Fourir transform Th Fourir rprsntation of signals plays an important rol in both continuous and discrt signal procssing In this sction w
More informationDynamic Modelling of Hoisting Steel Wire Rope. Da-zhi CAO, Wen-zheng DU, Bao-zhu MA *
17 nd Intrnational Confrnc on Mchanical Control and Automation (ICMCA 17) ISBN: 978-1-6595-46-8 Dynamic Modlling of Hoisting Stl Wir Rop Da-zhi CAO, Wn-zhng DU, Bao-zhu MA * and Su-bing LIU Xi an High
More informationDerivation of Electron-Electron Interaction Terms in the Multi-Electron Hamiltonian
Drivation of Elctron-Elctron Intraction Trms in th Multi-Elctron Hamiltonian Erica Smith Octobr 1, 010 1 Introduction Th Hamiltonian for a multi-lctron atom with n lctrons is drivd by Itoh (1965) by accounting
More informationMA 262, Spring 2018, Final exam Version 01 (Green)
MA 262, Spring 218, Final xam Vrsion 1 (Grn) INSTRUCTIONS 1. Switch off your phon upon ntring th xam room. 2. Do not opn th xam booklt until you ar instructd to do so. 3. Bfor you opn th booklt, fill in
More information1 General boundary conditions in diffusion
Gnral boundary conditions in diffusion Πρόβλημα 4.8 : Δίνεται μονοδιάτατη πλάκα πάχους, που το ένα άκρο της κρατιέται ε θερμοκραία T t και το άλλο ε θερμοκραία T 2 t. Αν η αρχική θερμοκραία της πλάκας
More informationSection 11.6: Directional Derivatives and the Gradient Vector
Sction.6: Dirctional Drivativs and th Gradint Vctor Practic HW rom Stwart Ttbook not to hand in p. 778 # -4 p. 799 # 4-5 7 9 9 35 37 odd Th Dirctional Drivativ Rcall that a b Slop o th tangnt lin to th
More informationIntroduction to Condensed Matter Physics
Introduction to Condnsd Mattr Physics pcific hat M.P. Vaughan Ovrviw Ovrviw of spcific hat Hat capacity Dulong-Ptit Law Einstin modl Dby modl h Hat Capacity Hat capacity h hat capacity of a systm hld at
More informationCHAPTER 9. Interpolation functions for 2D elements. Numerical Integration. Modeling Considerations
HAPTER 9 Intrpolation functions for D lmnts Numrical Intgration Modling onsidrations Pascal s triangl Dgr of Numbr of Elmnt with th complt trms in th nods polnomial polnomial Triangular Elmnts (Figur not
More informationCollisions between electrons and ions
DRAFT 1 Collisions btwn lctrons and ions Flix I. Parra Rudolf Pirls Cntr for Thortical Physics, Unirsity of Oxford, Oxford OX1 NP, UK This rsion is of 8 May 217 1. Introduction Th Fokkr-Planck collision
More informationDIFFERENTIAL EQUATION
MD DIFFERENTIAL EQUATION Sllabus : Ordinar diffrntial quations, thir ordr and dgr. Formation of diffrntial quations. Solution of diffrntial quations b th mthod of sparation of variabls, solution of homognous
More informationNONLINEAR ANALYSIS OF PLATE BENDING
NONLINEAR ANALYSIS OF PLATE BENDING CONTENTS Govrning Equations of th First-Ordr Shar Dformation thor (FSDT) Finit lmnt modls of FSDT Shar and mmbran locking Computr implmntation Strss calculation Numrical
More informationLecture 37 (Schrödinger Equation) Physics Spring 2018 Douglas Fields
Lctur 37 (Schrödingr Equation) Physics 6-01 Spring 018 Douglas Filds Rducd Mass OK, so th Bohr modl of th atom givs nrgy lvls: E n 1 k m n 4 But, this has on problm it was dvlopd assuming th acclration
More informationAddition of angular momentum
Addition of angular momntum April, 07 Oftn w nd to combin diffrnt sourcs of angular momntum to charactriz th total angular momntum of a systm, or to divid th total angular momntum into parts to valuat
More informationMor Tutorial at www.dumblittldoctor.com Work th problms without a calculator, but us a calculator to chck rsults. And try diffrntiating your answrs in part III as a usful chck. I. Applications of Intgration
More informationBackground: We have discussed the PIB, HO, and the energy of the RR model. In this chapter, the H-atom, and atomic orbitals.
Chaptr 7 Th Hydrogn Atom Background: W hav discussd th PIB HO and th nrgy of th RR modl. In this chaptr th H-atom and atomic orbitals. * A singl particl moving undr a cntral forc adoptd from Scott Kirby
More informationCOHORT MBA. Exponential function. MATH review (part2) by Lucian Mitroiu. The LOG and EXP functions. Properties: e e. lim.
MTH rviw part b Lucian Mitroiu Th LOG and EXP functions Th ponntial function p : R, dfind as Proprtis: lim > lim p Eponntial function Y 8 6 - -8-6 - - X Th natural logarithm function ln in US- log: function
More informationModule 8 Non equilibrium Thermodynamics
Modul 8 Non quilibrium hrmodynamics ctur 8.1 Basic Postulats NON-EQUIIRIBIUM HERMODYNAMICS Stady Stat procsss. (Stationary) Concpt of ocal thrmodynamic qlbm Extnsiv proprty Hat conducting bar dfin proprtis
More informationLecture 28 Title: Diatomic Molecule : Vibrational and Rotational spectra
Lctur 8 Titl: Diatomic Molcul : Vibrational and otational spctra Pag- In this lctur w will undrstand th molcular vibrational and rotational spctra of diatomic molcul W will start with th Hamiltonian for
More information5.80 Small-Molecule Spectroscopy and Dynamics
MIT OpnCoursWar http://ocw.mit.du 5.80 Small-Molcul Spctroscopy and Dynamics Fall 008 For information about citing ths matrials or our Trms of Us, visit: http://ocw.mit.du/trms. Lctur # 3 Supplmnt Contnts
More informationThat is, we start with a general matrix: And end with a simpler matrix:
DIAGON ALIZATION OF THE STR ESS TEN SOR INTRO DUCTIO N By th us of Cauchy s thorm w ar abl to rduc th numbr of strss componnts in th strss tnsor to only nin valus. An additional simplification of th strss
More informationMATH 1080 Test 2-SOLUTIONS Spring
MATH Tst -SOLUTIONS Spring 5. Considr th curv dfind by x = ln( 3y + 7) on th intrval y. a. (5 points) St up but do not simplify or valuat an intgral rprsnting th lngth of th curv on th givn intrval. =
More informationDivision of Mechanics Lund University MULTIBODY DYNAMICS. Examination Name (write in block letters):.
Division of Mchanics Lund Univrsity MULTIBODY DYNMICS Examination 7033 Nam (writ in block lttrs):. Id.-numbr: Writtn xamination with fiv tasks. Plas chck that all tasks ar includd. clan copy of th solutions
More informationProbability and Stochastic Processes: A Friendly Introduction for Electrical and Computer Engineers Roy D. Yates and David J.
Probability and Stochastic Procsss: A Frindly Introduction for Elctrical and Computr Enginrs Roy D. Yats and David J. Goodman Problm Solutions : Yats and Goodman,4.3. 4.3.4 4.3. 4.4. 4.4.4 4.4.6 4.. 4..7
More information1973 AP Calculus AB: Section I
97 AP Calculus AB: Sction I 9 Minuts No Calculator Not: In this amination, ln dnots th natural logarithm of (that is, logarithm to th bas ).. ( ) d= + C 6 + C + C + C + C. If f ( ) = + + + and ( ), g=
More informationOne Dimensional State Space Approach to Thermoelastic Interactions with Viscosity
7 IJSRST Volum 3 Issu 8 Print ISSN: 395-6 Onlin ISSN: 395-6X Thmd Sction: Scincand Tchnology On Dimnsional Stat Spac Approach to Thrmolastic Intractions with Viscosity Kavita Jain Rnu Yadav Dpartmnt of
More informationFourier Transforms and the Wave Equation. Key Mathematics: More Fourier transform theory, especially as applied to solving the wave equation.
Lur 7 Fourir Transforms and th Wav Euation Ovrviw and Motivation: W first discuss a fw faturs of th Fourir transform (FT), and thn w solv th initial-valu problm for th wav uation using th Fourir transform
More informationINC 693, 481 Dynamics System and Modelling: Linear Graph Modeling II Dr.-Ing. Sudchai Boonto Assistant Professor
INC 69, 48 Dynamics Systm and Modlling: Linar Graph Modling II Dr.-Ing. Sudchai Boonto Assistant Profssor Dpartmnt of Control Systm and Instrumntation Enginring King Mongkut s Unnivrsity of Tchnology Thonuri
More informationCalculus II (MAC )
Calculus II (MAC232-2) Tst 2 (25/6/25) Nam (PRINT): Plas show your work. An answr with no work rcivs no crdit. You may us th back of a pag if you nd mor spac for a problm. You may not us any calculators.
More informationSupplementary Materials
6 Supplmntary Matrials APPENDIX A PHYSICAL INTERPRETATION OF FUEL-RATE-SPEED FUNCTION A truck running on a road with grad/slop θ positiv if moving up and ngativ if moving down facs thr rsistancs: arodynamic
More informationChemical Physics II. More Stat. Thermo Kinetics Protein Folding...
Chmical Physics II Mor Stat. Thrmo Kintics Protin Folding... http://www.nmc.ctc.com/imags/projct/proj15thumb.jpg http://nuclarwaponarchiv.org/usa/tsts/ukgrabl2.jpg http://www.photolib.noaa.gov/corps/imags/big/corp1417.jpg
More informationDifferential Equations
UNIT I Diffrntial Equations.0 INTRODUCTION W li in a world of intrrlatd changing ntitis. Th locit of a falling bod changs with distanc, th position of th arth changs with tim, th ara of a circl changs
More informationExam 1. It is important that you clearly show your work and mark the final answer clearly, closed book, closed notes, no calculator.
Exam N a m : _ S O L U T I O N P U I D : I n s t r u c t i o n s : It is important that you clarly show your work and mark th final answr clarly, closd book, closd nots, no calculator. T i m : h o u r
More informationSME 3033 FINITE ELEMENT METHOD. Bending of Prismatic Beams (Initial notes designed by Dr. Nazri Kamsah)
Bnding of Prismatic Bams (Initia nots dsignd by Dr. Nazri Kamsah) St I-bams usd in a roof construction. 5- Gnra Loading Conditions For our anaysis, w wi considr thr typs of oading, as iustratd bow. Not:
More informationAnother view for a posteriori error estimates for variational inequalities of the second kind
Accptd by Applid Numrical Mathmatics in July 2013 Anothr viw for a postriori rror stimats for variational inqualitis of th scond kind Fi Wang 1 and Wimin Han 2 Abstract. In this papr, w giv anothr viw
More informationIntegration by Parts
Intgration by Parts Intgration by parts is a tchniqu primarily for valuating intgrals whos intgrand is th product of two functions whr substitution dosn t work. For ampl, sin d or d. Th rul is: u ( ) v'(
More informationElastic Analysis of Functionally Graded Variable Thickness Rotating Disk by Element Based Material Grading
Journal of Solid Mchanics ol. 9, No. 3 (017) pp. 650-66 Elastic Analysis of Functionally Gradd ariabl hicknss Rotating Disk by Elmnt Basd Matrial Grading A.K. hawait 1,*, L. Sondhi 1, Sh. Sanyal, Sh. Bhowmick
More informationFerroelectrics 342:73-82, 2006 Computational Modeling of Ferromagnetic Shape Memory Thin Films
Frrolctrics 4:7-8 6 Computational Modling of Frromagntic Shap Mmory Thin Films J. Liakhova M. Luskin and T. Zhang School of Mathmatics 6 Church St. SE Univrsity of Minnsota Minnapolis MN 55455 USA Email:
More informationSCALING OF SYNCHROTRON RADIATION WITH MULTIPOLE ORDER. J. C. Sprott
SCALING OF SYNCHROTRON RADIATION WITH MULTIPOLE ORDER J. C. Sprott PLP 821 Novmbr 1979 Plasma Studis Univrsity of Wisconsin Ths PLP Rports ar informal and prliminary and as such may contain rrors not yt
More informationHomotopy perturbation technique
Comput. Mthods Appl. Mch. Engrg. 178 (1999) 257±262 www.lsvir.com/locat/cma Homotopy prturbation tchniqu Ji-Huan H 1 Shanghai Univrsity, Shanghai Institut of Applid Mathmatics and Mchanics, Shanghai 272,
More informationBifurcation Theory. , a stationary point, depends on the value of α. At certain values
Dnamic Macroconomic Thor Prof. Thomas Lux Bifurcation Thor Bifurcation: qualitativ chang in th natur of th solution occurs if a paramtr passs through a critical point bifurcation or branch valu. Local
More informationMATH 319, WEEK 15: The Fundamental Matrix, Non-Homogeneous Systems of Differential Equations
MATH 39, WEEK 5: Th Fundamntal Matrix, Non-Homognous Systms of Diffrntial Equations Fundamntal Matrics Considr th problm of dtrmining th particular solution for an nsmbl of initial conditions For instanc,
More informationDiscontinuous Galerkin and Mimetic Finite Difference Methods for Coupled Stokes-Darcy Flows on Polygonal and Polyhedral Grids
Discontinuous Galrkin and Mimtic Finit Diffrnc Mthods for Coupld Stoks-Darcy Flows on Polygonal and Polyhdral Grids Konstantin Lipnikov Danail Vassilv Ivan Yotov Abstract W study locally mass consrvativ
More informationSECTION where P (cos θ, sin θ) and Q(cos θ, sin θ) are polynomials in cos θ and sin θ, provided Q is never equal to zero.
SETION 6. 57 6. Evaluation of Dfinit Intgrals Exampl 6.6 W hav usd dfinit intgrals to valuat contour intgrals. It may com as a surpris to larn that contour intgrals and rsidus can b usd to valuat crtain
More informationME 321 Kinematics and Dynamics of Machines S. Lambert Winter 2002
3.4 Forc Analysis of Linkas An undrstandin of forc analysis of linkas is rquird to: Dtrmin th raction forcs on pins, tc. as a consqunc of a spcifid motion (don t undrstimat th sinificanc of dynamic or
More informationNote If the candidate believes that e x = 0 solves to x = 0 or gives an extra solution of x = 0, then withhold the final accuracy mark.
. (a) Eithr y = or ( 0, ) (b) Whn =, y = ( 0 + ) = 0 = 0 ( + ) = 0 ( )( ) = 0 Eithr = (for possibly abov) or = A 3. Not If th candidat blivs that = 0 solvs to = 0 or givs an tra solution of = 0, thn withhold
More informationAs the matrix of operator B is Hermitian so its eigenvalues must be real. It only remains to diagonalize the minor M 11 of matrix B.
7636S ADVANCED QUANTUM MECHANICS Solutions Spring. Considr a thr dimnsional kt spac. If a crtain st of orthonormal kts, say, and 3 ar usd as th bas kts, thn th oprators A and B ar rprsntd by a b A a and
More information( ) Abstract. 2 FEDSS method basic relationships. 1 Introduction. 2.1 Tensorial formulation
Displacmnt basd continuous strss rcovry procdur, Mijuca D, Brkoviæ M. and Draškoviæ Z., Advancs in Finit Elmnt Tchnology, ISBN 0 948749 4, Ed. B.H.V.Topping, Civil Comp Prss, 7-34 (996). Abstract In this
More informationThe Transmission Line Wave Equation
1//5 Th Transmission Lin Wav Equation.doc 1/6 Th Transmission Lin Wav Equation Q: So, what functions I (z) and V (z) do satisfy both tlgraphr s quations?? A: To mak this asir, w will combin th tlgraphr
More informationThe Relativistic Stern-Gerlach Force C. Tschalär 1. Introduction
Th Rlativistic Strn-Grlach Forc C. Tschalär. Introduction For ovr a dcad, various formulations of th Strn-Grlach (SG) forc acting on a particl with spin moving at a rlativistic vlocity in an lctromagntic
More informationThe Implementation of Finite Element Method for Poisson Equation
Th Implmntation of Finit Elmnt Mthod for Poisson Equation Wnqiang Fng Abstract This is my MATH 574 cours projct rport. In this rport, I giv som dtails for implmnting th Finit Elmnt Mthod (FEM) via Matlab
More informationElectromagnetic scattering. Graduate Course Electrical Engineering (Communications) 1 st Semester, Sharif University of Technology
Elctromagntic scattring Graduat Cours Elctrical Enginring (Communications) 1 st Smstr, 1388-1389 Sharif Univrsity of Tchnology Contnts of lctur 8 Contnts of lctur 8: Scattring from small dilctric objcts
More information