Best ODE algorithm to use for time/velocity independent potentials?

  • Thread starter Signifier
  • Start date
  • #1
76
0

Main Question or Discussion Point

Hello, I am not sure if this is the right place to ask this, but I don't readily see a "numerical methods" forum here so I assumed this would be the place to go. Sorry if I overlooked another place to post this!

Anyway, I have some points interacting via potentials that are dependent only on distances (IE, not time or velocity), and I am wondering what would be the best ODE algorithm to use to simulate the motions of these points. The number of interacting points is fairly large, from 150 to 300. Let's say they're interacting via an additive pairwise potential similar to Newton's law of gravitation.

What would be a decent ODE algorithm to use that is both fairly accurate (and not subject to huge error propagation) and not too complicated for solving the equations of motion?

Any help would be greatly appreciated.
 

Answers and Replies

  • #2
Dr Transport
Science Advisor
Gold Member
2,330
446
Runge-Kutta is a standard ODE algorithm for problems like yours. Try a fourth order for reasonably decent convergence.
 
  • #3
76
0
Is RK4 noticeably better than, for example, Beeman or velocity-Verlet? I've never used RK4, though I know it's the word for many people. I don't know much about it. Right now, I'm using Beeman for energy conservation, though what I'm most interested in is accuracy of trajectories. Would RK give me a bump up in this category?

Thank you in advance!
 

Related Threads on Best ODE algorithm to use for time/velocity independent potentials?

Replies
18
Views
3K
  • Last Post
Replies
7
Views
2K
  • Last Post
Replies
2
Views
14K
Replies
25
Views
2K
Replies
4
Views
205
Replies
2
Views
1K
Replies
13
Views
6K
Replies
4
Views
1K
Top