Example - 48 - Crew Module Atmospheric Re-entry Experiment

[1]:
from AMAT.planet import Planet
from AMAT.vehicle import Vehicle
[2]:
import numpy as np
import matplotlib.pyplot as plt

This notebook simulates the atmospheric entry of Crew Module Atmospheric Re-entry Experiment. https://www.isro.gov.in/sites/default/files/flipping_book/LVM-3/CARE/files/assets/common/downloads/LVM3-brochure.pdf

[3]:
# Set up the planet and atmosphere model.
planet=Planet("EARTH")
planet.h_skip = 126.0E3
planet.h_trap = 0.0E3
planet.loadAtmosphereModel('../atmdata/Earth/earth-gram-avg.dat', 0, 1, 2, 3)
[4]:
# Set up the vehicle
vehicle1=Vehicle('care', 3735.0, 350.0, 0.0, np.pi*3.1**2.0*0.25, 0.0, 3.0, planet)
[5]:
# Set up entry parameters
vehicle1.setInitialState(126.0,0.0,0.0,5.326,0.0,0.0,0.0,0.0)
[6]:
# Set up solver
vehicle1.setSolverParams(1E-6)
[7]:
# Propogate vehicle entry trajectory
vehicle1.propogateEntry (120.0*60.0,0.1,0.0)
[8]:
# import rcParams to set figure font type
from matplotlib import rcParams
[9]:
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/crew-module-atmospheric-re-entry-experiment-earth.png',bbox_inches='tight')
plt.savefig('../plots/crew-module-atmospheric-re-entry-experiment-earth.pdf', dpi=300,bbox_inches='tight')
plt.savefig('../plots/crew-module-atmospheric-re-entry-experiment-earth.eps', dpi=300,bbox_inches='tight')

plt.show()
../_images/examples_example-48-crew-module-atmospheric-re-entry-experiment-earth_10_0.png
[10]:
plt.figure(figsize=(6,6))
plt.plot(vehicle1.stag_pres_atm[0:2700], vehicle1.q_stag_total[0:2700], '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/crew-module-atmospheric-re-entry-experiment-earth-thermal.png',bbox_inches='tight')
plt.savefig('../plots/crew-module-atmospheric-re-entry-experiment-earth-thermal.pdf', dpi=300,bbox_inches='tight')
plt.savefig('../plots/crew-module-atmospheric-re-entry-experiment-earth-thermal.eps', dpi=300,bbox_inches='tight')

plt.show()
../_images/examples_example-48-crew-module-atmospheric-re-entry-experiment-earth_11_0.png