Skip to content

Commit de12cac

Browse files
Fix SRP page
1 parent 3390f2a commit de12cac

File tree

1 file changed

+7
-8
lines changed
  • docs/nyxspace/MathSpec/models

1 file changed

+7
-8
lines changed

docs/nyxspace/MathSpec/models/srp.md

Lines changed: 7 additions & 8 deletions
Original file line numberDiff line numberDiff line change
@@ -1,5 +1,3 @@
1-
# Solar radiation pressure
2-
31
SRP is implemented as a `ForceModel`, meaning that the implementation computes a force from which the instantaneous mass of the spacecraft is subsequently divided to apply an acceleration. Therefore, as opposed to most formulations in the literature, we will use a force formulation.
42

53
Currently, Nyx only supports spherical SRP model.
@@ -20,18 +18,19 @@ Note that $C_r$ and $\mathcal{A}$ are stored in the "context" passed to the EOM
2018

2119
First, we compute the position of the Sun as seen from the spacecraft, and its unit vector, respectively $\mathbf{r_\odot}$ and $\mathbf{\hat r_\odot}$. Then, we compute the shadowing factor, $k$, using the [eclipse model](../celestial/eclipse.md).
2220

23-
Compute the norm of Sun vector in AU, $ ||\mathbf{r_\odot}||_{\text{AU}}$ by dividing the $\mathbf{r_\odot} $ vector by 1 AU.
21+
Compute the norm of Sun vector in AU, $||\mathbf{r_\odot}||$ by dividing the $\mathbf{r_\odot}$ vector by 1 AU.
2422

25-
Compute the flux pressure as follows:
23+
Compute the flux pressure ($\Phi$) as follows:
2624

27-
$$\Phi_{\text{SRP}} = \frac{k\phi}{c} \left(\frac{1.0}{||\mathbf{r_\odot}||_{\text{AU}}} \right)^2$$
25+
$$\Phi = \frac{k\phi}{c} \left(\frac{1.0}{||\mathbf{r_\odot}||} \right)^2$$
2826

2927
Finally, return the SRP force [^2]:
3028

31-
$$ \mathbf{F}_{\text{SRP}} = C_r \mathcal{A} \Phi_{\text{SRP}} \mathbf{\hat r_\odot}$$
29+
$$ \mathbf{F} = C_r \mathcal{A} \Phi \mathbf{\hat r_\odot}$$
30+
31+
## Estimation
3232

33-
!!! note
34-
Although the above derivation mentions the Sun, Nyx trivially supports any other light source regardless of the integration frame.
33+
Nyx can estimate the coefficient of reflectivity $C_r$ in orbit determination scenarios, as shown in [the Lunar Reconnaissance Orbiter](../../showcase/04_lro_od/index.md) example.
3534

3635
??? check "Validation"
3736
Nyx has three validation scenarios for the SRP computation to ensure that we test full illumination (`srp_earth_full_vis`), long penumbra passages (`srp_earth_meo_ecc_inc`), and very short penumbra passages (`srp_earth_penumbra`). In all of the test cases, we propagate a spacecraft for 24 **days** to ensure that a high amount of error can accumulate if the modeling is incorrect. The worst absolute position error is _high_ compared to GMAT: 287 meters.

0 commit comments

Comments
 (0)