Scalar and Vector Field Functionality

Introduction

Vectors and Scalars

In physics, we deal with two kinds of quantities – scalars and vectors.

A scalar is an entity which only has a magnitude – no direction. Examples of scalar quantities include mass, electric charge, temperature, distance, etc.

A vector, on the other hand, is an entity that is characterized by a magnitude and a direction. Examples of vector quantities are displacement, velocity, magnetic field, etc.

A scalar can be depicted just by a number, for e.g. a temperature of 300 K. On the other hand, vectorial quantities like acceleration are usually denoted by a vector. Given a vector \(\mathbf{V}\), the magnitude of the corresponding quantity can be calculated as the magnitude of the vector itself \(\Vert \mathbf{V} \Vert\), while the direction would be specified by a unit vector in the direction of the original vector, \(\mathbf{\hat{V}} = \frac{\mathbf{V}}{\Vert \mathbf{V} \Vert}\).

For example, consider a displacement of \((3\mathbf{\hat{i}} + 4\mathbf{\hat{j}} + 5\mathbf{\hat{k}})\) m, where , as per standard convention, \(\mathbf{\hat{i}}\), \(\mathbf{\hat{j}}\) and \(\mathbf{\hat{k}}\) represent unit vectors in the \(\mathbf{X}\), \(\mathbf{Y}\) and \(\mathbf{Z}\) directions respectively. Therefore, it can be concluded that the distance traveled is \(\Vert 3\mathbf{\hat{i}} + 4\mathbf{\hat{j}} + 5\mathbf{\hat{k}} \Vert\) m = \(5\sqrt{2}\) m. The direction of travel is given by the unit vector \(\frac{3}{5\sqrt{2}}\mathbf{\hat{i}} + \frac{4}{5\sqrt{2}}\mathbf{\hat{j}} + \frac{5}{5\sqrt{2}}\mathbf{\hat{k}}\).

Fields

In general, a \(field\) is a vector or scalar quantity that can be specified everywhere in space as a function of position (Note that in general a field may also be dependent on time and other custom variables). In this module, we deal with 3-dimensional spaces only. Hence, a field is defined as a function of the \(x\), \(y\) and \(z\) coordinates corresponding to a location in 3D space.

For example, temperate in 3 dimensional space (a temperature field) can be written as \(T(x, y, z)\) – a scalar function of the position. An example of a scalar field in electromagnetism is the electric potential.

In a similar manner, a vector field can be defined as a vectorial function of the location \((x, y, z)\) of any point in space.

For instance, every point on the earth may be considered to be in the gravitational force field of the earth. We may specify the field by the magnitude and the direction of acceleration due to gravity (i.e. force per unit mass ) \(g(x, y, z)\) at every point in space.

To give an example from electromagnetism, consider an electric potential of form \(2{x}^{2}y\), a scalar field in 3D space. The corresponding conservative electric field can be computed as the gradient of the electric potential function, and expressed as \(4xy\mathbf{\hat{i}} + 2{x}^{2}\mathbf{\hat{j}}\). The magnitude of this electric field can in turn be expressed as a scalar field of the form \(\sqrt{4{x}^{4} + 16{x}^{2}{y}^{2}}\).

Implementation of fields in sympy.physics.vector

In sympy.physics.vector, every ReferenceFrame instance is assigned basis vectors corresponding to the \(X\), \(Y\) and \(Z\) directions. These can be accessed using the attributes named x, y and z respectively. Hence, to define a vector \(\mathbf{v}\) of the form \(3\mathbf{\hat{i}} + 4\mathbf{\hat{j}} + 5\mathbf{\hat{k}}\) with respect to a given frame \(\mathbf{R}\), you would do

>>> from sympy.physics.vector import ReferenceFrame
>>> R = ReferenceFrame('R')
>>> v = 3*R.x + 4*R.y + 5*R.z

Vector math and basic calculus operations with respect to vectors have already been elaborated upon in other sections of this module’s documentation.

On the other hand, base scalars (or coordinate variables) are implemented as special SymPy Symbol s assigned to every frame, one for each direction from \(X\), \(Y\) and \(Z\). For a frame R, the \(X\), \(Y\) and \(Z\) base scalar Symbol s can be accessed using the R[0], R[1] and R[2] expressions respectively.

Therefore, to generate the expression for the aforementioned electric potential field \(2{x}^{2}y\), you would have to do

>>> from sympy.physics.vector import ReferenceFrame
>>> R = ReferenceFrame('R')
>>> electric_potential = 2*R[0]**2*R[1]
>>> electric_potential
2*R_x**2*R_y

In string representation, R_x denotes the \(X\) base scalar assigned to ReferenceFrame R. Essentially, R_x is the string representation of R[0].

Scalar fields can be treated just as any other SymPy expression, for any math/calculus functionality. Hence, to differentiate the above electric potential with respect to \(x\) (i.e. R[0]), you would have to use the diff method.

>>> from sympy.physics.vector import ReferenceFrame
>>> R = ReferenceFrame('R')
>>> electric_potential = 2*R[0]**2*R[1]
>>> from sympy import diff
>>> diff(electric_potential, R[0])
4*R_x*R_y

Like vectors (and vector fields), scalar fields can also be re-expressed in other frames of reference, apart from the one they were defined in – assuming that an orientation relationship exists between the concerned frames. This can be done using the express method, in a way similar to vectors - but with the variables parameter set to True.

>>> from sympy.physics.vector import ReferenceFrame
>>> R = ReferenceFrame('R')
>>> electric_potential = 2*R[0]**2*R[1]
>>> from sympy.physics.vector import dynamicsymbols, express
>>> q = dynamicsymbols('q')
>>> R1 = R.orientnew('R1', rot_type = 'Axis', amounts = [q, R.z])
>>> express(electric_potential, R1, variables=True)
2*(R1_x*sin(q(t)) + R1_y*cos(q(t)))*(R1_x*cos(q(t)) - R1_y*sin(q(t)))**2

Moreover, considering scalars can also be functions of time just as vectors, differentiation with respect to time is also possible. Depending on the Symbol s present in the expression and the frame with respect to which the time differentiation is being done, the output will change/remain the same.

>>> from sympy.physics.vector import ReferenceFrame
>>> R = ReferenceFrame('R')
>>> electric_potential = 2*R[0]**2*R[1]
>>> q = dynamicsymbols('q')
>>> R1 = R.orientnew('R1', rot_type = 'Axis', amounts = [q, R.z])
>>> from sympy.physics.vector import time_derivative
>>> time_derivative(electric_potential, R)
0
>>> time_derivative(electric_potential, R1).simplify()
2*(R1_x*cos(q(t)) - R1_y*sin(q(t)))*(3*R1_x**2*cos(2*q(t))/2 -
R1_x**2/2 - 3*R1_x*R1_y*sin(2*q(t)) - 3*R1_y**2*cos(2*q(t))/2 -
R1_y**2/2)*Derivative(q(t), t)