Research Article  Open Access
Thermal Protection System Design of a Reusable Launch Vehicle Using Integral Soft Objects
Abstract
In the present paper, a modelling procedure of the thermal protection system designed for a conceptual Reusable Launch Vehicle is presented. A special parametric model, featuring a scalar field irradiated by a set of bidimensional soft objects, is developed and used to assign an almost arbitrary distribution of insulating materials over the vehicle surface. The model fully exploits the autoblending capability of soft objects and allows a rational distribution of thermal coating materials using a limited number of parameters. Applications to different conceptual vehicle configurations of an assigned thickness map, and material layout show the flexibility of the model. The model is finally integrated in the framework of a multidisciplinary analysis to perform a trajectorybased TPS sizing, subjected to fixed thermal constraints.
1. Introduction
In the last decade, a growing number of projects have been focused on the development of fully Reusable Launch Vehicles (RLV), designed for a crewrescue mission. Successful demonstrative flights of private companies, like SpaceX, Virgin Galactic, and Sierra Nevada Corporation, and the activity promoted by the European Space Agency are aimed at improving operability of RLV [1, 2]. Consequently, a great deal of research effort has been put to design RLV as blended wingbodies, because of the promising tradeoff between aerodynamic efficiency, crossrange, and aeroheating performances during the reentry [3]. The EXPERT (European eXPErimental Reentry Testbed) program and the IXV (Intermediate eXperimental Vehicle), which performed an atmospheric lifting reentry from orbital speed, are just examples of such demonstrators developed to predict performances of a fullscale vehicle. Besides, the X37B, an unmanned lifting body developed by Boeing, has been put into orbit by an Atlas V rocket, performing a successful liftingguided reentry. Furthermore, the foreseeable opportunity for space tourism represented by experimental flights of Virgin Galactic’s SpaceShipTwo and SpaceX has also emerged [4]. The requirements currently considered for RLV design are (i) to allow a very lowg (nearly 1.5 g) reentry, with a landing on a conventional runway; (ii) to adopt a lightweight (passive), fully reusable thermal protection system to withstand several flights without any replacement; and (iii) to provide vehicle autonomy to land at predefined locations for rescue issues [2, 3]. In order to fulfill all those requirements, the duration of reentry flight increases and consequently the integrated heat load absorbed by the structure [2].
The above circumstances may conflict with the adoption of a fully reusable TPS, eventually restricting the choice of insulating materials and penalizing the total mass [5]. However, antithetical requirements between room for the payload, weight, and vehicle operability demand a tradeoff between vehicle shape and TPS sizing. In preliminary design practice, thousands of design configurations are typically evaluated by an optimization algorithm to find the best fit [5–10]. As a consequence, a preliminary appraisal of vehicle performances is commonly performed using highefficiency, loworder fidelity methods that give support to a multidisciplinary analysis performed with a computational effort which fit the typical timeline of the conceptual design phase [11].
The aerothermal environment is a basic design criterion for either TPS sizing or choice of materials [12, 13]. TPS sizing is generally performed once the reentry trajectory is assigned, having computed the peak heating flux and the time integrated heat load [14]. External rocket propulsion systems allow an RLV a less severe heating due to their different ascent trajectories.
Therefore, an RLV operates an unpowered crew rescue with the TPS sized mainly by its aerothermal reentry corridor. On the other hand, an RLV that performs an ascent phase with airbreathing propulsion systems is sized considering more severe heating conditions [15]. Several works dealing with TPS sizing have been published in the literature. Lobbia [8] determined the sizing of a TPS in the framework of a multidisciplinary optimization. Material densities and maximum reuse temperature were computed. TPS mass was estimated assuming the category of materials used for the space shuttle and thickness distribution assigned on a review of HL20 materials for each component. Trajectorybased TPS sizing has been proposed by Olynick [13] for a winged vehicle concept. The peak heating temperature was determined considering an X33 trajectory, discretized in a number of fixed waypoints. The resulting aerothermal database was used as an input for a onedimensional conduction analysis, and several onedimensional stackups of different materials representative of TPS were consequently dimensioned. Bradford [15] developed an engineering software tool for aeroheating analysis and TPS sizing. The tool is applicable in the conceptual design phase, for reusable nonablative thermal protection systems. The thermal model was based on a onedimensional analysis, and TPS was modeled considering a stackup of ten different material layers. Mazzaracchio [14] proposed a method to identify the optimal combination between the ablative and the reusable part of TPS. This was dictated by a tradeoff between the reentry trajectory and the TPS sizing. In Wurster and Stone [16], a selection of design trajectories, and wind tunnel data were integrated with CFD simulations and engineering prediction to perform TPS sizing for an HL20 reentry vehicle. In the present work, following an approach proposed in Ref. [17], a soft objectderived representation for TPS thickness and material assignation is introduced. According to the legacy formulation of this technique, originally developed in computer graphics (C.G.) for the rendering of complex organic shapes [18], threedimensional object surfaces are (implicitly) obtained by defining a set of source points (or even more complex varieties) irradiating a potential field that is subsequently rendered as a chosen isosurface. A major feature offered by this methodology is that the topology of the regions represented by these objects can be arbitrary and easily controlled just by altering the mutual spatial position of the primitives. These very desirable properties have made implicit modelling techniques extremely popular in the field of photorealistic rendering of organic shapes whose topology changes dynamically, and many advances and developments have been proposed over time to improve their capabilities and/or reduce some unwanted sideeffects [19, 20]. In the present paper, following a quite different paradigm developed in [21], the full potential field irradiated by a set of bidimensional soft objects is congruently mapped on a discretized paneled RLV shape. The methodology is able to create arbitrary TPS distributions seamlessly increasing the thickness where critical heat loads are experienced and dropping out elsewhere. A similar, slightly modified procedure can also be applied to create an arbitrary binary map representing membership functions for TPS materials. This binary map can be operated independently on thickness distribution (or locally synchronized with different thickness maps). The present formulation is formalized in the framework of a parametric model which exploits simple variations of parameters to perform the soft object mapping over a discretized surface. Applications of the developed procedure to different vehicle configurations show the flexibility of the method. In addition, the procedure can be easily embedded in a multidisciplinary computation, to perform a trajectorybased TPS sizing on a RLV subjected to fixed thermal constraint.
2. Methodology
2.1. Soft Object Definition
Soft objects constitute a modelling technique introduced in computer graphics to represent threedimensional objects having complex and organic shapes. Blinn first developed a soft object model to display the appearance of electron density clouds in a covalent bond of a molecule [18]. According to the model formulation, curved (closed) surfaces can be modeled defining potential fields , namely, blobs. Several blobs can be connected smoothly by the selfblending property, by performing an algebraic summation of their potential fields [22]:
The commonly adopted notation separates the distance field from the field function . The strength accounts for the value of in a generic point , at distance from a key point .
The potential field generated by blobs can be generally taken into account computing a specific isosurface (subsequently processed) by a raster conversion algorithm: where the threshold in equation (4) selects an isosurface of . Blinn originally proposed the “blobby molecule,” an isotropically decaying Gaussian function modulated by strength, and radius [18]. The blobby molecule is a potential function with infinite support. This aspect affects the computational effort in a practical implementation, because it has to be evaluated in all points of the space. However, in the literature, several finite support potential functions have been proposed for different modelling purposes [22]. The field function used in the present work has a finite support and assumes normalized values in the range between 0 and 1 [22]:
The parameter defines the hardness of blob and controls the level of blending between two soft objects.
2.2. TwoDimensional Integral Soft Object Field
Twodimensional soft objects preserve the selfblending property. However, it is not always easy to create a rational distribution of a set of independent point source objects for application purposes. Therefore, blobs are conveniently and easily arranged in macroaggregates. Figure 1 shows that a potential field is created superposing discrete blobs having radius and centers mutually placed on a line segment of length over an equally spaced twodimensional grid of step . If two or more blobs superpose, the strength of the potential field is obtained summing up the strengths of each blob (yellow colored region). This procedure relies on a similar idea to the one developed in [21] to generate selfstiffened structural panels. Specifically, the full (integral) potential field irradiated by a set of discrete twodimensional point source blobs generates a seamless potential field. This approach is quite different from isocontour tracking commonly adopted to represent soft objects.
(a)
(b)
2.3. Modelling of TwoDimensional Stick Primitives
A superposition of point source blobs with key points placed on a geometric segment (straight or curved) is denoted from now on as a “stick.” However, as shown in Figure 1(b), equation (1) creates a stick with a support having “bulges.” Increasing the number of sticks, the shape of the support becomes more regular. However, the strength is not bounded (see Figure 1(b)). The above drawback is overcome by modifying the definition of the potential field given by equation (1). A bounded potential field, regardless of the number of the blobs used on a stick, is obtained with the relation
Equation (6), where , expresses the global potential field irradiated by a set of blobs at a generic point of space placed at a distance from the key points, as the max between the previous potentials accounted for by the assembly layer and the current potential over the plane disk of radius : where is the total number of blobs present on the Bgrid. Figures 2(a) and 2(b) show the strength field of a twodimensional stick primitive obtained using and 20, respectively, computed with equation (6). By increasing the number of blobs on a stick, the strength of is still bounded to a maximum unit value. Figures 2(c) and 2(d) show the same behavior for a tapered primitive having a linear variation of the blob radius along the axis of sticks.
(a)
(b)
(c)
(d)
3. RLV Shape Modelling and Thermal Protection System Sizing Criteria
In the present work, we assume that the generic shape of an RLV is represented by a grid formed by a quadrangular and/or by a degenerated triangular panel grid. Grid points are obtained using a proprietary procedure that authors fully detailed in [23, 24]. Without going into details of the shape model, we remark that the mesh arrangement over the RLV surface is obtained with no NURBS support surface: a threedimensional parametric wireframe is created using cubic rational Bsplines and used to reconstruct the computational surface grid. The control parameter allows a wide range of shape variations to handle different design objectives (thermal or dynamical) for a reentry mission. Grid topology is equivalent to a spherical surface with no singularities (open poles) and allows a mapping of the points in UV coordinates over an equivalent cylindrical surface. The above considerations ensure a topologically invariant shape. In previous papers proposed by the authors [23, 24], a multidisciplinary shape optimization for an RLV comprising a trajectorybased TPS sizing procedure was developed. The TPS was modeled using two insulating materials placed at different locations along the vehicle surface. A different mapping (thickness distribution and longitudinal location) of the two materials with different operational temperature was adopted. The sizing of insulating materials required the computation of aerothermal loads across the reentry trajectory from a LEO orbit up to , which is considered the limit below which thermal heating can be neglected. TPS thickness was parametrically sized according to thermal requirements assumed in the optimization; a simple (but very rough), bilinear distribution of the TPS thickness along the longitudinal axis of the vehicle and a linear distribution across each cross section were adopted. The maximum allowable temperature values (depending on the adopted material) for the interior and exterior surfaces of TPS outlined the thermal constraint to be fulfilled by the sizing procedure.
4. Soft Object Design of TPS
4.1. Rationale
The modelling procedure for the TPS is defined starting from the definition of a set of soft objects which are represented on the topological map associated with the current morphology of the object, as shown in Figure 3. Consequently, the supports of the sticks are adjusted according to the normalized dimensions relative to this map. The topological map is emulated introducing a twodimensional grid (from now on denoted as Bgrid) having the same topology tree as the vehicle open grid (number of points, panels, and connectivity) except for the unit size. A geometric mapping between the Bgrid and the vehicle grid is established, and elements of the Bgrid are biunivocally mapped onto corresponding elements of the vehicle surface (see Figure 3). Several stick primitives are emulated on the Bgrid placing a number of equally spaced isotropic blobs, with radius and length , in normalized units. Stick emulation is performed by overlapping blobs using the special formulation reported in [21] that ensures a convergent envelope of the finite support and a limited value of the blob strength. An exemplificative spatial distribution of sticks on the Bgrid is shown in Figure 3. Position and orientation of each stick is determined by assigning coordinates of centers and precession angles , respectively, with respect to a Cartesian frame of reference Oxz oriented as in Figure 3. Therefore, a generic distribution of sticks created on the vehicle grid is equally mapped on the vehicle surface whatever the morphological map considered. In the present case, gray colored regions (1) denote points of the Bgrid mapped on the windward side of the RLV shape (see Figure 3), while white regions (2) relate to leeward regions of the vehicle. Regions of the vehicle surface mainly subjected to heating peaks during the reentry maneuver are the (i) nose, (ii) leading edge, and (iii) tail. The global potential field generated by the sticks onto the Bgrid is adjusted in a suitable dimensional scale and subsequently mapped on the mesh panels of the vehicle surface grid to obtain an easy and powerful control of the thickness distribution. The proposed methodology is able to create virtually arbitrary TPS distributions and can be easily tuned up to locally increase the thickness where critical heat loads are expected and dropping out elsewhere. A similar, slightly modified procedure is also applied to create an arbitrary binary map distribution of different TPS materials that may be operated independently of the thickness distribution.
5. Parametric Model of a Thermal Protection System
5.1. Thickness Modelling
As a demonstrative example, a parametric representation of a thermal protection system is obtained using a limited set of stick primitives (), oriented as shown in Figure 4. Skin sticks characterized by a large radius and limited strength are spread over the skin surface in the longitudinal direction in order to provide a thickness graded baseline. A constant minimum thickness is superposed in all remaining points of the Bgrid, ensuring a nonzero value in any point of the grid. Furthermore, additional parametric sticks, specifically positioned and oriented to affect thickness in critical regions, such as the nose, leading edge, and trailing edge, complete the support for the TPS and create a rational distribution of insulating material suitable with a reentry mission. The parametric position of sticks and the axis of orientation are defined by assigning centroid coordinates and and angle θ_{th}, measured with respect to the system of reference reported in Figure 4. Length () and strength () are expressed with the parametric relations:
Skin () and nose () sticks have a tapered support obtained by imposing a linear variation of the point source blob radius. Conversely, a constant radius is adopted for the leading edge () and trailing edge () sticks.
5.2. Material Modelling
A similar but completely independent stickbased parameterization has been also defined to model a dynamic distribution map of different insulating materials, denoted here generically as material “0” and material “1.” We assume that material “1” outperforms material “0.” Therefore, material “1” is a natural candidate to insulate the nose, leading edge, and trailing edge. TPS materials are assigned according to a discrete distribution. Differently than sticks used for thickness distribution, this additional set of primitives returns just binary values used to define specific materials. In this case, the field function (see relations in (9)) assumes a constant value equal to one inside the finite support of a stick and zero elsewhere. Centroid coordinates ( and ) and length of these additional sticks are given by parametric relations in (9) with normalized parameters reported in Table 1:

6. Additional Considerations about Integral Soft Objects for TPS Modelling
In order to better clarify the rationale underlying the proposed methodology, some additional considerations are provided next. As a general premise, it should be emphasized that implicit modelling techniques canonically used in C.G. differ in many respects from those reformulated and used in the present context. About this, a brief description of the approaches followed in C.G. is preliminarily given, to highlight both common points and main differences with respect to the illustrated methods. In C.G., implicit modelling is commonly used for the rendering of complex organic shapes. In these methods, some objects (usually referred to as blobs or primitives) of appropriate dimensionality (2D, 3D), typically represented by their own morphological skeletons, are conceived as “emitters” of suitable finite support field functions, expressed as distance laws in an appropriate norm (usually Euclidean). These primitives are allowed to mutually interact with each other by simply overlapping their finite supports, cumulating in that way the field intensities where superposition takes place. According to an implicit representation, a specific instance (isosurface) of the global field associated to an assigned isovalue can be finally visualized. Depending on the chosen rendering method (e.g., raytracing), the isosurface can preserve its implicit formulation or be evaluated through suitable progressive sampling algorithms (i.e., octree) to be translated into discrete polygonal elements, typically triangular meshes. In general, the mathematical structure of the isosurfaces will be characterized by smooth curvatures, highlighting an intrinsic capability to generate automatic fillets in those spatial regions where primitives overlap. Fillets can be controlled locally by introducing a hardness parameter in the field functions, which makes the primitives “harder” or “softer” in blending. As previously stated, in the present work, a methodology just roughly inspired by the aforementioned techniques has been introduced to arbitrarily control thickness and material assignment using a limited number of control parameters. The proposed approach uses the same finite support field functions employed in implicit modelling; therefore, the desirable capability to generate arbitrary topologies by simply controlling the mutual spatial position of primitives is maintained. However, differently from what is being implemented in C.G., where the discrete resolution of the tracked isosurface is often adaptive and conditioned by the expected rendering quality and the amount of local curvatures, these primitives are now just projected onto a twodimensional flat grid with a fixed resolution that acts as an extended finite support; the global field resulting from primitive interaction is then used integrally, not just represented with single contours. Therefore, no implicit formulation is used at all. These circumstances explain the definition of “integral soft objects” (opposed to “implicit modelling”) used in the present context to describe the “agents” that operate to assemble the maps of thickness and material. In this new framework, the legacy control parameters for field functions adopted in C.G. are still maintained, although with a different semantic connotation. Specifically, the strength parameter has been reinterpreted either as the maximum thickness value of TPS locally transferred onto the topological map or the maximum integer value of the set of pointers; instead, the hardness parameter becomes the thickness gradient with respect to the distance in the finite support. A major aspect that affects the implementation of the method is given by the low spatial resolution of the topological map where TPS thickness and materials are transferred. Due to the specific methodologies (panel methods described next) used in the resolution of the thermal and fluid dynamic fields on the body at different speed regimes, this grid is extremely coarse if compared to the typical resolutions used in C.G. (only a few hundred quadrangular elements are used). In addition, the resulting thickness and material distributions are sampled only at the centroids of the panels and assumed constant for each panel. These intrinsic limitations have determined the above described implementation choices; those rationales are briefly given here: (i) a stick primitive has been preferred and systematically adopted in this context, because effective TPS thickness distributions likely happen along the UV parametric directions of the topological map. Effectiveness and flexibility of these distributions have been increased by also allowing tapered and steerable sticks. (ii) A surrogate representation of a skeletonbased stick is simply obtained by distributing a finite number of onepoint skeleton primitives (circular blobs), mixing each other with the MAX blending operator outlined by equation (6). As mentioned earlier, the rationale of this choice is that the field function of the assembled stick primitive still remains limited, and at the same time, only a small number of circular blobs (slightly more than a dozen) are required to obtain a satisfactory envelope of the finite support without bulges because of the low intrinsic resolution of the grid combined with the actual ranges assigned to the finite supports of primitives. Moreover, this simplified approach allows in perspective the definition of some other primitives based on different skeletons (homeomorphic with the line segment, for example, based on splines) that might potentially ensure more sensitivity in the TPS definition problem with no major modification required. (iii) Although it is virtually possible to use any suitable blending function to merge the fields of different preassembled primitives (for example, the sum of the local magnitudes), the MAX blending function is also still preferred for this purpose, because the intrinsic capability offered by this operator to generate “zigguratstyle” step fields when the involved primitives exhibit different strength values has been considered desirable for an effective TPS thickness distribution.
7. An Example of TPS Modelling Capabilities
The previously introduced modelling procedure has been applied on a conceptual RLV shape created with the model described in Section 4 and detailed in [23, 24]. The applicability of the procedure is shown for the arbitrarily chosen distribution of stick primitives that creates a morphologically adaptive TPS on two RLV shapes with different dimensions: RLV1 with length m, wingspan m, and cabin height m and RLV2 with length m, wingspan m, and cabin height m.
The parameters characterizing the distribution of thickness and the materials are reported in Table 1. Figures 5(a) and 5(b) show the application of TPS modelling over the first configuration (RLV1), on leeward (a) and windward (b) surfaces. Different colors denote different values of thickness and are represented in a dimensional scale. It can be observed that the thickness map can be easily tuned up for best covering of regions where maximum heat loads occurs (i.e., the nose and leading edge). Figure 5 shows the capability to create arbitrary seamless thickness distributions up to the value of the baseline thickness which has been arbitrarily set equal to m (denoted in blue color). This corresponds to a region of the leeward surface not covered by the skin stick. Figures 5(c) and 5(d) show the map of two different insulating materials created with equation (9). Red colors indicate material 1, which is placed on regions of the vehicle subjected to higher heat loads. Comparisons between Figures 6(a) and 6(b) and Figures 6(c) and 6(d) also exhibit the capability of the model to handle independently both the thickness and material distribution. Finally, Figures 6(a) and 6(b) and Figures 6(c) and 6(d) show the same blob distribution adopted either for thickness or for material modelling applied on a different RLV configuration (RLV2). The procedure creates, as expected, the same TPS distribution both for thickness or materials on two different shapes and is completely independent by their morphology.
(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)
8. Application of the TPS Sizing Procedure for a Conceptual RLV Configuration
The TPS modelling procedure developed in the previous sections has been implemented in the ANSYS^{®} Parametric Design Language [25], to perform a trajectorybased sizing of a TPS for an RLV designed for a LEO reentry mission. A multidisciplinary analysis comprising aerodynamics, heating analysis, trajectory estimation, and mass estimation is implemented to determine the aerothermal loads on the vehicle. The entire flowchart of the procedure has been discussed and detailed in [23, 24]. For the sake of brevity, here, it will be addressed with reference to the specific application, specifying the assumptions adopted. TPS is designed to withstand heating for an unpowered reentry maneuver performed from an altitude of km and down to . The reentry is finalized with a conventional landing performed at the prescribed speed. The mission is based on the following stages during reentry: (1) hypersonic phase, (2) supersonic phase, (3) subsonic phase, and (4) landing. As we want to address the applicability of the developed sizing methodology for a conceptual design configuration, loworder fidelity methods [3, 5] are used for each subdiscipline to reduce the computational time [11].
8.1. Aerodynamics
Aerodynamic coefficients for the hypersonic phase of the reentry are computed with a public domain panel flow solver using Newtonian flow theory available in Ref. [26], in the waypoints reported in Table 2. According to Newtonian theory, hypersonic flows are modeled as an ensemble of particles impacting the surface of a body approximated to a flat plate at incidence. Therefore, considering a discrete panelization of the body surface, the pressure force acting on the th panel having an area and directed along the outward unit normal vector to the panel is then computed as . The nondimensional force coefficients in a body frame of reference are given by computed as the summation of pressure forces over the total number of mesh panels. Newtonian approximation is acceptable in the range of altitudes and Mach numbers up to , as viscousforce contribution in the axial direction on hypersonic bodies decreases with an increasing angle of attack [3, 8, 27]. Validation of the present computations on aerodynamic coefficients has been previously performed in [24], together verifying the mesh independence of aerodynamic coefficients. Different choices of impact methods are adopted in the computation and are indicated in Table 3. The rationale behind the partitioning of the nose and wingbody region is related to the different flow incidence occurring on the wingbody panels requiring the use of appropriate compression methods adopted along the windward region (see Table 3). A similarity between the flow on the wingbody panels and the flow downstream of an attached shock can be assumed, and the tangent cone method can be used [8]. Aerodynamic coefficients of the vehicle at an incompressible Mach number are computed by adopting the integral formulation of the potential generated at a point by a distribution of singularities (sources and doublets) with strength and , respectively, being and being the body and the wake surface, respectively. The freely distributed panel code available in Ref. [28] is adopted for the current computation of this aerodynamic regime. The solver requires a discretization of the geometry by quadrilateral panels with sources and doublets of constant strength. An approximation of equation (11) on the computational mesh is


The unknown values of doublet strength are obtained by solving the linear system of equations obtained by discretizing equation (11). To make the solution uniquely defined, tangency condition is applied, i.e., , and the source strength is set equal to on each control node. Furthermore, a wake surface, still based on quadrilateral panels, has been added to apply the Kutta condition on each panel at the trailing edge. The intensity of the wake panels is set equal to the difference between the doublet strength of the upper and lower doublet panels at the trailing edge. Having computed the doublet strengths, the determination of velocity, pressure, and therefore aerodynamic coefficients is performed.
8.2. Mass Estimation
The mass of the thermal protection system has been estimated modelling the TPS with two of the materials adopted for the space shuttle thermal insulation system whose thermal properties are reported in [29, 30]: the Reinforced CarbonCarbon (RCC) composite and HighTemperature Reusable Surface Insulation (HRSI) tiles, made of coated LI900 silica ceramics. RCC material can be selected for the thermal insulation of the nose, leading edge, or trailing edge, where the peak wall temperature is expected. A mass decoupling model, based on the set of relations reported in [31], is adopted in the current procedure to compute the mass of the vehicle: being excluding the contribution of the thermal protection system, where , , , , and are reported in Table 4. The mass of the thermal protection system , with , , and being the density of the insulating material, area of the discrete th panel of TPS, and thickness of the th TPS panel, respectively, computed with the procedure developed in Section 5.

8.3. Flight Dynamics
The calculation of the trajectory is performed in a nonrotating, inertial, EarthCentered, EarthFixed (ECEF) frame of reference, according to the general system of equation of planetary flight reported in Ref. [32]. The vehicle is assumed as a mass point, describing a nonplanar reentry trajectory with a constant bank angle , assigned to ensure a crossrange performance. A negative value of initial reentry flightpath angle is given to avoid the skip phenomenon [32]. The angle of attack changes following an implicitly defined modulation law, reported in Table 2. The initial value of longitude , latitude , and flight azimuth are all set to zero. The vehicle dynamics is described by a fourdegreeoffreedom point mass model. Trajectory equations have been integrated by using an implicit NewtonRaphson method to reduce the computational time of the overall procedure, performing a time step sensitivity analysis to ensure convergence of the numerical integration.
8.4. Heating Analysis
The thermal state of the TPS is efficiently determined computing its exterior () and interior wall temperatures () which determines the choice of the material. The thermal analysis is only performed on the windward side as this is the maximum heated region of the TPS; the leeward side of the vehicle is supposed to be the fixed wall temperature region. Heat radiated from the external TPS wall provides the cooling of the vehicle. Catalytic recombination, low density effects, and thermal radiation from nonconvex surfaces are neglected. The above simplifying assumptions are considered reasonable for a conceptual design phase and determine a conservative overestimation of the TPS mass. A onedimensional model of the TPS thickness [6] is used to determine TPS interior wall temperatures exploiting the exterior TPS surface temperature. Radiative equilibrium is assumed on wall [5]: where is the surface emissivity and is the StefanBoltzmann constant. The convective heating at the wall at the stagnation point is approximated using the following correlation: where is the specific heat at the wall, the freestream density, and the vehicle velocity. The windward convective heating on the vehicle is evaluated using a correlation for spheres, cylinders, and flat plates [33]: where the constant is specifically characterized for a laminar or turbulent boundary layer. TPS wall temperatures are determined using kinematic trajectory data, in the range where peak values were expected to occur. Equation (18) is solved using a NewtonRaphson method to determine the wall temperature on each panel of the TPS thickness discretized as shown Figure 7. The interior TPS wall temperature is computed at each node of the discrete th panel of the TPS integrating in time a onedimensional unsteady heatdiffusion model [12, 23]. Both the initial and boundary conditions assigned as provide a wellposed heatdiffusion problem, numerically solved with a finite difference method. The thermal state of TPS is globally defined assuming, for each computational element , the following scalar relations: where the variable vanishes for negative values of the difference where it has been supposed that K is the value of the maximum allowable TPS temperature which adheres to structural elements of the vehicle and is the reuse temperature limit depending on the material [30].
8.5. Results
The multidisciplinary analysis previously described has been applied to perform sizing and material distribution of a TPS for an RLV with dimensions chosen equal to length m, wingspan m, and cabin height m. Thickness distribution of Figures 8(a) and 8(b) is obtained by setting .
(a)
(b)
(c)
(d)
Position, axis of orientation, and strength of stick primitives are manually and iteratively assigned, rationally choosing stick parameters, to cover the most heated regions of the TPS surface. A minimum value of the thickness set equal to m allows the convergence of the heat conduction problem (see equation (19)). A maximum value of thickness is chosen equal to m and takes into account the conservative assumptions made in Section 8.5; therefore, it is expected to be oversized but in agreement with the approximations used in the conceptual design phase [5]. It is observed in Figures 9(a) and 9(b) that the seamless distribution of thickness is obtained along the wingspan directions. The nose thickness is increased only on the windward side, because the leeward side is assumed to be at constant temperature during the reentry. Figures 8(c) and 8(d) show the binary material distribution between the RCC and Li900 tiles. The selected configuration fulfills the thermal constraint K on each panel of the discrete interior TPS wall as shown in Figures 9(a) and 9(b), and the computed mass of the thermal protection system is kg.
(a)
(b)
9. Conclusions
In the present paper, a special modelling procedure of the thermal protection system designed for a conceptual Reusable Launch Vehicle has been developed. A set of macroaggregates of point source blobs organized in envelopes of finite supports and with a bounded strength has been successfully created on the topological map associated with the computational grid. Applications of the modelling procedure to different design configurations highlighted the sensitivity and powerful control to radically change the TPS using a limited number of parameters. This feature has been illustrated executing a multidisciplinary analysis regarding a conceptual configuration where a simple iterative and manual sizing of the thermal protection system has been successfully performed, without adopting any optimization procedure.
Data Availability
The data supporting the analysis are from previously reported studies and datasets, which have been cited in Ref. [20, 21]. Specifically, the multidisciplinary analysis discussed in Section 7 has been performed on an RLV shape formulating a multiobjective optimization problem. The problem setup, including the range of design variables, has been detailed in [21]. Thermal properties of materials used in the TPS sizing procedure are available in [27].
Conflicts of Interest
The authors declare that there is no conflict of interest regarding the publication of this paper.
Acknowledgments
This work was supported by Università degli Studi della Campania Luigi Vanvitelli.
References
 D. Sziroczak and H. Smith, “A review of design issues specific to hypersonic flight vehicles,” Progress in Aerospace Sciences, vol. 84, pp. 1–28, 2016. View at: Publisher Site  Google Scholar
 E. Hirschel and C. Weiland, Selected Aerothermodynamic Design Problems of Hypersonic Flight Vehicles, SpringerVerlag/AIAA, 2009. View at: Publisher Site
 D. Dirkx and E. Mooij, Conceptual Shape Optimization of Entry Vehicles Applied to Capsules and Winged Fuselage Vehicles, Springer Aerospace Technology, Springer International Publishing Switzerland, 2017. View at: Publisher Site
 D. DePasquale, A. C. Charania, and J. R. Olds, “The emerging orbital space tourism industry: new insight into demand and prospect for success,” in AIAA Paper 20067478, Atlanta, Georgia, 2006. View at: Google Scholar
 D. Dirkx and E. Mooij, “Optimization of entryvehicle shapes during conceptual design,” Acta Astronautica, vol. 94, no. 1, pp. 198–214, 2014. View at: Publisher Site  Google Scholar
 G. Ridolfi, E. Mooji, D. Dirkx, and S. Corpino, “Robust multidisciplinary optimization of unmanned entry capsules,” in AIAA Paper 20125006, Minneapolis, Minnesota, 2012. View at: Google Scholar
 L. F. Rowell, R. D. Braun, J. R. Olds, and R. Unal, “Multidisciplinary conceptual design optimization of space transportation systems,” Journal of Aircraft, vol. 36, no. 1, pp. 218–226, 1999. View at: Publisher Site  Google Scholar
 M. A. Lobbia, “Multidisciplinary design optimization of waveriderderived crew reentry vehicles,” Journal of Spacecraft and Rockets, vol. 54, no. 1, pp. 233–245, 2017. View at: Publisher Site  Google Scholar
 N. F. Foster and G. S. Dulikravich, “Threedimensional aerodynamic shape optimization using genetic and gradient search algorithms,” Journal of Spacecraft and Rockets, vol. 34, no. 1, pp. 36–42, 1997. View at: Publisher Site  Google Scholar
 T.t. Zhang, Z.g. Wang, W. Huang, and L. Yan, “Parameterization and optimization of hypersonicgliding vehicle configurations during conceptual design,” Aerospace Science and Technology, vol. 58, pp. 225–234, 2016. View at: Publisher Site  Google Scholar
 W. S. Hinman and C. T. Johansen, “Rapid prediction of hypersonic blunt body flows for parametric design studies,” Aerospace Science and Technology, vol. 58, pp. 48–59, 2016. View at: Publisher Site  Google Scholar
 S. Kumar and S. P. Mahulikar, “Design of thermal protection system for reusable hypersonic vehicle using inverse approach,” Journal of Spacecraft and Rockets, vol. 54, no. 2, pp. 436–446, 2017. View at: Publisher Site  Google Scholar
 D. Olynick, “Trajectorybased thermal protection system sizing for an X33 winged vehicle concept,” Journal of Spacecraft and Rockets, vol. 35, no. 3, pp. 249–257, 1998. View at: Publisher Site  Google Scholar
 A. Mazzaracchio, “Thermal protection system and trajectory optimization for orbital plane change aeroassisted maneuver,” Journal of Aerospace Technology and Management, vol. 5, no. 1, pp. 49–64, 2013. View at: Publisher Site  Google Scholar
 J. E. Bradford, “Thermal protection system sizing and selection for RLVs using the sentry code,” in AIAA Paper 20064605, 2006. View at: Google Scholar
 K. E. Wurster and H. W. Stone, “Aerodynamic heating environment definition/thermal protection system selection for the HL20,” Journal of Spacecraft and Rockets, vol. 30, no. 5, pp. 549–557, 1993. View at: Publisher Site  Google Scholar
 A. Aprovitola, L. Iuspa, C. Rainone, and A. Viviani, “Soft objects based procedure for thermal protection system modelling of reusable launch vehicles,” in HiSST: International Conference on HighSpeed Vehicle & Science Technology, Moscow, Russia, November 2018. View at: Google Scholar
 J. F. Blinn, “A generalization of algebraic surface drawing,” ACM Transaction on Graphics, vol. 68, no. 3, pp. 123–134, 1982. View at: Google Scholar
 O. Gourmel, L. Barthe, M. P. Cani et al., “A gradientbased implicit blend,” ACM Transactions on Graphics, vol. 32, no. 2, pp. 1–12, 2013. View at: Publisher Site  Google Scholar
 B. Angles, M. Tarini, B. Wyvill, L. Barthe, and A. Tagliasacchi, “Sketchbased implicit blending,” ACM Transactions on Graphics, vol. 36, no. 6, pp. 1–13, 2017. View at: Publisher Site  Google Scholar
 L. Iuspa, “Free topology generation of selfstiffened panels using skeletonbased integral soft objects,” Computers & Structures, vol. 158, pp. 184–210, 2015. View at: Publisher Site  Google Scholar
 C. Blanc and C. Schlick, “Extended field functions for soft objects,” in Proc. Implicit Surface, pp. 21–32, Grenoble, France, 1995. View at: Google Scholar
 A. Viviani, L. Iuspa, and A. Aprovitola, “An optimizationbased procedure for selfgeneration of reentry vehicles shape,” Aerospace Science and Technology, vol. 68, pp. 123–134, 2017. View at: Publisher Site  Google Scholar
 A. Viviani, L. Iuspa, and A. Aprovitola, “Multiobjective optimization for reentry spacecraft conceptual design using a freeform shape generator,” Aerospace Science and Technology, vol. 71, pp. 312–324, 2017. View at: Publisher Site  Google Scholar
 ANSYS, “Workbench Documentation.”. View at: Google Scholar
 R. L. Carmichael, “Public domain computer programs for the aeronautical engineer, version 1. CDROM,” Tech. Rep., Public Domain Aeronautical Software, 2015, http://www.pdas.com. View at: Google Scholar
 M. J. Grant and R. D. Braun, “Analytic hypersonic aerodynamics for conceptual design of entry vehicles,” in AIAA Paper 20101212, Orlando, Florida, 2010. View at: Google Scholar
 D. Filkovic, “APAME 3D panel method,” http://www.3dpanelmethod.com. View at: Google Scholar
 D. E. Myers, C. J. Martin, and M. L. Blosser, “Parametric weight comparison of advanced metallic, ceramic tile, and ceramic blanket thermal protection systems,” Tech. Rep., NASA/TM210289, 2000. View at: Google Scholar
 TPSX Material Properties Database, “National aeronautics and space administration,” http://www.tpsx.arc.nasa.gov. View at: Google Scholar
 R. R. Rohrschneider, Development of a mass estimating relationship database for launch vehicle conceptual design masters project, School of Aerospace Engineering, Georgia Institute of Technology, 2002.
 C. Weiland, Computational Space Flight Mechanics, Springer, 2010. View at: Publisher Site
 M. E. Tauber and H. G. Adelman, “Thermal environment of transatmospheric vehicles,” Journal of Aircraft, vol. 25, no. 4, pp. 355–363, 1988. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2019 Andrea Aprovitola et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.