I Force fields in curvilinear coordinate systems

AI Thread Summary
The discussion focuses on calculating work done by a force along paths in cylindrical and spherical coordinates, highlighting difficulties in parametrizing curves compared to Cartesian coordinates. The user successfully applies the Cartesian method but struggles with cylindrical parametrization, particularly when transitioning between different paths like quarter circles. Clarification is provided on the distinction between linear and curvilinear coordinates, emphasizing the need to express positions and forces in terms of standard or normalized bases. The importance of considering changing basis vectors when dealing with parametric curves in curvilinear coordinates is also discussed. Understanding these concepts is crucial for correctly applying integrals in non-linear coordinate systems.
calculo2718
Messages
25
Reaction score
1
I am trying to solve problems where I calculate work do to force along paths in cylindrical and spherical coordinates.

I can do almost by rote the problems in Cartesian: given a force ##\vec{F} = f(x,y,z)\hat{x} + g(x,y,z)\hat{y}+ h(x,y,z)\hat{z}## I can parametricize my some curve ##\gamma: \vec{r}(t) = x(t)\hat{x} + y(t)\hat{y} + z(t)\hat{z}## then take ##\int_\gamma \vec{F}\cdot\vec{r}'(t)dt ##

I know it should work similarly in cylindrical, if I have ##\vec{F}(R,\phi) = f(R,\phi)\hat{R} + g(R,\phi)\hat{\phi}##, The problem is I am having trouble thinking through the curve parametrization in cylindrical and similarly in spherical. For example, If I want to integrate along the path where ##\phi = 0##, from some radius ##r_1## to ##r_2## I would say
##\gamma_1: \vec{r}(t) = t\hat{R} + 0\hat{\phi}##

If I then want to go from along the quarter circle path from ##\phi = 0## to ##\phi = \frac{\pi}{2}## then

##\gamma_2: \vec{r}(t) = r_2\hat{R} + r_2\cdot t\hat{\phi}##

If I want to go back down to ##r_1##
##\gamma_3: \vec{r}(t) = t\hat{R} + \frac{\pi}{2}\hat{\phi}##

This is where my confusion gets in the way, the ##\gamma_2## curve makes it look like I can in general parametricize a curve in cylindrical as ##\vec{r}(t) = R(t)\hat{R} + \phi(t)R(t)\hat{\phi}## where on this curve ##R(t) = r_2## and ##\phi(t) = t##, however the ##\gamma_3## curve is not of this form and if I try to parametricize it that way(##\vec{r}(t) = t\hat{R} + \frac{\pi}{2}t\hat{\phi}##) it is incorrect(I tested it with a Force only in the ##\hat{\phi}## direction which should give me 0, but I get a non-0 integral).

How can I parametricize a general curve in cylindrical and spherical the way I have done above for cartesian? What am I missing when trying to do these problems?

In addition how can I derive from my ##\vec{r}(t)## in cylindrical the line element in cylindrical ##d\vec{r} = dR\hat{R} + rd\phi\hat{\phi}##. I feel like understand how this is done will help me see more clearly why I am confused
 
Physics news on Phys.org
I am going to give you a quick general overview which I think will answer most of your questions:

To start there are two things going on when we speak in terms of coordinate systems. There are the linear coordinate systems, including Cartesian coordinates where the coordinates are the components of a position vector in a specific coordinate system, \vec{r}=x\hat{\imath} + y\hat{\jmath}+z\hat{k}.

You then have non-linear (curvilinear) coordinates (polar, cylindrical, spherical, etc). When using general curvilinear coordinates you should express positions either as coordinate positions (coordinate tuples which are not considered a vector) or using the standard basis which is effectively giving the cartesian coordinate as function of these non-linear coordinates. Example:
\vec{r} = R\cos(\phi)\hat{\imath} + R\sin(\phi)\hat{\jmath}+z\hat{k}
in cylindrical coordinates (R,\phi,z)

Now if you are considering other vectors besides the position vector (such as velocity or force or differential position for a parameterized curve) then you can express it either in the standard basis but with components still a function of the curvilinear coordinates, or in a coordinate basis, or in a normalized coordinate basis.

Example: In a (2-dimensional) central force problem where the force is toward the origin and with magnitude k/r^2 where r is distance from the origin. (I am going to use (r,\theta) polar coordinates in the plain since that's the familiar convention to me.)

I can express my force vector in purely rectangular coordinates and in the standard basis as:
\vec{F} = \frac{-k x}{(x^2+y^2)^{3/2}}\hat{\imath}+\frac{-k y}{(x^2+y^2)^{3/2}}\hat{\imath}
of course this is rather awkward and we can instead express the force in terms of the polar coordinates but still in the standard basis:
\vec{F} = \frac{-k \cos(\theta)}{r^2}\hat{\imath}+\frac{-k \sin(\theta)}{r^2}\hat{\imath}
even better we might decide to use a normalized coordinate basis. By a coordinate basis we mean that we use the gradient of each coordinate as our basis. These are not necessarily unit vectors so if we want a normalized basis we would then make each into a unit vector to get a normalized coordinate basis (which is then not the coordinate basis as such!).
In polar coordinates:
r = \sqrt{x^2 + y^2} , \nabla r = \frac{x}{r} \hat{\imath}+ {y}{r}\hat{\jmath} = \hat{r}
it happens to already be a unit vector. In fully polar coordinates it is \nabla r = \hat{r} = \cos(\theta)\hat{\imath} + \sin(\theta)\hat{\jmath}.
\theta = \tan^{-1}(y/x) + \text{const}, \nabla \theta = \frac{-y}{r^2}\hat{\imath} + \frac{x}{r^2}\hat{\jmath}
it is not a unit vector but r\nabla \theta is so we can use:
\hat{\theta} = r\nabla\theta = \frac{-y}{r}\hat{\imath} + \frac{x}{r}\hat{\jmath} = -\sin(\theta)\hat{\imath} + \cos(\theta)\hat{\jmath}

We can now express the force above in both polar coordinates and a normalized coordinate basis:
\vec{F} = \frac{-k}{r}\hat{r}

Note that the curvilinear coordinate basis and normalize coordinate basis can in general depend on the coordinates which can complicate mightily the differential calculus. That's a whole chapter of a textbook by itself. But for the parametric curve case it is not so bad.

If you are dealing with a parametric curve you are giving your position coordinates (in whatever coordinate system) as functions of the parameter. But do remember that for curvilinear case where you are using a coordinate or normalized coordinate basis, you must remember that the basis vectors are also changing with the parameter. You simply apply the product rule. E.g.

\dot{\vec{F}} = \frac{d}{dt}\left(\frac{-k}{r^2}\right)\hat{r} + \frac{-k}{r^2}\frac{d}{dt}\hat{r} = ...

I hope this has clarified the subject for you a bit.
 
thank you! this was helpful.
 
Thread 'Gauss' law seems to imply instantaneous electric field propagation'
Imagine a charged sphere at the origin connected through an open switch to a vertical grounded wire. We wish to find an expression for the horizontal component of the electric field at a distance ##\mathbf{r}## from the sphere as it discharges. By using the Lorenz gauge condition: $$\nabla \cdot \mathbf{A} + \frac{1}{c^2}\frac{\partial \phi}{\partial t}=0\tag{1}$$ we find the following retarded solutions to the Maxwell equations If we assume that...
Maxwell’s equations imply the following wave equation for the electric field $$\nabla^2\mathbf{E}-\frac{1}{c^2}\frac{\partial^2\mathbf{E}}{\partial t^2} = \frac{1}{\varepsilon_0}\nabla\rho+\mu_0\frac{\partial\mathbf J}{\partial t}.\tag{1}$$ I wonder if eqn.##(1)## can be split into the following transverse part $$\nabla^2\mathbf{E}_T-\frac{1}{c^2}\frac{\partial^2\mathbf{E}_T}{\partial t^2} = \mu_0\frac{\partial\mathbf{J}_T}{\partial t}\tag{2}$$ and longitudinal part...
Thread 'Recovering Hamilton's Equations from Poisson brackets'
The issue : Let me start by copying and pasting the relevant passage from the text, thanks to modern day methods of computing. The trouble is, in equation (4.79), it completely ignores the partial derivative of ##q_i## with respect to time, i.e. it puts ##\partial q_i/\partial t=0##. But ##q_i## is a dynamical variable of ##t##, or ##q_i(t)##. In the derivation of Hamilton's equations from the Hamiltonian, viz. ##H = p_i \dot q_i-L##, nowhere did we assume that ##\partial q_i/\partial...
Back
Top