Numerical Methods for Partial Differential Equations
|
|
- Kathryn Perry
- 6 years ago
- Views:
Transcription
1 Numerical Methods for Partial Differential Equations Eric de Sturler University of Illinois at Urbana-Champaign The calculus of variations deals with maxima, minima, and stationary values of (definite) integrals. Consider the following minimization problem. Find the minimum over all (sufficiently smooth and bounded) functions v(x) of t f (v, v t, x)dx (or f (v, v )dx as in the book), such that v() = v and v() = v. Note that this is (again) a minimization problem over a space of functions. It turns out that this problem can be replaced by an equivalent partial differential equation. Also f(r, s, x) must be sufficiently smooth. Typically, the integral measures the total energy of some system. Eric de Sturler 2-2 /27/
2 The basic idea is that we add an arbitrary other function (a variation) to the unknown minimizer (maximizer, stationarizer(?)) and derive the conditions that this should always increase (for minimizer,...) the value of the integral. So assume u(x) with u() = u and u() = u is the minimizer and add v(x) to u(x) for an arbitary function v(x) with ævæ =, v() = and v() =. If u(x) solves our problem, we must have the following d t d f (u + v, u + v, x)dx = =. The directional derivative in any direction must vanish at u ( = ). The conditions this places on u(x) allow us to compute this function by deriving an equivalent partial differential equation. Eric de Sturler 2 Given the regularity assumptions on u and v we can differentiate under the integral to get d t d f (u + v, u + v, x)dx = = v + v dx, Integration by parts of the second term gives v dx = dv = [f u v] This gives for the whole integral vd = d dx vdx d t d f (u + v, u + v, x)dx = = v + v dx = d dx vdx =. Eric de Sturler /27/
3 Now the equation d dx vdx =, for all v, means that d dx = for all x. Otherwise, assume d dx > (or < ) over some small interval. We can choose v such that it is nonzero over (part of) that interval and zero everywhere else. Clearly, this would violate the equation above. Hence, u is given by the partial differential equation d dx = for x c (, ) and u() = u and. u() = u This is called the Euler-Lagrange equation for the minimization problem. Eric de Sturler 2 Now consider the minimization problem, find u(x) with u() = u, u () = u, u() = u, u () = u, that minimizes F (u, u, u, x)dx. d We consider again F (u + v, u + v, u + v, x)dx =. d This gives (note v = v = at the end points) v + v + v dx = v dx = g v d dx v + v dx = dv = v d dx v dx = d dx = dv = d 2 dx 2 vdx Substituting this gives d dx + d2 for all : dx 2 vdx = v d dx + d2 dx 2 for x c (, ). Eric de Sturler /27/
4 We also have the concept of essential and natural boundary conditions: In the case of fixed boundary values (Dirichlet) we see that that the boundary values must be enforced by restricting the trial (and test) space. Now consider we minimize f (u, u, x)dx without requiring specific boundary values. The analogous procedure gives d dx vdx +, for all. v = v Now taking functions v such that v() = v() = implies d dx =. This leaves the condition v, which must hold for arbitrary values = of v() and v() since there are no (such) restrictions on v. This implies = at the end points (natural boundary conditions). Eric de Sturler 2 Of course, more general boundary conditions can also be included. Note that the procedure described gives stationary points. In many cases (especially static problems) the functional is quadratic or, more generally, convex and the stationary point (function) gives a minimum. For example, the minimum of the total energy of some mechanical system. In other cases (especially dynamic problems) the functional is not minimized but we find a stationary point (function) of a functional. For example a stationary point of the difference between kinetic and potential energy. Eric de Sturler /27/
5 The first example is the total energy of an elastic bar under tangential load ; see section Find (with ) that minimizes f(x) v(x) v() = v() = the total energy given by I (v, v )dx = [ 2 av2 fv]dx The integral represents the internal elastic energy and the load potential. ØI We have Øv = f ØI and Øv = av, so the Euler-Lagrange equation gives ØI Øv ØI Øv dx = f av dx = (which is the weak formulation) Integration by parts then gives [ f ( av ) ] dx = for all : f (av ) = with v() = v() = Eric de Sturler 2 The second example concerns the movement of a system consisting of a mass connected to a spring moving frictionless over some (hor.) surface. Let r(t) denote the displacement, such that r = means no force applied. Kinetic energy: T = 2 mṙ2 = 2 mṙ (t) 2 Potential energy: U = 2 kr2 = 2 kr (t) 2 Total energy: E = T + U = 2 mṙ2 + 2 kr2 Langrangian: L = T U = 2 mr.2 2 kr2 We assume as initial conditions r() = r and r () = r. Now Hamilton s principle states that we look for a stationary point of the integral over the Lagrangian. Eric de Sturler 2 9- /27/
6 Find r such that r() = r and ṙ() = ṙ and r is a stationary point of L(r, ṙ) = 2mṙ 2 2kr 2 dt Using variations with () = () = the Euler-Lagrange equation gives. kr + mṙ dt = e mr + kr = Note that for this system the total energy is constant in time: d dt E = mṙr + krṙ = ṙ (mr + kr) = Excluding the less useful choice ṙ = we find again mr + kr =. Eric de Sturler 2 Clearly these two derivations are related. Hamilton showed that the following coupled system of first order equations is equivalent to the Euler-Lagrange equations (in general). Let H(p, r, t) = pṙ L(r, ṙ, t) (the Legendre transformation of L(ṙ) ), where p = ØL Øṙ and ṗ = ØL Ør. Then the coupled system ṗ = ØH Ør and q. = ØH Øp is equivalent to the EL equation. For our problen ṗ = ØH Ør = ØL Ør = kr is the force exerted by the spring, and p = ØL Øṙ = mṙ is the momentum of the system. Now if the Lagrangian represents a mechanical system with L = T U, the difference between kinetic and potential energy, and the kinetic energy is a quadratic function of ṙ, then the Hamiltonian H represents the total energy: H = ØL Ør. ṙ + L = mṙ 2 L = 2T (T U) = T + U = E. Eric de Sturler 2-2 /27/
7 These characterizations of static and dynamic systems are the subject of analytical mechanics (introduced by Euler and Lagrange). For more details see (two excellent books): The Variational Principles of Mechanics, C. Lanczos, Dover Mathematical Methods of Classical Mechanics, V.I. Arnold, Springer Another good book discussing these minimizations and their relation to differential equations in weak form and finite element solution is Finite Element Solution of Boundary Value Problems (Theory and Computation), O. Axelsson and V.A. Baker, SIAM (Classics in Applied Mathematics) Eric de Sturler 2 We now briefly consider the variational method in higher dimensions. F(u, u x, u y, x, y)dxdy for x c u = for x c Ø Adding a function with = on the boundary and taking the directional derivative gives + x x + y y dxdy = Applying Gauß theorem gives and (weak formulation) Ø Ø x x dxdy = Ø x n ds Øx x dxdy = Øx x dxdy Ø Ø y y dxdy = Ø y n 2 ds Øy y dxdy = Øy y dxdy Ø Øx x Ø Øy y dxdy = Eric de Sturler /27/
8 The following theorems are very useful Gauss: $ u d = Ø u $ n ds Where d is the unit surface/volume element, and ds is the unit line/surface element. u If we apply Gauss theorem to the function v we get u 2 v 2 v Øx + 2v 2 Øx 2 d = Ø (u v )n + (u 2 v 2 )n 2 ds Now assume v 2 = (constant): v Øx d = Øx v + u Øv Øx d = Ø (u v )n ds uv x d = Ø (uv)n ds u x vd, or applying the same idea for y uv y d = Ø (uv)n 2 ds u y vd Eric de Sturler 2 Clearly, exactly the same derivation in three dimensions will give Ø Øx x Ø Øy y Ø Øz z dxdydz = So the Euler-Lagrange equation becomes d i= Ø Øx i x i = The boundary condition u = for x c Ø is an essential boundary condition. Removing the constraint on the boundary value leads to the natural boundary condition(s) (in 2D) Ø x n + y n 2 ds = for all implies x n + y n 2 =. Eric de Sturler /27/
9 We can also use Gauss theorem to derive the weak form of the Poisson equation. u = f for x c The weak form is given by v u d = fv d for all v (of some space) Now apply Gauss to $ (v u) d = Ø (v u) $ n ds w $ v Øx + v Øx 2 d = Ø Øx v Øx + Ø Øx 2 v Øx 2 d = Øv Øx Øx + Øv Øx 2 Øx 2 d + v Ø2 u Øx 2 + v Ø2 u Øx2 d = Ø (v u) $ n ds 2 Finally, this gives v $ u d = Ø (v u) $ n ds vßu d and since v is a scalar function v $ u d = Ø v Øn ds vßu d Eric de Sturler 2 Poisson equation: u = f for x c The weak form (after partial applying Gauss theorem) is then given by vßu d = v $ u d + Ø v Øn ds = fv d Now taking Dirichlet boundary conditions, u given and v = on the boundary, or taking Neumann or mixed boundary conditions, leads to the usual essential or natural boundary conditions. Essential boundary conditions: u = g for x c Ø ; so v = for x c Ø and the boundary integral vanishes. Natural boundary conditions: Øn + u = for x c Ø 2 ; so Ø 2 v u $ n ds = Ø 2 v Øn ds = Ø 2 uv ds Ø 2 v ds for x c Ø 2 Eric de Sturler /27/
10 Stokes theorem: S ( % u) $ n ds = u $ d, where n is the normal to the surface S, ds is the unit surface element, is the tangent to the boundary curve (counterclockwise), and d is the unit line element (on the boundary curve). Green s theorem can be derived from either Gauss or Stokes theorem: u $ v d = Ø v( u $ n) d vßu d = Ø v Øn d vßu d is easily derived from Gauss theorem. Using Stokes theorem and taking x 3 = ( S is part of the (x, x 2 ) plane): S 2 Øx Øx 2 dx dx 2 = (u 2 n u n 2 )d Eric de Sturler /27/
Numerical Methods for Partial Differential Equations
Numerical Methods for Partial Differential Equations Eric de Sturler University of Illinois at Urbana-Champaign Read section 8. to see where equations of type (au x ) x = f show up and their (exact) solution
More informationLegendre Transforms, Calculus of Varations, and Mechanics Principles
page 437 Appendix C Legendre Transforms, Calculus of Varations, and Mechanics Principles C.1 Legendre Transforms Legendre transforms map functions in a vector space to functions in the dual space. From
More informationPhysics 5153 Classical Mechanics. Canonical Transformations-1
1 Introduction Physics 5153 Classical Mechanics Canonical Transformations The choice of generalized coordinates used to describe a physical system is completely arbitrary, but the Lagrangian is invariant
More informationMassachusetts Institute of Technology Department of Physics. Final Examination December 17, 2004
Massachusetts Institute of Technology Department of Physics Course: 8.09 Classical Mechanics Term: Fall 004 Final Examination December 17, 004 Instructions Do not start until you are told to do so. Solve
More informationHamilton-Jacobi theory
Hamilton-Jacobi theory November 9, 04 We conclude with the crowning theorem of Hamiltonian dynamics: a proof that for any Hamiltonian dynamical system there exists a canonical transformation to a set of
More informationMath Tune-Up Louisiana State University August, Lectures on Partial Differential Equations and Hilbert Space
Math Tune-Up Louisiana State University August, 2008 Lectures on Partial Differential Equations and Hilbert Space 1. A linear partial differential equation of physics We begin by considering the simplest
More informationDuality, Dual Variational Principles
Duality, Dual Variational Principles April 5, 2013 Contents 1 Duality 1 1.1 Legendre and Young-Fenchel transforms.............. 1 1.2 Second conjugate and convexification................ 4 1.3 Hamiltonian
More informationPhysics 451/551 Theoretical Mechanics. G. A. Krafft Old Dominion University Jefferson Lab Lecture 18
Physics 451/551 Theoretical Mechanics G. A. Krafft Old Dominion University Jefferson Lab Lecture 18 Theoretical Mechanics Fall 018 Properties of Sound Sound Waves Requires medium for propagation Mainly
More informationChapter 12. Partial di erential equations Di erential operators in R n. The gradient and Jacobian. Divergence and rotation
Chapter 12 Partial di erential equations 12.1 Di erential operators in R n The gradient and Jacobian We recall the definition of the gradient of a scalar function f : R n! R, as @f grad f = rf =,..., @f
More informationComparison of Models for Finite Plasticity
Comparison of Models for Finite Plasticity A numerical study Patrizio Neff and Christian Wieners California Institute of Technology (Universität Darmstadt) Universität Augsburg (Universität Heidelberg)
More informationChapter 1: The Finite Element Method
Chapter 1: The Finite Element Method Michael Hanke Read: Strang, p 428 436 A Model Problem Mathematical Models, Analysis and Simulation, Part Applications: u = fx), < x < 1 u) = u1) = D) axial deformation
More informationHamiltonian. March 30, 2013
Hamiltonian March 3, 213 Contents 1 Variational problem as a constrained problem 1 1.1 Differential constaint......................... 1 1.2 Canonic form............................. 2 1.3 Hamiltonian..............................
More informationThe Particle-Field Hamiltonian
The Particle-Field Hamiltonian For a fundamental understanding of the interaction of a particle with the electromagnetic field we need to know the total energy of the system consisting of particle and
More informationVariational Principles
Part IB Variational Principles Year 218 217 216 215 214 213 212 211 21 218 Paper 1, Section I 4B 46 Variational Principles Find, using a Lagrange multiplier, the four stationary points in R 3 of the function
More informationPhysics 141 Energy 1 Page 1. Energy 1
Physics 4 Energy Page Energy What I tell you three times is true. Lewis Carroll The interplay of mathematics and physics The mathematization of physics in ancient times is attributed to the Pythagoreans,
More informationScientific Computing WS 2018/2019. Lecture 15. Jürgen Fuhrmann Lecture 15 Slide 1
Scientific Computing WS 2018/2019 Lecture 15 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 15 Slide 1 Lecture 15 Slide 2 Problems with strong formulation Writing the PDE with divergence and gradient
More informationGeometric Mechanics and Global Nonlinear Control for Multi-Body Dynamics
Geometric Mechanics and Global Nonlinear Control for Multi-Body Dynamics Harris McClamroch Aerospace Engineering, University of Michigan Joint work with Taeyoung Lee (George Washington University) Melvin
More informationVirtual Work and Variational Principles
Virtual Work and Principles Mathematically, the structural analysis problem is a boundary value problem (BVP). Forces, displacements, stresses, and strains are connected and computed within the framework
More informationCOMPLETE ALL ROUGH WORKINGS IN THE ANSWER BOOK AND CROSS THROUGH ANY WORK WHICH IS NOT TO BE ASSESSED.
BSc/MSci EXAMINATION PHY-304 Time Allowed: Physical Dynamics 2 hours 30 minutes Date: 28 th May 2009 Time: 10:00 Instructions: Answer ALL questions in section A. Answer ONLY TWO questions from section
More informationFALL TERM EXAM, PHYS 1211, INTRODUCTORY PHYSICS I Saturday, 14 December 2013, 1PM to 4 PM, AT 1003
FALL TERM EXAM, PHYS 1211, INTRODUCTORY PHYSICS I Saturday, 14 December 2013, 1PM to 4 PM, AT 1003 NAME: STUDENT ID: INSTRUCTION 1. This exam booklet has 14 pages. Make sure none are missing 2. There is
More informationREVIEW. Hamilton s principle. based on FW-18. Variational statement of mechanics: (for conservative forces) action Equivalent to Newton s laws!
Hamilton s principle Variational statement of mechanics: (for conservative forces) action Equivalent to Newton s laws! based on FW-18 REVIEW the particle takes the path that minimizes the integrated difference
More informationSimple Examples on Rectangular Domains
84 Chapter 5 Simple Examples on Rectangular Domains In this chapter we consider simple elliptic boundary value problems in rectangular domains in R 2 or R 3 ; our prototype example is the Poisson equation
More informationIntroduction to the Calculus of Variations
236861 Numerical Geometry of Images Tutorial 1 Introduction to the Calculus of Variations Alex Bronstein c 2005 1 Calculus Calculus of variations 1. Function Functional f : R n R Example: f(x, y) =x 2
More informationThe Simple Double Pendulum
The Simple Double Pendulum Austin Graf December 13, 2013 Abstract The double pendulum is a dynamic system that exhibits sensitive dependence upon initial conditions. This project explores the motion of
More informationVariational Formulations
Chapter 2 Variational Formulations In this chapter we will derive a variational (or weak) formulation of the elliptic boundary value problem (1.4). We will discuss all fundamental theoretical results that
More informationCurves in the configuration space Q or in the velocity phase space Ω satisfying the Euler-Lagrange (EL) equations,
Physics 6010, Fall 2010 Hamiltonian Formalism: Hamilton s equations. Conservation laws. Reduction. Poisson Brackets. Relevant Sections in Text: 8.1 8.3, 9.5 The Hamiltonian Formalism We now return to formal
More informationPoincaré (non-holonomic Lagrange) Equations
Department of Theoretical Physics Comenius University Bratislava fecko@fmph.uniba.sk Student Colloqium and School on Mathematical Physics, Stará Lesná, Slovakia, August 23-29, 2010 We will learn: In which
More informationVariational principles and Hamiltonian Mechanics
A Primer on Geometric Mechanics Variational principles and Hamiltonian Mechanics Alex L. Castro, PUC Rio de Janeiro Henry O. Jacobs, CMS, Caltech Christian Lessig, CMS, Caltech Alex L. Castro (PUC-Rio)
More informationScalar and vector fields
Scalar and vector fields What is a field in mathematics? Roughly speaking a field defines how a scalar-valued or vectorvalued quantity varies through space. We usually work with scalar and vector fields.
More informationA Short Essay on Variational Calculus
A Short Essay on Variational Calculus Keonwook Kang, Chris Weinberger and Wei Cai Department of Mechanical Engineering, Stanford University Stanford, CA 94305-4040 May 3, 2006 Contents 1 Definition of
More informationG : Statistical Mechanics
G25.2651: Statistical Mechanics Notes for Lecture 1 Defining statistical mechanics: Statistical Mechanics provies the connection between microscopic motion of individual atoms of matter and macroscopically
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationMultivariable Calculus
2 Multivariable Calculus 2.1 Limits and Continuity Problem 2.1.1 (Fa94) Let the function f : R n R n satisfy the following two conditions: (i) f (K ) is compact whenever K is a compact subset of R n. (ii)
More informationTHE CALCULUS OF VARIATIONS
THE CALCULUS OF VARIATIONS LAGRANGIAN MECHANICS AND OPTIMAL CONTROL AARON NELSON Abstract. The calculus of variations addresses the need to optimize certain quantities over sets of functions. We begin
More informationSuggested Solution to Assignment 7
MATH 422 (25-6) partial diferential equations Suggested Solution to Assignment 7 Exercise 7.. Suppose there exists one non-constant harmonic function u in, which attains its maximum M at x. Then by the
More informationSketchy Notes on Lagrangian and Hamiltonian Mechanics
Sketchy Notes on Lagrangian and Hamiltonian Mechanics Robert Jones Generalized Coordinates Suppose we have some physical system, like a free particle, a pendulum suspended from another pendulum, or a field
More informationVariational Principles for Equilibrium Physical Systems
Variational Principles for Equilibrium Physical Systems 1. Variational Principles One way of deriving the governing equations for a physical system is the express the relevant conservation statements and
More information[#1] R 3 bracket for the spherical pendulum
.. Holm Tuesday 11 January 2011 Solutions to MSc Enhanced Coursework for MA16 1 M3/4A16 MSc Enhanced Coursework arryl Holm Solutions Tuesday 11 January 2011 [#1] R 3 bracket for the spherical pendulum
More informationAdvanced Dynamics. - Lecture 4 Lagrange Equations. Paolo Tiso Spring Semester 2017 ETH Zürich
Advanced Dynamics - Lecture 4 Lagrange Equations Paolo Tiso Spring Semester 2017 ETH Zürich LECTURE OBJECTIVES 1. Derive the Lagrange equations of a system of particles; 2. Show that the equation of motion
More informationFunctions of Several Variables
Functions of Several Variables The Unconstrained Minimization Problem where In n dimensions the unconstrained problem is stated as f() x variables. minimize f()x x, is a scalar objective function of vector
More informationA very short introduction to the Finite Element Method
A very short introduction to the Finite Element Method Till Mathis Wagner Technical University of Munich JASS 2004, St Petersburg May 4, 2004 1 Introduction This is a short introduction to the finite element
More informationMath background. Physics. Simulation. Related phenomena. Frontiers in graphics. Rigid fluids
Fluid dynamics Math background Physics Simulation Related phenomena Frontiers in graphics Rigid fluids Fields Domain Ω R2 Scalar field f :Ω R Vector field f : Ω R2 Types of derivatives Derivatives measure
More informationA TASTE OF CALCULUS OF VARIATIONS
A TASTE OF CALCULUS OF VARIATIONS LIVIU I. NICOLAESCU ABSTRACT. These are notes for a talk at Arlo Caine s MPPM Lectures Series, University of Notre Dame. We survey through examples some of the basic problems,
More informationMATHEMATICAL PHYSICS
MATHEMATICAL PHYSICS Third Year SEMESTER 1 015 016 Classical Mechanics MP350 Prof. S. J. Hands, Prof. D. M. Heffernan, Dr. J.-I. Skullerud and Dr. M. Fremling Time allowed: 1 1 hours Answer two questions
More informationPhysics 451/551 Theoretical Mechanics. G. A. Krafft Old Dominion University Jefferson Lab Lecture 18
Physics 451/551 Theoretical Mechanics G. A. Krafft Old Dominion University Jefferson Lab Lecture 18 Theoretical Mechanics Fall 18 Properties of Sound Sound Waves Requires medium for propagation Mainly
More informationMath 342 Partial Differential Equations «Viktor Grigoryan
Math 342 Partial ifferential Equations «Viktor Grigoryan 3 Green s first identity Having studied Laplace s equation in regions with simple geometry, we now start developing some tools, which will lead
More informationfor changing independent variables. Most simply for a function f(x) the Legendre transformation f(x) B(s) takes the form B(s) = xs f(x) with s = df
Physics 106a, Caltech 1 November, 2018 Lecture 10: Hamiltonian Mechanics I The Hamiltonian In the Hamiltonian formulation of dynamics each second order ODE given by the Euler- Lagrange equation in terms
More information1. Let a(x) > 0, and assume that u and u h are the solutions of the Dirichlet problem:
Mathematics Chalmers & GU TMA37/MMG800: Partial Differential Equations, 011 08 4; kl 8.30-13.30. Telephone: Ida Säfström: 0703-088304 Calculators, formula notes and other subject related material are not
More informationQUALIFYING EXAMINATION, Part 1. Solutions. Problem 1: Mathematical Methods. r r r 2 r r2 = 0 r 2. d 3 r. t 0 e t dt = e t
QUALIFYING EXAMINATION, Part 1 Solutions Problem 1: Mathematical Methods (a) For r > we find 2 ( 1 r ) = 1 ( ) 1 r 2 r r2 = 1 ( 1 ) r r r 2 r r2 = r 2 However for r = we get 1 because of the factor in
More informationProblem set 2: Solutions Math 207B, Winter 2016
Problem set : Solutions Math 07B, Winter 016 1. A particle of mass m with position x(t) at time t has potential energy V ( x) an kinetic energy T = 1 m x t. The action of the particle over times t t 1
More informationLagrangian for Central Potentials
Physics 411 Lecture 2 Lagrangian for Central Potentials Lecture 2 Physics 411 Classical Mechanics II August 29th 2007 Here we will review the Lagrange formulation in preparation for the study of the central
More information06. Lagrangian Mechanics II
University of Rhode Island DigitalCommons@URI Classical Dynamics Physics Course Materials 2015 06. Lagrangian Mechanics II Gerhard Müller University of Rhode Island, gmuller@uri.edu Creative Commons License
More informationMathematical Analysis II, 2018/19 First semester
Mathematical Analysis II, 208/9 First semester Yoh Tanimoto Dipartimento di Matematica, Università di Roma Tor Vergata Via della Ricerca Scientifica, I-0033 Roma, Italy email: hoyt@mat.uniroma2.it We basically
More informationLecture Notes for PHY 405 Classical Mechanics
Lecture Notes for PHY 405 Classical Mechanics From Thorton & Marion s Classical Mechanics Prepared by Dr. Joseph M. Hahn Saint Mary s University Department of Astronomy & Physics September 1, 2005 Chapter
More information7.2 Calculus of Variations
7 CALCULUS OF VARIATIONS c 006 Gilbert Strang 7 Calculus of Variations One theme of this book is the relation of equations to minimum principles To minimize P is to solve P = 0 There may be more to it,
More informationPHY 5246: Theoretical Dynamics, Fall November 16 th, 2015 Assignment # 11, Solutions. p θ = L θ = mr2 θ, p φ = L θ = mr2 sin 2 θ φ.
PHY 5246: Theoretical Dynamics, Fall 215 November 16 th, 215 Assignment # 11, Solutions 1 Graded problems Problem 1 1.a) The Lagrangian is L = 1 2 m(ṙ2 +r 2 θ2 +r 2 sin 2 θ φ 2 ) V(r), (1) and the conjugate
More informationV (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)
IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common
More informationThe first order formalism and the transition to the
Problem 1. Hamiltonian The first order formalism and the transition to the This problem uses the notion of Lagrange multipliers and Legendre transforms to understand the action in the Hamiltonian formalism.
More informationPHYS 705: Classical Mechanics. Hamiltonian Formulation & Canonical Transformation
1 PHYS 705: Classical Mechanics Hamiltonian Formulation & Canonical Transformation Legendre Transform Let consider the simple case with ust a real value function: F x F x expresses a relationship between
More informationNoether s Theorem. 4.1 Ignorable Coordinates
4 Noether s Theorem 4.1 Ignorable Coordinates A central recurring theme in mathematical physics is the connection between symmetries and conservation laws, in particular the connection between the symmetries
More informationMath 265H: Calculus III Practice Midterm II: Fall 2014
Name: Section #: Math 65H: alculus III Practice Midterm II: Fall 14 Instructions: This exam has 7 problems. The number of points awarded for each question is indicated in the problem. Answer each question
More information4.1 Analysis of functions I: Increase, decrease and concavity
4.1 Analysis of functions I: Increase, decrease and concavity Definition Let f be defined on an interval and let x 1 and x 2 denote points in that interval. a) f is said to be increasing on the interval
More informationA Very General Hamilton Jacobi Theorem
University of Salerno 9th AIMS ICDSDEA Orlando (Florida) USA, July 1 5, 2012 Introduction: Generalized Hamilton-Jacobi Theorem Hamilton-Jacobi (HJ) Theorem: Consider a Hamiltonian system with Hamiltonian
More informationName: Instructor: Lecture time: TA: Section time:
Math 222 Final May 11, 29 Name: Instructor: Lecture time: TA: Section time: INSTRUCTIONS READ THIS NOW This test has 1 problems on 16 pages worth a total of 2 points. Look over your test package right
More informationDynamics and Stability application to submerged bodies, vortex streets and vortex-body systems
Dynamics and Stability application to submerged bodies, vortex streets and vortex-body systems Eva Kanso University of Southern California CDS 140B Introduction to Dynamics February 5 and 7, 2008 Fish
More informationC2 Differential Equations : Computational Modeling and Simulation Instructor: Linwei Wang
C2 Differential Equations 4040-849-03: Computational Modeling and Simulation Instructor: Linwei Wang Part II Variational Principle Calculus Revisited Partial Derivatives Function of one variable df dx
More informationPhysical Dynamics (PHY-304)
Physical Dynamics (PHY-304) Gabriele Travaglini March 31, 2012 1 Review of Newtonian Mechanics 1.1 One particle Lectures 1-2. Frame, velocity, acceleration, number of degrees of freedom, generalised coordinates.
More informationMATH 332: Vector Analysis Summer 2005 Homework
MATH 332, (Vector Analysis), Summer 2005: Homework 1 Instructor: Ivan Avramidi MATH 332: Vector Analysis Summer 2005 Homework Set 1. (Scalar Product, Equation of a Plane, Vector Product) Sections: 1.9,
More informationContents. MATH 32B-2 (18W) (L) G. Liu / (TA) A. Zhou Calculus of Several Variables. 1 Multiple Integrals 3. 2 Vector Fields 9
MATH 32B-2 (8W) (L) G. Liu / (TA) A. Zhou Calculus of Several Variables Contents Multiple Integrals 3 2 Vector Fields 9 3 Line and Surface Integrals 5 4 The Classical Integral Theorems 9 MATH 32B-2 (8W)
More informationAnalytical Dynamics: Lagrange s Equation and its Application A Brief Introduction
Analytical Dynamics: Lagrange s Equation and its Application A Brief Introduction D. S. Stutts, Ph.D. Associate Professor of Mechanical Engineering Missouri University of Science and Technology Rolla,
More informationChapter 1. Principles of Motion in Invariantive Mechanics
Chapter 1 Principles of Motion in Invariantive Mechanics 1.1. The Euler-Lagrange and Hamilton s equations obtained by means of exterior forms Let L = L(q 1,q 2,...,q n, q 1, q 2,..., q n,t) L(q, q,t) (1.1)
More informationLecture Notes for PHY 405 Classical Mechanics
Lecture Notes for PHY 45 Classical Mechanics From Thorton & Marion s Classical Mechanics Prepared by Dr. Joseph M. Hahn Saint Mary s University Department of Astronomy & Physics October 11, 25 Chapter
More informationMarion and Thornton. Tyler Shendruk October 1, Hamilton s Principle - Lagrangian and Hamiltonian dynamics.
Marion and Thornton Tyler Shendruk October 1, 2010 1 Marion and Thornton Chapter 7 Hamilton s Principle - Lagrangian and Hamiltonian dynamics. 1.1 Problem 6.4 s r z θ Figure 1: Geodesic on circular cylinder
More informationHamiltonian Mechanics
Alain J. Brizard Saint Michael's College Hamiltonian Mechanics 1 Hamiltonian The k second-order Euler-Lagrange equations on con guration space q =(q 1 ; :::; q k ) d @ _q j = @q j ; (1) can be written
More informationMathematical Methods of Physics I ChaosBook.org/ predrag/courses/phys Homework 1
PHYS 6124 Handout 6 23 August 2012 Mathematical Methods of Physics I ChaosBook.org/ predrag/courses/phys-6124-12 Homework 1 Prof. P. Goldbart School of Physics Georgia Tech Homework assignments are posted
More informationClassical mechanics of particles and fields
Classical mechanics of particles and fields D.V. Skryabin Department of Physics, University of Bath PACS numbers: The concise and transparent exposition of many topics covered in this unit can be found
More informationWe first review various rules for easy differentiation of common functions: The same procedure works for a larger number of terms.
1 Math 182 Lecture Notes 1. Review of Differentiation To differentiate a function y = f(x) is to find its derivative f '(x). Another standard notation for the derivative is Dx(f(x)). Recall the meanings
More informationl1, l2, l3, ln l1 + l2 + l3 + ln
Work done by a constant force: Consider an object undergoes a displacement S along a straight line while acted on a force F that makes an angle θ with S as shown The work done W by the agent is the product
More informationThe Principle of Least Action
The Principle of Least Action Anders Svensson Abstract In this paper, the principle of least action in classical mechanics is studied. The term is used in several different contexts, mainly for Hamilton
More information15. Hamiltonian Mechanics
University of Rhode Island DigitalCommons@URI Classical Dynamics Physics Course Materials 2015 15. Hamiltonian Mechanics Gerhard Müller University of Rhode Island, gmuller@uri.edu Creative Commons License
More informationCalculus of Variations Summer Term 2015
Calculus of Variations Summer Term 2015 Lecture 12 Universität des Saarlandes 17. Juni 2015 c Daria Apushkinskaya (UdS) Calculus of variations lecture 12 17. Juni 2015 1 / 31 Purpose of Lesson Purpose
More information2 A Model, Harmonic Map, Problem
ELLIPTIC SYSTEMS JOHN E. HUTCHINSON Department of Mathematics School of Mathematical Sciences, A.N.U. 1 Introduction Elliptic equations model the behaviour of scalar quantities u, such as temperature or
More informationSOLUTIONS TO THE FINAL EXAM. December 14, 2010, 9:00am-12:00 (3 hours)
SOLUTIONS TO THE 18.02 FINAL EXAM BJORN POONEN December 14, 2010, 9:00am-12:00 (3 hours) 1) For each of (a)-(e) below: If the statement is true, write TRUE. If the statement is false, write FALSE. (Please
More informationGeneralized Coordinates, Lagrangians
Generalized Coordinates, Lagrangians Sourendu Gupta TIFR, Mumbai, India Classical Mechanics 2012 August 10, 2012 Generalized coordinates Consider again the motion of a simple pendulum. Since it is one
More informationAMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends
AMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends Lecture 3: Finite Elements in 2-D Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Finite Element Methods 1 / 18 Outline 1 Boundary
More information8 Rotational motion of solid objects
8 Rotational motion of solid objects Kinematics of rotations PHY166 Fall 005 In this Lecture we call solid objects such extended objects that are rigid (nondeformable) and thus retain their shape. In contrast
More informationAnnouncements. Topics: Homework:
Announcements Topics: - sections 7.3 (the definite integral +area), 7.4 (FTC), 7.5 (additional techniques of integration) * Read these sections and study solved examples in your textbook! Homework: - review
More informationJohann Bernoulli ( )he established the principle of the so-called virtual work for static systems.
HISTORICAL SURVEY Wilhelm von Leibniz (1646-1716): Leibniz s approach to mechanics was based on the use of mathematical operations with the scalar quantities of energy, as opposed to the vector quantities
More informationRobotics. Islam S. M. Khalil. November 15, German University in Cairo
Robotics German University in Cairo November 15, 2016 Fundamental concepts In optimal control problems the objective is to determine a function that minimizes a specified functional, i.e., the performance
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Sobolev Embedding Theorems Embedding Operators and the Sobolev Embedding Theorem
More informationClassical Mechanics in Hamiltonian Form
Classical Mechanics in Hamiltonian Form We consider a point particle of mass m, position x moving in a potential V (x). It moves according to Newton s law, mẍ + V (x) = 0 (1) This is the usual and simplest
More informationLecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.
Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference
More informationg(x,y) = c. For instance (see Figure 1 on the right), consider the optimization problem maximize subject to
1 of 11 11/29/2010 10:39 AM From Wikipedia, the free encyclopedia In mathematical optimization, the method of Lagrange multipliers (named after Joseph Louis Lagrange) provides a strategy for finding the
More informationGetting Some Big Air
Getting Some Big Air Control #10499 February 14, 2011 Abstract In this paper we address the problem of optimizing a ramp for a snowboarder to travel. Our approach is two sided. We first address the forward
More informationCalculus of Variations
ECE 68 Midterm Exam Solution April 1, 8 1 Calculus of Variations This exam is open book and open notes You may consult additional references You may even discuss the problems (with anyone), but you must
More informationPRACTICE PROBLEMS. Please let me know if you find any mistakes in the text so that i can fix them. 1. Mixed partial derivatives.
PRACTICE PROBLEMS Please let me know if you find any mistakes in the text so that i can fix them. 1.1. Let Show that f is C 1 and yet How is that possible? 1. Mixed partial derivatives f(x, y) = {xy x
More informationFigure 25:Differentials of surface.
2.5. Change of variables and Jacobians In the previous example we saw that, once we have identified the type of coordinates which is best to use for solving a particular problem, the next step is to do
More informationChapter Two: Numerical Methods for Elliptic PDEs. 1 Finite Difference Methods for Elliptic PDEs
Chapter Two: Numerical Methods for Elliptic PDEs Finite Difference Methods for Elliptic PDEs.. Finite difference scheme. We consider a simple example u := subject to Dirichlet boundary conditions ( ) u
More informationWRT in 2D: Poisson Example
WRT in 2D: Poisson Example Consider 2 u f on [, L x [, L y with u. WRT: For all v X N, find u X N a(v, u) such that v u dv v f dv. Follows from strong form plus integration by parts: ( ) 2 u v + 2 u dx
More informationCIV-E1060 Engineering Computation and Simulation Examination, December 12, 2017 / Niiranen
CIV-E16 Engineering Computation and Simulation Examination, December 12, 217 / Niiranen This examination consists of 3 problems rated by the standard scale 1...6. Problem 1 Let us consider a long and tall
More information