NKTg Law: Simulating Neptune's Orbit Using NASA Data

### 🧠 Overview

This project demonstrates how to use the **NKTg Law**, a physics-based model that incorporates mass variation, to simulate and predict the motion of Neptune using only NASA's 2023 data and Python.

By modeling **atmospheric gas loss**, we successfully predicted Neptune’s **2024 orbital parameters** — and achieved nearly **zero error** in position and velocity compared to actual NASA data.

The project shows how simple math and real data can be combined using Python to simulate real-world planetary motion — no advanced numerical solvers required.

---

### 🎯 Goal

• Apply the NKTg Law to Neptune
• Use only 2023 NASA data
• Simulate Neptune’s 2024 orbit
• Compare predicted results with NASA’s published 2024 data

💡 Inspiration

markdown

Copy

Edit

• NASA publishes reliable orbital and planetary data.
• Traditional motion models assume constant mass.
• The **NKTg Law** allows us to simulate motion even if the body is **losing mass** (e.g., gas giants).
• Inspired by applying physical modeling to real data using open tools like **Python and Linux**.

🔬 How It Works

markdown

Copy

Edit

The NKTg Law defines motion through two conserved quantities:

• `NKTg₁ = x × p` (Position × Momentum)
• `NKTg₂ = (dm/dt) × p` (Mass change × Momentum)

Where:

• `x`: position
• `v`: velocity
• `m`: mass
• `p = m × v`
• `dm/dt`: mass loss rate

If `NKTg₁` and `NKTg₂` are constant, and you know the new `m`, you can calculate `x` and `v` **in reverse** — ideal for modeling time-forward orbital behavior.

This simulation assumes:

• Constant `v = 5.43 km/s`
• Mass loss: `–0.00002000 kg/s`

Then we solve:

```python

p = m × v

x = NKTg₁ / p

yaml

Copy

Edit

---

## 📊 Data Sources

```markdown

• [NASA JPL Horizons](https://ssd.jpl.nasa.gov/horizons) – Position and velocity of Neptune
• [NASA Planetary Fact Sheet](https://nssdc.gsfc.nasa.gov/planetary/factsheet/neptunefact.html) – Mass of Neptune
• [Nature: Hydrogen Escape](https://www.nature.com/articles/35036049) – Assumed gas loss rate

🐍 Code Snippet (Python Prototype)

# Given:

NKTg1 = 2.503e36 # from 2023 data

v = 5.43

mass_2024 = 1.02429900e26

# Calculate momentum and position

p = mass_2024 * v

x = NKTg1 / p

print("Predicted position (km):", x)

✅ Results

Date x Error (km) v Error (km/s) Mass Error (%)

2024-01-01 0 0 ~0.000020%

2024-07-01 0 0 ~0.000020%

Position and velocity match NASA exactly.

Minor mass deviation is consistent with gas loss model.

💻 Tools Used

Python 3.10

NumPy + Pandas for basic math and data

LibreOffice Calc to view and format result tables

Linux Mint 21 for the entire environment

🚀 What's Next?

Build a Qt GUI application to allow anyone to simulate motion

Extend simulation to Uranus, Jupiter, or even comets

Visualize orbits using Matplotlib or 3D rendering

🤝 Collaboration

Interested in:

Orbital mechanics?

Qt or PyQt development?

Planetary simulations?

Feel free to connect or suggest improvements. The model is lightweight and easy to port to embedded or scientific platforms.