Bipolar cylindrical coordinates: Difference between revisions

From formulasearchengine
Jump to navigation Jump to search
en>Addbot
m Bot: Migrating 1 interwiki links, now provided by Wikidata on d:q4915481
en>Jonesey95
m Fixing citation DOI or LCCN parameter error using AutoEd
 
Line 1: Line 1:
[[File:Oblate spheroidal coordinates half hyperboloid.png|thumb|right|350px|Figure 1: Coordinate isosurfaces for a point ''P'' (shown as a black sphere) in oblate spheroidal coordinates (''μ'', ''ν'', ''φ'').  The ''z''-axis is vertical, and the foci are at ±2.  The red oblate spheroid (flattened sphere) corresponds to ''μ'' = 1, whereas the blue half-hyperboloid corresponds to ''ν'' = 45°.  The azimuth ''φ'' = −60° measures the [[dihedral angle]] between the green ''x''-''z'' half-plane and the yellow half-plane that includes the point '''P'''.  The [[Cartesian coordinate system|Cartesian coordinates]] of '''P''' are roughly (1.09, −1.89, 1.66).]]
The writer is called Irwin Wunder but it's not the most masucline name out there. One of the things she enjoys most is to do aerobics and now she is trying to earn cash with it. Hiring has been my profession for some time but I've already applied for another 1. North Dakota is where me and my husband live.<br><br>Also visit my webpage [http://rivoli.enaiponline.com/user/view.php?id=438251&course=1 rivoli.enaiponline.com]
 
'''Oblate spheroidal coordinates''' are a three-dimensional [[orthogonal coordinates|orthogonal]] [[coordinate system]] that results from rotating the two-dimensional [[elliptic coordinates|elliptic coordinate system]] about the non-focal axis of the ellipse, i.e., the symmetry axis that separates the foci. Thus, the two foci are transformed into a ring of radius <math>a</math> in the ''x''-''y'' plane.  (Rotation about the other axis produces [[prolate spheroidal coordinates]].)  Oblate spheroidal coordinates can also be considered as a [[limiting case]] of [[ellipsoidal coordinates]] in which the two largest [[semi-axis|semi-axes]] are equal in length.
 
Oblate spheroidal coordinates are often useful in solving [[partial differential equation]]s when the boundary conditions are defined on an [[oblate spheroid]] or a [[hyperboloid|hyperboloid of revolution]]. For example, they played an important role in the calculation of the [[Perrin friction factors]], which contributed to the awarding of the 1926 [[Nobel Prize in Physics]] to [[Jean Baptiste Perrin]].  These friction factors determine the [[rotational diffusion]] of molecules, which affects the feasibility of many techniques such as [[protein NMR]] and from which the hydrodynamic volume and shape of molecules can be inferred.  Oblate spheroidal coordinates are also useful in problems of electromagnetism (e.g., dielectric constant of charged oblate molecules), acoustics (e.g., scattering of sound through a circular hole), fluid dynamics  (e.g., the flow of water through a firehose nozzle) and the diffusion of materials and heat (e.g., cooling of a red-hot coin in a water bath)
 
==Definition (μ, ν, φ)==
[[File:OblateSpheroidCoord.png|thumb|350 px|right|Figure 2: Plot of the oblate spheroidal coordinates &mu; and &nu; in the ''x''-''z'' plane, where φ is zero and ''a'' equals one. The curves of constant ''&mu;'' form red ellipses, whereas those of constant ''&nu;'' form cyan half-hyperbolae in this plane. The ''z''-axis runs vertically and separates the foci; the coordinates ''z'' and ν always have the same sign. The surfaces of constant &mu; and &nu; in three dimensions are obtained by rotation about the ''z''-axis, and are the red and blue surfaces, respectively, in Figure 1.]]
The most common definition of oblate spheroidal coordinates (μ, ν, φ) is
 
:<math>
x = a \ \cosh \mu \ \cos \nu \ \cos \phi
</math>
 
:<math>
y = a \ \cosh \mu \ \cos \nu \ \sin \phi
</math>
 
:<math>
z = a \ \sinh \mu \ \sin \nu
</math>
 
where μ is a nonnegative real number and the angle ν lies between ±90°.  The azimuthal angle φ can fall anywhere on a full circle, between ±180°.  These coordinates are favored over the alternatives below because they are not degenerate; the set of coordinates (μ, ν, φ) describes a unique point in Cartesian coordinates (''x'', ''y'', ''z''). The reverse is also true, except on the ''z''-axis and the disk in the ''x-y'' plane inside the focal ring.
 
===Coordinate surfaces===
 
The surfaces of constant μ form [[Oblate spheroid|oblate]] [[spheroids]], by the trigonometric identity
 
:<math>
\frac{x^{2} + y^{2}}{a^{2} \cosh^{2} \mu} +
\frac{z^{2}}{a^{2} \sinh^{2} \mu} = \cos^{2} \nu + \sin^{2} \nu = 1
</math>
 
since they are [[ellipse]]s rotated about the ''z''-axis, which separates their foci.   An ellipse in the ''x''-''z'' plane (Figure 2) has a [[semi-major axis|major semiaxis]] of length ''a'' cosh μ along the ''x''-axis, whereas its [[semi-minor axis|minor semiaxis]] has length ''a'' sinh μ along the ''z''-axis.  The foci of all the ellipses in the ''x''-''z'' plane are located on the ''x''-axis at ±''a''.
 
Similarly, the surfaces of constant ν form one-sheet half [[hyperboloid]]s of revolution by the hyperbolic trigonometric identity
 
:<math>
\frac{x^{2} + y^{2}}{a^{2} \cos^{2} \nu} -
\frac{z^{2}}{a^{2} \sin^{2} \nu} = \cosh^{2} \mu - \sinh^{2} \mu = 1
</math>
 
For positive ν, the half-hyperboloid is above the ''x''-''y'' plane (i.e., has positive ''z'') whereas for negative ν, the half-hyperboloid is below the ''x''-''y'' plane (i.e., has negative ''z'').  Geometrically, the angle ν corresponds to the angle of the [[asymptote]]s of the hyperbola.  The foci of all the hyperbolae are likewise located on the ''x''-axis at ±''a''.
 
===Inverse transformation===
 
The (μ, ν, φ) coordinates may be calculated from the Cartesian coordinates (''x'', ''y'', ''z'') as follows.  The azimuthal angle φ is given by the formula
 
:<math>
\tan \phi = \frac{y}{x}
</math>
 
The cylindrical radius ρ of the point P is given by
 
:<math>
\rho^{2} = x^{2} + y^{2}
</math>
 
and its distances to the foci in the plane defined by φ is given by
 
:<math>
d_{1}^{2} = (\rho + a)^{2} + z^{2}
</math>
 
:<math>
d_{2}^{2} = (\rho - a)^{2} + z^{2}
</math>
 
The remaining coordinates μ and ν can be calculated from the equations
 
:<math>
\cosh \mu = \frac{d_{1} + d_{2}}{2a}
</math>
 
:<math>
\cos \nu = \frac{d_{1} - d_{2}}{2a}
</math>
 
where the sign of μ is always non-negative, and the sign of ν is the same as that of ''z''.
 
Another method to compute the inverse transform is
 
:<math>
\mu = \operatorname{Re} \operatorname{arccosh} \frac{\rho + z i}{a}
</math>
 
:<math>
\nu = \operatorname{Im} \operatorname{arccosh} \frac{\rho + z i}{a}
</math>
 
:<math>
\phi = \arctan \frac{y}{x}
</math>
 
where
 
:<math>
\rho = \sqrt{x^2 + y^2}
</math>
 
===Scale factors===
 
The scale factors for the coordinates μ and ν are equal
 
:<math>
h_{\mu} = h_{\nu} = a\sqrt{\sinh^{2}\mu + \sin^{2}\nu}
</math>
 
whereas the azimuthal scale factor equals
 
:<math>
h_{\phi} = a \cosh\mu \ \cos\nu
</math>
 
Consequently, an infinitesimal volume element equals
 
:<math>
dV = a^{3} \cosh\mu \ \cos\nu \
\left( \sinh^{2}\mu + \sin^{2}\nu \right) d\mu d\nu d\phi
</math>
 
and the Laplacian can be written
 
:<math>
\nabla^{2} \Phi =
\frac{1}{a^{2} \left( \sinh^{2}\mu + \sin^{2}\nu \right)}
\left[
\frac{1}{\cosh \mu} \frac{\partial}{\partial \mu}
\left( \cosh \mu \frac{\partial \Phi}{\partial \mu} \right) +
\frac{1}{\cos \nu} \frac{\partial}{\partial \nu}
\left( \cos \nu \frac{\partial \Phi}{\partial \nu} \right)
\right] +
\frac{1}{a^{2} \left( \cosh^{2}\mu+\cos^{2}\nu \right)}
\frac{\partial^{2} \Phi}{\partial \phi^{2}}
</math>
 
Other differential operators such as <math>\nabla \cdot \mathbf{F}</math> and <math>\nabla \times \mathbf{F}</math> can be expressed in the coordinates (μ, ν, φ) by substituting the scale factors into the general formulae found in [[orthogonal coordinates]].
 
===Basis Vectors===
 
The orthonormal basis vectors for the <math>\mu,\nu,\phi</math> coordinate system can be expressed in Cartesian coordinates as
 
:<math>
\hat{e}_{\mu} = \frac{1}{\sqrt{\sinh^2 \mu + \sin^2 \nu}}
\left( \sinh \mu \cos \nu \cos \phi \boldsymbol{\hat{i}} + \sinh \mu \cos \nu \sin \phi \boldsymbol{\hat{j}} + \cosh \mu \sin \nu \boldsymbol{\hat{k}}\right)
</math>
 
:<math>
\hat{e}_{\nu} = \frac{1}{\sqrt{\sinh^2 \mu + \sin^2 \nu}}
\left(
-\cosh \mu \sin \nu \cos \phi \boldsymbol{\hat{i}} - \cosh \mu \sin \nu \sin \phi \boldsymbol{\hat{j}} + \sinh \mu \cos \nu \boldsymbol{\hat{k}}
\right)
</math>
 
:<math>
\hat{e}_{\phi} = -\sin \phi \boldsymbol{\hat{i}} + \cos \phi \boldsymbol{\hat{j}}
</math>
 
where <math>\boldsymbol{\hat{i}}, \boldsymbol{\hat{j}}, \boldsymbol{\hat{k}}</math> are the Cartesian unit vectors. Here, <math>\hat{e}_{\mu}</math> is the outward normal vector to the oblate spheroidal surface of constant <math>\mu</math>, <math>\hat{e}_{\phi}</math> is the same azimuthal unit vector from spherical coordinates, and <math>\hat{e}_{\nu}</math> lies in the tangent plane to the oblate spheroid surface and completes the right-handed basis set.
 
==Definition (ζ, ξ, φ)==
 
Another set of oblate spheroidal coordinates <math>(\zeta,\xi,\phi)</math> are sometimes used where <math>\zeta = \sinh \mu</math> and <math>\xi = \sin \nu</math> (Smythe 1968).  The curves of constant <math>\zeta</math> are oblate spheroids and the curves of constant <math>\xi</math> are the hyperboloids of revolution. The coordinate <math>\zeta</math> is restricted by <math>0 \le \zeta < \infty</math> and <math>\xi</math> is restricted by <math>-1 \le \xi < 1</math>.
 
The relationship to [[Cartesian coordinates]] is
 
:<math>x = a\sqrt{(1+\zeta^2)(1-\xi^2)}\,\cos \phi\,</math>
 
:<math>y = a\sqrt{(1+\zeta^2)(1-\xi^2)}\,\sin \phi\,</math>
 
:<math>z  = a \zeta \xi\,</math>
 
===Scale factors===
 
The scale factors for <math>(\zeta, \xi, \phi)</math> are:
 
:<math>
h_{\zeta} = a\sqrt{\frac{\zeta^2 + \xi^2}{1+\zeta^2}}
</math>
 
:<math>
h_{\xi} = a\sqrt{\frac{\zeta^2 + \xi^2}{1 - \xi^2}}
</math>
 
:<math>
h_{\phi} = a\sqrt{(1+\zeta^2)(1 - \xi^2)}
</math>
 
Knowing the scale factors, various functions of the coordinates can be calculated by the general method outlined in the [[orthogonal coordinates]] article. The infinitesimal volume element is:
 
:<math>
dV = a^{3} (\zeta^2+\xi^2)\,d\zeta\,d\xi\,d\phi
</math>
 
The gradient is:
:<math>
\nabla V =
\frac{1}{h_{\zeta}} \frac{\partial V}{\partial \zeta} \,\hat{\zeta}+
\frac{1}  {h_{\xi}} \frac{\partial V}{\partial \xi}  \,\hat{\xi}+
\frac{1} {h_{\phi}} \frac{\partial V}{\partial \phi}  \,\hat{\phi}
</math>
 
The divergence is:
 
:<math>
\nabla \mathbf{F} = \frac{1}{a(\zeta^2+\xi^2)}
\left\{
\frac{\partial}{\partial \zeta} \left(\sqrt{1+\zeta^2}\sqrt{\zeta^2+\xi^2}F_\zeta\right)        +
\frac{\partial}  {\partial \xi} \left(\sqrt{1-\xi^2}\sqrt{\zeta^2+\xi^2}F_\xi\right)
\right\}
+\frac{1}{\sqrt{1+\zeta^2}\sqrt{1-\xi^2}} \frac{\partial F_\phi}{\partial \phi}
</math>
 
and the Laplacian equals
 
:<math>
\nabla^{2} V =
\frac{1}{a^2 \left( \zeta^2 + \xi^2 \right)}
\left\{
\frac{\partial}{\partial \zeta} \left[
\left(1+\zeta^2\right) \frac{\partial V}{\partial \zeta}
\right] +
\frac{\partial}{\partial \xi} \left[
\left( 1 - \xi^2 \right) \frac{\partial V}{\partial \xi}
\right]
\right\}
+ \frac{1}{a^2 \left( 1+\zeta^2 \right) \left( 1 - \xi^{2} \right)}
\frac{\partial^2 V}{\partial \phi^{2}}
</math>
 
=== Oblate spheroidal harmonics ===
 
: ''See also [[Oblate spheroidal wave function]].''
 
As is the case with [[spherical coordinates]] and [[spherical harmonics]], Laplace's equation may be solved by the method of [[separation of variables]] to yield solutions in the form of '''oblate spheroidal harmonics''', which are convenient to use when boundary conditions are defined on a surface with a constant oblate spheroidal coordinate.  
 
Following the technique of [[separation of variables]], a solution to Laplace's equation is written:
 
:<math>V=Z(\zeta)\,\Xi(\xi)\,\Phi(\phi)\,</math>
 
This yields three separate differential equations in each of the variables:
 
:<math>\frac{d}{d\zeta}\left[(1+\zeta^2)\frac{dZ  }{d\zeta}\right]+\frac{m^2Z  }{1+\zeta^2}-n(n+1)Z  =0</math>
 
:<math>\frac{d}{d\xi  }\left[(1-\xi^2  )\frac{d\Xi}{d\xi  }\right]-\frac{m^2\Xi}{1-\xi^2  }+n(n+1)\Xi=0</math>
 
:<math>\frac{d^2\Phi}{d\phi^2}=-m^2\Phi</math>
 
where ''m'' is a constant which is an integer because the &phi; variable is periodic with period 2&pi;. ''n'' will then be an integer. The solution to these equations are:
 
:<math>Z_{mn}  =A_1P_n^m(i\zeta)+A_2Q_n^m(i\zeta)</math>
 
:<math>\Xi_{mn} =A_3P_n^m(\xi)+A_4Q_n^m(\xi)</math>
 
:<math>\Phi_m  =A_5e^{im\phi}+A_6e^{-im\phi}\,</math>
 
where the <math>A_i</math> are constants and <math>P_n^m(z)</math> and <math>Q_n^m(z)</math> are [[associated Legendre polynomial]]s of the first and second kind respectively. The product of the three solutions is called an ''oblate spheroidal harmonic'' and the general solution to Laplace's equation is written:
 
:<math>V=\sum_{n=0}^\infty\sum_{m=0}^\infty\,Z_{mn}(\zeta)\,\Xi_{mn}(\xi)\,\Phi_m(\phi)</math>
 
The constants will combine to yield only four independent constants for each harmonic.
 
==Definition (σ, τ, φ)==
 
[[File:Oblate spheroidal coordinates full hyperboloid.png|thumb|right|350px|Figure 3: Coordinate isosurfaces for a point P (shown as a black sphere) in the alternative oblate spheroidal coordinates (σ, τ, φ). As before, the oblate spheroid corresponding to σ is shown in red, and φ measures the azimuthal angle between the green and yellow half-planes. However, the surface of constant τ is a full one-sheet hyperboloid, shown in blue. This produces a two-fold degeneracy, shown by the two black spheres located at (''x'', ''y'', ±''z'').]]
 
An alternative and geometrically intuitive set of oblate spheroidal coordinates (σ, τ, φ) are sometimes used, where σ = [[cosh]] μ and τ = [[cosine|cos]] ν.<ref>Abramowitz and Stegun, p. 752.</ref>  Therefore, the coordinate σ must be greater than or equal to one, whereas τ must lie between ±1, inclusive.  The surfaces of constant σ are oblate spheroids, as were those of constant μ, whereas the curves of constant τ are full hyperboloids of revolution, including the half-hyperboloids corresponding to ±ν.  Thus, these coordinates are degenerate; ''two'' points in Cartesian coordinates (''x'', ''y'', ±''z'') map to ''one'' set of coordinates (σ, τ, φ).  This two-fold degeneracy in the sign of ''z'' is evident from the equations transforming from oblate spheroidal coordinates to the [[Cartesian coordinates]]
 
:<math>
x = a\sigma\tau \cos \phi\,
</math>
 
:<math>
y = a\sigma\tau \sin \phi\,
</math>
 
:<math>
z^{2}  = a^{2} \left( \sigma^{2} - 1 \right) \left(1 - \tau^{2} \right)
</math>
 
The coordinates <math>\sigma</math> and <math>\tau</math> have a simple relation to the distances to the focal ring.  For any point, the ''sum'' <math>d_{1}+d_{2}</math> of its distances to the focal ring equals <math>2a\sigma</math>, whereas their ''difference'' <math>d_{1}-d_{2}</math> equals <math>2a\tau</math>.  Thus, the "far" distance to the focal ring is <math>a(\sigma+\tau)</math>, whereas the "near" distance is <math>a(\sigma-\tau)</math>.
 
===Coordinate surfaces===
 
Similar to its counterpart μ, the surfaces of constant σ form [[Oblate spheroid|oblate]] [[spheroids]]
 
:<math>
\frac{x^{2} + y^{2}}{a^{2} \sigma^{2}} +
\frac{z^{2}}{a^{2} \left(\sigma^{2} -1\right)} = 1
</math>
 
Similarly, the surfaces of constant τ form full one-sheet [[hyperboloid]]s of revolution
 
:<math>
\frac{x^{2} + y^{2}}{a^{2} \tau^{2}} -
\frac{z^{2}}{a^{2} \left( 1 - \tau^{2} \right)} = 1
</math>
 
===Scale factors===
 
The scale factors for the alternative oblate spheroidal coordinates <math>(\sigma, \tau, \phi)</math> are
 
:<math>
h_\sigma = a\sqrt{\frac{\sigma^2 - \tau^2}{\sigma^2 - 1}}
</math>
 
:<math>
h_\tau = a\sqrt{\frac{\sigma^2 - \tau^2}{1 - \tau^{2}}}
</math>
 
whereas the azimuthal scale factor is <math>h_{\phi} = a \sigma \tau</math>.
 
Hence, the infinitesimal volume element can be written
 
:<math>
dV = a^3 \sigma \tau \frac{\sigma^2 - \tau^2}{\sqrt{\left( \sigma^2 - 1 \right) \left( 1 - \tau^2 \right)}} \, d\sigma \, d\tau \, d\phi
</math>
 
and the Laplacian equals
 
:<math>
\nabla^2 \Phi =
\frac{1}{a^2 \left( \sigma^2 - \tau^2 \right)}
\left\{
\frac{\sqrt{\sigma^2 -1}}{\sigma}
\frac{\partial}{\partial \sigma} \left[
\left( \sigma\sqrt{\sigma^2 - 1} \right) \frac{\partial \Phi}{\partial \sigma}
\right] +
\frac{\sqrt{1 - \tau^2}}{\tau}
\frac{\partial}{\partial \tau} \left[
\left( \tau\sqrt{1 - \tau^2} \right) \frac{\partial \Phi}{\partial \tau}
\right]
\right\}
+ \frac{1}{a^2 \sigma^2 \tau^2 }
\frac{\partial^2 \Phi}{\partial \phi^2}
</math>
 
Other differential operators such as <math>\nabla \cdot \mathbf{F}</math>  and <math>\nabla \times \mathbf{F}</math> can be expressed in the coordinates <math>(\sigma, \tau)</math> by substituting  the scale factors into the general formulae  found in [[orthogonal coordinates]].
 
As is the case with [[spherical coordinates]], Laplaces equation may be solved by the method of [[separation of variables]] to yield solutions in the form of '''oblate spheroidal harmonics''', which are convenient to use when boundary conditions are defined on a surface with a constant oblate spheroidal coordinate (See Smythe, 1968).
 
==References==
{{reflist}}
 
==Bibliography==
===No angles convention===
*{{cite book | author = Morse PM, Feshbach H | year = 1953 | title = Methods of Theoretical Physics, Part I | publisher = McGraw-Hill | location = New York | page = 662}}  Uses ξ<sub>1</sub> = a sinh μ,  ξ<sub>2</sub> = sin ν, and ξ<sub>3</sub> = cos φ.
*{{cite book | author = Zwillinger D | year = 1992 | title = Handbook of Integration | publisher = Jones and Bartlett | location = Boston, MA | isbn = 0-86720-293-9 | page = 115}}  Same as Morse & Feshbach (1953), substituting ''u''<sub>''k''</sub> for ξ<sub>''k''</sub>.
*{{cite book | last = Smythe | first = WR| title = Static and Dynamic Electricity |edition = 3rd | publisher = McGraw-Hill | location = New York | year = 1968}}
*{{cite book | author = Sauer R, Szabó I | year = 1967 | title = Mathematische Hilfsmittel des Ingenieurs | publisher = Springer Verlag | location = New York | page = 98 | lccn = 6725285}}  Uses hybrid coordinates ξ = sinh μ,  η = sin ν, and φ.
 
===Angle convention===
*{{cite book | author = Korn GA, Korn TM |year = 1961 | title = Mathematical Handbook for Scientists and Engineers | publisher = McGraw-Hill | location = New York | page = 177 | lccn = 5914456}}  Korn and Korn use the (μ, ν, φ) coordinates, but also introduce the degenerate (σ, τ, φ) coordinates.
*{{cite book | author = Margenau H, Murphy GM | year = 1956 | title = The Mathematics of Physics and Chemistry | publisher = D. van Nostrand | location = New York| page = 182 | lccn = 5510911 }}  Like Korn and Korn (1961), but uses [[colatitude]] θ = 90° - ν instead of [[latitude]] ν.
*{{cite book | author = Moon PH, Spencer DE | year = 1988 | chapter = Oblate spheroidal coordinates (η, θ, ψ) | title = Field Theory Handbook, Including Coordinate Systems, Differential Equations, and Their Solutions | edition = corrected 2nd ed., 3rd print | publisher = Springer Verlag | location = New York | isbn = 0-387-02732-7 | pages = 31&ndash;34 (Table 1.07)}}  Moon and Spencer use the colatitude convention θ = 90° - ν, and rename φ as ψ.
 
===Unusual convention===
*{{cite book | author = Landau LD, Lifshitz EM, Pitaevskii LP | year = 1984 | title = Electrodynamics of Continuous Media (Volume 8 of the [[Course of Theoretical Physics]]) | edition = 2nd | publisher = Pergamon Press | location = New York | isbn = 978-0-7506-2634-7 | pages = 19&ndash;29 }}  Treats the oblate spheroidal coordinates as a limiting case of the general [[ellipsoidal coordinates]].  Uses (ξ, η, ζ) coordinates that have the units of distance squared.
 
==External links==
*[http://mathworld.wolfram.com/OblateSpheroidalCoordinates.html MathWorld description of oblate spheroidal coordinates]
 
{{Orthogonal coordinate systems}}
 
[[Category:Coordinate systems]]

Latest revision as of 05:33, 25 August 2014

The writer is called Irwin Wunder but it's not the most masucline name out there. One of the things she enjoys most is to do aerobics and now she is trying to earn cash with it. Hiring has been my profession for some time but I've already applied for another 1. North Dakota is where me and my husband live.

Also visit my webpage rivoli.enaiponline.com