[Date Index]
[Thread Index]
[Author Index]
Re: cross-product in cylindrical problem
*To*: mathgroup at smc.vnet.net
*Subject*: [mg51263] Re: [mg51258] cross-product in cylindrical problem
*From*: Andrzej Kozlowski <akoz at mimuw.edu.pl>
*Date*: Mon, 11 Oct 2004 01:25:22 -0400 (EDT)
*References*: <200410100952.FAA25275@smc.vnet.net>
*Sender*: owner-wri-mathgroup at wolfram.com
On 10 Oct 2004, at 18:52, news wrote:
>
> I'm really puzzled by this behavior of Mathematica, I have two vectors
> in
> cylindrical coordinates and would like to take their cross-product in
> cylindrical, but it seems to give me incorrect answer, see below:
>
> define parametric path {r,phi,z}
>
> In[110]:=
> f[\[Rho]_, \[Phi]_] = {\[Rho], \[Phi], 0}
> Out[110]=
> {\[Rho], \[Phi], 0}
>
> take derivates of path w.r.t. r then w.r.t phi, get {1,0,0}, and
> {0,1,0}
>
> In[113]:=
> v1 = D[f[\[Rho], \[Phi]], \[Rho]]
> v2 = D[f[\[Rho], \[Phi]], \[Phi]]
> Out[113]=
> {1, 0, 0}
> Out[114]=
> {0, 1, 0}
>
> then cross them in cylindrical coords, and should get {0,0,1}, but
> instead
> get wrong answer below
>
> In[117]:=
> n = CrossProduct[v1, v2, Cylindrical[\[Rho], \[Phi], z]] //
> FullSimplify
> Out[117]=
> {0, 0, 0}
>
> As you can see, when I cross {1,0,0} with {0,1,0} in cylindrical
> coords, I
> get {0,0,0}, when I should be getting {0,0,1}.
>
> Can anyone help?
>
>
The problem is with the meaning of "two vectors in cylindrical
coordinates". Even if your default coordinate system is cylindrical,
Mathematica still represents tangent vectors using Cartesian
coordinates. In other words, the meaning of a vector {a,b,c} in
cylindrical coordinates is the vector given by:
<< Calculus`VectorAnalysis`
CoordinatesToCartesian[{a,b,c},Cylindrical]
Thus, the vectors {1,0,0} and {0,1,0} are turned to:
CoordinatesToCartesian[{1, 0, 0}, Cylindrical]
{1, 0, 0}
and
CoordinatesToCartesian[{0,1,0},Cylindrical]
{0,0,0}
So you get the zero vector (since {0,1,0} in cylindrical coordinates
represents just the origin) and hence the result you get.
I am not hundred percent sure what you meant. But I suppose thatt by
{1,0,0} and {0,1,0} in cylindrical coordinates you meant the basic
vector fields in the tangent bundle of R^3 (in other words vectors
changing from point to point). Presumably (I am guessing) the vector
(field) {1,0,0} is
JacobianMatrix[] . {1, 0, 0}
{Cos[Ï?], Sin[Ï?], 0}
and the vector (field) {0,1,0} is
JacobianMatrix[] . {0, 1, 0}
{(-Ï?)*Sin[Ï?], Ï?*Cos[Ï?], 0}
in Cartesian coordinates.
In that case
Simplify[Cross[JacobianMatrix[] . {1, 0, 0}, JacobianMatrix[] . {0, 1,
0}]]
{0, 0, Ï?}
which is the vector field sometimes denoted by {0,0,1} in cylindrical
coordinates.
Andrzej Kozlowski
Chiba, Japan
http://www.akikoz.net/~andrzej/
http://www.mimuw.edu.pl/~akoz/
Prev by Date:
**derivatives in cylindrical coord's, need help**
Next by Date:
**Re: plotting groups of polynomial roots**
Previous by thread:
**cross-product in cylindrical problem**
Next by thread:
**Re: cross-product in cylindrical problem**
| |