Abstract

The reverse-flow combustor is extensively applied in small engines due to its compact structure. The cold-state flow characteristics, the ignition process, and flame chemiluminescence characteristics of the reserve-flow combustor were investigated experimentally. Using the particle imaging velocimeter, the cold-state time-averaged flow fields at three different total pressure loss coefficients are examined. The ignition processes of the combustor, as well as flame characteristics during stable combustion at different equivalent ratios , are analyzed based on flame chemiluminescence imaging when is 3%. As increases from 1% to 5%, the flow field of the combustor remains almost unchanged while the flow velocity increases. For =3%, the ignition time drops from 26 ms to 16 ms with an increase of . In addition, the ignition performance is enhanced. Proportions of the T2 phase (quasi-stable phase) are 73%, 68%, and 66%, respectively, suggesting that the quasi-stable phase is the key to successful ignition. As increases, the flame becomes thinner and shorter. Furthermore, standard deviations of continuous image areas of and drop from 0.0143 and 0.0132 to 0.0115 and 0.0109, respectively, which indicates the enhancement of combustion stability. According to the brightness distributions of and along the axial direction of the combustor at different equivalent ratios, the production of is significantly affected by the temperature. It can be adopted as an important sign of heat release during combustion.

1. Introduction

The combustor of an aero-engine enhances the airflow temperature with the combustion; it should have characteristics such as reliable ignition, stable combustion, and favorable outlet temperature [1]. The flow field of the combustor directly affects the fuel atomization and mixing between fuel and the airflow, which in turn determines ignition and combustion performance. So, it has been the focus of research.

A spray device composed of a matching fuel nozzle and a multi-stage swirl device, usually called a swirl cup, is widely used in aero-engines. A lot of research has been done on the single-head swirl cup model. Mongia et al. [2] studied the influence of pressure and temperature on the time-averaged flow field, as well as the influence of the size of the swirl cup and the structure of the venturi on the flow field on the CFM-56 aero-engine. Li et al. [3, 4] carried out numerical simulation of the influence of inlet and outlet boundary conditions on the flow characteristics of a three annular cyclone (TARS) and carried out a study on the flow characteristics of the central recirculation zone and the shear layer between the inner and outer flows. Vashahi et al. [5] used the PIV technique to study multi-swirl flow in a rectangular model, comparing the distribution of axial and radial velocities to determine the confinement effects of a rectangular combustion chamber. Most of these studies are limited to straight-flow combustor. However, the airflow in the upper circuit and airflow in the flame tube of the reverse-flow combustor flow in opposite directions. These features make the combustion chamber structure more compact and the length of the engine shortened. It is mainly used in turboshaft engines. The difference in the structure of the reverse-flow combustor and the straight-flow combustor will inevitably lead to differences in the characteristics of flow field and combustion.

Some studies have been carried out on the reverse-flow combustor, such as Reddy et al. experimentally investigated cold-state flow characteristics in the reverse-flow combustor with the pitot tube probe. Moreover, the authors compared the experimental data with numerical calculation results. It was concluded that an increase in the Reynolds number at the inlet can lead to an increase in axial velocity and divergence of radial velocity in the combustor [68]. Bharani et al. employed a direction-sensitive pressure probe and standard hot-wire probe to measure the flow field in the reverse-flow combustor. In addition, the authors clarified the flow characteristics of the inlet airflow with different swirling intensities in the outer circular channel of the reserve-flow combustor. Lastly, they investigated the effects on jets from different jet holes [9, 10]. Currently, on account of compact structures, researchers mainly adopt probes and numerical simulations to examine the flow field characteristics and combustion performance in the reverse-flow combustor. However, the mainstream 3D strong rotational flow field test method employs particle imaging velocimeter (PIV) to examine the flow field structure. This method can measure 2D planar flow field based on true condition in the combustor, thereby obtaining the entire flow field structure in the combustor [1114]. Accordingly, numerical simulation results can be validated by the experimental data. This study reasonably arranged visual optical windows in the reserve-flow combustor model that can satisfy optical path requirements of some measurement instruments such as PIV. The flow field of the reverse-flow combustor is investigated with PIV. For the ignition performance of the combustion chamber, a reasonable flow field is the prerequisite for successful ignition. Favorable ignition performance is also an important index for measuring combustor performance. Lefebvre et al. conducted a significant number of investigations on the ignition performance of a combustor [1517]. Based on the model of the premixed combustion chamber, the authors investigated both ignition and flame propagation mechanisms and extracted the characteristics at different development phases of ignition [18]. Philip et al. validated the test results via large-eddy simulation (LES) calculation and further analyzed different ignition stages in detail [19]. Prieur et al. focused on the effects of equivalent ratio and flow velocity on the ignition time. The authors found that fuel distribution plays an important role in both the propagation and stabilization of flames [20]. Chemiluminescence combustion refers to the luminescence at a particular wavelength when free radicals excited in combustion chemical reaction transit from high- to low-energy levels; It is easy to use under severe environments without high-power laser and complex light paths [21, 22]. In hydrocarbon fuel, leading chemiluminescence behaviors are obtained by and at 430 nm and 307 nm, respectively. The main reaction mechanisms are described as follows:

and show different characteristics in the chemical reaction. These elements are extensively applied when determining chemical zones, identification of heat release, and evaluation of combustion equivalent ratio [2325]. Higgins et al. investigated the relations of chemiluminescence of with pressure and equivalent ratio, which expanded the application range of chemiluminescence [26]. Hardalupas et al. measured chemiluminescence imaging and examined the characteristics of premixed flames of natural gas at different equivalent ratios. The authors concluded that and are important markers of heat release rate [27]. Garcia et al. indicated application prospects of chemiluminescence when evaluating the distribution of the equivalent ratio of the combustor [28]. Orain et al. focused on the effects of equivalent ratios and different fuel compositions on the heat release rate. Moreover, the authors measured , , and via chemiluminescence imaging [29]. According to their results, different fuel compositions differ in the optimal marker of heat release rate.

The previous studies exhibit important guidance for exploring flow, ignition, and combustion characteristics. In this study, based on the reserve-flow combustor model, the inner flow field structure was measured with PIV, and the ignition process, as well as flame characteristics, was examined using flame chemiluminescence imaging. The present study can deepen the understanding of flow, ignition, and combustion process in a reverse-flow combustor.

2. Combustor Model and Experiment System

2.1. Reverse-Flow Combustor Model

The fan-shaped reserve-flow combustor is shown in Figure 1(a). Figure 1(b) shows a typical flow structure diagram of reverse-flow combustor. More details of the model combustion chamber can be found in [30], but the axial position and layout of the upper and lower jet holes are different from the structure used in [30]. Table 1 shows the parameters of jet holes.

2.2. Experiment Parameters

The experiment system is shown in Figure 2; it mainly consists of the reverse-flow combustor, the air supply system, the fuel system, and the measurement system. Further, the measurement system includes PIV measurement system and spontaneous radiation measurement system. Sampling rakes were arranged to measure inlet total pressure () and outlet total pressure () with a maximum error of 0.2%, and inlet total temperature () with a maximum error of 0.4%. The inlet mass flow rate () was measured by the LUGB/E vortex shedding flowmeter with a maximum error of 1%. The fuel pump was used to support the liquid fuel of RP-3 kerosene. The fuel mass flow rate () is measured through the pressure drop of the atomizer. The relationship between the flow rate and the pressure drop is fitted as:

The ignition electronic nozzle was arranged on the sidewall of the single-head combustor, which was 2 mm away from the sidewall. The PIV system was made up of a Bobcat B2041 digital camera. An Nd:YAG double-pulse laser with a laser wavelength of 532 nm, maximum single-pulse energy of 200 mJ, and a maximum operating frequency of 15 Hz was used. A MicroPulse 725 synchronous controller was used to control the camera, the laser, and the image capturing board [31]. The essence of PIV measurement is to characterize the flow structure by acquiring the flow characteristics of the tracer particles. Stokes number () is introduced as the criterion for selecting tracer particles [32]. is the ratio of particle relaxation time to fluid characteristic time, which represents the ratio of particle inertia and diffusion. The calculation formula is: where is the reaction time of the tracer particles, and is the characteristic time of the flow. When , the diffusion of particles is more dominant than the inertial force, and the PIV test results are more credible. So, the tracer particles were introduced at a position 1 m upstream of the intake section, and MgO particles with a particle size of 5 μm were selected as tracer particles in this study to ensure .

The total pressure loss directly affects the performance of the engine. The total pressure loss coefficient is one of the important parameters to characterize the state of the combustor. It can be calculated by:

The parameters of the experiment are shown in Table 2.

Due to the fact that the laser plane and the CCD camera viewing plane are perpendicular to each other in the PIV experiment, the actual measurement area of flow field at the center section P0 of the combustor is shown in Figure 3(a). Here, , , and denote the position lines at different axial distances from the outlet of the swirler ( =10 mm,  = 20 mm, and  = 40 mm). The origin of the coordinate system is located at the lower left corner of the measurement area, and the positive direction of axis is the same as the direction of the inlet velocity of the combustor, and opposite to the direction of the main flow inside the reverse-flow combustor.

The spontaneous radiation system adopts the combination of Phantom V2012 high-speed camera and Invisible version UVi image intensifier, equipped with Nikon 50 mm fixed focus lens, to capture high frame rate images of ignition process and image during stable combustion process. The filter with a wavelength of 430 nm and a bandwidth of 5 nm as well as an filter with a wavelength of 307 nm and a bandwidth of 5 nm was arranged in front of the camera. Since are imaged via the filters at particular wavelengths with fairly low light transmittance outside the bandwidth range, the effects of some particles, such as soot, on the measurement of can be neglected. In view of the structure of the combustor model and the framing characteristics of the camera, the actual measurement area of the flame at the central section P0 is shown in Figure 3(b). The sampling frequency was set to 2 000 Hz. 8-bit grayscale images were obtained. The pixels were calibrated before image acquisition. Finally, images with a resolution of were obtained [33].

2.3. Experimental Error Analysis

Errors generated in the test process can be calculated according to the following formula [34, 35]: where is the uncertainty of the derived physical quantity; is the uncertainty corresponding to each measurement parameter.

Taking the uncertainty of the equivalence ratio as an example, the calculation formula of is:

Here, is the stoichiometric ratio of RP-3 kerosene.  =0.067. According to formula (6) and (7), the uncertainty of is calculated as: where and are the uncertainties of and , respectively.

Table 3 shows the relevant values of uncertainties in experiments.

3. Results and Discussions

3.1. The Flow Field

The experiment was conducted at inlet static pressure 101 kPa and the temperature 300 K. The results and a fitting curve are shown in Figure 4. The pressure loss increases with the reference velocity () increasing. The largest pressure loss under experiment operations which is fixed at the maximum combustor inlet velocity (16.7 m/s) is 5%. The fitting curve is defined as: where the -Square is 0.999, which means the fitting curve is in accord with the experiment results ideally. So the higher the reference velocity is, the more drastic the total pressure loss changes.

The cold-state time-averaged flow fields at section P0 of the reverse-flow combustor with different values were measured by PIV. Figure 5 shows the cold-state time-averaged flow fields at section P0 when =1%, =3%, and =5%. The flow fields at section P0 under different values show similar structures with almost identical vortex positions, jet penetration depths, and flow lines. The vortex center positions of the vortex in the primary recirculation zone (primary vortex) and the vortex formed after the collision of the primary hole jets (secondary vortex) are recorded as R and S, respectively. The positions of R and S under different values are shown in Table 4. It means has little effect on the flow field structure.

Figure 6 displays the distribution of and at different axial distances ( =10, 20, and 40 mm). The values of and can be normalized by ; is the when is 5%. The effect of imposes consistent effects on flow structure at any axial positions, suggesting the irrelevance of the flow field structure to the inlet condition of the combustor. However, the value at the inlet can lead to an increase in the flow velocity. At  =10 mm, radial height in the central recirculation zone at the axial distance can be reflected by the axial velocity component at  ≈6-20 mm. Moreover, when  <20 mm, both velocity components and show low-velocity and high-velocity flow alternation along the radial direction. This can increase both and gradients and strengthen the flow mixing action. When  >20 mm, and change gently along the radial direction. At  = 20 mm, the radial height of the central recirculation zone increases, and it is accompanied by a large-gradient flow of and . When  >15 mm, the gradient of v further increases. However, basic flow characteristics remain almost unchanged. At  = 40 mm and  ≈4-36 mm, i.e., the primary zone of the combustor, flow moves smoothly to the combustor downstream, while the flow velocity near both upper and lower wall sizes varies significantly. This is mainly because a large adverse pressure gradient exists near the downstream wall of upper and lower primary holes, thereby producing local vortexes.

3.2. Ignition

It can be found that has a slight effect on the flow field of the reverse-flow combustor through the research in the previous section. The spontaneous radiation experiments were conducted at =3%.

The ignition process of the combustor was recorded. The equivalence ratio of the combustor was adjusted by adjusting the valves in the fuel supply pipeline to supply fuel with different mass flow rates. In order to highlight the details of the image and improve the anti-interference capability of the image, the grayscale images of flame obtained by the high-speed camera were converted into pseudo-color images through Matlab software. The typical flame development diagrams during the ignition process of the reverse-flow combustor are shown in Figure 7.

It can be assumed that the flame’s illumination intensity during the ignition process can reflect the instantaneous heat release rate of the combustion. Furthermore, by taking integral of the illumination intensities of various pixels in the gray-level images during the ignition process, variation of the illumination intensity per unit pixel with time can be obtained, as shown in Figure 8.

The ignition process of the reverse-flow combustor can be divided into three phases [36]. In the first turbulent dissipation phase, the high-energy spark plug discharges, evaporating and igniting the mixed gas. Flame intensity drops gradually to form the initial flame nucleus. At the same time, the followability of small-sized droplet particles is better, and better mixing is formed in PRZ, so the initial flame nucleus near the PRZ. Next, the second phase is the quasi-stable phase. A low-velocity recirculation zone is conducive for the formation and stabilization of the initial flame. At this phase, fresh mixed air entertained by the recirculating vortex absorbs energy. Simultaneously, the initial flame nucleus gradually develops along with the recirculating vortex, thereby leading to steady heat release. During the entire process, heat release is first slightly lower than heat loss. Then, with the development of the initial fire nucleus, heat release slightly exceeds the heat loss, thereby resulting in gradual energy accumulation. Finally, flame develops to stable combustion state, and the combustor ignites successfully at the third combustion-enhanced phase. The flame expansion first occurs in PRZ. The chemical reaction is enhanced with the continuous flame expansion, and the temperature of the unburned air-fuel mixture near flame and the fuel evaporation rate are greatly increased. At the same time, the rapid increase of unburned air-fuel mixture temperature and fuel evaporation rate will further promote the enhancement of chemical reaction. Flame expands rapidly under the action of this pair of positive feedbacks, and the flame develops from the upper PRZ to the lower PRZ. Finally, the combustor is successfully ignited. The combustion chemical reaction continues to be transmitted downstream and ignites the entire combustion chamber rapidly to achieve stable combustion with the action of turbulent dissipation; the flame intensity increases to a peak value and tends to be stable. The total ignition can be calculated as , where , , and denote the durations of three phases.

Table 5 lists the distributions of , , and during the ignition process at different equivalent ratios when =3%. The combustor can be successfully ignited at three different equivalent ratios. Moreover, the total ignition delay time decreases with the enhancement of equivalent ratio, accompanied by the enhancement of ignition performance. To be specific, remains almost unchanged since the first phase is mainly affected by the energy of the high-energy igniter itself. As the equivalent ratio increases, drops by 31.6% and 44.7%, while is reduced by 20% and 30%. The second and third phases are closely correlated with fuel/gas spatial distribution. Accordingly, the increasing equivalent ratio increases the fuel mist concentration in the combustor, which can significantly change the ignition times in the latter two phases. Under three operating conditions, the second phase occupies a dominant role in the entire ignition process. Proportions of in three phases are 73%, 68%, and 66%, while the duration of the quasi-stable phase decreases slightly with an increase in the equivalent ratio. Therefore, the second phase imposes a significant effect on successful combustor ignition. By combining the cool-state time-averaged flow field in Figure 4, the flame mainly stays in the central recirculation zone. Furthermore, flow velocity in the recirculation zone is low, and fresh air can be constantly supplemented, thereby leading to constant energy accumulation. As energy is accumulated to a certain degree, combustion can be gradually enhanced. Afterward, flame develops steadily in the central recirculation zone, becomes more intensive, and moves downstream towards the combustor. Finally, the combustor can be successfully ignited. Therefore, flow structure significantly affects the ignition performance of the combustor. It can be concluded that favorable flow structure design is of great significance.

3.3. Chemiluminescence

is tightly related to the combustion reaction and can be used to obtain heat release characteristics. Distributions of intermediate products and during stable combustion in the combustor at different equivalent ratios were measured via spontaneous radiation imaging. After cropping, filtering, and normalization on the original flame images via spontaneous radiation imaging, the continuous image series taken under each operating condition were averaged because of strong turbulent flame pulses. The averaged image represents a stable combustion image under the corresponding operating conditions. Finally, the image boundary was extracted by taking pseudo-color processing and OSTU algorithm [31], as shown in Figure 9.

During the combustion of hydrocarbon fuel, a tight correlation between and combustion reaction can be confirmed. Therefore, the distribution of can be used for reflecting the flame structure. Figure 10 shows the distributions of at three different equivalent ratios as well as the comparison of boundary contours. It can be observed that flames at different equivalent ratios show almost unchanged contour features. Based on the contour, the flame image can be divided into the flame development zone (Zone 1), complete combustion zone (Zone 2), and combustion wake zone (Zone 3). In Zone 1, the flame front develops downstream towards the combustor in a V-shaped pattern with an expansion angle of approximately 80°. The flame can be radially extended to the height of the flame tube, exhibiting typical non-premixed diffusion combustion characteristics. This suggests that the change of equivalent ratio imposes a minor effect on the flame’s initial expansion angle. The expansion angle mainly depends on the spray cone angle of the centrifugal nozzle. In Zone 2, the flame undergoes complete combustion and spreads to the entire flame tube. High-brightness zones near the upper and lower wall surface can be observed. However, the upper and lower high-brightness zones show different contour characteristics that correspond to the edges of the upper and the lower recirculation zones in the flow field. This can be attributed to the asymmetry of the jets from upper and lower primary holes due to the different axial positions of these two primary holes in the reverse-flow combustor. The jets from primary holes truncate the primary flow and interact with the rotational flow, thereby forming asymmetrical upper and lower primary recirculation zones in the head part. The recirculating vortex promotes the mixing between the fuel and the fresh air. On the other hand, due to the low velocity of the recirculating vortex, the retention time of fuel can also be prolonged with violent combustion. Moreover, the flame is truncated due to the truncation of jets from primary holes. Due to upper and lower recirculating vortexes, a small part of fuel and gas still spreads downstream. This reduces the local equivalent ratio and weakens the flame intensity. In the combustion wake region, combustion flames move downstream near the wall surfaces, where flame intensity is significantly weakened. This is because a small portion of unburned fuel moves downstream with the primary flow after the combustion in the primary zone. Furthermore, flow pressure around the high-velocity jet column of upper and lower main holes drops, which leads to vortex flowing downstream near the wall of the main hole. This can provide necessary flow conditions for stable combustion, which is consistent with the flow field structure shown in Figure 5. In addition, as the equivalent ratio increases, the fuel distribution concentration increases, the combustion reaction can be accelerated, and the flame front moves forward. Furthermore, the reaction intensity in the combustion wake zone increases and the flame stabilizes. Accordingly, the flame in the Zone 3 becomes thinner and shorter.

Based on the extracted boundary contours of averaged and images, the areas of and at three different equivalent ratios can be calculated and normalized according to the area of the entire image:

Standard deviation () of the areas of two intermediate products in continuous image series at different operating conditions (Figure 11) is calculated as . Overall, as increases (the combustor is still under lean oil combustion at that moment), the reaction can be accelerated, and fuel undergoes complete combustion. This leads to a decline in flame length and the thinning of the flame surface. Accordingly, the proportion of areas of two combustion intermediate products decreases, i.e., the areas of and gradually drop with an increase in the equivalent ratio. As the equivalent ratio increases, the combustion is enhanced, both reaction and heat release effect become more compact and stable, and the ratios of area to area at three equivalent ratios are increased by 3.67%, 3.02%, and 2.13%, respectively. The standard difference of the intermediate product area proportion in continuous image series, denoted as S, can be used for reflecting the difference between combustion stability and various operating conditions at different equivalent ratios. Since the turbulent combustion in the combustor is a non-steady process, random vortexes are produced under the interaction between the jest from primary holes and the rotational flows in the head. This can impose disturbance on the combustor flames, thereby affecting both reaction and heat release zones. As the equivalent ratio increases, the standard deviations of continuous image areas of and drop from 0.0143 and 0.0132 to 0.0115 and 0.0109, respectively. This suggests that the combustion stability is enhanced with an increase in the equivalent ratio. Heat release hysteresis is present compared with combustion chemical reaction, and it is subjected to multiple factors including energy, turbulent dissipation, and thermal radiation. Therefore, the standard deviation of the heat release zone is smaller than the one of the reaction zone. Thus, it can be concluded that equivalent ratio change affects the sizes of the combustion reaction zone, heat release zone, and combustion stability under the same airflow conditions.

By applying the normalization process on the brightness of all pixels in time-averaged images at different equivalent ratios, the brightness distributions of and along different axial directions of the combustor at different equivalent ratios are obtained, as shown in Figure 12. Since the combustion in the primary recirculation zone is most intense and it weakens downstream along with the combustion, the mean brightness of and first increases and then drops from the outlet of the swirler to the downstream direction of the combustor. Brightness distributions of and also fluctuate. In the upstream of the fluctuation zone (i.e., the primary combustion flame zone), the mean brightness of and at different axial positions increases with the equivalent ratio. This may be because the increasing equivalent ratio can intensify the chemical reaction, thereby leading to an increase in the spontaneous radiation intensities of and . Since the fluorescence signal of the spontaneous radiation of is higher than that of under the same conditions, the intensity of exceeds that of at the same axial position. In addition, two among three curves in the fluctuation zone are intersectant. The intersection points along the flow direction in the combustor are denoted as P1, P2, and P3. By taking Point P2 as the reference point, the distance between the axial position corresponding to the curve peak and the axial distance between the axial position corresponding to the curve peak and P2 is denoted as . It can be observed that is larger than . On the one hand, the position of peak lags behind that of peak. Moreover, as the flame propagation intensity weakens towards the downstream part of the combustor, the corresponding axial position of in the fluctuation zone moves more forward than the position of . This suggests that the production of is significantly affected by the temperature. Accordingly, can be adopted as an important sign of combustion heat release.

4. Conclusions

This study experimentally investigated the flow field, ignition, and flame structure characteristics of a reverse-flow combustor. The following conclusions can be drawn: (1)As total pressure loss coefficient varies from 1% to 5%, the flow in the combustor is under a self-modeling state. The time-averaged flow field structure remains almost unchanged. However, the flow velocity in the combustor is increased(2)When =3%, the ignition time decreases from 26 ms to 16 ms as the equivalent ratio increased, suggesting the enhancement of ignition performance. The entire ignition process can be divided into three phases: turbulent dissipation phase, quasi-stable phase, and combustion-enhanced phase. The change of mainly affects the times consumed in the quasi-stable phase and the combustion-enhanced phase (T2 and T3), respectively. In particular, the quasi-stable phase is an important factor that determines whether the combustor can be successfully ignited. As increases from 0.6 to 0.66 and 0.72, T2 drops by 31.6% and 44.7%, respectively. On the other hand, T3 can be reduced by 20% and 30%, respectively. And it can be founded that flow field of the combustor is tightly correlated with ignition performance and flame propagation behavior(3)The distribution of , with close relation to the combustion reaction, can be used for characterizing the flame front position. An increase in the equivalent ratio can intensify the chemical reaction and enhance combustion stability, resulting in a decrease of and areas. Standard deviations of the continuous image areas of and drop from 0.0143 and 0.0132 to 0.0115 and 0.0109, respectively. For the same equivalent ratio, the axial position of peak point moves backward further than the axial position of the peak point. Moreover, the axial position corresponding to in the fluctuation zone moves more forward than the axial position of . This suggests that the distribution of is significantly affected by temperature. As such, it is an important sign of heat release in combustion

Nomenclature
PIV:Particle image velocimetry
:Equivalence ratio
Wfuel:Fuel mass flow rate
FAR:Fuel-air ratio
Pt3-Pt4:Inlet/outlet total pressure
:Liner pressure loss
PRZ:Primary recirculation zone
RP-3:One kind of liquid kerosene
Wair:Air mass flow rate
FARST:Stoichiometric ratio
umax:Inlet velocity at with 5%.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors would like to acknowledge support from the National Science and Technology Major Project (Grant Number 2017-III-0002-0026), National Science and Technology Major Project (Grant Number 2019-III-0014-0058), and the preprint as per the following link: https://papers.ssrn.com/sol3/papers.cfm?abstract_id=4086466.