Magnus effect, what is the correct formula?

Click For Summary

Discussion Overview

The discussion revolves around the Magnus effect and the search for the correct formula to calculate the lift force on a spinning golf ball. Participants explore various formulas, their completeness, and the implications of different variables involved in the calculations, including the lift coefficient and the role of spin.

Discussion Character

  • Exploratory
  • Technical explanation
  • Debate/contested
  • Mathematical reasoning

Main Points Raised

  • One participant expresses doubts about the completeness of several formulas for the Magnus effect, particularly regarding the omission of the lift coefficient in some equations.
  • Another participant clarifies that the vector ω is not a unit vector and discusses the directionality of the cross product (ω x v) in relation to the lift force.
  • Concerns are raised about dimensional analysis of the formula Fm = ½ * CL * ρ * A * v² (ω x v), questioning the dependency of the lift force on the magnitude of ω.
  • Some participants note that the lift coefficient (CL) is often an empirical representation that may depend on the rotation rate of the ball.
  • One participant mentions that the first formula discussed could be a valid approximation, ignoring viscosity, but acknowledges that it does not account for separation effects relevant to dimpled golf balls.
  • A later reply introduces a code repository for simulating lift and drag coefficients based on empirical data, suggesting a method for deriving these coefficients as functions of velocity and spin.

Areas of Agreement / Disagreement

Participants do not reach a consensus on a single correct formula for the Magnus effect. Multiple competing views and interpretations of the formulas remain, particularly regarding the role of the lift coefficient and the dimensional analysis of the equations.

Contextual Notes

Limitations in the discussion include unresolved assumptions about the definitions of variables, the complexity of the physics involved, and the challenges in finding a simple formula that accurately represents all factors influencing the Magnus effect.

Joe28965
Messages
2
Reaction score
0
So I've been investigating the magnus effect for I guess 2 weeks, but I've not been able to find a formula that would seem to be the correct formula. I did find some formula, but I have my doubts about the completeness/correctness of all of these formulas. Note that I am trying to simulate what would happen to a golf ball, so the formulas I found that were for calculating cylinders I have not included.

Fm = π2 * ρ * ω * v * r3

My main issue with this formula is that the lift coefficient is not included in this formula, meaning it doesn't matter if the ball is smooth or (as in the case with golf balls) dimpled.

Fm = S*(ω x v)

I found this formula here http://farside.ph.utexas.edu/teaching/329/lectures/node43.html but it was not noted what S was, only a way to calculate it for the example of baseballs.

Fm = ½ * CL * ρ * A * v2 (ω x v)

I found this formula a lot, but the thing that bugs me with this one is that the cross product of ω and v is not a cross product of the vectors, but of the unit vectors, meaning that the last part is just to calculate the direction of the force. This would not be a problem if that was not the only place ω is used in that formula. Meaning that according to that formula ω is only used to calculate the direction, not the actual size of the force. I find it highly unlikely that the speed of which the ball is rotating does not affect the force, only the fact it is or is not rotating.

I could be wrong of course about these formulas. Or I could be interpreting them wrong. If any of you actually know something about the magnus force that could help me, I'd love for you to reply to this topic and maybe help me out, since I'm stuck at the moment.
 
Physics news on Phys.org
ω is not a unit vector. Assume the ball is moving left to right with back spin, ω is a vector pointed towards the observer, v is a vector pointed to the right, and (ω x v) is a vector pointed up with a magnitude of |ω| |v|.
 
but then the units wouldn't match up again. Right now it's:

Kg*m/s2 = kg/m3 * m2 * m2/s2 = kg*m4/m3*s2 = Kg*m/s2

If you would add ω and v as normal vectors, not unit vectors.
 
So that should have been ##(\hat{\omega} \ x \ \hat{v})## ?

Apparently the spin affects coefficient of lift:

http://spiff.rit.edu/richmond/baseball/traj/traj.html

That article was about baseballs. The spin effect is more evident in the case of ping pong balls.
 
Last edited:
As with all topics in aerodynamics, things are typically complicated here. For example, lift (and drag) coefficients are typically an empirical (or model) representation of a much more complicated phenomenon. In this case, the rotation rate would be wrapped up into the lift coefficient. On the other hand, as with the first equation, the lift coefficient need not be used. The first equation represents the lift on a ball due to the Magnus effect ignoring viscosity (if my memory serves me) so it is perfectly valid as an approximation. It won't account for separation effects, and therefore for golf ball dimples, though.

If you want to include all the physics then you are just out of luck in finding a simple formula because one simply doesn't exist. You would have to flat out solve the Navier-Stokes equations, and any simplifications you make from there will be removing at the very least a little bit of the physics for the sake of practicality.
 
mm6c5w.jpg
23sviwx.png
mm6c5w.jpg
23sviwx.png
Joe28965 said:
but then the units wouldn't match up again. Right now it's:

Kg*m/s2 = kg/m3 * m2 * m2/s2 = kg*m4/m3*s2 = Kg*m/s2

If you would add ω and v as normal vectors, not unit vectors.

Well folks, I was going to open a new thread to ask this but I found this one.

I ask here because the question has not been properly answered, that is what I think.

The dimensional analysis of the equation is what bothers me.

Fm = ½ * CL * ρ * A * v2 (ω x v)

The term Fm= ½ * Cl * ρ * A * v2 is already dimensionally correct. We have:
Force= (Mass * Lenght)/Time2 = Density (Mass/Length3) * Area (Lenght2)* Speed(Lenght2/Time2) = Mass * Lenght / Time2

And this result raises the question on this term:

## \vec w \times \vec v ##

If it is dimensionless we are saying the lift force is not dependent on w modulus which is absurd.

If ##\vec w## and ##\vec v## have their dimensions we get:

## \vec w \times \vec v ## = IwI * IvI* sin (90) * unit vector.

Since w= (1/time) and v = lenght/time I think w * v = Lenght/Time 2

Force = (Mass * Lenght)/ (Time2) * Lenght/Time 2 = Mass Lenght2 / Time4

Which is not correct.

So what am I doing wrong?

In this article:
http://spiff.rit.edu/richmond/baseball/traj/traj.html#spin
The force is given by:
mm6c5w.jpg


I found the same formula in another article:
https://www.researchgate.net/publication/328816286_Model_trajectories_for_a_spinning_tennis_ball_I_The_service_stroke
Guess what, now it is different:

23sviwx.png


Can anybody clarify what is the correct dimensional analysis of this equation?

Thanks.
 

Attachments

  • mm6c5w.jpg
    mm6c5w.jpg
    6.5 KB · Views: 3,254
  • 23sviwx.png
    23sviwx.png
    3.6 KB · Views: 3,661
Last edited by a moderator:
I found on another book FL does not depend explicitly on w, it does just because CL is a function of w.

I hope this helps others.
 
Joe28965 said:
So I've been investigating the magnus effect for I guess 2 weeks, but I've not been able to find a formula that would seem to be the correct formula. I did find some formula, but I have my doubts about the completeness/correctness of all of these formulas. Note that I am trying to simulate what would happen to a golf ball, so the formulas I found that were for calculating cylinders I have not included.

Fm = π2 * ρ * ω * v * r3

My main issue with this formula is that the lift coefficient is not included in this formula, meaning it doesn't matter if the ball is smooth or (as in the case with golf balls) dimpled.

Fm = S*(ω x v)

I found this formula here http://farside.ph.utexas.edu/teaching/329/lectures/node43.html but it was not noted what S was, only a way to calculate it for the example of baseballs.

Fm = ½ * CL * ρ * A * v2 (ω x v)

I found this formula a lot, but the thing that bugs me with this one is that the cross product of ω and v is not a cross product of the vectors, but of the unit vectors, meaning that the last part is just to calculate the direction of the force. This would not be a problem if that was not the only place ω is used in that formula. Meaning that according to that formula ω is only used to calculate the direction, not the actual size of the force. I find it highly unlikely that the speed of which the ball is rotating does not affect the force, only the fact it is or is not rotating.

I could be wrong of course about these formulas. Or I could be interpreting them wrong. If any of you actually know something about the magnus force that could help me, I'd love for you to reply to this topic and maybe help me out, since I'm stuck at the moment.
I just recently ran into the same problem as you, and I've been testing the different formulas I've been finding on a shot simulator I'm developing, (for fun purposes) and the formula I found was working. was this one.

Fm = ½ * CL * ρ * A * v2 (ω x v)

where ω and v are the unit vectors.

The think here is as somebody else mentioned CL is in function of the spin magnitude.

What I did to find that function I found a code repository that has lift and drag coefficients (constant) for each of the different clubs a golfer can use, https://github.com/x00126177/GolfBall-Simulator
I then went and looked up the average ball speeds for those clubs for Tour players https://blog.trackmangolf.com/trackman-average-tour-stats/

made a plot in spreadsheets with the data and used a polynomial regression to fit the drag, and lift coefficients to functions of velocity and spin respectively.

here are the functions I found for the drag and lift coefficients from the data I gathered.
CD(v) = 0.0606 + 0.0129v - 2.57E-04v2 + 1.49E-06v3

CD(ω) = 0.176 - 4x10-4ω + 1.36x10-6ω2 - 1.58x10-9ω3 + 6.41x10-13ω4

The results I get are close to what Trackman calculates, but are a bit longer. my guess at the moment is that I'm missing the wind force calculation adjustment.
 

Similar threads

  • · Replies 9 ·
Replies
9
Views
4K
  • · Replies 25 ·
Replies
25
Views
4K
  • · Replies 1 ·
Replies
1
Views
3K
  • · Replies 2 ·
Replies
2
Views
3K
  • · Replies 3 ·
Replies
3
Views
3K
  • · Replies 15 ·
Replies
15
Views
3K
  • · Replies 4 ·
Replies
4
Views
4K
  • · Replies 7 ·
Replies
7
Views
3K
  • · Replies 7 ·
Replies
7
Views
4K
  • · Replies 1 ·
Replies
1
Views
3K