esa / LADDS

Large-scale Deterministic Debris Simulation - Codebase for the ARIADNA Study between TU Munich and ESA's Advanced Concepts Team.
GNU General Public License v3.0
5 stars 3 forks source link

Run new conjunction determination with dt = 1s and dt = 5s #67

Closed gomezzz closed 2 years ago

gomezzz commented 2 years ago

Feature

Desired behavior

Same # of conjunction

What needs to be done

Run with different dt in cfg

How it can be tested

Run

gomezzz commented 2 years ago

Soooo, these results are from 3 different runs with identical simulation times and settings but different timesteps.

For completeness' sake, this is the used config (for dt = 5s, for 1s and 10s iteration are adjusted accordingly.)

sim:
  logLevel: info    # Available levels are off, critical, err, warn, info, debug, trace
  iterations: 6307200 # Number of simulation iterations
  maxAltitude: 10000 # Maximum satellite altitude [km]
  deltaT: 5.0 # [s]

  prop: # Which propagation model components should be applied
    useKEPComponent: true # Keplerian propagation
    useJ2Component: true # J2 spherical harmonic approximation
    useC22Component: true # C22 spherical harmonic approximation
    useS22Component: true # S22 spherical harmonic approximation
    useSOLComponent: true # Solar gravitational pull
    useLUNComponent: true # Lunar gravitational pull
    useSRPComponent: true # Solar radiation pressure
    useDRAGComponent: true # Atmospheric drag

io:
  output_file: output.csv # Output file name
  posFileName: pos.csv # Name of the CSV file with satellite positions [km]
  velFileName: v.csv # Name of the CSV file with satellite velocities [km/s]
  #progressOutputFrequency : 50 # Frequency of writing iteration number and encounter info to std::out.
  vtkWriteFrequency: 4000 # Frequency of writing to vtk file [iterations]
  #constellationList: StarlinkGen1,100,7200 # ';'-seperated constellations consisting of path (to constellation directory), start time, duration
  #constellationFrequency: 10 # Frequency of adding satellites to simulation [iterations]

autopas:
  logLevel: off
  cutoff: 80.0 # Cutoff for autopas force interaction
  skin: 160.0 # twice the maximum distance a particle can travel between two list rebuilds
  rebuildFrequency: 20
  desiredCellsPerDimension: 32 # Desired number of cells per dimension
  tuningMode: false
#  Newton3: "enabled, disabled"
#  DataLayout: "AoS, SoA"
#  Container: "varVerletListsAsBuilt, verletListCells"
#  Traversal: "vvl_as_build, vlc_c01, vlc_c18, vlc_sliced_balanced, vlc_slice_c02"

Let's start with looking at conjunction plots.

dt = 10

image

image

dt = 5

image

image

dt = 1

image

image

Orbital elements over time.

dt = 10

image

image

image

image

image

image

image

image

image

image

image

image

dt = 5

image

image

image

image

image

image

image

image

image

image

image

image

dt = 1

image

image

image

image

image

image

image

image

image

image

image

image

Closest particle distribution

dt = 10

100km

image

25km

image

10km

image

dt = 5

100km

image

25km

image

10km

image

dt = 1

100km

image

25km

image

10km

image

Waiting for #68 to analyze orbital elements vs conjunctions (we still have escaping particles atm)