moon_shine

ssapy.compute.moon_shine(r_moon, r_sat, r_earth, r_sun, radius, albedo, albedo_moon, albedo_back, albedo_front, area_panels)[source][source]

Calculate the fractional flux of sunlight reflected from the Moon to the satellite.

This function computes the flux of sunlight reflected from the Moon to the satellite, including contributions from both the front and back surfaces of the satellite’s solar panels.

Parameters:

r_moon(n, 3) numpy.ndarray: Array of coordinates representing the position of the Moon.
r_sat(n, 3) numpy.ndarray: Array of coordinates representing the position of the satellite.
r_earth(n, 3) numpy.ndarray: Array of coordinates representing the position of the Earth.
r_sun(n, 3) numpy.ndarray: Array of coordinates representing the position of the Sun.
radiusfloat: Radius of the satellite in meters.
albedofloat: Albedo of the satellite’s surface.
albedo_moonfloat: Albedo of the Moon.
albedo_backfloat: Albedo of the back surface of the satellite’s solar panels.
albedo_frontfloat: Albedo of the front surface of the satellite’s solar panels.
area_panelsfloat: Area of the satellite’s solar panels in square meters.

Returns:

dict: Dictionary containing the flux contributions from the Moon to the satellite: - ‘moon_bus’: Fraction of light reflected off the satellite’s bus from the Moon. - ‘moon_panels’: Fraction of light reflected off the satellite’s panels from the Moon.

Notes:

The function assumes that the solar panels are always facing the Sun and calculates flux based on the phase angles.
Flux contributions from both the front and back surfaces of the solar panels are computed.

Example:

>>> r_moon = np.array([[1e8, 1e8, 1e8]])
>>> r_sat = np.array([[1e7, 1e7, 1e7]])
>>> r_earth = np.array([[0, 0, 0]])
>>> r_sun = np.array([[1e11, 0, 0]])
>>> moon_shine(r_moon, r_sat, r_earth, r_sun, radius=0.4, albedo=0.20, albedo_moon=0.12, albedo_back=0.50, albedo_front=0.05, area_panels=100)
{'moon_bus': array([...]), 'moon_panels': array([...])}