*Notes to a video lecture on http://www.unizor.com*

__Divergence__

In this lecture we continue studying the characteristics of

*fields*and, in particular, discuss another usage of symbol

**∇**(

*nabla*) that represents a set (or a triplet, or a pseudo-vector, or simply "vector") of three operators of partial differentiation, performed for each of the three dimensions of space we live in, in given Cartesian coordinates.

If needed, review the topic of partial derivatives in the "Calculus" part of the course "Math 4 Teens" on UNIZOR.COM.

We use the symbol

**∇**as a shorthand for a triplet of operators of partial differentiation

**∇**= {

*∂*,

**/**∂**x***∂*,

**/**∂**y***∂*}

**/**∂**z**Previous lecture about

*gradient*of a

*scalar field*

*defined*

**F(x,y,z)****∇**as a "multiplication" of a "vector"

*F(x,y,z)***∇**by "scalar"

*, resulting in a "vector", as traditional multiplication of a vector by a constant should:*

**F(x,y,z)****∇**

*F(x,y,z) =*

={=

*∂*,

**F()/**∂**x***∂*,

**F()/**∂**y***∂*}

**F()/**∂**z**The physical meaning of the

*gradient*at any point {

*} of a scalar field*

**x,y,z***was a*

**F(x,y,z)***vector field*(direction and magnitude at each point) of the greatest change of this field, thereby defining from a

*scalar field*

*a*

**F(x,y,z)***vector field*

**∇**.

*F(x,y,z)*As an application of this concept, if

*represents an air density at point {*

**F(x,y,z)***}, which is a scalar, the*

**x,y,z****∇**represents the direction and magnitude of the most significant potential movement of the air, a vector, the main wind direction and its strength at each point.

*F(x,y,z)*The

*divergence*is the result of a different operation with "vector"

**∇**.

Consider, again, an example of air distribution with different air density (and, therefore, pressure) in different areas of space. As a result, we can observe the direction and strength of a wind at any point - a

*vector field*of air velocities (speed and direction of a wind) at each point.

The air goes from zones of higher density towards zones of lower density. While this process goes, the amount of air in a zone of higher density diminishes, while the amount of air in a zone of lower density increases until the density (and, therefore, pressure) equalizes. During this process the amount of air in a zone of average density in between the above two zones will remain steady, what air comes from a high density zone will go to a lower density.

The

*divergence*of a vector field (for example, of a vector field of velocities of the wind at each point in space) quantitatively reflects the situation with air zones described above, showing the areas, from which air flows from and those it flows into.

*Divergence*is a scalar defined for each point of a

*vector field*. It's positive for those areas, from which something (like air, or electrons, or energy) is flowing out, negative for areas receiving that something and zero if the amounts of coming into and going out are the same.

Let's approach this mathematically.

For each point {

*} in space there is vector*

**x,y,z***with projections on coordinate axes*

**V(x,y,z)***,*

**V**_{x}(x,y,z)*and*

**V**_{y}(x,y,z)*so that the following vector equality is held:*

**V**_{z}(x,y,z)

**V(x,y,z) = V**

+ V_{x}(x,y,z)·i ++ V

_{y}(x,y,z)·j + V_{z}(x,y,z)·kwhere

*,*

**i***and*

**j***are unit vectors along X-, Y- and Z-axis.*

**k**Consider a parallelepiped

*centered at point*

**ABCDEFGH***and infinitesimal dimensions Δ*

**P(x,y,z)***⨯Δ*

**x***⨯Δ*

**y***.*

**z**Then coordinates of vertices are

*{*

**A***Δ*

**x−½***Δ*

**x,y−½***Δ*

**y,z−½***}*

**z***{*

**B***Δ*

**x+½***Δ*

**x,y−½***Δ*

**y,z−½***}*

**z***{*

**C***Δ*

**x+½***Δ*

**x,y+½***Δ*

**y,z−½***}*

**z***{*

**D***Δ*

**x−½***Δ*

**x,y+½***Δ*

**y,z−½***}*

**z***{*

**E***Δ*

**x−½***Δ*

**x,y−½***Δ*

**y,z+½***}*

**z***{*

**F***Δ*

**x+½***Δ*

**x,y−½***Δ*

**y,z+½***}*

**z***{*

**G***Δ*

**x+½***Δ*

**x,y+½***Δ*

**y,z+½***}*

**z***{*

**H***Δ*

**x−½***Δ*

**x,y+½***Δ*

**y,z+½***}*

**z**To bring some flavor into mathematics, let's use our analogy of the

*vector field*

*as a direction and magnitude of the wind at point*

**V(x,y,z)***.*

**P(x,y,z)**To evaluate the air contribution by this wind into point

*, we will calculate how much air coming in and going out of the parallelepiped*

**P(x,y,z)***per unit of time and take its limit, as the dimensions of this parallelepiped are infinitesimal.*

**ABCDEFGH**We will consider three major directions along the three axes of the coordinate system and coordinate representation of the vector

*.*

**V(x,y,z)**Projection of the vector

*on the X-axis*

**V(x,y,z)***brings air into our parallelepiped through its left side and takes the air out through its right side.*

**V**_{x}(x,y,z)·iIn the middle of the left side this projection is

*Δ*

**V**_{x}(x−½*.*

**x,y,z)·i**In the middle of the right side this projection is

*Δ*

**V**_{x}(x+½*.*

**x,y,z)·i**Since we assume that our parallelepiped has infinitesimal dimensions, we can assume that the projection of the vector

*on the X-axis at each point of the left side*

**V(x,y,z)***is the same as in the middle of this side at point {*

**ADHE***Δ*

**x−½***}.*

**x,y,z**Similarly, we assume that the projection of the vector

*on the X-axis at each point of the right side*

**V(x,y,z)***is the same as in the middle of this side at point {*

**BCGF***Δ*

**x+½***}.*

**x,y,z**The amount of air moving out through the right side of a parallelepiped is proportional to the average speed of wind at the right side

*Δ*

**V**_{x}(x+½*times the area of the right side Δ*

**x,y,z)***·Δ*

**y***.*

**z**Similarly, the amount of air moving in through the left side of a parallelepiped is proportional to the average speed of wind at the left side

*Δ*

**V**_{x}(x−½*times the area of the left side Δ*

**x,y,z)***·Δ*

**y***.*

**z**The net result of air moving in and out of the parallelepiped along the X-axis is proportional to

*[*

**A**_{x}(x,y,z) =*Δ*

**V**_{x}(x+½*Δ*

**x,y,z) −**

− V− V

_{x}(x−½*]·Δ*

**x,y,z)***·Δ*

**y**

**z**Recall from the calculus the following property of differentiable function

*:*

**f(x)**

**f(b) − f(a) =**d**f(c)/**d**x · (b−a)**where point

*is somewhere between*

**c***and*

**a***.*

**b**If

*→*

**a***(that is,*

**b***is an infinitesimal variable),*

**b−a***also tends to*

**c***(or*

**a***) and*

**b***d*.

**f(c)/**d**x ≅**d**f(a)/**d**x ≅**

≅d≅

**f(b)/**d**x ≅**d**f(½(a+b))/**d**x**Using this property in our case with

*Δ*

**a=x−***and*

**x***Δ*

**b=x+***, with Δ*

**x***infinitesimal, we derive*

**x***Δ*

**V**_{x}(x+½*Δ*

**x,y,z) −**

− V− V

_{x}(x−½*Δ*

**x,y,z) ≅**

≅∂≅

**V**∂_{x}(x,y,z)/**x ·**

**x**Now the expression for an amount of air accumulated in the vicinity of point {

*} is*

**x,y,z***Δ*

**A**

≅∂_{x}(x,y,z) ≅≅

**V**∂_{x}(x,y,z)/**x ·***Δ*

**x·***Δ*

**y·**

**z**Since Δ

*Δ*

**x·***Δ*

**y·***is the volume of our parallelepiped, dividing both sides of the equation above by it will give a change (positive for increment or negative for decrement) of air density in the infinitesimal neighborhood of point {*

**z***}, contributed by a component of the vector field in the X-axis direction.*

**x,y,z**

**D**∂_{x}(x,y,z) ≅**V**∂_{x}(x,y,z)/**x**Analogously, contribution to air density by Y-component of the vector field

*is*

**V(x,y,z)**

**D**∂_{y}(x,y,z) ≅**V**∂_{y}(x,y,z)/**y**And Z-component of this vector field adds

**D**∂_{z}(x,y,z) ≅**V**∂_{z}(x,y,z)/**z**The combined contribution of air density will be a sum of these individual contributions, a

**divergence**(

*div*) of the vector field

*:*

**V(x,y,z)***div*

**V(x,y,z) =**∂**V**∂_{x}(x,y,z)/**x +**

+∂+

**V**∂_{y}(x,y,z)/**y +**∂**V**∂_{z}(x,y,z)/**z**To simplify this expression,

**∇**to the rescue.

Considering

**∇**as a "vector"

**∇**= {

*∂*}

**/**∂**x,**∂**/**∂**y,**∂**/**∂**z**and

*being a vector in coordinate representation*

**V(x,y,z)**{

*}*

**V**_{x}(x,y,z),V_{y}(x,y,z),V_{z}(x,y,z)the expression for a

**divergence**above can be represented as a scalar (dot) product of these two "vectors":

*div V(x,y,z) = *

**∇**

*· V(x,y,z)*Other examples of the usage of a concept

*divergence*are: divergence of the electric field (not zero where there are electric charges, see the First Maxwell equation),

divergence of the magnetic field (always zero, see the Second Maxwell equation),

divergence of the gravitational field (non-zero inside any material object, zero in vacuum),

and any others, where it makes sense to talk about

*vector field*.

## No comments:

Post a Comment