|
|
Line 95: |
Line 95: |
|
\left(-\frac{1}{r}E^r_{lm}+\frac{\mathrm{d}E^{(1)}_{lm}}{\mathrm{d}r}+\frac{1}{r}E^{(1)}_{lm}\right)\mathbf{\Phi}_{lm}\right)</math> |
|
\left(-\frac{1}{r}E^r_{lm}+\frac{\mathrm{d}E^{(1)}_{lm}}{\mathrm{d}r}+\frac{1}{r}E^{(1)}_{lm}\right)\mathbf{\Phi}_{lm}\right)</math> |
|
|
|
|
|
|
==Examples== |
|
==Example: application to electrodynamics== |
|
===Application to electrodynamics== |
|
The VHS are especially useful in the study of multipole radiation fields. For instance, a magnetic multipole is due to an oscillating current with angular frequency <math>\omega\,</math> and complex amplitude |
|
The VHS are especially useful in the study of multipole radiation fields. For instance, a magnetic multipole is due to an oscillating current with angular frequency <math>\omega\,</math> and complex amplitude |
|
|
|
|
Line 126: |
Line 127: |
|
In this way, the partial differential equations have been transformed in a set of ordinary differential equations. |
|
In this way, the partial differential equations have been transformed in a set of ordinary differential equations. |
|
|
|
|
|
==Example: application to fluid dynamics== |
|
===Application to fluid dynamics=== |
|
|
|
|
|
In the calculation of the [[Stokes' law]] for the drag than a viscous fluid exerts on a small spherical particle, the velocity distribution obeys [[Navier-Stokes equations]] neglecting inertia, i.e. |
|
In the calculation of the [[Stokes' law]] for the drag than a viscous fluid exerts on a small spherical particle, the velocity distribution obeys [[Navier-Stokes equations]] neglecting inertia, i.e. |
Definition
The Vector Spherical Harmonics (VSH) are an extension of the scalar Spherical harmonics for the use with vector fields.
Several conventions have been used to define the VSH[1][2][3][4]. We follow that of Barrera et al. Given a scalar spherical harmonic
we define three VSH:
![{\displaystyle \mathbf {Y} _{lm}=Y_{lm}{\hat {\mathbf {r} }}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2939e0959a607f7c2b51fbc83aba73a2f7cb3cb0)
![{\displaystyle \mathbf {\Psi } _{lm}=r\nabla Y_{lm}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e5d0d4637bd38b7a4ddb397c2f8d43d6b133c655)
![{\displaystyle \mathbf {\Phi } _{lm}=\mathbf {r} \times \nabla Y_{lm}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/776814b9163f018a6c15aef9fbaf470147a142c3)
being
the unitary vector along the radial direction. The radial factors are included to guarantee that the dimensions of the VSH are the same as the ordinary spherical harmonics and that the VSH do not depend on the radial spherical coordinate.
The interest of these new vector fields is to separate the radial dependence from the angular one when using spherical ordinates, so that a vector field admits a multipole expansion
![{\displaystyle \mathbf {E} =\sum _{l=0}^{\infty }\sum _{m=-l}^{l}\left(E_{lm}^{r}(r)\mathbf {Y} _{lm}+E_{lm}^{(1)}(r)\mathbf {\Psi } _{lm}+E_{lm}^{(2)}(r)\mathbf {\Phi } _{lm}\right)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c9833d256febd27550f3ba31dc3088448577d377)
The labels on the components reflect that
is the radial component of the vector field, while
and
are transverse components.
Main Properties
Orthogonality
The VSH are orthogonal in the usual three-dimensional way
![{\displaystyle \mathbf {Y} _{lm}\cdot \mathbf {\Psi } _{lm}=0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/93c70f72fe4ee4d10495e9b1742fc55b98c8e917)
![{\displaystyle \mathbf {Y} _{lm}\cdot \mathbf {\Phi } _{lm}=0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/17371ac0480a34337c092f1d51808da0d316da9a)
![{\displaystyle \mathbf {\Psi } _{lm}\cdot \mathbf {\Phi } _{lm}=0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/aca1cce442a623a1771c2fdce39dd9eb16cdc7ea)
but also in the Hilbert space
![{\displaystyle \int \mathbf {Y} _{lm}\cdot \mathbf {Y} _{l'm'}^{*}\,\mathrm {d} \Omega =\delta _{ll'}\delta _{mm'}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/54640c27918c5decf186257837567a09dd4f9220)
![{\displaystyle \int \mathbf {\Psi } _{lm}\cdot \mathbf {\Psi } _{l'm'}^{*}\,\mathrm {d} \Omega =l(l+1)\delta _{ll'}\delta _{mm'}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4b727e6bea71db9eed320a374902fee8242266ca)
![{\displaystyle \int \mathbf {\Phi } _{lm}\cdot \mathbf {\Phi } _{l'm'}^{*}\,\mathrm {d} \Omega =l(l+1)\delta _{ll'}\delta _{mm'}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/604afe65de48e21abca4891b1decbb22e1ad5553)
![{\displaystyle \int \mathbf {Y} _{lm}\cdot \mathbf {\Psi } _{l'm'}^{*}\,\mathrm {d} \Omega =0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d602f1d669439752f0ece21ff1acfcbcbd80a21a)
![{\displaystyle \int \mathbf {Y} _{lm}\cdot \mathbf {\Phi } _{l'm'}^{*}\,\mathrm {d} \Omega =0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/813bcdee5a3a112d1b86588b52046f4b4748c946)
![{\displaystyle \int \mathbf {\Psi } _{lm}\cdot \mathbf {\Phi } _{l'm'}^{*}\,\mathrm {d} \Omega =0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/262967190d5d0f0d7115f87ef23313d306be7592)
Vector multipole moments
The orthogonality relatios allow to compute the spherical multipole moments of a vector field as
![{\displaystyle E_{lm}^{r}=\int \mathbf {E} \cdot \mathbf {Y} _{lm}^{*}\,\mathrm {d} \Omega }](https://wikimedia.org/api/rest_v1/media/math/render/svg/7e2f1d3786c2ed63106d22efd02e340d3daa5645)
![{\displaystyle E_{lm}^{(1)}={\frac {1}{l(l+1)}}\int \mathbf {E} \cdot \mathbf {\Psi } _{lm}^{*}\,\mathrm {d} \Omega }](https://wikimedia.org/api/rest_v1/media/math/render/svg/dbcfe86d54b98008748fbbc05b1962ece0055570)
![{\displaystyle E_{lm}^{(2)}={\frac {1}{l(l+1)}}\int \mathbf {E} \cdot \mathbf {\Phi } _{lm}^{*}\,\mathrm {d} \Omega }](https://wikimedia.org/api/rest_v1/media/math/render/svg/55755d8053054cb2114f8f393a923714fc7cfb12)
The gradient of a scalar field
Given the multipole expansion of a scalar field
![{\displaystyle \phi =\sum _{l=0}^{\infty }\sum _{m=-l}^{l}\phi _{lm}(r)Y_{lm}(\theta ,\phi )}](https://wikimedia.org/api/rest_v1/media/math/render/svg/475ee6c1a1ea8dbf6203ee07df3ae0c4521f214f)
we can express its gradient in terms of the VSH as
![{\displaystyle \nabla \phi =\sum _{l=0}^{\infty }\sum _{m=-l}^{l}\left({\frac {\mathrm {d} \phi _{lm}}{\mathrm {d} r}}\mathbf {Y} _{lm}+{\frac {\phi _{lm}}{r}}\mathbf {\Psi } _{lm}\right)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/58cbacb73c1f2429fb93f6f6b38142c8391b0c4b)
Divergence
For any multipole field we have
![{\displaystyle \nabla \cdot \left(f(r)\mathbf {Y} _{lm}\right)=\left({\frac {\mathrm {d} f}{\mathrm {d} r}}+{\frac {2}{r}}f\right)Y_{lm}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ea375aea2717a694d0793e4ab8b863d1b708f2ae)
![{\displaystyle \nabla \cdot \left(f(r)\mathbf {\Psi } _{lm}\right)=-{\frac {l(l+1)}{r}}fY_{lm}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/053504d374ee7d0054f0306f3f300022748b14ac)
![{\displaystyle \nabla \cdot \left(f(r)\mathbf {\Phi } _{lm}\right)=0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/61494f370773f8887f131f2d62b22dced17d6eaf)
By superposition we obtain the divergence of any vector field
![{\displaystyle \nabla \cdot \mathbf {E} =\sum _{l=0}^{\infty }\sum _{m=-l}^{l}\left({\frac {\mathrm {d} E_{lm}^{r}}{\mathrm {d} r}}+{\frac {2}{r}}E_{lm}^{r}-{\frac {l(l+1)}{r}}E_{lm}^{(1)}\right)Y_{lm}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/55d4de2254b3a875a5a076c7fa777e43956f6593)
we see that the component on
is always solenoidal.
Curl
For any multipole field we have
![{\displaystyle \nabla \times \left(f(r)\mathbf {Y} _{lm}\right)=-{\frac {1}{r}}f\mathbf {\Phi } _{lm}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/020c1b4f2e05a38baf76d20676f2fe3820ebd302)
![{\displaystyle \nabla \times \left(f(r)\mathbf {\Psi } _{lm}\right)=\left({\frac {\mathrm {d} f}{\mathrm {d} r}}+{\frac {1}{r}}f\right)\mathbf {\Phi } _{lm}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d24b8eb273205d2ce539d3739ce19e206fcf087f)
![{\displaystyle \nabla \times \left(f(r)\mathbf {\Phi } _{lm}\right)=-{\frac {l(l+1)}{r}}f\mathbf {Y} _{lm}-\left({\frac {\mathrm {d} f}{\mathrm {d} r}}+{\frac {1}{r}}f\right)\mathbf {\Psi } _{lm}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/56d01cdea4338bf785a1508712850486c423a68a)
By superposition we obtain the curl of any vector field
![{\displaystyle \nabla \times \mathbf {E} =\sum _{l=0}^{\infty }\sum _{m=-l}^{l}\left(-{\frac {l(l+1)}{r}}E_{lm}^{(2)}\mathbf {Y} _{lm}-\left({\frac {\mathrm {d} E_{lm}^{(2)}}{\mathrm {d} r}}+{\frac {1}{r}}E_{lm}^{(2)}\right)\mathbf {\Psi } _{lm}+\left(-{\frac {1}{r}}E_{lm}^{r}+{\frac {\mathrm {d} E_{lm}^{(1)}}{\mathrm {d} r}}+{\frac {1}{r}}E_{lm}^{(1)}\right)\mathbf {\Phi } _{lm}\right)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d257f8e8578b9738593f3b1f614470f107970e35)
Examples
=Application to electrodynamics
The VHS are especially useful in the study of multipole radiation fields. For instance, a magnetic multipole is due to an oscillating current with angular frequency
and complex amplitude
![{\displaystyle {\hat {\mathbf {J} }}=J(r)\mathbf {\Phi } _{lm}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f2f9ac8d6503f53533424004b3a0d7f08c0ff1a8)
and the corresponding electric and magnetic fields can be written as
![{\displaystyle {\hat {\mathbf {E} }}=E(r)\mathbf {\Phi } _{lm}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ec5c92c42c1ab8ad0940fd35231f5be3f7fabd6f)
![{\displaystyle {\hat {\mathbf {B} }}=B^{r}(r)\mathbf {Y} _{lm}+B^{(1)}(r)\mathbf {\Psi } _{lm}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3a73b4b1e6824af3907bd3a1319f3e7c87206f98)
Substituting into Maxwell equations, Gauss' law is automathically satisfied
![{\displaystyle \nabla \cdot {\hat {\mathbf {E} }}=0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f586d4096e6abb463fa92e624fe5a5ff23a51399)
while Faraday's law decouples in
![{\displaystyle \nabla \times {\hat {\mathbf {E} }}=-\mathrm {j} \omega {\hat {\mathbf {B} }}\qquad \Rightarrow \qquad \left\{{\begin{array}{l}\displaystyle {\frac {l(l+1)}{r}}E=\mathrm {j} \omega B^{r}\\\ \\\displaystyle {\frac {\mathrm {d} E}{\mathrm {d} r}}+{\frac {E}{r}}=\mathrm {j} \omega B^{(1)}\end{array}}\right.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/6a63158500633f50ed7bc48e7eab402eafee7b65)
Gauss' law for the magnetic field implies
![{\displaystyle \nabla \cdot {\hat {\mathbf {B} }}=0\quad \Rightarrow \quad {\frac {\mathrm {d} B^{r}}{r}}+{\frac {2}{r}}B^{r}-{\frac {l(l+1)}{r}}B^{(1)}=0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5d035a21337200bd2fb12b00b7b4cfba1e0b8e83)
and Ampère-Maxwell's equation gives
![{\displaystyle \nabla \times {\hat {\mathbf {B} }}=\mu _{0}{\hat {\mathbf {J} }}+\mathrm {j} \mu _{0}\varepsilon _{0}\omega {\hat {\mathbf {E} }}\quad \Rightarrow \quad -{\frac {B^{r}}{r}}+{\frac {\mathrm {d} B^{(1)}}{\mathrm {d} r}}+{\frac {B^{(1)}}{r}}=\mu _{0}J+\mathrm {j} \omega \mu _{0}\varepsilon _{0}E}](https://wikimedia.org/api/rest_v1/media/math/render/svg/adb36bdd8012d1989e80e8f261160c35173b1444)
In this way, the partial differential equations have been transformed in a set of ordinary differential equations.
Application to fluid dynamics
In the calculation of the Stokes' law for the drag than a viscous fluid exerts on a small spherical particle, the velocity distribution obeys Navier-Stokes equations neglecting inertia, i.e.
![{\displaystyle \nabla \cdot \mathbf {v} =0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/760025de0202db29f86a1c2c4bef18d54ccc0f64)
![{\displaystyle \mathbf {0} =-\nabla p+\eta \nabla ^{2}\mathbf {v} }](https://wikimedia.org/api/rest_v1/media/math/render/svg/6542457581a611661422f280a8a2ae3a70c77ade)
with the boundary conditions
![{\displaystyle \mathbf {v} =\mathbf {0} \quad (r=a)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5b31a142330b4b9ea874468dfc0030961d70654e)
![{\displaystyle \mathbf {v} =-\mathbf {U} _{0}\quad (r\to \infty )}](https://wikimedia.org/api/rest_v1/media/math/render/svg/27418386f995be5e17b153fe4f973f78f7e17556)
being
the relative velocity of the particle to the fluid far from the particle. In spherical coordinates this velocity at infinity can be written as
![{\displaystyle \mathbf {U} _{0}=U_{0}\left(\cos \theta \,{\hat {\mathbf {r} }}-\sin \theta \,{\hat {\mathbf {\theta } }}\right)=U_{0}\left(\mathbf {Y} _{10}+\mathbf {\Psi } _{10}\right)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ff82cf5bad7c3741a0493360d05ec26b71cffb89)
The last expression suggest a expansion on spherical harmonics for the liquid velocity and the pressure
![{\displaystyle p=P(r)Y_{10}\,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/44e66cdcf7151118ba616465a76ff07301956ec0)
![{\displaystyle \mathbf {v} =v^{r}(r)\mathbf {Y} _{10}+v^{(1)}(r)\mathbf {\Psi } _{10}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e71cee172c94ab9de87ceeb7939501a2c9b6ee7b)
Substitution of Stokes equation produces a set of ordinary differential equations for the coefficients.
See also
External links
[1] Vector Spherical Harmonics at Eric Weisstein's Mathworld
References
- ^
R.G. Barrera, G.A. Estévez and J. Giraldo, Vector spherical harmonics and their application to magnetostatic, Eur. J. Phys. 6 287-294 (1985)
- ^ B. Carrascal, G.A. Estevez, P. Lee and V. Lorenzo Vector spherical harmonics and their application to classical electrodynamics, Eur. J. Phys., 12, 184-191 (1991)
- ^
E. L. Hill, The theory of Vector Spherical Harmonics, Am. J. Phys. 22, 211-214 (1954)
- ^
E. J. Weinberg, Monopole vector spherical harmonics, Phys. Rev. D. 49, 1086-1092 (1994)