2.13 Variation and Linearisation of Kinematic Tensors

Size: px
Start display at page:

Download "2.13 Variation and Linearisation of Kinematic Tensors"

Transcription

1 Section.3.3 Variation an Linearisation of Kinematic ensors.3. he Variation of Kinematic ensors he Variation In this section is reviewe the concept of the variation, introce in Part I, 8.5. he variation is efine as follows: consier a fnction (, with ( a secon fnction which is at most infinitesimally ifferent from ( at every point x, Fig..3. ( ( ( Figre.3.: the variation x hen efine ( ( he Variation (.3.) he operator is calle the variation symbol an is calle the variation of (. he variation of ( is nerstoo to represent an infinitesimal change in the fnction at x. Note from the figre that a variation of a fnction is ifferent to a ifferential. he orinary ifferentiation gives a measre of the change of a fnction reslting from a specifie change in the inepenent variable (in this case. Also, note that the inepenent variable oes not participate in the variation process; the variation operator imparts an infinitesimal change to the fnction at some fixe x formally, one can write this as x 0. he Commtative Properties of the variation operator () (.3.) 37

2 Section.3 Proof: ( ) ( x x () ( ( (.3.3) Proof: x ( x ( x ( x x ( ( ( Variation of a Fnction Consier A, a scalar-, vector-, or tensor-vale fnction of, variation to,, A changes to A A. When we apply a. he variation of A is then efine as A (, ) A ( ) A ( ) (.3.4) (in the limit as 0 ). his can be expresse sing the concept of the irectional A, so that erivative in the sal way (see.6.): consier the fnction A A an 0 A A0 A/, or 0 Setting then gives Eqn..3.4; ths A A. A aylor expansion gives A A A (.3.5) 0 A( ) A( ) A[ ] (.3.6) where A[ ] is the irectional erivative of A in the irection ; the irectional erivative in this context is the variation of A: A, A[ ] A 0 (.3.7) 38

3 Section.3 For example, consier the scalar fnction P : E, where P an E are secon orer tensors. hen E, E E[ E] P : E E P : E 0 (.3.8) he secon variation is efine as A A A[ ] A 0 (.3.9) For example, for a scalar fnction of a vector, the chain rle an Eqn..3. give, 0 0 (.3.0) Variation of Fnctions of the Displacement In what follows is iscsse the change (variation) in fnctions A () when the isplacement (or velocity) fiels nergo a variation. hese ieas are sefl in formlating variational principles of mechanics (see, for example, 3.9). Shown in Fig..3. is the crrent configration frozen at some instant in time. he isplacement fiel is then allowe to nergo a variation. his change to the isplacement fiel eviently changes kinematic tensors, an these changes are now investigate. Note that this variation to the isplacement inces a variation to x, x, bt X remains nchange, X 0. X reference configration ( x crrent configration Figre.3.: a variation of the isplacement 39

4 Section.3 o evalate the variation of the eformation graient F, F x X an Eqn...43, F Gra I isplacement fiel, note that the efinition.3.7,,, where is the. One has, from F F F, 0 0 Gra F I Gra (.3.) Noting the first commtative property of the variation,.3., this can also be expresse as, Gra Note that is completely inepenent of the fnction. F (.3.) Here are some other examples, involving the inverse eformation graient, the Green- Lagrange strain, the inverse right Cachy-Green strain an the spatial line element: { Problem -3} F C F E F εf F gra εf (.3.3) where ε is the small strain tensor, Eqn One also has, sing the chain rle for the irectional erivative, Eqn..5.8, the irectional erivative for the eterminant, Eqn..5.3, the trace relation.0.0e, Eqn...8b, J, et F, et F F et F F et FGra F et FF : Gra (.3.4) JtrGra F Jtrgra J iv 330

5 Section.3 Example o pt some of the above concepts into a simple an less abstract setting, consier the following scenario: a bar over 0 X is extene, as illstrate in Fig..3.3, accoring to: x X 3 X x 3 (.3.5) he eformation graient is F Gra 4X (.3.6) So, for example, in the initial configration (A), an infinitesimal line element at X 0 oes not stretch ( F 0 ) whereas a line element at X stretches by 4. he inverse eformation graient is F grax 8 x 3 (.3.7) his implies that, in the crrent configration (B), an infinitesimal line element at x 3 is the same size as its conterpart in the initial configration ( F 0 ) whereas a line element at x 5 shrinks by a factor of 4 when retrning to the initial configration F ( ) 4 X A A B X F 4 C C 4 X Figre.3.3: a motion an a variation 33

6 Section.3 Now introce a variation, which moves the bar from configration B to configration C: x X 3 (.3.8) he point at 3 moves to 3 3 an the point at 5 moves to 5 5. (his variation happens to be a simple linear fnction of x, bt it can be anything for or prposes here.) he variation is plotte below as a fnction of X an x. 5 3 x 5 3 Figre.3.4: the variation as a fnction of x an X Differentiating Eqns..3.9, the graients of the variations are X X Gra 4 gra (.3.0) which are the slopes in Figre.3.4. o calclate the F associate with the new variation configration, i.e. F points X have now move to: an so 3 3, note that X X (.3.) 33

7 Section.3 F Gra X 3 4X4X (.3.) his says that an infinitesimal line element at X 0 oes not stretch when moving to configration C ( F 0 ) whereas a line element at X stretches by 4 4. Sbtracting Eqn..3.7 form Eqn..3.: 4 From Eqn..3.0, this verifies Eqn..3., that F F F X (.3.3) F Gra (.3.4) We col also calclate the variation of F by moving irectly from configration B to configration C. he movement of the particles from B to C is given by Eqn..3.9: 3 F Gra X 3 4X. X an so, base on this motion, o calclate the that the new crrent position x is (Eqn..3.): F associate with the new variation configration, i.e. F, note x X X C 3 3 xc X 3 (.3.5) his means that the point 3 3 in configration C correspons to X 0 an the point 5 5 correspons to X. hen, F xc 3 C 8 3 xc (.3.6) So an element at the point 3 3 in configration C oes not change in size as it is mappe back to the initial configration, whereas an element at the point 5 5 shrinks back to the / 4 4, as inicate in Fig initial configration by a factor of Alternatively, since xc xx, this can be written as F 8 x 3 (.3.7) 333

8 Section.3 Sbtracting Eqn..3.8 from Eqn..3.7, the variation of the inverse eformation graient is then F F F x 3 x x (.3.8) Using a series expansion, ),, for small (neglecting terms of orer F 8 x 3 (.3.9) From Eqns..3.8 an.3.0, this verifies the relation.3.3: A formla for the inverse eformation graient is F I/ x, bt that F F gra δ (.3.30) F I x. F Igra. However, note that / C he Lie Variation he Lie-variation is efine for spatial vectors an tensors as a variation holing the eforme basis constant. For example, i j b L a aijg g (.3.3) he object is first plle-back, the variation is then taken an finally a psh-forwar is carrie ot. For example, analogos to..66, a a (.3.3) L, For example, consier the Lie-variation of the Eler-Almansi strain e. First, from..56b, b b e E e E F. From..40a,. hen.3.3b gives εf b e b b F εf ε L e, (.3.33) 334

9 Section.3.3. Linearisation of Kinematic Fnctions Linearisation of a Fnction As for the variation, consier A, a scalar-, vector-, or tensor-vale fnction of. If nergoes an increment, then, analogos to.3.4, A A[ ] he irectional erivative A[ ] linearization of A with respect to is efine to be A (.3.34) in this context is also enote by A, A A,,. he L A (.3.35) Using exactly the same metho of calclation as was se for the variations above, the linearization of F an E, for example, are L F L E where ε gra gra, F F F Gra, E E E F εf is the linearise small strain tensor ε. (.3.36) Linearisation of Variations of a Fnction One can also linearise the variation of a fnction. For example,, A, A, he secon term here is the irectional erivative his leas to an expression similar to vector, L A (.3.37) A[, ] A[ ] A (.3.38) A 0. For example, for a scalar fnction of a (.3.39) 335

10 Section.3 Consier now the virtal Green-Lagrange strain,.3.b, E F εf. o carry ot the linearization of E, it is convenient to first write it in the form E F εf F gra Gra gra F F F Gra (.3.40) hen E E Gra F F Gra (.3.4) Recall that the variation is inepenent of ; this eqation is being linearise with respect to, an is naffecte by the linearization (see Fig..3.3 below). However, the motion, an in particlar F, are affecte by the increment in. hs { Problem 4} E sym Gra Gra (.3.4) reference configration crrent configration Figre.3.3: linearisation As with the variational operator, one can efine the linearization of a spatial tensor as involving a pll back, followe by the irectional erivative, an finally the psh forwar operation. hs, a a (.3.43).3.3 Problems 336

11 Section.3. Use Eqn..., E F F I GravF, Eqn..3., F Gra,, an Eqn...8b, grav, to show that E F εf, where ε is the small strain tensor, Eqn Use.3. to show that the variation of the inverse eformation graient F is F F gra. [Hint: ifferente the relation F F I by the proct rle an then se the relation vf grav Gra for vector v.] 3. Use the efinition C F F to show that C F εf. sym A A A to show that 4. Use the relation Gra F F Gra sym Gra E Gra 5. Use e ε gra gra to show that the e e sym Gra sym gra Gra gra 337

THE DISPLACEMENT GRADIENT AND THE LAGRANGIAN STRAIN TENSOR Revision B

THE DISPLACEMENT GRADIENT AND THE LAGRANGIAN STRAIN TENSOR Revision B HE DISPLACEMEN GRADIEN AND HE LAGRANGIAN SRAIN ENSOR Revision B By om Irvine Email: tom@irvinemail.org Febrary, 05 Displacement Graient Sppose a boy having a particlar configration at some reference time

More information

Logarithmic, Exponential and Other Transcendental Functions

Logarithmic, Exponential and Other Transcendental Functions Logarithmic, Eponential an Other Transcenental Fnctions 5: The Natral Logarithmic Fnction: Differentiation The Definition First, yo mst know the real efinition of the natral logarithm: ln= t (where > 0)

More information

Theorem (Change of Variables Theorem):

Theorem (Change of Variables Theorem): Avance Higher Notes (Unit ) Prereqisites: Integrating (a + b) n, sin (a + b) an cos (a + b); erivatives of tan, sec, cosec, cot, e an ln ; sm/ifference rles; areas ner an between crves. Maths Applications:

More information

n s n Z 0 on complex-valued functions on the circle. Then sgn n + 1 ) n + 1 2) s

n s n Z 0 on complex-valued functions on the circle. Then sgn n + 1 ) n + 1 2) s . What is the eta invariant? The eta invariant was introce in the famos paper of Atiyah, Patoi, an Singer see [], in orer to proce an inex theorem for manifols with bonary. The eta invariant of a linear

More information

II. First variation of functionals

II. First variation of functionals II. First variation of functionals The erivative of a function being zero is a necessary conition for the etremum of that function in orinary calculus. Let us now tackle the question of the equivalent

More information

Euler equations for multiple integrals

Euler equations for multiple integrals Euler equations for multiple integrals January 22, 2013 Contents 1 Reminer of multivariable calculus 2 1.1 Vector ifferentiation......................... 2 1.2 Matrix ifferentiation........................

More information

The total derivative. Chapter Lagrangian and Eulerian approaches

The total derivative. Chapter Lagrangian and Eulerian approaches Chapter 5 The total erivative 51 Lagrangian an Eulerian approaches The representation of a flui through scalar or vector fiels means that each physical quantity uner consieration is escribe as a function

More information

Chapter 2 Lagrangian Modeling

Chapter 2 Lagrangian Modeling Chapter 2 Lagrangian Moeling The basic laws of physics are use to moel every system whether it is electrical, mechanical, hyraulic, or any other energy omain. In mechanics, Newton s laws of motion provie

More information

Math Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors

Math Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors Math 18.02 Notes on ifferentials, the Chain Rule, graients, irectional erivative, an normal vectors Tangent plane an linear approximation We efine the partial erivatives of f( xy, ) as follows: f f( x+

More information

Implicit Differentiation

Implicit Differentiation Implicit Differentiation Thus far, the functions we have been concerne with have been efine explicitly. A function is efine explicitly if the output is given irectly in terms of the input. For instance,

More information

UNIT IV BOUNDARY LAYER AND FLOW THROUGH PIPES Definition of bonary layer Thickness an classification Displacement an momentm Thickness Development of laminar an trblent flows in circlar pipes Major an

More information

Calculus of Variations

Calculus of Variations Calculus of Variations Lagrangian formalism is the main tool of theoretical classical mechanics. Calculus of Variations is a part of Mathematics which Lagrangian formalism is base on. In this section,

More information

Topic 2.3: The Geometry of Derivatives of Vector Functions

Topic 2.3: The Geometry of Derivatives of Vector Functions BSU Math 275 Notes Topic 2.3: The Geometry of Derivatives of Vector Functions Textbook Sections: 13.2 From the Toolbox (what you nee from previous classes): Be able to compute erivatives scalar-value functions

More information

Desert Mountain H. S. Math Department Summer Work Packet

Desert Mountain H. S. Math Department Summer Work Packet Corse #50-51 Desert Montain H. S. Math Department Smmer Work Packet Honors/AP/IB level math corses at Desert Montain are for stents who are enthsiastic learners of mathematics an whose work ethic is of

More information

Derivative of a Constant Multiple of a Function Theorem: If f is a differentiable function and if c is a constant, then

Derivative of a Constant Multiple of a Function Theorem: If f is a differentiable function and if c is a constant, then Bob Brown Math 51 Calculus 1 Chapter 3, Section Complete 1 Review of the Limit Definition of the Derivative Write the it efinition of the erivative function: f () Derivative of a Constant Multiple of a

More information

Applications of the Wronskian to ordinary linear differential equations

Applications of the Wronskian to ordinary linear differential equations Physics 116C Fall 2011 Applications of the Wronskian to orinary linear ifferential equations Consier a of n continuous functions y i (x) [i = 1,2,3,...,n], each of which is ifferentiable at least n times.

More information

Table of Common Derivatives By David Abraham

Table of Common Derivatives By David Abraham Prouct an Quotient Rules: Table of Common Derivatives By Davi Abraham [ f ( g( ] = [ f ( ] g( + f ( [ g( ] f ( = g( [ f ( ] g( g( f ( [ g( ] Trigonometric Functions: sin( = cos( cos( = sin( tan( = sec

More information

Calculus in the AP Physics C Course The Derivative

Calculus in the AP Physics C Course The Derivative Limits an Derivatives Calculus in the AP Physics C Course The Derivative In physics, the ieas of the rate change of a quantity (along with the slope of a tangent line) an the area uner a curve are essential.

More information

HOMEWORK 2 SOLUTIONS

HOMEWORK 2 SOLUTIONS HOMEWORK 2 SOLUTIONS PHIL SAAD 1. Carroll 1.4 1.1. A qasar, a istance D from an observer on Earth, emits a jet of gas at a spee v an an angle θ from the line of sight of the observer. The apparent spee

More information

Optimal Operation by Controlling the Gradient to Zero

Optimal Operation by Controlling the Gradient to Zero Optimal Operation by Controlling the Graient to Zero Johannes Jäschke Sigr Skogesta Department of Chemical Engineering, Norwegian University of Science an Technology, NTNU, Tronheim, Norway (e-mail: {jaschke}{skoge}@chemeng.ntn.no)

More information

Mehmet Pakdemirli* Precession of a Planet with the Multiple Scales Lindstedt Poincare Technique (2)

Mehmet Pakdemirli* Precession of a Planet with the Multiple Scales Lindstedt Poincare Technique (2) Z. Natrforsch. 05; aop Mehmet Pakemirli* Precession of a Planet with the Mltiple Scales Linstet Poincare Techniqe DOI 0.55/zna-05-03 Receive May, 05; accepte Jly 5, 05 Abstract: The recently evelope mltiple

More information

JUST THE MATHS UNIT NUMBER DIFFERENTIATION 2 (Rates of change) A.J.Hobson

JUST THE MATHS UNIT NUMBER DIFFERENTIATION 2 (Rates of change) A.J.Hobson JUST THE MATHS UNIT NUMBER 10.2 DIFFERENTIATION 2 (Rates of change) by A.J.Hobson 10.2.1 Introuction 10.2.2 Average rates of change 10.2.3 Instantaneous rates of change 10.2.4 Derivatives 10.2.5 Exercises

More information

Calculus of Variations

Calculus of Variations 16.323 Lecture 5 Calculus of Variations Calculus of Variations Most books cover this material well, but Kirk Chapter 4 oes a particularly nice job. x(t) x* x*+ αδx (1) x*- αδx (1) αδx (1) αδx (1) t f t

More information

Conservation Laws. Chapter Conservation of Energy

Conservation Laws. Chapter Conservation of Energy 20 Chapter 3 Conservation Laws In orer to check the physical consistency of the above set of equations governing Maxwell-Lorentz electroynamics [(2.10) an (2.12) or (1.65) an (1.68)], we examine the action

More information

Modeling of a Self-Oscillating Cantilever

Modeling of a Self-Oscillating Cantilever Moeling of a Self-Oscillating Cantilever James Blanchar, Hi Li, Amit Lal, an Doglass Henerson University of Wisconsin-Maison 15 Engineering Drive Maison, Wisconsin 576 Abstract A raioisotope-powere, self-oscillating

More information

Math Implicit Differentiation. We have discovered (and proved) formulas for finding derivatives of functions like

Math Implicit Differentiation. We have discovered (and proved) formulas for finding derivatives of functions like Math 400 3.5 Implicit Differentiation Name We have iscovere (an prove) formulas for fining erivatives of functions like f x x 3x 4x. 3 This amounts to fining y for 3 y x 3x 4x. Notice that in this case,

More information

Lie symmetry and Mei conservation law of continuum system

Lie symmetry and Mei conservation law of continuum system Chin. Phys. B Vol. 20, No. 2 20 020 Lie symmetry an Mei conservation law of continuum system Shi Shen-Yang an Fu Jing-Li Department of Physics, Zhejiang Sci-Tech University, Hangzhou 3008, China Receive

More information

2.20 Marine Hydrodynamics Lecture 3

2.20 Marine Hydrodynamics Lecture 3 2.20 Marine Hyroynamics, Fall 2018 Lecture 3 Copyright c 2018 MIT - Department of Mechanical Engineering, All rights reserve. 1.7 Stress Tensor 2.20 Marine Hyroynamics Lecture 3 1.7.1 Stress Tensor τ ij

More information

Lecture XII. where Φ is called the potential function. Let us introduce spherical coordinates defined through the relations

Lecture XII. where Φ is called the potential function. Let us introduce spherical coordinates defined through the relations Lecture XII Abstract We introuce the Laplace equation in spherical coorinates an apply the metho of separation of variables to solve it. This will generate three linear orinary secon orer ifferential equations:

More information

Lecture 1b. Differential operators and orthogonal coordinates. Partial derivatives. Divergence and divergence theorem. Gradient. A y. + A y y dy. 1b.

Lecture 1b. Differential operators and orthogonal coordinates. Partial derivatives. Divergence and divergence theorem. Gradient. A y. + A y y dy. 1b. b. Partial erivatives Lecture b Differential operators an orthogonal coorinates Recall from our calculus courses that the erivative of a function can be efine as f ()=lim 0 or using the central ifference

More information

Lecture 2 Lagrangian formulation of classical mechanics Mechanics

Lecture 2 Lagrangian formulation of classical mechanics Mechanics Lecture Lagrangian formulation of classical mechanics 70.00 Mechanics Principle of stationary action MATH-GA To specify a motion uniquely in classical mechanics, it suffices to give, at some time t 0,

More information

Math 11 Fall 2016 Section 1 Monday, September 19, Definition: A vector parametric equation for the line parallel to vector v = x v, y v, z v

Math 11 Fall 2016 Section 1 Monday, September 19, Definition: A vector parametric equation for the line parallel to vector v = x v, y v, z v Math Fall 06 Section Monay, September 9, 06 First, some important points from the last class: Definition: A vector parametric equation for the line parallel to vector v = x v, y v, z v passing through

More information

Lectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs

Lectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs Lectures - Week 10 Introuction to Orinary Differential Equations (ODES) First Orer Linear ODEs When stuying ODEs we are consiering functions of one inepenent variable, e.g., f(x), where x is the inepenent

More information

Basic Thermoelasticity

Basic Thermoelasticity Basic hermoelasticity Biswajit Banerjee November 15, 2006 Contents 1 Governing Equations 1 1.1 Balance Laws.............................................. 2 1.2 he Clausius-Duhem Inequality....................................

More information

Rainer Friedrich

Rainer Friedrich Rainer Frierich et al Rainer Frierich rfrierich@lrztme Holger Foysi Joern Sesterhenn FG Stroemngsmechanik Technical University Menchen Boltzmannstr 5 D-85748 Garching, Germany Trblent Momentm an Passive

More information

The derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x)

The derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x) Y. D. Chong (2016) MH2801: Complex Methos for the Sciences 1. Derivatives The erivative of a function f(x) is another function, efine in terms of a limiting expression: f (x) f (x) lim x δx 0 f(x + δx)

More information

G j dq i + G j. q i. = a jt. and

G j dq i + G j. q i. = a jt. and Lagrange Multipliers Wenesay, 8 September 011 Sometimes it is convenient to use reunant coorinates, an to effect the variation of the action consistent with the constraints via the metho of Lagrange unetermine

More information

Integration Review. May 11, 2013

Integration Review. May 11, 2013 Integration Review May 11, 2013 Goals: Review the funamental theorem of calculus. Review u-substitution. Review integration by parts. Do lots of integration eamples. 1 Funamental Theorem of Calculus In

More information

5.5 U-substitution. Solution. Z

5.5 U-substitution. Solution. Z CHAPTER 5. THE DEFINITE INTEGRAL 22 5.5 U-sbstittion Eample. (a) Fin the erivative of sin( 2 ). (b) Fin the anti-erivative cos( 2 ). Soltion. (a) We se the chain rle: sin(2 )=cos( 2 )( 2 ) 0 =cos( 2 )2

More information

Switching Time Optimization in Discretized Hybrid Dynamical Systems

Switching Time Optimization in Discretized Hybrid Dynamical Systems Switching Time Optimization in Discretize Hybri Dynamical Systems Kathrin Flaßkamp, To Murphey, an Sina Ober-Blöbaum Abstract Switching time optimization (STO) arises in systems that have a finite set

More information

Further Differentiation and Applications

Further Differentiation and Applications Avance Higher Notes (Unit ) Prerequisites: Inverse function property; prouct, quotient an chain rules; inflexion points. Maths Applications: Concavity; ifferentiability. Real-Worl Applications: Particle

More information

SEG Houston 2009 International Exposition and Annual Meeting

SEG Houston 2009 International Exposition and Annual Meeting Fonations of the metho of M fiel separation into pgoing an ongoing parts an its application to MCSM ata Michael S. Zhanov an Shming Wang*, Universit of Utah Smmar The renee interest in the methos of electromagnetic

More information

Math 273b: Calculus of Variations

Math 273b: Calculus of Variations Math 273b: Calcls of Variations Yacob Kreh Homework #3 [1] Consier the 1D length fnctional minimization problem min F 1 1 L, or min 1 + 2, for twice ifferentiable fnctions : [, 1] R with bonary conitions,

More information

APPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France

APPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France APPROXIMAE SOLUION FOR RANSIEN HEA RANSFER IN SAIC URBULEN HE II B. Bauouy CEA/Saclay, DSM/DAPNIA/SCM 91191 Gif-sur-Yvette Ceex, France ABSRAC Analytical solution in one imension of the heat iffusion equation

More information

Lagrangian and Hamiltonian Mechanics

Lagrangian and Hamiltonian Mechanics Lagrangian an Hamiltonian Mechanics.G. Simpson, Ph.. epartment of Physical Sciences an Engineering Prince George s Community College ecember 5, 007 Introuction In this course we have been stuying classical

More information

θ x = f ( x,t) could be written as

θ x = f ( x,t) could be written as 9. Higher orer PDEs as systems of first-orer PDEs. Hyperbolic systems. For PDEs, as for ODEs, we may reuce the orer by efining new epenent variables. For example, in the case of the wave equation, (1)

More information

1 Lecture 13: The derivative as a function.

1 Lecture 13: The derivative as a function. 1 Lecture 13: Te erivative as a function. 1.1 Outline Definition of te erivative as a function. efinitions of ifferentiability. Power rule, erivative te exponential function Derivative of a sum an a multiple

More information

18-660: Numerical Methods for Engineering Design and Optimization

18-660: Numerical Methods for Engineering Design and Optimization 8-66: Nmerical Methos for Engineering Design an Optimization in Li Department of ECE Carnegie Mellon University Pittsbrgh, PA 53 Slie Overview Geometric Problems Maximm inscribe ellipsoi Minimm circmscribe

More information

1 Lecture 20: Implicit differentiation

1 Lecture 20: Implicit differentiation Lecture 20: Implicit ifferentiation. Outline The technique of implicit ifferentiation Tangent lines to a circle Derivatives of inverse functions by implicit ifferentiation Examples.2 Implicit ifferentiation

More information

Lecture 3: Development of the Truss Equations.

Lecture 3: Development of the Truss Equations. 3.1 Derivation of the Stiffness Matrix for a Bar in Local Coorinates. In 3.1 we will perform Steps 1-4 of Logan s FEM. Derive the truss element equations. 1. Set the element type. 2. Select a isplacement

More information

DIFFERENTIAL GEOMETRY, LECTURE 15, JULY 10

DIFFERENTIAL GEOMETRY, LECTURE 15, JULY 10 DIFFERENTIAL GEOMETRY, LECTURE 15, JULY 10 5. Levi-Civita connection From now on we are intereste in connections on the tangent bunle T X of a Riemanninam manifol (X, g). Out main result will be a construction

More information

6 General properties of an autonomous system of two first order ODE

6 General properties of an autonomous system of two first order ODE 6 General properties of an autonomous system of two first orer ODE Here we embark on stuying the autonomous system of two first orer ifferential equations of the form ẋ 1 = f 1 (, x 2 ), ẋ 2 = f 2 (, x

More information

Calculus and optimization

Calculus and optimization Calculus an optimization These notes essentially correspon to mathematical appenix 2 in the text. 1 Functions of a single variable Now that we have e ne functions we turn our attention to calculus. A function

More information

Complementing the Lagrangian Density of the E. M. Field and the Surface Integral of the p-v Vector Product

Complementing the Lagrangian Density of the E. M. Field and the Surface Integral of the p-v Vector Product Applie Mathematics,,, 5-9 oi:.436/am..4 Pblishe Online Febrary (http://www.scirp.org/jornal/am) Complementing the Lagrangian Density of the E. M. Fiel an the Srface Integral of the p- Vector Proct Abstract

More information

The Principle of Least Action

The Principle of Least Action Chapter 7. The Principle of Least Action 7.1 Force Methos vs. Energy Methos We have so far stuie two istinct ways of analyzing physics problems: force methos, basically consisting of the application of

More information

Differentiation ( , 9.5)

Differentiation ( , 9.5) Chapter 2 Differentiation (8.1 8.3, 9.5) 2.1 Rate of Change (8.2.1 5) Recall that the equation of a straight line can be written as y = mx + c, where m is the slope or graient of the line, an c is the

More information

The continuity equation

The continuity equation Chapter 6 The continuity equation 61 The equation of continuity It is evient that in a certain region of space the matter entering it must be equal to the matter leaving it Let us consier an infinitesimal

More information

and from it produce the action integral whose variation we set to zero:

and from it produce the action integral whose variation we set to zero: Lagrange Multipliers Monay, 6 September 01 Sometimes it is convenient to use reunant coorinates, an to effect the variation of the action consistent with the constraints via the metho of Lagrange unetermine

More information

Designing Information Devices and Systems II Fall 2017 Note Theorem: Existence and Uniqueness of Solutions to Differential Equations

Designing Information Devices and Systems II Fall 2017 Note Theorem: Existence and Uniqueness of Solutions to Differential Equations EECS 6B Designing Information Devices an Systems II Fall 07 Note 3 Secon Orer Differential Equations Secon orer ifferential equations appear everywhere in the real worl. In this note, we will walk through

More information

Implicit Differentiation

Implicit Differentiation Implicit Differentiation Implicit Differentiation Using the Chain Rule In the previous section we focuse on the erivatives of composites an saw that THEOREM 20 (Chain Rule) Suppose that u = g(x) is ifferentiable

More information

The Ehrenfest Theorems

The Ehrenfest Theorems The Ehrenfest Theorems Robert Gilmore Classical Preliminaries A classical system with n egrees of freeom is escribe by n secon orer orinary ifferential equations on the configuration space (n inepenent

More information

Linear and quadratic approximation

Linear and quadratic approximation Linear an quaratic approximation November 11, 2013 Definition: Suppose f is a function that is ifferentiable on an interval I containing the point a. The linear approximation to f at a is the linear function

More information

u t v t v t c a u t b a v t u t v t b a

u t v t v t c a u t b a v t u t v t b a Nonlinear Dynamical Systems In orer to iscuss nonlinear ynamical systems, we must first consier linear ynamical systems. Linear ynamical systems are just systems of linear equations like we have been stuying

More information

A note on the Mooney-Rivlin material model

A note on the Mooney-Rivlin material model A note on the Mooney-Rivlin material moel I-Shih Liu Instituto e Matemática Universiae Feeral o Rio e Janeiro 2945-97, Rio e Janeiro, Brasil Abstract In finite elasticity, the Mooney-Rivlin material moel

More information

Section The Chain Rule and Implicit Differentiation with Application on Derivative of Logarithm Functions

Section The Chain Rule and Implicit Differentiation with Application on Derivative of Logarithm Functions Section 3.4-3.6 The Chain Rule an Implicit Differentiation with Application on Derivative of Logarithm Functions Ruipeng Shen September 3r, 5th Ruipeng Shen MATH 1ZA3 September 3r, 5th 1 / 3 The Chain

More information

Lecture 4 : General Logarithms and Exponentials. a x = e x ln a, a > 0.

Lecture 4 : General Logarithms and Exponentials. a x = e x ln a, a > 0. For a > 0 an x any real number, we efine Lecture 4 : General Logarithms an Exponentials. a x = e x ln a, a > 0. The function a x is calle the exponential function with base a. Note that ln(a x ) = x ln

More information

Diagonalization of Matrices Dr. E. Jacobs

Diagonalization of Matrices Dr. E. Jacobs Diagonalization of Matrices Dr. E. Jacobs One of the very interesting lessons in this course is how certain algebraic techniques can be use to solve ifferential equations. The purpose of these notes is

More information

Mathematical Model for Vibrations Analysis of The Tram Wheelset

Mathematical Model for Vibrations Analysis of The Tram Wheelset TH ARCHIVS OF TRANSPORT VOL. XXV-XXVI NO - 03 Mathematical Moel for Vibrations Analysis of The Tram Wheelset Anrzej Grzyb Karol Bryk ** Marek Dzik *** Receive March 03 Abstract This paper presents a mathematical

More information

A simple model for the small-strain behaviour of soils

A simple model for the small-strain behaviour of soils A simple moel for the small-strain behaviour of soils José Jorge Naer Department of Structural an Geotechnical ngineering, Polytechnic School, University of São Paulo 05508-900, São Paulo, Brazil, e-mail:

More information

Exam 2 Review Solutions

Exam 2 Review Solutions Exam Review Solutions 1. True or False, an explain: (a) There exists a function f with continuous secon partial erivatives such that f x (x, y) = x + y f y = x y False. If the function has continuous secon

More information

3.7 Implicit Differentiation -- A Brief Introduction -- Student Notes

3.7 Implicit Differentiation -- A Brief Introduction -- Student Notes Fin these erivatives of these functions: y.7 Implicit Differentiation -- A Brief Introuction -- Stuent Notes tan y sin tan = sin y e = e = Write the inverses of these functions: y tan y sin How woul we

More information

Problem set 2: Solutions Math 207B, Winter 2016

Problem 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 information

5-4 Electrostatic Boundary Value Problems

5-4 Electrostatic Boundary Value Problems 11/8/4 Section 54 Electrostatic Bounary Value Problems blank 1/ 5-4 Electrostatic Bounary Value Problems Reaing Assignment: pp. 149-157 Q: A: We must solve ifferential equations, an apply bounary conitions

More information

Thermal conductivity of graded composites: Numerical simulations and an effective medium approximation

Thermal conductivity of graded composites: Numerical simulations and an effective medium approximation JOURNAL OF MATERIALS SCIENCE 34 (999)5497 5503 Thermal conuctivity of grae composites: Numerical simulations an an effective meium approximation P. M. HUI Department of Physics, The Chinese University

More information

Shape functions in 1D

Shape functions in 1D MAE 44 & CIV 44 Introuction to Finite Elements Reaing assignment: ecture notes, ogan.,. Summary: Prof. Suvranu De Shape functions in D inear shape functions in D Quaratic an higher orer shape functions

More information

x f(x) x f(x) approaching 1 approaching 0.5 approaching 1 approaching 0.

x f(x) x f(x) approaching 1 approaching 0.5 approaching 1 approaching 0. Engineering Mathematics 2 26 February 2014 Limits of functions Consier the function 1 f() = 1. The omain of this function is R + \ {1}. The function is not efine at 1. What happens when is close to 1?

More information

Vectors in two dimensions

Vectors in two dimensions Vectors in two imensions Until now, we have been working in one imension only The main reason for this is to become familiar with the main physical ieas like Newton s secon law, without the aitional complication

More information

The Sokhotski-Plemelj Formula

The Sokhotski-Plemelj Formula hysics 24 Winter 207 The Sokhotski-lemelj Formula. The Sokhotski-lemelj formula The Sokhotski-lemelj formula is a relation between the following generalize functions (also calle istributions), ±iǫ = iπ(),

More information

Relation between the propagator matrix of geodesic deviation and the second-order derivatives of the characteristic function

Relation between the propagator matrix of geodesic deviation and the second-order derivatives of the characteristic function Journal of Electromagnetic Waves an Applications 203 Vol. 27 No. 3 589 60 http://x.oi.org/0.080/0920507.203.808595 Relation between the propagator matrix of geoesic eviation an the secon-orer erivatives

More information

23 Implicit differentiation

23 Implicit differentiation 23 Implicit ifferentiation 23.1 Statement The equation y = x 2 + 3x + 1 expresses a relationship between the quantities x an y. If a value of x is given, then a corresponing value of y is etermine. For

More information

CHAPTER 1 : DIFFERENTIABLE MANIFOLDS. 1.1 The definition of a differentiable manifold

CHAPTER 1 : DIFFERENTIABLE MANIFOLDS. 1.1 The definition of a differentiable manifold CHAPTER 1 : DIFFERENTIABLE MANIFOLDS 1.1 The efinition of a ifferentiable manifol Let M be a topological space. This means that we have a family Ω of open sets efine on M. These satisfy (1), M Ω (2) the

More information

model considered before, but the prey obey logistic growth in the absence of predators. In

model considered before, but the prey obey logistic growth in the absence of predators. In 5.2. First Orer Systems of Differential Equations. Phase Portraits an Linearity. Section Objective(s): Moifie Preator-Prey Moel. Graphical Representations of Solutions. Phase Portraits. Vector Fiels an

More information

Lecture 3Section 7.3 The Logarithm Function, Part II

Lecture 3Section 7.3 The Logarithm Function, Part II Lectre 3Section 7.3 The Logarithm Fnction, Part II Jiwen He Section 7.2: Highlights 2 Properties of the Log Fnction ln = t t, ln = 0, ln e =. (ln ) = > 0. ln(y) = ln + ln y, ln(/y) = ln ln y. ln ( r) =

More information

Sturm-Liouville Theory

Sturm-Liouville Theory LECTURE 5 Sturm-Liouville Theory In the three preceing lectures I emonstrate the utility of Fourier series in solving PDE/BVPs. As we ll now see, Fourier series are just the tip of the iceberg of the theory

More information

x f(x) x f(x) approaching 1 approaching 0.5 approaching 1 approaching 0.

x f(x) x f(x) approaching 1 approaching 0.5 approaching 1 approaching 0. Engineering Mathematics 2 26 February 2014 Limits of functions Consier the function 1 f() = 1. The omain of this function is R + \ {1}. The function is not efine at 1. What happens when is close to 1?

More information

QF101: Quantitative Finance September 5, Week 3: Derivatives. Facilitator: Christopher Ting AY 2017/2018. f ( x + ) f(x) f(x) = lim

QF101: Quantitative Finance September 5, Week 3: Derivatives. Facilitator: Christopher Ting AY 2017/2018. f ( x + ) f(x) f(x) = lim QF101: Quantitative Finance September 5, 2017 Week 3: Derivatives Facilitator: Christopher Ting AY 2017/2018 I recoil with ismay an horror at this lamentable plague of functions which o not have erivatives.

More information

Mark J. Machina CARDINAL PROPERTIES OF "LOCAL UTILITY FUNCTIONS"

Mark J. Machina CARDINAL PROPERTIES OF LOCAL UTILITY FUNCTIONS Mark J. Machina CARDINAL PROPERTIES OF "LOCAL UTILITY FUNCTIONS" This paper outlines the carinal properties of "local utility functions" of the type use by Allen [1985], Chew [1983], Chew an MacCrimmon

More information

Mechanics Physics 151

Mechanics Physics 151 Mechanics Physics 151 Lecture 3 Continuous Systems an Fiels (Chapter 13) Where Are We Now? We ve finishe all the essentials Final will cover Lectures 1 through Last two lectures: Classical Fiel Theory

More information

Math Chapter 2 Essentials of Calculus by James Stewart Prepared by Jason Gaddis

Math Chapter 2 Essentials of Calculus by James Stewart Prepared by Jason Gaddis Math 231 - Chapter 2 Essentials of Calculus by James Stewart Prepare by Jason Gais Chapter 2 - Derivatives 21 - Derivatives an Rates of Change Definition A tangent to a curve is a line that intersects

More information

Light flavor asymmetry of polarized quark distributions in thermodynamical bag model

Light flavor asymmetry of polarized quark distributions in thermodynamical bag model Inian Jornal of Pre & Applie Physics Vol. 5, April 014, pp. 19-3 Light flavor asymmetry of polarize qark istribtions in thermoynamical bag moel K Ganesamrthy a & S Mrganantham b* a Department of Physics,

More information

TMA 4195 Matematisk modellering Exam Tuesday December 16, :00 13:00 Problems and solution with additional comments

TMA 4195 Matematisk modellering Exam Tuesday December 16, :00 13:00 Problems and solution with additional comments Problem F U L W D g m 3 2 s 2 0 0 0 0 2 kg 0 0 0 0 0 0 Table : Dimension matrix TMA 495 Matematisk moellering Exam Tuesay December 6, 2008 09:00 3:00 Problems an solution with aitional comments The necessary

More information

Linear First-Order Equations

Linear First-Order Equations 5 Linear First-Orer Equations Linear first-orer ifferential equations make up another important class of ifferential equations that commonly arise in applications an are relatively easy to solve (in theory)

More information

Rank, Trace, Determinant, Transpose an Inverse of a Matrix Let A be an n n square matrix: A = a11 a1 a1n a1 a an a n1 a n a nn nn where is the jth col

Rank, Trace, Determinant, Transpose an Inverse of a Matrix Let A be an n n square matrix: A = a11 a1 a1n a1 a an a n1 a n a nn nn where is the jth col Review of Linear Algebra { E18 Hanout Vectors an Their Inner Proucts Let X an Y be two vectors: an Their inner prouct is ene as X =[x1; ;x n ] T Y =[y1; ;y n ] T (X; Y ) = X T Y = x k y k k=1 where T an

More information

Analytic Scaling Formulas for Crossed Laser Acceleration in Vacuum

Analytic Scaling Formulas for Crossed Laser Acceleration in Vacuum October 6, 4 ARDB Note Analytic Scaling Formulas for Crosse Laser Acceleration in Vacuum Robert J. Noble Stanfor Linear Accelerator Center, Stanfor University 575 San Hill Roa, Menlo Park, California 945

More information

SYDE 112, LECTURE 1: Review & Antidifferentiation

SYDE 112, LECTURE 1: Review & Antidifferentiation SYDE 112, LECTURE 1: Review & Antiifferentiation 1 Course Information For a etaile breakown of the course content an available resources, see the Course Outline. Other relevant information for this section

More information

GLOBAL PHASE PORTRAITS OF SOME REVERSIBLE CUBIC CENTERS WITH NONCOLLINEAR SINGULARITIES

GLOBAL PHASE PORTRAITS OF SOME REVERSIBLE CUBIC CENTERS WITH NONCOLLINEAR SINGULARITIES This is a preprint of: Global phase portraits of some reversible cbic centers with noncollinear singlarities, Magalena Cabergh, Joan Torregrosa, Internat. J. Bifr. Chaos Appl. Sci. Engrg., vol. 23(9),

More information

Math 1271 Solutions for Fall 2005 Final Exam

Math 1271 Solutions for Fall 2005 Final Exam Math 7 Solutions for Fall 5 Final Eam ) Since the equation + y = e y cannot be rearrange algebraically in orer to write y as an eplicit function of, we must instea ifferentiate this relation implicitly

More information

Chapter 9 Method of Weighted Residuals

Chapter 9 Method of Weighted Residuals Chapter 9 Metho of Weighte Resiuals 9- Introuction Metho of Weighte Resiuals (MWR) is an approimate technique for solving bounary value problems. It utilizes a trial functions satisfying the prescribe

More information

ON THE RIEMANN EXTENSION OF THE SCHWARZSCHILD METRICS

ON THE RIEMANN EXTENSION OF THE SCHWARZSCHILD METRICS ON THE RIEANN EXTENSION OF THE SCHWARZSCHILD ETRICS Valerii Dryuma arxiv:gr-qc/040415v1 30 Apr 004 Institute of athematics an Informatics, AS R, 5 Acaemiei Street, 08 Chisinau, olova, e-mail: valery@ryuma.com;

More information

05 The Continuum Limit and the Wave Equation

05 The Continuum Limit and the Wave Equation Utah State University DigitalCommons@USU Founations of Wave Phenomena Physics, Department of 1-1-2004 05 The Continuum Limit an the Wave Equation Charles G. Torre Department of Physics, Utah State University,

More information