Unit Normal Vector in Acceleration

[Philosophaie]

What is the Unit Normal Vector in Acceleration.$$\vec{u}_n$$ How would I go about calculating the normal. Wiki and Wolfram Alpha were very vague about calculating the normal of an orbit. I would like to know the normal to the orbit of a mass orbiting a much larger mass not counting any other gravitating masses in proximity.

Magnitude$$v = \sqrt{vx^2 + vy^2 + vz^2}$$
Unit Tangent Vector$$\vec{u}_v =(\frac{vx}{v} ,\frac{vy}{v} ,\frac{vz}{v})$$
Time derivative of velocity$$\frac{dv}{dt} =(\frac{(\frac{d(vx^2)}{dt} + \frac{d(vy^2)}{dt} + \frac{d(vz^2)}{dt})}{v})$$
Acceleration$$\vec{a} = \frac{dv}{dt}*\vec{u}_v + \frac{v^2}{r}*\vec{u}_n$$
How do you calculate$$\vec{u}_n$$

 

[Chestermiller]

$\vec{u_n}$ is equal to $\frac{d\vec{u_v}}{dt}$ divided by the magnitude of $\frac{d\vec{u_v}}{dt}$.

 

[ZapperZ]

If this is an orbit in a plane (which most orbits and conic section motion are), why aren't making your life easier by doing this in plane polar coordinates?

Zz.

 

[Philosophaie]

Now to acceleration, is the equation correct. My data is not coming out right.$$\vec{a} = \frac{dv}{dt} * \vec{u}_v + \frac{v^2}{r} * \frac{\frac{d(\vec{v})}{dt}}{|\frac{d(\vec{v})}{dt}|}$$I am interested in the centripetal acceleration!

 

[D H]

Now to acceleration, is the equation correct my data is not coming out right.$$\vec{a} = \frac{dv}{dt} * \vec{u}_v + \frac{v^2}{r} * \frac{\frac{d(\vec{v})}{dt}}{|\frac{d(\vec{v})}{dt}|}$$

That is not correct. You are writing velocity as a scalar times the unit vector in the direction of the velocity vector, $\vec v = v\hat u_v$. Here $\hat u_v \equiv \vec v / ||\vec v||$. Differentiating with respect to time yields $\vec a = \dot v \vec u _v + v\dot{\hat u}_v$. Your second term is not the same as $v\dot{\hat u}_v$.

I am interested in the centripetal acceleration!

Why? It's pretty much meaningless. You want the central force acceleration, and that's given by Newton's law of gravitation. You can compute centripetal acceleration by subtracting the tangential component ($v\hat u_v$) from the gravitational acceleration, but why bother?

 

[Astrum]

If this is an orbit in a plane (which most orbits and conic section motion are), why aren't making your life easier by doing this in plane polar coordinates?

Zz.

I suggested the same thing on his/her last two or three posts, for some reason, he/she refuses to listen.

Let me try again, maybe he/she will hear me this time.

USE POLAR COORDINATES, THEY WERE CREATED EXACTLY FOR PROBLEMS LIKE THIS

 

[Chestermiller]

There is a much easier way of getting $\vec{u_n}$. If $\vec{a}$ is the acceleration vector, then the component of $\vec{a}$ in the direction of the velocity vector is $\vec{a}\centerdot \vec{u_v}$. So the vector component of the acceleration vector in the direction perpendicular to the velocity vector is $\vec{a}-(\vec{a}\centerdot \vec{u_v})\vec{u_v}$. The unit normal vector $\vec{u_n}$ is equal to this vector divided by its own magnitude:
$$\vec{u_n}=\frac{\vec{a}-(\vec{a}\centerdot \vec{u_v})\vec{u_v}}{\vert\vec{a}-(\vec{a}\centerdot \vec{u_v})\vec{u_v}\vert}$$
Now, the denominator is given by
$$\vert\vec{a}-(\vec{a}\centerdot \vec{u_v})\vec{u_v}\vert=\sqrt{a^2-(\vec{a}\centerdot \vec{u_v})^2}$$
So,
$$\vec{u_n}=\frac{\vec{a}-(\vec{a}\centerdot \vec{u_v})\vec{u_v}}{\sqrt{a^2-(\vec{a}\centerdot \vec{u_v})^2}}$$
Finally, since $\vec{u_v}=\frac{\vec{v}}{v}$
$$\vec{u_n}=\frac{v^2\vec{a}-(\vec{a}\centerdot \vec{v})\vec{v}}{v\sqrt{(av)^2-(\vec{a}\centerdot \vec{v})^2}}$$
This provides a simple expression for getting the unit normal to the trajectory even when one is using Cartesian coordinates.

 

[D H]

This provides a simple expression for getting the unit normal to the trajectory even when one is using Cartesian coordinates.

Even easier, use the triple product $\frac {\vec v \times \vec a \times \vec v}{v^2}$.

However, this vector has zero physical meaning. So why bother?

 

[Chestermiller]

Even easier, use the triple product $\frac {\vec v \times \vec a \times \vec v}{v^2}$.

However, this vector has zero physical meaning. So why bother?

I like your thinking. Very clever and very cute. Cudos!

There is a bit of confusion, however. Your final result is not the unit normal. The denominator is not equal to the magnitude of the vector in the numerator. I don't know whether you intended this or not. Actually, the denominator in my expression is equal to the magnitude of the vector in the numerator of your expression. That must mean that our two numerators must be mathematically equivalent. I have not tried to verify this, but it should be possible by expressing the results in component form. But, there must be a more elegant way than this.

I agree that the unit normal does not have much physical meaning. What really is important is the curvature vector:

$$\frac{d\vec{u_v}}{ds}=\frac{\vec{u_n}}{r}$$
where r is the local radius of curvature of the trajectory and ds is differential position along the trajectory: ds=vdt.

 

[D H]

You're right. My expression is the component of acceleration that is normal to the velocity vector. Easy enough to normalize, for all that's worth.

Even the curvature isn't all that useful for orbits. At least I've never found it to be.

 

[jhae2.718]

For an orbit, as others have said, the easiest thing to use is polar coordinates. First construct a reference frame fixed to the large body that you'll use to approximate an inertial frame, and then construct a rotating reference frame pointing toward the orbiting satellite. Suppose the inertial frame is defined by the three orthonormal vectors $\hat{\boldsymbol \imath}_1$, $\hat{\boldsymbol \imath}_2$, and $\hat{\boldsymbol \imath}_3$, and that the rotating frame is defined by $\hat{\boldsymbol e}_1$, $\hat{\boldsymbol e}_2$, and $\hat{\boldsymbol e}_3$.

If we say that $\hat{\boldsymbol e}_1$ points to the mass center of the satellite, we can define the inertial position vector as $${\boldsymbol p} = r \hat{\boldsymbol e}_1.$$ Now, to make calculating the inertial velocity and acceleration in the ${\boldsymbol e}^+$ frame, we need a result called the kinematic transport theorem. Consider two reference frames ${\boldsymbol a}^+$ and ${\boldsymbol b}^+$, where the ${\boldsymbol b}^+$ frame is rotating with respect to the ${\boldsymbol a}^+$ frame with angular velocity ${\boldsymbol \omega}_\text{a/b}$. The derivative of a vector ${\boldsymbol p}$ in the ${\boldsymbol a}^+$ frame, coordinatized in the ${\boldsymbol b}^+$ frame, comes from two components:

1. The rate of change of ${\boldsymbol p}$ in the ${\boldsymbol b}^+$ frame
2. The angular velocity of ${\boldsymbol b}^+$ relative to ${\boldsymbol a}^+$

This gives us the kinematic transport theorem, which states that $$\frac{{}^a\text{d}}{\text{d}t}\left({\boldsymbol p}\right) = \frac{{}^b\text{d}}{\text{d}t}\left({\boldsymbol p}\right) + {\boldsymbol \omega}_\text{a/b} \times {\boldsymbol p}.$$ This is an incredibly powerful result, as it let's us easily use as many reference frames as we want while still being able to find how observers in different frames see vectors changing.

Why do we want to do this? The ultimate goal is to find the inertial acceleration. Recall that it is an axiom of Newtonian mechanics that inertial reference frames exist and that Newton's laws hold in inertial reference frames. The KTT makes it almost trivial to use non-inertial frames and still find the inertial velocity and acceleration.

Let's apply the KTT to find the inertial velocity of our ${\boldsymbol p}$ vector. Let's assume we are in the plane spanned by $\hat{\boldsymbol e}_1$, $\hat{\boldsymbol e}_2$ and that $\hat{\boldsymbol \imath}_3 = \hat{\boldsymbol e}_3$. (That is, we are considering planar motion, which we can do since orbits are planar.) Suppose the angular velocity vector is ${\boldsymbol \omega}_\text{e/i} = \dot{\theta}\hat{\boldsymbol e}_3$
$$\begin{align*} {\boldsymbol v} = \dot{\boldsymbol p} &= \dot{r} \hat{\boldsymbol e}_1 + \dot{\theta}\hat{\boldsymbol e}_3 \times r \hat{\boldsymbol e}_1 \\ {\boldsymbol v} &= \dot{r} \hat{\boldsymbol e}_1 + r\dot{\theta}\hat{\boldsymbol e}_2 \end{align*}$$
If we repeat this to find the inertial acceleration, we find that
$${\boldsymbol a} = \dot{\boldsymbol v} =\ddot{\boldsymbol p} = \left(\ddot{r} - r\dot{\theta}^2\right) \hat{\boldsymbol e}_1 + \left(r\ddot{\theta} + 2\dot{r}\dot{\theta}\right)\hat{\boldsymbol e}_2$$(Verifying this expression is left as an exercise for the reader.)

Now, we can apply Newton's second law, ${\boldsymbol F} = m{\boldsymbol a}$. Using Newton's law of gravitation, we find that
$${\boldsymbol F} = -\frac{GMm}{r^2}\hat{\boldsymbol e}_1 = -\frac{\mu m}{r^2}\hat{\boldsymbol e}_1$$where $\mu=GM$ is the gravitational parameter. This results in the equations of motion for an orbiting satellite[1],$$\begin{align*} \ddot{r} - r\dot{\theta}^2 &= -\frac{\mu}{r^2} \\ r\ddot{\theta} + 2\dot{r}\dot{\theta} &= 0 \end{align*}$$

[1] More specifically, this is the position of the satellite's mass center. For a rigid body, you can attach a third reference frame to the body and use various attitude coordinatizations to fully define the pos and attitude of the satellite.

 

[UltrafastPED]

See Frenet–Serret formulas: http://en.wikipedia.org/wiki/Frenet–Serret_formulas

They provide the apparatus for the kinematic properties of any curved path, including orbits.

 

[Philosophaie]

I have had some luck with adding a simple Centripetal Acceleration to the mix:$$a_{centripetal} = \frac{v^2}{r}$$Making the Acceleration equations:


$$a_x = \frac{d||v||}{dt} * uv_x + ||v|| * \frac{d(uv_x)}{dt} + \frac{vx ^ 2 }{ r}$$
$$a_y = \frac{d||v||}{dt} * uv_y + ||v|| * \frac{d(uv_y)}{dt} + \frac{vy ^ 2 }{ r}$$
$$a_z = \frac{d||v||}{dt} * uv_z + ||v|| * \frac{d(uv_z)}{dt} + \frac{vz ^ 2}{ r}$$

It does not sit well with me that the beautiful derivative of time of the magnitude yields velocity then the derivative of velocity does not yield the acceleration and needs an additional centripetal factor for the total acceleration.

 

[D H]

Now, we can apply Newton's second law, ${\boldsymbol F} = m{\boldsymbol a}$. Using Newton's law of gravitation, we find that
$${\boldsymbol F} = -\frac{GMm}{r^2}\hat{\boldsymbol e}_1 = -\frac{\mu m}{r^2}\hat{\boldsymbol e}_1$$where $\mu=GM$ is the gravitational parameter. This results in the equations of motion for an orbiting satellite[1],$$\begin{align*} \ddot{r} - r\dot{\theta}^2 &= -\frac{\mu}{r^2} \\ r\ddot{\theta} + 2\dot{r}\dot{\theta} &= 0 \end{align*}$$

You've ignored the motion of the primary, jhae. Working in a frame with a fixed central body is a non-inertial (accelerating) frame. You need to account for this in your equations of motion. Do that and you'll find that the gravitational parameter μ becomes G(M+m) rather than just GM.

See Frenet–Serret formulas: http://en.wikipedia.org/wiki/Frenet–Serret_formulas

They provide the apparatus for the kinematic properties of any curved path, including orbits.

The problem with those formulae is that the fundamental theorem of space curves has little applicability to orbits. I've found but one paper (and it's a technical report rather than a journal paper) that attempts to use the Frenet-Serret frame to describe orbits, and that one paper has zero citations. That paper: Livieratos, Evangelos. "Intrinsic Parameters and Satellite Orbital Elements." ifp (2003): 295.

Even the plethora of mechanisms used to describe orbits in the two body problem become a lot less useful as soon as one goes beyond a pair of point masses (e.g., the N body problem, non-spherical gravity, atmospheric perturbations, etc.). One has little choice but to revert to numerical integration.

 

[UltrafastPED]

Frenet-Serret is not for finding the orbits; you use it after you have the path. The original question was how to find the unit normal vector in acceleration.

 

[D H]

I have had some luck with adding a simple Centripetal Acceleration to the mix:$$a_{centripetal} = \frac{v^2}{r}$$

If this r is the distance to the primary, that is not the centripetal acceleration except in the case of a circular orbit. For any other type of orbit, the center of curvature is not located at the central body. Gravitation is a central force, not a centripetal force.

Rather than trying to reinvent a 400 year old wheel, I strongly suggest you go to a library. Pick up or buy a copy of "Fundamentals of Astrodynamics and Applications" by David Vallado or "Fundamentals of Astrodynamics" by Bate, Mueller, and White. Vallado is very extensive but can be ridiculously expensive. It costs $150 for a new third edition in paperback. A used version of the older second edition is much more affordable. Bate, Mueller, and White is cheap but it's much less extensive in coverage.

 

[D H]

Frenet-Serret is not for finding the orbits; you use it after you have the path. The original question was how to find the unit normal vector in acceleration.

There is a problem here: Why?

This thread is a prototypical example of an X-Y problem. In an X-Y problem, someone tries to solve problem X, and thinks approach Y would work. The person then starts running into problems with this approach Y and asks for help on this non-viable approach Y rather than on the original problem X.

 

[jhae2.718]

You've ignored the motion of the primary, jhae. Working in a frame with a fixed central body is a non-inertial (accelerating) frame. You need to account for this in your equations of motion. Do that and you'll find that the gravitational parameter μ becomes G(M+m) rather than just GM.

You're correct, of course. For the proper two-body equations of motion I should have constructed a relative position vector ${\boldsymbol r} = {\boldsymbol r}_2 - {\boldsymbol r}_1$ and worked with that. I ignored the motion of the primary object (approximating it as a fixed point in space) and focused on the motion of the satellite in order to simplify things to try to show the method more clearly. Perhaps that was a mistake to do so (and I certainly should have stated those assumptions in my earlier post!).