# Design Methodology and Performance Evaluation of New Generation Sounding Rockets.

1. IntroductionSounding rockets are specialized missiles generally used to investigate the region of space between 50 and 700 km where it is difficult to enter with the traditional atmospheric balloon or low orbit satellites. They can be used in several missions, such as detection of the solar activity and anomalies, analysis of the constituents of the upper atmosphere, thermal analysis on new materials, and generally, measurements of the space surrounding the Earth [1-3]. The most important characteristic of a sounding rocket is the capability of achieving microgravity conditions for the payload, without the need of a manned mission in the ISS to perform experiments and then reducing the costs of the mission. The quality of microgravity depends on the absence of drag and other gyroscopic forces, so the guidance of the rocket assumes the utmost importance. Another good characteristic of the sounding rocket in contrast to the satellite missions is the sounding rocket's ability to retrieve the payload with the help of a heat shield and a parachute. For all of these reasons, a vast number of sounding rocket families were developed, such as the American "Black Brant" and the European "Maxus." The Black Brant family includes sounding rockets that range from one through three stages. The Black Brant VC [4] (Figure 1) can be launched from a steerable launch tower with different kick angles to achieve different apogee altitudes. It is controlled via spin-up motors and canard wings, which provide a fast response at the cost of small instability. Four rear tail surfaces guarantee an adequate static margin. The Maxus [4, 5] is a very powerful European single-stage sounding rocket launched in Kiruna (Sweden) and can reach a maximum altitude of 700 km, while the range does not exceed 100 km. It has a microgravity duration of about 12 min, and the microgravity level [DELTA]g is near [10.sup.-4]m/[sec.sup.2]. TVC control system is highly sophisticated in order to guarantee strict margins for the range and to avoid damages to both people and goods. Table 1 portrays some useful data to compare the performance of the two rockets.

Both of these sounding rockets are propelled with a solid state motor, due to their simplicity and light weights with respect to liquid propellant. A crucial point is that a sounding rocket has to be cheaper than a high-tech rocket, and in order to achieve this objective, new design methodologies are needed. Many authors (e.g., Chowdhury et al. [6, 7] and Casalino and Pastrone [8, 9]) are interested in designing new sounding rockets to explore the potentiality of new hybrid motors for hypersonic flight. Nonaka et al. [10] and Woo et al. [11] investigated the reusability of the sounding rocket, and Cremaschi et al. [12] investigated the optimization of a trajectory with the use of a commercial trajectory simulation tool (ASTOS). However, a very limited number of studies are concerned with fast sizing and aerodynamic losses during the ascent path. The purpose of this work is to find a novel and simple methodology to size and design for future sounding rockets, considering accurate aerodynamics and three guidance algorithms. In particular, the technique explained in this work is applied to a single-stage rocket whose three-dimensional trajectory is simulated using the following guidance schemes:

(i) Constant attitude guidance

(ii) Near radial guidance

(iii) Sun-pointing guidance

This work purpose is to find a procedure that can provide the minimal volume rocket that can guarantee a fixed time of microgravity with different scientific payloads inside. Definitely, the methodology at hand yields configurations that can provide a well-developed preliminary design for future sounding rockets, with an accurate aerodynamic model and using three distinct guidance schemes.

2. Rocket Geometry and Sizing

2.1. Basic Geometry. With the purpose of maintaining the configuration as simple as possible, the sounding rocket is composed of one stage and a module hosting the avionics, experiments, and recovery system (parachute and the heat shield).

The basic geometry of the rocket is composed of the following:

(i) A nose cone of parabolic shape

(ii) A fuselage of cylindrical shape

(iii) Four canard wings in cruciform configuration

(iv) Four tail surfaces of trapezoidal shape

This configuration is chosen with the purpose of creating an algorithm as general as possible. Incidentally, this geometry is adopted by a lot of sounding rockets in the market, like the "Black Brant VC" in Figure 1.

Every rocket configuration in this work is defined by the total length l and the maximum diameter d, while all the other parameters like the propulsive ones are obtained through scaling factors. The grid used for the computation is

(i) l from 3 m to 15 m with step of 1 m

(ii) d from 0.2 m to 1.5 m with step of 0.05 m.

Thus, a total of 243 different configurations were analyzed. The design methodology that is being presented relates geometry, weight, aerodynamics, and propulsive performances.

2.2. Main Body. For each configuration, the length of the ogive is proportional to the total length of the sounding rocket

[l.sub.ogive] = [C.sub.ogive] l, (1)

where [C.sub.ogive] = 0.17 is a dimensionless coefficient. The dimensional ogive point distribution is then obtained by multiplying the nondimensional points in the ([xi], [eta]) frame in Figure 2 with [l.sub.ogive] and d. The main body, which contains the propulsive system and the avionics, is cylindrical with diameter d, while its length is found through

[l.sub.body] = l - [l.sub.ogive]. (2)

2.3. Aerodynamic Surfaces. The aerodynamic surfaces are sized according to the following ratios:

[mathematical expression not reproducible]. (3)

Then the following quantities are computed as follows:

[h.sub.T] = [square root of (4/3 [A.sub.C])], [h.sub.C] = [square root of (2 [A.sub.T])], (4)

while the distances with respect to the nose tip are as follows:

[d.sub.C] = [l.sub.ogive] + 0.4m, [d.sub.T] = l - [h.sub.T]. (5)

The main concept under this procedure is to maintain the shape of the wings, while selecting their size in relation to the overall size of the main body. This is in order to avoid wings with nonconventional shapes that would result in undesired aerodynamic effects. Figure 2 portrays the wing geometry with the profile adopted for both wing systems. The nonvariable distance in (5) 0.4 m is meant to provide a fixed volume needed for the canard avionic and the recovery system. Figure 2 portrays an example of a configuration with the disposition of the leading edges along the longitudinal axis of the rocket.

2.4. Mass Scaling. In order to describe properly a real sounding rocket mission, a fixed mass including parachute, heat shield, and metal case is used for every configuration. Table 2 portrays in detail these auxiliary masses, which form [m.sub.aux]

With regard to the initial mass at launch, a common initial density is used in order to find a reasonable and scaled mass distribution for every rocket configuration. This initial density is [mathematical expression not reproducible], which is taken from the Black Brant VC and is considered to have a reasonable value. This value includes propellant, structural, and auxiliary masses. So the initial mass for every configuration without the ogive mass is

[mathematical expression not reproducible]. (6)

For every configuration. 50, 100, and 200 kg of scientific payload are used, so the total initial mass at launch is

[m.sub.0] = [m.sub.0body] + [m.sub.U] + [m.sub.aux]. (7)

The propulsive technology of this single-stage rocket is a solid propellant one, allowing for the following structural coefficient:

[epsilon] = [m.sub.S]/[m.sub.P] + [m.sub.S] = 0.11. (8)

The structural and propellant masses are obtained through the following expressions:

[mathematical expression not reproducible]. (9)

In the end, (6), (7), (8), and (9) yield the overall rocket mass at launch, portrayed in Figure 3 as a function of l and d, together with the corresponding propellant mass.

2.5. Propulsion System Scaling. To scale the thrust for every sounding rocket size, the initial launch acceleration and the specific impulse are maintained constant:

T = [n.sub.0] g[m.sub.0], (10)

where [n.sub.0] = 3g, while the propellant mass consumption is

[??] = -T/[I.sub.sp]g. (11)

The value for [n.sub.0] is chosen in order to constrain the rocket dynamical stress at launch. Equation (11) leads to the burnout time:

[t.sub.B] = [m.sub.P] / [??]. (12)

Figure 4 portrays the thrust and the burnout time scaling with the rocket size.

3. Rocket Aerodynamic Modeling

Aerodynamic forces affect the rocket motion, which partially occurs in the atmosphere. Thus, the rocket aerodynamics is to be modeled appropriately. In this work, four quantities are assumed to determine the aerodynamic forces: Mach number M, Reynolds number Re, angle of attack [alpha], and sideslip angle [beta]. To provide an accurate model for each configuration, two fundamental steps are needed:

(1) Derivation of the aerodynamic coefficients [C.sub.D], [C.sub.L], and [C.sub.Q] from DATCOM [13] at a relevant number of Mach, Reynolds, angles of attack, and sideslip angles

(2) Cubic interpolation of the aerodynamic coefficients during the integration of the dynamics equations

Following the approach presented by Mangiacasale [14], Pallone et al. [15], and Pontani [16], the aerodynamics of the rocket was modeled through the Missile DATCOM software [13]. This software provides an accurate esteem of the aerodynamic forces without using a lot of computational resources; thus, it can provide a complete aerodynamic database for each configuration in a short time. The following discrete values of M, Re, [alpha], and [beta] are used:

(i) M = [0.2, 0.5, 0.6, 0.8, 0.9, 1.1, 1.2, 1.5, 2, 2.5, 3, 3.5, 4, 4.5, 5, 5.5, 6, 6.5, 7, 7.5, 8, 10]

(ii) Re = [0.5, 1, 1.5, 2] x [10.sup.5]

(iii) [alpha] = [0, 1, 2, 3, 4, 5] deg

(iv) [beta] = [0, 1, 2, 3, 4, 5] deg

This means that DATCOM was run 3024 times for each (l, d) combination. The aerodynamic surfaces used in the computation are the cross surfaces of each configuration and correspond to the circular section, that is,

[S.sup.(i)] = [pi] [(d/2).sup.2]. (13)

Figure 5 portrays an example of the aerodynamic coefficients for a single configuration.

4. Rocket Dynamics

The problem investigated in this work consists in finding the sounding rocket with the minimum total volume that can achieve microgravity durations not less than 10 or 15 minutes. In the absence of significant attitude angular velocities, a rocket is considered in microgravity conditions when it is flying at an altitude above 100 km [5].

This approximation is based on real evidence and calculations, allowing for it to be used as a practical rule to analyze the rocket performance. Thus, the problem at hand reduces to finding the minimal rocket volume (l, d) that can guarantee [DELTA][T.sub.mg] of 10 and 15 min, fixing the propulsive technology ([I.sub.sp], [epsilon]) and the geometrical ratios between its parts ([mathematical expression not reproducible]).

4.1. Reference Frames. A description of the reference frames used in this work is useful to understand the guidance implementation. In the following, the notation [R.sub.j] ([- or +] [xi]) (j = 1, 2, 3) denotes the elementary rotation about axis j by angle [xi] while the [- or +] sign means, respectively, clockwise and counterclockwise.

(1) The Earth-centered frame (ECI) [mathematical expression not reproducible] is a Cartesian inertial reference frame defined as follows. Its origin O is the center of the Earth. The unit vector [[??].sub.3] is aligned with the Earth axis of rotation and is positive northward, whereas [[??].sub.1] is aligned with the vernal axis, which corresponds to the intersection of the Earth equatorial plane with the ecliptic plane.

(2) The Earth-centered Earth-fixed frame (ECEF) [mathematical expression not reproducible] represents a reference system that is rigidly attached to the Earth. Under the assumption that the planet rotates with a constant angular rate [[omega].sub.E], also the ECEF frame rotates with angular rate [[omega].sub.E]. The origin of the coordinate system is the center of the Earth, while the unit vector [??] intersects the Greenwich reference meridian at all times, and [??] is aligned with the planet axis of rotation and is positive northward. [[omega].sub.E][??] is the vector rotation rate of the ECEF frame with respect to the ECI frame. The reference meridian is the Greenwich side real-time [[theta].sub.G], given by [[theta].sub.G](t) = [[theta].sub.G]([bar.t]) + [[omega].sub.E](t - [bar.t]) where [bar.t] refers to a generic time instant. The ECEF frame is obtained through a counterclockwise rotation by angle [[theta].sub.G] around the third axis, that is, [mathematical expression not reproducible].

(3) The local horizontal frame (LH) [mathematical expression not reproducible] is a reference system with the origin in the rocket center of mass. Specifically, [??] is the unit vector aligned with the instantaneous position [mathematical expression not reproducible] is directed along the local east direction, and [??] is aligned with the local north direction. LH frame is obtained from the ECEF frame through the following rotations: [mathematical expression not reproducible].

(4) The auxiliary orbital frame (AO) [mathematical expression not reproducible] is defined making reference to the rocket relative velocity [v.sub.R] = [v.sub.I] - [[omega].sub.E] x r. The unit vector [[??].sub.R] is aligned with the projection of the velocity [v.sub.r] onto the horizontal plane [mathematical expression not reproducible]. The AO frame is obtained from the LH frame with the simple rotation [mathematical expression not reproducible].

(5) The relative velocity frame (RV) [mathematical expression not reproducible] is obtained from the AO frame through the following clockwise rotation: [mathematical expression not reproducible].

(6) The wind axis frame (WA) [mathematical expression not reproducible] is defined with reference to the rocket velocity with respect to the local velocity of the atmosphere [v.sub.a] = [[omega].sub.E] x r so [v.sub.a] = [v.sub.R] and just a counterclockwise rotation about axis 2 by the bank angle [sigma]: [mathematical expression not reproducible].

(7) The auxiliary body axis frame (ABA) [mathematical expression not reproducible] is a reference system aligned with the rocket inertia axes. The unit vector [[??].sub.B] is orthogonal to the rocket plane of symmetry, while the two remaining unit vectors lie on it. The ABA frame is obtained from the WA frame through a sequence of two rotations: [mathematical expression not reproducible].

(8) The body axis (BA) [mathematical expression not reproducible] is, like the ABA frame, also attached to the rocket inertial axes only with a different orientation. Specifically, [[??].sub.B] is aligned with the rocket longitudinal axis [mathematical expression not reproducible] is directed downward so [mathematical expression not reproducible] and as consequence [mathematical expression not reproducible]. These relations lead to the following elementary rotation:

[mathematical expression not reproducible]. (14)

The knowledge of the instantaneous orientation of [mathematical expression not reproducible] is equivalent to identifying the instantaneous attitude of the rocket.

(9) The launch inertial frame (LI) [mathematical expression not reproducible] is an inertial system centered in the launch location at the launch instant. The unit vector [[??].sub.0] is aligned with the local north direction, [[??].sub.0] is directed toward the local east direction, and [[??].sub.0] is aligned with the downward direction. The LI frame is obtained from the ECI frame with a sequence of three rotations: the first two of them take the ECI frame to an intermediary reference frame, named auxiliary launch inertial frame (ALI) [mathematical expression not reproducible] using the geographical latitude and longitude of the launch base at the launch time, while the third one relates the axes of the ALI frame to the LI frame. The complete sequence is as follows:

[mathematical expression not reproducible]. (15)

The BA frame is also related to the LI frame through a sequence of three rotations using the attitude angles ([phi], [theta], [psi]). The complete sequence is [mathematical expression not reproducible].

In order to avoid ambiguities in the definition of the vectors of interest, the angles introduced in this section have to be constrained. Table 3 shows the limitation for the angles used in this work.

Any vector can be obtained with respect to any reference frame using the rotations defined in this section. These rotations can be used to define the rocket kinematics and dynamics through the motion, and Figures 6 and 7 portray all the reference frames used in this work, together with the related rotations.

4.2. Forces. The main term of the gravitational force G is given in the RV frame simply by

[mathematical expression not reproducible]. (16)

The components of the aerodynamic force can be referred to the WA frame and then projected in the RV frame with a simple rotation. The three components are the lift force, the drag force, and the side force.

[mathematical expression not reproducible], (17)

where

[mathematical expression not reproducible]. (18)

The aerodynamic forces have to be projected in the RV frame with a rotation about the bank angle around the second axis:

[mathematical expression not reproducible]. (19)

The aerodynamic forces are taken into account up to 100 km of altitude where the atmospheric density becomes negligible.

The thrust can be projected along the RV frame with the in-plane [[alpha].sub.T] and the out-of-plane [[beta].sub.T]:

[mathematical expression not reproducible], (20)

and to avoid ambiguities, the following constraints must hold:

-[pi] < [[alpha].sub.T] [less than or equal to] [pi], -[pi]/2 [less than or equal to] [[beta].sub.T] [less than or equal to] [pi]/2. (21)

The thrust vector T is considered aligned with the rocket longitudinal axis so two angles, [[alpha].sub.T] and [[beta].sub.T] suffices to describe its direction.

4.3. Equations of Motion. The equations of motion that govern the three-dimensional rocket dynamics can be conveniently written in terms of {r, [[lambda].sub.G], La, [[gamma].sub.R], [v.sub.R], [[zeta].sub.R]}. These variables refer to the relative motion in the RV frame. They form the state vector [x.sub.R] of the launch vehicle (in rotating coordinates).

[mathematical expression not reproducible]. (22)

The initial values for (22) are reported in Table 4.

Equation (22) is valid for both the propelled arc and the ballistic flight of the ogive. The values for La and [[lambda].sub.G] in Table 4 are referred to the Malindi launch base, while the launch is in the east direction ([[zeta].sub.R]). The control laws for ([[alpha].sub.T], [[beta].sub.T]), are determined using three distinct guidance scheme, namely, the constant attitude guidance, the near radial guidance, and the sun-pointing guidance, which are described in the following section. The state [x.sub.R] is continuos across the first stage separation, which occurs at time [t.sub.B]. After that, the payload continues its motion in a ballistic trajectory.

5. Guidance Algorithms

5.1. Constant Attitude Guidance. With this guidance, the rocket orientation does not change during the flight and was actually implemented in the Black Brant X (cf. [17]). This is ideal to achieve nominal (perfect) microgravity conditions ([DELTA]g = 0), due to rotation rates equal to 0. The attitude angles {[phi], [theta], [psi]} remain constant; therefore, in this type of guidance, the following expressions can be assumed as follows:

[phi] = [[phi].sub.0] [theta] = [[theta].sub.0] [psi] = [[psi].sub.0]. (23)

The remaining angles {[[alpha].sub.T], [[beta].sub.T], [alpha], [beta], [sigma]} can be found using analytical formulas at any times using the rotation matrices portrayed in Figure 6.

5.2. Near Radial Guidance. This guidance is called near radial because it is defined in the AO reference frame and the thrust direction is defined by

[[alpha].sub.T] = [[gamma].sub.R0] + t/[t.sub.B] ([pi]/2 - [[gamma].sub.R0]) - [[gamma].sub.R], [[beta].sub.T] = 0 (24)

whereas [[gamma].sub.R0] = [[theta].sub.0] is the pitch angle at the initial time and is equal to the initial flight path angle. With [[beta].sub.T] = 0, the angle that the rocket forms with the local horizontal axis is given by ([[alpha].sub.T] + [[gamma].sub.R]). At the initial time, this angle corresponds to the initial flight path angle, while at the burnout time, it corresponds to ([[alpha].sub.T] + [[gamma].sub.R]) = [pi]/2. The roll angle is can be chosen arbitrarily and is set to 0 deg. For the remaining angles {[theta], [psi], [alpha], [beta], [sigma]}, they are obtained through analytical expressions through the rotation matrices in Figure 6. The thrust direction tends to be aligned with the radial direction. Thus, in principle, this guidance leads to obtain the mean maximal altitude for the rocket at hand. Moreover, with this guidance, the usual values of [DELTA]g are in the order of [10.sup.-4] m/[sec.sup.2].

5.3. Sun-Pointing Guidance. The sun-pointing guidance is based on maintaining constant attitude of the rocket with respect to the sun position. The attitude angles {[phi], [theta], [psi]} are computed in an analytical fashion using the rotations in Figure 6. The first step to obtain the attitude of the sounding rocket is to locate the sun position with respect to the ECI frame:

[mathematical expression not reproducible], (25)

where [[lambda].sub.S] and [[xi].sub.S] are the sun declination and the right ascension, respectively. Then, projection of this unit vector into the ALI frame yields as follows:

[mathematical expression not reproducible]. (26)

So the thrust direction is obtained by aligning the body axis of the launcher with the sun direction in the ALI frame, as shown in Figure 8. To avoid that the rocket will fly with a [[gamma].sub.r] near to zero while trying to catch the sun slightly above the local horizon, the thrust direction is constrained to have a displacement from [[??].sub.[dot encircle]] not exceeding the maximal value [eta] = 10deg (cf. Figure 8). So the thrust direction is defined imposing [[xi].sub.T] = [[xi].sub.[dot encircle]] and [[lambda].sub.T] = max ([[lambda].sub.[dot encircle]], [pi]/2 - [eta]). The angles [[alpha].sub.T] and [[beta].sub.T] are then obtained through rotation of [??] from the ALI frame to the RV frame using the rotations reported in Figure 6 and considering that the thrust is aligned with the longitudinal axis of the rocket [[??].sub.B]. The pitch and yaw angles are then obtained with the following expressions:

[mathematical expression not reproducible]. (27)

As roll angle is arbitrary, it can be selected after choosing the angle [[iota].sub.S] related to the sun position in the body frame. The roll angle [phi] is obtained by comparing [[??].sub.[dot encircle]] in the ALI frame with its projection in the RV frame.

From Figure 8(c), it is straightforward to obtain the following angles:

[mathematical expression not reproducible]. (28)

They univocally define the vector [mathematical expression not reproducible] in the body frame. Comparing the coordinates of [[??].sub.[dot encircle]] in the ALI frame to the ones in the BA frame leads the following expression:

[mathematical expression not reproducible]. (29)

Equation (29) leads to the roll angle [phi].

Then, the aerodynamic angles [alpha], [beta], [sigma] can be obtained through the rotation matrices in Figure 6, which are not reported for the sake of brevity. In short, the sun-pointing guidance requires specifying the following parameters:

(i) [[theta].sub.S] which is the sun position in the ECI frame

(ii) [eta] which constraints the possible thrust direction

(iii) [[theta].sub.G0], defined by the date and hour of the launch

(iv) [l.sub.S], needed to obtain the roll angle [phi]

As the attitude is also in this case constant, the sun-pointing guidance is a particular case of the constant attitude guidance. Hence, it maintains its main advantages like the ideally perfect microgravity condition ([DELTA]g = 0) and can provide the perfect conditions for thermal experiments in the space.

6. Numerical Results

Figures 9-11 portray the results of the computation, which are, respectively, for the constant attitude guidance (Figure 9), the near radial guidance (Figure 10), and the sun-pointing guidance (Figure 11). From Figure 9 to Figures 11(b), 11(d), and 11(f)), if the desired [DELTA][T.sub.mg] is not chosen, it is possible to obtain a maximum [DELTA][T.sub.mg] for each length. This result is particularly visible for the shorter lengths while is less evident for the bigger configurations. Panels (a), (c), and (e) of Figures 9-11 show every (l, d) combination which guarantee 600 and 900 sec of [DELTA][T.sub.mg]. The optimal configuration corresponds to the point in which each curve is tangent to an equilateral isovolume hyperbola in the (l, [d.sup.2]) graph. As the tangence point always occurs in the lowest right part of the graph, two constraints were inserted, for the sake of avoiding too slender or bulky configurations. In general, the constraint for the l/d is

5 [less than or equal to] l/d [less than or equal to] 20. (30)

It should be noticed that that the optimal configurations for all the three guidance are nearly the same, and this means that the preliminary design of the sounding rocket can be done regardless of the guidance type. Tables 5 and 6 show the optimal rocket configuration for 600 and 900 sec of [DELTA] [T.sub.mg] and for all the payload masses. The configuration with l = 7.29m and d = 0.36m in Table 5 is particularly similar to the Black Brant VC, demonstrating that with similar [m.sub.U] and [DELTA][T.sub.mg] the algorithm can be effective to optimize the size of the rocket.

Figures 12-14 portray the complete state vector elements and other variables for the optimal configuration with [m.sub.U] = 50 kg and [DELTA][T.sub.mg] = 10 min, using the near radial guidance.

For the rocket structural integrity, the constraint Q[[alpha].sub.tot] < 150000Pa * deg must be satisfied, where

Q = 1/2 [rho][v.sup.2.sub.R] (31)

and

[[alpha].sub.tot] = arccos (cos [alpha] cos [beta]). (32)

From Figure 13(f), it is apparent that the constraint is satisfied, and the same occurs for all the cases.

7. Conclusions

Designing a sounding rocket is a difficult task, and this work proposes a method to simplify the early stage modeling process and performance evaluation. The new methodology at hand has four principal phases: (i) mass and geometry sizing based on reasonable assumptions, (ii) scaled propulsion modeling, (iii) accurate aerodynamic modeling, and (iv) the performance evaluation using different guidance schemes. With regard to the guidance algorithms, three schemes where presented which are the constant attitude, near radial, and sun-pointing guidances. While the constant attitude guidance is a well-known scheme, the other two are novel guidances. In particular, the benefits of the sun-pointing guidance are presented, which are the nearly perfect microgravity condition and the possibility to orient the rocket with a defined attitude with respect to the sun. The results clearly show that the algorithm is capable of providing a fast and reasonable sizing in accordance with the mission requirements. The results obtained can be used as a reference point for a montecarlo campaign to analyze the impact points in the presence of local winds or as a first guess solution for more sophisticated optimization algorithms.

Nomenclature g: Earth gravitational acceleration, m/[sec.sup.2] [[mu].sub.E]: Earth gravitational parameter, [km.sup.2]/[sec.sup.3] d: Rocket diameter, m [l.sub.ogive]: Ogive length, m [mathematical expression Dimensionless total initial mass not reproducible] [m.sub.0]: Initial mass, kg [m.sub.P] Propellant mass, kg [V.sub.body]: Rocket volume, [m.sup.3] u: Mass ratio T: Thrust, N [I.sub.SP]: Specific impulse, sec [t.sub.B]: Burnout time, sec [C.sub.L]: Lift coefficient [C.sub.Q]: Side force coefficient D: Drag force, N M: Mach number [alpha]: Angle of attack, deg S: Reference surface, [m.sup.2] [DELTA][T.sub.mg]: Microgravity duration, sec [[theta].sub.G]: Greenwich sidereal time, deg [L.sub.L]: Launch base latitude, deg [[lambda].sub.G]: Longitude, deg [v.sub.R]: Relative velocity m/[sec.sub.2] [sigma]: Bank angle, deg [theta]: Pitch angle, deg [[alpha].sub.T]: In-plane control angle, deg m: Mass, kg [[lambda].sub.[encircle]]: Sun declination ALI frame, deg [[lambda].sub.S: Sun declination ECI frame, deg [[lambda].sub.T: Thrust declination ALI frame, deg [L.sub.S]: Sun declination BA frame, deg [[omega].sub.E]: Earth rotation rate, deg/[sec.sub.2] l: Rocket total length, m t: Time, sec [l.sub.body]: Rocket body length, m [C.sub.ogive]: Dimensionless ogive length [m.sub.aux]: Auxiliary mass, kg [m.sub.S]: Structural mass, kg [epsilon]: Structural coefficient [m.sub.U]: Payload mass, kg [n.sub.0]: Acceleration at launch, g [??]: Mass rate, kg/sec [m.sub.0body]: Cylinder initial mass, kg [C.sub.D]: Drag coefficient L: Lift force, N Q: Side force, N Re: Reynolds number [beta]: Sideslip angle, deg [rho]: Atmospheric density, kg/[m.sup.3] r: Radius, m La: Latitude, deg [[lambda].sub.aL]: Launch base longitude, deg [[gamma].sub.R]: Relative flight path angle, deg [[zeta].sub.R]: Heading angle, deg [phi]: Roll angle, deg [psi]: Yaw angle, deg [[beta].sub.T]: Out-of-plane control angle, deg Q: Dynamic pressure, Pa [[xi].sub.[encircle]]: Sun right ascension ALI frame, deg [[xi].sub.S]: Sun right ascension ECI frame, deg [[xi].sub.T]: Thrust right ascension ALI frame, deg [[chi].sub.S]: Sun right ascension BA frame, deg.

https://doi.org/10.1155/2018/1678709

Conflicts of Interest

The authors declare that they have no conflicts of interest.

References

[1] A. Asensio Ramos and R. Manso Sainz, "Signal detection for spectroscopy and polarimetry," Astronomy and Astrophysics, vol. 547, p. A113, 2012.

[2] E. Luvsandamdin, S. Spiefiberger, M. Schiemangk et al., "Development of narrow linewidth, micro-integrated extended cavity diode lasers for quantum optics experiments in space," Applied Physics B, vol. 111, no. 2, pp. 255-260, 2013.

[3] S. Wieman, L. Didkovsky, T. Woods, A. Jones, and C. Moore, "Sounding rocket observations of active region soft X-ray spectra between 0.5 and 2.5 nm using a modified SDO/EVE instrument," Solar Physics, vol. 291, no. 12, pp. 3567-3582, 2016.

[4] Sounding Rocket Program Handbook, National Aeronautics and Space Administration, Goddard Space Flight Center, 2001.

[5] G. Seibert, The History of Sounding Rockets and Their Controibution to European Space Research, ESA Publications Division, 2006.

[6] S. M. Chowdhury, J. P. de la Beaujardiere, M. Brooks, and L. Roberts, "An integrated six degree-of-freedom trajectory simulator for hybrid sounding rockets," in 49th AIAA Aerospace Sciences Meeting including the New Horizons Forum and Aerospace Exposition, Orlando, Florida, January 2011.

[7] S. M. Chowdhury, Design and Performance Simulation of a Hybrid Sounding Rocket, Master's thesis, University of KwaZulu-Natal, 2012.

[8] L. Casalino and D. Pastrone, "Optimization of hybrid sounding rockets for hypersonic testing," in 45th AIAA/ASME/SAE/ASEE Joint Propulsion Conference & Exhibit, Denver, Colorado, August 2009.

[9] L. Casalino and D. Pastrone, "Optimization of hybrid sounding rockets for hypersonic testing," Journal of Propulsion and Power, vol. 28, no. 2, pp. 405-411, 2012.

[10] S. Nonaka, H. Ogawa, Y. Naruo, and Y. Inatani, "System design and tecnical demonstrations for reusable sounding rocket," in 20th ESA Symposium on European Rocket and Balloon Programmes and Related Research, Hyere, France, October 2011.

[11] D. L. Y. Woo and J. A. Martin, "Reusable sounding-rocket design," Journal of Spacecraft and Rockets, vol. 32, no. 2, pp. 376-378, 1995.

[12] F. Cremaschi, S. Weikert, A. Wiegand, W. Jung, and F. Scheuerpflug, "Sounding rocket trajectory simulation and optimization with astos," in 19th ESA Symposium on European Rocket and Balloon Programmes and Related Research, 2009.

[13] W. B. Blake, Missile DATCOM User's Manual--1997 Fortran 90 Revision, United States Air Force, 1998.

[14] L. Mangiacasale, Meccanica del volo atmosferico. Terne di riferimento, equazioni di moto, linearizzazione, stabilita, Ingegneria, 2000.

[15] M. Pallone, M. Pontani, and P. Teofilatto, "Accurate modeling and near optimal ascent trajectory of microsatellite launch vehicles via firework algorithm," in 67th International Astro-nautical Congress, 2016.

[16] M. Pontani, "Particle swarm optimization of ascent trajectories of multistage launch vehicles," Acta Astronautica, vol. 94, no. 2, pp. 852-864, 2014.

[17] L. Ljunge and L. Hall, "S19 guidance of the Black Brant X sounding rocket," Journal of Guidance, Control, and Dynamics, vol. 7, no. 2, pp. 156-160, 1984.

Marco Pallone (iD), Mauro Pontani (iD), Paolo Teofilatto, and Angelo Minotti (iD)

University of Rome "La Sapienza", Via Salaria 851/831, 00138 Rome, Italy

Correspondence should be addressed to Marco Pallone; pallone.1420138@studenti.uniroma1.it

Received 24 March 2017; Accepted 25 September 2017; Published 6 February 2018

Academic Editor: Corin Segal

Caption: Figure 1: Black Brant VC.

Caption: Figure 2: Main surfaces and reference quantities.

Caption: Figure 3: Example of scaled [m.sub.0] and [m.sub.P] with m[U.sub.] = 50kg.

Caption: Figure 4: Example of scaled T and [t.sub.B] with [m.sub.U] = 50 kg.

Caption: Figure 5: Example of the drag force coefficients provided by the DATCOM software for the configuration l = 6 m and d = 0.8 m.

Caption: Figure 6: Reference frames and related rotation matrices.

Caption: Figure 7: Main reference frames and related angles.

Caption: Figure 8: Sun-pointing guidance principal angles and axes.

Caption: Figure 9: Optimal results for the constant attitude guidance.

Caption: Figure 10: Optimal results for the near radial guidance.

Caption: Figure 11: Optimal results for the sun-pointing guidance.

Caption: Figure 12: State variables for the configuration with l = 7.29 m and d = 0.36 m (near radial guidance).

Caption: Figure 13: Control law, attitude angles, and dynamic pressure in the thrust arc for the configuration with l = 7.29 m and d = 0.36 m (near radial guidance).

Caption: Figure 14: Aerodynamic angles and coefficients in the thrust arc for the configuration with l = 7.29 m and d = 0.36 m (near radial guidance).

Table 1: "Maxus" and "Black Brant VC" data. BBVC Maxus Length (m) 8.3 16.2 Max diameter (m) 0.44 1 Mass liftoff (kg) 1544 12298 Propellant mass (kg) 1017 10042 Av. thrust (N) 69374 372500 Burnout time (sec) 32 63 Table 2: Auxiliary mass [m.sub.aux] of the ogive. Parachute mass (kg) 10 Heat shield mass (kg) 35 Metal case mass (kg) 5 Total [m.sub.aux] (kg) 50 Table 3: Angle constraints. -[pi] < [[theta].sub.G] [less -[pi] < [[lambda].sub.G] [less than or equal to] [pi] than or equal to] [pi] [pi]/2 < La [less than or -[pi] < [[zeta].sub.R] [less equal to] [pi]/2 than or equal to] [pi] [pi]/2 < [[gamma].sub.R] -[pi] < [sigma] [less than or [less than or equal to] [pi]/2 equal to] [pi] [pi]/2 < [beta] [less than or -[pi] < [alpha] [less than or equal to] [pi]/2 equal to] [pi] -[pi] < [psi] [less than or -[pi]/2 < [theta] [less than equal to] [pi] or equal to] [pi]/2 -[pi] < [phi] [less than or -- equal to] [pi] Table 4: Initial condition for (22). r (km) [R.sub.E] + 0.001 [v.sub.R] (km/s) 0.008 [[gamma].sub.R] (deg) 86 La (deg) -2.938 [[zeta].sub.R] (deg) 0 [[lambda].sub.G] (deg) 40.213 Table 5: Optimal configurations for each payload mass, [DELTA][T.sub.mg] = 10 min. [m.sub.U] (kg) 50 100 200 l (m) 7.29 8.25 9.69 d (m) 0.36 0.41 0.48 Table 6: Optimal configurations for each payload mass, [DELTA][T.sub.mg] = 15 min. [m.sub.U] (kg) 50 100 200 l (m) 10.98 12.42 14.53 d (m) 0.55 0.62 0.73

Printer friendly Cite/link Email Feedback | |

Title Annotation: | Research Article |
---|---|

Author: | Mauro, Marco Pallone; Pontani; Teofilatto, Paolo; Minotti, Angelo |

Publication: | International Journal of Aerospace Engineering |

Article Type: | Report |

Geographic Code: | 1USA |

Date: | Jan 1, 2018 |

Words: | 5978 |

Previous Article: | The Impact of Concrete Pavement Field Floor to Vehicle Missile Launching Process. |

Next Article: | Dynamics Analysis of Deployable Structures considering a Two-Dimensional Coupled Thermo-Structural Effect. |

Topics: |