Spillage-Adaptive Fixed-Geometry Bump Inlet of Wide Speed Range

: In this work, a new spillage-adaptive bump inlet concept is proposed to widen the speed range for hypersonic air-breathing ﬂight vehicles. Various approaches to improve the inlet start-ability are summarized and compared, among which the bump-inlet pattern holds the merits of high lift-to-drag ratio, boundary layer diversion, and ﬂexible integration ability. The proposed spillage-adaptive concept ensures the inlet starting performance by spilling extra mass ﬂow away at low speed number conditions. The inlet presetting position is determined by synthetically evaluating the ﬂow uniformity and the low-kinetic-energy ﬂuid proportion. The numerical results show that the ﬂow spillage of the inlet increases with the inﬂow speed decrease, which makes the inlet easier to start at low speed conditions (M 2.5–6.0). The effects of the boundary layer on spillage are also studied in this work. The new integration pattern releases the ﬂow spillage potentials of three-dimensional inward-turning inlets by reasonably arranging the inlet compression on the bump surface. Future work will focus on the spillage-controllable design method.


Introduction
Hypersonic inlet design is crucial to achieving the target of rapid global arrival [1][2][3][4]. The contribution of inlets to the propulsion system is estimated [5] as follows: According to the maximum temperature limit which the material can bear at present (the temperature limit is approximately 2000 K), the converted energy by inlet compression occupies 12% of the added energy by combustion at M 1.8, while it increases to 66.7% at M 3.4. This energy reaches approximately 230% at M 4.5. The trend reveals that a tiny increase of the energy loss in the inlets will induce considerable deterioration in the whole characteristics of the propulsion system at high-speed conditions [6][7][8].
The inlet start is the essential requirement to ensure its efficient working at off-design conditions. Abundant studies have been conducted to clarify the unstart mechanism and widen the inlet speed range. The direct cause of inlet unstart is that the swallowed inflow exceeds the flow capacity of the throat. Except for the backpressure fluctuation from the combustor [9][10][11], the internal contraction ratio (ICR) is the key parameter to determine the inlet start. Kantrowitz [12] first proposed the inviscid ICR calculating equation based on 1D steady flow assumption. The corresponding ICR limit is obtained when the ending shock is located at the entrance, and M throat approaches 1.0. Xie [13,14] developed a threelayer model to predict the ICR limit at viscous conditions. The improvement of inlet start performance is carried out in two directions: one is to decrease the internal compression ratio and the other is to remove the boundary layer. Wang [15] developed an image detection method to quantitively analyze the correlation between shock oscillation and the dynamic wall pressure in hypersonic experiments. Yu [16] proposed a high external compression ratio basic flow field. The speed range is widened by approximately 3.5% in M 3.5 through decreasing the external/internal contraction ratio by 38%. Xiong [17] developed and the incident shock strength. The former is related to the flow capacity decrease and the latter is related to low spillage ability.

SES Effects Caused by Overlarge Lateral Compression
The SES is observed in the authors' previous work [5] on a typical central-spillage bump inlet, which is demonstrated in Figure 1. The inflow is expanded over the bump and compressed in the inlet. The z-direction flow deflection of this inflow is larger than that of the parallel inflow. Accordingly, the lateral compression is larger. The compression is equivalent to the oblique shock wave of M∞ = 5, β = 12.31°, which forms a SES near the inlet bottom surface. The shock profile is clearly observed by the boundary between yellow and red colors. The SES effects will deteriorate the inlet performance in two aspects: (1) The internal shock system of the inlet distorts owing to the SES, which brings the additional energy loss; (2) the effective circulation area decreases, which forms an aerodynamic throat. Therefore, the flow capacity decreases and the inlet unstart can be easily induced.

Modular Assessment Approach for the Spillage Pattern Analysis
The direct cause for flow spillage is the detachment of the inlet incident shock. This section employs a modular assessment approach to estimate the spillage window area at low speed conditions. As shown in Figure 2, the lower-left mode stands for the typical central spillage pattern, and the upper-right mode stands for the side spillage pattern. The flat plane ABCD at the bottom represents the original incident shock of the inlet. At low speed conditions, the incident shock will detach from the compression center (e.g., points O1 and A1/B1). The spillage window then forms between the leading edge CD and the compression center (shown as the red-colored part in the figure). AOPC and BOPD are two unit-squares for simplicity. The detailed analysis is as follows: (1) the bottom line CD is identical for two patterns. This notion indicates that the areas of spillage windows are determined by the track length of compression centers Lcc (e.g., curves OO1 and AA1/BB1). (2) The hypotenuse of the spillage window (e.g., O1C and A1C) represents the streamwise span between the leading edge and the compression center. Specifically, the longer hypotenuse is led from a larger contraction ratio (i.e., a stronger incident shock). (3) In these two patterns with an identical inlet width space, the shock strength of central spillage pattern is larger than that of the side spillage pattern (O1C > A1C). The Lcc of central spillage pattern is longer than that of the side spillage pattern, which leads to better spillage ability of the latter at the wide low speed range. As analyzed above, the side-spillage pattern holds larger potentials in flow spillage ability at low speed conditions.

Modular Assessment Approach for the Spillage Pattern Analysis
The direct cause for flow spillage is the detachment of the inlet incident shock. This section employs a modular assessment approach to estimate the spillage window area at low speed conditions. As shown in Figure 2, the lower-left mode stands for the typical central spillage pattern, and the upper-right mode stands for the side spillage pattern. The flat plane ABCD at the bottom represents the original incident shock of the inlet. At low speed conditions, the incident shock will detach from the compression center (e.g., points O 1 and A 1 /B 1 ). The spillage window then forms between the leading edge CD and the compression center (shown as the red-colored part in the figure). AOPC and BOPD are two unit-squares for simplicity. The detailed analysis is as follows: (1) the bottom line CD is identical for two patterns. This notion indicates that the areas of spillage windows are determined by the track length of compression centers L cc (e.g., curves OO 1 and AA 1 /BB 1 ). (2) The hypotenuse of the spillage window (e.g., O 1 C and A 1 C) represents the streamwise span between the leading edge and the compression center. Specifically, the longer hypotenuse is led from a larger contraction ratio (i.e., a stronger incident shock).
(3) In these two patterns with an identical inlet width space, the shock strength of central spillage pattern is larger than that of the side spillage pattern (O 1 C > A 1 C). The L cc of central spillage pattern is longer than that of the side spillage pattern, which leads to better spillage ability of the latter at the wide low speed range. As analyzed above, the side-spillage pattern holds larger potentials in flow spillage ability at low speed conditions.

Novel Spillage-Adaptive Integration Pattern
The novel concept of the spillage-adaptive integration pattern is proposed based on the analysis in Sections 2.1 and 2.2. Figure 3 demonstrates that the pattern consists of the pressure-controllable bump and two symmetrical side-located inlets. The pressure-controllable bump is the transition part between the airframe and the inlets. The parabolic profile of bump leading edge is flexible to be integrated with various airframes. In contrast with the conventional streamline-tracing bump, the pressure-controllable bump is inversely generated by the prescribed pressure distribution. Its pressure distribution controls the boundary layer diversion and the flow uniformity. Two inward-turning inlets are set back-to-back along the central plane of the bump. The region after the red-ribbon area in Figure 3 is the inlet presetting position. This region is designed to be occupied with high-momentum fluid. The flow uniformity is adjusted from expanded to parallel.

Integrated Design of the Spillage-Adaptive Bump Inlet
The integrated design method of the spillage-adaptive bump inlet is presented in this section. The pressure-controllable bump is inversely designed by the prescribed surface distribution. The inlets are designed on the basis of the classical osculating plane theory. The framework of the design is shown in Figure 4.

Novel Spillage-Adaptive Integration Pattern
The novel concept of the spillage-adaptive integration pattern is proposed based on the analysis in Sections 2.1 and 2.2. Figure 3 demonstrates that the pattern consists of the pressure-controllable bump and two symmetrical side-located inlets. The pressurecontrollable bump is the transition part between the airframe and the inlets. The parabolic profile of bump leading edge is flexible to be integrated with various airframes. In contrast with the conventional streamline-tracing bump, the pressure-controllable bump is inversely generated by the prescribed pressure distribution. Its pressure distribution controls the boundary layer diversion and the flow uniformity. Two inward-turning inlets are set back-to-back along the central plane of the bump. The region after the red-ribbon area in Figure 3 is the inlet presetting position. This region is designed to be occupied with high-momentum fluid. The flow uniformity is adjusted from expanded to parallel.

Novel Spillage-Adaptive Integration Pattern
The novel concept of the spillage-adaptive integration pattern is proposed based on the analysis in Sections 2.1 and 2.2. Figure 3 demonstrates that the pattern consists of the pressure-controllable bump and two symmetrical side-located inlets. The pressure-controllable bump is the transition part between the airframe and the inlets. The parabolic profile of bump leading edge is flexible to be integrated with various airframes. In contrast with the conventional streamline-tracing bump, the pressure-controllable bump is inversely generated by the prescribed pressure distribution. Its pressure distribution controls the boundary layer diversion and the flow uniformity. Two inward-turning inlets are set back-to-back along the central plane of the bump. The region after the red-ribbon area in Figure 3 is the inlet presetting position. This region is designed to be occupied with high-momentum fluid. The flow uniformity is adjusted from expanded to parallel.

Integrated Design of the Spillage-Adaptive Bump Inlet
The integrated design method of the spillage-adaptive bump inlet is presented in this section. The pressure-controllable bump is inversely designed by the prescribed surface distribution. The inlets are designed on the basis of the classical osculating plane theory. The framework of the design is shown in Figure 4.

Integrated Design of the Spillage-Adaptive Bump Inlet
The integrated design method of the spillage-adaptive bump inlet is presented in this section. The pressure-controllable bump is inversely designed by the prescribed surface distribution. The inlets are designed on the basis of the classical osculating plane theory. The framework of the design is shown in Figure 4.  The design of the pressure-controllable bump includes the leading edge and the surface pressure distribution. The four design parameters of the bump leading edge are the cone-shock angle δshock, the intercepting height Hintercept, the bump width Wbump, and the bump length Lbump. The δshock is determined by the freestream speed and the derived cone angle, which can be calculated by the axisymmetric conical shock equations. The leading edge is obtained by the intersection of the conical shock wave and intercept plane. In this study, the planes parallel to xoy are taken as the intercept plane. The Hintercept is the height between the intercepting plane and the central line. The Wbump and Lbump are two related parameters. Herein, the bump leading edge is uniquely determined when δshock, Hintercept, and Wbump (or Lbump) are fixed.
The PR mechanism is used for the bump surface pressure distribution. The PR concept originated from the mountain ridge, which blocks the wind near the ground from blowing to the other side of the mountain. A relatively high-pressure area is also set on the bump surface. The pressure is higher than that of both sides, thus making it like a ridge of pressure. The boundary layer is diverted by the outward pressure gradient ▽out. The flow direction is adjusted from expanded to parallel by the inward pressure gradient ▽in. The πpeak is the peak value of PR, and the Wpeak is the width of PR.
The RBF-based (RBF: radius-based function) permeable boundary method [34] is previously developed to generate a bump surface with the PR pressure distribution. The classical permeable boundary method requires strict mesh topology, which is commonly applied in 2D aeronautical impeller design. The calculating efficiency and accuracy are seriously restricted in 3D surface with arbitrary pressure distribution. The previous study of authors coupled the permeable boundary method with the RBF network. The aerodynamic values of every single point on the bump surface are obtained by the interpolation of the reference points. The interpolation uses radius weighting to determine the values of solution points, which ensures the flow continuity and has better physical sense than other mathematical interpolation methods.
The inlet presetting location is determined on the basis of the bump flow field. The HKE flow ratio and the flow uniformity are two main aspects to be evaluated. As previously mentioned, the boundary layer diversion depends on ▽out and the flow uniformity depends on▽in. Variables ▽out and ▽in are designed for the corresponding boundary layer thickness and the flow expansion ratio. The schematic of the inlet presetting location is shown in Figure 5 with the translucent red part. The upper boundary is the bump shock, the lower boundary is the bump surface, and the left/right boundaries are the HKE flow regions that are within the PR. The inlet lip (i.e., the inlet entrance) should be set in this region and capture the certain amount of mass flow for the engine. The area of inlet exit The design of the pressure-controllable bump includes the leading edge and the surface pressure distribution. The four design parameters of the bump leading edge are the cone-shock angle δ shock , the intercepting height H intercept , the bump width W bump , and the bump length L bump . The δ shock is determined by the freestream speed and the derived cone angle, which can be calculated by the axisymmetric conical shock equations. The leading edge is obtained by the intersection of the conical shock wave and intercept plane. In this study, the planes parallel to xoy are taken as the intercept plane. The H intercept is the height between the intercepting plane and the central line. The W bump and L bump are two related parameters. Herein, the bump leading edge is uniquely determined when δ shock , H intercept , and W bump (or L bump ) are fixed.
The PR mechanism is used for the bump surface pressure distribution. The PR concept originated from the mountain ridge, which blocks the wind near the ground from blowing to the other side of the mountain. A relatively high-pressure area is also set on the bump surface. The pressure is higher than that of both sides, thus making it like a ridge of pressure. The boundary layer is diverted by the outward pressure gradient out . The flow direction is adjusted from expanded to parallel by the inward pressure gradient in . The π peak is the peak value of PR, and the W peak is the width of PR.
The RBF-based (RBF: radius-based function) permeable boundary method [34] is previously developed to generate a bump surface with the PR pressure distribution. The classical permeable boundary method requires strict mesh topology, which is commonly applied in 2D aeronautical impeller design. The calculating efficiency and accuracy are seriously restricted in 3D surface with arbitrary pressure distribution. The previous study of authors coupled the permeable boundary method with the RBF network. The aerodynamic values of every single point on the bump surface are obtained by the interpolation of the reference points. The interpolation uses radius weighting to determine the values of solution points, which ensures the flow continuity and has better physical sense than other mathematical interpolation methods.
The inlet presetting location is determined on the basis of the bump flow field. The HKE flow ratio and the flow uniformity are two main aspects to be evaluated. As previously mentioned, the boundary layer diversion depends on out and the flow uniformity depends on in . Variables out and in are designed for the corresponding boundary layer thickness and the flow expansion ratio. The schematic of the inlet presetting location is shown in Figure 5 with the translucent red part. The upper boundary is the bump shock, the lower boundary is the bump surface, and the left/right boundaries are the HKE flow regions that are within the PR. The inlet lip (i.e., the inlet entrance) should be set in this region and capture the certain amount of mass flow for the engine. The area of inlet exit is is calculated by the streamwise projected area of inlet lip and the contraction ratio. While the shape of inlet exit is assigned by the combustion chamber. A coefficient RLKE to evaluate the low kinetic energy (LKE) flow ratio of the reference inlet entrance is proposed in this study. This factor is the equivalent LKE flow area divided by the area between shock and bump, that indicates the usage of the inlet presetting region. The definition is as follows: where W is the prescribed width of inlet, yshock(z) is the y-coordinate of shock position at a specific z-coordinate, ybump(z) is the y-coordinate of bump surface at a specific z-coordinate, HKE is the equivalent thickness of the kinetic energy loss, which converts the kinetic energy loss to an equivalent thickness in the inviscid condition. This factor is defined as: where ρ, u are density and x-direction speed and ρ∞, u∞ are the parameters of freestream. The ηKE is the kinetic energy efficiency and is defined as: where γ is the specific heat ratio, T * is the temperature, σ is the total pressure recovery. RLKE distribution at the inlet presetting region is shown in Figure 6. The abscissa is the half width of the inlet. The bump width is 0.698 m. The RLKE is lower than 0.0113 when W < 0.2 m, which indicates that the boundary layer is well diverted, and the HKE occupies most of the region. The RLKE is larger than 0.0113 when W > 0.2 m, which indicates that the boundary layer diversion ratio is small. The distribution of RLKE is flat when W < 0.2 m, but it rapidly increases when W > 0.2 m. The inlet width is finally determined as 0.4 m. A coefficient R LKE to evaluate the low kinetic energy (LKE) flow ratio of the reference inlet entrance is proposed in this study. This factor is the equivalent LKE flow area divided by the area between shock and bump, that indicates the usage of the inlet presetting region. The definition is as follows: where W is the prescribed width of inlet, y shock (z) is the y-coordinate of shock position at a specific z-coordinate, y bump (z) is the y-coordinate of bump surface at a specific z-coordinate, H KE is the equivalent thickness of the kinetic energy loss, which converts the kinetic energy loss to an equivalent thickness in the inviscid condition. This factor is defined as: where ρ, u are density and x-direction speed and ρ ∞ , u ∞ are the parameters of freestream. The η KE is the kinetic energy efficiency and is defined as: where γ is the specific heat ratio, T * is the temperature, σ is the total pressure recovery. R LKE distribution at the inlet presetting region is shown in Figure 6. The abscissa is the half width of the inlet. The bump width is 0.698 m. The R LKE is lower than 0.0113 when W < 0.2 m, which indicates that the boundary layer is well diverted, and the HKE occupies most of the region. The R LKE is larger than 0.0113 when W > 0.2 m, which indicates that the boundary layer diversion ratio is small. The distribution of R LKE is flat when W < 0.2 m, but it rapidly increases when W > 0.2 m. The inlet width is finally determined as 0.4 m. With regard to the inlet subjected to the uniform inflow, the arrangement of osculating planes is the key factor to the inlet spillage at lower speed. The pre-compression strength of the inlet integrated with the forebody is another crucial aspect to be concerned about. The vertical compressing component is larger than the horizontal compressing component due to the elliptical shape of the bump shock. The incident shock detaches from the inlet lip and further distorts while the inflow speed decreases. The detached shock will vertically interact with the bump shock first and then horizontally. Hence, the side spillage pattern of the identical inlet contraction ratio is preferable to the bump inlet integration, which confirms better spillage ability in the low speeds. The spillage ratio and the mass capturing ratio can be balanced by different side spillage locations on the inlet lip.
The projection of the entrance and exit for the new spillage-adaptive inlet is shown in Figure 7. Line AB is the connecting part of two inlets. The solid black line is the projection of the inlet entrance on the x-direction. Meanwhile, the solid red line is the projection of the inlet exit. Point O is the compression center of the inlet. The ratio of projection length between AB and BC is 1.02, which confirms that the inlet lip is within the presetting location ( Figure 5). The contraction ratio of the projected area of the entrance and exit is 3.65, of which the contraction ratio from point A to A′ is 2.80, the contraction ratio from point B to B′ is 2.07, and the contraction ratio from point C to C′ is 2.17. The contraction ratio of typical positions AA′, BB′, and CC′ is relatively uniform, which can weaken the lateral compression wave intensity between the rear section of the compression section of the inlet to a certain extent. The projection of exit profile is circular, and the AOC is a spline. The entrance curve AOC is tangent to the outlet curve A′OC′ at point O. Thus, the inlet will not produce additional external drag. With regard to the inlet subjected to the uniform inflow, the arrangement of osculating planes is the key factor to the inlet spillage at lower speed. The pre-compression strength of the inlet integrated with the forebody is another crucial aspect to be concerned about. The vertical compressing component is larger than the horizontal compressing component due to the elliptical shape of the bump shock. The incident shock detaches from the inlet lip and further distorts while the inflow speed decreases. The detached shock will vertically interact with the bump shock first and then horizontally. Hence, the side spillage pattern of the identical inlet contraction ratio is preferable to the bump inlet integration, which confirms better spillage ability in the low speeds. The spillage ratio and the mass capturing ratio can be balanced by different side spillage locations on the inlet lip.
The projection of the entrance and exit for the new spillage-adaptive inlet is shown in Figure 7. Line AB is the connecting part of two inlets. The solid black line is the projection of the inlet entrance on the x-direction. Meanwhile, the solid red line is the projection of the inlet exit. Point O is the compression center of the inlet. The ratio of projection length between AB and BC is 1.02, which confirms that the inlet lip is within the presetting location ( Figure 5). The contraction ratio of the projected area of the entrance and exit is 3.65, of which the contraction ratio from point A to A is 2.80, the contraction ratio from point B to B is 2.07, and the contraction ratio from point C to C is 2.17. The contraction ratio of typical positions AA , BB , and CC is relatively uniform, which can weaken the lateral compression wave intensity between the rear section of the compression section of the inlet to a certain extent. The projection of exit profile is circular, and the AOC is a spline. The entrance curve AOC is tangent to the outlet curve A OC at point O. Thus, the inlet will not produce additional external drag.  The high external compression ratio basic flow field for designing the IWR inlet is developed in the authors' previous study [16], which has proved its efficiency for M 3.5 inflow condition. This study uses the same approach to improve the flow field in M 6.0 condition. The typical internal conical flow field (ICFC) is shown in the left subplot of Figure 8. A large expansion zone can be observed at the entrance of the compression section. The Mach number increases from 5.2 to 5.3 and then decreases to 5.2, which brings an additional oblique shock wave at x = 4 in the basic flow field. On the one hand, the interaction between shocks is introduced to increase the energy loss. On the other hand, the expanded flow reduces the compression efficiency. Therefore, the flow field length required to obtain the same Mach number at the exit is larger. In contrast with the basic flow field of ICFC, the basic flow field of the new ICFC+ achieves the following improvements: (1) the ICFC+ flow field improves the intensity of the upstream initial shock wave. Accordingly, the incident point of the initial shock wave is moved from 3.6 to 3.4, and the Mach number after the wave is reduced from the original 5.3 to 5.1; (2) in the initial section of the ICFC+ flow field, the expansion zone has been eliminated, and the initial section of the basic flow field is from the original two channels. The incident shock is simplified as an incident shock. The Mach number distribution along the path demonstrates that the Mach number of the wall corresponding to x coordinate of the ICFC+ flow field is at least 0.1 lower than that of the ICFC flow field, maintaining a good compression efficiency; (3) the length of the flow field is five for the ICFC flow field to compress the flow from M 6.0 to 3.2, while the length of the flow field is 4.6 for the ICFC+ flow field to compress the flow from M 6.0 to 3.1. The basic flow field is shortened by 8%, which indicates that the new ICFC+ flow field greatly improves the compression efficiency by eliminating the upstream expansion zone and shortens the flow field length required to achieve the same compression performance. The configuration of the new spillage-adaptive bump inlet is shown in Figure 9. The two inlet channels of the double channel lateral overflow scheme are symmetrically The high external compression ratio basic flow field for designing the IWR inlet is developed in the authors' previous study [16], which has proved its efficiency for M 3.5 inflow condition. This study uses the same approach to improve the flow field in M 6.0 condition. The typical internal conical flow field (ICFC) is shown in the left subplot of Figure 8. A large expansion zone can be observed at the entrance of the compression section. The Mach number increases from 5.2 to 5.3 and then decreases to 5.2, which brings an additional oblique shock wave at x = 4 in the basic flow field. On the one hand, the interaction between shocks is introduced to increase the energy loss. On the other hand, the expanded flow reduces the compression efficiency. Therefore, the flow field length required to obtain the same Mach number at the exit is larger. In contrast with the basic flow field of ICFC, the basic flow field of the new ICFC+ achieves the following improvements: (1) the ICFC+ flow field improves the intensity of the upstream initial shock wave. Accordingly, the incident point of the initial shock wave is moved from 3.6 to 3.4, and the Mach number after the wave is reduced from the original 5.3 to 5.1; (2) in the initial section of the ICFC+ flow field, the expansion zone has been eliminated, and the initial section of the basic flow field is from the original two channels. The incident shock is simplified as an incident shock. The Mach number distribution along the path demonstrates that the Mach number of the wall corresponding to x coordinate of the ICFC+ flow field is at least 0.1 lower than that of the ICFC flow field, maintaining a good compression efficiency; (3) the length of the flow field is five for the ICFC flow field to compress the flow from M 6.0 to 3.2, while the length of the flow field is 4.6 for the ICFC+ flow field to compress the flow from M 6.0 to 3.1. The basic flow field is shortened by 8%, which indicates that the new ICFC+ flow field greatly improves the compression efficiency by eliminating the upstream expansion zone and shortens the flow field length required to achieve the same compression performance.  The high external compression ratio basic flow field for designing the IWR inlet is developed in the authors' previous study [16], which has proved its efficiency for M 3.5 inflow condition. This study uses the same approach to improve the flow field in M 6.0 condition. The typical internal conical flow field (ICFC) is shown in the left subplot of Figure 8. A large expansion zone can be observed at the entrance of the compression section. The Mach number increases from 5.2 to 5.3 and then decreases to 5.2, which brings an additional oblique shock wave at x = 4 in the basic flow field. On the one hand, the interaction between shocks is introduced to increase the energy loss. On the other hand, the expanded flow reduces the compression efficiency. Therefore, the flow field length required to obtain the same Mach number at the exit is larger. In contrast with the basic flow field of ICFC, the basic flow field of the new ICFC+ achieves the following improvements: (1) the ICFC+ flow field improves the intensity of the upstream initial shock wave. Accordingly, the incident point of the initial shock wave is moved from 3.6 to 3.4, and the Mach number after the wave is reduced from the original 5.3 to 5.1; (2) in the initial section of the ICFC+ flow field, the expansion zone has been eliminated, and the initial section of the basic flow field is from the original two channels. The incident shock is simplified as an incident shock. The Mach number distribution along the path demonstrates that the Mach number of the wall corresponding to x coordinate of the ICFC+ flow field is at least 0.1 lower than that of the ICFC flow field, maintaining a good compression efficiency; (3) the length of the flow field is five for the ICFC flow field to compress the flow from M 6.0 to 3.2, while the length of the flow field is 4.6 for the ICFC+ flow field to compress the flow from M 6.0 to 3.1. The basic flow field is shortened by 8%, which indicates that the new ICFC+ flow field greatly improves the compression efficiency by eliminating the upstream expansion zone and shortens the flow field length required to achieve the same compression performance. The configuration of the new spillage-adaptive bump inlet is shown in Figure 9. The two inlet channels of the double channel lateral overflow scheme are symmetrically The configuration of the new spillage-adaptive bump inlet is shown in Figure 9. The two inlet channels of the double channel lateral overflow scheme are symmetrically placed back-to-back, and the inlet profile is completely located within the PR region. The placed back-to-back, and the inlet profile is completely located within the PR region. The dimensions of the bump inlet are 2371 mm × 1161 mm × 265 mm. The shape of the bump inlet is connected with the flat plate in this work to simplify the flow field analysis with low energy flow and study the performance of the inlet with different boundary layer thicknesses.

Preparation of Numerical Simulation
The design condition is M 6.0 at altitude H = 24 km. The aerodynamic performance of the bump inlet was examined by viscous numerical simulations. The flow was modeled as the compressible steady flow using Reynolds-Averaged Navier-Stokes (RANS) equations simulated by the ANSYS CFX commercial software. The governing equations are demonstrated in the following form: where U is the conservative vector, Fc is the convective flux vector, Fv is the viscous flux vector, n = nxi + nyj + nzk is the unit normal vector, t is the time, Ω is the control volume, and S is the control volume surface. The vectors U, Fc, and Fv are defined as: In the above equations, the ρ is the density. The u, v, and w are respectively the x/y/z partial velocities on Cartesian coordinate system. The E, p, and τ are respectively the total energy per unit mass, the pressure, the viscous stress tensor. The Πx~Πz are given by: where the T is the temperature and the λ is the thermal conductivity.

Preparation of Numerical Simulation
The design condition is M 6.0 at altitude H = 24 km. The aerodynamic performance of the bump inlet was examined by viscous numerical simulations. The flow was modeled as the compressible steady flow using Reynolds-Averaged Navier-Stokes (RANS) equations simulated by the ANSYS CFX commercial software. The governing equations are demonstrated in the following form: where U is the conservative vector, F c is the convective flux vector, F v is the viscous flux vector, n = n x i + n y j + n z k is the unit normal vector, t is the time, Ω is the control volume, and S is the control volume surface. The vectors U, F c , and F v are defined as: In the above equations, the ρ is the density. The u, v, and w are respectively the x/y/z partial velocities on Cartesian coordinate system. The E, p, and τ are respectively the total energy per unit mass, the pressure, the viscous stress tensor. The Π x~Πz are given by: where the T is the temperature and the λ is the thermal conductivity. The high-resolution scheme based on MUSCL (Monotonic Upstream-centered Scheme for Conservation Laws) reconstruction was employed to calculate the flux terms. That is, all results are calculated by CFX commercial software 2019 R3 based on high resolution in numeric turbulence and advection scheme. The calculation converges when the residuals are less than 1 × 10 −5 .
The simulation for mixed external-internal flow field is a special case to the turbulence model determination. Zhang et al. [35] studied the impacts of turbulence models on the numerical results. Both k-ε and Spalart-Allmaras turbulence models match well with experiment data of a 2D supersonic inlet, which contains external and internal flow field. The k-ε turbulence model was used in this study while considering both accuracy and efficiency. Besides, a previous research [36] obtained that in the simulation of hypersonic inlets, the calculation results branch between using variable and constant specific heat ratio at high speeds (approximately M 4.0), as shown in Figure 10. The difference between variable/fixed specific heat ratio cases can be ignored when M < 4.0. When M ≥ 4.0, the difference increases with M. The total temperature of fixed specific heat ratio is approximately 18% higher than that of variable specific heat ratio. Therefore, this study uses the variable specific heat ratio in the CFD simulation. Thus, the piecewise-polynomial method is chosen to calculate the specific heat and viscosity by the Sutherland law.
The high-resolution scheme based on MUSCL (Monotonic Upstream-centered Scheme for Conservation Laws) reconstruction was employed to calculate the flux terms. That is, all results are calculated by CFX commercial software 2019 R3 based on high resolution in numeric turbulence and advection scheme. The calculation converges when the residuals are less than 1 × 10 −5 .
The simulation for mixed external-internal flow field is a special case to the turbulence model determination. Zhang et al. [35] studied the impacts of turbulence models on the numerical results. Both k-ε and Spalart-Allmaras turbulence models match well with experiment data of a 2D supersonic inlet, which contains external and internal flow field. The k-ε turbulence model was used in this study while considering both accuracy and efficiency. Besides, a previous research [36] obtained that in the simulation of hypersonic inlets, the calculation results branch between using variable and constant specific heat ratio at high speeds (approximately M 4.0), as shown in Figure 10. The difference between variable/fixed specific heat ratio cases can be ignored when M < 4.0. When M ≥ 4.0, the difference increases with M. The total temperature of fixed specific heat ratio is approximately 18% higher than that of variable specific heat ratio. Therefore, this study uses the variable specific heat ratio in the CFD simulation. Thus, the piecewise-polynomial method is chosen to calculate the specific heat and viscosity by the Sutherland law. The 3D structured mesh is established for the numerical simulation, which is shown in Figure 11. A 3.5 m flat plate is set in front of the bump surface for boundary layer development, which is divided by four parts to create boundary layer with different thicknesses. The length of the four parts is successively 0.5, 1.0, 1.0, 1.0 m. That is, if the front two (purple and red) plates are set as free-slip wall condition and the rear two (yellow and green) plates are set as no-slip wall condition, the actual length for boundary layer development (LBL) is 2.0 m. Thereby, it can create inflow conditions with various nominal boundary layer thicknesses as 9.00, 21.05, 32.37, 41.72 mm at M 6.0, which are successively 5.1%, 11.8%, 18.2%, 23.4% of the inlet entrance height. To obtain the requested standard k-ε turbulence model, the minimum mesh height of near-wall grid is 0.15 mm, and the grids in the y-direction are stretched with the increasing ratio of 1.2 refined with geometric proportion rule. The y+ is kept as 20~110 to meet the standard k-ε turbulence model. A Cgrid is set to surround the bump, which accurately fits the bump leading edge. An O-grid is set in the internal flow field of inlet, which accurately fits the inlet wall. A grid independence study is conducted for the bump inlet configuration. The mesh ranging between 3.5 and 6.0 million is generated. The results are shown in Table 1. The σ is the total pressure recovery, the π is the static pressure ratio, which is calculated by the exit pressure divided by the freestream pressure. The φ is the mass capturing ratio. All the aerodynamic The 3D structured mesh is established for the numerical simulation, which is shown in Figure 11. A 3.5 m flat plate is set in front of the bump surface for boundary layer development, which is divided by four parts to create boundary layer with different thicknesses. The length of the four parts is successively 0.5, 1.0, 1.0, 1.0 m. That is, if the front two (purple and red) plates are set as free-slip wall condition and the rear two (yellow and green) plates are set as no-slip wall condition, the actual length for boundary layer development (L BL ) is 2.0 m. Thereby, it can create inflow conditions with various nominal boundary layer thicknesses as 9.00, 21.05, 32.37, 41.72 mm at M 6.0, which are successively 5.1%, 11.8%, 18.2%, 23.4% of the inlet entrance height. To obtain the requested standard k-ε turbulence model, the minimum mesh height of near-wall grid is 0.15 mm, and the grids in the y-direction are stretched with the increasing ratio of 1.2 refined with geometric proportion rule. The y+ is kept as 20~110 to meet the standard k-ε turbulence model. A C-grid is set to surround the bump, which accurately fits the bump leading edge. An O-grid is set in the internal flow field of inlet, which accurately fits the inlet wall. A grid independence study is conducted for the bump inlet configuration. The mesh ranging between 3.5 and 6.0 million is generated. The results are shown in Table 1. The σ is the total pressure recovery, the π is the static pressure ratio, which is calculated by the exit pressure divided by the freestream pressure. The ϕ is the mass capturing ratio. All the aerodynamic characteristics in the current work are obtained by mass-averaged method on corresponding position (e.g., the freestream, the throat, the inlet exit, etc.,). Due to the aerodynamic characteristics discrepancy between different meshes, the reasonable mesh in the domain is approximately 4.3 million.  The boundary settings for calculation are as follows. The no-slip and adiabatic wall conditions are adopted for entire walls, which include the bump inlet and the flat plate. The symmetry condition is adopted for the central plane and the freestream condition is adopted for the outer part of the whole calculating domain. The supersonic outlet condition is adopted for exit of the inlet. The numerical methods in this study have been effectively validated by You and Liang [37], who used wind tunnel experiments and showed that the method can calculate reasonable and reliable results. The numerical method has also been validated in the previous study [5] through a comparison with an experimental study of bump-inlet integration shown in Figure 12. The results show that the difference of static pressure between the numerical and the experimental results is less than 2.7, which indicates that the two results agree with each other. As validated above, the numerical methods effectively calculate the true flow field.  The boundary settings for calculation are as follows. The no-slip and adiabatic wall conditions are adopted for entire walls, which include the bump inlet and the flat plate. The symmetry condition is adopted for the central plane and the freestream condition is adopted for the outer part of the whole calculating domain. The supersonic outlet condition is adopted for exit of the inlet. The numerical methods in this study have been effectively validated by You and Liang [37], who used wind tunnel experiments and showed that the method can calculate reasonable and reliable results. The numerical method has also been validated in the previous study [5] through a comparison with an experimental study of bump-inlet integration shown in Figure 12. The results show that the difference of static pressure between the numerical and the experimental results is less than 2.7, which indicates that the two results agree with each other. As validated above, the numerical methods effectively calculate the true flow field.

Aerodynamic Characteristics of the Spillage-adaptive Integration at the Design Condition
The aerodynamic characteristics of the new integration (pattern 1) are shown in Table  2. The results of the typical central spillage pattern (pattern 2) are also presented for comparison. The mass capturing ratio φ1 is 0.4% higher than φ2. The Mexit,1 is 3.6% lower than Mexit,2. The total pressure recovery σ1 is 2.0% higher than σ2. The exit static pressure π1 is 7.1% lower than π2. The RM is the ratio of the deceleration, which is Mexit divided by M∞. The kinetic energy efficiency ηKE,1 is almost the same with ηKE,2. The two integration pat-

Aerodynamic Characteristics of the Spillage-adaptive Integration at the Design Condition
The aerodynamic characteristics of the new integration (pattern 1) are shown in Table 2. The results of the typical central spillage pattern (pattern 2) are also presented for comparison. The mass capturing ratio ϕ 1 is 0.4% higher than ϕ 2 . The M exit,1 is 3.6% lower than M exit,2 . The total pressure recovery σ 1 is 2.0% higher than σ 2 . The exit static pressure π 1 is 7.1% lower than π 2 . The R M is the ratio of the deceleration, which is M exit divided by M ∞ . The kinetic energy efficiency η KE,1 is almost the same with η KE,2 . The two integration patterns maintain similar performance in the design condition. The flow structure of the new spillage-adaptive bump inlet is shown in Figure 13. The pressure contours of the bump, cowl lip, and inlet are presented. The M contours of various crosswise-sections along the x-direction are also obtained. The detailed analysis is as follows. (1) The boundary layer develops along the x-direction, and its thickness increases. Once the layer is encountered with the bump surface, it is diverted by out and accumulates at the corner between the bump and the flat plate. (2) The inlets are located after the bump shock; thus, no interaction exists between the bump shock and the incident shock at design condition. Meanwhile, the LKE flow affects the two shocks. The bump shock and the incident shock distort by the influence of the LKE flow. The bump shock has not attached well to the leading edge. The incident shock distorts upstream, which decreases the mass capturing ratio. (3) The PR mechanism and the side spillage pattern demonstrate that the streamlines at the inlet bottom wall are uniform. No concentrating trend can be observed for the streamlines. This notion indicates that the flow field of inlet is uniform to avoid overlarge side compression. Specifically, the SES effects (Section 2.1) are avoided in the new bump inlet.

Start-Ability of the New Spillage-Adaptive Bump Inlet
The new bump inlet is numerically studied at M 2.0-6.0 to investigate the speed range. In Table 3, the swallowed mass has not exceeded the inlet flow capacity. Meanwhile, the compression ratio of M 2.0 is 3.6, which is 56% lower than that of M 2.5 condition. This value is particularly low to bear the pressure fluctuation from the downstream combustor. Thus, the inlet speed range is M 2.5-6.0. The incident shock distorts toward the lip, and the flow spillage increases with the decrease in inflow speed. The compression efficiency of the inlet decreases; hence, the πexit decreases. At low speeds, the actual capture area of the inlet decreases with the increase in spillage. The actual contraction ratio of the inlet decreases which results in reduction in πexit. The shock system is affected by the flow spillage, and variables σ and ηKE decrease.

Start-Ability of the New Spillage-Adaptive Bump Inlet
The new bump inlet is numerically studied at M 2.0-6.0 to investigate the speed range. In Table 3, the swallowed mass has not exceeded the inlet flow capacity. Meanwhile, the compression ratio of M 2.0 is 3.6, which is 56% lower than that of M 2.5 condition. This value is particularly low to bear the pressure fluctuation from the downstream combustor. Thus, the inlet speed range is M 2.5-6.0. The incident shock distorts toward the lip, and the flow spillage increases with the decrease in inflow speed. The compression efficiency of the inlet decreases; hence, the π exit decreases. At low speeds, the actual capture area of the inlet decreases with the increase in spillage. The actual contraction ratio of the inlet decreases which results in reduction in π exit . The shock system is affected by the flow spillage, and variables σ and η KE decrease. The spillage ratio ϕ spillage is used in this work to better analyze the influence of flow spillage on the inlet start-ability. This concept is defined as the spilled mass flow of the low speed condition ( The ϕ spillage -M ∞ distribution is shown in Figure 14. The analysis is as follows: (1) the relationship between two parameters is quasi-parabolic.   Figure 15, which shows the interaction between the incident shock and the thick boundary layer (LBL = 3.5 m). The detailed analysis is as follows: 1) the shock system of M∞ = 5.0 is similar to that of M∞ = 6.0 (shown in Figure 13). Flow spillage (9.5%) occurs at the upper part of the inlet lip. 2) When M∞ decreases to 4.0, the spillage increases by 26.8% than the  Figure 15, which shows the interaction between the incident shock and the thick boundary layer (L BL = 3.5 m). The detailed analysis is as follows: 1) the shock system of M ∞ = 5.0 is similar to that of M ∞ = 6.0 (shown in Figure 13). Flow spillage (9.5%) occurs at the upper part of the inlet lip. 2) When M ∞ decreases to 4.0, the spillage increases by 26.8% than the design condition. The incident shock distorts the upstream to make a spillage window at the compression center and the side wall of the inlet. The LKE fluid accumulates around the side wall. This phenomenon restricts the flow capacity by forming an aerodynamic throat that is smaller in area than the geometry. 3) The decrease in M ∞ leads to the main flow spilling from the adjacent region of the compression center. Two impact patterns of spillage to the flow field can be observed. In the speed range of M 3.0 to 5.0, the spillage ratio is lower than 40%. The spilled flow has not encountered bump shock. In the speed range of M 2.0 to 3.0, the spillage ratio is over 40%. The interaction between the inlet incident shock and the bump shock leads to the change in the external flow field, which can be observed on the middle section (colored in red).  Figure 15, which shows the interaction between the incident shock and the thick boundary layer (LBL = 3.5 m). The detailed analysis is as follows: 1) the shock system of M∞ = 5.0 is similar to that of M∞ = 6.0 (shown in Figure 13). Flow spillage (9.5%) occurs at the upper part of the inlet lip. 2) When M∞ decreases to 4.0, the spillage increases by 26.8% than the design condition. The incident shock distorts the upstream to make a spillage window at the compression center and the side wall of the inlet. The LKE fluid accumulates around the side wall. This phenomenon restricts the flow capacity by forming an aerodynamic throat that is smaller in area than the geometry. 3) The decrease in M∞ leads to the main flow spilling from the adjacent region of the compression center. Two impact patterns of spillage to the flow field can be observed. In the speed range of M 3.0 to 5.0, the spillage ratio is lower than 40%. The spilled flow has not encountered bump shock. In the speed range of M 2.0 to 3.0, the spillage ratio is over 40%. The interaction between the inlet incident shock and the bump shock leads to the change in the external flow field, which can be observed on the middle section (colored in red).

Effects of Boundary Layer Thickness to Aerodynamic Characteristics
The aerodynamic characteristics of the inlet under different boundary layer thicknesses are shown in Figure 16. The LBL is the length for boundary layer development. When the LBL decreases from 3.5 to 1.0 m, the LKE ratio of the captured inflow decreases. Accordingly, the Mexit and σ increase, and the gradient of these two parameters decrease with the LBL. This phenomenon can be explained on the basis of the boundary layer thickness. The nominal boundary layer thickness δ (i.e., the height where 99% of mainstream speed locates) is positively correlated to the development length of boundary layer, which is shown as,

Effects of Boundary Layer Thickness to Aerodynamic Characteristics
The aerodynamic characteristics of the inlet under different boundary layer thicknesses are shown in Figure 16. The L BL is the length for boundary layer development. When the L BL decreases from 3.5 to 1.0 m, the LKE ratio of the captured inflow decreases. Accordingly, the M exit and σ increase, and the gradient of these two parameters decrease with the L BL . This phenomenon can be explained on the basis of the boundary layer thickness. The nominal boundary layer thickness δ (i.e., the height where 99% of mainstream speed locates) is positively correlated to the development length of boundary layer, which is shown as, where ν is the kinetic viscosity coefficient, x is the length for boundary layer development, U is the mainstream speed. When ν and U are fixed, x and δ are in parabolic relation. Variable x stands for L BL in this study, and δ denotes the boundary layer thickness. Considering 0 < δ < 1, the decrease in x will lead to the decrease in δ. When L BL decreases with the identical ∆x, the σ will increase, but the ∆σ will decrease. This condition explains the trend in Figure 16b. The Mach number will increase with the decrease in L BL , but ∆M will decrease, which can explain the trend in Figure 16a.

Effects of Boundary Layer Thickness to Aerodynamic Characteristics
The aerodynamic characteristics of the inlet under different boundary layer thicknesses are shown in Figure 16. The LBL is the length for boundary layer development. When the LBL decreases from 3.5 to 1.0 m, the LKE ratio of the captured inflow decreases. Accordingly, the Mexit and σ increase, and the gradient of these two parameters decrease with the LBL. This phenomenon can be explained on the basis of the boundary layer thickness. The nominal boundary layer thickness δ (i.e., the height where 99% of mainstream speed locates) is positively correlated to the development length of boundary layer, which is shown as, where ν is the kinetic viscosity coefficient, x is the length for boundary layer development, U is the mainstream speed. When ν and U are fixed, x and δ are in parabolic relation. Variable x stands for LBL in this study, and δ denotes the boundary layer thickness. Considering 0 < δ < 1, the decrease in x will lead to the decrease in δ. When LBL decreases with the identical ∆x, the σ will increase, but the ∆σ will decrease. This condition explains the trend in Figure 16b. The Mach number will increase with the decrease in LBL, but ∆M will decrease, which can explain the trend in Figure 16a. According to Figure 16c,d, σ and ηKE increase to a certain extent. When LBL decreases from 3.5 to 2.0 m, the mass capturing ratio φ only increases by 2.6%, while it increases by 13.7% when LBL decreases from 2.0 to 1.0 m. This notion indicates that the incident shock is distorted by the disturbance of the boundary layer at LBL > 2.0 m. Meanwhile, the disturbance at LBL < 2.0 m is not strong enough to induce shock distortion. The ηKE is positively related to the LBL. According to the definition of ηKE (Equation (3)), a highly nonlinear relationship exists between ηKE and LBL. A clear trend in the ηKE-LBL distribution is difficult to obtain. The ηKE-RM distribution is presented in the following analysis to evaluate the overall compressing efficiency of inlet. The inlet efficiency is shown in Figure 17. The hypersonic inlet expert Van Wie [38] proposed to use the kinetic efficiency ηKE for the inlet performance evaluation. The definition of ηKE is shown in Equation (3). This expression includes the kinetic energy of the inlet exit and the transferred internal energy caused by the flow expansion. The ηKE-RM distribution is used to show the inlet efficiency. Van wie [38][39][40][41] proposed the ηKE-RM relation below, which is based on 2D inlet experiments.
Kouichiro [42] proposed another ηKE-RM relation, which is based on side-compression According to Figure 16c,d, σ and η KE increase to a certain extent. When L BL decreases from 3.5 to 2.0 m, the mass capturing ratio ϕ only increases by 2.6%, while it increases by 13.7% when L BL decreases from 2.0 to 1.0 m. This notion indicates that the incident shock is distorted by the disturbance of the boundary layer at L BL > 2.0 m. Meanwhile, the disturbance at L BL < 2.0 m is not strong enough to induce shock distortion. The η KE is positively related to the L BL . According to the definition of η KE (Equation (3)), a highly nonlinear relationship exists between η KE and L BL . A clear trend in the η KE -L BL distribution is difficult to obtain. The η KE -R M distribution is presented in the following analysis to evaluate the overall compressing efficiency of inlet. The inlet efficiency is shown in Figure 17. The hypersonic inlet expert Van Wie [38] proposed to use the kinetic efficiency η KE for the inlet performance evaluation. The definition of η KE is shown in Equation (3). This expression includes the kinetic energy of the inlet exit and the transferred internal energy caused by the flow expansion. The η KE -R M distribution is used to show the inlet efficiency. Van wie [38][39][40][41] proposed the η KE -R M relation below, which is based on 2D inlet experiments.

Conclusions
A novel side-spillage pattern for hypersonic bump inlet integration was proposed in this work to improve the start-ability at low-speed conditions. The conclusions are presented as follows.
The side-embedded shock effects were analyzed and a modular assessment was employed. It proved that the side-spillage pattern holds better flow spillage potentials than the typical central spillage pattern. The novel spillage-adaptive integration pattern was proposed based on the above analysis. Two side-located inlets are set back-to-back on a bump surface. The compression centers of inlets are located at two sides over the bump. The inlet presetting region is determined by four parameters: the bump configuration, the bump shock profile, the low-kinetic-energy (LKE) flow ratio, and the lateral flow deflection.
The numerical results show that the new bump inlet has similar aerodynamic characteristics to the typical central spillage bump inlet. The speed range (starting at M 2.5-6.0) of the inlet is significantly wider than the typical integration pattern (starting at M 5.0-6.0). The relationship between the spillage ratio and the inflow speed is quasi-parabolic. The spillage ability increases with the decrease of the inflow speed.
The aerodynamic characteristics of inlets under various boundary layer ratios are obtained. The total pressure recovery (σ) and mass capturing ratio (φ) are sensitive to the thick boundary layer. When LBL > 2.0 m, the σ is rapidly increased with the decrease in LBL while at LBL < 2.0 m, the φ is rapidly increased with the decrease in LBL. The analysis on the inlet efficiency indicated that the performance of the new spillage-adaptive bump inlet (results with thick boundary layer) is close to that of the side compression inlet (results without LKE inflow).  Kouichiro [42] proposed another η KE -R M relation, which is based on side-compression inlet experiments. The relation is shown as follows: This work uses the above two relations as the standards to evaluate inlet efficiency. In Figure 17, the solid blue line stands for the η KE (1) relation, and the red dashed line indicates the η KE (2) relation. The higher inlet efficiency locates at the upper-right side of the figure. The lower inlet efficiency locates at the lower-left side of the figure. When the L BL decreases from 3.5 to 2.0 m, the η KE -R M distribution is almost parallel to the η KE (1) and η KE (2) curves, which indicates that the inlet efficiency remains steady in this L BL range. When the L BL decreases from 2.0 to 1.0 m, the η KE -R M distribution approaches the η KE (2) curve. This notion indicates that the inlet efficiency starts to increase in this L BL range.

Conclusions
A novel side-spillage pattern for hypersonic bump inlet integration was proposed in this work to improve the start-ability at low-speed conditions. The conclusions are presented as follows.
The side-embedded shock effects were analyzed and a modular assessment was employed. It proved that the side-spillage pattern holds better flow spillage potentials than the typical central spillage pattern. The novel spillage-adaptive integration pattern was proposed based on the above analysis. Two side-located inlets are set back-to-back on a bump surface. The compression centers of inlets are located at two sides over the bump. The inlet presetting region is determined by four parameters: the bump configuration, the bump shock profile, the low-kinetic-energy (LKE) flow ratio, and the lateral flow deflection.
The numerical results show that the new bump inlet has similar aerodynamic characteristics to the typical central spillage bump inlet. The speed range (starting at M 2.5-6.0) of the inlet is significantly wider than the typical integration pattern (starting at M 5.0-6.0). The relationship between the spillage ratio and the inflow speed is quasi-parabolic. The spillage ability increases with the decrease of the inflow speed.
The aerodynamic characteristics of inlets under various boundary layer ratios are obtained. The total pressure recovery (σ) and mass capturing ratio (ϕ) are sensitive to the thick boundary layer. When L BL > 2.0 m, the σ is rapidly increased with the decrease in L BL while at L BL < 2.0 m, the ϕ is rapidly increased with the decrease in L BL. The analysis on the inlet efficiency indicated that the performance of the new spillage-adaptive bump inlet (results with thick boundary layer) is close to that of the side compression inlet (results without LKE inflow).