Example - 40 - Oceanus - Uranus (Concept)

from AMAT.planet import Planet
from AMAT.vehicle import Vehicle

import numpy as np
import matplotlib.pyplot as plt

This notebook simulates the atmospheric entry of the Uranus entry probe which was part of the Oceanus mission study https://doi.org/10.1016/j.asr.2017.02.012

# Set up the planet and atmosphere model.
planet=Planet("URANUS")
planet.h_skip = 1000.0E3
planet.h_trap = -100.0E3
planet.loadAtmosphereModel('../atmdata/Uranus/uranus-ames.dat', 0 , 1 , 2, 3)

# Set up the vehicle
vehicle=Vehicle('Oceanus-Uranus-Probe', 120.0, 170.0, 0.0, np.pi*0.8**2.0*0.25, 0.0, 0.28, planet)

# Set up entry parameters
vehicle.setInitialState(1000.0,0.0,0.0,21.47,0.0,-31.5,0.0,0.0)

# Set up solver
vehicle.setSolverParams(1E-6)

# Propogate vehicle entry trajectory
vehicle.propogateEntry (100*60.0,0.1,0.0)

# import rcParams to set figure font type
from matplotlib import rcParams

fig = plt.figure(figsize=(12,8))
plt.rc('font',family='Times New Roman')
params = {'mathtext.default': 'regular' }
plt.rcParams.update(params)

plt.subplot(2, 2, 1)
plt.plot(vehicle.v_kmsc, vehicle.h_kmc, 'k-', linewidth=2.0)
plt.xlabel('Speed, km/s',fontsize=14)
plt.ylabel('Altitude, km', fontsize=14)
ax=plt.gca()
ax.tick_params(direction='in')
ax.yaxis.set_ticks_position('both')
ax.xaxis.set_ticks_position('both')
ax.tick_params(axis='x',labelsize=14)
ax.tick_params(axis='y',labelsize=14)

plt.subplot(2, 2, 2)
plt.plot(vehicle.acc_net_g, vehicle.h_kmc, 'b-', linewidth=2.0)
plt.xlabel('Deceleration, Earth g',fontsize=14)
plt.ylabel('Altitude, km', fontsize=14)
ax=plt.gca()
ax.tick_params(direction='in')
ax.yaxis.set_ticks_position('both')
ax.xaxis.set_ticks_position('both')
ax.tick_params(axis='x',labelsize=14)
ax.tick_params(axis='y',labelsize=14)

plt.subplot(2, 2, 3)
plt.plot(vehicle.q_stag_total, vehicle.h_kmc,'r-', linewidth=2.0)
plt.xlabel('Stagnation point heat-rate, '+r'$W/cm^2$',fontsize=14)
plt.ylabel('Altitude, km', fontsize=14)
ax=plt.gca()
ax.tick_params(direction='in')
ax.yaxis.set_ticks_position('both')
ax.xaxis.set_ticks_position('both')
ax.tick_params(axis='x',labelsize=14)
ax.tick_params(axis='y',labelsize=14)


plt.subplot(2, 2, 4)
plt.plot(vehicle.heatload/1.0E3, vehicle.h_kmc, 'm-', linewidth=2.0)
plt.xlabel('Stagnation point heat-load, '+r'$kJ/cm^2$',fontsize=14)
plt.ylabel('Altitude, km', fontsize=14)
ax=plt.gca()
ax.tick_params(direction='in')
ax.yaxis.set_ticks_position('both')
ax.xaxis.set_ticks_position('both')
ax.tick_params(axis='x',labelsize=14)
ax.tick_params(axis='y',labelsize=14)

plt.savefig('../plots/oceanus-uranus-probe.png',bbox_inches='tight')
plt.savefig('../plots/oceanus-uranus-probe.pdf', dpi=300,bbox_inches='tight')
plt.savefig('../plots/oceanus-uranus-probe.eps', dpi=300,bbox_inches='tight')

plt.show()