arxiv: v3 [math-ph] 24 Aug 2010

Transcription

1 Inverse Vector Operators Shaon Sahoo Department of Physics, Indian Institute of Science, Bangalore 5600, India. arxiv:0804.9v [math-ph] 4 Aug 00 Abstract In different branches of physics, we frequently deal with vector del operator ( ). This del operator is generally used to find curl or divergence of a vector function or gradient of a scalar function. In many important cases, we need to know the parent vector whose curl or divergence is known or require to find the parent scalar function whose gradient is known. But the task is not very easy, especially in case of finding vector potential whose curl is known. Here, inverse curl, inverse divergence and inverse gradient operators are defined to solve those problems easily. Keywords del operator, curl, gradient, divergence, solenoidal vector, vector potentials, curvilinear coordinate system. Introduction In physical science, usefulness of the concepts of vector and scalar potentials is indisputable. They can be used to replace more conventional concepts of magnetic and electric fields. In fact, in modern era of science, these concepts of potentials became so useful and popular that their presence can be seen in virtually all disciplines related somehow to electromagnetism! The operator which relates those potentials to the conventional fields is the so-called vector del operator ( ). This operator plays a very important role in a wide range of physics besides being indispensable in electromagnetism. As for example, it is now very common to see the operator in fluid dynamics, quantum mechanics, statistical mechanics, classical mechanics, theory of relativity, thermodynamics, etc. There are three basic operations, namely, curl, divergence and gradient, which can be performed by the operator to relate different quantities of importance. For the ease of discussion, let us divide the functions in two categories, (a) potential functions (this can be vector or scalar) and (b) field functions (again can be vector or scalar). Above mentioned operations are performed on potential functions to obtain corresponding field functions. While these operations are easy to perform, a back operation to get a potential function from a given field function is not always trivial. Since in some cases, it is convenient to replace one type of functions with others, one should be able to interchange these two types of functions comfortably. Unfortunately, a potential function is not unique for a given field function and there is no general and unique procedure to obtain a potential function from a given field function. In fact, due to this non-uniqueness, one may wonder whether it is at all possible to find any general procedure for this purpose. It is exactly the concern addressed in this article. Here it is shown that, it is possible to find some general procedures in the form of inverse vector operators, which can be applied easily to the field functions to obtain the corresponding potential functions. shaon@physics.iisc.ernet.in

2 All the inverse operators are defined here in the orthogonal curvilinear co-ordinate system []. Let u, u, u be the three mutually perpendicular families of co-ordinate surfaces and ê, ê, ê as unit vectors normal to respective surfaces. Inverse Curl Operator At first we mention some symbols, used later. a. ( ) indicates inverse curl operator. b. ( + + )(u ) does integration w.r.t. u, while other two variables of the integrand are to be treated fixed. If the integrand is the rd component (along ê ) of a vector, then + (u ) acts only on the part of the integrand having the variable u. Similarly (u ) acts only on the part of the integrand does not contain the variable u. For example, let φ(u,u,u ) Ψ (u,u )+Ψ (u,u,u ) be the component of a vector along ê, then, (p + +p )φu p + Ψ u + p Ψ u, with p and p being any two numbers (u,u are to be treated as constants). Now take two vectors AA ê + A ê + A ê and BB ê + B ê + B ê, such that, B A () Using inverse curl operator, this vector potential A can be expressed as, From eqn () it is obvious that, A( ) B () B ( A)0 [ h B )+ h B )+ h B )]0 h h h u u u u [c (u + )+c (u )]+ u [c (u + )+c (u )]+ u [c (u + )+c (u )]0 c (u + u )+ c (u + u )+ c (u + )0 () u where c (u + ) and c (u ) are two parts of the term h h B, containing and not containing u respectively, such that, h h B c (u + )+c (u ). Other terms can be defined similarly. From eqn () we get, B A h h h h ê h ê h ê u u u h A h A h A (4)

3 B ê + B ê + B ê ê [ A ) A )] h h u u + ê + ê [ A ) A )] h h u u [ A ) A )] (5) h h u u We already have defined components of B in terms of c s. Now, we have to choose h A, h A and h A in terms of c s in such a way that, r.h.s. of the eqn (5) becomes same as l.h.s. Let (by inspection), h A k c (u + )u + k c (u + )u + k c (u )u + k 4 c (u )u Similarly, h A k 5 c (u + )u + k 6 c (u + )u + k 7 c (u )u + k 8 c (u )u and, h A k 9 c (u + )u + k 0 c (u + )u + k c (u )u + k c (u )u where, k,k,,k are some constants to be determined. (Note that, c terms are not included for h A, since, h A does not appear in the expression of st component B in eqn (5). Also note that, terms in the expression of h A are in the integral form w.r.t. variables u and u but not u, as in eqn (5), h A is never differentiated w.r.t. u.) Now, see the coefficient of ê (expression of B ) on the r.h.s. of eqn (5), h h [ u A ) u A )] h h [ u {k 5 c (u + )u + k 6 c (u + )u + k 7 c (u )u + k 8 c (u )u } u {k c (u + )u + k c (u + )u + k c (u )u + k 4 c (u )u }] h h [k 5 c (u + )+k 6 u c (u + )u + k 7 c (u )+k 8 0 k u c (u + )u k c (u + ) k 0 k 4 c (u )] h h [k 5 c (u + ) k { u c (u + )+ u c (u + )}u + k 7 c (u )] [Taking k 5 k, k 7 k 4 and k 6 k ] h h [k 5 c (u + ) k { u c (u + )}u + k 7 c (u )] [Using eqn ()] h h [k 5 c (u + )+k c (u + )+k 7c (u )] h h [k 5 c (u + )+k 7c (u )] [Taking k k 5 ]

4 h h [ c (u + )+ c (u )] [Considering k 5 / and k 7 /] h h [c (u + )+c (u )] h h h B ) B Similarly we can show that, if, k k 5 k 9 k k 6 k 0 / and, k k 7 k k 4 k 8 k /, then, h h [ u A ) u A )]B and h h [ u A ) u A )]B So, using these values of k,k,,k, we get, A h [ c (u + )u c (u + )u ]+ h [ c (u )u c (u )u ] A h [ c (u + )u c (u + )u ]+ h [ c (u )u c (u )u ] A h [ c (u + )u c (u + )u ]+ h [ c (u )u c (u )u ] Now from eqn (), we get, ( ) B AA ê + A ê + A ê ê h [ c (u + )u c (u + )u ]+ ê h [ c (u )u c (u )u ] + ê h [ c (u + )u c (u + )u ]+ ê h [ c (u )u c (u )u ] + ê h [ c (u + )u c (u + )u ]+ ê h [ c (u )u c (u )u ] ê ê (u ) (u ) h c (u + ) c (u + ) + ê (u ) (u ) h c (u + ) c (u + ) ê (u ) (u ) h c (u + ) c (u + ) (u ) (u ) h c (u ) c (u ) + ê (u ) (u ) h c (u ) c (u ) ê (u ) (u ) h c (u ) c (u ) ê ê ê h h h ê ê ê h h h (u ) (u ) (u ) (u ) (u ) (u ) c (u + ) c (u + ) c (u + ) c (u ) c (u ) c (u ) 4

5 ê ê ê h h h (u ) (u ) (u ) c (u + )+ c (u ) c (u + )+ c (u ) c (u + )+ c (u ) ê ê ê h h h ( + + )(u ) ( + + )(u ) ( + + )(u ) h h B h h B h h B This is the desired expression for our inverse curl operator. Since, ( φ)0, so, it should be noted that, for a given solenoidal vector, corresponding vector potential is not unique, which can be in general expressed as, - A( ) B+ φ where φ is an arbitrary scalar function.. An Example in Cartesian System One example in Cartesian system would clarify the procedure of getting vector potential for a solenoidal vector. In this system, ê, ê, ê, u, u and u are respectively replaced by î, ĵ, ˆk, x, y, and z. Here h h h. Now take a solenoidal vector B(xyz+y )î+(xz+y) ĵ z(+yz/)ˆk. (note B0) So, vector potentials of ˆB are,- A( ) B+ φ î ĵ ˆk ( + + )(x) ( + + )(y) ( + + )(z) xyz+y xz+y z yz / + φ î[ ( z yz /)y+ 0y yz xzz] ĵ[ xyzz+ y z ( z yz /)x 0x] ˆk[ yx+ xzx xyzy y y]+ φ î[(zy+z y /4)/+(yz/+xz /4)] ĵ[(zx+yxz /)/+(xyz /6+y z/)] + ˆk[ (yx/+x z/4)+xy z/6+y /6]+ φ 5

6 .. Verification of the Result î ĵ ˆk A x y z + ( φ) (zy+z y /4)/ (zx+xyz /)/ (yx/+x z/4) +(yz/+xz /4) (xyz /6+y z/) +(xy z/6+y /6) î(xyz+y )+ ĵ(xz+y) ˆkz(+yz/)+0 B Inverse Divergence Operator Here we indicate this operator by ( ). Now if A and φ be a vector and a scalar function respectively, such that φ A h h h [ u h A )+ u h A )+ u h A )], then A can be expressed as A( ) φ. We can now define the operator as (by inspection), ( ) k ê h h h h h (u )+k ê h h h h h (u )+k ê h h h h h (u ) where, k + k + k.. Verification For a given φ, A( ) φ k ê h h h h h φ(u )+k ê h h h h h φ(u )+k ê h h h h h φ(u ) Now note that, divergence of the obtained vector gives the scalar function φ, A h h h [k h h h φ+k h h h φ+k h h h φ] (k + k + k )φφ Since, ( B) 0, so, it should be noted that for a given scalar function, corresponding potential function is not unique, which can be in general expressed as, A( ) φ+ B where B is any vector. 6

7 4 Inverse Gradient Operator We indicate this operator by ( ). Let AA ê +A ê +A ê and φ(u,u,u ) be a vector and a scalar function, such that A φ. Since, the vector field represented by A is conservative, thus, line integral of A is path independent. Choose (a,b,c) as initial point, such that φ(a,b,c) exists. Now φ(u,u,u ) (u,u,u ) dφ (a,b,c) + c 0 [c 0 is integration constant] (u,u,u ) ( φ) dr (a,b,c) (u,u,u ) A dr (a,b,c) + c 0 + c 0 (a,b,u ) A h du (a,b,c) + (a,u,u ) A h du (a,b,u ) + (u,u,u ) A h du (a,u,u ) + c 0 So, φ( ) A u a u u A h u + A h u + A h u + c 0 b c u u a a u b Here in the first integration u is variable but other two (u,u ) are fixed, and so on for other two integrations. 7

8 Note: The original work was done in 00 and published in 006 in Science and Culture []. This write-up is almost same as published one. References [] Any text on vector del operator and curvilinear co-ordinate system. For example: by George B. Arfken and Hans J. Weber, Mathematical Methods For Physicists, Academic Press: 5th edition (000); st and nd chapters. [] Shaon Sahoo, Sci. & Cult., 7, (-), 89 (006). 8