CL_fo_srpAcc

Description

• Acceleration due to solar radiation pressure (SRP).

  The acceleration is directed from the Sun to the considered position and is inversely proportional to the distance squared.

  The SRP coefficient (coefp) is defined by: coefp = cp * area / mass (with cp: reflectivity coefficient, between 1 and 2).

  Eclipses can be taken into account or not. If they are, the acceleration is multiplied by a factor equal to 1 out of the eclipse region and less than 1 otherwise.

• Notes:

  - The coordinates frame can be any frame. The origin of the frame does not matter, except if eclipses are taken into account, in which case the origin of the frame must be the center of the eclipsing body (usually the same as the central body).

  - The calculation of the eclipse uses the radius of the eclipsing body (er), and the radius of the Sun (ersun). If ersun is empty ([]) then an internal value is used.

  - p0 is the solar radius pressure at 1 AU. If p0 is empty ([]) then the internal value is used (equal to the total solar irradiance divided by the speed of light).

• See Force models for more details.

Parameters

pos:

Position vector (from any position) [m]. (3xN or 3x1)

pos_sun:

Sun position (same origin as pos) [m]. (3xN or 3x1)

coefp:

SRP coefficient (cp*area/mass) [m^2/kg]. (1xN or 1x1)

ecl:

(optional, boolean) %t if eclipses are taken into account; %f otherwise. Default is %t. (1x1)

er:

(optional) Equatorial radius of eclipsing body. Default is %CL_eqRad. [m] (1x1)

ersun:

(optional) Equatorial radius of the Sun. Default is [] (internal value is used). [m] (1x1)

p0:

(optional) Solar radiation pressure at 1 AU. Default is [] (internal value is used). [N/m^2] (1x1)

acc:

Acceleration [m/s^2]. (3xN)

Authors

• CNES - DCT/SB

See also

• CL_fo_srpPot

• CL_gm_eclipseCheck

Examples

pos = [7000.e3; 0; 0]; // ECI pos_sun = CL_eph_sun(CL_dat_cal2cjd(2000,3,21)); // ECI coefp = 1.5 * 10 / 1000; CL_fo_srpAcc(pos, pos_sun, coefp, %t)