search

groops-devs / groops

A software toolkit for gravity field recovery and GNSS processing
GNU General Public License v3.0
185 stars 107 forks source link

Non-gravitational forces from ACT1B when simulating an orbit #96

Closed ReubenLi1999 closed 2 years ago

ReubenLi1999 commented 2 years ago

Hello developers,

I realised that you guys have done a great job of simulating a pure dynamic orbit from a given force function which includes non-gravitational forces.

It seems that the non-gravitational forces used here are calculated based on models, instead of measurements by the accelerometers. Have you guys considered calculating the non-conservative forces by the measurements of accelerometers?

It would be great if you add this feature.

Best regard,

Reuben

akvas commented 2 years ago

Hi @ReubenLi1999,

thanks for this enhancement proposal! What would be your use case?

If you already have an approximate orbit, GROOPS can integrate dynamic orbits with accelerometer data. The program you are looking for is PreprocessingVariationalEquation. You can find a detailed description of how to use it in our GRACE gravity field recovery cookbook. Would this be applicable to your use case?

Cheers, Andreas

ReubenLi1999 commented 2 years ago

Hi @akvas,

Thank you for the helpful guidance. I am sorry that I did not interpret my use case clearly. I want to propagate dynamic orbits of GRACE-FO from a given force model (static gravity and tides included) and ACT1B data, just like what you have done in SimulateOrbit. The only difference is that ACT1B is required in my computation. Therefore, I am sorry that the content of GRACE gravity filed recovery cookbook does not seem to be what I want.

Additionally, as you mentioned in PreprocessingVariationalEquation, non-conservative forces like solar radiation pressure need the orientation of the satellite (inputfileStarCamera). However, I did not find any interface to input the orientation of the satellite (SCA) when SimulateOrbit is calculating the solar radiation pressure. Would you mind enlightening me on how to use the orientation of the satellite to calculate the non-gravitational forces?

Thanks for your time.

Best regard,

Reuben

akvas commented 2 years ago

Hi @ReubenLi1999,

I put together a minimal example which integrates an orbit for GRACE-C over a day. Note that you will probably need to adapt the file paths for the L1B files and gravity field model. Let me know if that solves your problem.

Cheers, Andreas