7. Tensor Calculus in Euclidean Spaces#
This section deals with tensor calculus in Euclidean space or on manifolds embedded in Euclidean spaces, focusing on \(d\)-dimensional spaces with \(d \le 3\), with inner product.
This section may rely on results of differential geometry.
7.1. Coordinates#
A set of parameters \(\{q^a\}_{a=1:d}\) to represent vector (or point) in space,
if \(\vec{r} \in E^{d}\), \(a=1:d\).
In \(E^3\),
Coordinate lines, 2-parameter family of lines, keeping 2 coordinates constant. As an example, coordinate lines with constant \(q^2, \, q^3\)
\[\vec{r}_1(q^1) = \vec{r}(q^1, \bar{q}^2, \bar{q}^3) \ .\]Coordinate surfaces, 1-parameter family of surfaces, keeping 1 coordinate constant. As an example, coordinate surfaces with constant \(q^1\),
\[\vec{r}_{23}(q^2, q^3) = \vec{r}(\bar{q}^1, q^2, q^3) \ .\]
Definition 7.1 (Regular parametrization)
If \(\frac{\partial \vec{r}}{\partial q^a} \ne 0\).
7.1.1. Natural basis#
Definition 7.2 (Natural basis)
Vectors of natural basis
Definition 7.3 (Reciprocal basis (todo move to Tensor Algebra))
Given a basis \(\{ \vec{b}_a \}_{a}\), its reciprocal basis the set of vector \(\{ \vec{b}^b \}_b\) defined as
being \(\delta_a^b\) Kronecker delta.
Definition 7.4 (Christoffel symbols)
Christoffel symbols (of the \(2^{nd}\) kind) are defined as the components of the partial derivatives of the vectors of a natural basis w.r.t. the coordinates referred to the natural basis itself,
Properties of Christoffel symbols
Exploiting the definition of reciprocal basis, Christoffel symbols can be written as
Symmetry. Symmetry os the lower indices
readily follows Schwartz theorem about partial derivatives
7.2. Fields#
Function of the points in space \(F: E^d \rightarrow V^r\), being \(V^r\) a space of tensors of order \(r\).
7.3. Differential operators#
7.3.1. Directional derivative#
The directional derivative for an arbitrary vector \(\vec{v} \in V\)
is evaluated as the limit for \(\alpha \rightarrow 0\) of the incremental ratio
7.3.2. Gradient#
The gradient is the differential operator is the first-order differential operator appearing in the definition of the directional derivative, \(\nabla F(\vec{r})\). It takes a tensor field \(F(\vec{r})\) of order \(r\) and gives a tensor field \(\nabla F(\vec{r})\) of order \(r+1\). Given a set of coordinates \(\{q^a\}_{a=1:d}\), the gradient can be written using the reciprocal basis of the natural basis as
Examples. …
Example 7.1 (Gradient of a scalar field - with general coordinates \(q^{a}\))
Applying the definition (7.2) of gradient operator, it readily follows
Example 7.2 (Gradient of a vector field - with general coordinates \(q^{a}\))
Applying the definition (7.2) of gradient operator, rule for the derivative of a product and the definition (7.1) of Christoffel symbols to write derivatives of base vectors,
Example 7.3 (Gradient of a \(2^{nd}\)-order tensor field - with general coordinates \(q^{a}\))
Applying the definition (7.2) of gradient operator, rule for the derivative of a product and the definition (7.1) of Christoffel symbols to write derivatives of base vectors,
7.3.3. Divergence#
Divergence opearator is a first-order differential operator that can be defined as the contraction of the first two indices of the gradient,
It takes a tensor field \(F(\vec{r})\) of order \(r \ge 1\) and gives a tensor field \(\nabla \cdot F(\vec{r})\) of order \(r-1 \ge 0\).
Example 7.4 (Divergence of a vector field - with general coordiantes \(q^{a}\))
Applying contraction to the gradient of a vector field, it readily follows,
Example 7.5 (Divergence of a \(2^{nd}\)-order tensor field - with general coordiantes \(q^{a}\))
Applying contraction to the gradient of a vector field, it readily follows,
7.3.4. Laplacian#
Laplacian operator is second-order differential operator that can be defined as the divergence of the gradient,
Example 7.6 (Laplacian of a scalar field - with general coordinates \(q^a\))
Example 7.7 (Laplacian of a vector field - with general coordinates \(q^a\))
7.3.5. Curl#
7.4. Integrals in \(E^d\), \(d \le 3\)#
7.4.1. Line integrals#
7.4.1.1. Density#
Integrals
represent the summation of contributions \(F(\vec{r})\) over elementary segments of path \(\gamma\), whose dimension is \(|d \vec{r}|\), i.e. implicitly means
Given a regular parametrization of the curve \(\vec{r}(q^1)\) (with increasing \(q^1\) so that \(|dq^1| = dq^1\)), and the differential \(d \vec{r} = \vec{r}'(q^1) \, d q^1\), the integral can be written as an integral in the parameter \(q^1\)
with \(\vec{r}(q^1_a)\), \(\vec{r}(q^1_b)\) the extreme points of path \(\gamma\).
7.4.1.2. Work#
Integrals
implicitly mean
as \(\hat{t} = \frac{d \vec{r}}{|d \vec{r}|}\). Given a regular parametrization of the curve \(\vec{r}(q^1)\) (with increasing \(q^1\) so that \(|dq^1| = dq^1\)), and the differential \(d \vec{r} = \vec{r}'(q^1) \, d q^1\), the integral can be written as an integral in the parameter \(q^1\)
7.4.2. Surface integrals#
Given two coordinates \(q^1, \, q^2\) describing a surface, \(\vec{r}(q^1, q^2)\) the elementary surface with unit normal reads
and the elementary surface thus reads
7.4.2.1. Density#
Integrals
implicitly mean
Given regular parametrization of the surface, \(\vec{r}(q^1, \, q^2), \ (q^1, q^2) \in Q^{12}\), the integral can be written as the multi-dimensional integral in coordinates \(q^1, \ q^2\),
7.4.2.2. Flux#
Integrals
implicitly mean
Given regular parametrization of the surface, \(\vec{r}(q^1, \, q^2), \ (q^1, q^2) \in Q^{12}\), the integral can be written as the multi-dimensional integral in coordinates \(q^1, \ q^2\),
7.4.3. Volume#
7.4.3.1. Density#
Integrals
implicitly mean
Given regular parametrization of the volume, \(\vec{r}(q^1, \, q^2, \, q^3), \ (q^1, q^2, q^3) \in Q\), the integral can be written as the multi-dimensional integral in coordinates \(q^1, \, q^2, \, q^3\),
7.4.4. Theorems#
7.4.4.1. Two useful lemmas#
The next lemma forms the foundation of the well-known divergence theorem and gradient theorem:
the proof of these two theorems is based on a straightforward repeated use of this lemma.
Given how simple this lemma is and how frequently it is applied in writing balances and, more generally, in integration by parts, it is very convenient to remember this simple result.
Theorem 7.1 (Lemma 1.)
Under the assumptions of Green’s lemma in the plane,
Proof
The reasoning closely follows the one used for the proof of Green’s lemma in the plane.
For \({\partial A}/{\partial z}\):
The most complex step is transitioning from the integral over \((x,y) \in R\) to the surface integral over \(S\), the boundary of volume \(V\): the infinitesimal area element \(dR\) in the xy-plane is equal to \(dR = dx dy\); the drawing and the proof refer to a simple volume (as in the case of Green’s lemma in the plane, the results can be generalized to domains of arbitrary shape).
It is possible to divide the surface \(S\) into two “halves” \(S^+: z = f_2(x,y)\) and \(S^-: z = f_1(x,y)\) such that
\(S^+ \cup S^- = S\), and the outward normal has positive and negative z-components respectively (\(S^+: \mathbf{\hat{n}}\cdot\mathbf{\hat{z}}>0\), \(S^-: \mathbf{\hat{n}} \cdot\mathbf{\hat{z}}<0\)).
The surface element \(dR\) is also the projection of the surface element \(dS\) onto the xy-plane: in general, \(dS\) will not be parallel to the xy-plane and thus will be larger than \(dR\). It’s not difficult to show that:
We can now continue the proof:
Just as the previous lemma forms the basis for the proof of the gradient theorem and the divergence theorem,
the following lemma forms the basis for the proof of the curl theorem.
Theorem 7.2 (Lemma 2.)
Under the assumptions of Green’s lemma in the plane,
Proof
For \(A\mathbf{\hat{e}_x}\), we have \({\nabla} \times (A \mathbf{\hat{e}_x}) = {\partial A}/{\partial z} \mathbf{\hat{e}_y} - {\partial A}/{\partial y} \mathbf{\hat{e}_z}\). The surface \(S\) is written in parametric form as: \(\mathbf{r} = x\mathbf{\hat{e}_x} + y\mathbf{\hat{e}_y} + z(x,y)\mathbf{\hat{e}_z}\). The vector \({\partial \mathbf{r}}/{\partial y} = \mathbf{\hat{e}_y} + {\partial z}/{\partial y} \mathbf{\hat{e}_z} \) is tangent to the surface \(S\) and hence perpendicular to the normal \(\mathbf{\hat{n}}\):
Now writing \([\mathbf{\nabla} \times (A \mathbf{\hat{e}_x})] \cdot \mathbf{\hat{n}}\):
Recognizing that \(\frac{\partial A(x,y,z(x,y))}{\partial y} = \frac{\partial A}{\partial z} \frac{\partial z}{\partial y} + \frac{\partial A}{\partial y}\), we can write:
7.4.4.2. Gradient theorem#
…todo assumptions…
Proof for simple domains \(\ V\)
This result immediately follows from Lemma 1 Theorem 7.1
having (1) exploited the freedom to put unit vectors of the Cartesian basis inside the integrals, as they’re unifrom - constant in space -, and (2) recognized the expression of the gradient of a scalar field expressed using Cartesian coordinates, as shown in Example 7.8.
7.4.4.3. Divergence theorem#
…todo assumptions…
Proof
This result immediately follows from Lemma 1 Theorem 7.1
having (1) exploited the freedom to put unit vectors of the Cartesian basis inside the integrals, as they’re unifrom - constant in space -, and (2) recognized the expression of the divergence of a vector field expressed using Cartesian coordinates, as shown in Example 7.11.
7.4.4.4. Curl theorem#
…todo assumptions…
Proof
This proof seamlessly follows from Lemma Theorem 7.2, applied to all the Cartesian contributions of the vector field
as