ENH: Generic Surfaces and Generic Linear Surfaces

[MateusStano]

Pull request type

• [x] Code changes (bugfix, features)

Checklist

• [ ] Tests for the changes have been added (if needed)
• [ ] Docs have been reviewed and added / updated
• [ ] Lint (black rocketpy/ tests/) has passed locally
• [x] All tests (pytest tests -m slow --runslow) have passed locally
• [ ] CHANGELOG.md has been updated (if relevant)

Description

• Added the GenericSurface class that takes in coefficients defined by either a function or a .csv

• Added the LinearGenericSurface class that takes in coefficients derivatives and calculates the total coefficients assuming linearity

• Changed all add_surfaces position arguments to be 3 dimensional

• Improved Flight class structure to be indifferent to the type of AeroSurface

• There is an example simulation in coeff_testing.ipynb of calisto using coefficients that were extracted directly from rocketpy's standard model (Barrowman's). It uses the .csvs that have been committed. These .csvs have a LOT of points, this is due to the shepard interpolation in the Function class, which makes the interpolated values used in the simulation extremely noisy. The files are this big for validation purposes. We should reduce them and create tests with them

Breaking change

• [ ] Yes
• [x] No

---

Remaining Tasks

• [x] 1. Add a representation of GenericSurface in the Drawing.
• [x] 2. Create Plots and Prints classes for both GenericSurface and LinearGenericSurface
• [x] 3. Add sideslip and partial angle of attack to flight plots and/or prints
• [ ] 4. Create a method in the Rocket class to define the coefficients of the entire rocket. The method should receive a GenericSurface or LinearGenericSurface and use its coefficients for simulation, overwriting any other added aerodynamic surface. The aerodynamic surfaces should still be used for the drawing though. The received generic surface should be positioned in the center of the dry mass of the rocket for the desired behavior
• [ ] 5. Currently, we receive cL, cD, cQ, however, user's might prefer to input the coefficients relating to the body axes directly (cX, cY, cZ). This should be an optional parameter in the instantiation of the class.
• [x] 6. Pylint
• [x] 7. Tests
• [x] 8. Doc pages

It would be truly great to get some help on these minor tasks

Future Tasks

1. Create an axisymmetric surface class that needs only cD, cL, either center of pressure or cm, and cl (roll moment coefficient)
2. The Rocket class should have total coefficient attributes: cL, cQ, cD, cm, cn, cl. They should be included in Plots or Prints
3. Add lateral stability analysis. The current center of pressure, static, and stability margins only works for axisymmetric configurations. GenericSurfaces and LinearGenericSurfaces are currently not included in those analysis.

---

Contributors

Special thanks to @kevin-alcaniz and Faraday Rocketry UPV for their contributions!!

[codecov[bot]]

Codecov Report

Attention: Patch coverage is 85.51587% with 73 lines in your changes missing coverage. Please review.

Project coverage is 76.02%. Comparing base (bae7cae) to head (b9b44dd). Report is 101 commits behind head on develop.

Files with missing lines	Patch %	Lines
rocketpy/simulation/flight.py	66.19%	24 Missing :warning:
rocketpy/prints/aero_surface_prints.py	33.33%	14 Missing :warning:
rocketpy/plots/rocket_plots.py	45.00%	11 Missing :warning:
rocketpy/rocket/rocket.py	76.19%	10 Missing :warning:
...ocketpy/rocket/aero_surface/fins/free_form_fins.py	94.61%	7 Missing :warning:
rocketpy/rocket/aero_surface/generic_surface.py	95.60%	4 Missing :warning:
rocketpy/plots/aero_surface_plots.py	92.00%	2 Missing :warning:
rocketpy/prints/rocket_prints.py	75.00%	1 Missing :warning:

Additional details and impacted files

```diff @@ Coverage Diff @@ ## develop #680 +/- ## =========================================== + Coverage 75.45% 76.02% +0.56% =========================================== Files 96 99 +3 Lines 10887 11278 +391 =========================================== + Hits 8215 8574 +359 - Misses 2672 2704 +32 ```

:umbrella: View full report in Codecov by Sentry.
:loudspeaker: Have feedback on the report? Share it here.

[Lucas-Prates]

I get the following error from Scipy when trying to load the coefficients using a csv file with insufficient data. Maybe it's worth checking the number of rows on the csv?