Curl, Divergence, Gradient, and Laplacian in Cylindrical and Spherical Coordinate Systems

Similar documents
u 3 = u 3 (x 1, x 2, x 3 )

Section 11.6: Directional Derivatives and the Gradient Vector

MATH 1080 Test 2-SOLUTIONS Spring

FEM FOR HEAT TRANSFER PROBLEMS دانشگاه صنعتي اصفهان- دانشكده مكانيك

Objective Mathematics

First derivative analysis

Electromagnetic scattering. Graduate Course Electrical Engineering (Communications) 1 st Semester, Sharif University of Technology

That is, we start with a general matrix: And end with a simpler matrix:

Finite element discretization of Laplace and Poisson equations

BSc Engineering Sciences A. Y. 2017/18 Written exam of the course Mathematical Analysis 2 August 30, x n, ) n 2

Direct Approach for Discrete Systems One-Dimensional Elements

1973 AP Calculus AB: Section I

MATHEMATICS PAPER IIB COORDINATE GEOMETRY AND CALCULUS. Note: This question paper consists of three sections A, B and C.

Mathematical Notation Math Calculus & Analytic Geometry III

2. Background Material

Background: We have discussed the PIB, HO, and the energy of the RR model. In this chapter, the H-atom, and atomic orbitals.

MAXIMA-MINIMA EXERCISE - 01 CHECK YOUR GRASP

Mock Exam 2 Section A

2F1120 Spektrala transformer för Media Solutions to Steiglitz, Chapter 1

Function Spaces. a x 3. (Letting x = 1 =)) a(0) + b + c (1) = 0. Row reducing the matrix. b 1. e 4 3. e 9. >: (x = 1 =)) a(0) + b + c (1) = 0

Lecture 2: Review of Vector Calculus

Announce. ECE 2026 Summer LECTURE OBJECTIVES READING. LECTURE #3 Complex View of Sinusoids May 21, Complex Number Review

MATHEMATICS PAPER IB COORDINATE GEOMETRY(2D &3D) AND CALCULUS. Note: This question paper consists of three sections A,B and C.

Mathematical Notation Math Calculus & Analytic Geometry III

EELE 3331 Electromagnetic I Chapter 3. Vector Calculus. Islamic University of Gaza Electrical Engineering Department Dr.

10. The Discrete-Time Fourier Transform (DTFT)

Appendix: Orthogonal Curvilinear Coordinates. We define the infinitesimal spatial displacement vector dx in a given orthogonal coordinate system with

Derivation of Electron-Electron Interaction Terms in the Multi-Electron Hamiltonian

2008 AP Calculus BC Multiple Choice Exam

Engineering Mathematics I. MCQ for Phase-I

The Matrix Exponential

PHYS-333: Problem set #2 Solutions

The Matrix Exponential

Solution: APPM 1360 Final (150 pts) Spring (60 pts total) The following parts are not related, justify your answers:

Deepak Rajput

1997 AP Calculus AB: Section I, Part A

Hydrogen Atom and One Electron Ions

Gradient, Divergence and Curl in Curvilinear Coordinates

1 1 1 p q p q. 2ln x x. in simplest form. in simplest form in terms of x and h.

2.3 Matrix Formulation

surface of a dielectric-metal interface. It is commonly used today for discovering the ways in

THE SPINOR FIELD THEORY OF THE PHOTON

8. THEOREM If the partial derivatives f x. and f y exist near (a, b) and are continuous at (a, b), then f is differentiable at (a, b).

Calculus II (MAC )

Mathematical Concepts & Notation

Math 102. Rumbos Spring Solutions to Assignment #8. Solution: The matrix, A, corresponding to the system in (1) is

Mat 267 Engineering Calculus III Updated on 9/19/2010

SAFE HANDS & IIT-ian's PACE EDT-15 (JEE) SOLUTIONS

NEW APPLICATIONS OF THE ABEL-LIOUVILLE FORMULA

Introduction to Condensed Matter Physics

Elements of Statistical Thermodynamics

COHORT MBA. Exponential function. MATH review (part2) by Lucian Mitroiu. The LOG and EXP functions. Properties: e e. lim.

Division of Mechanics Lund University MULTIBODY DYNAMICS. Examination Name (write in block letters):.

SME 3033 FINITE ELEMENT METHOD. Bending of Prismatic Beams (Initial notes designed by Dr. Nazri Kamsah)

Chapter 3 - Vector Calculus

1997 AP Calculus AB: Section I, Part A

Addition of angular momentum

u r du = ur+1 r + 1 du = ln u + C u sin u du = cos u + C cos u du = sin u + C sec u tan u du = sec u + C e u du = e u + C

Partial Derivatives: Suppose that z = f(x, y) is a function of two variables.

Exponential Functions

EAcos θ, where θ is the angle between the electric field and

Addition of angular momentum

MATHEMATICS (B) 2 log (D) ( 1) = where z =

AS 5850 Finite Element Analysis

Middle East Technical University Department of Mechanical Engineering ME 413 Introduction to Finite Element Analysis

ANSWER KEY 1. [A] 2. [C] 3. [B] 4. [B] 5. [C] 6. [A] 7. [B] 8. [C] 9. [A] 10. [A] 11. [D] 12. [A] 13. [D] 14. [C] 15. [B] 16. [C] 17. [D] 18.

INC 693, 481 Dynamics System and Modelling: The Language of Bound Graphs Dr.-Ing. Sudchai Boonto Assistant Professor

Einstein Equations for Tetrad Fields

Why is a E&M nature of light not sufficient to explain experiments?

A Spatial Version of the Theorem of the Angle of Circumference

Trusses. Introduction. Plane Trusses Local and global coordinate systems. Plane Trusses Local and global coordinate systems

Vector Calculus. A primer

ENGR 323 BHW 15 Van Bonn 1/7

Chapter 5. Introduction. Introduction. Introduction. Finite Element Modelling. Finite Element Modelling

Introduction to Vector Calculus (29) SOLVED EXAMPLES. (d) B. C A. (f) a unit vector perpendicular to both B. = ˆ 2k = = 8 = = 8

CHAPTER 1. Introductory Concepts Elements of Vector Analysis Newton s Laws Units The basis of Newtonian Mechanics D Alembert s Principle

Ch. 24 Molecular Reaction Dynamics 1. Collision Theory

Class 8 Multiplication of Polynomials

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.

S.Y. B.Sc. (IT) : Sem. III. Applied Mathematics. Q.1 Attempt the following (any THREE) [15]

MCE503: Modeling and Simulation of Mechatronic Systems Discussion on Bond Graph Sign Conventions for Electrical Systems

a b ixā + y b ixb + ay

Contents. MATH 32B-2 (18W) (L) G. Liu / (TA) A. Zhou Calculus of Several Variables. 1 Multiple Integrals 3. 2 Vector Fields 9

Higher-Order Discrete Calculus Methods

DSP-First, 2/e. LECTURE # CH2-3 Complex Exponentials & Complex Numbers TLH MODIFIED. Aug , JH McClellan & RW Schafer

cycle that does not cross any edges (including its own), then it has at least

Vectors and Fields. Vectors versus scalars

송석호 ( 물리학과 )

Calculus Revision A2 Level

Answers & Solutions. for MHT CET-2018 Paper-I (Mathematics) Instruction for Candidates

nd the particular orthogonal trajectory from the family of orthogonal trajectories passing through point (0; 1).

PHYSICS 489/1489 LECTURE 7: QUANTUM ELECTRODYNAMICS

5. B To determine all the holes and asymptotes of the equation: y = bdc dced f gbd

송석호 ( 물리학과 )

Chapter 4: Partial differentiation

Unfired pressure vessels- Part 3: Design

Chapter 10. The singular integral Introducing S(n) and J(n)

MAT 203 Vector Calculus [Kalkulus Vektor]

Basic Polyhedral theory

Review of Linear Algebra

Transcription:

A P P E N D I X B Curl, Divrgnc, Gradint, and Laplacian in Cylindrical and Sphrical Coordinat Systms 788 In Chaptr 3, w introducd th curl, divrgnc, gradint, and Laplacian and drivd th xprssions or thm in th Cartsian coordinat systm. In this appndix, w driv th corrsponding xprssions in th cylindrical and sphrical coordinat systms. Considring irst th cylindrical coordinat systm, w rcall rom Sction.3 that th ininitsimal box dind by th thr orthogonal suracs intrscting at point Pr,, z2 and th thr orthogonal suracs intrscting at point Qr + dr, + d, z + dz2 is as shown in Fig. B.. From th basic dinition o th curl o a vctor introducd in Sction 3.3 and givn by A c AC d a (B.) S: S n max w ind th componnts o Aas ollows, with th aid o Fig. B.: A2 r = d: dz: [A z ] r, + d2 - [A z ] r,2 d: r d Aabcda lim d: ara abcd dz: A z [A ] r,z2 r d + [A z ] r, + d2 dz - [A ] r,z + dz2 r d - [A z ] r,2 dz - A z r d dz dz: [A ] r,z2 - [A ] r,z + dz2 dz (B.2a)

Appndix B 789 c d h Q(r dr, d, z dz) dz b a P(r,, z) dr r d g (r dr) d FIGURE B. Ininitsimal box ormd by incrmnting th coordinats in th cylindrical coordinat systm. Aada A2 dz: ara ad dz: [A z] r,2 dz + [A r ], z + dz2 dr -[A z ] r + dr,2 dz - [A r ],z2 dr dr dz (B.2b) dz: [A r ],z + dz2 - [A r ],z2 dz [A z ] r,2 - [A z ] r + dr,2 dr = A r z - A z r Aagba A2 z ara agb d: [A r],z2 dr + [A ] r + dr,z2 r + dr2 d - [A r ] + d,z2 dr - [A ] r,z2 r d r dr d d: (B.2c) [ra ] r + dr,z2 - [ra ] r,z2 r dr r ra 2 - A r r [A r ],z2 - [A r ] + d,z2 d: r d

79 Appndix B Curl, Divrgnc, Gradint, and Laplacian Combining (B.2a), (B.2b), and (B.2c), w obtain th xprssion or th curl o a vctor in cylindrical coordinats as A = c r A z - A z da r + c A r z + r c r ra 2 - A r da z a r a z a r r = 5 5 r z - A z r da (B.3) A r ra A z To ind th xprssion or th divrgnc, w us th basic dinition o th divrgnc o a vctor, introducd in Sction 3.3 and givn by # A = AS A # ds lim v: v (B.4) Evaluating th right sid o (B.4) or th box o Fig. B., w obtain # [A r] r + dr r + dr2 d dz - [A r ] r r d dz + [A ] + d dr dz - [A ] dr dz + [A z ] z + dz r dr d - [A z ] z r dr d A r dr d dz d: dz: [ra r ] r + dr - [ra r ] r r dr [A z ] z + dz - [A z ] z dz: dz A r ra r2 + r + A z z [A ] + d - [A ] d: r d (B.5) To obtain th xprssion or th gradint o a scalar, w rcall rom Sction.3 that in cylindrical coordinats, Thror, dl = dr a r + r d a + dz a z (B.6) d = r dr + d + z dz = a r a r + r a + z a zb # dr ar + r d a + dz a z 2 = # dl (B.7)

Appndix B 79 Thus, = r a r + r a + z a z (B.8) To driv th xprssion or th Laplacian o a scalar, w rcall rom Sction 5. that 2 = # (B.9) Thn using (B.5) and (B.8), w obtain 2 r ar r b + r a r b + z a z b r ar r b + 2 r 2 2 + 2 z 2 (B.) Turning now to th sphrical coordinat systm, w rcall rom Sction.3 that th ininitsimal box dind by th thr orthogonal suracs intrscting at Pr, u, 2 and th thr orthogonal suracs intrscting at Qr + dr, u + du, + d2 is as shown in Fig. B.2. From th basic dinition o th curl o a vctor givn by (B.), w thn ind th componnts o Aas ollows, with th aid o Fig. B.2: d (r dr) sin u d d a P(r, u, ) dr c (r dr) du h Q(r dr, u du, d) r du r sin (u du) d FIGURE B.2 b g Ininitsimal box ormd by incrmnting th coordinats in th sphrical coordinat systm.

792 Appndix B Curl, Divrgnc, Gradint, and Laplacian Aabcda A2 r du: ara abcd d: du: d: A # Aada dl A2 u d: ara ad d: [A u] r,2 r du + [A ] r,u + du2 r sinu + du2 d - [Au] r, + d2r du - [A ] r,u2 d [A r ] u, + d2 - [A r ] u, 2 d: d r 2 sin u du d [A sin u] r,u + du2 - [A sin u] r,u2 du: du [A u ] r,2 - [A u ] r, + d2 d: d = u A sin u2 - A u [A ] r,u2 d + [A r ] u, + d2 dr -[A ] r + dr,u2 r + dr2 sin u d - [A r ] u,2 dr dr d (B.a) (B.b) = [ra ] r,u2 - [ra ] r + dr, u2 A r - r r dr r ra 2 Aagba A2 ara agb du: du: [A r] u,2 dr + [A u ] r + dr,2 r + dr2 du -[A r ] u + du,2 dr - [A u ] r,2 r du r dr du [ra u ] r + dr,2 - [ra u ] r,2 r dr [A r ] u,2 dr - [A r ] u + du,2 dr du: r du = r r ra u2 - A r r u (B.c)

Combining (B.a), (B.b), and (B.c), w obtain th xprssion or th curl o a vctor in sphrical coordinats as Appndix B 793 A = c u A sin u2 - A u da r + r c A r sin u - r ra 2da u + r c r ra u2 - A r u da (B.2) a r r 2 sin u = 5 r a u u a r 5 A r ra u A To ind th xprssion or th divrgnc, w us th basic dinition o th divrgnc o a vctor givn by (B.4), and by valuating its right sid or th box o Fig. B.2, w obtain # A du: d: [A r ] r + dr r + dr2 2 sin u du d - [A r ] r r 2 sin u du d c + [A u ] u + du r sinu + du2 dr d - [Au]u dr d s + [A ] + d r dr du - [A ] r dr du [r 2 A r ] r + dr - [r 2 A r ] r r 2 dr [A ] + d - [A ] d: d 2 r r2 A r 2 + r 2 sin u dr du d du: [A u sin u] u + du - [A u sin u] u du u A u sin u2 + A (B.3) To obtain th xprssion or th gradint o a scalar, w rcall rom Sction.3 that in sphrical coordinats, dl = dr a r + r du a u + d a (B.4)

794 Appndix B Curl, Divrgnc, Gradint, and Laplacian Thror, d = r = # dl dr + u du + d = a r a r + r u a u + a b # dr ar + r du a u + d a 2 (B.5) Thus, = r a r + r u a u + a (B.6) To driv th xprssion or th Laplacian o a scalar, w us (B.9), in conjunction with (B.3) and (B.6). Thus, w obtain 2 = r 2 ar2 r r b + u a r + a = r 2 ar2 r r b + 2 + r 2 sin 2 u 2 b r 2 sin u u sin ub u asin u u b (B.7)