Gravitational Sphere of Influence | Orbital Mechanics with Python

The two-body problem involving particles with masses M1 and M2 can be reduced to a one-body problem using reduced mass μ = m1m2/(m1 + m2). The Lagrangian equation L = T - V helps derive this reduction. Kepler's Laws describe planetary motion: First Law states planets move in elliptical orbits with the Sun at one focus; Second Law (Law of Areas) states equal areas are swept in equal times, meaning planets move faster when closer to the Sun; Third Law states T² ∝ r³, meaning orbital period squared is proportional to orbital radius cubed.

Kepler's three laws of planetary motion describe that planets move in elliptical orbits with the Sun at one focus (First Law), sweep equal areas in equal times (Second Law, proven through angular momentum conservation where dA/dt = L/(2m) is constant), and have orbital periods related to their semi-major axes by T² ∝ a³ (Third Law). For the two-body problem, the system can be reduced to a single body of reduced mass μ = M₁M₂/(M₁+M₂) moving under central force, with the center of mass remaining stationary; when one mass dominates (like the Sun over Earth), the reduced mass approximates the smaller mass and the center of mass approaches the larger mass, justifying treating the larger body as stationary.

Kepler's First Law states that planets move in elliptical orbits with the Sun at one focus. The derivation begins by recognizing that planetary motion under gravitational attraction follows the inverse square law. The two-body problem (Sun-planet interaction) is reduced to a one-body problem using reduced mass μ = (m₁m₂)/(m₁+m₂). For the Sun-planet system, this reduced mass replaces the actual planetary mass in equations of motion. The general orbit equation for inverse square law motion is 1/r = (μc/J²)[1 + e cosθ], where μ is reduced mass, c = GMsMp, J is angular momentum, and e is eccentricity.

Kepler's three laws describe planetary motion: (1) Planets move in elliptical orbits with the Sun at one focus; (2) The radius vector sweeps equal areas in equal times; (3) The square of the orbital period is proportional to the cube of the semi-major axis. From Newton's laws, the two-body problem reduces to motion under central force with conserved energy and angular momentum. The trajectories are conic sections: ellipses (bound, E<0), parabolas (E=0), and hyperbolas (E>0). The first cosmic velocity is v₁ = √(GM/r) for circular orbit, and the second cosmic velocity is v₂ = √(2GM/r) for escape.

The two-body problem describes how two celestial bodies move under mutual gravitational attraction. In the solar system, the Sun dominates because it is vastly more massive than all other bodies. Johannes Kepler discovered that planets move in elliptical orbits with the Sun at one focus, not circular orbits as previously believed. Isaac Newton formulated the law of universal gravitation, stating that every particle attracts every other particle with a force proportional to the product of their masses and inversely proportional to the square of the distance. This force causes bodies to move along conic sections—ellipses, parabolas, or hyperbolas—depending on their energy. Kepler's Second Law states that a line joining a planet and the Sun sweeps equal areas in equal times, reflecting conservation of angular momentum. Kepler's Third Law establishes that orbital period squared is proportional to average distance cubed. In two-body systems, both bodies actually move around their common center of mass, but by fixing coordinates to one body, we can observe relative motion as if it were moving around a stationary point.

Matplotlib is a Python library for creating various types of graphs, while NumPy provides mathematical calculation tools; to generate a graph, import both libraries using 'import matplotlib.pyplot as plt' and 'import numpy as np', create data arrays using functions like 'np.linspace()' for time values and 'np.cos()' for function values, then plot using 'plt.plot(t, y)' and display with 'plt.show()'; increasing the number of data points results in smoother, more fluid graphs, and graphs can be customized with titles and axis labels using 'plt.title()', 'plt.xlabel()', and 'plt.ylabel()'.

Two essential libraries extend Python's capabilities: (1) NumPy provides efficient numerical computing with arrays, offering functions like np.array() for creation, np.mean() for averages, np.sum() for totals, np.max() for maximum, and np.std() for standard deviation; (2) Matplotlib enables data visualization with functions like plt.plot() for line charts, plt.title() for titles, and plt.xlabel()/plt.ylabel() for axis labels. Together, these libraries enable data analysts to perform numerical calculations and create visualizations that reveal patterns and insights in data, completing the data analysis workflow from raw data to actionable insights.

Python's scientific computing ecosystem relies on importing packages like NumPy and Matplotlib. NumPy provides efficient array operations for numerical computing, with functions like 'np.linspace()' creating evenly spaced arrays. Matplotlib enables data visualization through the 'plt.plot()' function. Packages can be imported with nicknames (e.g., 'import numpy as np') for convenience. Standard nicknames like 'np' for NumPy and 'plt' for matplotlib.pyplot are widely recognized. These packages work together to create numerical arrays and generate plots, forming the foundation for scientific computing in Python.

This section covers essential Python libraries for numerical computing and data analysis. NumPy is imported using: import numpy as np. Arrays are created using np.array() and support element-wise operations. The shape attribute returns dimensions. Boolean indexing enables filtering: array[array > 5]. Element-wise operations between arrays require compatible shapes. The linspace function generates evenly spaced values. Zero matrices are created using np.zeros((rows, columns)). Matrix operations like multiplication (np.dot()) and inversion (np.linalg.inv()) are supported. Pandas is imported using: import pandas as pd. DataFrames are created from dictionaries. The info() method displays column names, data types, and non-null counts. The describe() method shows statistical summaries. Missing data appears as NaN values and can be handled using dropna() or fillna(). Duplicate rows are detected using duplicated() and removed using drop_duplicates(). Correlation between columns is calculated using corr(), returning values from -1 to 1. Scatter plots use plot(kind='scatter', x='x_col', y='y_col'). Histograms use plot(kind='hist', column='col'). Matplotlib is imported using: import matplotlib.pyplot as plt. Basic plots use plt.plot(x, y). Color customization uses the color parameter. Markers can be customized. Axis labels use plt.xlabel() and plt.ylabel(). Titles use plt.title(). Legends use plt.legend(). Multiple plots can be combined. The show() function displays plots. Interactive mode (plt.ion()) should be disabled during evaluations.

Matplotlib is a Python library for creating static, animated, and interactive visualizations. To create basic plots, import matplotlib.pyplot as plt, then use plt.plot() for line plots, plt.scatter() for scatter plots, and plt.hist() for histograms. Customize plots by adding titles with plt.title(), axis labels with plt.xlabel() and plt.ylabel(), and legends with plt.legend(). Display the plot using plt.show().

The Cartesian coordinate system uses three orthonormal unit vectors (ex, ey, ez) that are fixed in space. For a point M with coordinates (x, y, z), the position vector is OM = x·ex + y·ey + z·ez. The displacement vector is dOM = dx·ex + dy·ey + dz·ez. Velocity is the first derivative: v = ẋ·ex + ẏ·ey + ż·ez. Acceleration is the second derivative: a = ẍ·ex + ÿ·ey + z̈·ez. Since the basis vectors are constant, only the coordinate components are differentiated. This system is simplest for linear motion without rotation.

To study motion, coordinate systems with defined axes are essential. For 3D space, use x, y, z axes; for 2D planes, use x and y; for 1D lines, use only x. The position vector is expressed as r = xi + yj + zk, where components are found by projecting onto axes. Vectors are added by combining corresponding components. The magnitude (distance from origin) is calculated using |r| = √(x² + y² + z²), derived from the Pythagorean theorem. Motion complexity depends on dimensionality: ant crawling on ground needs 2D (x,y), car on straight road needs 1D (x).

This lesson covers position vectors in 2D and 3D space. A position vector has its tail at the origin (0,0) or (0,0,0) and its head at a point (x,y) or (x,y,z). Vectors can be expressed as column vectors (written top to bottom), row vectors (written left to right), or in basis vector form using unit vectors i, j, k. Components are separated by spaces, not commas. To find coordinates in 3D, draw perpendiculars from the point to each axis. To draw a position vector, first find the coordinates from the vector expression, then plot the point and draw a line from the origin to that point. For negative components, move in the negative direction along each axis.

In Cartesian coordinates, any point is represented by position vector R = xî + yĵ + zk̂, where unit vectors are constant. Velocity is V = dR/dt = ẋî + ẏĵ + żk̂, and velocity squared is V² = ẋ² + ẏ² + ż². The position differential is dR = dxî + dyĵ + dzk̂. This simple form makes Cartesian coordinates ideal for basic mechanics problems, though more complex systems require generalized coordinates.

This section covers the fundamental mathematical framework for describing motion in two and three dimensions. The 2D Cartesian coordinate system uses X and Y axes with an origin, where any point is represented by coordinates (x, y). Velocity is the first derivative of position: v = (dx/dt)i + (dy/dt)j. Acceleration is the second derivative of position: a = (d²x/dt²)i + (d²y/dt²)j. In 3D motion, the system extends to include the Z-axis with coordinates (x, y, z), and velocity and acceleration gain corresponding Z-components. These vector representations allow precise mathematical description of particle motion in space.