Numerical Differentiation
|
|
- Gloria Haynes
- 5 years ago
- Views:
Transcription
1 Part 5 Capter 19 Numercal Derentaton PowerPonts organzed by Dr. Mcael R. Gustason II, Duke Unversty Revsed by Pro. Jang, CAU All mages copyrgt Te McGraw-Hll Companes, Inc. Permsson requred or reproducton or dsplay.
2 Capter Objectves Understandng d te applcaton o g-accuracy numercal derentaton ormulas or equspaced data. Knowng ow to evaluate dervatves or unequally spaced data. Understandng ow Rcardson etrapolaton s appled or numercal derentaton. Recognzng te senstvty o numercal derentaton to data error. Knowng ow to evaluate dervatves n MATLAB wt te d and gradent unctons. Knowng ow to generate contour plots and vector elds wt MATLAB.
3 Introducton to Derentaton Te one dmensonal orms o some consttutve laws commonly used Law Equaton Pyscal Area Gradent Flu Fourer s law dt qk d Heat conducton Temperature Heat Proportonal constant Termal conductvty t dc Fck s law J D d Mass duson Concentraton Mass Dusvty D Arcy s law d qk d Flow troug porous meda Head Flow Hydraulc conductvty Om s dv J Electrcal l Current low Voltage Current law d conductvty Newton s du vscosty t Fluds Vl Velocty Sear stress d law Hooke s law L E L Elastcty Deormaton Stress Dynamc vscosty Young s modulus
4 Derentaton Te matematcal denton o a dervatve begns wt a derence appromaton: y and as : ndependent varable s allowed to approac zero, te derence o or y: dependent varable becomes a dervatve: dy d lm 0
5 Hg-Accuracy Derentaton Formulas Taylor seres epanson can be used to generate g-accuracy ormulas or dervatves by usng lnear algebra to combne te epanson around several ponts. l l d Tree categores or te ormula nclude orward nte-derence, backward ntederence, and centered nte-derence.
6 Derentaton Tere are also backward derence and centered derence appromatons, dependng on te ponts used: Forward: ' 1 Backward: ' 1 1 Centered: O O ' O 2
7 Derentaton cont te rst order dervatve Forward derence appromaton o O O te rst order dervatve Backward derence appromaton o O O te rst order dervatve Centered derence appromaton o O O O Hger-Accuracy
8 Forward Fnte-Derence Hger-Accuracy
9 Backward Fnte-Derence
10 Centered Fnte-Derence
11 Eample /2 Q. Recall tat n E. 4.4 we estmated te dervatve o at =0.5 usng orward derences and a step sze o =0.25. Te results are summarzed n te table below. Te eact value o 0.5= Repeat te computaton wt g accuracy ormulas. Backward O Centered O 2 Forward O Estmate t 21.7% -2.4% -26.5%
12 Sol. 2 2 Eample / Forward derence o O 2 s computed as t =5.82 % Backward derence o O 2 s computed as t =3.77 % Centered derence o O 4 s computed as t =0 % Backward Centered Forward O O 2 O t 21.7% -2.4% -26.5%
13 Rcardson Etrapolaton As wt ntegraton, te Rcardson etrapolaton can be used to combne two lower-accuracy estmates o te dervatve to produce a ger-accuracy estmate. For te cases were tere are two O 2 estmates and te nterval s alved = /2, an mproved O estmate may be ormed usng: D 4 3 D D 1 For te cases were tere are two O 4 estmates and te nterval s alved 2 = 1 /2, an mproved O 6 estmate may be ormed usng: D D D 1 For te cases were tere are two O 6 estmates and te nterval s alved 2 = 1 /2, an mproved O 8 estmate may be ormed usng: D 64 D D 1
14 Eample 19.2 Q. Usng te same uncton as n E.19.1, estmate te rst dervatve at =0.5 or a step sze o 1 =0.5, and 2 =0.25. Use te Rcardson etrapolaton to compute mproved estmate. Te eact soluton s Sol. Te rst dervatve wt centered derence D t = 9.6% D t = 2.4% Usng te Rcardson etrapolaton, te mproved Estmate s 4 1 D D 2 D D
15 Unequally Spaced Data One way to calculated dervatves o unequally spaced data s to determne a polynomal l t and take ts dervatve at a pont. As an eample, usng a second-order Lagrange polynomal to t tree ponts and takng ts dervatve yelds:
16 Eample 19.3 Q. A temperature s measured nsde te sol as sown below. Compute te eat lu nto te ground at te ar-sol nterace. q z 0 dt k dz z 0 were q=eat lu W/m2, k=termal conductvty or sol =0.5 W/m K, T=TemperatureK, z=dstance measured rom te surace nto te sol K / m W W W qz mk m m 2
17 Dervatves and Integrals or Data wt Errors A sortcomng o numercal derentaton s tat t tends to amply errors n data, wereas ntegraton tends to smoot data errors. One approac or takng dervatves o data wt errors s to t a smoot, derentable uncton to te data and take te dervatve o te uncton. a Data wt no error b Resultng numercal derentaton o curve a c Data moded slgtly d Resultng numercal derentaton o curve a > Small data errors are ampled -> Small data errors are ampled by numercal derentaton.
18 Numercal Derentaton wt MATLAB MATLAB as two bult-n unctons to elp take dervatves, d and gradent: d Returns te derence between adjacent elements n dy./d Returns te derence between adjacent values n y dvded by te correspondng derence n adjacent values o
19 Numercal Derentaton wt MATLAB = gradent, Determnes te dervatve o te data n at eac o te ponts. Te program uses orward derence or te rst pont, backward derence or te last pont, and centered derence or te nteror ponts. s te spacng between ponts; omtted =1. Te major advantage o gradent over d s gradent s result s te same sze as te orgnal data. Gradent can also be used to nd partal dervatves or matrces: [, y] = gradent,
20 Vsualzaton MATLAB can generate contour plots o unctons as well as vector elds. Assumng and y represent a mesgrd o and y values and z represents a uncton o and y, contour, y, z can be used to generate a contour plot [, y]=gradentz, can be used to generate partal dervatves and quver, y,, y can be used to generate vector elds
ORDINARY DIFFERENTIAL EQUATIONS EULER S METHOD
Numercal Analss or Engneers German Jordanan Unverst ORDINARY DIFFERENTIAL EQUATIONS We wll eplore several metods o solvng rst order ordnar derental equatons (ODEs and we wll sow ow tese metods can be appled
More informationShuai Dong. Isaac Newton. Gottfried Leibniz
Computatonal pyscs Sua Dong Isaac Newton Gottred Lebnz Numercal calculus poston dervatve ntegral v velocty dervatve ntegral a acceleraton Numercal calculus Numercal derentaton Numercal ntegraton Roots
More 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 informationLecture 26 Finite Differences and Boundary Value Problems
4//3 Leture 6 Fnte erenes and Boundar Value Problems Numeral derentaton A nte derene s an appromaton o a dervatve - eample erved rom Talor seres 3 O! Negletng all terms ger tan rst order O O Tat s te orward
More informationCHAPTER 4d. ROOTS OF EQUATIONS
CHAPTER 4d. ROOTS OF EQUATIONS A. J. Clark School o Engneerng Department o Cvl and Envronmental Engneerng by Dr. Ibrahm A. Assakka Sprng 00 ENCE 03 - Computaton Methods n Cvl Engneerng II Department o
More informationOPTIMISATION. Introduction Single Variable Unconstrained Optimisation Multivariable Unconstrained Optimisation Linear Programming
OPTIMIATION Introducton ngle Varable Unconstraned Optmsaton Multvarable Unconstraned Optmsaton Lnear Programmng Chapter Optmsaton /. Introducton In an engneerng analss, sometmes etremtes, ether mnmum or
More informationNUMERICAL DIFFERENTIATION
NUMERICAL DIFFERENTIATION 1 Introducton Dfferentaton s a method to compute the rate at whch a dependent output y changes wth respect to the change n the ndependent nput x. Ths rate of change s called the
More informationNumerical Methods School of Mechanical Engineering Chung-Ang University
Part 5 Chapter 19 Numerical Differentiation Prof. Hae-Jin Choi hjchoi@cau.ac.kr 1 Chapter Objectives l Understanding the application of high-accuracy numerical differentiation formulas for equispaced data.
More informationSolution for singularly perturbed problems via cubic spline in tension
ISSN 76-769 England UK Journal of Informaton and Computng Scence Vol. No. 06 pp.6-69 Soluton for sngularly perturbed problems va cubc splne n tenson K. Aruna A. S. V. Rav Kant Flud Dynamcs Dvson Scool
More informationFinite Difference Method
7/0/07 Instructor r. Ramond Rump (9) 747 698 rcrump@utep.edu EE 337 Computatonal Electromagnetcs (CEM) Lecture #0 Fnte erence Method Lecture 0 These notes ma contan coprghted materal obtaned under ar use
More 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 informationSingle Variable Optimization
8/4/07 Course Instructor Dr. Raymond C. Rump Oce: A 337 Phone: (95) 747 6958 E Mal: rcrump@utep.edu Topc 8b Sngle Varable Optmzaton EE 4386/530 Computatonal Methods n EE Outlne Mathematcal Prelmnares Sngle
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 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 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 informationA New Recursive Method for Solving State Equations Using Taylor Series
I J E E E C Internatonal Journal of Electrcal, Electroncs ISSN No. (Onlne) : 77-66 and Computer Engneerng 1(): -7(01) Specal Edton for Best Papers of Mcael Faraday IET Inda Summt-01, MFIIS-1 A New Recursve
More informationUNIVERSITY OF BOLTON RAK ACADEMIC CENTRE BENG(HONS) MECHANICAL ENGINEERING SEMESTER TWO EXAMINATION 2017/2018 FINITE ELEMENT AND DIFFERENCE SOLUTIONS
OCD0 UNIVERSITY OF BOLTON RAK ACADEMIC CENTRE BENG(HONS) MECHANICAL ENGINEERING SEMESTER TWO EXAMINATION 07/08 FINITE ELEMENT AND DIFFERENCE SOLUTIONS MODULE NO. AME6006 Date: Wednesda 0 Ma 08 Tme: 0:00
More informationMathematical Economics MEMF e ME. Filomena Garcia. Topic 2 Calculus
Mathematcal Economcs MEMF e ME Flomena Garca Topc 2 Calculus Mathematcal Economcs - www.seg.utl.pt/~garca/economa_matematca . Unvarate Calculus Calculus Functons : X Y y ( gves or each element X one element
More information, rst we solve te PDE's L ad L ad n g g (x) = ; = ; ; ; n () (x) = () Ten, we nd te uncton (x), te lnearzng eedbac and coordnates transormaton are gve
Freedom n Coordnates Transormaton or Exact Lnearzaton and ts Applcaton to Transent Beavor Improvement Kenj Fujmoto and Tosaru Suge Dvson o Appled Systems Scence, Kyoto Unversty, Uj, Kyoto, Japan suge@robotuassyoto-uacjp
More informationMultivariate Ratio Estimator of the Population Total under Stratified Random Sampling
Open Journal of Statstcs, 0,, 300-304 ttp://dx.do.org/0.436/ojs.0.3036 Publsed Onlne July 0 (ttp://www.scrp.org/journal/ojs) Multvarate Rato Estmator of te Populaton Total under Stratfed Random Samplng
More 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 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 informationComplex Variables. Chapter 18 Integration in the Complex Plane. March 12, 2013 Lecturer: Shih-Yuan Chen
omplex Varables hapter 8 Integraton n the omplex Plane March, Lecturer: Shh-Yuan hen Except where otherwse noted, content s lcensed under a BY-N-SA. TW Lcense. ontents ontour ntegrals auchy-goursat theorem
More informationSummary with Examples for Root finding Methods -Bisection -Newton Raphson -Secant
Summary wth Eamples or Root ndng Methods -Bsecton -Newton Raphson -Secant Nonlnear Equaton Solvers Bracketng Graphcal Open Methods Bsecton False Poston (Regula-Fals) Newton Raphson Secant All Iteratve
More informationPHYS 1443 Section 004 Lecture #12 Thursday, Oct. 2, 2014
PHYS 1443 Secton 004 Lecture #1 Thursday, Oct., 014 Work-Knetc Energy Theorem Work under rcton Potental Energy and the Conservatve Force Gravtatonal Potental Energy Elastc Potental Energy Conservaton o
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 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 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 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 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 informationEE 330 Lecture 24. Small Signal Analysis Small Signal Analysis of BJT Amplifier
EE 0 Lecture 4 Small Sgnal Analss Small Sgnal Analss o BJT Ampler Eam Frda March 9 Eam Frda Aprl Revew Sesson or Eam : 6:00 p.m. on Thursda March 8 n Room Sweene 6 Revew rom Last Lecture Comparson o Gans
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 informationAerodynamics. Finite Wings Lifting line theory Glauert s method
α ( y) l Γ( y) r ( y) V c( y) β b 4 V Glauert s method b ( y) + r dy dγ y y dy Soluton procedure that transforms the lftng lne ntegro-dfferental equaton nto a system of algebrac equatons - Restrcted to
More informationCISE301: Numerical Methods Topic 2: Solution of Nonlinear Equations
CISE3: Numercal Methods Topc : Soluton o Nonlnear Equatons Dr. Amar Khoukh Term Read Chapters 5 and 6 o the tetbook CISE3_Topc c Khoukh_ Lecture 5 Soluton o Nonlnear Equatons Root ndng Problems Dentons
More information4 Round-Off and Truncation Errors
HK Km Slgtly moded 3//9, /8/6 Frstly wrtte at Marc 5 4 Roud-O ad Trucato Errors Errors Roud-o Errors Trucato Errors Total Numercal Errors Bluders, Model Errors, ad Data Ucertaty Recallg, dv dt Δv v t Δt
More informationBasic concept of reactive flows. Basic concept of reactive flows Combustion Mixing and reaction in high viscous fluid Application of Chaos
Introducton to Toshhsa Ueda School of Scence for Open and Envronmental Systems Keo Unversty, Japan Combuston Mxng and reacton n hgh vscous flud Applcaton of Chaos Keo Unversty 1 Keo Unversty 2 What s reactve
More informationInstituto Tecnológico de Aeronáutica FINITE ELEMENTS I. Class notes AE-245
Insttuto Tecnológco de Aeronáutca FIITE ELEMETS I Class notes AE-5 Insttuto Tecnológco de Aeronáutca 5. Isoparametrc Elements AE-5 Insttuto Tecnológco de Aeronáutca ISOPARAMETRIC ELEMETS Introducton What
More informationThe Finite Element Method: A Short Introduction
Te Fnte Element Metod: A Sort ntroducton Wat s FEM? Te Fnte Element Metod (FEM) ntroduced by engneers n late 50 s and 60 s s a numercal tecnque for solvng problems wc are descrbed by Ordnary Dfferental
More informationA Tale of Friction Basic Rollercoaster Physics. Fahrenheit Rollercoaster, Hershey, PA max height = 121 ft max speed = 58 mph
A Tale o Frcton Basc Rollercoaster Physcs Fahrenhet Rollercoaster, Hershey, PA max heght = 11 t max speed = 58 mph PLAY PLAY PLAY PLAY Rotatonal Movement Knematcs Smlar to how lnear velocty s dened, angular
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 informationTR/95 February Splines G. H. BEHFOROOZ* & N. PAPAMICHAEL
TR/9 February 980 End Condtons for Interpolatory Quntc Splnes by G. H. BEHFOROOZ* & N. PAPAMICHAEL *Present address: Dept of Matematcs Unversty of Tabrz Tabrz Iran. W9609 A B S T R A C T Accurate end condtons
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 10. Numerical Solution Methods for Engineering Analysis
Appled Engneerng Analss - sldes or class teachng* Chapter Numercal Soluton Methods or Engneerng Analss * Based on the tetbook on Appled Engneerng Analss, b Ta-Ran Hsu, publshed b John Wle & Sons, 8 (ISBN
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 informationTR/28. OCTOBER CUBIC SPLINE INTERPOLATION OF HARMONIC FUNCTIONS BY N. PAPAMICHAEL and J.R. WHITEMAN.
TR/8. OCTOBER 97. CUBIC SPLINE INTERPOLATION OF HARMONIC FUNCTIONS BY N. PAPAMICHAEL and J.R. WHITEMAN. W960748 ABSTRACT It s sown tat for te two dmensonal Laplace equaton a unvarate cubc splne approxmaton
More informationNumerical Solution of Singular Perturbation Problems Via Deviating Argument and Exponential Fitting
Amercan Journal of Computatonal and Appled Matematcs 0, (): 49-54 DOI: 0.593/j.ajcam.000.09 umercal Soluton of Sngular Perturbaton Problems Va Devatng Argument and Eponental Fttng GBSL. Soujanya, Y.. Reddy,
More informationAbsorbing Markov Chain Models to Determine Optimum Process Target Levels in Production Systems with Rework and Scrapping
Archve o SID Journal o Industral Engneerng 6(00) -6 Absorbng Markov Chan Models to Determne Optmum Process Target evels n Producton Systems wth Rework and Scrappng Mohammad Saber Fallah Nezhad a, Seyed
More informationOn a nonlinear compactness lemma in L p (0, T ; B).
On a nonlnear compactness lemma n L p (, T ; B). Emmanuel Matre Laboratore de Matématques et Applcatons Unversté de Haute-Alsace 4, rue des Frères Lumère 6893 Mulouse E.Matre@ua.fr 3t February 22 Abstract
More informationLesson 16: Basic Control Modes
0/8/05 Lesson 6: Basc Control Modes ET 438a Automatc Control Systems Technology lesson6et438a.tx Learnng Objectves Ater ths resentaton you wll be able to: Descrbe the common control modes used n analog
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 informationSolving Fractional Nonlinear Fredholm Integro-differential Equations via Hybrid of Rationalized Haar Functions
ISSN 746-7659 England UK Journal of Informaton and Computng Scence Vol. 9 No. 3 4 pp. 69-8 Solvng Fractonal Nonlnear Fredholm Integro-dfferental Equatons va Hybrd of Ratonalzed Haar Functons Yadollah Ordokhan
More informationSpring Force and Power
Lecture 13 Chapter 9 Sprng Force and Power Yeah, energy s better than orces. What s net? Course webste: http://aculty.uml.edu/andry_danylov/teachng/physcsi IN THIS CHAPTER, you wll learn how to solve problems
More informationCSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography
CSc 6974 and ECSE 6966 Math. Tech. for Vson, Graphcs and Robotcs Lecture 21, Aprl 17, 2006 Estmatng A Plane Homography Overvew We contnue wth a dscusson of the major ssues, usng estmaton of plane projectve
More information158 Calculus and Structures
58 Calculus and Structures CHAPTER PROPERTIES OF DERIVATIVES AND DIFFERENTIATION BY THE EASY WAY. Calculus and Structures 59 Copyrigt Capter PROPERTIES OF DERIVATIVES. INTRODUCTION In te last capter you
More informationCOSC 6374 Parallel Computation
COSC 67 Parallel Comptaton Partal Derental Eqatons Edgar Gabrel Fall 0 Nmercal derentaton orward derence ormla From te denton o dervatves one can derve an appromaton or te st dervatve Te same ormla can
More informationMain components of the above cycle are: 1) Boiler (steam generator) heat exchanger 2) Turbine generates work 3) Condenser heat exchanger 4) Pump
Introducton to Terodynacs, Lecture -5 Pro. G. Cccarell (0 Applcaton o Control olue Energy Analyss Most terodynac devces consst o a seres o coponents operatng n a cycle, e.g., stea power plant Man coponents
More 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 informationNON-PARABOLIC INTERFACE MOTION FOR THE 1-D STEFAN PROBLEM Dirichlet Boundary Conditions
Hernandez, E. M., et al.: Non-Parabolc Interace Moton or the -D Stean Problem... THERMAL SCIENCE: Year 07, Vol., No. 6A, pp. 37-336 37 NON-PARABOLIC INTERFACE MOTION FOR THE -D STEFAN PROBLEM Drchlet Boundary
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 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 informationParameter estimation class 5
Parameter estmaton class 5 Multple Ve Geometr Comp 9-89 Marc Pollefes Content Background: Projectve geometr (D, 3D), Parameter estmaton, Algortm evaluaton. Sngle Ve: Camera model, Calbraton, Sngle Ve Geometr.
More informationElectrical Circuits 2.1 INTRODUCTION CHAPTER
CHAPTE Electrcal Crcuts. INTODUCTION In ths chapter, we brefly revew the three types of basc passve electrcal elements: resstor, nductor and capactor. esstance Elements: Ohm s Law: The voltage drop across
More informationSENSITIVITY APPROACH TO OPTIMAL CONTROL FOR AFFINE NONLINEAR DISCRETE-TIME SYSTEMS
448 Asan Journal o Control Vol. 7 No. 4 pp. 448-454 December 25 Bre Paper SENSIIVIY APPROACH O OPIMAL CONROL FOR AFFINE NONLINEAR DISCREE-IME SYSEMS Gong-You ang Nan Xe Peng Lu ABSRAC hs paper deals wth
More informationECEN 5005 Crystals, Nanocrystals and Device Applications Class 19 Group Theory For Crystals
ECEN 5005 Crystals, Nanocrystals and Devce Applcatons Class 9 Group Theory For Crystals Dee Dagram Radatve Transton Probablty Wgner-Ecart Theorem Selecton Rule Dee Dagram Expermentally determned energy
More informationReport on Image warping
Report on Image warpng Xuan Ne, Dec. 20, 2004 Ths document summarzed the algorthms of our mage warpng soluton for further study, and there s a detaled descrpton about the mplementaton of these algorthms.
More information2.29 Numerical Fluid Mechanics Fall 2011 Lecture 6
REVIEW of Lecture 5 2.29 Numercal Flud Mechancs Fall 2011 Lecture 6 Contnuum Hypothess and conservaton laws Macroscopc Propertes Materal covered n class: Dfferental forms of conservaton laws Materal Dervatve
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 informationUniversity of Washington Department of Chemistry Chemistry 452/456 Summer Quarter 2014
Lecture 16 8/4/14 Unversty o Washngton Department o Chemstry Chemstry 452/456 Summer Quarter 214. Real Vapors and Fugacty Henry s Law accounts or the propertes o extremely dlute soluton. s shown n Fgure
More information6.3.4 Modified Euler s method of integration
6.3.4 Modfed Euler s method of ntegraton Before dscussng the applcaton of Euler s method for solvng the swng equatons, let us frst revew the basc Euler s method of numercal ntegraton. Let the general from
More informationThe Karush-Kuhn-Tucker. Nuno Vasconcelos ECE Department, UCSD
e Karus-Kun-ucker condtons and dualt Nuno Vasconcelos ECE Department, UCSD Optmzaton goal: nd mamum or mnmum o a uncton Denton: gven unctons, g, 1,...,k and, 1,...m dened on some doman Ω R n mn w, w Ω
More 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 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 informationWork is the change in energy of a system (neglecting heat transfer). To examine what could
Work Work s the change n energy o a system (neglectng heat transer). To eamne what could cause work, let s look at the dmensons o energy: L ML E M L F L so T T dmensonally energy s equal to a orce tmes
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 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 informationPubH 7405: REGRESSION ANALYSIS. SLR: INFERENCES, Part II
PubH 7405: REGRESSION ANALSIS SLR: INFERENCES, Part II We cover te topc of nference n two sessons; te frst sesson focused on nferences concernng te slope and te ntercept; ts s a contnuaton on estmatng
More informationUniversity of Washington Department of Chemistry Chemistry 452/456 Summer Quarter 2014
Lecture 12 7/25/14 ERD: 7.1-7.5 Devoe: 8.1.1-8.1.2, 8.2.1-8.2.3, 8.4.1-8.4.3 Unversty o Washngton Department o Chemstry Chemstry 452/456 Summer Quarter 2014 A. Free Energy and Changes n Composton: The
More informationChapter 20 The First Law of Thermodynamics
Chapter he Frst aw o hermodynamcs. developng the concept o heat. etendng our concept o work to thermal processes 3. ntroducng the rst law o thermodynamcs. Heat and Internal Energy Internal energy: s the
More informationPowerPoints organized by Dr. Michael R. Gustafson II, Duke University
Part 5 Chapter 19 Numerical Differentiation PowerPoints organized by Dr. Michael R. Gustafson II, Duke University All images copyright The McGraw-Hill Companies, Inc. Permission required for reproduction
More informationNumerical Methods Solution of Nonlinear Equations
umercal Methods Soluton o onlnear Equatons Lecture Soluton o onlnear Equatons Root Fndng Prolems Dentons Classcaton o Methods Analytcal Solutons Graphcal Methods umercal Methods Bracketng Methods Open
More informationPHYS 1101 Practice problem set 12, Chapter 32: 21, 22, 24, 57, 61, 83 Chapter 33: 7, 12, 32, 38, 44, 49, 76
PHYS 1101 Practce problem set 1, Chapter 3: 1,, 4, 57, 61, 83 Chapter 33: 7, 1, 3, 38, 44, 49, 76 3.1. Vsualze: Please reer to Fgure Ex3.1. Solve: Because B s n the same drecton as the ntegraton path s
More informationHongyi Miao, College of Science, Nanjing Forestry University, Nanjing ,China. (Received 20 June 2013, accepted 11 March 2014) I)ϕ (k)
ISSN 1749-3889 (prnt), 1749-3897 (onlne) Internatonal Journal of Nonlnear Scence Vol.17(2014) No.2,pp.188-192 Modfed Block Jacob-Davdson Method for Solvng Large Sparse Egenproblems Hongy Mao, College of
More informationThermal-Fluids I. Chapter 18 Transient heat conduction. Dr. Primal Fernando Ph: (850)
hermal-fluds I Chapter 18 ransent heat conducton Dr. Prmal Fernando prmal@eng.fsu.edu Ph: (850) 410-6323 1 ransent heat conducton In general, he temperature of a body vares wth tme as well as poston. In
More informationPolynomial Interpolation
Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximatinga function fx, wose values at a set of distinct points x, x, x,, x n are known, by a polynomial P x suc
More informationJournal of Universal Computer Science, vol. 1, no. 7 (1995), submitted: 15/12/94, accepted: 26/6/95, appeared: 28/7/95 Springer Pub. Co.
Journal of Unversal Computer Scence, vol. 1, no. 7 (1995), 469-483 submtted: 15/12/94, accepted: 26/6/95, appeared: 28/7/95 Sprnger Pub. Co. Round-o error propagaton n the soluton of the heat equaton by
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 informationChapter Runge-Kutta 2nd Order Method for Ordinary Differential Equations
Cter. Runge-Kutt nd Order Metod or Ordnr Derentl Eutons Ater redng ts cter ou sould be ble to:. understnd te Runge-Kutt nd order metod or ordnr derentl eutons nd ow to use t to solve roblems. Wt s te Runge-Kutt
More informationTEST 5 (phy 240) 2. Show that the volume coefficient of thermal expansion for an ideal gas at constant pressure is temperature dependent and given by
ES 5 (phy 40). a) Wrte the zeroth law o thermodynamcs. b) What s thermal conductvty? c) Identyng all es, draw schematcally a P dagram o the arnot cycle. d) What s the ecency o an engne and what s the coecent
More informationrisk and uncertainty assessment
Optmal forecastng of atmospherc qualty n ndustral regons: rsk and uncertanty assessment Vladmr Penenko Insttute of Computatonal Mathematcs and Mathematcal Geophyscs SD RAS Goal Development of theoretcal
More informationONE DIMENSIONAL TRIANGULAR FIN EXPERIMENT. Technical Advisor: Dr. D.C. Look, Jr. Version: 11/03/00
ONE IMENSIONAL TRIANGULAR FIN EXPERIMENT Techncal Advsor: r..c. Look, Jr. Verson: /3/ 7. GENERAL OJECTIVES a) To understand a one-dmensonal epermental appromaton. b) To understand the art of epermental
More informationFor secondary motion, physical effects Simulation (physically based animation) For character animation. Motion capturing / motion synthesis
3 4 Introducton to Anmaton Technques For character anmaton Keyramng Moton capturng / moton synthess For secondary moton physcal eects Procedural Smulaton (physcally based anmaton) Motvaton Anmaton rom
More informationInitial-value Technique For Singularly Perturbed Two Point Boundary Value Problems Via Cubic Spline
Unversty of Central Florda Electronc Teses and Dssertatons Masters Tess (Open Access) Intal-value Tecnque For Sngularly Perturbed Two Pont Boundary Value Problems Va Cubc Splne 010 Lus G. Negron Unversty
More informationGenerative and Discriminative Models. Jie Tang Department of Computer Science & Technology Tsinghua University 2012
Generatve and Dscrmnatve Models Je Tang Department o Computer Scence & Technolog Tsnghua Unverst 202 ML as Searchng Hpotheses Space ML Methodologes are ncreasngl statstcal Rule-based epert sstems beng
More informationMultigrid Methods and Applications in CFD
Multgrd Metods and Applcatons n CFD Mcael Wurst 0 May 009 Contents Introducton Typcal desgn of CFD solvers 3 Basc metods and ter propertes for solvng lnear systems of equatons 4 Geometrc Multgrd 3 5 Algebrac
More informationSolving Singularly Perturbed Differential Difference Equations via Fitted Method
Avalable at ttp://pvamu.edu/aam Appl. Appl. Mat. ISSN: 193-9466 Vol. 8, Issue 1 (June 013), pp. 318-33 Applcatons and Appled Matematcs: An Internatonal Journal (AAM) Solvng Sngularly Perturbed Dfferental
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 informationORDINARY DIFFERENTIAL EQUATIONS
6 ORDINARY DIFFERENTIAL EQUATIONS Introducton Runge-Kutt Metods Mult-step Metods Sstem o Equtons Boundr Vlue Problems Crcterstc Vlue Problems Cpter 6 Ordnr Derentl Equtons / 6. Introducton In mn engneerng
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 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 information5 The Laplace Equation in a convex polygon
5 Te Laplace Equaton n a convex polygon Te most mportant ellptc PDEs are te Laplace, te modfed Helmoltz and te Helmoltz equatons. Te Laplace equaton s u xx + u yy =. (5.) Te real and magnary parts of an
More information