Past exam of the mathematics course of the University of Cambridge / 2009 / ia / Paper 3 / 4B (Vector Calculus)

Table of contents
- Solution 4B

Solution

For Cartesian components,

\partial_{j} x_{i} = δ_{ij}

, the Kronecker delta. Differentiating

r^{2} = x_{i} x_{i}

gives

\partial_{j} r = \frac{x _{j}}{r} (r > 0) .

Put

s = k_{ℓ} x_{ℓ}

. Using the product rule and

\partial_{j} r^{- q} = - q x_{j} r^{- q - 2}

gives the velocity gradient

d_{ij} = - \frac{k _{i} x _{j}}{r ^{3}} + \frac{k _{j} x _{i}}{r ^{3}} + \frac{s δ _{ij}}{r ^{3}} - \frac{3 s x _{i} x _{j}}{r ^{5}} .

This is the rate of strain and vorticity of a Stokeslet, apart from the overall physical normalization. Its symmetric and antisymmetric parts of a matrix are

E_{ij} = \frac{d _{ij} + d _{ji}}{2} = \frac{s δ _{ij}}{r ^{3}} - \frac{3 s x _{i} x _{j}}{r ^{5}}, W_{ij} = \frac{d _{ij} - d _{ji}}{2} = \frac{k _{j} x _{i} - k _{i} x _{j}}{r ^{3}} .

Here

E_{ji} = E_{ij}

and

W_{ji} = - W_{ij}

, and

d = E + W

. Taking the trace gives the divergence

\nabla \cdot u = d_{ii} = \frac{3 s}{r ^{3}} - \frac{3 s r ^{2}}{r ^{5}} = 0.

For the curl, use the Levi-Civita symbol. Because

d_{kj} = \partial_{j} u_{k}

, the symmetric part makes no contribution, and

(\nabla \times u)_{i} = ϵ_{ijk} d_{kj} = \frac{ϵ _{ijk} ( k _{j} x _{k} - k _{k} x _{j} )}{r ^{3}} = \frac{2 ϵ _{ijk} k _{j} x _{k}}{r ^{3}} .

Thus

\nabla \times u = 2 k \times x / r^{3}

. All derivative identities hold away from the singular origin.

 Ancestors (9)

 View article source

 Discussion (0) New discussion

There are no discussions about this article yet.

 Articles by others on the same topic (0)

There are currently no matching articles.

  See all articles in the same topic Create my own version