[Math] Calculate surface integral of point charge located outside the surface

calculusmultivariable-calculus

Problem

Given:

$$\vec r = r(\sin \theta \cos \phi, \sin \theta \sin \phi, \cos \theta)$$

$S: | \vec r | = a$, with $\hat n$ outwards

$$\vec r_0 = \frac{3a}{5}(\hat x + \hat y – \hat z) = \frac{3a}{5}(1,1,-1)$$

$$\vec F = k \frac{\vec r – \vec r_0}{ {| \vec r – \vec r_0 |}^3 }$$

Calculate:

$$\int_S \vec F \cdot \mathrm d \vec S $$

Solution (added after accepted answer)

Thanks to the accepted answer which confirmed that the explicit calculation on/in $S$ is messy, and the suggestion of using translation invariance, I have noted down the solution below (let me know if you have further suggestions).

An argument is used which concludes that the integrand is zero inside the sphere $S$.

Using Gauss' theorem: $\int_S \vec F \cdot \mathrm d \vec S = \int_V \nabla \cdot \vec F \mathrm dV $.

Since the calculation is not easy to explicitly calculate in (or on) $S$, instead it is shown that $\nabla \cdot \vec F = 0$ exterior to $S_{\epsilon}$, which is a sphere containing the singularity, and since $S$ is exterior to $S_{\epsilon}$ the integrand must be $0$ and so $\int_V \nabla \cdot \vec F \mathrm dV = 0$.

$$S_{\epsilon}: |\vec r – \vec r_0 | = \epsilon, \epsilon > 0$$

translate the coordinate system so that $S_{\epsilon}$ is the origin in the translated system:

$$\vec R = \vec r – \vec r_0 = R (\sin \theta \cos \phi, \sin \theta \sin \phi, \cos \theta)$$

$$F(\vec R) = k \frac{\vec R}{|\vec R|^3} = k \frac{1}{R^2} \hat R$$

$$S_{\epsilon}: |\vec R| = \epsilon$$

$$\nabla \cdot \vec F (\vec R) = \frac{1}{R^2} \frac{\partial}{\partial R} \left(R^2 \vec F_R \right) = \frac{1}{R^2} \frac{\partial}{\partial R} \left(R^2 k \frac{1}{R^2} \right) = 0$$, ($k$ is a constant)

This is true for any value of $R$ except at the singularity ($\epsilon$ can be made as small as required), since $S$ does not contain the singularity, the integrand ($\nabla \cdot \vec F(\vec r)$ is 0 and the integral is therefore zero.

My question (updated)

I am having difficulties in explicitly calculating the value of this integral.
Specifically, the divergence of the field becomes messy, I am not able to see how I can use the symmetry of $S$ due to $\vec r_0$. (When applying the Gauss' theorem, I am stuck in evaluating the divergence of the field).

I am able to argue for that this integral is indeed zero (the field only has one singularity and it is exterior to the sphere $S$). However, I am not able to explicitly show (by calculation) that this integral is zero. Thankful for any help.

Please note that this is not homework, I am studying for an exam.

My question (original)

What is an easy way to calculate this integral? Any suggestions on the approaches below?

My apologies if this is due to lack of some basic knowledge (I am back studying after 2.5 years)

Intuitively I understand the integral is zero (the point charge is located outside the sphere, anything flowing into the sphere will also flow out), however, I have issues with the calculation.

Calculate directly

$$\int_S \vec F \cdot \mathrm d \vec S = \int_0^{\pi} \mathrm d \theta \int_0^{2\pi} \mathrm d \phi r^2 \sin {\theta} \vec F_r \cdot \hat r = \int_0^{\pi} \mathrm d \theta \int_0^{2\pi} \mathrm d \phi r^2 \sin \theta k \frac{a – \frac{3a}{5} \left( \sin \theta \cos \phi + \sin \theta \sin \phi – \cos \theta) \right)}{ {\left| \vec r – \vec r_0 \right|}^3 }$$

Using Gauss theorem

I am tempted to use Gauss theorem, $\int_S \vec F \cdot \mathrm d \vec S = \int_V \nabla \cdot \vec F \mathrm dV $, in spherical coordinates. However $\vec F$ has components also in $\hat \theta$ and $\hat \phi$ (due to $\vec r_0$) and ${\left| \vec r – \vec r_0 \right|}^3$ is not that nice to derivate.

$$\vec F_r = k \frac{\vec r – \vec r_0}{ {\left| \vec r – \vec r_0 \right|}^3 } \cdot \hat r = k \frac{r – \frac{3a}{5} \left( \sin \theta \cos \phi + \sin \theta \sin \phi – \cos \theta) \right)}{ {\left| \vec r – \vec r_0 \right|}^3 }$$

$$\vec F_{\theta} = k \frac{\vec r – \vec r_0}{ {\left| \vec r – \vec r_0 \right|}^3 } \cdot \hat \theta = k \frac{ \frac{3a}{5} \left( \cos \theta \cos \phi + \cos \theta \sin \phi – \sin \theta) \right)}{ {\left| \vec r – \vec r_0 \right|}^3 }$$

$$\vec F_{\phi} = k \frac{\vec r – \vec r_0}{ {\left| \vec r – \vec r_0 \right|}^3 } \cdot \hat \phi = k \frac{ \frac{3a}{5} \left( – \sin \phi + \cos \phi ) \right)}{ {\left| \vec r – \vec r_0 \right|}^3 }$$

$$\nabla \cdot \vec F = \frac{1}{r^2} \frac{\partial}{\partial r} \left(r^2 \vec F_r \right) + \frac{1}{r \sin \theta} \frac{\partial}{\partial \theta} \left( \sin \theta F_{\theta} \right) + \frac{1}{r \sin \theta} \frac{\partial F_{\phi}}{\partial \phi}$$

Best Answer

(I don't know what you mean by $\hat x+\hat y -\hat z$, so I disregard the information about $r_0$.) Your field $\vec F$ is the gravitational field produced by a point mass at $r_0$. The divergence of this field is $\equiv0$ away from $r_0$. If the point $r_0$ is in the exterior of the sphere $S$ then it follows by Gauss' theorem that the integral in question is $0$. If $r_0$ lies in the interior of $S$ then you should draw a tiny sphere $S_\epsilon$ with center $r_0$ and apply Gauss' theorem to the region between $S_\epsilon$ and $S$. You are left with an integral over $S_\epsilon$ which you can almost do in your head.

Related Question