Example - 47 - Venus Flagship - Decadal Study - 2010 (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 Venus Flagship Mission, 2010 Decadal Mission Concept Study (page 159). https://vfm.jpl.nasa.gov/files/Venus+Flagship+Mission+Study+090501-compressed.pdf

# Set up the planet and atmosphere model.
planet=Planet("VENUS")
planet.h_skip = 220.0E3
planet.h_trap = 10.0E3
planet.loadAtmosphereModel('../atmdata/Venus/venus-gram-avg.dat', 0, 1, 2, 3)

# Set up the vehicle
vehicle1=Vehicle('venus-flagship', 2020.0, 349.0, 0.0, np.pi*2.65**2.0*0.25, 0.0, 0.30, planet)

# Set up entry parameters
vehicle1.setInitialState(220.0,0.0,0.0,11.1,0.0,-13.8,0.0,0.0)

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

# Propogate vehicle entry trajectory
vehicle1.propogateEntry (120.0*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(vehicle1.v_kmsc, vehicle1.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(vehicle1.acc_net_g, vehicle1.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(vehicle1.q_stag_total, vehicle1.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(vehicle1.heatload/1.0E3, vehicle1.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/decadal-2010-venus-flagship.png',bbox_inches='tight')
plt.savefig('../plots/decadal-2010-venus-flagship.pdf', dpi=300,bbox_inches='tight')
plt.savefig('../plots/decadal-2010-venus-flagship.eps', dpi=300,bbox_inches='tight')

plt.show()

plt.figure(figsize=(6,6))
plt.plot(vehicle1.stag_pres_atm[0:1700], vehicle1.q_stag_total[0:1700], 'r-', linewidth=2.0)
plt.xlabel('Stagnation pressure, atm',fontsize=14)
plt.ylabel('Stagnation point heat-rate, '+r'$W/cm^2$',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/decadal-2010-venus-flagship-thermal.png',bbox_inches='tight')
plt.savefig('../plots/decadal-2010-venus-flagship-thermal.pdf', dpi=300,bbox_inches='tight')
plt.savefig('../plots/decadal-2010-venus-flagship-thermal.eps', dpi=300,bbox_inches='tight')

plt.show()