Drone Launched Short Range Rockets

A concept of drone launched short range rockets (DLSRR) is presented. A drone or an aircraft rises DLSRR to a release altitude of up to 20 $km$. At the release altitude, the drone or an aircraft is moving at a velocity of up to 700 $m/s$ and a steep angle of up to 68$^o$ to the horizontal. After DLSRRs are released, their motors start firing. DLSRRs use slow burning motors to gain altitude and velocity. At the apogee of their flight DLSRRs release projectiles which fly to the target and strike it at high impact velocity. The projectiles reach a target at ranges of up to 442 $km$ and impact velocities up to 1.88 $km/s$. We show that a rocket launched at high altitude and high initial velocity does not need expensive thermal protection to survive ascent. Delivery of munitions to target by DLSRRs should be much less expensive than delivery by a conventional rocket. %% Even though delivery of munitions by bomber aircraft is even less expensive, a bomber needs to fly close to the target, while a DLSRR carrier releases the rockets from a distance of at least 200 $km$ from the target. %% All parameters of DLSRRs, and their trajectories are calculated based on theoretical (mechanical and thermodynamical) analysis and on several MatLab programs.


Introduction
In this work we discuss a possibility of a short range rocket launched from a fighter or a drone. The carrier aircraft or drone is called Rocket Launcher Aircraft (RLA). RLA can either be a modified fighter like F16, or an aircraft or drone constructed specifically for a rocket-launching mission. The rocket or rockets are called Drone Launched Short Range Rockets (DLSRR). In this work, we present the models and calculate performance of DLSRR, but not RLA. The author plans to write another work on RLA.
RLA raises DLSRRs to an altitude of 12 km -20 km and fires them. At the firing point, RLA and hence DLSRRs have a velocity of 500 m/s -700 m/s in the ground system of coordinates. Their motion is at an angle of 53 o to 68 o with respect to the horizontal. DLSRR has a slow-burning motor which burns 60 s -80 s. After the DLSRR motor burns out, it follows a ballistic trajectory. As the DLSRR reaches the apogee of 77 km -111 km, it releases one or more warheads which fly to the target. The DLSRR itself lands a few km away from the target.
Very little work has been done on military rockets with reusable first stage. Between 2010 and 2012, US Air Force had a project called "Reusable Booster System" to assist long-range rockets [1]. Jeffery Becker's concept is Theater Guided Missile Carrier (T-CVG), which is similar to our RLA. In Becker's concept, a modified version of reusable first stage of Falcon I rocket is used as RLA. This stage rises almost vertically, firing the second stage rockets which have a range of 1,600 km to 8,000 km [2]. Our proposal is more modest. It is likely, that multiple forms of military rocket systems with reusable first stage will exist.
Delivery of payload to target should be less expensive than conventional short range rocket systems. It may be noticed that conventional fighters and bombers also deliver munitions to a relatively long distance at a relatively low price. Unlike RLA, conventional fighters and bombers expose themselves to antiaircraft fire, which may make their operations inefficient. RLA never flies into hostile territory.
The current work consists of 6 Sections. In Section 1 we describe the state of art technology for long range precision strikes. In Subsection 1.1, we describe the state of the art rocket systems capable of making precise strikes to the range of 80 km to 600 km also mentioning some longer range missiles. In Tables 1 and 2, we present the cost of delivery per kg of payload. For ranges of 156 km to 600 km, the costs are from $ 2,700 per kg to $ 6,700 per kg depending on the missile. Thus we establish the need for a less expensive system. In Subsection 1.2, we describe historic and state of art guided munitions. In Subsection 1.3, we describe some antiaircraft missiles. The main reason the RLA -DLSRR system has not been developed so far is that bomber aircraft had capability of long range precision strike. The new generation of antiaircraft missiles may make bombers obsolete in any sizable confrontation.
In Section 2, we introduce the physics of rocket flight. We present the equations describing the rocket trajectory. We also introduce the concepts of Specific Impulse, gravity loss, and drag loss. These concepts are important for rocket development. The material in Section 2 is used to write programs and perform calculations with results presented in subsequent sections.
In Section 3, we discuss solid propellant rockets. In Subsection 3.1, we describe the rockets we need for DLSRRs In Subsection 3.2, we present properties of several fuel compositions. In Subsection 3.3, we discuss possible fuel grain shapes. We also discuss different rocket container materials along with their corresponding advantages and disadvantages.
In Section 4, we calculate rocket and warhead performance. In Subsection 4.1 we describe the performance of DLSRRs fired from RLA. Ranges obtained for these rockets are 263 km to 442 km depending on launch conditions and propellant burning time. In Subsection 4.2 we calculate the impact velocity of several projectiles. The Hypervelocity Projectile is described in Subsection 4.4.1. Europrojectile is a small tungsten projectile described in Subsubsection 4.2.2. Europrojectile has the best performance. In Section 5, we describe the aerodynamic heating and thermal protection. In Subsection 5.1 we present the general expressions for aerodynamic heating of different parts of the rocket. In Subsubsection 5.1.1, we present the equations for the stagnation point heat flux. In Subsubsection 5.1.2, we present the equations for general turbulent flow heat flux. In Subsubsection 5.1.3, we present the equations for the heat flux on rocket cylinder. In Subsubsection 5.1.4, we present the equations for the nose cone heat flux. The material in Subsection 5.1 is used to write programs and perform calculations with the results presented in Subsection 5.3. In Subsection 5.2 we discuss the heat shields. In Subsubsection 5.2.1, we describe radiation heat shields. In Subsubsection 5.2.2, we describe heat sink shields. In Subsection 5.3 we calculate the heating rates and temperatures of different parts of DLSRR during flight. In Subsubsection 5.3.1, we calculate stagnation point temperature as a function of time. The temperature never exceeds 1,200 o C -well within the capability of thermal protection. In Subsubsection 5.3.2, we calculate temperatures of several points on the nose cone as a function of time. The temperature never exceeds 960 o C. In Subsubsection 5.3.2, we calculate temperature DLSRR cylinder base as a function of time. The cylinder base temperature never exceeds 94 o C -at this temperature, aluminum does not soften.
In Section 6, we present a conclusion. In Appendix, we describe the MatLab programs we have written and used in this work.
1 Artillery Rockets and Warheads

Artillery Rockets
Short range ballistic missiles can make precise strikes at ranges from 150 km to 1,000 km. Some of the modern systems are extremely expensive. As we see from Table 1 and Table 2, the cost per 1 kg of munitions delivered to the target varies between $1,500 and $6,500 depending on the range.
Scud missile was the original Soviet short range ballistic missile. Scud-A introduced in 1957 delivered a 950 kg payload to a range of 180 km. Scud-B introduced in 1964 delivered a 985 kg payload to a range of 300 km. Over 7,000 Scud-B missiles have been built. Scud-C introduced in 1965 delivered a 600 kg payload to a range of 600 km. Scud-D introduced in 1989 delivered a 985 kg payload to a range of 700 km. Scud missiles have been used extensively in Iran-Iraq War. Between October 1988 and February 1992 approximately 1,700 to 2,000 Scud missiles were launched in Afghanistan's civil war [3]. Original Scud missiles had poor guidance, but in modern times, relatively inexpensive and precise guidance systems are available. Russian short range rocket Iskander delivers a 480 kg payload to a distance of 500 km or a 700 kg unit to a distance of 400 km. The missile has very precise guidance. Russia has only 136 such missiles. The price is kept secret, but should be at least $4 Million per unit [4, p.125].
American short range rocket ATACMS delivers a 230 kg payload to 300 km or a 560 kg payload to 165 km. The missile has very precise guidance. USA has 3,700 units [5]. The unit cost of ATACMS is $1.5 Million [6, p.19]. Israel has developed LORA ("Long Range Attack") missile, which delivers a 570 kg payload to a distance of 400 km. The missile has very precise guidance. The unit cost is not disclosed [7].
China has Dong-Feng ballistic missiles. Some of these missiles are intercontinental ballistic missiles. Main short-range ballistic missiles are DF-11 and DF-15. DF-11 launches an 800 kg payload to a range of 280 km or a 500 kg payload to a range of 350 km. DF-11A launches an 500 kg payload to a range of 350 km. DF-15 launches an 500 kg payload to a range of 600 km. DF-15A launches an 600 kg payload to a range of 600 km. DF-15B launches an 600 kg payload to a range of 700 km [8, p.48]. China possesses about 1,200 ballistic rockets with ranges 300 km to 1,000 km [9, p.31].
Artillery rockets can deliver guided submunitions much better than artillery projectiles, since they are not subject to extreme accelerations. Guided Multiple Launch Rocket System (GMLRS) is used by USA. M30 and M31 rockets deliver a 90 kg payload to 84 km. The missile has very precise guidance [10, p.112 The unit price is undisclosed [11]. China has several powerful rocket artillery systems. AR3 is a 370 mm caliber rocket which delivers a 200 kg warhead to a range of 220 km with great accuracy. Each launcher has 8 rockets [12].
USA did use the rockets described above in recent local conflicts. Out of 17,184 GMLRS Unitary rockets produced, 3 [13].
The measure of expense for a missile system at a given range is the cost per unit delivered weight. Below we tabulate the costs for the few solid propellant rockets for which the information is not classified: Liquid propellant rockets are less expensive, but they require about an hour to be prepared for firing. Some of them are listed in Table 2 below. The first two prices in 1990s are given in [15]. The prices listed below are double the old prices. 2,000 km 1,000 kg $ 5,100,000 $ 5,100 As can be seen from the above table, the cost of delivering munitions by short range ballistic missiles is very high. The costs of long range missiles are relatively low for their range.

Guided Munitions
One of the most important factors enabling long-range missiles and artillery is appearance and development of guided munitions. Precise guidance allows projectiles to strike targets tens or hundreds kilometers away with an accuracy of a few meters.
The most common guided munition is the air bomb. Guided bombs and glide bombs were first used by Germany in 1943 [17, p.4]. Since that time, relatively inexpensive guided bombs came into common use -during operations Enduring Freedom in 2001 and Iraqi Freedom in 2003, over 50% of all bombs used were guided [17, p. 20]. Joint Direct Attack Munition (JDAM) is a guidance kit which can be attached to any unguided bomb in order to transfer it into a guided bomb [17, p.213]. "Through FY 2005, 105,286 JDAM kits had been procured for an average unit-procurement cost of $21,379 each" [17, p.217]. Recently the Boeing corporation has produced 300,000 th JDAM [18]. Warheads used by the system we are proposing in this paper are similar to bombs dropped from an apogee of the second-stage rocket flight. They are much lighter than conventional bombs, but they will experience strong heating due to high velocity.
Guided projectiles and projectile guidance kits are much more hardy than bomb or missile guidance -they must survive acceleration of thousands of g's inside a cannon. Low cost guidance systems for cannon launched projectiles are becoming available [19]. The present generation of guided tube artillery munitions uses precision guidance kits attached to non-guided projectiles. The kit currently used by US Army is "M1156 Precision Guidance Kit (PGK)". PGK consists of a fuse placed in the nose of a non-guided projectile. M1156 uses GPS guidance. The projectile equipped with a fuse was suppose to have a circular error probability (CEP) of 30 m [23, p. 21]. Almost always new technology performs below expected parameters with M1156 being a nice exception. "During its first lot acceptance testing in April 2015, the M1156 demonstrated median accuracy of less than 10 m with a reliability of 97% when fired from the 155 mm/39-calibre M109A6 Paladin." [21, p.8] A German test had even better result. "From a distance of 27 km, 90 percent of the PGKequipped German shells landed within 5 meters of the target" [20]. PGK M1156 have a unit price of $10,000 [22].   [24]. One of the main advantages of M1156 is its simplicity. Indeed, (see Fig. 1) The kit has a single moving part, the canard assembly,which can only rotate along the longitudinal axis, the wings having a fixed cant; two couples of opposite wings have the same direction and thus provide lift, while the two despin wings provide counter-rotation [23, p. 21].
100,000 PGKs are ordered [21, p.8]. Israel is also developing a guidance kit called "Silver Bullet". The kit can be fitted to the nose of any conventional or rocket assisted projectile. This PGK does have four moving wings [23, p. 21-22]. Another guidance kit developed by Israel is TopGun [25]. Guidance kits used to guide projectiles should be able to guide warheads released by second-stage rockets.
The Hypervelocity Projectile (HPV) has been produced and tested. It is likely the next generation of guided munitions [26, p.12]. Its guidance should easily survive extreme acceleration within a launcher and heating within atmosphere. The projectile's guidance system should work even if it is shot from an electromagnetic cannon to a range of 400 km. HPV has a conical shape. It is made mostly out of tungsten. HPV's weight is 11.4 kg, length 64.2 cm, and caliber 7.83 cm. It is shown in Figure 2 below:

Antiaircraft Missiles
One of the main reasons the DLSRR technology has not been developed this far is that technology for precision strike at long range existed for many decades. As we have discussed in Subsection 1.2, bomber aircraft made it possible to bring guided munition to the target at a relatively low cost. Thus, there was no need for any unconventional delivery system.
In large-scale modern warfare, bomber aircraft may become a lost technology due to proliferation of very powerful anti-aircraft missiles. The main Russian surface to air missile in 2020 is S400. S400 missile with GRAU index 40N6 has a range of 400 km, ceiling of 185 km, and maximum velocity 2,000 m/s. The unit cost of the missile is $3.6 Million [27]. A new SAM called S500 should exceed the capabilities of S400. A new generation of SAMs may also employ Multiple Kill Vehicle (MKV). One rocket releases several guided submunitions each of which is capable of destroying the targeted aircraft [28, p.42].
USA also has very powerful antiaircraft missiles. SM3 Block 1B missile can destroy any aircraft or an incoming missile within 700 km. The warhead weighs 29 kg. The missile costs $12 million in 2018 dollars [29]. SM6 missile can destroy any aircraft within 300 km, it has a 64 kg warhead and a price tag of $4 million in 2018 dollars [30].
These missiles would make it impossible for any conventional bomber to operate in any sizable conflict. RLA does not approach the target closer then about 200 km. Moreover, it is possible to design an RLA which stays in the air only for a short period of time not exceeding 10 minutes.

The Physics of Rocket Flight
Consider a rocket that starts ascent at an altitude h 0 , which is generally 12 km to 20 km, with the initial speed is v 0 , which is generally 500 m/s to 700 m/s, with the angle with respect to horizontal line α 0 , generally between 51 o and 68 o . Figure 3: Forces acting on a firing rocket Three forces acting on a rocket during flight are presented in Figure 3. The gravitational force is M(t)g pointed in −ŷ direction. Notice, that while the rocket fuel is burning, the rocket mass is decreasing. The air resistance or drag force F d is acting in the opposite direction of velocity. The magnitude of F d is ρ is the air density, and A is the base area of the rocket. Let M be the rocket mass. Then thrust force F t acting in the direction of velocity is given by whereṀ is the fuel burning rate, and v e is the equivalent exhaust velocity. The exhaust velocity increases as the rocket rises out of dense atmospheric layers. For the DLSRR launched at a high altitude, v e barely changes. For some propellant grains, the fuel burning rate changes during the burn time. In our case we assume an almost steady burning ratė where t b is the burn time and f p is the propellant mass fraction given by Combined mass of propellant, rocket, and payload .

(2.4)
Since we know the forces acting on the rocket, we calculate its trajectory via the MatLab programs FirstStage.m and Rocket.m , which are described in the Appendix.
The change in velocity produced by the rocket engine is where M(t) is the rocket mass at time t. According to Tsialkovski Rocket Equation [31], where v e is the average exhaust velocity. In our case, v e ≈ 2, 100 m/s. For expensive rockets with flame temperatures in excess of 2,800 o C, v e ≈ 2, 600 m/s.
where M(t) is the rocket mass at time t, t 0 is the time it takes the rocket to reach the apogee, and At this point we define the effective loss of rocket velocity due to gravity. It is called gravitational loss and denoted v g . We define this loss in terms of the rocket's speed and altitude at the apogee. First, assume that the rocket is given its impulse v r instantaneously, and the drag loss is negligible. Such assumption is an abstract limit for a rocket, but it may be reality for a projectile fired at high altitude. Then the total kinetic and potential energy per unit rocket (projectile) mass at the beginning of trajectory is where h 0 is the launch altitude, v 0 is the launch velocity, and v r is the velocity gain due to the action of the rocket engine. Second, we incorporate the aerodynamic drag loss into (2.8) to obtain where v d is the drag loss. Third, we incorporate the gravitational loss into (2.9) to obtain where v g is the gravitational loss. In the above expressions, v 0 , v r , v d , and h 0 are known, while both E and v g are unknown. In order to calculate v g , we calculate E using the parameters of rocket motion at the apogee. The maximum ascent of the rocket, or the height the rocket attains at apogee is denoted h A . The rocket speed at the apogee is v A . The vector of this speed is in the forward direction. The total kinetic and potential energy per unit rocket mass at the apogee is which is the same as the rocket energy earlier in the path given in (2.10). Equating (2.10) and

Solid Propellant Rockets
In the present and the two following sections we consider a solid propellant rocket. It is schematically presented in Figure 4 below.

Rocket tube
Rocket tube The DLSRR should consist of a rocket with one or more guided warheads. These rockets would be much less expensive than conventional rockets with similar range and payload. First, the fact that the DLSRR starts at an altitude of over 12 km to 20 km with a velocity 500 m/s to 700 m/s would reduce △v needed for the rocket to achieve a given range by at least 750 m/s. The reduction in needed △v is greater than the initial velocity due to the fact that the rocket starts out at high altitude. First, it decreases the drag loss. Second, initial elevation itself reduces △v needed to achieve a given range. Second, as we show in the following sections, a rocket starting at high altitude with high initial velocity does not need to experience high firing acceleration. While most artillery rockets burn their fuel within 2-3 seconds, the DLSRR will be able to achieve good result even if it takes 40-60 seconds to burn its fuel. A solid-propellant rocket with slow burning fuel is much less expensive to manufacture. Third, slow-burning fuels have much lower flame temperature [32], which greatly reduces the cost of exhaust and initial guidance. There is a general pattern that propellants with higher combustion temperatures have higher burning rates than propellants with lower combustion temperatures at the same pressure. For example, ammonium nitrate propellants have low burning rates. Even though exceptions exist, they are rare. This phenomenon has a physical explanation. Lower temperature flames send lower heat flux back to the propellant. Lower heat flux causes the propellant to pyrolyse at lower rate. Thus, the overall regression and burning rate of the propellant is lower.
The initial guidance of DLSRR should be simple. It must make sure that DLSSR does not sway off the course by more than a kilometer. Probably inertial guidance with fins within the nozzle can accomplish the task.

Propellant Composition
Most modern artillery rockets as well as space rocket boosters use a solid propellant containing 70% ammonium perchlorate (AP), 15% aluminum powder, and 15% HTPB binder [33, p.479]. Such propellant is expensive to manufacture due to the fact that AP is highly explosive. Other propellants contain high explosives and highly nitrated nitrocellulose. All of them are expensive and dangerous to handle. All the aforementioned propellants have flame temperatures of 2,500 o C to 3,500 o C. Such exhaust temperatures require expensive nozzles.
DLSRR have a choice of fuels burning at lower temperature and rate. Some formulations contain ammonium nitrate (AN). The chemical formula for ammonium nitrate is NH 4 NO 3 . Its heat of formation is 367 kJ/mol. Ammonium Nitrate density is 1.725 g/cm 3 [34, p. 256]. Ammonium Nitrate undergoes phase transitions at -18 o C and +32 o C. These phase transitions are accompanied by volume change of about 4%, thus they must be avoided in order to avoid fuel grain destruction [34, p. 263]. Ammonium Nitrate can be phase-stabilized by addition of 10% potassium nitrate or 2% potassium fluoride [34, p. 267-268].
Other formulations have mildly nitrated nitrocellulose. The chemical formula for Nitrocellulose is C 6 H 10−x O 5 NO 2 x . Its atomic mass is 162 + 45x. Its binding energy is 961.5 − 103.6 · x kJ/mol. The number 0 < x < 3 determines the oxidizer content of nitrocellulose. Nitroglycerin (NG) is used in double base propellants with nitrocellulose [35]. Nitroglycerin formula is H 5 C 3 N 3 O 9 . Its atomic mass is 227 amu. Its heat of formation is 380 kJ/mol. The propellants we are interested in should have flame temperatures slightly below 1,500 o C.
In Table 3 below, we list several propellants. In the first column, we list propellant composition. In the second column, we list the flame temperature. In the third column, we list exhaust velocity into vacuum given an initial pressure of 40 atm and expansion of 15. Both temperature and the exhaust velocity are calculated using the program called Rocket Propulsion Analsis [36], which is available online. The forth column lists the burning rate at 40 atm. The fifth column lists power coefficient. For almost all propellants, the burning rate is approximated by where r b is the burning rate, P is pressure, P 0 is the reference pressure, and n is the power coefficient.

Propellant
Temperature  Propellants using ammonium nitrate oxidizer may experience burn rate instability for low and moderate pressures. More research is needed in order to determine feasibility and cost of using such propellants. Ammonium Dinitramide and Hydrazinium Nitroformate have been proposed as solid propellant oxidizers. Propellant formulations using these oxidizers have excellent performance. The problem with these oxidizers is their cost of hundreds of dollars per kilogram.
Finding a propellant optimal in terms of both cost and properties remains an open problem. For now we assume, that our propellant has flame temperature 1,450 o C, burning rate of 2.0 mm/s, and exhaust velocity of 2,100 m/s. Extrapolation of the data in [42, p.151] also predicts the aforementioned combination of burning rate and exhaust velocity. Using the densities of energetic materials in [42, p.37], the density of our propellant should be 1.6 g/cm 3 .

Grain Shape
The grain shape is the shape of the fuel within the rocket tube. Different rockets use a wide variety of grain shapes illustrated in Figure 5 below: The shape may change over the length of the rocket tube. We mention that different types of grain shapes provide different thrust curves. In the present paper, we focus on grain shape providing constant thrust. Propellant loading is the fraction of rocket tube volume occupied by the propellant. The best propellant grain shape for the DLSRR satisfies three requirements. First, the burning propellant surface area should experience minimal change as the propellant is burning. This would ensure steady thrust. Second, the propellant loading should be as high as possible. This would maximise fuel mass ratio within the rocket. Third, the design should be as simple as possible. This would minimize cost.
A diagram of a rocket motor cross-section is below: Aluminum 6061 T6 is one of the best candidates for the rocket wall. Aluminum 6061 T6 pipe of diameter up to 25 cm and wall thickness 1 cm is sold for $14.30 per kg pipe weight [43]. It has density of 2.7 g/cm 3 and tensile yield strength of 2,720 atm [44]. The mass of solid propellant rocket wall for a given volume and pressure is inversely proportional to the wall's specific strength, which is yield strength divided by density. Specific strength has units of Aluminum 6061 T6 has a specific strength of 1.15·10 5 J/kg. Duralumin and titanium have specific strength about twice as high as Aluminum 6061 T6, but they are tens of times more expensive than Aluminum 6061 T6. Composite materials with specific strength of 1.7 · 10 5 J/kg and price of $60 per kg in 2009 are available. These materials are widely used in compressed natural gas tanks [45]. Currently, these tanks are for sale at about $60 per kg wall material.  The rocket has diameter of 30 cm, total mass of 300 kg, total length of 4.3 m of which 1.0 m is the ogive warhead compartment. The propellant mass fraction is f p = 0.50. We took typical values for a 30 cm rocket, that can be found in the literature. The exhaust velocity is v e = 2, 100 m/s. Using RASAero software [46] we have calculated the drag coefficient for aforementioned artillery rocket at different Mach numbers and presented in Table 4 below. The RASAero software takes in rocket radius, rocket cylinder length, rocket nose length, and rocket nose shape. The software calculates the drag coefficient. The warhead compartment consists of a 10 kg container which ejects up to 60 kg of projectiles at the apogee. The projectiles are guided. They are made of high density material like tungsten or depleted uranium.
Performance of DLSRR for several release altitudes, release velocities, and firing times is present in Table 5      In this subsection we calculate the impact velocity for the warhead consisting of the Hypervelocity Projectile (HPV). HPV is described in Subsection 1.2 and presented in Fig. 2. Recall, that HPV's weight is 11.4 kg, length 64.2 cm, and caliber 7.83 cm. A 70 kg warhead compartment can carry up to 5 such projectiles. Using RASAero software [46] we have calculated the drag coefficient for the HPV at different Mach numbers and presented in Table 6 below. As we mentioned earlier, the RASAero software takes in rocket radius, rocket cylinder length, rocket nose length, and rocket nose shape. The software calculates the drag coefficient.  In Rocket.m simulations described in Subsection 4.1, we have calculated DLSRR apogee altitude, DLSRR apogee velocity, and DLSRR apogee horizontal distance. These calculations were performed for release altitudes of 12 km, 16 km, and 20 km. Corresponding release velocities were 500 m/s, 600 m/s, and 700 m/s. We input the results of Rocket.m simulations and the drag coefficient tabulated in Table 6 above into Impact.m. The program Impact.m calculates the projectiles' ranges, impact velocities, and descent drag loss. The results of calculations are presented in Table  7 below. The first row is the initial velocity. The second row is the firing altitude. The third row is the rocket firing time. The fourth row is the apogee altitude. The fifth row is the apogee velocity. The sixth row is the warhead horizontal range from the firing point. The seventh row is the impact velocity. The eighth row is the descent drag loss.  As we see from the above table, the HPV has very high aerodynamic drag loss. Even though minimization of drag loss has been an important goal in the design of HPV, the main design goal for HPV was an ability to survive the acceleration of a cannon launch. Rocket launched projectiles do not require the ability to survive extreme acceleration, hence HPV design may be suboptimal for rocket-launched projectiles.

European Hypersonic Projectile
The Europrojectile has not been built yet and exists only as a concept. Based on paper [47], we present a sketch below.
The volume of the cylindrical body of height 4.2 d and radius 0.35 d is The volume for LV-Haack nose is derived from the cone equations [48] as 0.56πLR 2 , where L is the cone length and R is the cone radius. Thus, in our case the cone volume is 0.50d 3 . The total volume of the projectile is 3.55d 3 . If the average density of the projectile is ρ, then its total mass is The projectile of our choice has caliber d = 7.5 cm, density ρ = 11 g/cm 3 , length 67.5 cm, and weight of 16.5 kg. A 70 kg warhead compartment can carry up to 3 such projectiles. Using RASAero software [46] we have calculated the drag coefficient for the HPV at different Mach numbers and presented in Table 8   We input the results of Rocket.m simulations and the drag coefficient tabulated in Table 6 above into Impact.m. The program Impact.m calculates the projectiles' ranges, impact velocities, and descent drag loss. The results of calculations are presented in Table 9 below. The first row is the initial velocity. The second row is the firing altitude. The third row is the rocket firing time. The fourth row is the apogee altitude. The fifth row is the apogee velocity. The sixth row is the warhead horizontal range from the firing point. The seventh row is the impact velocity. The eighth row is the descent drag loss.  From the above table we see that the Europrojectile has high impact velocity. At such velocity, the projectile can inflict considerable damage due to kinetic energy. One type of kinetic energy projectile which inflicts high area damage is a flechette projectile. This projectile explodes at a distance of about a hundred meters from the target releasing thousands of tungsten flechettes weighing 1 g to 2 g. For a projectile with high impact velocity, the flechettes cause extensive damage [49].

Aerodynamic Heating and Thermal Protection
An object moving through atmosphere at supersonic velocity experiences two types of aerodynamic heating. The first type is shock wave heating. The supersonic object produces a shock wave in front of it. Air temperature behind a shock wave can be as high as 20,000 o C, which has been experienced by the Galileo probe when entering the Jupiter's atmosphere [60, p.9]. For the rockets we are considering the temperature of the shocked air will never exceed "mere" 1,850 o C, but this temperature is still high.
The second type is skin friction heating. This is due to the high-speed friction of the rocket nose and sides against supersonic or hypersonic stream. Generally, this friction contributes more to overall dynamic heating than the shock wave.
In this work, we use simplified expressions to find an upper bound of heat flux. Exact calculation of heat flux requires extensive calculations using Aerodynamic Theory. Designing an actual thermal protection requires not only advanced theoretical and computational work, but also many laboratory experiments. Most rockets take thousands of expert-years to design [51].
The expressions themselves do not differentiate between heat flux caused by the shock wave and heat flux caused by friction. The heat transfer equations fall into two different classes. These classes represent two types of flow -laminar and turbulent. In this Section we present the expressions for heat flux generated by laminar and turbulent flows.
At the nose of the vehicle, the flow is always laminar. At a short distance from the nose, it becomes turbulent and stays turbulent for the rest of the vehicle length. Generally, turbulent flow transfers a higher heat flux to the wall then a laminar flow.
The expressions used to estimate laminar and turbulent heat fluxes on different types of surfaces have been obtained over decades as mathematical fits to experimental data. These expressions may or may not have deeper physical meanings. Some data patterns may have more than one expression providing a good fit [52,53,54,55,56,57].
One of the main advantages of DLSRR over conventional long-range ballistic missiles is the relatively low heating of DLSRR during ascent. The DLSRR reaches high velocities only in significantly rarified atmosphere. Thus, the heat flux and heat load are considerably lower than those experienced by a conventional long-range rocket.

General Expressions for Heat Flux
As a rocket flies through the air at high speed, it experiences aerodynamic heating. The rocket nose experiences laminar flow heating, while the nose cone and rocket cylinder experience turbulent flow heating. Even though some parts of the nose cone may experience laminar flow at high altitudes, turbulent flow heating rate still provides a valid upper bound. In this Subsection we present general expressions for the heat flux experienced by different points on a rocket surface during ascent.
Recovery temperature is the temperature the surface would have if it lost zero energy by thermal radiation or inward conduction. Recovery temperature is such that air at that temperature has enthalpy h aw . Adiabatic wall enthalpy, h aw is given by where h a is the ambient air enthalpy, v is the vehicle velocity, and r f is the recovery factor [52, p.9]. The recovery factor is a dimensionless number, which is 1 at a stagnation point, and 0.9 for a turbulent boundary flow [53, p.100]. A crude approximation for the recovery temperature is given by [53, p.100]: where T r is the recovery temperature, T a is the ambient air temperature.

Stagnation Point Heat Flux
The stagnation point is a point in a flow field where the local velocity of the fluid is zero [54, p.17]. The nose tip is the only stagnation point in the flow around a rocket. The laminar heat flux to the stagnation point can be approximated by [52,p.6,Eq. (40)]: where Q S denotes the total thermal energy absorbed by the stagnation point, andQ S is its time derivative. In Eq. (5.3) above, ρ is the air density measured in kg/m 3 , r is the nose radius in meters, h w is the specific enthalpy of the air at the wall, h aw is the specific enthalpy of the air at a hypothetical adiabatic wall, h a is the specific enthalpy of the ambient air, v is the rocket velocity in m/s with respect to air andQ is a constant. Sutton and Graves give a similar expression [55] witḣ According to Chapman [56],Q Detra and Hidalgo give a velocity dependent expression [57], where V kps is the rocket velocity in kilometers per second. In this work, we usė Q K = 1.83 · 10 −4 W /m 2 , which is the highest heating rate. By including the blackbody radiation emanating from the rocket, we obtain the total heat flux passing through the stagnation point: where σ B = 5.67 · 10 −8 W /m 2 K 4 is the Stefan-Boltzmann constant, and e w < 1 is the wall emissivity.

General Turbulent Flow Heat Flux
Turbulent flow heat flux is generally 3 to 6 times greater than laminar flow heat flux. The hypersonic flow around any rocket almost always starts out laminar and then turns turbulent. The transition from laminar to turbulent flow occurs at Reynolds number 10 5 ≤ Re ≤ 10 6 [58, p.19]. The Reynolds number is defined in Equations (5.11) and (5.13) below. The flow is turbulent over almost all surface area of the hypersonic vehicles we are considering in this paper.
Below we estimate the skin friction heat flux for general wall temperature and inclination. According to [52, p.iii], the skin friction heat transfer is defined in terms of Stanton number, The enthalpy supplied to the air by the shock wave travelling in front of the moving rocket is equal to v 2 /2. Thus, Substituting (5.10) into (5.9), we obtaiṅ The Stanton number at distance ry from the leading edge can be expressed in terms of skin friction coefficient [59, p.305]: The above equation is valid for ambient air temperature for most maneuvers. The effect of heating of air by friction and shock wave is described below. The skin friction coefficient is very well approximated by [62, p.4]: where T a is the ambient air temperature, T * is the air temperature at the boundary, µ a is the ambient air viscosity, and µ * is the viscosity at the boundary. (5.20) The range of validity for Eq. (5.20) above is wide enough to describe the conditions for all ascending rockets. In Equation (5.20) above, log represents log base 10. We express (5.20) in terms of natural logarithm: Combining (5.21) and (5.14), we obtaiṅ For the rockets and velocities we are considering, the Reynolds number varies between 1 · 10 6 and 1 · 10 8 . The quotient T * /T a varies between 1 and 8. For all values of Re and T * /T a considered, C f is strongly and strictly decreasing function of T * /T a . Hence, skin friction coefficient and turbulent heating are strongly decreasing with rising wall temperature.

Rocket Cylinder
For the rocket cylinder of DLSRR, the turbulent flow heat flux is calculated by (5.22) with the Reynolds number Re given by (5.17). At the base of the rocket cylinder, the distance from the stagnation point is x = 1 m. Substituting this distance into (5.17), we obtain the Reynolds number for the leading edge: For further parts of the rocket cylinder, the Reynolds number is higher and the heat flux is a little lower.

Nose Cone
Conical shape surface experiences 1.3 times the heat flux experienced by a flat plate [63, p. 11].
Another factor leading to an increase in turbulent flow heating rate on a cone relative to cylinder wall is the fact that the air pressure on the cone will exceed ambient atmospheric pressure. Any shape moving through the air stream with Mach number M having tangential angle θ to the air stream will incur a dynamic pressure of where P is the ambient atmospheric pressure. We calculate the total pressure exerted on an inclined surface in an air stream by adding ambient pressure and dynamic pressure. The total pressure will exceed the ambient pressure by a factor of where v is the rocket velocity in m/s, and θ is the angle between the cone surface and the oncoming air stream . The heat transfer rate given in (5.22) will be multiplied by the same amount. Thus, where T * is the temperature for which air enthalpy is given by (5.19), T a is the air temperature, T w is the rocket wall temperature, and Re is the Reynolds number given by (5.17).
Using the above information about the aerodynamic heating we have written a MatLab program ThermalAnalysis.m described in Appendix. This program calculates the temperatures of DLSRR stagnation point, several points on DLSRR nose cone, and rocket cylinder base during the powered portion of DLSRR ascent. The program takes in the time series data containing ambient air temperature, ambient air density, and rocket velocity. These time series are produced by the program Rocket.m, which is also described in Appendix.
Equation (5.8) is used to calculate stagnation point heat flux. Equation (5.26) is used to calculate the heat flux on the rocket nose cone caused by turbulent flow. All of the aforementioned fluxes are functions of the local rocket wall temperature. The program finds wall temperatures at the stagnation point and several points on rocket nose such that the absorbed heat flux is equal to the heat flux radiated away. These temperatures are determined for all time intervals during the rocket's powered ascent.
Equation (5.22) is used to calculate the heat flux on the rocket cylinder base caused by the turbulent flow. The cylinder does not radiate away any of absorbed heat. Being a heat sink shield, it absorbs the heat and gets hotter. The program calculates the cylinder base temperature at the end of the powered flight.
The program ThermalAnalysis.m performs calculations for rocket heating. The results of these calculations are presented in Subsection 5.3.

Thermal Protection Systems
There are three categories of heat shields. These are radiative shields, ablative shields, and heat sink shields. Radiative shields, which radiate away almost all of the heat flux they receive. They conduct almost no heat. DLSRR nose cone uses a radiative heat shield. Ablative shields absorb the heat flux. The heat is dissipated by pyrolysis and sublimation of the shield material [50, p.80]. DLSRR does not use ablative heat shields anywhere. Heat sink shields absorb heat flux. These shields generally consist of light metals which can absorb considerable thermal energy without becoming extremely hot. DLSRR cylinder rocket uses Aluminum 6061 T6 cylinder wall as a "inherent" heat shield. During the rocket flight the wall experience heating, but not enough to become dysfunctional.

Radiative Heat Shields
A typical radiative heat shield consists of three layers. The outermost layer is composed of a thin metal sheet. The metal should be resistant to oxidation and have high emissivity. Emissivity is the "ratio of the radiant flux emitted per unit area to that of an ideal black body at the same temperature" [41, p.2-45]. The next layer is flexible thermal insulation composed of ceramic fiber. The innermost layer is the underlying rigid structure. This is the shield used on DLSRR nose cone. In some works, the whole three layer shield is called a metal heat shield [64].

Outermost Layer -Oxidation Resistant Metals
The outermost layer of the metal heat shield is subject to both high heat flux and hypersonic oncoming wind. This layer must have high melting point and high oxidation resistance. Oxidation resistance of material is defined in terms of its recession rate under the influence of high temperature oxidizing environment. The recession rate depends on material, temperature, air density, and oncoming wind velocity. In some cases, recession follows a linear law: d r (t) = K linear t, (5.27) where d r (t) is the recession distance, t is the exposure time, and K linear is a linear recession constant. The constant K linear has units of m/s. Some oxidation-resistant metals are protected from oxidation by an outer oxide layer [65]. This layer is protective for temperatures beyond the oxide melting point. For instance, refractory metals such as tungsten and molybdenum have poor oxidation resistance due to the low melting points of their corresponding oxides. In cases where metal oxide does provide protection, metal recession follows a parabolic law: where d r (t) is the recession distance, t is the exposure time, and K parabolic is a parabolic recession constant. The constant K parabolic has units of m 2 /s. For temperatures up to 1,300 o C, Nickel has good oxidation resistance. In static air at atmospheric pressure and 1,300 o C temperature, it has parabolic oxidation rate with constant K parabolic = 2.5 · 10 −15 m 2 /s [65]. Under these conditions, it takes 10,000 seconds to oxidize a layer of nickel 0.005 mm thick. Inconel's oxidation resistance is slightly lower. In static air at atmospheric pressure and 1,300 o C temperature, it has parabolic oxidation rate with constant Chromium is oxidation-resistant up to 1,550 o C [70]. A layer of chromium oxide protects chromium from rapid oxidation at high temperatures. Pure chromium should have been protected by its oxide up to its melting point of about 1,900 o C. Unfortunately, the Cr 2 O 3 forms eutectic with lower oxides and melts away at 1,650 o C [70]. In air at atmospheric pressure and 1,550 o C temperature, chromium has parabolic oxidation rate with constant K parabolic = 1 · 10 −12 m 2 /s [71].
More expensive metals and alloys are needed to provide oxidation protection at higher temperatures [72]. For DLSRR nose cone, however, Nickel, Inconel, and Chromium are sufficient.
Second Layer -High Temperature Thermal Insulation Relatively inexpensive ceramic felt and cloth thermal isolation is available. CeramFiber thermal insulation works up to 1,260 o C. A sheet weighing 3 kg/m 2 costs $58 per m 2 [73]. Superwool HT Paper 2 mm thick has areal density of 0.45 kg/m 2 , thermal conductivity below 0.25 W /(m o K), and good properties up to 1,300 o C [74, p.46]. Zirconia Grade ceramic fiber paper works up to 1,420 o C. A sheet weighing 3 kg/m 2 costs $147 per m 2 [75]. All the aforementioned fibers can resist temperatures 100 o C higher than their classification for about a minute.

DLSRR Thermal Insulation
DLSRR nose cone has a simple radiative heat shield. The stagnation point and the first 10 cm of the cone are covered by a 1 mm thick chromium sheet. The second layer at that area consists of a 1 cm thick Zirconia Grade ceramic fiber paper. This combination can withstand temperatures of up to 1,520 o C which is 1,793 o K for duration of SRS ascent. The rest of the nose cone is covered by a 0.5 mm thick nickel or Inconel sheet. The second layer consists of CeramFiber blanket 1 cm thick. This shield section can withstand temperatures of up to 1,260 o C which is 1,533 o K.

Heat Sink Shields
The heat sink shield consists of three layers similar to metal radiative heat shield. The outermost layer is composed of metal. The next layer is flexible thermal insulation composed of ceramic fiber. The innermost layer is the underlying rigid structure. The difference between heat sink shield and radiative heat shield is the function of the outer layer. Whereas the outer layer of the radiative heat shield radiates away the heat flux, the outer layer of heat sink shield absorbs the thermal energy. Heat sink shields are most useful for absorbing short and powerful heat fluxes. These shields have been used for suborbital vehicles since 1950s [76, p.23].
Inherent heat sink shield is a part of a rocket or aircraft which has another primary function and also fulfils an extra function as a heat shield. For instance, a military aircraft may have armor which can act as a heat sink shield during short-time supersonic maneuvers. A solid rocket's metal shell mainly serves to contain the pressure of fuel grain combustion product. This shell may serve as a inherent heat sink shield as long as it does not overheat during hypersonic flight.
DLSRR has such inherent heat sink shield. As we mentioned in Section 4.1, the rocket cylinder consists of a 0.5 cm thick aluminum wall. Below we calculate the heat load such wall can withstand.
Rocket cylinder wall can be heated only so far as it retains tensile strength to hold the solid rocket motor. Aluminum 6061 T6 is the material composing the artillery rocket cylinder wall. Its Long-range rockets fired from the ground sustain an order of magnitude higher heat load. These rockets must have extensive outer heat shields. This is one of the reasons such rockets are expensive [53, p. 103-106].

DLSRR Heating
In this subsection we estimate the maximal heat load on various parts of the artillery rocket (DL-SRR) during its' ascent. The 30 cm DLSRR is described in great detail in Section 3.4. In this Subsection we describe the heating of three parts of DLSRR: the stagnation point, the nose cone and the rocket cylinder. The nose cone experiences greater heat flux than the rocket cylinder. The steps in heat flux calculations for the three parts are somewhat different.
The three parts of DLSRR have different heat shields. Both DLSRR stagnation point and DLSRR nose cone use a radiative heat shield described in Subsection 6.3. The temperature of this shield at any moment depends on the local heat flux. DLSRR cylinder uses an inherent heat sink shield described in Subsection 6.3. The temperature of this shield grows slowly almost in proportion to the absorbed heat load. DLSRR outer shell containing the three aforementioned parts is presented below: where σ B = 5.67 · 10 −8 W /m 2 K 4 is the Stefan-Boltzmann constant. The upper bound for stagnation point temperature is calculated by setting the maximum heat flux to zero: Notice, that in the above Eq. (5.31), the air enthalpy at wall temperature h w is a function of the wall temperature T w . The program ThermalAnalysis.m uses Eq. (5.31) to calculate stagnation point temperature. This calculation is done for every point on DLSRR trajectory.

DLSRR Nose Cone
During the rocket ascent, DLSRR nose cone experiences turbulent flow heating. This heat is radiated away by oxidized Inconel covering the stagnation point. The emissivity of oxidized Inconel surface is at least 0.75 [67, p.E-387]. The total heat flux to the nose cone is obtained by subtracting the radiative heat flux from the convective heat flux given in (5.26). For reader's convenience, we provide the formula for total heat flux with the notations described below formula (5.26). (5.32) The angle of inclination, θ , towards the air stream is different for different parts of nose cone. The upper bound for nose cone temperature is calculated by setting the maximum heat flux to zero: Notice, that in the above Eq. (5.33), both the air enthalpy at wall temperature h w and boundary air temperature T * are functions of the wall temperature T w . In order to solve Eq. (5.33), we must find the angle of inclination θ for different points on the nose cone. We visualize the nose cone below: Figure 11: Nose cone Any point A on the cone can be located in terms of an angular coordinate φ ∈ [0, 2 π) and a nonangular coordinate Y . If B is the projection of A onto Y axis then Y is the distance between A and the stagnation point. Heat flux and heat shield temperature are independent of φ due to the angular symmetry. Below, we calculate the heat flux and the heat shield temperature for several values of Y . For an ogive cone with radius 0.15 m and length 1.0 m, the inclination of any plane segment at coordinate Y is The program ThermalAnalysis.m uses Eq. (5.33) to calculate nose cone temperatures at different points on the nose cone as functions of time. This calculation is done for every point on DLSRR trajectory. The heat shield temperatures for DLSRR rocket for different firing altitudes, initial velocities, and rocket firing times are shown in Figure 12 below.   The heat flux at the stagnation point is proportional to √ ρ. The heat flux on the nose cone is proportional to ρ, where ρ is the ambient air density. As we see on the last plot of the last row, the heat flux 3 cm from stagnation point approaches the heat flux at the stagnation point as the rocket ascends.

DLSRR Cylinder
The turbulent flow heating rate is given by (5.26). All the heat is absorbed by the 0.5 cm thick aluminum rocket body. Almost no heat is radiated away. Assuming DLSRR body starts out at 20 o C, the program ThermalAnalysis.m calculates how much the rocket body heats up during the motor firing. After firing, the rocket body may absorb more heat, but it will have no effect on rocket performance. The temperature of the rocket cylinder base after firing is tabulated in Table  10.

Conclusion
In this work we have demonstrated the feasibility of drone launched short range rockets launched from a rocket launcher aircraft. DLSRR deliver warheads to targets at a distance of 263 km to 442 km with impact velocity of 1,400 m/s to 1,880 m/s. RLA never approaches the target closer than 220 km making it much less vulnerable to anti-aircraft weapons than modern fighters and bombers. The cost of delivery of munitions by the DLSRR should be much lower than the cost of delivery of munitions by modern short range ballistic missiles (SRBM). This is the case due to the greatly reduced requirements placed on the DLSRR compared to TBM. First, DLSRR require much lower △v than SRBM. Second, DLSRR can burn fuel over 30 s to 80 s which is much slower than most modern TBMs. Third, as we have shown in Subsection 5.3, the DLSRR reach high velocity in high stratosphere and experience only moderate heating. Stagnation point temperature never exceeds 1,200 o C, and nose cone temperature never exceeds 960 o C. Such temperatures require relatively inexpensive heat protection.
A DLSRR having an initial mass of 300 kg can deliver a 60 kg charge to the target. The charge can consist of one or more projectiles. The ranges and impact velocities for a warhead containing four Europrojectiles are listed in Table 9. High release altitude and velocity combined with shorter burning times increase the warhead's range and impact velocity. The same factors increase launch cost.
Rocket.m performs calculations for the DLSRR. The user inputs DLSRR's mass, diameter, drag coefficient, propellant mass fraction, exhaust velocity, initial inclination and propellant burning time. Effective frontal area of the rocket is While the rocket engine is burning the thrust is obtained by combining (2.2) and (2.3): where M is the rocket mass, f p is the propellant mass fraction, t b is the engine burning time, and v e is the exhaust velocity. Using this data, the program starts to perform iterative calculations on the rocket.
At time t = 0, Rocket.m has the following: rocket's mass M; rocket's coordinates in xy -plane, where x is the horizontal distance and y is the altitude; rocket's velocity v -which has horizontal component v x and vertical component v y . Rocket.m uses the rocket's altitude to determine air density ρ and sonic velocity v s at time t = 0. For air density as a function of altitude, we use US Standard Atmosphere 1976. Rocket.m calculates the Mach number as M = v/v s . The drag coefficient C d is determined from the Mach number. The drag force is given by by (2.1): The drag force acts in the direction opposite of rocket motion. After the rocket's engine burns out, F d = 0. Rocket acceleration is derived from all forces acting on the rocket as shown in Figure 3: wherev is a unit vector in the direction of rocket motion, andŷ is a unit vector in y or upward direction.
Rocket.m uses the information given in the above paragraph to calculate similar information at a time dt. If the rocket's engines are still burning, then the rocket's mass has decreased by the mass of propellent consumed: otherwise the rocket mass stays constant. The rocket's velocity has changed by v(t + dt) = v(t) + a(t) dt.
The rocket's position has changed by r(t + dt) = r(t) + v(t) dt.
Based on the new mass, position, and velocity, Rocket.m calculates the rocket's drag force and acceleration at time dt.
Rocket.m performs another iteration obtaining the rocket's drag force and acceleration at time 2 dt. Rocket.m performs many iterations, until the rocket reaches its apogee -the maximal ascent. Generally, dt = 0.1 s and it takes 140 s to 160 s for the rocket to reach the apogee. Rocket.m performs 1,400 to 1,600 iterations within a few seconds.
Rocket.m calculates DLSRR's flight apogee altitude y, flight apogee horizontal velocity v x , and flight apogee horizontal coordinate x. Recall, that apogee vertical velocity is v y = 0. Rocket.m also calculates DLSRR's drag loss by Eq. (2.7) and gravity loss by Eq. (2.12).
Rocket.m produces a time series of DLSRR's coordinates (x, y), and velocity v = v x , v y with time interval dt. The time series of ambient air densities and temperatures determined by altitude y as well as the Mach numbers are also recorded.
Impact.m performs calculations on the projectiles released at the apogee. The user inputs the projectile's mass, diameter, drag coefficient, apogee altitude, apogee, velocity, and apogee horizontal coordinate. The program calculates impact distance and velocity.
Mathematically, Impact.m is a simplified version of Rocket.m. As the rocket reaches the apogee, it releases one or more warhead. The warhead(s) keep the rocket's horizontal and vertical coordinates (x, y). They keep the rocket's horizontal velocity v x , and vertical velocity, which is 0 at the apogee. The warhead(s) have a new caliber d 1 ≪ d, and hence a new frontal area They also have a new mass M 1 ≪ M, which stays the same throughout the descent. The warhead(s) have new drag coefficient dependance on Mach number C d1 (M ).
As the warhead(s) descend from the apogee to the target, Impact.m calculates their acceleration, velocity, and coordinate for each time increment dt. Drag force is calculated by Eq. (A.1) with the new frontal area and drag coefficient function. Acceleration is calculated by wherev is a unit vector in the direction of rocket motion, andŷ is a unit vector in y or upward direction.
Impact.m calculates warhead's descent time, impact velocity, and horizontal distance of the impact. Impact.m also calculates the warhead's drag loss by Eq. (2.7). The warhead experiences no gravity loss during the descent. The combination of Rocket.m and Impact.m calculate the warheads horizontal range from the point of rocket release, and the flight time from rocket firing to impact.
ThermalAnalysis.m calculates the temperatures of SSR stagnation point, several points on SSR nose cone, and rocket cylinder base during the powered portion of SSR ascent. The program takes in the time series data containing ambient air temperature, ambient air density, and rocket velocity. These time series are produced by Rocket.m.
For every time interval dt, ThermalAnalysis.m solves Eq. (5.31) for the stagnation point wall temperature. Likewise, for every time interval dt, ThermalAnalysis.m solves Eq. (5.33) for wall temperature at several points on the rocket's nose. Notice, that the temperatures are found for each time independently of temperatures in the previous times.
ThermalAnalysis.m uses Eq. (5.22) to calculate the heat flux on the rocket cylinder base caused by the turbulent flow. The cylinder does not radiate away any of absorbed heat. Being a heat sink shield, it absorbs the heat and gets hotter. ThermalAnalysis.m assumes that the rocket cylinder base starts at 20 o C, and calculates the time series of rocket base temperatures as the rocket base absorbs heat flux. ThermalAnalysis.m displays the cylinder base temperature at the end of rocket firing.
Rocket wall temperature time series at stagnation point as well as distances of 3 cm, 10 cm, and 30 cm from the stagnation point are plotted in Figure 12. Cylinder base temperature at the end of rocket firing are presented in Table 10.