ENGI Duffing s Equation Page 4.65
|
|
- Jason Hood
- 5 years ago
- Views:
Transcription
1 ENGI Duffing s Equation Page Duffing s Equation Among the simplest models of damped non-linear forced oscillations of a mechanical or electrical system with a cubic stiffness term is Duffing s equation: d x dx 3 + a + bx + cx = d cosωt () In section 4.0, we considered the simple undamped pendulum: d x g + sin x = 0 () L When x is very small, sin x x and () reduces to the ODE for simple harmonic motion. 3 x The next order approximation is sin x x, so that () becomes 6 3 d x g g x + x = 0 (3) L L 6 dx If we add a damping term a and a forcing function d cos ωt, then (3) becomes Duffing s equation ().
2 ENGI Duffing s Equation Page 4.66 Special Case : Conduct a stability analysis for the undamped unforced Duffing s equation d x 3 + ω x + cx = 0 (4) The equivalent first order system is dx dy = y = ω x cx 3 (5) Critical points: ( y = 0) and x= 0 or x = c Near (0, 0) the linear approximation is x 0 x = y ω 0 y The characteristic equation is det(a - λi) = 0 λ + ω = 0 The eigenvalues are a pure imaginary pair (0, 0) is a centre. It is stable but not asymptotically stable. If c > 0, then this is the only critical point of (4). (6) If c < 0, then there are two other critical points, at ±,0. c
3 ENGI Duffing s Equation Page 4.67 Special Case : (continued) Near ±,0 c dx P P = ( y) x + y near ( ± ω / c, 0) near ( ± ω / c, 0) x ( / c, 0) c y ± ω ( ± ω/ c, 0) = ( 0) x + () y y c = dy 3 Q Q = ( ω x cx ) x + y near ( ) ( ± ω/ c, 0) near ± ω / c, 0 x ( / c, 0) c y ± ω ( ± ω/ c, 0) = ω 3c x + ( 0) y = ω x c c c The linear approximation to (5) near ±,0 is c x 0 x = c y ω 0 y The characteristic equation is det(a - λi) = 0 λ ω = 0 The eigenvalues are a distinct real pair with opposite sign ±,0 are saddle points. They are unstable. c (7)
4 ENGI Duffing s Equation Page 4.68 Exact Solution of Special Case : The system (5), dx = y dy 3 = ω x cx 3 dy dy dx ω x c x = = dx y 4 3 y A ω x cx ydy = ( ω x cx ) dx = 4 4 cx Therefore the orbits in the phase space (x, y) are y = A ω x, where A is an arbitrary constant. If c > 0, then all orbits are closed, about the centre at (0, 0).
5 ENGI Duffing s Equation Page 4.69 Special Case : (continued) If c < 0, then those orbits far enough away from the centre are open, due to the influence of the saddle points at ±,0. c The part of the phase space between the two saddle points resembles that for the undamped pendulum on page 4.07: The orbits passing through the saddle points separates closed orbits from open orbits and is called the separatrix. The positive y axis intercept of each orbit is just the value of A for that orbit. The separatrix has x axis intercepts at the saddle points. Therefore, for the separatrix, cx c ω A = ω x + = ω + c x / c = c c =± ω The equation of the separatrix is 4 4 ω cx y = ω x +, ( c < 0) c
6 ENGI Duffing s Equation Page 4.70 Special Case : Conduct a stability analysis for the damped unforced Duffing s equation d x dx 3 + a + ω x + cx = 0 (8) The equivalent first order system is dx dy = y 3 = ω x cx ay (9) The critical points are the same as in special case : ( y = 0) and x= 0 or x = c Near (0, 0) the linear approximation is x 0 x = y ω a y The characteristic equation is det(a - λi) = 0 λ + aλ + ω = 0 a ± a 4ω λ = The critical point is stable if a > 0 and unstable if a < 0. It is a focus if a 4ω < 0 and a node otherwise. (0) If c > 0, then this is the only critical point of (8). If c < 0, then there are two other critical points, at ±,0. c
7 ENGI Duffing s Equation Page 4.7 Special Case : (continued) Near ±,0 c dx P P = ( y) x y near ( ± ω/ c, 0) near ( / c, 0) x + ± ω ( / c, 0) c y ± ω ( ± ω/ c, 0) = ( 0) x + () y = y c dy 3 Q Q = ( ω x cx ay) x y near ( ) ( ± ω/ c, 0) near / c, 0 x + ± ω ( / c, 0) c y ± ω ( ± ω/ c, 0) = ω 3c x + ( a) y = ω x ay c c c The linear approximation to (9) near ±,0 is c x 0 x c = y ω a y The characteristic equation is det(a - λi) = 0 λ + aλ ω = 0 a ± a + 8ω λ = The eigenvalues are a distinct real pair with opposite sign ±,0 are saddle points. They are unstable. c () The presence of the damping term changes the centre into a stable focus (for physically reasonable values of a, ω and c, or, for particularly strong damping, a stable node). The form of the separatrix is more complicated, as trajectories leaving either saddle point in the direction of the origin are swept by the damping term into the focus (or node) instead of moving around the centre to the other saddle point. There are no closed orbits; just orbits that terminate at the origin or a saddle point and orbits that retreat to infinity.
8 ENGI Duffing s Equation Page 4.7 Special Case : (continued) An enhanced sample phase portrait plot from Maple is shown here: and, zooming in,
9 ENGI More Examples Page More Examples Example 4.. Examine the stability of the linear second order differential equation d x dx + + ( 4π + ) x = 0 and find the complete solution for the initial conditions x 0 = 0, y 0 = x 0 = π. ( ) ( ) ( ) The system can be rewritten as the first order system x 0 x = y ( 4π + ) y The only critical point is at the origin. ( ) ( ) ( ) ( ) ( ) D = a d + 4bc = π + = 6π D < 0 and (a + d) < 0 the critical point is an asymptotically stable focus. ( a+ d) ± D ± 6π λ = = = ± π j Using the formula on page 4.30, u =, v π, u d, c ( 4π ) The general solution is t x t = e c cos π t πsin πt + c π cos πt + sin πt = = + = = +. () ( 3( ) 4( )) t () = ( 4π + )( cosπ + sinπ ) y t e c t c t 3 4 [and one can check that dx = y is indeed true.] π ( x( 0, ) y( 0) ) = ( 0,π) ( c3+ πc4, c3) = 0, 4 π + π ( c3, c4) =, 4π + 4π +
10 ENGI More Examples Page 4.74 Example 4.. (continued) The complete solution is t e ( xt (), yt ()) = 4π + ( π ( πsinπt cosπt) + ( π cosπt + sin πt), ( 4π + )( π cosπt sinπt) ) t ( x( t), y ( t) ) = e ( sin πt, ( π cos πt sin πt) ) As t, both functions x(t) and y(t) tend to zero. The resulting phase space diagram is
11 ENGI More Examples Page 4.75 Example 4.. Examine the stability of the linear second order differential equation d x dx + + ( 4π + ) x = ( 3π + ) cosπt πsinπ t The complete solution for the initial conditions x( 0) =, y( 0) = x ( 0) = π. can be obtained by building upon the solution to Example 4.. and is ( ) t t ( x() t, y () t ) = cosπt + e sin πt, πsinπt + e ( π cos πt sin πt) In this case, the steady state solution (after the transient terms have vanished) is lim x t, y t = cos πt, π sinπt t ( ( ) ( )) ( ) so that the orbit in the phase space approaches the ellipse π x + y = π This ellipse must therefore be the limit cycle for the system. A plot of the orbit in the phase space is shown here. Note how the solution curve in the phase space can wander inside and outside the limit cycle more than once, before finally settling down to its asymptotic approach as the transient terms become negligible.
12 ENGI More Examples Page 4.76 Example 4.. (continued) Different sets of initial conditions can generate orbits that look very different at first, before they settle down into their steady-state configuration near the limit cycle.
13 ENGI Liénard s Theorem Page Liénard s Theorem If f (x) is an even function for all x and g(x) is an odd function for all x and g(x) > 0 for all x > 0 and F( x) ( ) then the system x = f t is such that F (x) = 0 has exactly one positive root, γ, and 0 F (x) < 0 for 0 < x < γ and F (x) > 0 and non-decreasing for x > γ, or, equivalently, ( ) ( ) x = y, y = f x y g x d x dx + f ( x) + g( x) = 0 has a unique limit cycle enclosing the origin and that limit cycle is asymptotically stable. When all of the conditions of Liénard s theorem are satisfied, the system has exactly one periodic solution, towards which all other trajectories spiral as t. Example 4.3. Let f (x) = μ ( x ) and g(x) = x, (with μ > 0), then Liénard s ODE becomes d x dx μ ( x ) + x = 0 which is Van der Pol s equation (section 4.08). Checking the conditions of Liénard s theorem: f (x) = μ ( x ) is an even function. g(x) = x is an odd function, positive for all x > 0. x x 3 3 t x x F( x) = μ( t ) = μ t = + μ x = μx F (x) = 0 has only one positive root, γ = 3. F (x) < 0 for 0 < x < 3 and F (x) > 0 and increasing for x > 3. 0 Therefore Van der Pol s equation possesses a unique and asymptotically stable limit cycle.
14 ENGI Gradient, Divergence, Curl Page The Gradient Operator A brief review is provided here for the gradient operator in both Cartesian and orthogonal non-cartesian coordinate systems. Sections in this Chapter: 5.0 Gradient, Divergence, Curl and Laplacian (Cartesian) 5.0 Differentiation in Orthogonal Curvilinear Coordinate Systems 5.03 Summary Table for the Gradient Operator 5.04 Derivatives of Basis Vectors 5.0 Gradient, Divergence, Curl and Laplacian (Cartesian) Let z be a function of two independent variables (x, y), so that z = f (x, y). The function z = f (x, y) defines a surface in 3. At any point (x, y) in the x-y plane, the direction in which one must travel in order to experience the greatest possible rate of increase in z at that point is the direction of the gradient vector, f ˆ f f i ˆj x y The magnitude of the gradient vector is that greatest possible rate of increase in z at that point. The gradient vector is not constant everywhere, unless the surface is a plane. (The symbol is usually pronounced del ). The concept of the gradient vector can be extended to functions of any number of variables. If u = f (x, y, z, t), then f f f f f x y z t. If v is a function of position r and time t, while position is in turn a function of time, then by the chain rule of differentiation, d v v dx v dy v dz v dr v v x y z t t d v v v v t which is of use in the study of fluid dynamics. T
15 ENGI Gradient, Divergence, Curl Page 5.0 The gradient operator can also be applied to vectors via the scalar (dot) and vector (cross) products: The divergence of a vector field F(x, y, z) is div F T F F F F F F F3 x y z x y z T 3 A region free of sources and sinks will have zero divergence: the total flux into any region is balanced by the total flux out from that region. The curl of a vector field F(x, y, z) is ˆi curl F F ˆj F kˆ F3 F F x y z F F3 y z x F F F3 z x y In an irrotational field, curl Whenever or Proof: F 0. F for some twice differentiable potential function, curl F 0 curl grad T T F F F3 F x y z 0 ˆ i 0 x x y z z y curl ˆ j 0 y y z x x z ˆ k 0 z z x y y x
16 ENGI Gradient, Divergence, Curl Page 5.03 Among many identities involving the gradient operator is div curl F F 0 for all twice-differentiable vector functions F Proof: div curl F F3 F F F 3 F F x y z y z x z x y F3 F F F 3 F F 0 x y x z y z y x z x z y The divergence of the gradient of a scalar function is the Laplacian: div grad f f f f f f x y z for all twice-differentiable scalar functions f. In orthogonal non-cartesian coordinate systems, the expressions for the gradient operator are not as simple.
17 ENGI Curvilinear Gradient Page Differentiation in Orthogonal Curvilinear Coordinate Systems For any orthogonal curvilinear coordinate system (u, u, u 3 ) in the unit tangent vectors along the curvilinear axes are ˆ r e T, h u where the scale factors h i r u i. 3, ˆi i i i The displacement vector r can then be written as r ueˆ ˆ ˆ ue u3e 3, where the unit vectors ê i form an orthonormal basis for 3. 0 i j eˆi e ˆ j i j i j The differential displacement vector dr is (by the Chain Rule) r r r dr du du du h du eˆ h du eˆ h du e ˆ u u u3 and the differential arc length ds is given by ds dr dr h du h du h3 du3 The element of volume dv is x, y, z dv h h h du du du du du du 3 3 u 3, u, u3 Jacobian x y u u x y u u x y u u 3 3 z u z u z u Example 5.0.: Find the scale factor h θ for the spherical polar coordinate system x, y, z r sin cos, r sin sin, r cos : T r x y z r cos cos r cos sin r sin r h r r r cos cos cos sin sin r cos cos sin sin cos sin 3 du du du r r r 3 T
18 ENGI Gradient Summary Page Summary Table for the Gradient Operator eˆ eˆ ˆ e3 Gradient operator h u h u h u 3 3 eˆ V eˆ ˆ V e3 V Gradient V h u h u h u 3 3 Divergence h h F h h F h h F F h h h 3 u u u Curl F h h h h eˆ h F u h eˆ h F 3 u h eˆ h F u3 Laplacian V = h h h 3 h h u h 3 V u u h3 h V h u u 3 h h h3 V u 3 Scale factors: Cartesian: h x = h y = h z =. Cylindrical polar: h = h z =, h =. Spherical polar: h r =, h = r, h = r sin. Example 5.03.: The Laplacian of V in spherical polars is sin V V r sin V V r sin r r sin or V V V cot V V V r r r r r sin
19 ENGI Gradient Summary Page 5.06 Example A potential function V r is spherically symmetric, (that is, its value depends only on the distance r from the origin), due solely to a point source at the origin. There are no other sources or sinks anywhere in V r. 3. Deduce the functional form of V r is spherically symmetric V r,, f r In any regions not containing any sources of the vector field, the divergence of the vector field F V (and therefore the Laplacian of the associated potential function V) must be zero. Therefore, for all r 0, div F V V 0 But sin V r sin V r sin r r sin d V dv sin r r sin dr dr d dv dv dv r 0 r B B r dr dr dr dr V Br A, where A, B are arbitrary constants of integration. Therefore the potential function must be of the form B V r,, A r This is the standard form of the potential function associated with a force that obeys the inverse square law F. r
20 ENGI Basis Vectors Page Derivatives of Basis Vectors d Cartesian: ˆ d ˆ d i j kˆ 0 r x ˆi y ˆj z kˆ v x ˆi y ˆj z kˆ Cylindrical Polar Coordinates: x cos, y sin, z z d d ˆ ˆ r ˆ z kˆ d ˆ d ˆ v ˆ ˆ z kˆ d kˆ 0 [radial and transverse components of v ] Spherical Polar Coordinates. The declination angle θ is the angle between the positive z axis and the radius vector r. 0 < θ < π. The azimuth angle is the angle on the x-y plane, measured anticlockwise from the positive x axis, of the shadow of the radius vector. 0 < < π. z = r cos θ. The shadow of the radius vector on the x-y plane has length r sin θ. It then follows that x = r sin θ cos and y = r sin θ sin. d d ˆ ˆ d r sin ˆ r r r ˆ d ˆ d d rˆ cos ˆ v r rˆ r ˆ r sin ˆ d ˆ d sin rˆ cos ˆ
21 ENGI Basis Vectors Page 5.08 Example Find the velocity and acceleration in cylindrical polar coordinates for a particle travelling along the helix x = 3 cos t, y = 3 sin t, z = t. Cylindrical polar coordinates: x cos, y sin, z z y x y, tan x 9cos t 9sin t sin t tan tan t t 3cost r 3 ˆ z kˆ z t z dr v ˆ ˆ zkˆ 0 ˆ 3 ˆ kˆ 6 ˆ kˆ [The velocity has no radial component the helix remains the same distance from the z axis at all times.] d v a 6 ˆ kˆ 6 ˆ 0 ˆ [The acceleration vector points directly at the z axis at all times.] Other examples are in the problem sets. END OF CHAPTER 5
ENGI Gradient, Divergence, Curl Page 5.01
ENGI 94 5. - Gradient, Divergence, Curl Page 5. 5. The Gradient Operator A brief review is provided here for the gradient operator in both Cartesian and orthogonal non-cartesian coordinate systems. Sections
More informationENGI Gradient, Divergence, Curl Page 5.01
ENGI 940 5.0 - Gradient, Divergence, Curl Page 5.0 5. e Gradient Operator A brief review is provided ere for te gradient operator in bot Cartesian and ortogonal non-cartesian coordinate systems. Sections
More informationENGI 9420 Lecture Notes 4 - Stability Analysis Page Stability Analysis for Non-linear Ordinary Differential Equations
ENGI 940 Lecture Notes 4 - Stability Analysis Page 4.01 4. Stability Analysis for Non-linear Ordinary Differential Equations A pair of simultaneous first order homogeneous linear ordinary differential
More information154 Chapter 9 Hints, Answers, and Solutions The particular trajectories are highlighted in the phase portraits below.
54 Chapter 9 Hints, Answers, and Solutions 9. The Phase Plane 9.. 4. The particular trajectories are highlighted in the phase portraits below... 3. 4. 9..5. Shown below is one possibility with x(t) and
More informationENGI 4430 Line Integrals; Green s Theorem Page 8.01
ENGI 443 Line Integrals; Green s Theorem Page 8. 8. Line Integrals Two applications of line integrals are treated here: the evaluation of work done on a particle as it travels along a curve in the presence
More informationSTABILITY. Phase portraits and local stability
MAS271 Methods for differential equations Dr. R. Jain STABILITY Phase portraits and local stability We are interested in system of ordinary differential equations of the form ẋ = f(x, y), ẏ = g(x, y),
More informationMath 266: Phase Plane Portrait
Math 266: Phase Plane Portrait Long Jin Purdue, Spring 2018 Review: Phase line for an autonomous equation For a single autonomous equation y = f (y) we used a phase line to illustrate the equilibrium solutions
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 informationENGI Linear Approximation (2) Page Linear Approximation to a System of Non-Linear ODEs (2)
ENGI 940 4.06 - Linear Approximation () Page 4. 4.06 Linear Approximation to a System of Non-Linear ODEs () From sections 4.0 and 4.0, the non-linear system dx dy = x = P( x, y), = y = Q( x, y) () with
More informationr t t x t y t z t, y t are zero, then construct a table for all four functions. dy dx 0 and 0 dt dt horizontal tangent vertical tangent
3. uggestions for the Formula heets Below are some suggestions for many more formulae than can be placed easily on both sides of the two standard 8½"" sheets of paper for the final examination. It is strongly
More informationA plane autonomous system is a pair of simultaneous first-order differential equations,
Chapter 11 Phase-Plane Techniques 11.1 Plane Autonomous Systems A plane autonomous system is a pair of simultaneous first-order differential equations, ẋ = f(x, y), ẏ = g(x, y). This system has an equilibrium
More informationENGI 9420 Lecture Notes 4 - Stability Analysis Page Stability Analysis for Non-linear Ordinary Differential Equations
ENGI 940 Lecture Notes 4 - Stability Analysis Page 4.0 4. Stability Analysis for Non-linear Ordinary Differential Equations A pair of simultaneous first order homogeneous linear ordinary differential equations
More informationCurvilinear coordinates
C Curvilinear coordinates The distance between two points Euclidean space takes the simplest form (2-4) in Cartesian coordinates. The geometry of concrete physical problems may make non-cartesian coordinates
More informationNonlinear Autonomous Systems of Differential
Chapter 4 Nonlinear Autonomous Systems of Differential Equations 4.0 The Phase Plane: Linear Systems 4.0.1 Introduction Consider a system of the form x = A(x), (4.0.1) where A is independent of t. Such
More informationENGI 4430 Parametric Vector Functions Page dt dt dt
ENGI 4430 Parametric Vector Functions Page 2-01 2. Parametric Vector Functions (continued) Any non-zero vector r can be decomposed into its magnitude r and its direction: r rrˆ, where r r 0 Tangent Vector:
More information1. < 0: the eigenvalues are real and have opposite signs; the fixed point is a saddle point
Solving a Linear System τ = trace(a) = a + d = λ 1 + λ 2 λ 1,2 = τ± = det(a) = ad bc = λ 1 λ 2 Classification of Fixed Points τ 2 4 1. < 0: the eigenvalues are real and have opposite signs; the fixed point
More informationIntroduction to Vector Calculus (29) SOLVED EXAMPLES. (d) B. C A. (f) a unit vector perpendicular to both B. = ˆ 2k = = 8 = = 8
Introduction to Vector Calculus (9) SOLVED EXAMPLES Q. If vector A i ˆ ˆj k, ˆ B i ˆ ˆj, C i ˆ 3j ˆ kˆ (a) A B (e) A B C (g) Solution: (b) A B (c) A. B C (d) B. C A then find (f) a unit vector perpendicular
More informationENGI 4430 Line Integrals; Green s Theorem Page 8.01
ENGI 4430 Line Integrals; Green s Theorem Page 8.01 8. Line Integrals Two applications of line integrals are treated here: the evaluation of work done on a particle as it travels along a curve in the presence
More informationSample Solutions of Assignment 10 for MAT3270B
Sample Solutions of Assignment 1 for MAT327B 1. For the following ODEs, (a) determine all critical points; (b) find the corresponding linear system near each critical point; (c) find the eigenvalues of
More informationChapter 6 Nonlinear Systems and Phenomena. Friday, November 2, 12
Chapter 6 Nonlinear Systems and Phenomena 6.1 Stability and the Phase Plane We now move to nonlinear systems Begin with the first-order system for x(t) d dt x = f(x,t), x(0) = x 0 In particular, consider
More informationAE301 Aerodynamics I UNIT B: Theory of Aerodynamics
AE301 Aerodynamics I UNIT B: Theory of Aerodynamics ROAD MAP... B-1: Mathematics for Aerodynamics B-: Flow Field Representations B-3: Potential Flow Analysis B-4: Applications of Potential Flow Analysis
More informationLECTURE 8: DYNAMICAL SYSTEMS 7
15-382 COLLECTIVE INTELLIGENCE S18 LECTURE 8: DYNAMICAL SYSTEMS 7 INSTRUCTOR: GIANNI A. DI CARO GEOMETRIES IN THE PHASE SPACE Damped pendulum One cp in the region between two separatrix Separatrix Basin
More informationNotes 3 Review of Vector Calculus
ECE 3317 Applied Electromagnetic Waves Prof. David R. Jackson Fall 2018 A ˆ Notes 3 Review of Vector Calculus y ya ˆ y x xa V = x y ˆ x Adapted from notes by Prof. Stuart A. Long 1 Overview Here we present
More informationMathematical Tripos Part IA Lent Term Example Sheet 1. Calculate its tangent vector dr/du at each point and hence find its total length.
Mathematical Tripos Part IA Lent Term 205 ector Calculus Prof B C Allanach Example Sheet Sketch the curve in the plane given parametrically by r(u) = ( x(u), y(u) ) = ( a cos 3 u, a sin 3 u ) with 0 u
More informationENGI 4430 Surface Integrals Page and 0 2 r
ENGI 4430 Surface Integrals Page 9.01 9. Surface Integrals - Projection Method Surfaces in 3 In 3 a surface can be represented by a vector parametric equation r x u, v ˆi y u, v ˆj z u, v k ˆ where u,
More informationLecture 38. Almost Linear Systems
Math 245 - Mathematics of Physics and Engineering I Lecture 38. Almost Linear Systems April 20, 2012 Konstantin Zuev (USC) Math 245, Lecture 38 April 20, 2012 1 / 11 Agenda Stability Properties of Linear
More informationENGI Parametric Vector Functions Page 5-01
ENGI 3425 5. Parametric Vector Functions Page 5-01 5. Parametric Vector Functions Contents: 5.1 Arc Length (Cartesian parametric and plane polar) 5.2 Surfaces of Revolution 5.3 Area under a Parametric
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 information2.20 Fall 2018 Math Review
2.20 Fall 2018 Math Review September 10, 2018 These notes are to help you through the math used in this class. This is just a refresher, so if you never learned one of these topics you should look more
More informationGeometry and Motion Selected answers to Sections A and C Dwight Barkley 2016
MA34 Geometry and Motion Selected answers to Sections A and C Dwight Barkley 26 Example Sheet d n+ = d n cot θ n r θ n r = Θθ n i. 2. 3. 4. Possible answers include: and with opposite orientation: 5..
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 informationCourse Notes Math 275 Boise State University. Shari Ultman
Course Notes Math 275 Boise State University Shari Ultman Fall 2017 Contents 1 Vectors 1 1.1 Introduction to 3-Space & Vectors.............. 3 1.2 Working With Vectors.................... 7 1.3 Introduction
More informationENGI Partial Differentiation Page y f x
ENGI 3424 4 Partial Differentiation Page 4-01 4. Partial Differentiation For functions of one variable, be found unambiguously by differentiation: y f x, the rate of change of the dependent variable can
More informationSummary of various integrals
ummary of various integrals Here s an arbitrary compilation of information about integrals Moisés made on a cold ecember night. 1 General things o not mix scalars and vectors! In particular ome integrals
More informationCBE 6333, R. Levicky 1. Orthogonal Curvilinear Coordinates
CBE 6333, R. Levicky 1 Orthogonal Curvilinear Coordinates Introduction. Rectangular Cartesian coordinates are convenient when solving problems in which the geometry of a problem is well described by the
More informationDisclaimer: This Final Exam Study Guide is meant to help you start studying. It is not necessarily a complete list of everything you need to know.
Disclaimer: This is meant to help you start studying. It is not necessarily a complete list of everything you need to know. The MTH 234 final exam mainly consists of standard response questions where students
More informationEELE 3331 Electromagnetic I Chapter 3. Vector Calculus. Islamic University of Gaza Electrical Engineering Department Dr.
EELE 3331 Electromagnetic I Chapter 3 Vector Calculus Islamic University of Gaza Electrical Engineering Department Dr. Talal Skaik 2012 1 Differential Length, Area, and Volume This chapter deals with integration
More informationexample consider flow of water in a pipe. At each point in the pipe, the water molecule has a velocity
Module 1: A Crash Course in Vectors Lecture 1: Scalar and Vector Fields Objectives In this lecture you will learn the following Learn about the concept of field Know the difference between a scalar field
More informationMATH H53 : Final exam
MATH H53 : Final exam 11 May, 18 Name: You have 18 minutes to answer the questions. Use of calculators or any electronic items is not permitted. Answer the questions in the space provided. If you run out
More informationChapter 3 Vectors. 3.1 Vector Analysis
Chapter 3 Vectors 3.1 Vector nalysis... 1 3.1.1 Introduction to Vectors... 1 3.1.2 Properties of Vectors... 1 3.2 Coordinate Systems... 6 3.2.1 Cartesian Coordinate System... 6 3.2.2 Cylindrical Coordinate
More informationProblem Set Number 5, j/2.036j MIT (Fall 2014)
Problem Set Number 5, 18.385j/2.036j MIT (Fall 2014) Rodolfo R. Rosales (MIT, Math. Dept.,Cambridge, MA 02139) Due Fri., October 24, 2014. October 17, 2014 1 Large µ limit for Liénard system #03 Statement:
More informationMAC2313 Final A. (5 pts) 1. How many of the following are necessarily true? i. The vector field F = 2x + 3y, 3x 5y is conservative.
MAC2313 Final A (5 pts) 1. How many of the following are necessarily true? i. The vector field F = 2x + 3y, 3x 5y is conservative. ii. The vector field F = 5(x 2 + y 2 ) 3/2 x, y is radial. iii. All constant
More informationDivergence Theorem December 2013
Divergence Theorem 17.3 11 December 2013 Fundamental Theorem, Four Ways. b F (x) dx = F (b) F (a) a [a, b] F (x) on boundary of If C path from P to Q, ( φ) ds = φ(q) φ(p) C φ on boundary of C Green s Theorem:
More informationChapter 1. Vector Algebra and Vector Space
1. Vector Algebra 1.1. Scalars and vectors Chapter 1. Vector Algebra and Vector Space The simplest kind of physical quantity is one that can be completely specified by its magnitude, a single number, together
More information5. Suggestions for the Formula Sheets
5. uggestions for the Formula heets Below are some suggestions for many more formulae than can be placed easily on both sides of the two standard 8½" " sheets of paper for the final examination. It is
More informationMath 312 Lecture Notes Linear Two-dimensional Systems of Differential Equations
Math 2 Lecture Notes Linear Two-dimensional Systems of Differential Equations Warren Weckesser Department of Mathematics Colgate University February 2005 In these notes, we consider the linear system of
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 informationEE2007: Engineering Mathematics II Vector Calculus
EE2007: Engineering Mathematics II Vector Calculus Ling KV School of EEE, NTU ekvling@ntu.edu.sg Rm: S2-B2b-22 Ver 1.1: Ling KV, October 22, 2006 Ver 1.0: Ling KV, Jul 2005 EE2007/Ling KV/Aug 2006 EE2007:
More informationDifferential Operators and the Divergence Theorem
1 of 6 1/15/2007 6:31 PM Differential Operators and the Divergence Theorem One of the most important and useful mathematical constructs is the "del operator", usually denoted by the symbol Ñ (which is
More informationDivergence Theorem Fundamental Theorem, Four Ways. 3D Fundamental Theorem. Divergence Theorem
Divergence Theorem 17.3 11 December 213 Fundamental Theorem, Four Ways. b F (x) dx = F (b) F (a) a [a, b] F (x) on boundary of If C path from P to Q, ( φ) ds = φ(q) φ(p) C φ on boundary of C Green s Theorem:
More informationSample Solutions of Assignment 9 for MAT3270B
Sample Solutions of Assignment 9 for MAT370B. For the following ODEs, find the eigenvalues and eigenvectors, and classify the critical point 0,0 type and determine whether it is stable, asymptotically
More informationTwo dimensional oscillator and central forces
Two dimensional oscillator and central forces September 4, 04 Hooke s law in two dimensions Consider a radial Hooke s law force in -dimensions, F = kr where the force is along the radial unit vector and
More informationProblem Solving 1: Line Integrals and Surface Integrals
A. Line Integrals MASSACHUSETTS INSTITUTE OF TECHNOLOY Department of Physics Problem Solving 1: Line Integrals and Surface Integrals The line integral of a scalar function f ( xyz),, along a path C is
More informationDef. (a, b) is a critical point of the autonomous system. 1 Proper node (stable or unstable) 2 Improper node (stable or unstable)
Types of critical points Def. (a, b) is a critical point of the autonomous system Math 216 Differential Equations Kenneth Harris kaharri@umich.edu Department of Mathematics University of Michigan November
More informationMultiple Integrals and Vector Calculus (Oxford Physics) Synopsis and Problem Sets; Hilary 2015
Multiple Integrals and Vector Calculus (Oxford Physics) Ramin Golestanian Synopsis and Problem Sets; Hilary 215 The outline of the material, which will be covered in 14 lectures, is as follows: 1. Introduction
More information1 The pendulum equation
Math 270 Honors ODE I Fall, 2008 Class notes # 5 A longer than usual homework assignment is at the end. The pendulum equation We now come to a particularly important example, the equation for an oscillating
More informationMathematical Preliminaries
Mathematical Preliminaries Introductory Course on Multiphysics Modelling TOMAZ G. ZIELIŃKI bluebox.ippt.pan.pl/ tzielins/ Table of Contents Vectors, tensors, and index notation. Generalization of the concept
More informationGRE Math Subject Test #5 Solutions.
GRE Math Subject Test #5 Solutions. 1. E (Calculus) Apply L Hôpital s Rule two times: cos(3x) 1 3 sin(3x) 9 cos(3x) lim x 0 = lim x 2 x 0 = lim 2x x 0 = 9. 2 2 2. C (Geometry) Note that a line segment
More informationReview of Vector Analysis in Cartesian Coordinates
Review of Vector Analysis in Cartesian Coordinates 1 Scalar: A quantity that has magnitude, but no direction. Examples are mass, temperature, pressure, time, distance, and real numbers. Scalars are usually
More informationVector Calculus, Maths II
Section A Vector Calculus, Maths II REVISION (VECTORS) 1. Position vector of a point P(x, y, z) is given as + y and its magnitude by 2. The scalar components of a vector are its direction ratios, and represent
More informationDifferential Equations 2280 Sample Midterm Exam 3 with Solutions Exam Date: 24 April 2015 at 12:50pm
Differential Equations 228 Sample Midterm Exam 3 with Solutions Exam Date: 24 April 25 at 2:5pm Instructions: This in-class exam is 5 minutes. No calculators, notes, tables or books. No answer check is
More informationENGI 9420 Lecture Notes 1 - ODEs Page 1.01
ENGI 940 Lecture Notes - ODEs Page.0. Ordinary Differential Equations An equation involving a function of one independent variable and the derivative(s) of that function is an ordinary differential equation
More informationCURRENT MATERIAL: Vector Calculus.
Math 275, section 002 (Ultman) Fall 2011 FINAL EXAM REVIEW The final exam will be held on Wednesday 14 December from 10:30am 12:30pm in our regular classroom. You will be allowed both sides of an 8.5 11
More informationMath Review Night: Work and the Dot Product
Math Review Night: Work and the Dot Product Dot Product A scalar quantity Magnitude: A B = A B cosθ The dot product can be positive, zero, or negative Two types of projections: the dot product is the parallel
More informationLecture Notes Introduction to Vector Analysis MATH 332
Lecture Notes Introduction to Vector Analysis MATH 332 Instructor: Ivan Avramidi Textbook: H. F. Davis and A. D. Snider, (WCB Publishers, 1995) New Mexico Institute of Mining and Technology Socorro, NM
More information1 Curvilinear Coordinates
MATHEMATICA PHYSICS PHYS-2106/3 Course Summary Gabor Kunstatter, University of Winnipeg April 2014 1 Curvilinear Coordinates 1. General curvilinear coordinates 3-D: given or conversely u i = u i (x, y,
More informationDirectional Derivative and the Gradient Operator
Chapter 4 Directional Derivative and the Gradient Operator The equation z = f(x, y) defines a surface in 3 dimensions. We can write this as z f(x, y) = 0, or g(x, y, z) = 0, where g(x, y, z) = z f(x, y).
More informationCHAPTER 4 ELECTROMAGNETIC WAVES IN CYLINDRICAL SYSTEMS
CHAPTER 4 ELECTROMAGNETIC WAVES IN CYLINDRICAL SYSTEMS The vector Helmholtz equations satisfied by the phasor) electric and magnetic fields are where. In low-loss media and for a high frequency, i.e.,
More informationSOME PROBLEMS YOU SHOULD BE ABLE TO DO
OME PROBLEM YOU HOULD BE ABLE TO DO I ve attempted to make a list of the main calculations you should be ready for on the exam, and included a handful of the more important formulas. There are no examples
More informationChapter 3 - Vector Calculus
Chapter 3 - Vector Calculus Gradient in Cartesian coordinate system f ( x, y, z,...) dr ( dx, dy, dz,...) Then, f f f f,,,... x y z f f f df dx dy dz... f dr x y z df 0 (constant f contour) f dr 0 or f
More information(1) Recap of Differential Calculus and Integral Calculus (2) Preview of Calculus in three dimensional space (3) Tools for Calculus 3
Math 127 Introduction and Review (1) Recap of Differential Calculus and Integral Calculus (2) Preview of Calculus in three dimensional space (3) Tools for Calculus 3 MATH 127 Introduction to Calculus III
More informationIntroduction to the Phase Plane
Introduction to the Phase Plane June, 06 The Phase Line A single first order differential equation of the form = f(y) () makes no mention of t in the function f. Such a differential equation is called
More informationMathematical Concepts & Notation
Mathematical Concepts & Notation Appendix A: Notation x, δx: a small change in x t : the partial derivative with respect to t holding the other variables fixed d : the time derivative of a quantity that
More informationProblem set 7 Math 207A, Fall 2011 Solutions
Problem set 7 Math 207A, Fall 2011 s 1. Classify the equilibrium (x, y) = (0, 0) of the system x t = x, y t = y + x 2. Is the equilibrium hyperbolic? Find an equation for the trajectories in (x, y)- phase
More informationProblem Solving 1: The Mathematics of 8.02 Part I. Coordinate Systems
Problem Solving 1: The Mathematics of 8.02 Part I. Coordinate Systems In 8.02 we regularly use three different coordinate systems: rectangular (Cartesian), cylindrical and spherical. In order to become
More informationStudy Guide for Exam #2
Physical Mechanics METR103 November, 000 Study Guide for Exam # The information even below is meant to serve as a guide to help you to prepare for the second hour exam. The absence of a topic or point
More informationPractice Problems for Final Exam
Math 1280 Spring 2016 Practice Problems for Final Exam Part 2 (Sections 6.6, 6.7, 6.8, and chapter 7) S o l u t i o n s 1. Show that the given system has a nonlinear center at the origin. ẋ = 9y 5y 5,
More informationCURRENT MATERIAL: Vector Calculus.
Math 275, section 002 (Ultman) Spring 2012 FINAL EXAM REVIEW The final exam will be held on Wednesday 9 May from 8:00 10:00am in our regular classroom. You will be allowed both sides of two 8.5 11 sheets
More informationSolutions for the Practice Final - Math 23B, 2016
olutions for the Practice Final - Math B, 6 a. True. The area of a surface is given by the expression d, and since we have a parametrization φ x, y x, y, f x, y with φ, this expands as d T x T y da xy
More informationChapter 2. Vector Calculus. 2.1 Directional Derivatives and Gradients. [Bourne, pp ] & [Anton, pp ]
Chapter 2 Vector Calculus 2.1 Directional Derivatives and Gradients [Bourne, pp. 97 104] & [Anton, pp. 974 991] Definition 2.1. Let f : Ω R be a continuously differentiable scalar field on a region Ω R
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 informationMultiple Integrals and Vector Calculus: Synopsis
Multiple Integrals and Vector Calculus: Synopsis Hilary Term 28: 14 lectures. Steve Rawlings. 1. Vectors - recap of basic principles. Things which are (and are not) vectors. Differentiation and integration
More informationCDS 101 Precourse Phase Plane Analysis and Stability
CDS 101 Precourse Phase Plane Analysis and Stability Melvin Leok Control and Dynamical Systems California Institute of Technology Pasadena, CA, 26 September, 2002. mleok@cds.caltech.edu http://www.cds.caltech.edu/
More informationBrief Review of Vector Algebra
APPENDIX Brief Review of Vector Algebra A.0 Introduction Vector algebra is used extensively in computational mechanics. The student must thus understand the concepts associated with this subject. The current
More informationTraveling Harmonic Waves
Traveling Harmonic Waves 6 January 2016 PHYC 1290 Department of Physics and Atmospheric Science Functional Form for Traveling Waves We can show that traveling waves whose shape does not change with time
More informationStability of Dynamical systems
Stability of Dynamical systems Stability Isolated equilibria Classification of Isolated Equilibria Attractor and Repeller Almost linear systems Jacobian Matrix Stability Consider an autonomous system u
More informationDO NOT BEGIN THIS TEST UNTIL INSTRUCTED TO START
Math 265 Student name: KEY Final Exam Fall 23 Instructor & Section: This test is closed book and closed notes. A (graphing) calculator is allowed for this test but cannot also be a communication device
More informationEven-Numbered Homework Solutions
-6 Even-Numbered Homework Solutions Suppose that the matric B has λ = + 5i as an eigenvalue with eigenvector Y 0 = solution to dy = BY Using Euler s formula, we can write the complex-valued solution Y
More informationPractice Midterm Exam 1. Instructions. You have 60 minutes. No calculators allowed. Show all your work in order to receive full credit.
MATH202X-F01/UX1 Spring 2015 Practice Midterm Exam 1 Name: Answer Key Instructions You have 60 minutes No calculators allowed Show all your work in order to receive full credit 1 Consider the points P
More informationMath 210, Final Exam, Practice Fall 2009 Problem 1 Solution AB AC AB. cosθ = AB BC AB (0)(1)+( 4)( 2)+(3)(2)
Math 2, Final Exam, Practice Fall 29 Problem Solution. A triangle has vertices at the points A (,,), B (, 3,4), and C (2,,3) (a) Find the cosine of the angle between the vectors AB and AC. (b) Find an
More information1 + f 2 x + f 2 y dy dx, where f(x, y) = 2 + 3x + 4y, is
1. The value of the double integral (a) 15 26 (b) 15 8 (c) 75 (d) 105 26 5 4 0 1 1 + f 2 x + f 2 y dy dx, where f(x, y) = 2 + 3x + 4y, is 2. What is the value of the double integral interchange the order
More informationChapter 7. Nonlinear Systems. 7.1 Introduction
Nonlinear Systems Chapter 7 The scientist does not study nature because it is useful; he studies it because he delights in it, and he delights in it because it is beautiful. - Jules Henri Poincaré (1854-1912)
More information18.325: Vortex Dynamics
8.35: Vortex Dynamics Problem Sheet. Fluid is barotropic which means p = p(. The Euler equation, in presence of a conservative body force, is Du Dt = p χ. This can be written, on use of a vector identity,
More informationDepartment of Mathematics IIT Guwahati
Stability of Linear Systems in R 2 Department of Mathematics IIT Guwahati A system of first order differential equations is called autonomous if the system can be written in the form dx 1 dt = g 1(x 1,
More informationIntroduction and Vectors Lecture 1
1 Introduction Introduction and Vectors Lecture 1 This is a course on classical Electromagnetism. It is the foundation for more advanced courses in modern physics. All physics of the modern era, from quantum
More informationEE2007: Engineering Mathematics II Vector Calculus
EE2007: Engineering Mathematics II Vector Calculus Ling KV School of EEE, NTU ekvling@ntu.edu.sg Rm: S2-B2a-22 Ver: August 28, 2010 Ver 1.6: Martin Adams, Sep 2009 Ver 1.5: Martin Adams, August 2008 Ver
More informationASTR 320: Solutions to Problem Set 2
ASTR 320: Solutions to Problem Set 2 Problem 1: Streamlines A streamline is defined as a curve that is instantaneously tangent to the velocity vector of a flow. Streamlines show the direction a massless
More informationFaculty of Engineering, Mathematics and Science School of Mathematics
Faculty of Engineering, Mathematics and Science School of Mathematics SF Engineers SF MSISS SF MEMS MAE: Engineering Mathematics IV Trinity Term 18 May,??? Sports Centre??? 9.3 11.3??? Prof. Sergey Frolov
More information송석호 ( 물리학과 )
http://optics.hanyang.ac.kr/~shsong 송석호 ( 물리학과 ) Introduction to Electrodynamics, David J. Griffiths Review: 1. Vector analysis 2. Electrostatics 3. Special techniques 4. Electric fields in mater 5. Magnetostatics
More informationChapter 9 Uniform Circular Motion
9.1 Introduction Chapter 9 Uniform Circular Motion Special cases often dominate our study of physics, and circular motion is certainly no exception. We see circular motion in many instances in the world;
More information