# A Novel Method for Blunting the Leading Edge of Waverider with Specified Curvature

Authors: Xiang Fu; Jun Liu; Feng Ding; Shengye Wang; Tianheng Dou
Journal: International Journal of Aerospace Engineering
 (2020)

---

## Abstract

An ideal waverider has an infinite sharp leading edge, which causes difficulty for manufacture and aerothermal protection. Therefore, the leading edge of the waverider must be blunted. For this purpose, a parametric method for blunting the leading edge of the waverider is proposed here, which can fulfill the goals of setting a leading-edge blunt radius, achieving geometric continuity, and realizing the parametric design. First is the blunting procedure of the proposed method incorporating the construction of two-dimensional blunt curves and the integration of these curves on a three-dimensional waverider configuration. Second, waveriders blunted with different geometric continuities are built with corresponding computing grids generated. Numerical methods are then introduced and validated by the benchmark cases. Finally, results from these blunted configurations are presented and compared in terms of their geometric and flow characteristics. It shows that the proposed method has a better performance in the head region of the waverider and is thereby more suitable for the practical design.


---

## Body

## 1. Introduction

The waverider constitutes one of the most promising configurations for hypersonic flight owing to its superior aerodynamic performance and thus finds applications in various fields. Such a configuration in general is analytically designed via streamline tracing in a pregenerated or preselected supersonic/hypersonic flow field that contains the shock wave, thereby incorporating an infinitely sharp leading edge with the shock wave attached. Owing to this, the high-pressure stream behind the shock wave is restricted under the waverider and the high lift-to-drag ratio is consequently achieved, which is crucial in the hypersonic flight.Since the waverider was first invented in 1959 [1], efforts has been made to improve its performance and adaptability. As a result, various design methods and strategies have been proposed. These investigations are related to three major issues, that is, the design of basic flow fields, the parameterization and optimization, and the remodel design. The first concerns designing a single flow field or assembling multiple flow fields. Research objects of the former include flows over wedges [2, 3], cones [4–8], constricting ducts [9], and other three-dimensional base bodies [10, 11]. For the latter, parts from similar flow fields are assembled to achieve a larger design space. In the osculating methods [12–15], parts of flow fields are combined in the spanwise direction. And the multistage compression waverider [16] is built on multiple conical flow fields, which are in the longitudinal arrangement. As for the second issue, it is common to parameterize the waverider by the characteristic profile curves, which include the upper/lower surface profile curve (USPC/LSPC), the leading edge profile curve (LEPC) or its projection, the shock wave profile curve (SWPC), and the hybrid design of these curves [17, 18]. Based on the parameterization, we can proceed with the optimization of the waverider. The majority of research focuses on the improvement of the lift-to-drag ratio [19, 20], while the interests of the community also involve the mission-oriented optimization [21, 22] and the multiobject optimization [23]. The third theme focuses on modifying the geometry of the ideal waverider for realistic purposes. For instance, two waveriders are assembled to get a better wide-speed-range performance [24]; the outer wings are installed on the waverider to improve the low-speed performance [25]. Furthermore, a creative example is to remodel the design principle of characteristic profile curves to construct an airframe/inlet integrated waverider [26]. As the waverider has been extensively researched to this date, it is chosen for several practical hypersonic flight projects [27–29]. Nevertheless, some of the practical considerations proposed in Ref. [30] still needs close attention to further the application of the waverider in practice.When compared with blunt vehicles like the space shuttle, the sharp leading edge of the ideal waverider makes it unqualified for hypersonic flight. Therefore, the sharp leading edge must be blunted in practice for manufacturability improvement, aerothermal protection, and some handling concerns. Two major aspects have been actively pursued in various blunting methods. In one aspect, the sharp leading-edge geometry is modified by adding or removing the material [30, 31] with a special focus on the volume of the waverider. Specifically, while the former lifts the upper surface by separating it from the lower one to preserve the planform of the waverider and increase its volume, the latter truncates the sharp leading edge until a wanted blunt geometry is achieved. The other aspect concerns the choice of blunt curves that include simple ones like circular and elliptical arcs and parametric ones like Bézier curves [32] and rational Bézier curves [33]. The employment of either circular or elliptical arcs can explicitly specify the leading-edge curvature but causes geometric discontinuity at the connection of the blunted part and the original geometry. In contrast, when a single parametric curve is adopted, the higher-order geometric continuity is realized, yet without the blunt radius at the leading-edge point being directly specified. Moreover, in addition to these two aspects, the effects of different bluntness on the performance of the waverider have also been studied [34]. It shows that the heat-transfer rate at the stagnation point varies inversely with the square root of the head blunt radius, while the lift-to-drag ratio inversely with the head blunt radius [35, 36]. To accommodate the variation of thermal load at different positions and to reduce aerodynamic losses, a method with the nonuniform blunt radius is proposed and studied [37].To pursue this line of research, this paper is aimed at developing an efficient method based on parametric curves to generate the blunted waverider configuration. The proposed method mainly concerns the design of the blunt curve and involves two cubic Bézier curves to construct the blunted leading edge. Being tentatively viewed as an extension of the method using a single Bézier curve, this method can fulfill the goals of setting the leading-edge blunt radius, achieving second-order geometric continuity, and realizing parametric design, which cannot be fully accomplished by its counterparts. The design approach in detail is described in Section2, and physical models and numerical methods are presented in Section 3. The discussion over the influence of different blunting methods on geometric and flow characteristics of waverider is discussed in Section 4.


## 2. Method Description

A typical waverider configuration consists of an upper surface, a compression surface, and a base surface. The longitudinal cross-section of an ideal waverider configuration can be viewed as a forward wedge-like geometry, which allows for parametric curves to be used to generate a curved blunt leading-edge shape [33]. An outline of the proposed blunting method in this paper can be divided into two phases: first, use two cubic Bézier curves to fix the shape of the leading edge for the two-dimensional case; secondly, extending the aforementioned process alongside the spanwise direction to realize the integration on a three-dimensional waverider configuration. The adding-material strategy is adopted in this paper, while it does not influence the applicability of this method with the removing-material strategy.


### 2.1. Construction of 2D Leading Edge

Two cubic Bézier curves are used in the blunting procedure at the longitudinal section of the waverider. As shown in Figure1, the upper surface profile is vertically separated from the lower for a distance of h. The original leading point is split into two new endpoints, namely, P0 and P6. The original wedge-like geometry is blunted and replaced with two cubic Bézier curves P0P3 and P3P6. The new leading-edge point at this section is reset to P3, which has a relative horizontal distance lLE and a relative vertical distance hLE to P0, the original leading point.Figure 1
Illustration of key points on the blunt curves.Solving a cubic Bézier curve from the interpolation conditions is referred to as the geometric Hermite interpolation (GHI) [38] problem. The interpolation conditions here are the unit tangent vector and the signed curvature (see Figure 2 for all the involved conditions). The unit tangent vectors v0 and v2 depend on the geometry of the original waverider geometry, so it is for the signed curvatures k0 and k2. Nevertheless, the unit tangent vector v1 is fixed to vertical upward, while the signed curvature at P3 is valued artificially.Figure 2
Illustration of interpolation conditions on the blunt curves.To be specific, the curveP0P3 is chosen to illustrate the solution procedure of the cubic Bézier curve. P3P6 can be solved in the same manner. Note that the endpoints P0 and P3 are already specified, while the control points P1 and P2 are to be solved.The positions ofP1 and P2 are defined by
(1)P1=P0+δ0v0,P2=P3−δ1v1,and a system of quadratic equations is constructed for solving δ [38], which is written as
(2)v0×v1δ0=d×v1−32k1δ12,v0×v1δ1=v0×d−32k0δ02,in which d=P3−P0.When blunting the leading edge of waverider, the results ofv0×v1, d×v1, and v0×d are apparently equal to zero in no case, while k0 and k2 are likely to be equal to or proximate to zero. In this case, Equation (2) can be simplified by introducing the new parameters ρ and S, which are defined by
(3)δ0=ρ0d×v1v0×v1,δ1=ρ1v0×dv0×v1,(4)S0=32k0d×v12v0×dv0×v12,S1=32k1v0×d2d×v1v0×v12.Substituting Equations (3) and (4) into Equation (2) yields an equivalent system
(5)ρ0=1−S1ρ12,ρ1=1−S0ρ02.By solving Equation (5), the parameters δ are obtained, and the cubic Bézier curve P0P3 is constructed. The value of S0 determines the solution of Equation (5). When S0 is nonzero, Equation (5) is transformed into a quartic polynomial equation of ρ0 as Equation (6), the roots of which are equal to the eigenvalues of its companion matrix. As there would be multiple real roots, the affine invariants λ and μ [39] are used to check the shape of the cubic Bézier curve.
(6)ρ04−2S0ρ02+2S02S1ρ0+S1−1S02S1=0.


### 2.2. Integration on 3D Waverider Geometry

When integrating the two-dimensional blunt curve mentioned earlier in the three-dimensional configuration, factors to be considered include the adding/removing material method, the direction of blunting curves, the variation of leading-edge radius, etc. The preferred methods are presented in this subsection.As the adding material method is adopted here, the upper surface of the waverider is separated vertically from the lower one for a distance ofh. Nevertheless, the transition distance varies alongside the leading edge when the variable blunt radius is considered. Here, h is set to be the function of the spanwise coordinate, namely, hz. If the method of removing material is chosen, the truncation length varies alongside the leading edge [33], too.The blunting curve is sketched on the reference plane that is parallel to the symmetry plane and the free stream. The merit lies in the blunt radius front straight to the free stream being directly specified. Besides, the reference plane can also be perpendicular to the leading edge or parallel to the osculating plane [33].It is important and meaningful to control how the radius of curvature is distributed along the spanwise direction, from the head to the tip. In Ref. [32], Equation (7) is proposed for the off-centerline leading edge radius, RSW, which is written as a function of the centerline radius, RCT, and the local leading-edge sweep, λSW. With this equation, the heating rate at the off-centerline leading-edge position would be the same as the centerline, which makes the most use of the thermal protection system.
(7)RSW=RCTcosλSW2.2.The blunt leading edge of the three-dimensional waverider is constructed by lofting all the two-dimensional leading edges, which is implemented through computer-aided design (CAD).


## 2.1. Construction of 2D Leading Edge

Two cubic Bézier curves are used in the blunting procedure at the longitudinal section of the waverider. As shown in Figure1, the upper surface profile is vertically separated from the lower for a distance of h. The original leading point is split into two new endpoints, namely, P0 and P6. The original wedge-like geometry is blunted and replaced with two cubic Bézier curves P0P3 and P3P6. The new leading-edge point at this section is reset to P3, which has a relative horizontal distance lLE and a relative vertical distance hLE to P0, the original leading point.Figure 1
Illustration of key points on the blunt curves.Solving a cubic Bézier curve from the interpolation conditions is referred to as the geometric Hermite interpolation (GHI) [38] problem. The interpolation conditions here are the unit tangent vector and the signed curvature (see Figure 2 for all the involved conditions). The unit tangent vectors v0 and v2 depend on the geometry of the original waverider geometry, so it is for the signed curvatures k0 and k2. Nevertheless, the unit tangent vector v1 is fixed to vertical upward, while the signed curvature at P3 is valued artificially.Figure 2
Illustration of interpolation conditions on the blunt curves.To be specific, the curveP0P3 is chosen to illustrate the solution procedure of the cubic Bézier curve. P3P6 can be solved in the same manner. Note that the endpoints P0 and P3 are already specified, while the control points P1 and P2 are to be solved.The positions ofP1 and P2 are defined by
(1)P1=P0+δ0v0,P2=P3−δ1v1,and a system of quadratic equations is constructed for solving δ [38], which is written as
(2)v0×v1δ0=d×v1−32k1δ12,v0×v1δ1=v0×d−32k0δ02,in which d=P3−P0.When blunting the leading edge of waverider, the results ofv0×v1, d×v1, and v0×d are apparently equal to zero in no case, while k0 and k2 are likely to be equal to or proximate to zero. In this case, Equation (2) can be simplified by introducing the new parameters ρ and S, which are defined by
(3)δ0=ρ0d×v1v0×v1,δ1=ρ1v0×dv0×v1,(4)S0=32k0d×v12v0×dv0×v12,S1=32k1v0×d2d×v1v0×v12.Substituting Equations (3) and (4) into Equation (2) yields an equivalent system
(5)ρ0=1−S1ρ12,ρ1=1−S0ρ02.By solving Equation (5), the parameters δ are obtained, and the cubic Bézier curve P0P3 is constructed. The value of S0 determines the solution of Equation (5). When S0 is nonzero, Equation (5) is transformed into a quartic polynomial equation of ρ0 as Equation (6), the roots of which are equal to the eigenvalues of its companion matrix. As there would be multiple real roots, the affine invariants λ and μ [39] are used to check the shape of the cubic Bézier curve.
(6)ρ04−2S0ρ02+2S02S1ρ0+S1−1S02S1=0.


## 2.2. Integration on 3D Waverider Geometry

When integrating the two-dimensional blunt curve mentioned earlier in the three-dimensional configuration, factors to be considered include the adding/removing material method, the direction of blunting curves, the variation of leading-edge radius, etc. The preferred methods are presented in this subsection.As the adding material method is adopted here, the upper surface of the waverider is separated vertically from the lower one for a distance ofh. Nevertheless, the transition distance varies alongside the leading edge when the variable blunt radius is considered. Here, h is set to be the function of the spanwise coordinate, namely, hz. If the method of removing material is chosen, the truncation length varies alongside the leading edge [33], too.The blunting curve is sketched on the reference plane that is parallel to the symmetry plane and the free stream. The merit lies in the blunt radius front straight to the free stream being directly specified. Besides, the reference plane can also be perpendicular to the leading edge or parallel to the osculating plane [33].It is important and meaningful to control how the radius of curvature is distributed along the spanwise direction, from the head to the tip. In Ref. [32], Equation (7) is proposed for the off-centerline leading edge radius, RSW, which is written as a function of the centerline radius, RCT, and the local leading-edge sweep, λSW. With this equation, the heating rate at the off-centerline leading-edge position would be the same as the centerline, which makes the most use of the thermal protection system.
(7)RSW=RCTcosλSW2.2.The blunt leading edge of the three-dimensional waverider is constructed by lofting all the two-dimensional leading edges, which is implemented through computer-aided design (CAD).


## 3. Physical Models and Numerical Methods

The generation of two sharp waveriders based on the osculating cone-derived (OCD) method [12] is presented in this section. The parameters of the proposed blunting method are described, accompanied by the introduction of the comparable blunting method that uses the elliptical blunt curve. In the second subsection, the numerical methods are illustrated and validated via two benchmark problems.


### 3.1. Physical Models

When generating an ideal waverider, the compression surface should be given the priority, as it decisively influences the aerodynamic performance of the waverider. Here, the classical OCD method is adopted to generate the lower surface, while the upper surface profile curve (USPC) and the shock wave profile curve (SWPC) are designed as the control curves (see Figure3). The design parameters include the length L, the width W, the design Mach number Ma, and the shock angle β, which are listed in Table 1. It should be noted that all the waveriders presented in this paper have the same characteristic sizes, which are 2100 mm in length and 960 mm in width.Figure 3
Control curves of an OCD waverider on the base plane.Table 1
Design parameters of the compression surface of sharp waveriders.
ParametersL (mm)W (mm)Maβ (°)Values2100.0960.08.09.0To better investigate the influence of different blunting methods, two waveriders that have the same compression surface yet different upper surfaces are generated. The first waverider whose configuration is referred to as WR1 has a free-stream upper surface. The second waverider named WR2 has a reshaped upper surface to realize a better volumetric characteristic and a higher lift-to-drag ratio. To be specific, while its volume is concentrated to the middle of the waverider, its sides are thinner than the ideal one, just like the HiFIRE-4 vehicle [28]. This kind of design is more demanding on the blunting method than the free-stream surface, as it has a variable geometry characteristic on the leading edge. The two sharp waveriders are shown in Figure 4.Figure 4
Isometric views of two sharp waveriders with different upper surfaces: (a) WR1; (b) WR2.
(a)(b)For the sake of comparison, a blunting method that provides zero-order geometric continuity is introduced. It uses a half elliptical arc to connect the upper and lower surface profiles of the waverider. To be brief, the new blunting method proposed in this paper is referred to as the G2 method, while the method concerning the elliptical curve is referred to as the G0 method. As for the G0 method, its leading-edge point is exactly the vertex of the elliptical arc, and its position is related to the leading-edge radius and the transition distance of the upper surface. The new leading point of the G2 method is placed on the bisector of the tangent lines passing through the two endpoints, as shown in Figure 5.Figure 5
Difference in the position of the new leading-edge point: (a) G0 method (black); (b) G2 method (red).The blunt radius and the vertical transition distance are set to the same value when blunting waveriders through the two methods. In Figure6, the blunt radius calculated from Equation (7) and the design variation are demonstrated together. The former is taken as a reference line for the design of the blunt radius. It shows that the design variation ranges from 10 mm at the center plane to 0.63 mm at the wingtip.Figure 6
Variations of blunt radius alongside the spanwise direction.The WR1 and WR2 configurations are blunted by the G0 and G2 methods, respectively. Consequently, four blunted waveriders named WR1-G0, WR1-G2, WR2-G0, and WR2-G2 are obtained. It is noted that all four blunted configurations have the same planform and the same distribution of blunt radius. Figure 7 gives the isometric view of the WR2-G2 configuration.Figure 7
Isometric view of the WR2-G2 configuration.


### 3.2. Numerical Methods

The waveriders in this paper are all designed based on the inviscid flow field, and no viscous correction is implemented on these configurations. Therefore, both the inviscid and the viscous calculation are considered when evaluating the performance of waverider by the computational fluid dynamic (CFD) method.To simulate the viscous flow around the waverider, the three-dimensional Reynolds averaged Navier-Stokes (RANS) equations are solved by the implicit finite volume method (FVM). The governing equations of the mean flow are closed by the shear-stress transport (SST)k-ω turbulence model, a two-equation model that is popular in both industry and academics. The second-order spatially accurate upwind scheme with the advection upstream (AUSM) splitting approach is employed to split the flux vector. The Courant-Friedrichs-Lewy (CFL) number is maintained at 0.5 to ensure stability. The convergence is reached if the residuals fall by more than four orders of magnitude, and the difference between the computed inflow and outflow mass flux falls below 0.1%. The boundary condition of the wall is assumed to be isothermal (Tw=300K) and no-slip. The air is assumed to be a thermally and calorically perfect gas with the viscosity modelled according to the well-known Sutherland law. For the inviscid calculation, the Euler equations are solved in the similar manner as the RANS equations.Two hypersonic experimental cases [40] are chosen to validate the accuracy of the numerical methods applied in this paper. The first validation case VC1 focuses on measuring the surface pressure of the double ellipsoid, while the second one VC2 concentrates on the measurement of the surface heat flux. The experimental conditions are given in Table 2. y+ of the first cell height and the cell Reynolds number Regrid are used as two criteria to guide the generation of the waverider grid: y+ should be set to smaller than 1 for the simulation of force [41]; Regrid should be set to smaller than 10 for the simulation of heat flux [42]. Thus, a set of multiblock structured grids is generated for the numerical simulation. As depicted in Figure 8, it has 2,078,848 cells and a first cell height of 0.005 mm. For the VC1 case, the maximum value of y+ is 1.10, while Regrid for the VC2 case is 1.45.
(8)y+=ρwuτΔyμw,Regrid=ρ∞v∞Δyμ∞.Table 2
Experimental conditions of the two validation cases.
CasesMaRe (106 × m-1)α (°)P0 (MPa)Τ0 (K)VC18.0219.808.5720VC210.022.206.91457Figure 8
Illustration of the structured grid for the two validation cases.In this paper, the pressure coefficientCP and the Stanton number St are employed as normalizations of the surface pressure and the surface heat flux, respectively, for further comparison. The definitions of CP and St are expressed in
(9)CP=p−p∞1/2ρ∞v∞2,St=qTw−T∞ρ∞v∞cp∞.The results of the numerical simulation are compared with the experimental data (see Figure9 for details). While there exists a discrepancy, the changing trends of numerical results are consistent with the experimental data, which validates the numerical methods applied in this paper.Figure 9
Comparisons between numerical results and experimental data: (a)CP results from VC1; (b) St results from VC2.
(a)(b)When it comes to the blunted waveriders, the simulated flight is executed atMa=8.0 and H=30km with an angle of attack of 0°, which is shown in Table 3. The computational grids for the four blunted waveriders have 1,529,616 cells, with its first cell height of grids set to 0.001 mm, which corresponds to y+=0.55 and Regrid=3.68. Figure 10 gives the overall view and the close-up view of the WR2-G2 mesh.Table 3
Computational conditions for the blunt waveriders.
CasesMaRe (106 × m-1)α (°)P0 (MPa)Τ0 (K)Blunt waveriders8.02.95011.443127.77Figure 10
Illustration of the structured grid for WR2-G2: (a) overall view; (b) close-up view on the head.
(a)(b)


## 3.1. Physical Models

When generating an ideal waverider, the compression surface should be given the priority, as it decisively influences the aerodynamic performance of the waverider. Here, the classical OCD method is adopted to generate the lower surface, while the upper surface profile curve (USPC) and the shock wave profile curve (SWPC) are designed as the control curves (see Figure3). The design parameters include the length L, the width W, the design Mach number Ma, and the shock angle β, which are listed in Table 1. It should be noted that all the waveriders presented in this paper have the same characteristic sizes, which are 2100 mm in length and 960 mm in width.Figure 3
Control curves of an OCD waverider on the base plane.Table 1
Design parameters of the compression surface of sharp waveriders.
ParametersL (mm)W (mm)Maβ (°)Values2100.0960.08.09.0To better investigate the influence of different blunting methods, two waveriders that have the same compression surface yet different upper surfaces are generated. The first waverider whose configuration is referred to as WR1 has a free-stream upper surface. The second waverider named WR2 has a reshaped upper surface to realize a better volumetric characteristic and a higher lift-to-drag ratio. To be specific, while its volume is concentrated to the middle of the waverider, its sides are thinner than the ideal one, just like the HiFIRE-4 vehicle [28]. This kind of design is more demanding on the blunting method than the free-stream surface, as it has a variable geometry characteristic on the leading edge. The two sharp waveriders are shown in Figure 4.Figure 4
Isometric views of two sharp waveriders with different upper surfaces: (a) WR1; (b) WR2.
(a)(b)For the sake of comparison, a blunting method that provides zero-order geometric continuity is introduced. It uses a half elliptical arc to connect the upper and lower surface profiles of the waverider. To be brief, the new blunting method proposed in this paper is referred to as the G2 method, while the method concerning the elliptical curve is referred to as the G0 method. As for the G0 method, its leading-edge point is exactly the vertex of the elliptical arc, and its position is related to the leading-edge radius and the transition distance of the upper surface. The new leading point of the G2 method is placed on the bisector of the tangent lines passing through the two endpoints, as shown in Figure 5.Figure 5
Difference in the position of the new leading-edge point: (a) G0 method (black); (b) G2 method (red).The blunt radius and the vertical transition distance are set to the same value when blunting waveriders through the two methods. In Figure6, the blunt radius calculated from Equation (7) and the design variation are demonstrated together. The former is taken as a reference line for the design of the blunt radius. It shows that the design variation ranges from 10 mm at the center plane to 0.63 mm at the wingtip.Figure 6
Variations of blunt radius alongside the spanwise direction.The WR1 and WR2 configurations are blunted by the G0 and G2 methods, respectively. Consequently, four blunted waveriders named WR1-G0, WR1-G2, WR2-G0, and WR2-G2 are obtained. It is noted that all four blunted configurations have the same planform and the same distribution of blunt radius. Figure 7 gives the isometric view of the WR2-G2 configuration.Figure 7
Isometric view of the WR2-G2 configuration.


## 3.2. Numerical Methods

The waveriders in this paper are all designed based on the inviscid flow field, and no viscous correction is implemented on these configurations. Therefore, both the inviscid and the viscous calculation are considered when evaluating the performance of waverider by the computational fluid dynamic (CFD) method.To simulate the viscous flow around the waverider, the three-dimensional Reynolds averaged Navier-Stokes (RANS) equations are solved by the implicit finite volume method (FVM). The governing equations of the mean flow are closed by the shear-stress transport (SST)k-ω turbulence model, a two-equation model that is popular in both industry and academics. The second-order spatially accurate upwind scheme with the advection upstream (AUSM) splitting approach is employed to split the flux vector. The Courant-Friedrichs-Lewy (CFL) number is maintained at 0.5 to ensure stability. The convergence is reached if the residuals fall by more than four orders of magnitude, and the difference between the computed inflow and outflow mass flux falls below 0.1%. The boundary condition of the wall is assumed to be isothermal (Tw=300K) and no-slip. The air is assumed to be a thermally and calorically perfect gas with the viscosity modelled according to the well-known Sutherland law. For the inviscid calculation, the Euler equations are solved in the similar manner as the RANS equations.Two hypersonic experimental cases [40] are chosen to validate the accuracy of the numerical methods applied in this paper. The first validation case VC1 focuses on measuring the surface pressure of the double ellipsoid, while the second one VC2 concentrates on the measurement of the surface heat flux. The experimental conditions are given in Table 2. y+ of the first cell height and the cell Reynolds number Regrid are used as two criteria to guide the generation of the waverider grid: y+ should be set to smaller than 1 for the simulation of force [41]; Regrid should be set to smaller than 10 for the simulation of heat flux [42]. Thus, a set of multiblock structured grids is generated for the numerical simulation. As depicted in Figure 8, it has 2,078,848 cells and a first cell height of 0.005 mm. For the VC1 case, the maximum value of y+ is 1.10, while Regrid for the VC2 case is 1.45.
(8)y+=ρwuτΔyμw,Regrid=ρ∞v∞Δyμ∞.Table 2
Experimental conditions of the two validation cases.
CasesMaRe (106 × m-1)α (°)P0 (MPa)Τ0 (K)VC18.0219.808.5720VC210.022.206.91457Figure 8
Illustration of the structured grid for the two validation cases.In this paper, the pressure coefficientCP and the Stanton number St are employed as normalizations of the surface pressure and the surface heat flux, respectively, for further comparison. The definitions of CP and St are expressed in
(9)CP=p−p∞1/2ρ∞v∞2,St=qTw−T∞ρ∞v∞cp∞.The results of the numerical simulation are compared with the experimental data (see Figure9 for details). While there exists a discrepancy, the changing trends of numerical results are consistent with the experimental data, which validates the numerical methods applied in this paper.Figure 9
Comparisons between numerical results and experimental data: (a)CP results from VC1; (b) St results from VC2.
(a)(b)When it comes to the blunted waveriders, the simulated flight is executed atMa=8.0 and H=30km with an angle of attack of 0°, which is shown in Table 3. The computational grids for the four blunted waveriders have 1,529,616 cells, with its first cell height of grids set to 0.001 mm, which corresponds to y+=0.55 and Regrid=3.68. Figure 10 gives the overall view and the close-up view of the WR2-G2 mesh.Table 3
Computational conditions for the blunt waveriders.
CasesMaRe (106 × m-1)α (°)P0 (MPa)Τ0 (K)Blunt waveriders8.02.95011.443127.77Figure 10
Illustration of the structured grid for WR2-G2: (a) overall view; (b) close-up view on the head.
(a)(b)


## 4. Results and Discussion

In this section, the influence of different blunting methods on the characteristics of waveriders is studied and discussed, including geometric characteristics and flow characteristics.


### 4.1. Comparison of Geometric Characteristics

Table4 gives an overall picture of geometric characteristics of the two original sharp waveriders and the four blunt waveriders, including the total surface area S, the area of the base surface SB, the total volume V, and the available volume VA. The available volume corresponds to the inner region that has a distance larger than 25 mm from the nearest upper and lower surfaces. It shows that the change of the upper surface has a bigger impact on the area of the base surface and the total volume than the total surface area, while the impact on the available volume can be reduced to a negligible level by the design of the upper surface. Besides, owing to the control of blunting parameters, there is nearly no difference on the geometric characteristics between the waveriders blunted by two blunting methods.Table 4
Overall geometric characteristics of the blunted waveriders.
ConfigurationsS (m2)SB (m2)V (m3)VA (m3)WR12.650.07090.05640.0161WR22.690.02450.04190.0157WR1-G02.780.08080.07040.0229WR1-G22.770.08090.07050.0229WR2-G02.840.03440.05590.0213WR2-G22.830.03440.05590.0213The blunt curve on the central plane is observed for further comparison. For the WR1 configurations, the blunting curves from the G0 and G2 methods are similar; a relatively large deviation is also found on the upper curves, which is depicted in Figure 11(a). When it comes to the WR2 configurations, there is an apparent discrepancy between the two blunting curves. In Figure 11(b), the curve generated by the G2 method is lower than its counterpart and shows a smoother connection with the base geometry.Figure 11
Blunting curves on the central plane of waveriders: (a) blunting curves from WR1-G0 and WR1-G2; (b) blunting curves from WR2-G0 and WR2-G2.
(a)(b)


### 4.2. Comparison of Flow Characteristics

Table5 gives the overall aerodynamic characteristics of all six waveriders, including the drag coefficient CD, the lift coefficient CL, and the lift-to-drag ratio L/D. Note that both the viscous and inviscid results are listed in Table 5. It is found that the viscosity has large influence on the aerodynamic characteristics of waverider, as the drag coefficient increases greatly after adding the viscosity. Meanwhile, the impact of bluntness is also considerable. However, the aerodynamic difference of two blunting methods is negligible, as the four blunted waveriders have the same distribution of bluntness alongside the leading edge.Table 5
Overall aerodynamic performances of the blunted waveriders.
CasesInviscid calculationViscous calculationCD×10−2CL×10−2L/DCD×10−2CL×10−2L/DWR10.06901.1716.960.2701.204.44WR20.1131.5413.630.2961.535.16WR1-G00.1711.136.610.3891.142.92WR1-G20.1721.156.680.3911.152.95WR2-G00.2081.456.970.4051.433.53WR2-G20.2111.476.970.4061.453.58Figure12 compares the surface characteristics on the central blunting curve between two blunting methods. For the G0 method, the distributions of CP and St are nearly identical for the upper part and the lower part of the central blunting curves, which corresponds to the line of x≤0mm in Figure 12. This is because of the symmetry of the elliptical arc. However, an adverse pressure gradient appears around x=0mm for the WR2-G0 configuration, which corresponds to the geometric discontinuity in Figure 11(b). Owing to the geometric simplicity of the two WR1 configurations, there is almost no fluctuation of surface characteristics for both the WR1-G0 and WR1-G2 configurations. As there is asymmetry between the upper and the lower part of the central blunting curves, the distributions of CP and St are separated for the WR1-G2 and WR2-G2 configurations. And all the distributions on the WR1-G2 and WR2-G2 configurations show the smooth transition from the blunted part to the base geometry.Figure 12
Surface characteristics on the central blunt curve: (a) comparison between WR1-G0 and WR1-G2; (b) comparison between WR2-G0 and WR2-G2.
(a)(b)Figure13 illustrates the contour of CP on the central plane, and no significant difference is found between the two methods. In Figure 14, the contours of St around the head of the blunted waveriders are presented. For the configurations blunted by the G0 method, there exist second peaks of St besides the stagnation point, which inflects the geometric protrusion corresponding to these peaks. On the contrary, the distributions of heat flux are more concentrated for the WR1-G2 and WR2-G2 configurations.Figure 13
Contour ofCP on the central plane of the blunted waveriders: (a) WR1-G0; (b) WR1-G2; (c) WR2-G0; (d) WR2-G2.
(a)(b)(c)(d)Figure 14
Contour ofSt around the head of blunted waveriders: (a) WR1-G0; (b) WR1-G2; (c) WR2-G0; (d) WR2-G2.
(a)(b)(c)(d)Figures15 and 16 illustrate the contours of surface characteristics alongside the leading edge. It can be observed that both CP and St decrease alongside the leading edge to the downstream. However, there is no significant difference between geometries blunted by the two methods. When compared to the contours of CP, St shows a smoother decrease because of the design of the blunt radius.Figure 15
Contours of surface characteristics alongside leading edges of WR1-G0 and WR1-G2: (a) contours of CP; (b) contours of St.
(a)(b)Figure 16
Contours of surface characteristics alongside leading edges of WR2-G0 and WR2-G2: (a) contours of CP; (b) contours of St.
(a)(b)Figure17 presents both the inviscid and viscous contours of CP at the base plane of waveriders. As shown in Figures 17(a) and 17(b), the high-pressure flow is restricted under the lower surface and no flow leaks in the inviscid calculation, which indicates that the design procedure of the lower surface is valid under the design condition. However, the high-pressure air leaks to the upper surface when the viscosity is considered. When it comes the leading-edge bluntness, it has a more significant influence to the flow field of waverider than the viscosity, as the leak is much larger than in Figures 17(a) and 17(b). The influence of the blunting method is relatively small when compared with the bluntness and the viscosity. The waveriders blunted by the two methods have similar shape of shock wave and leak of flow.Figure 17
Contours ofCP on the base plane of waveriders: (a) WR1; (b) WR2; (c) WR1-G0; (d) WR1-G0; (e) WR1-G2; (f) WR2-G2.
(a)(b)(c)(d)(e)(f)


## 4.1. Comparison of Geometric Characteristics

Table4 gives an overall picture of geometric characteristics of the two original sharp waveriders and the four blunt waveriders, including the total surface area S, the area of the base surface SB, the total volume V, and the available volume VA. The available volume corresponds to the inner region that has a distance larger than 25 mm from the nearest upper and lower surfaces. It shows that the change of the upper surface has a bigger impact on the area of the base surface and the total volume than the total surface area, while the impact on the available volume can be reduced to a negligible level by the design of the upper surface. Besides, owing to the control of blunting parameters, there is nearly no difference on the geometric characteristics between the waveriders blunted by two blunting methods.Table 4
Overall geometric characteristics of the blunted waveriders.
ConfigurationsS (m2)SB (m2)V (m3)VA (m3)WR12.650.07090.05640.0161WR22.690.02450.04190.0157WR1-G02.780.08080.07040.0229WR1-G22.770.08090.07050.0229WR2-G02.840.03440.05590.0213WR2-G22.830.03440.05590.0213The blunt curve on the central plane is observed for further comparison. For the WR1 configurations, the blunting curves from the G0 and G2 methods are similar; a relatively large deviation is also found on the upper curves, which is depicted in Figure 11(a). When it comes to the WR2 configurations, there is an apparent discrepancy between the two blunting curves. In Figure 11(b), the curve generated by the G2 method is lower than its counterpart and shows a smoother connection with the base geometry.Figure 11
Blunting curves on the central plane of waveriders: (a) blunting curves from WR1-G0 and WR1-G2; (b) blunting curves from WR2-G0 and WR2-G2.
(a)(b)


## 4.2. Comparison of Flow Characteristics

Table5 gives the overall aerodynamic characteristics of all six waveriders, including the drag coefficient CD, the lift coefficient CL, and the lift-to-drag ratio L/D. Note that both the viscous and inviscid results are listed in Table 5. It is found that the viscosity has large influence on the aerodynamic characteristics of waverider, as the drag coefficient increases greatly after adding the viscosity. Meanwhile, the impact of bluntness is also considerable. However, the aerodynamic difference of two blunting methods is negligible, as the four blunted waveriders have the same distribution of bluntness alongside the leading edge.Table 5
Overall aerodynamic performances of the blunted waveriders.
CasesInviscid calculationViscous calculationCD×10−2CL×10−2L/DCD×10−2CL×10−2L/DWR10.06901.1716.960.2701.204.44WR20.1131.5413.630.2961.535.16WR1-G00.1711.136.610.3891.142.92WR1-G20.1721.156.680.3911.152.95WR2-G00.2081.456.970.4051.433.53WR2-G20.2111.476.970.4061.453.58Figure12 compares the surface characteristics on the central blunting curve between two blunting methods. For the G0 method, the distributions of CP and St are nearly identical for the upper part and the lower part of the central blunting curves, which corresponds to the line of x≤0mm in Figure 12. This is because of the symmetry of the elliptical arc. However, an adverse pressure gradient appears around x=0mm for the WR2-G0 configuration, which corresponds to the geometric discontinuity in Figure 11(b). Owing to the geometric simplicity of the two WR1 configurations, there is almost no fluctuation of surface characteristics for both the WR1-G0 and WR1-G2 configurations. As there is asymmetry between the upper and the lower part of the central blunting curves, the distributions of CP and St are separated for the WR1-G2 and WR2-G2 configurations. And all the distributions on the WR1-G2 and WR2-G2 configurations show the smooth transition from the blunted part to the base geometry.Figure 12
Surface characteristics on the central blunt curve: (a) comparison between WR1-G0 and WR1-G2; (b) comparison between WR2-G0 and WR2-G2.
(a)(b)Figure13 illustrates the contour of CP on the central plane, and no significant difference is found between the two methods. In Figure 14, the contours of St around the head of the blunted waveriders are presented. For the configurations blunted by the G0 method, there exist second peaks of St besides the stagnation point, which inflects the geometric protrusion corresponding to these peaks. On the contrary, the distributions of heat flux are more concentrated for the WR1-G2 and WR2-G2 configurations.Figure 13
Contour ofCP on the central plane of the blunted waveriders: (a) WR1-G0; (b) WR1-G2; (c) WR2-G0; (d) WR2-G2.
(a)(b)(c)(d)Figure 14
Contour ofSt around the head of blunted waveriders: (a) WR1-G0; (b) WR1-G2; (c) WR2-G0; (d) WR2-G2.
(a)(b)(c)(d)Figures15 and 16 illustrate the contours of surface characteristics alongside the leading edge. It can be observed that both CP and St decrease alongside the leading edge to the downstream. However, there is no significant difference between geometries blunted by the two methods. When compared to the contours of CP, St shows a smoother decrease because of the design of the blunt radius.Figure 15
Contours of surface characteristics alongside leading edges of WR1-G0 and WR1-G2: (a) contours of CP; (b) contours of St.
(a)(b)Figure 16
Contours of surface characteristics alongside leading edges of WR2-G0 and WR2-G2: (a) contours of CP; (b) contours of St.
(a)(b)Figure17 presents both the inviscid and viscous contours of CP at the base plane of waveriders. As shown in Figures 17(a) and 17(b), the high-pressure flow is restricted under the lower surface and no flow leaks in the inviscid calculation, which indicates that the design procedure of the lower surface is valid under the design condition. However, the high-pressure air leaks to the upper surface when the viscosity is considered. When it comes the leading-edge bluntness, it has a more significant influence to the flow field of waverider than the viscosity, as the leak is much larger than in Figures 17(a) and 17(b). The influence of the blunting method is relatively small when compared with the bluntness and the viscosity. The waveriders blunted by the two methods have similar shape of shock wave and leak of flow.Figure 17
Contours ofCP on the base plane of waveriders: (a) WR1; (b) WR2; (c) WR1-G0; (d) WR1-G0; (e) WR1-G2; (f) WR2-G2.
(a)(b)(c)(d)(e)(f)


## 5. Conclusions

A parametric method for blunting the three-dimensional hypersonic waverider is proposed in this paper. At the theoretical level, the proposed method satisfies the needs of setting a leading-edge blunt radius, achieving second-order geometric continuity, and realizing the parametric design. From the numerical investigation, the new method shows no inferior performance compared with the comparable method, while having a better geometric continuity and a higher degree of freedom. The practical design of the waverider is demanding on this method, as the upper surface has variable geometric characteristic alongside the leading edge. Nevertheless, it is also found that the geometric continuity has little influence on the overall aerodynamic performances when the blunting parameters are fixed. And the impact of geometric continuity on the flow transition needs further study.



---

*Source: 8885062-2020-12-21.xml*