How to change CSpice initial velocity from an Asteroid SPK file?

Click For Summary
SUMMARY

This discussion focuses on modifying the initial velocity of an asteroid using the CSpice toolkit, specifically for asteroid 7482 (1994 PC1). The user successfully integrated CSpice with a modified C++ implementation of Newtonian mechanics, utilizing the four-order symplectic Candy-Rozmus integration method. Key tools mentioned include the de440.bsp and 2007482(1994 PC1).bsp SPK files, and the importance of unit conversion and frame specifications (ECLIPJ2000) is emphasized. The user recommends having APK Editor as an additional tool for this process.

PREREQUISITES
  • Familiarity with CSpice for planetary and asteroid position calculations
  • Understanding of Newtonian mechanics and symplectic integration methods
  • Experience with C++ programming for modifying existing code
  • Knowledge of SPK files, specifically de440.bsp and 2007482(1994 PC1).bsp
NEXT STEPS
  • Explore the CSpice toolkit documentation for advanced features
  • Study the implementation of the Candy-Rozmus integration method in C++
  • Learn about unit conversion in astrodynamics, particularly for ECLIPJ2000
  • Investigate the use of APK Editor for modifying software tools
USEFUL FOR

Astronomers, astrophysicists, and software developers working on asteroid trajectory calculations and simulations will benefit from this discussion.

forofo
Messages
4
Reaction score
1
TL;DR
Change spice spk file velocity to determinate the position of an asteroid
Hello, I've made a SPK file for asteroid 7482 (1994 PC1) with Horizon. I wan't to change the initial velocity with cspice, because I want to know where it will be in a future time at the speed changed. (now I'm using Newtons calc but is slow an error increases with time). This is for calculating NEO's asteroid position desviation.
 
Astronomy news on Phys.org
I found the answer: In this link https://github.com/fding/nbody (Newtonian mechanics and a four-order symplectic Candy-Rozmus integration), there is c++ code and phyton. I've modify the c++, so now I calc the inital positions and speeds with cspice and pass then to Newton candy with other asteroid speed. Then I run both to the future to see the asteroid distance to any planet. Just take care about units conversion, cspice frame(ECLIPJ2000) and obverser naif Solar System(0), Used de440.bsp and 2007482(1994 PC1).bsp
 
Reply
  • Informative
Likes   Reactions: berkeman
You might want have APK Editor in your tool kit.
 

Similar threads

Time required to change the velocity of an asteroid?
  • · Replies 23 ·
Replies
23
Views
3K
Undergrad Change in Momentum from mass contribution to the moving object
  • · Replies 61 ·
3
Replies
61
Views
4K
Undergrad How Can I Improve My 2D Newtonian Gravity Simulator?
  • · Replies 7 ·
Replies
7
Views
5K
Finding initial velocity of a projectile. Help with force and acceleration
  • · Replies 15 ·
Replies
15
Views
11K
How to Apply Friction and Gravity to a 3D Velocity Vector?
  • · Replies 32 ·
2
Replies
32
Views
9K
Undergrad Damped oscillator with changing mass
  • · Replies 131 ·
5
Replies
131
Views
8K
Undergrad Modeling a Bottle Rocket: Analyzing the Maximum Altitude
  • · Replies 4 ·
Replies
4
Views
5K
Undergrad Does G effectively change in an expanding metric?
  • · Replies 15 ·
Replies
15
Views
2K
Atmospheric Heating from Meteorite Impact
  • · Replies 2 ·
Replies
2
Views
3K
Graduate How is the observed redshift caused by a recessional velocity?
  • · Replies 9 ·
Replies
9
Views
5K