Next Article in Journal
Numerical Analysis Results of Debonding Damage Effects for an SHM System Application on a Typical Composite Beam
Next Article in Special Issue
Agglomerate Size Evolution in Solid Propellant Combustion under High Pressure
Previous Article in Journal
A Multilayered and Multifactorial Health Assessment Method for Launch Vehicle Engine under Vibration Conditions
Previous Article in Special Issue
Laser-Induced Ignition and Combustion of Single Micron-Sized Al-Li Alloy Particles in High Pressure Air/N2
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Flow Coefficient and Starting Performance Prediction of Variable Geometry Curved Axisymmetric Inlet

1
College of Aircraft Engineering, Nanchang Hangkong University, Nanchang 330063, China
2
College of Energy and Power Engineering, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China
*
Author to whom correspondence should be addressed.
Aerospace 2023, 10(6), 506; https://doi.org/10.3390/aerospace10060506
Submission received: 4 May 2023 / Revised: 24 May 2023 / Accepted: 25 May 2023 / Published: 27 May 2023
(This article belongs to the Special Issue Combustion Evaluation and Control of Solid Rocket Motors)

Abstract

:
With the development of combined cycle engines, it is urgent to estimate more quickly and accurately the flow capture capacity and starting performance of variable geometry inlets over a wide Mach number range. Based on the flow field and parameter fitting, two prediction methods for the curved axisymmetric inlet with lip translation scheme have been proposed. The method based on the flow field of the reference inlet is more efficient than the parameters-based prediction method, as it can accurately predict the lip translation distance and the corresponding flow coefficient over the entire working range of the inlet without additional numerical calculations. Moreover, the starting Mach number is accurately predicted by the fitting method based on the throat Mach number of the reference inlet, with a relative error of only 0.95% compared to the numerical simulation. The flow coefficient-based method is simple and accurate for predicting lip translation distances with a known starting Mach number, with a relative error of only 1.65% compared to numerical simulations. The prediction approaches can overcome the drawbacks of the standard iterative algorithms and significantly enhance computational accuracy and efficiency.

1. Introduction

With the continuous expansion of hypersonic flight speed and airspace, the combined-cycle engine has become a research hotspot, with examples including the Rocket Based Combined Cycle engine (RBCC) [1], the Air Turbo Rocket engine (ATR) [2,3], the Turbine Based Combined Cycle engine (TBCC) [4,5,6], the Pre-cooled Cycle engine [7,8,9,10], and so on. The inlet is one of the core components of the combined cycle engine, and its operating range needs to be further extended. In this case, a fixed geometry inlet can neither take into account the performance requirements of high Mach numbers and the starting requirements of low Mach numbers at the same time nor satisfy the engine flow requirements over a wide range of Mach numbers. Therefore, the adoption of variable geometry in the design of inlets is an inevitable trend [11,12,13,14].
The axisymmetric inlet has good engineering design advantages due to its simple structure, convenient manufacturing, and high utilization of windward areas. Since the 1950s, many studies have been carried out and the results have been successfully applied to hypersonic vehicles, such as the US Hypersonic Research Engine (HRE) [15], the Hyfly project [16], and the CIAM-NASA Scramjet Flight Program [17]. Xie et al. [18] systematically introduced the design method of the axisymmetric inlet and conducted wind tunnel experiments. Wang et al. [19] carried out a multi-objective optimization design of the axisymmetric inlet profile. Due to the narrow operating range of the engine at the early stage, the studies focused on the fixed geometry inlets, while for the later-developed engine with a wide operating range, especially combined cycle engines, the variable geometry inlet was required.
Most variable geometry inlets adopt the scheme of moving the center cone, such as the classic TBCC engine inlet of J58 [20], the RBCC engine inlet of the GTX plan [21], and the SABRE (Synergistic Air-Breathing Rocket Engine) single-cone multistage adjustable axisymmetric inlet [22]. Maru et al. [23] proposed a multi-row disk (MRD) inlet which can automatically change the distance between multi-row conical disks to ensure that the shock wave intersects at the lip under different inflow Mach numbers, but there are extra losses between the disks. Kojima et al. [24] used a movable cone to control the conical shock before the inlet in the air-turbo ramjet expander-cycle engine (ATREX) to keep the throat Mach number at about 1.3, thereby facilitating high compression performance. Cheng et al. [25] proposed a variable geometry inlet concept combined with local sub-flow circulation which can improve the flow capture capability of the inlet at off-design points and keep the total pressure recovery coefficient unchanged. Yuan et al. [26] proposed a three-dimensional aerodynamic design scheme of the inlet operating at Mai = 1.5–4.0 and control laws of the bleed slots are gained. Cai et al. [27] proposed an axisymmetric adjustable inlet that can meet both flow distribution and compression requirements of two channels in a wide working range for the precooled combined-cycle engine. The inlet is controlled by the coordinated translation of the central cone and the splitter plate. Nan et al. [28] performed geometric adjustment by moving the central cone to match the aerodynamic compression surface so as to achieve full flow capture.
Most of the above axisymmetric inlets are designed using the traditional straight-cone method, so the performance of the off-design points is significantly decreased. Recently, Zhang and co-workers [29] proposed a new design method of a curved compression inlet and carried out a series of experiments. The results show that the inlet performs well in the whole operating range. Under the condition of given Mach number distribution, Li et al. [30] used the method of characteristics to inversely design the axisymmetric inlet with good performance. In addition, the variable geometry inlet usually has a complex adjustment structure, which will induce sealing and cooling problems and add additional structural mass. A new variable geometry scheme with simple control and convenient adjustment, called the lip adjustment scheme, was developed to overcome the shortcomings. It can not only ensure the self-starting performance under a low Mach number but also can easily adjust the flow coefficient [30,31,32].
In this paper, the methods for fast and accurate prediction of the flow capture capability and starting performance of the inlet are studied for the curved compression axisymmetric inlet with lip translation over a wide range of Mach numbers. The structure of the flow field in the external pressure section of the traditional binary multi-wedge inlet is simple, and the inviscid calculation can be carried out through the shock relations, which yield lip-adjusted positions corresponding to the desired flow coefficients. The flow fields behind conical and curved shocks are too complex to be calculated by those methods, so the iterative trial calculation method is used [33]. In addition, there are barely published literature reports about the prediction of the starting performance of the curved compression axisymmetric inlet [34]. Therefore, the prediction methods based on the flow field and parameters fitting are proposed; these are verified by numerical simulation. It is shown that the prediction accuracy and computational efficiency are very high and can meet prediction requirements of the flow and starting performance of the variable geometry inlet over a wide range of Mach numbers.

2. Aerodynamic Design of the Reference Inlet

The reference inlet (named Inlet1), which is the reference configuration of the variable geometry axisymmetric inlet in the range of Mai from 3.0 to 7.0, was designed. The design Mach number is Mai = 6.0 and the external pressure section was inverse-designed with the arctangent Mach number distribution law. This distribution law can weaken the intensity of the leading edge shock and increase the isentropic compression ratio, thereby improving the overall performance of the inlet. Based on the parameter studies of the external pressure section and the internal contraction section, a set of optimized parameters was selected to complete the design of the reference inlet, as shown in Table 1. The structure and wave system of the reference inlet is shown in Figure 1, where a simple lip translation scheme is used and the green part can be moved. Compared with the central cone movement scheme, this scheme is simpler and needs a smaller driving force [30].

3. Numerical Calculation Method

3.1. Numerical Details

Numerical simulation of the steady and unsteady flow fields in an axisymmetric inlet with an angle of attack of 0 degree was performed with the Fluent solver. The effectiveness of this solver has been widely verified by calculating hypersonic flow, which can meet the steady and unsteady calculation accuracy of the hypersonic inlet [35,36,37]. The Reynolds averaged Navier-Stokes equations in two dimensions were solved using a finite volume spatial discretization method. During the computation, the inviscid flux scheme was Roe’s method, and the splitting format of the convective flux was the AUSM (Advection Upstream Splitting Method) format. The discrete scheme is the second-order upwind scheme, and the turbulence model is the SST (Shear Stress Transport) k - ω with some extensions for modeling high-speed, wall-bounded flows and separation regions [38,39]. The calculation process regards the gas as ideal and considers the variation of the specific heat capacity and thermal conductivity with temperature and solves the viscosity using Sutherland’s formula. In addition, the unsteady flow field was solved by a two-time-step implicit time-discrete method, where the physical time-step Δt is given according to the following equation [35]:
Δ t < δ λ
where δ and λ are the minimum mesh size and characteristic wave velocity of the unsteady computational domain, respectively. In the present study, the physical time step is set as 10−6 s with 50 inner iterations in each time step. In the calculation, the residual and throat pressure were monitored. The convergence criterion is that the residuals of each equation and the monitored parameters are reduced to 10−4 or they no longer change; additionally, the import and export flow remains stable. The inflow parameters are listed in Table 2, where Mai and Hi are the inflow Mach number and height; pi and Ti are the inflow static pressure and static temperature.

3.2. Turbulence Model Validation

The simulation accuracy of the SST k - ω model for the hypersonic inlet flow field was verified by the experimental data in the open-circuit wind tunnel at the Aachen Jet Propulsion Laboratory [40]. The experimental model is a 2D mixed compression inlet, as shown in Figure 2. In this original model geometry, two consecutive ramps deflect the incoming flow. However, due to the restrictions of the wind tunnel, the first ramp of the inlet was omitted in the experimental tests. The given angle of attack of the freestream, α = 10°, arises from the missing first ramp. The outline of the inlet and the definition of the main dimensions are given in Figure 2 and Table 3. The presented two-dimensional numerical simulations belong to a configuration with an isolator length of l = 79.3 mm at Mai = 2.41 without throttling. The total pressure and total temperature are 540 kPa and 305 K, respectively. The grids near the wall are refined to ensure y+ below 1.0, and the total grid is 141,000 structured cells.
Compared to the experimental schlieren image, the computational results give a good prediction of the flow structure of the hypersonic inlet, as shown in Figure 3. The shock pattern, the appearance and extension of the separation packages, and the development of the boundary layer are in good agreement with the experiment. In addition, the calculated static pressure distributions of ramp and cowl surface (see Figure 4 and Figure 5) agree well with the experimental results. Overall, the applied numerical simulation yields satisfying results, so the SST k - ω turbulence model can be used in the current study.

3.3. Unsteady Numerical Method Validation

In order to evaluate the ability of Fluent code to solve unsteady compressible flow, an experimental study of the shock-focusing reflection process was conducted by Izumi et al. [41]. The experimental device is shown in Figure 6. The section size of the shock tunnel is 44 mm × 44 mm. The expression of the shock reflector surface is X = CY2. The height of the parabolic reflector is 28 mm. In the unsteady numerical verification calculation, the conditions of C = 0.5 and the Mach number of the moving shock Mas = 1.5 are selected. The computational domain mesh and boundary conditions are shown in Figure 7 with about 60,000 cells.
Assuming the specific heat ratio γ is 1.4, the vacuum pressure and temperature of the shock tunnel are set to 101,000 Pa and 300.00 K, respectively. The velocity of the moving shock-front V1 is zero and the parameters after the shock are calculated by the shock relation (Equations (2)–(4)), as shown in Table 4.
p 2 p 1 = 1 + 2 γ γ + 1 ( M a s 2 1 )
T 2 T 1 = ( 2 γ γ + 1 M a s 2 γ 1 γ + 1 ) ( 1 + γ 1 2 M a s 2 ) γ + 1 2 M a s 2
V 2 V 1 a 1 = 2 γ + 1 ( M a s 1 M a s )
Numerical calculations were conducted by solving the two-dimensional unsteady Euler equations. The time step for unsteady numerical calculation is 10−8 s, with 50 iterations per time step. Figure 8 shows the comparison between the experiment and the computation at different times during the shock-focusing reflection process. The incident shock propagates from left to right and the reflected shock propagates from right to left. The time t is nondimensionalized with respect to γ1/2D/a1, where γ is a specific heat ratio, a1 is the speed of sound in air ahead of the shock, and the half height of the shock reflector D is 14 mm. The upper part comprises the experimental schlieren images at different times, and the lower part comprises the numerical schlieren images at corresponding times. At the dimensionless time t’ = 1.68, the incident shock disappears, and only the reflected shock remains. The numerical calculation results capture well the whole shock reflection process and the flow field structure is highly consistent, which verifies that the Fluent software is reliable for solving unsteady compressible flows.

3.4. Grid Analysis

Considering the existence of shock/boundary layer interaction in the hypersonic inlet, the mesh is refined near the wall and in regions that require focusing, such as near the lip and in the internal pressure section. The mesh near the wall is refined in the boundary layer to ensure that the wall y+ is less than 1.0. The mesh and boundary conditions for the steady and unsteady computations are shown in Figure 9. The inflow conditions of the inlet are set as the pressure far field. The solid wall boundary is set to be a no-slip adiabatic wall. The exit is set as a pressure outlet. The performance of a grid-independent analysis confirms that the grid resolution used in this paper is sufficient to capture the physically relevant features. In Figure 10, the static pressure distribution along the central cone surface is shown for three different grid-refinement levels: coarse (30,320 structured cells), medium (127,000 structured cells), and fine (434,400 structured cells), indicating the grid at medium and beyond is sufficient to obtain the accurate solution. Therefore, a medium grid was selected for further numerical analysis.

4. Results and Analysis of the Reference Inlet

4.1. Flow Features

The two-dimensional steady calculation was carried out for the axisymmetric reference inlet in the range of Mai = 3.0–7.0, and Figure 11a–e illustrates the flow fields of the reference inlet. At the design point of Mai = 6.0, as shown in Figure 11a, the curved leading edge shock just hits the lip position, and the flow coefficient is 1.00. After the leading edge shock, the Mach number decreases continuously in the radial direction. The lip-curved shock is weak, with almost no separation at the shoulder of the inlet. As the compression increases, the Mach number decreases with the streamline. The isentropic compression wave is located between the leading edge shock and the lip shock. As shown in Figure 11b, the leading edge shock is more curved and penetrates into the lip at Mai = 7.0. Some expansion waves are generated near the lip, and the flow is accelerated to Mach number 9.0. The curved lip shock becomes stronger, hitting the front of the throat and producing a small separation package. At Mai = 5.0 and 4.0, as shown in Figure 11c,d, the leading edge shock waves become straighter and the overflow windows become larger, the lip shock waves become weaker, and there is no separation near the shoulder. At Mai = 3.0, Figure 11e shows that the inlet is not starting, but the flow field is generally stable, and there is no periodic oscillation phenomenon. At this time, the separation package and the induced oblique shock appear near the lip. The upper part of the internal pressure section is mainly a supersonic zone and the actual flow area gradually decreases. There is a large, closed separation package near the cone which occupies more than half of the channel, and the exit Mach number is close to 1.0.
In order to perform a careful analysis of the unstarted flow field for an axisymmetric reference inlet with a low Mach number, the unsteady calculation was performed at Mai = 3.0. The flow field shown in Figure 12a was selected as the initial flow field at t = 0 ms. The calculation results are shown in Figure 12b–g. Since Mai = 3.0 is lower than the starting Mach number, the captured flow entering the inlet cannot all pass through the throat, which is choked at this time. Figure 12a shows that the supersonic flow cannot propagate disturbance forward and the captured flow accumulates continuously at the throat, which will increase the pressure continuously. The high-pressure perturbation propagates upstream through the subsonic region of the boundary layer, forcing the boundary layer to thicken and bend the streamlines. Since the Mach number here is low, the interaction between the shock boundary layer is not serious. There is no significant separation package, only a small subsonic region at the upper and lower walls.
As the pressure further increases due to accumulation, the pressure perturbation continues to propagate upstream with the boundary layer. The separation packages appear on both the cowl wall and the cone (Figure 12b), and they as well as the induced shock wave continue moving forward (Figure 12c). Due to the increasing Mach number of the wavefront, the induced shock becomes stronger and the interaction between the shock wave and the boundary layer is enhanced. Separation packages are getting larger and longer. The actual flow capacity was further weakened. In addition, the interaction between the induced shock waves further aggravates the complexity of the flow field, with the highest separation package at the incident point of the induced shock. On the one hand, the interaction between the shock waves and the separation packages makes the induced shock waves move forward and strengthen with the forward movement of the separation packages (Figure 12d) until they reach the lip (Figure 12e). On the other hand, the induced shock enhances the separation package, which generates a dynamical coupling property. It can also be seen from Figure 12e that the induced shock waves near the cowl wall move to the front of the lip and gradually become a normal shock. The induced shock on the cone intersects with the normal shock wave, forming a typical “λ” shock, whose slip layer makes the flow field more complex, and the flow near the throat approaches sound velocity. The induced shock and separation package significantly increase the total pressure loss, which further reduces the actual flow capacity of the throat. Therefore, the induced shock at the cone can only stop far from the lip. Figure 12f shows the farthest position of the separation package from the lip at t = 3.30 ms. Stronger separation-induced shock appears in front of the lip, further enhancing the overflow capability of the inlet. The disturbance inside the lip can overflow from the outside of the lip after the induced shock leaves the lip, so the pressure disturbance can only spread forward from the separation zone on the cone, resulting in a significant increase in the separation package, which occupies most of the intake of the internal pressure section. Then there is a flow-balancing process. The separation package on the cone starts to move backward and decrease, and the induced shock starts to become weaker. The unstarted flow field is generally stable and no longer changes after t = 5.05 ms (Figure 12g). There is a separation-induced shock in front of the lip of the stable unstarted flow field. The separation package at the incidence point of the lip-induced shock is the largest, and a local subsonic flow appears after the shock. The flow in the supersonic region behind the lip-induced shock accelerates to Mach 1.8, the channel in the supersonic region further back becomes smaller, and the flow continues to decelerate towards 1.0.
In summary, the unstarted process with the low Mach number is a complex, unsteady process. Over time, the separation package adapts to the increasing or decreasing back pressure by the simultaneous adjustment of shape, position, and induced separation shock, and then reaches a steady state, which is known as the dynamic self-sustaining stability of the separation package.

4.2. Performance Parameters

The performance parameters of the reference inlet are calculated using the mass-weighted method, as stated in Table 5. These parameters could be used to quantify the inlet performance, including the flow coefficient φ, total pressure recovery coefficient σ, pressure ratio p/pi, and so on. Subscripts: th is the throat, and e is the exit.
The definition of flow coefficient φ is defined as:
φ = m ˙ e m ˙ c
wherein m ˙ e indicates the mass flow at exit; m ˙ c indicates the mass flow captured by the inlet. The coefficient of σ is defined as:
σ = p t p ti
herein, p t signifies the mass-weighted total pressure at the section; p ti denotes the inflow total pressure. The specific flow density function q (Ma) is given below:
q ( M a ) = M a [ 2 γ + 1 ( 1 + γ 1 2 M a 2 ) ] γ + 1 2 ( γ 1 )
wherein γ is the specific heat ratio. The inlet shows good performance for Mai = 3.5–7.0, especially the high flow coefficients at low Mach numbers, with flow coefficients greater than 0.60 for Mai = 3.5. At Mai = 3.0, the inlet is unstarted and the flow coefficient and the total pressure recovery coefficient decrease significantly, with the flow coefficient only half of its value at Mai = 3.5.
Ref. [23] presents an optimized high-performance tri-cone inlet (named Inlet2), and the geometry and performance parameters of the curved compression reference inlet (Inlet1) proposed are compared with it. First, the external pressure section of the reference inlet is 5.20% shorter than that of Inlet2. As the lengths of the isolator are not equal, the throat parameters are chosen for comparison, as shown in Figure 13. At the design point with Mai = 6.0, the performance of the reference inlet is basically equivalent to that of the tri-cone inlet, but at the off-design point with Mai less than 5.2, the performance of the reference inlet, especially the flow coefficient, is higher than that of Inlet2. The flow coefficient and total pressure recovery coefficient of the reference inlet are 10.4% and 1.36% higher than those of Inlet2 at Mai = 4.0, respectively. Thus, the comprehensive performance of the reference inlet is better over a wide range of Mach numbers.

5. Prediction Method of Flow Coefficient

The inlet must satisfy the engine flow requirement over a wide range of Mach numbers, so two prediction methods have been proposed: one based on the flow field and the other on parameter fitting.

5.1. Prediction Method Based on Flow Field

There is a one-to-one mapping relationship between the lip translation distance and the flow coefficient in the reference inlet flow field at a specified Mach number. According to this characteristic, by knowing the flow fields in a wide Mach number range the lip translation distance (defined Lm) under the target flow coefficient (defined φt) can be calculated. Figure 11 shows the flow fields at the reference inlet, and an example is given to illustrate the specific prediction method. Compared with the flow coefficient of the reference inlet, the flow coefficient at the mode transition Mach number (Mai = 4.0) requires an increase of 10%, and the flow coefficient at other Mach numbers such as Mai = 5.0 requires an increase of 6%.
Step 1 Based on the flow coefficient requirement, the target flow coefficient corresponding to one of the inflow Mach numbers is firstly selected for the preliminary prediction of the lip translation distance. In this paper, the reference inlet flow field with modal conversion Mach number Mai = 4.0 (see Figure 14) was selected. AB is the curved leading edge shock, the flow before AB is the free flow, CD is the lip translation line, Rc is the capture radius, Rt is the radius of the streamline before AB, and EC is the lip translation distance. Equation (8) is derived from the definition of flow coefficient (Equation (5)), which shows the relationship between Rt and Rc. In this example, Rc = 200 mm, φ = φt = 0.769, Rt = 175.374 mm can be obtained according to Equation (8), and the corresponding streamline is made according to Rt. The intersection point of the streamline and CD is E, then EC is the lip translation distance corresponding to the target flow coefficient. In this example, Lm = EC = 11.28 mm.
R t = R c φ
Step 2 Besides the mode transition Mach number (Mai = 4.0), continue to calculate the flow coefficient of the inlet at other Mach numbers after translating the lip distance by the value of EC, then evaluate whether it meets the requirements. In this example, Mai = 5.0 is selected for illustration, as shown in Figure 15. In the flow field of the reference inlet at Mai = 5.0, the position of point E is determined from the lip translation line at EC = 11.28 mm. Then, a streamline is made through point E, resulting in a radius Rt = 195.35 mm. The corresponding flow coefficient φ = 0.953 can be calculated by Equation (8), which is higher than the target flow coefficient φt = 0.934 at Mai = 5.0, indicating that the lip translation distance EC = 11.28 mm meets the design requirements.
Step 3 follows Step 2 by traversing all the other Mach numbers, computing the actual flow coefficients corresponding to the lip translation distance EC, and evaluating whether they meet the requirements. If all the requirements are met, the obtained lip translation distance EC is the final value. If not, choose one of the Mach numbers and the corresponding target flow coefficient that does not satisfy the requirement, and repeat the above procedure until all flow coefficients meet the requirement.
This method can quickly obtain the lip translation distance that meets the entire flow coefficient requirement in a given Mach number range, avoiding massive numerical simulations and thus greatly saving compactional time and resources.

5.2. Prediction Method Based on Parameter Fitting

Step 1 The first two-dimensional viscous numerical calculations of the above variable geometry inlet were performed at the mode transition Mach number Mai = 4.0, and the specific performance parameters are given in Table 6, with a positive sign for Lm indicating forward motion and a negative sign for backward motion. It can be seen that as the lip moves forward, the total pressure recovery coefficient and the pressure ratio of the exit increase. This is because the increase of the internal contraction ratio improves the compression efficiency, which is exactly the opposite when the lip moves backward.
Step 2 Polynomial curve fitting is performed on the flow coefficient φ and the translation distance Lm in Table 6, and the fitting formula is obtained as Equation (9), the correlation coefficient of which is R2 = 1.000, as is shown in Figure 16. Given the target flow coefficient φt = φ = 0.769, the translation distance Lm = 11.79 mm can be obtained according to Equation (9). The moving distance of the lip increases from −20 mm to 20 mm, and the flow coefficient increases from 0.557 to 0.812, indicating that the actual capture flow of the inlet can be effectively adjusted by changing the lip position.
φ = 3.61 × 1 0 5 L m 2 + 6.36 × 1 0 3 L m + 0.699
Step 3 Except for Mai = 4.0, continue with numerical simulations to calculate the flow coefficients corresponding to other inflow Mach numbers with lip forward Lm = 11.79 mm to evaluate whether they can meet the requirements. If all the requirements are satisfied, then the Lm = 11.79 mm is the final value. If not, one of the inflow Mach numbers and the corresponding target flow coefficient that does not satisfy the requirement is selected and the above steps are repeated until the flow coefficients satisfy the requirements for all Mach conditions.

5.3. Numerical Verification

The accuracy of the above prediction methods is verified by numerical simulations. Two-dimensional numerical calculations were performed for an axisymmetric inlet with a forward lip of 11.28 mm. The performance parameters at Mai = 4.0 and 5.0 are shown in Table 7. The predicted values of the flow coefficients are exactly equal to the numerical values, which proves that the first prediction method is very accurate. The second method is computationally intensive, less accurate, and does not predict the flow coefficients at other Mach numbers. Using a desktop computer (4-core X64 processor, Intel i7-8700 CPU, 8 GB memory), the numerical simulation of a case in Table 6 takes 1890 s. As a result, the additional eight cases in the second method take a total of 15,120 s, while the time can be reduced to 180 s by using the first method.
Figure 16 shows the flow field of the inlet with Lm = 11.28 mm at Mai = 4.0 and 5.0. Compared with the reference inlet (Figure 11), the shape of the leading edge shock is unchanged. The forward translation of the lip reduces the overflow and further improves compression efficiency.
Figure 16. The flow fields of the inlet with the lip moving forward at Mai = 4.0 and 5.0; (a) Mai = 4.0; (b) Mai = 5.0.
Figure 16. The flow fields of the inlet with the lip moving forward at Mai = 4.0 and 5.0; (a) Mai = 4.0; (b) Mai = 5.0.
Aerospace 10 00506 g016

6. Prediction Method of Starting Mach Number

In this paper, the reference inlet is unstarted at Mai = 3.0. There are two important issues that need to be addressed: one is the prediction of the starting Mach number, and the other is the prediction of the critical translation distance of the lip if the inlet can be started at Mai = 3.0.

6.1. Prediction of Starting Mach Number of Reference Inlet

Step 1 The root reason for the reference inlet at low Mach numbers not starting is throat choking [34,42]; the flow state of the throat affects the starting performance of the inlet directly. Therefore, the throat Mach number is chosen to predict the starting Mach number of the reference inlet in this paper. The throat Mach numbers corresponding to the low Mach number of Mai = 3.5–4.5 are fitted by an exponential function (Equation (10)), and the correlation coefficient is R2 = 0.9999:
1 q ( M a th ) = 0.1774 e 0.585 M a i
where q(Math) is the specific flow density function of the throat.
Step 2 Assuming that the reference inlet is in a critical starting state, without considering the viscosity and the non-uniformity of the throat, the ideal throat Mach number Math = 1.0 is introduced into Equation (10) and the starting Mach number Mas = 2.95 is obtained. However, the numerical result shows that the inlet is unstarted at Mai = 3.0 (Figure 12), indicating that the predicted value is smaller, and 1/q(Math) = 1.0 needs to be corrected. The starting Mach number Mas = 3.15 is obtained by trial calculation of numerical simulation; the throat Mach number Math = 1.41 and the corresponding 1/q(Math) = 1.12, so 1.10 is chosen as the modified value of 1/q(Math). Substitute this modified value into Equation (10) and the starting Mach number Mas = 3.12. The relative error of predicted Mas with the numerical result is −0.95%, which meets the requirements of prediction accuracy.

6.2. Prediction of Lip Translation Distance at the Specified Starting Mach Number

Due to the low Mach number at which the reference inlet is unstarted, the starting of the inlet can be achieved by reducing the actual capture flow through a backward translation of the lip. At the specified starting Mach number Mas = 3.0, the key to accurately predicting the lip backward distance is to quickly determine the critical flow coefficient. Two prediction methods are proposed in this paper.

6.2.1. Prediction Method Based on Fitting of Flow Coefficient

The flow coefficients of low inflow Mach numbers (Mai = 4.5, 4.0, and 3.5) in Table 5 are selected for polynomial fitting (Equation (11)), and the correlation coefficient is R2 = 1.000:
φ = 0.054 M a i 2 + 0.583 M a i 0.769
The flow coefficient φ is obtained by fitting the inflow Mach numbers greater than the starting Mach number, so Equation (11) can be used to predict the critical flow coefficient corresponding to the starting Mach number. By taking Mai = 3.0 into Equation (11), φ = 0.494 can be obtained. According to the flow coefficient prediction method based on the flow field in Section 5.1, the lip translation distance Lm = −6.15 mm when φ = 0.494.

6.2.2. Prediction Method Based on Fitting of Total Pressure Recovery Coefficient

The total pressure recovery coefficients of the throat of Mai = 4.5, 4.0, and 3.5 in Table 3 are selected for polynomial fitting (Equation (12)), and the correlation coefficient is R2 = 1.000:
σ th = 0.01 M a i 2 + 0.049 M a i + 0.785
Therefore, the fitting value of σth is 0.842 at Mai = 3.0. According to the definition of the flow coefficient under zero-dimensional conditions, φc also can be calculated by:
φ c = σ th q ( M a th ) q ( M a i ) / R ct
where the values of parameters at Mai = 3.0 are given in Table 8. It can be seen from Section 6.1 that 1/q(Math) is not equal to 1.0, and the flow coefficient φc = 0.539 can be obtained by Equation (13).
According to the numerical result of the external pressure section of the reference inlet at Mai = 3.0, as shown in Figure 17, the flow coefficient φ is 0.526 calculated by Equation (5), thus the flow coefficient φc calculated by Equation (13) is unreasonable. Comparing the flow coefficient calculated by Equation (13) φc with φ in Table 5, it was found that φ/φc is less than 1.0. The distribution of throat parameters is not uniform, so the correction coefficient η is introduced to correct the flow coefficient φc by Equation (14).
φ = η φ c
wherein η corresponds to the known inflow Mach numbers shown in Table 9. According to the linear interpolation at Mai = 4.0 and 3.5, η = 0.911 is obtained at Mai = 3.0, so φ = 0.491 is calculated by Equation (14). According to the flow coefficient prediction method based on the flow field in Section 5.1, the lip translation distance Lm = −6.49 mm when φ = 0.491.

6.2.3. Numerical Simulation Verification

In Section 6.2.1 and Section 6.2.2, the critical flow coefficients φ at the starting Mach number Mas = 3.0 are obtained by two different prediction methods, which are 0.494 and 0.491, and the corresponding lip translation distances Lm are −6.15 mm and −6.49 mm, respectively. The numerical results are shown in Figure 18, where the inlet with lip translation distance Lm = −6.15 mm is started at Mai = 3.0, and its critical translation distance corresponding to Mas = 3.0 is −6.05 mm. The relative errors of the prediction methods based on the fitting of flow coefficients and the total pressure recovery coefficients are 1.65% and 7.27%, respectively, indicating the higher accuracy of the first method.

7. Conclusions

A curved compressible axisymmetric reference inlet with controllable Mach number distribution was designed, and then prediction methods were proposed for the flow coefficient and starting performance of the inlet lip translation scheme. The results were analyzed and verified by numerical simulation. The conclusions are as follows:
(1)
The overall performance of the reference inlet is good and the flow fields have almost no separation within the range of Mai = 3.5–7.0. The flow coefficient remains above 0.60 at Mai = 3.5. Compared to the optimized high-performance traditional tri-cone inlet, the reference inlet has better comprehensive performance over a wide range of Mach numbers, with 10.4% and 1.36% higher flow coefficient and total pressure recovery coefficient at Mai = 4.0, respectively. At Mai = 3.0, the inlet is unstarted, the flow field is basically stable after t ≥ 5.05 ms, and its separation package shows obvious self-sustaining stability.
(2)
The prediction method based on the flow field can accurately predict the flow coefficients of the inlet with lip translation in the whole working range only by using the flow field of the reference inlet, while the prediction method based on parameter fitting needs to calculate the performance of multiple lip positions, which is computationally expensive and has low accuracy.
(3)
The starting Mach number of the reference inlet can be accurately predicted based on the throat Mach number fitting method, and the relative error with the numerical simulation result is less than 1%. In the lip translation distance prediction for the specified starting Mach number, the prediction method based on the fitting of the flow coefficient is simple and accurate and the relative error with the numerical simulation result is 1.65%, while the prediction method based on the fitting of the throat total pressure recovery coefficient fitting is relatively complex with more undetermined parameters, and the relative error is 7.27%.
The present study shows that the prediction methods for flow coefficient and starting performance can overcome the shortcomings of the traditional iterative algorithms. These methods can be considerably more computationally efficient and can quickly satisfy the design requirements of variable geometry inlets.

Author Contributions

Y.L. and D.S. performed the design method, computational investigations, and data analysis and contributed to the paper draft. Z.W. contributed to the simulation methodologies and to the paper draft. K.Z. contributed to reviewing and editing the paper draft. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (No. 11702205), the Jiangxi Province Innovation Leading Talent Project, the Open Fund for the Key Laboratory of Aeroengine Thermal Environment and Thermal Structure (No. CEPE2022005), and the Natural Science Foundation of Chongqing, China (No. cstc2020jcyj-msxmX0823).

Data Availability Statement

The datasets generated or analyzed during this study are available from the corresponding author on reasonable request.

Conflicts of Interest

The authors declare there is no conflict of interest regarding the publication of this paper.

References

  1. Li, J.; Zhu, G.; Zhang, S.; Yan, O.; Li, W.; Qin, F. Two-stage fuel injection performance of RBCC under scramjet mode. Aerosp. Sci. Technol. 2020, 105, 106062. [Google Scholar] [CrossRef]
  2. Nan, X.; Wang, S.; Li, P. Investigation on status and prospect of air turbine rocket. J. Rocket. Propuls. 2008, 34, 31–35. [Google Scholar]
  3. Kanda, T.; Kudo, K. Conceptual study of a combined-cycle engine for an aerospace plane. J. Propuls. Power 2012, 23, 301–309. [Google Scholar] [CrossRef]
  4. Liu, J.; Yuan, H.; Hua, Z.; Chen, W.; Ge, N. Experimental and numerical investigation of smooth turbine-based combined-cycle inlet mode transition. Aerosp. Sci. Technol. 2017, 60, 124–130. [Google Scholar] [CrossRef]
  5. Yu, H.; Guo, Y.; Yuan, X.; Wang, J. Flight/propulsion integrated control of over-under TBCC ENGINE BASED on GA-LQR method. Aerospace 2022, 9, 621. [Google Scholar] [CrossRef]
  6. Snyder, L.; Escher, D. Turbine based combination cycle (TBCC) propulsion subsystem integration. In Proceedings of the 40th AIAA/ASME/SAE/ASEE Joint Propulsion Conference and Exhibit, Fort Lauderdale, FL, USA, 11–14 July 2004; p. 3649. [Google Scholar]
  7. Shi, H. Parametric research and aerodynamic characteristic of a two-stage transonic compressor for a turbine based combined cycle engine. Aerospace 2022, 9, 346. [Google Scholar] [CrossRef]
  8. Wang, Z.; Wang, Y.; Zhang, J.; Zhang, B. Overview of the key technologies of combined cycle engine precooling systems and the advanced applications of micro-channel heat transfer. Aerosp. Sci. Technol. 2014, 39, 31–39. [Google Scholar] [CrossRef]
  9. Zhang, M.; Nan, X.; Liu, D. Principles and realizing ways of combined power system for precooling air turbo rocket. J. Rocket. Propuls. 2016, 42, 6–12. [Google Scholar]
  10. Taguchi, H.; Harada, K.; Kobayashi, H.; Hongoh, M.; Masaki, D.; Nishida, S. Mach 4 Simulating experiment of pre-cooled turbojet engine using liquid hydrogen. Aerospace 2022, 9, 39. [Google Scholar] [CrossRef]
  11. Bouchez, M.; Falempin, F. Air breathing space launcher interest of a fully variable geometry propulsion system status. In Proceedings of the 36th AIAA/ASME/SAE/ASEE Joint Propulsion Conference and Exhibit, Las Vegas, NV, USA, 24–28 July 2000; p. 2376. [Google Scholar]
  12. Zhang, H.; Guo, R.; Li, B. Research status of TBCC inlet and its key technologies. Acta Aerodyn. Sin. 2010, 28, 613–621. [Google Scholar]
  13. Liu, J.; Chen, J.; Yuan, H. Reynolds number effect on aerodynamic and starting characteristics of a two-dimensional hypersonic inlet. Aerospace 2022, 9, 811. [Google Scholar] [CrossRef]
  14. Meng, P.; Yang, H.; Man, Y. Development of hypersonic inlet vehicle integrative design technology. Phys. Gases 2021, 34, 66–83. [Google Scholar]
  15. Andrews, E.; Mackley, E. ANALYSIS of Experimental Results of the Inlet for the NASA Hypersonic Research Engine Aerothermodynamics Integration Model; NASA/TM: Washington, DC, USA, 1976; p. 336.
  16. Foelsche, R.; Leylegian, J.; Betti, A.; Chue, R.; Marconi, F.; Beckel, S.; Tayll, J.; Charletta, R.; Bakos, R. Progress on the development of a free flight atmospheric scramjet test technique. In Proceedings of the AIAA/CIRA 13th International Space Planes and Hypersonics Systems and Technologies Conference, Capua, Italy, 16–20 May 2005; p. 3297. [Google Scholar]
  17. Voland, R.; Auslende, A.; Smart, M. CIAM/NASA Mach 6.5 scramjet flight and ground test. In Proceedings of the 9th International Space Planes and Hypersonic Systems and Technologies Conference, Norfolk, VA, USA, 1–5 November 1999; p. 4848. [Google Scholar]
  18. Xie, L.; Guo, R. Numerical simulation and experimental validation of flow in mixed-compression axisymmetric supersonic inlet with fixed-geometry. Acta Aeronaut. Et Astronaut. Sin. 2007, 28, 78–83. [Google Scholar]
  19. Wang, S.; Eriqitai, C.; Ding, W. Multi-objective optimization design of hypersonic axisymmetric inlet. J. Aerosp. Power 2020, 35, 1392–1401. [Google Scholar]
  20. Colville, J.; Lewis, M. An aerodynamic redesign of the SR-71 inlet with applications to turbine combined cycle engine. In Proceedings of the 40th AIAA/ASME/SAE/ASEE Joint Propulsion Conference and Exhibit, Fort Lauderdale, FL, USA, 11–14 July 2004; p. 3481. [Google Scholar]
  21. Thomas, S.; Palac, D.; Trefny, C.; Roche, J. Performance Evaluation of the NASA GTX RBCC Flow Path; NASA/TM: Washington, DC, USA, 2001; p. 210953.
  22. Roger, L.; Alan, B. The SKYLON project. In Proceedings of the 17th AIAA International Space Planes and Hypersonic Systems and Technologies Conference, San Francisco, CA, USA, 11–14 April 2011; p. 2244. [Google Scholar]
  23. Maru, Y.; Tanatsugu, N.; Sato, T.; Kobayashi, H.; Kojima, T.; Okai, K. Multi-row disk arrangement concept for spike of axisymmetric air inlet. In Proceedings of the 40th AIAA/ASME/SAE/ASEE Joint Propulsion Conference and Exhibit, Fort Lauderdale, FL, USA, 11–14 July 2004; p. 3407. [Google Scholar]
  24. Kojima, T.; Tanatsugu, N.; Sato, T.; Kanda, M.; Enomoto, Y. Development study on axisymmetric air inlet for ATREX engine. In Proceedings of the 10th AIAA/NAL-NASDA-ISAS International Space Planes and Hypersonic Systems and Technologies Conference, Kyoto, Japan, 24–27 April 2001; p. 1895. [Google Scholar]
  25. Cheng, D.; Zhang, Y.; Gao, W.; Xue, Y. Investigation of a variable axisymmetric inlet with local secondary flow recirculation. J. Propuls. Technol. 2019, 40, 2003–2011. [Google Scholar]
  26. Yuan, H.; Liu, F.; Wang, X.; Zhou, Z. Design and analysis of a supersonic axisymmetric inlet based on controllable bleed slots. Aerosp. Sci. Technol. 2020, 107, 106347. [Google Scholar] [CrossRef]
  27. Cai, Y.; Jin, Z.; Zhou, J.; Zhou, H. Design scheme of combined multiple thermodynamic cycle engine inlet. Acta Aeronaut. Et Astronaut. Sin. 2020, 41, 123745. [Google Scholar]
  28. Nan, X.; Zhang, H. Design of supersonic axisymmetric variable geometric inlet with zero additional drag. J. Aerosp. Power 2020, 35, 1376–1382. [Google Scholar]
  29. Zhang, K. Review on curved surface compression technology of hypersonic inlet. J. Propuls. Technol. 2018, 39, 2227–2235. [Google Scholar]
  30. Li, Y.; Li, G.; Liu, X.; Ma, Y. Design of a curved axisymmetric variable geometry inlet with Mach number from 1.5 to 4.5. J. Rocket. Propuls. 2018, 44, 30–35. [Google Scholar]
  31. Jin, Z.; Zhang, K. Concept of a varied geometry scramjet inlet with rotatable cowl. J. Propuls. Technol. 2008, 29, 43–48. [Google Scholar]
  32. Sun, B.; Shao, Q.; Ni, K.; Zhou, C.; Zhang, K. Parallel plans for TBCC variable geometry inlet. J. Aerosp. Power 2015, 32, 2989–2996. [Google Scholar]
  33. Wang, Y.; Yuan, H.; Guo, R. Research on mass flow controlling technology of TBCC axisymmetric inlet. J. Propuls. Technol. 2013, 34, 1606–1615. [Google Scholar]
  34. Molder, S.; Timofeev, E.; Tahir, R. Flow starting in high compression hypersonic air inlets by mass spillage. In Proceedings of the 40th AIAA/ASME/SAE/ASEE Joint Propulsion Conference and Exhibit, Fort Lauderdale, FL, USA, 11–14 July 2004; p. 4130. [Google Scholar]
  35. Wang, W.; Guo, R. Unsteady characteristics of hypersonic inlet below self-starting Mach number. J. Aerosp. Power 2012, 27, 2733–2741. [Google Scholar]
  36. Huang, H.; Tan, H.; Sun, S.; Ling, Y. Evolution of supersonic corner vortex in a hypersonic inlet/isolator model. Phys. Fluids 2016, 28, 126101. [Google Scholar] [CrossRef]
  37. Gao, W.; Chen, J.; Liu, C.; Li, Z.; Yang, J.; Zeng, Y. Effects of vortex generators on unsteady unstarted flows of an axisymmetric inlet with nose bluntness. Aerosp. Sci. Technol. 2020, 104, 106021. [Google Scholar] [CrossRef]
  38. Menter, F. Two-equation eddy-viscosity turbulence models for engineering applications. AIAA J. 1994, 32, 1598–1605. [Google Scholar] [CrossRef]
  39. Huang, G.; Coakley, T. Turbulence modeling for complex hypersonic flows. In Proceedings of the 31st Aerospace Sciences Meeting, Reno, NV, USA, 11–14 January 1993; p. 200. [Google Scholar]
  40. Reinartz, B.; Herrmann, C.; Ballmann, J. Aerodynamic performance analysis of a hypersonic inlet isolator using computation and experiment. J. Propuls. Power 2003, 19, 868–875. [Google Scholar] [CrossRef]
  41. Izumi, K.; Aso, S.; Nishida, M. Experimental and computational studies focusing processes of shock waves reflected from parabolic reflectors. Shock. Waves 1994, 3, 213–222. [Google Scholar] [CrossRef]
  42. Van Wie, D.; Kwok, F.; Walsh, R. Starting characteristics of supersonic inlets. In Proceedings of the 32nd Joint Propulsion Conference and Exhibit, Orlando, FL, USA, 1–3 July 1996; p. 2914. [Google Scholar]
Figure 1. Structure and wave system of the reference inlet (Inlet1).
Figure 1. Structure and wave system of the reference inlet (Inlet1).
Aerospace 10 00506 g001
Figure 2. Main dimension diagram of the inlet [40].
Figure 2. Main dimension diagram of the inlet [40].
Aerospace 10 00506 g002
Figure 3. Comparison of the Mach number contour lines of the computation (top) and the schlieren image (bottom) without throttling.
Figure 3. Comparison of the Mach number contour lines of the computation (top) and the schlieren image (bottom) without throttling.
Aerospace 10 00506 g003
Figure 4. Comparison of pressure distribution for the ramp without throttling.
Figure 4. Comparison of pressure distribution for the ramp without throttling.
Aerospace 10 00506 g004
Figure 5. Comparison of pressure distribution for the cowl surface without throttling.
Figure 5. Comparison of pressure distribution for the cowl surface without throttling.
Aerospace 10 00506 g005
Figure 6. Shock focusing from a parabolic reflector [41].
Figure 6. Shock focusing from a parabolic reflector [41].
Aerospace 10 00506 g006
Figure 7. Computational domain meshes and boundary conditions.
Figure 7. Computational domain meshes and boundary conditions.
Aerospace 10 00506 g007
Figure 8. Comparison of experiment (upper half) with computation (lower half).
Figure 8. Comparison of experiment (upper half) with computation (lower half).
Aerospace 10 00506 g008
Figure 9. Computational domain and boundary conditions of the axisymmetric inlet.
Figure 9. Computational domain and boundary conditions of the axisymmetric inlet.
Aerospace 10 00506 g009
Figure 10. Pressure distribution of the central cone surface for refined grids.
Figure 10. Pressure distribution of the central cone surface for refined grids.
Aerospace 10 00506 g010
Figure 11. Flow fields of the reference inlet: (a) Mai = 6.0; (b) Mai = 7.0; (c) Mai = 5.0; (d) Mai = 4.0; (e) Mai = 3.0.
Figure 11. Flow fields of the reference inlet: (a) Mai = 6.0; (b) Mai = 7.0; (c) Mai = 5.0; (d) Mai = 4.0; (e) Mai = 3.0.
Aerospace 10 00506 g011
Figure 12. Dynamic process of the unstarted flow field at Mai = 3.0: (a) t = 0 ms; (b) t = 0.45 ms; (c) t = 1.05 ms; (d) t = 1.80 ms; (e) t = 2.30 ms; (f) t = 3.30 ms; (g) t = 5.05 ms.
Figure 12. Dynamic process of the unstarted flow field at Mai = 3.0: (a) t = 0 ms; (b) t = 0.45 ms; (c) t = 1.05 ms; (d) t = 1.80 ms; (e) t = 2.30 ms; (f) t = 3.30 ms; (g) t = 5.05 ms.
Aerospace 10 00506 g012aAerospace 10 00506 g012b
Figure 13. Throat parameters comparison of Inlet1 and Inlet2: (a) Flow coefficient and total pressure recovery coefficient; (b) Pressure ratio and throat Mach number.
Figure 13. Throat parameters comparison of Inlet1 and Inlet2: (a) Flow coefficient and total pressure recovery coefficient; (b) Pressure ratio and throat Mach number.
Aerospace 10 00506 g013
Figure 14. Prediction diagram of lip translation distance at Mai = 4.0.
Figure 14. Prediction diagram of lip translation distance at Mai = 4.0.
Aerospace 10 00506 g014
Figure 15. Prediction diagram of flow coefficient at Mai = 5.0.
Figure 15. Prediction diagram of flow coefficient at Mai = 5.0.
Aerospace 10 00506 g015
Figure 17. Flow field of external pressure section of the reference inlet at Mai = 3.0.
Figure 17. Flow field of external pressure section of the reference inlet at Mai = 3.0.
Aerospace 10 00506 g017
Figure 18. The starting flow field with different lip translation distances at Mai = 3.0: (a) Lm = −6.15 mm; (b) Lm = −6.05 mm.
Figure 18. The starting flow field with different lip translation distances at Mai = 3.0: (a) Lm = −6.15 mm; (b) Lm = −6.05 mm.
Aerospace 10 00506 g018
Table 1. Design parameters of the reference inlet.
Table 1. Design parameters of the reference inlet.
ParametersValueUnit
Intake capture radius; Rc200.0mm
Semi-cone angle of the leading edge; δc7.6°
Total compression angle; δt21.0°
Inner angle of the lip; θi12.0°
Total length of the inlet; Lt978.0mm
Length of the external pressure section; Lc730.0mm
Height of the throat; Hth16.0mm
Length of the isolator; Liso112.0mm
Total contraction ratio; Rct6.01-
Internal contraction ratio; Rci1.71-
Table 2. Different inflow parameters.
Table 2. Different inflow parameters.
MaiHi (km)pi (Pa)Ti (K)
6.025.02549.22221.55
7.027.01879.98223.54
5.0/4.523.03466.87219.57
4.0/3.520.08849.73216.65
Table 3. Main dimension parameters of the inlet [40].
Table 3. Main dimension parameters of the inlet [40].
ParametersValueUnit
Inner width52.0mm
Overall length400.0mm
Isolator height15.0mm
Second ramp angle δ221.5°
Cowl angle δ39.5°
Divergence angle δ45.0°
Table 4. Boundary conditions of the numerical calculation.
Table 4. Boundary conditions of the numerical calculation.
Boundary TypeMaipi (Pa)Ti (K)
Pressure far field0.6044248,290396.06
Pressure outlet-101,000300.00
Table 5. Overall performance parameters of the reference inlet.
Table 5. Overall performance parameters of the reference inlet.
Maiφσthpth/piMath1/q(Math)σepe/piMae
7.01.0000.43125.43.426.3000.35927.53.21
6.01.0000.64122.23.124.7470.53323.12.93
5.00.8810.76317.12.683.1230.65018.72.50
4.50.7610.80314.42.432.4710.68716.22.25
4.00.6990.82113.22.101.8370.71513.81.97
3.50.6100.83411.41.741.3760.75012.11.64
3.00.3080.46411.00.821.0300.4387.631.08
Table 6. Performance parameters of the inlet with lip moving forward at Mai = 4.0.
Table 6. Performance parameters of the inlet with lip moving forward at Mai = 4.0.
Lm (mm)φσthpth/piMathσepe/piMae
−200.5570.7939.682.270.6789.862.13
−150.5960.80410.462.230.69310.772.09
−100.6320.80411.452.180.69811.762.04
−50.6660.81612.322.150.71012.742.01
00.6990.82113.172.100.71513.791.97
50.7300.81714.322.050.71314.991.92
100.7590.81915.292.010.71716.031.88
150.7860.82016.181.980.72016.961.85
200.8120.82017.001.940.72417.861.82
Table 7. Overall performance parameters of the inlet with lip moving forward.
Table 7. Overall performance parameters of the inlet with lip moving forward.
Maiφσthpth/piMathσepe/piMae
4.00.7690.82015.52.000.71816.21.87
5.00.9530.76719.22.580.65220.92.43
Table 8. The parameter values of the Equation (13) at Mai = 3.0.
Table 8. The parameter values of the Equation (13) at Mai = 3.0.
ParametersValue
Mai3.0
q(Mai)0.236
q(Math)1.10
σth0.842
Rct6.01
Table 9. The correction coefficient η corresponds to the inflow Mach numbers.
Table 9. The correction coefficient η corresponds to the inflow Mach numbers.
Maiη
3.50.891
4.00.877
4.50.850
5.00.866
6.00.837
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Li, Y.; Sun, D.; Wu, Z.; Zhang, K. Flow Coefficient and Starting Performance Prediction of Variable Geometry Curved Axisymmetric Inlet. Aerospace 2023, 10, 506. https://doi.org/10.3390/aerospace10060506

AMA Style

Li Y, Sun D, Wu Z, Zhang K. Flow Coefficient and Starting Performance Prediction of Variable Geometry Curved Axisymmetric Inlet. Aerospace. 2023; 10(6):506. https://doi.org/10.3390/aerospace10060506

Chicago/Turabian Style

Li, Yongzhou, Di Sun, Zejun Wu, and Kunyuan Zhang. 2023. "Flow Coefficient and Starting Performance Prediction of Variable Geometry Curved Axisymmetric Inlet" Aerospace 10, no. 6: 506. https://doi.org/10.3390/aerospace10060506

APA Style

Li, Y., Sun, D., Wu, Z., & Zhang, K. (2023). Flow Coefficient and Starting Performance Prediction of Variable Geometry Curved Axisymmetric Inlet. Aerospace, 10(6), 506. https://doi.org/10.3390/aerospace10060506

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop