• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Criteria for hypersonic airbreathing propulsion and its experimental verification

    2021-04-06 02:09:38ZonglinJIANGZijianZHANGYunfengLIUChunWANGChangtongLUO
    CHINESE JOURNAL OF AERONAUTICS 2021年3期

    Zonglin JIANG, Zijian ZHANG, Yunfeng LIU, Chun WANG,Changtong LUO

    a State Key Laboratory of High Temperature Gas Dynamics, Institute of Mechanics, Chinese Academy of Sciences, Beijing 100190, China

    b School of Engineering Science, University of Chinese Academy of Sciences, Beijing 100049, China

    KEYWORDS Concept demonstration;Hypersonic propulsion;Oblique detonation engine;Scramjet;Shock wave;Supersonic combustion;Wind tunnel test

    Abstract Hypersonic airbreathing propulsion is one of the top techniques for future aerospace flight, but there are still no practical engines after seventy years’ development. Two critical issues are identified to be the barriers for the ramjet-based engine that has been taken as the most potential concept of the hypersonic propulsion for decades. One issue is the upstream-traveling shock wave that develops from spontaneous waves resulting from continuous heat releases in combustors and can induce unsteady combustion that may lead to engine surging during scramjet engine operation.The other is the scramjet combustion mode that cannot satisfy thrust needs of hypersonic vehicles since its thermos-efficiency decreases as the flight Mach number increases.The two criteria are proposed for the ramjet-based hypersonic propulsion to identify combustion modes and avoid thermal choking. A standing oblique detonation ramjet (Sodramjet) engine concept is proposed based on the criteria by replacing diffusive combustion with an oblique detonation that is a unique pressure-gain phenomenon in nature. The Sodramjet engine model is developed with several flow control techniques, and tested successfully with the hypersonic flight-duplicated shock tunnel.The experimental data show that the Sodramjet engine model works steadily, and an oblique detonation can be made stationary in the engine combustor and is controllable.This research demonstrates the Sodramjet engine is a promising concept and can be operated stably with high thermal efficiency at hypersonic flow conditions.

    1. Introduction

    Many 21st century challenges exist in aviation science and technology, and one of these is the hypersonic vehicle from the dream for human beings to fly faster, higher and further than ever. With hypersonic airliners, we can arrive anywhere in the world within two hours. With reusable transatmospheric planes, we can take off horizontally from an airport runway,accelerate into orbit around the Earth,then reenter into the atmosphere, and finally land at an airport. In this way,space access will become reliable,routine and affordable.However,those exciting projects turn out to be very difficult to achieve.1,2It took about sixty years from Wright Brothers’first try in 1903 to the Concorde supersonic airliner taking off in 1969. During the period, we started from the very beginning of aviation, and entered successfully into the supersonic era.Now,seventy years have passed since the first object of human origin achieved a flight speed being five times faster than the speed of sound in 1949, and the practical hypersonic flight is still far ahead.3

    There are several critical techniques for developing hypersonic vehicles,and the most challenging one must be the hypersonic airbreathing engine. The turbofan engine works well for a supersonic flight of the Mach numbers less than 3, but cannot be operated efficiently for hypersonic flights.In 1935,Rene′Leduc of France was issued a patent on a piloted aircraft propelled by a ramjet concept. Later, the supersonic combustion ramjet (scramjet) engine concept was proposed for hypersonic propulsion. The research on the scramjet engine has spread rapidly over the world, and the word ‘supersonic combustion’is accepted to describe the flow physics in the scramjet combustor. The continuous exploration on the scramjet engine has been lasting for about 70 years since then,so far it still remains as one of the hottest research topics in gas dynamic realm.Many scientific research papers reported progresses on supersonic combustion and scramjet engines, including combustion physics, thermal-efficiency and engine design principals that were investigated with numerical simulations, experiments,and even flight tests.4-9However, there is still no practical scramjet engine in use for aerospace engineering.

    By recalling the research progress achieved in recent decades, several physical issues are considered to play an important role in scramjet development. Two of those issues are fundamental and critical.The first one is the wave propagation in scramjet flow passages due to spontaneous combustion waves arising from sudden heat release, which may develop into upstream traveling shock waves. This phenomenon can excite unstable combustion and results in inlet unstart and engine thermal choking.10-12The other is the low engine thrust that was demonstrated with flight tests,especially at high flight Mach numbers.It is possible to get positive thrusts for a small flight-testing vehicle, but the thrust is too small to power any practical hypersonic airliner for commercial applications.8,9.

    Another candidate for hypersonic propulsion is so-called the Standing Oblique Detonation (SOD) engine. For the SOD engine, an oblique detonation is applied to replace the diffusion-dominated combustion in scramjet combustors. The SOD engine had been theoretically demonstrated to be of the high-power density, the short combustor length and the simple engine structure.13,14.Ostrander et al.further calculated the specific impulse at the similar flow condition and found that it is higher than the scramjet engine.15Yuan and Huang did a parameter comparison of the total pressure, entropy and exergy at combustor outlets, and their theoretical results showed that the performance of the SOD mode under condition of the C-J detonation is the best among supersonic combustion modes and its entropy increase at the nozzle outlet is also the lowest one.16Other research results showed that the SOD mode can be operated in a wide flight Mach numbers ranging from 6-16.13In order to make an oblique detonation stationary,more papers dedicated to detailed detonation front structures and its evolution, and show the oblique detonation develops via three stages: initiation, transition and the fulldeveloped one where transverse shock waves exist like normal detonation waves.17,18.However,there is no extensive research carried out like scramjet engines in the world,especially,there is a lack of the experimental works on the SOD engine,even it is well known that the detonation-driven engine has potential advantages over the scramjet.

    Two problems may be responsible for such a situation.One is that the standing oblique detonation in combustors will take a complex structure developing from shock reflections, shock/boundary interaction and combustion instability due to the confined flow passage, therefore, it becomes a very difficult task to make an oblique detonation stationary at the required position in the SOD combustor.19The other is the lack of the proper wind tunnel being capable of reproducing true hypersonic flight conditions that are necessary for testing the SOD engine in its full scale.Therefore,the ground facility for hypersonic propulsion testing has been a challenge for Mach numbers higher than 8 for decades.20-22

    When the inlet flow Mach number is getting higher and higher, the engine combustion mode may transit from subsonic, supersonic combustion to detonations at a ramjetbased engine with proper flow control techniques. So, both the scramjet and the SOD engines belong to the same class of the thermal engines. Only difference between the two engines is the operation manner with which how the combustion is organized. Therefore, this paper is dedicated to the investigation into the following four issues. The first issue is the generation of the upstream-traveling shock waves developing during scramjet combustion and its shock-enhanced mechanism. The phenomenon is important because it is the key issue related to scramjet operation instability.The second one is on the criteria for hypersonic propulsion engineers to identify combustion modes and avoid unsteady combustion. The third is the feasibility of the Sodramjet engine that is operated by maintaining a standing oblique detonation wave in its combustor. The last is on experimental verifications to confirm the concept and criteria proposed in this paper.These may not be all the key physical issues that have puzzle us a lot for decades during the practical hypersonic propulsion development, but are real the significant and fundamental ones being of worthy of much attention.

    2. Generation of upstream-traveling shock waves

    In the 1950s and 1960s, a variety of experimental scramjet engines had been built, and tested with hypersonic propulsion test facilities in the USA.7,22The significant accomplishment on the Hypersonic Research Engine had been achieved under the support of the National Aerospace Plane program(NASP)since the X-30 was to use scramjet engine as its propulsion.The engine model with flight weight is re-generatively cooled and of the aero-thermodynamic integrated engine configuration.Unfortunately,the X-30 concept demonstrator was flown only in wind tunnels and the program was canceled in 1994.Low maturity of its propulsion technology is one of the many problems.

    A milestone was reached on the day of 27 March 2004 when a scramjet engine fueled with hydrogen was tested with the X-43A, a small hypersonic flight-testing vehicle, and the flight test was successfully completed. ‘‘The X-43A was a turning point,” says Jim Pittman, the principal investigator for hypersonic at NASA. ‘‘We learned two things: Scramjets really do work. You really can get positive thrust out of a scramjet and you really can integrate a scramjet with a vehicle that you can fly and control. And both of those things are huge”.There are also other two important issues that were revealed with the flight test data. The first issue is that the acceleration of the X-43A is small during the Mach number 7 flight and reduces to zero for Mach number 10.9The second one is the identification of the strong shock wave that exists in the scramjet inlet because a sharp jump was observed from the pressure distribution along the axial station through the scramjet engine. This is the upstream-traveling shock wave that had been observed frequently in scramjet ground testing and may result in inlet unstart and unstable engine operation.Pittman’s two things are the brief summary for the past hypersonic research, and these two issues are also huge for future hypersonic engine development.23

    When a diffusion-dominated combustion takes place in open air or in subsonic flows, spontaneous waves generated from sudden heat releases will propagate radially at the local sound speed, decay and disappear quickly without drawing anyone’s attention. However, when the combustion occurs in supersonic flows, shock waves are generated from the nonlinear propagation of compression waves and enhanced by the confined wall reflection.These acoustic waves were widely recognized in both numerical and experimental works.10,11Actually,only an exceedingly small fraction of the chemical energy released in the combustion process is required to generate large excursions of acoustic oscillations, and the relevant combustion instabilities plagued the development of air-breathing propulsion systems.12In scramjet engines, an upstreamtraveling shock wave generated from these waves may result in severe disruption in engine operation.

    In order to identify the so-called upstream-traveling shock wave,a scramjet test was completed with the Performance Test Engine (PTE) model in the hypersonic-flight-duplicated wind tunnel (JF-12 shock tunnel). Fig. 1 presents the PTE model and its installation in test section. The JF-12 shock tunnel operated with the backward detonation driver is 268 m in total length with a 2.5 m nozzle,and can be used for reproducing the hypersonic flight condition from Mach number 5-9 at altitudes of 25-50 km.20,21The PTE model is 2.2 m long and consists of an inlet, a combustor, and an expansion nozzle. The inlet is three-dimensional one with two compression ramp surfaces.The combustor is made of with two backward steps, two fuel injections are distributed within the steps and the fuelinjected mass is adjustable for reaction heat release control.The engine nozzle is two-dimensional and 1.5 m in length.The test flow Mach number is 7 with a total flow temperature of 2200 K and the maximum dynamic pressure of 50 kPa.Forty pressure transducers are distributed from the inlet to the engine nozzle to monitor pressure variations during combustion. Two high-speed cameras are mounted toward test windows to observe both combustion products from the nozzle and the shock wave structure around the inlet.

    Fig. 1 Photo of performance test engine model and its installation in JF-12 shock tunnel.

    Fig. 2 Pressure variations recorded along scramjet engine model.

    Fig. 3 Time-sequential photos taken with high-speed cameras during engine surging and six frames showing about one cycle.

    Two typical pressure(p)variations recorded along scramjet engine model are plotted in Fig. 2, from which two operation modes are found for the PTE model. From the upper half of the figure, it is observable that the pressure profile is flat with slight pressure perturbations. The experimental data indicate that the PTE model is working continuously and a stable combustion is maintained. From the lower half of Fig. 2,the periodical pressure oscillations are identified. To investigate further the pressure oscillations, the time-sequential photos are presented in Fig. 3 and these frames were taken with high-speed cameras during the PTE model operation. From this figure, it can be observed that the combustion dies out and is re-ignited periodically. The average frequency is about 220 Hz by counting the peaks of the oscillatory pressure. It seems that the combustion is dying out slower than its reignition by checking the frame time. The Equivalence Ratio (ER)is only the difference between these two tests. This important phenomenon could be defined as the engine surging that is a critical issue for scramjet operation control because there would be no practical engine that is allowed to operate in such an unstable mode.

    Generally speaking, there exists a critical operation point for any scramjet engine at certain inflow conditions.Once heat release is over the critical operation point, upstream-traveling compression waves will get stronger and stronger, and then a shock wave is generated and accelerated due to the temperature gradient distribution ahead.Therefore,a thermal choking in engines takes place once the upstream-traveling shock wave reaches to the inlet entrance, the inlet goes into unstarts and the combustion flame dies out. Once the flame dies out, the upstream-traveling shock wave will decay rapidly so that the inlet would restart and the PTE model operates again. This phenomenon is taken as the hypersonic engine surging that is very interesting to understand combustion physics, but a quite tough problem for practical engine development. There must be a balance point at which the inflow Mach number matches with the reaction heat release, and then the shock wave becomes stationary in combustors.

    3. Shock-enhanced mechanisms

    Acoustic waves generated in nature due to sudden heat releases will decay quickly in open space. However, the physical phenomenon will behave quite differently due to wave reflection when combustion reactions occur in a scramjet engine flow passage that is a finite-confined space. Three physical issues play important roles in the shock/reaction interaction and will be discussed by assuming a supersonically-moving gas flow in a straight flow passage with a constant heat-releasing source.The first physical issue is a spherical shock wave around the area where the chemical heat is released continuously like what happens in the scramjet combustor. The flow in upstream is compressed because it is entering,and expands in downstream because it is leaving.The flow physics is the same to the Mach cone generation and the strength of the spherical shock wave depends on the amount of the released reaction heat, but the flow temperature in the heat-releasing region is very high due to the reaction. The second issue is the effect of the flow passage wall. The spherical shock wave is reflected soon from the passage wall after it is generated,and then develops rapidly into a planar one.The flow expansion behind the planar shock wave is much weaker than the spherical one so that the planar shock wave would not decay so easily. The third one is the temperature gradient in front of the upstream-traveling shock wave.The coming flow is compressed by a series of the oblique shock waves originating from the confined wall and its temperature is increased gradually from the inlet entrance to the combustor.Therefore,the temperature gradient is generated in the inlet, resulting local sound speed variations. From the above discussion, it is understandable that the upstream-traveling shock wave is generated from a series of combustion waves,enhanced by shock reflections from the confined wall and accelerated nonlinearly in the inlet due to the negative temperature gradient ahead. As a result, the shock wave is easier to generate and less to decay in the scramjet flow passage than in open space.

    To demonstrate the shock-enhanced mechanism discussed above, numerical results are presented to show the generation and propagation of upstream-traveling shock waves. The test case is a two-dimensional straight flow passage with a chemically-reacting source in its middle where the heat is released continuously at the rate being equal to one occurring in a scramjet combustor,as shown in Figs.4(a)and(b),respectively. The computational domain is 20 mm in width and 200 mm in length. The flow field is computed by solving the Euler equations and its pressure distributions are presented in Fig.4.Fig.4(a)is the result for the inlet flow Mach number of 2.5 and Fig. 4(b) is for the inlet flow Mach number of 4.5.For most of the experimental scramjet engines, their designed inflow Mach numbers are between 2.5 and 3.5. From this figure,it is observable that the shock wave is generated originally in a spherical form, and evolving into a planar one after wall reflections, and then propagates upstream. From those time sequential frames, the shock wave motion is much slower in Fig.4(b)than in Fig.4(a),but it still propagates upstream even though the inlet flow Mach number is as high as 4.5. This is why the low equivalence ratio is used frequently for most of scramjet experiments to maintain a stable combustion because the low heat release could reduce effectively the Mach number of the upstream-traveling shock wave.Therefore,for a ramjetbased propulsion engine, the inlet flow Mach number should match with the upstream-traveling shock wave to avoid inlet unstart and engine surging.

    Fig. 4 Upstream-traveling shock propagation within engine flow passage.

    4. Criteria for ramjet-based hypersonic propulsions

    The inlet flow Mach number is a key parameter for the scramjet engine to maintain stable combustion, and the widelyaccepted values range from 2.5 to 3.5 for the current scramjet engine models.At such the inlet flow Mach numbers,the static temperature at combustor entrances could arise to the autoignition level under hypersonic flight conditions. However, to develop practical hypersonic engines, the engine operation instability becomes a severe problem that has to be considered according to the discussion in the above chapter.

    To demonstrate the physical mechanism behind the engine operation instability in ramjet-based hypersonic propulsion systems, the one-dimensional steady flow with a constant heat addition is taken as an example.It is well known that the heat addition makes a flow approach the sonic state whether it is originally supersonic or subsonic. The maximum heat, qmax,required to drive the flow into the sonic state is given by the following equation:

    Fig. 5 The maximum heat required to drive one-dimensional flow to sonic state.

    where Ma is the flow Mach number,γ is the specific heat ratio,Cpis the constant pressure ratio and T01is the initial flow temperature. The maximum heat calculated with Eq. (1) is presented in Fig.5.From this figure,it is observable that the maximum heat needed in subsonic cases is much higher than in supersonic or hypersonic cases. For the subsonic cases, the chemical reaction heat contributes to both the flow dynamic energy and the gas temperature increases. For the supersonic cases,the heat is dedicated mainly to increase the flow temperature, therefore, the local sound speed increases quickly and the local Mach number decreases rapidly. This is the fundamental mechanism behind the unsteady combustion in scramjet engines and the maximum heat is the first critical criterion that could be used to distinguish the combustion modes in the ramjet-based hypersonic propulsion,that is,supersonic or subsonic in the heat release area after combustion reactions. It is necessary to point out that the upstream flow of the heat release area will become subsonic since the flow is coming and the downstream one will become supersonic because the flow is expanding.

    The maximum heat addition imposes a big barrier on the scramjet thrust because thermal engines are not allowed to operate near the thermal choking state, so that a low equivalence ratio is often used in scramjet engine testing to avoid subsonic combustion. Using the low equivalence ratio is an effective means to stabilize the engine operation, but results in even low scramjet thrust. It is understandable that the scramjet engine operation stability is closely related to the inlet flow Mach number and the released reaction heat. The higher is the inlet flow Mach number, the more the heat addition could be. The physical mechanism is also correct for all the ramjet-based hypersonic propulsion, unfortunately, this problem has not been well recognized so far. There has been the arguing about ‘‘supersonic” or ‘‘subsonic” combustion in hypersonic propulsion realm for decades, and Fig. 5 indicates the answer inspiringly. In the ramjet-based propulsion engine,the supersonic combustion may exist as long as the chemical reaction heat released is less than the maximum heat required to drive the flow into the sonic state.

    Once the chemical reaction heat released in ramjet-based engines is higher than the maximum heat, the subsonic combustion takes place. This combustion mode transition results in a shock wave that is generated from combustion waves and can propagate upstream. There must be another criterion at which the upstream-traveling shock wave can compete with the inlet flow. In other words, what is the critical inlet flow Mach number for a ramjet-based hypersonic propulsion engine?At such the Mach number,the engine can be operated stably at a full equivalence ratio without the shock-induced engine surging.This criterion is critical because there is no airliner that can accept an engine that may work in unstable combustion mode. To obtain this critical Mach number, the flow process in the ramjet-based hypersonic engine is simplified as a steady one-dimensional flow with continuous heat addition at a given reaction rate. During the heat-releasing process,the upstream-traveling shock wave will be generated after a sonic state is reached, and its Mach number depends on the amount of the reaction heat, gaseous media and its thermal state. Assuming a perfect gas, the critical Mach number for stable engine operation can be given by following after the C-J detonation theory.

    where MaCriis the Mach number of the upstream-traveling shock wave, and qmaxdefines the maximum amount of the chemical heat that can be released for a given combustible gas mixture. Subscript ‘‘0′′stands for the flow state before the heat addition and subscript ‘‘1” is for that after the addition.γ is the specific heat ratio and c is the sound speed.MaCriand qmaxcan be taken as two key parameters for the ramjetbased engine design and its operation. If the inlet flow Mach number is less than this critical value, the upstream-traveling shock wave will propagate into the inlet, and may result in not only unsteady combustion but also engine surging as experimentally demonstrated above.

    Fig. 6 Mach number, MaCri, variations with equivalence ratio for different fuels.

    The critical inlet flow Mach number is the second criterion proposed for hypersonic airbreathing propulsion,as defined by Eq. (2), and the criterion depends on both the reaction heat and the detonable gas mixture. Fig. 6 presents results of the critical Mach number, MaCrivarying with the equivalence ratio for different fuels. It can be seen that the Mach number is higher than 3 even the equivalence ratio is around 0.3.Unfortunately, for most scramjet engines, their designed inlet flow Mach numbers are between 2.5 and 3.5.This is the reason why many research papers reported that subsonic combustions are observed in their combustors,the relevant combustions are unstable and the upstream-traveling shock wave is observed at inlet entrances. So, the low equivalence ratio was chosen frequently to maintain stable combustions. Actually, the low equivalence ratio for stable engine operation is not acceptable because the engine thrust is reduced significantly since the total pressure loss has already been incurred during inlet flow compression.It is also possible to lower this critical Mach number by distributing the heat release sources along the combustor and expanding the flow while the gas mixture is reacting,however, the problem still exists when trying the full equivalence ratio to increase the engine thrust.

    The word, ‘supersonic combustion’ may not be an appropriate word that can be used to describe the flow physics in the ramjet-based hypersonic engines,and the phrase,‘combustion in supersonic flow’ seems more suitable to identify the chemically-reacting flows. In supersonic flows, the diffusive combustion results in successive acoustic waves propagating radially once it takes place if the first criterion is established,and the upstream travelling shock wave will be generated when the heat addition is higher than the maximum heat required to drive the flow into the sonic state. The more the heat releases,the stronger the shock wave develops. The shock wave is enhanced by the confined flow passage wall and further accelerated due to temperature gradients within engine inlets.There exists a competing mechanism between the supersonic coming flow and chemical reactions and the mechanism can be measured with the second criterion. In a ramjet-based engine, the inlet flow Mach number must match with qmaxfor stable combustion without engine surging. Eq. (2) presents the relationship that can be used to define the second criterion for the ramjet-based hypersonic propulsion without the upstreamtraveling shock wave that can propagate into the engine inlet and results in thermal choking.

    5. Verification of Sodramjet engine concept

    The second criterion defined with Eq. (2) presents the minimum inlet flow Mach number for the ramjet-based hypersonic propulsion. Over this Mach number, there is no shock wave that could propagate upstream into inlets,but a standing oblique detonation may be generated if an initiation source is provided in combustors.The higher the Mach number is,the small the angle of the standing oblique detonation becomes.The detonation is a shock-induced chemically-reacting wave that is a unique pressure-gain combustion in nature.With the standing oblique detonation, a nearly constant volume thermal-cycle can be organized for air-breathing hypersonic engines. Therefore, keeping the inlet flow Mach number higher than the criterion is absolutely necessary to make a detonation stationary.From this viewpoint, the criterion is actually a design parameter for the Sodramjet engine. In the engine, the coming gas flow is compressed to an auto-ignition level by the leading shock wave, the shock-induced reaction provides enough energy to support the leading shock wave in return.Therefore,the standing oblique detonation is not only stable, but also self-sustainable.The leading shock wave works like an efficient compressor and the chemical reaction coupling closely with it works like a turbine in modern turbojet engines. The Sodramjet engine belongs to the class of the ramjet-based hypersonic propulsion, but behaves quite different from the scramjet engine due to its high thermo-efficiency, simple combustor structure, low inflow compression loss, and stable engine operation.

    For the Sodramjet engine, the criterion indicates a minimum flow Mach number for inlet flows because it is necessary to generate an oblique shock wave with the post-shock temperature,Tigreaching to the auto-ignition level of combustible gas mixtures. According to the oblique shock wave formula, the relation between the required oblique shock angle, θ, and the inlet flow Mach number can be determined. The results calculated with the relation are presented in Fig. 7 where the ignition temperature, Tigis set to be three different values,respectively.From Fig.7,it is observable that for a given inlet flow Mach number, the higher the ignition temperature, the bigger the oblique shock angle becomes. It is also known that the overdriven detonation will cause more entropy increase than the C-J detonation so that choosing well-matched θ for detonation ignition could maximize thermal efficiency of the Sodramjet engine.

    Fig. 7 Oblique shock angle variations with inflow Mach number.

    Fig.8 Numerical result showing development of oblique shockinduced oblique detonation.

    Producing a proper oblique shock wave for ignition is the first step to make an oblique detonation stand in ramjet engines, and some physical issues need to pay attention for realizing this goal.23-25Fig.8 presents a temperature(T)distribution of an oblique detonation wave from the computational test case where the combustible gas is a hydrogen/air mixture,the inflow Mach number is 9, θ is taken to be 25° and the Navier-Stokes equations are solved to get solutions.From this figure, four physical phenomena could be identified. The first one is the flow region from the wedge tip to the triple point on the wave-front structure, which is called the transition region where the reaction from shock ignition transits to detonation. The second one is also a transition region where the detonation transiting from an overdriven state to a freepropagating state. This feature can be observed more clearly for an oblique detonation wave over a finite-length wedge.The third part is the cellar detonation and its macroscopic characters in the normal direction of the detonation wave front could be predicted with the C-J theory. The last one is the boundary layer where the chemically-reacting gas flow in the subsonic state is dominated.

    In order to design a wind tunnel test model for demonstrating the Sodramjet engine concept, a series of numerical simulations were carried out to investigate three issues that are important to make an oblique detonation stationary.17,26-28.The first issue is hypersonic pre-mixing since the inlet flow Mach number is usually much higher than 5.The fuel injection device must be moved into the inlet to enhance the diffusiondominated mixing process. The injection device, as shown in Fig. 9, utilizes struts with sharp edges to avoid the bow shock generation that may result in chemical reactions earlier than required. Injection holes are arranged to be normal to streamlines so that stream-wise vortices generated behind the injection device can be used to enhance the mixing. A long distance is necessary because the mixing process is dominated mainly by the diffusion mechanism. The second issue is the chemical reaction occurring in the boundary layer behind the first transition region. Although the chemical reaction is weak in the boundary layer,acoustic waves generated from the combustion will propagate upstream due to the local subsonic state and interfere with the standing oblique detonation. Therefore,a lean mixture in the boundary layer is a better choice when arranging injection holes. The third one is the boundary layer separation near the position where the oblique detonation stands. The separation bubble can induce a bow shock wave that may excite chemical reactions which can disturb the designed oblique detonation. The boundary layer bleeding is one way to solve the problem,and computational results show the idea works well. There may be other physical issues that are effective more or less in the Sodramjet engine design, but these issues are significantly important.

    According to the proposed criteria and three key issues discussed in the above chapters, the schematic of the Sodramjet engine model designed at Mach number 9 is shown in Fig. 9,where the OSW stands for oblique shock wave and the ODW for oblique detonation wave. The Sodramjet engine model is composed of three strut-injectors,a single-stage compression inlet,a combustor,and a nozzle.The dimension of the Sodramjet engine model is 2.2 m in length and 0.55 m in height. The inlet is composed of a 15°-inclined ramp that is 1.6 m in length. The combustor is 0.41 m in length and 0.0765 m in height. Obviously, the Sodramjet combustor is rather short by comparing with scramjets even it is to operate at the high Mach number. This feature is helpful for reducing both frictional force and thermal loads.Following the combustor is a short nozzle with a 15°expansion angle at one side,and 0.4 m in length. A bleed device is also equipped for boundary layer control to make oblique detonation stationary.

    The configuration of the concept demonstration model of the Sodramjet engine is shown in Fig.10(a)and its installation in the JF-12 shock tunnel test section is shown in Fig. 10(b).The JF-12 shock tunnel can provide the engine test with 100 ms test duration and its 2.5-meter nozzle can accommodate fully the test model in its uniform core flow region.20-22,29,30The front half of the Sodramjet combustor sidewall is replaceable with glass windows to ensure that the standing oblique detonation wave can be recorded with a high-speed camera. Pressure transducers are also installed along the test model to show pressure variations during engine operation.

    Experiments are carried out at a nominal Mach number of 9, and the oblique detonation standing in the combustor is observed to maintain a stable state for as long as 50 ms. One of the video frames is presented in Fig.11 with the corresponding hydrogen concentration from numerical simulations that show the hydrogen fraction distribution and pressure isolines together.From Fig.11(a),it is clearly observed that the standing oblique detonation wave exists in the combustor and the hydrogen burns out behind the leading shock in a very short distance, as shown in Fig. 11(b). This short distance from the leading shock wave to the line where the hydrogen is mostly consumed is the so-called chemically-reacting zone that appears to be coupled together with the oblique shock wave.The standing oblique detonation front is composed of two parts,the transition region and the fully-developed detonation.The transition region is from the originating point of the oblique shock wave to the position where the reacting front catches up with the oblique shock wave. This wind tunnel test confirms that the Sodramjet engine can be realized and is also controllable with three key techniques as discussed above.The concept of the oblique detonation engine had been proposed for many decades, but is demonstrated successfully with wind tunnels for the first time.The success benefits from the criteria proposed for the Sodramjet engine design and the flow control techniques also play a very important role in maintaining a standing oblique detonation at the required position.

    Fig. 9 Schematic of designed Sodramjet engine concept.

    Fig. 10 Concept demonstration model of Sodramjet engine and its installation in JF-12 shock tunnel.

    Fig.11 Experimental photo of standing oblique detonation in combustor of Sodramjet engine model(left)and corresponding hydrogen concentration from numerical simulations (right).

    The Sodramjet engine has several advantages over the ramjet-based propulsion engines. The first advantage is thermal efficiency. The thermodynamic process taking place in the Sodramjet engine can be taken as a nearly constant volume cycle, and the process in the ramjet-based engine is often considered to be a constant pressure one. Theoretical analysis shows that the thermal efficiency of the constant volume cycle is about 50%higher than the constant pressure cycle.11,14The second one is the combustor size. The combustor length is 410 mm for the current Sodramjet engine model but about 2000 mm or more for scramjet engines. This length reduction will cut down significantly the drag force acting on the combustor wall. The third one is the inflow compression loss.The designed inlet flow Mach number is between 2.5 and 3.5 for most of the scramjet engines, but the Sodramjet engine requires Mach numbers to be higher than 5 according to the proposed criterion. Therefore, the inlet flow compression loss of the Sodramjet engine is reduced greatly and the low total pressure loss contributes a lot to the engine efficiency. The fourth one is the thermal load reduction on the combustor.It is obvious that the short combustor will not only reduce the thermal load,but also prevent the heat loss from the combustor so that the engine efficiency would benefit.The last one is about combustion stability.The Sodramjet engine can work with the full equivalence ratio if the inlet flow Mach number meets the criterion. The oblique shock angle, θ will be smaller if the inflow Mach number increases, therefore, the standing oblique detonation could be automatically self-adjusted to be sustainable.Moreover,it is possible for us to control the wedge angle to make an oblique shock wave just strong enough to initiate an oblique detonation. The higher is the inlet flow Mach number, the smaller the θ becomes. The Sodramjet engine operating at near C-J detonation mode has the lowest entropy increase and the highest exergy at combustor outlets, which has never been reached by other kind engines before.14

    It is interesting to point out that the exhaust gas flow after the standing oblique detonation front is supersonic and the fuel is consumed at supersonic speeds. The combustion in the Sodramjet engine is really the supersonic combustion.Eq. (1) indicates the critical criterion at which the detonation takes place because the upstream-travelling shock is getting strong enough to excite chemical reactions. Therefore, a new type of hypersonic propulsion could be defined as ‘‘standing oblique detonation ramjet engine”, short for the Sodramjet engine. The engine belongs to the ramjet-based hypersonic propulsion, but its combustion mode is totally different from the scramjet engines. Theoretical analysis indicates that the Soramjet engine working as an air-breathing engine can offer a great potential to extend the flight Mach numbers from 6 to 16.11

    By recalling the development of the modern aviation industry,the reciprocating piston internal combustion engine is used for subsonic airplanes,and the turbojet engine is developed for supersonic flights. What kind of the engine is capable of powering hypersonic vehicles? 70 years’ exploration on hypersonic propulsion indicates that the revolutionary concept is really in need for hypersonic air-breathing engine development. The Sodramjet engine concept can be a very promising choice and the work presented here supports strongly this idea.

    6. Conclusions

    In this paper, the upstream-traveling shock wave and its enhanced mechanism are discussed, and two criteria proposed for developing hypersonic airbreathing propulsion and the Sodramjet engine concept are verified with wind tunnel tests.Progresses achieved are summarized as follows:

    (1) The upstream-traveling shock wave is identified to be the intrinsic shock wave of ramjet-based hypersonic propulsion engines. The shock wave occurs naturally due to continuous heat release in supersonic flows, and is enhanced by the engine flow passage walls and accelerated when propagating upstream through the inlet where the temperature gradient exists. The more the reaction heat releases, the higher the shock Mach number becomes. The engine surging induced by this shock wave is disclosed with a scramjet model experiment in the JF-12 shock tunnel. The problem is not well recognized before, but must be solved when designing any practical hypersonic engine that is required to operate stably at a full equivalence ratio.

    (2) By considering upstream-travelling shock waves, two criteria are proposed for developing hypersonic airbreathing propulsion under the assumption of the perfect gas, instantaneous heat release and onedimensional steady flow. The first criterion can be used to identify the combustion mode, that is, subsonic or supersonic combustion in traditional saying.The second criterion is critical not only for scramjet engines to maintain stable combustion, but also indicates a minimum inlet flow Mach number for developing Sodramjet engines.

    (3) A concept demonstration model of the Sodramjet engine is designed based on the second criterion with the help of several flow control techniques and the experiment is carried out in the JF-12 shock tunnel.The experimental data show that the standing oblique detonation could be made stable and it is automatically self-adjusted to be sustainable. It is the first time to verify the oblique detonation engine concept successfully with wind tunnel tests and the standing oblique detonation undergoes the initiation,transition and fully-developed stages. The experiments also demonstrated well the criteria proposed in this paper.According to research progresses on oblique detonation engines,the Sodramjet engine is very promising for hypersonic airbreathing propulsion because that the engine system is simple in its structure,efficient in the thermo-cycle,broad in its operation Mach number range and low in aerothermodynamic heating load.

    Declaration of Competing Interest

    The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

    Acknowledgment

    The authors would like to thank all the members in Research Group of Shock Wave and Detonation Physics who help a lot in the JF-12 shock tunnel tests. The research is supported by the National Natural Science Foundation of China (No.11532014).

    一边摸一边抽搐一进一小说| 成年女人看的毛片在线观看| 熟妇人妻久久中文字幕3abv| 97超碰精品成人国产| 国产av一区在线观看免费| 免费黄网站久久成人精品| 别揉我奶头 嗯啊视频| 老师上课跳d突然被开到最大视频| 欧美另类亚洲清纯唯美| videos熟女内射| 秋霞伦理黄片| 国产日韩欧美在线精品| 国产日韩欧美在线精品| 大又大粗又爽又黄少妇毛片口| 久久人妻av系列| 亚洲国产成人一精品久久久| 国产一级毛片在线| 久久人人爽人人片av| 国产精品伦人一区二区| 久久人人爽人人片av| 女人被狂操c到高潮| 成人二区视频| 如何舔出高潮| 欧美成人精品欧美一级黄| 99热这里只有精品一区| 免费在线观看成人毛片| 成年免费大片在线观看| 午夜精品国产一区二区电影 | 亚洲美女搞黄在线观看| 亚洲av不卡在线观看| 18禁动态无遮挡网站| 精品久久久久久成人av| 国产亚洲av片在线观看秒播厂 | 国产精品一二三区在线看| 晚上一个人看的免费电影| 久久久久久久久久久丰满| 十八禁国产超污无遮挡网站| 精品久久久噜噜| 一级爰片在线观看| 精品人妻一区二区三区麻豆| 亚洲天堂国产精品一区在线| 久久久欧美国产精品| 成人午夜精彩视频在线观看| 国产在线男女| 国内精品一区二区在线观看| 高清午夜精品一区二区三区| 久久久亚洲精品成人影院| 2022亚洲国产成人精品| 久久精品91蜜桃| 天堂中文最新版在线下载 | 国产精品久久视频播放| 亚洲最大成人手机在线| 成人av在线播放网站| 欧美xxxx性猛交bbbb| 人妻系列 视频| 99热这里只有精品一区| 十八禁国产超污无遮挡网站| 91精品一卡2卡3卡4卡| 岛国毛片在线播放| 久久久久国产网址| 国产精品国产三级专区第一集| 亚洲中文字幕一区二区三区有码在线看| 久久99精品国语久久久| 91aial.com中文字幕在线观看| 亚洲国产精品专区欧美| 色综合亚洲欧美另类图片| 久久久久久伊人网av| 免费电影在线观看免费观看| 一级毛片电影观看 | 亚洲成人久久爱视频| 午夜激情福利司机影院| 亚洲一区高清亚洲精品| 亚洲av电影不卡..在线观看| 最近最新中文字幕大全电影3| 麻豆成人av视频| 三级经典国产精品| 九九在线视频观看精品| 全区人妻精品视频| av在线蜜桃| 波多野结衣高清无吗| 日本色播在线视频| 97超视频在线观看视频| 午夜亚洲福利在线播放| 成人无遮挡网站| 日韩亚洲欧美综合| 韩国av在线不卡| 日本三级黄在线观看| 18+在线观看网站| 不卡视频在线观看欧美| 成年版毛片免费区| 午夜福利在线观看免费完整高清在| 国产三级中文精品| 久久久久久久久大av| 国产一区二区在线av高清观看| 亚洲av电影在线观看一区二区三区 | 成人二区视频| 亚洲国产精品国产精品| 人人妻人人看人人澡| 国产精品国产高清国产av| 国产老妇女一区| 成年女人看的毛片在线观看| 少妇丰满av| 久久人人爽人人爽人人片va| 天美传媒精品一区二区| 国模一区二区三区四区视频| 久久久国产成人免费| 日韩一区二区三区影片| 亚洲人成网站在线播| 一级毛片aaaaaa免费看小| 日韩,欧美,国产一区二区三区 | 久久精品久久精品一区二区三区| 18禁在线无遮挡免费观看视频| 国产片特级美女逼逼视频| 久久草成人影院| 一级毛片aaaaaa免费看小| 色综合色国产| 两个人的视频大全免费| 国产精品福利在线免费观看| 亚洲精品乱久久久久久| 69人妻影院| 国产精品久久久久久精品电影小说 | 日韩强制内射视频| 色播亚洲综合网| 午夜福利视频1000在线观看| 嘟嘟电影网在线观看| 国产色婷婷99| 亚洲最大成人中文| 国产午夜福利久久久久久| 国产一区二区三区av在线| 日韩大片免费观看网站 | 青春草国产在线视频| 黄色欧美视频在线观看| .国产精品久久| 日韩一本色道免费dvd| 国产精品久久久久久精品电影| 国产精品日韩av在线免费观看| 免费观看性生交大片5| 大香蕉97超碰在线| 成人性生交大片免费视频hd| 九九爱精品视频在线观看| 看黄色毛片网站| 九九在线视频观看精品| 日本色播在线视频| 丰满人妻一区二区三区视频av| 色哟哟·www| 天天躁日日操中文字幕| 久久久国产成人精品二区| 丝袜喷水一区| 亚洲国产日韩欧美精品在线观看| 亚洲欧美精品自产自拍| 99久久精品热视频| 99久久无色码亚洲精品果冻| 国产视频内射| 国产欧美日韩精品一区二区| 三级经典国产精品| 亚洲成人中文字幕在线播放| 国产成人福利小说| 亚洲四区av| 国产高清三级在线| 一级黄色大片毛片| 搡女人真爽免费视频火全软件| 一级二级三级毛片免费看| 国产69精品久久久久777片| 国产在线男女| 久久久亚洲精品成人影院| 99热网站在线观看| 日韩欧美精品免费久久| 日韩高清综合在线| 中文字幕制服av| 国产精品99久久久久久久久| 在线天堂最新版资源| 久久久久精品久久久久真实原创| 偷拍熟女少妇极品色| 国产中年淑女户外野战色| 大又大粗又爽又黄少妇毛片口| 欧美日韩一区二区视频在线观看视频在线 | 亚洲乱码一区二区免费版| 麻豆国产97在线/欧美| 免费观看的影片在线观看| 国产在线一区二区三区精 | 超碰av人人做人人爽久久| 尤物成人国产欧美一区二区三区| 国内少妇人妻偷人精品xxx网站| 99九九线精品视频在线观看视频| 最近最新中文字幕大全电影3| 久久久午夜欧美精品| 国产免费视频播放在线视频 | 高清午夜精品一区二区三区| 久久精品久久精品一区二区三区| 一级毛片aaaaaa免费看小| 久久99蜜桃精品久久| 毛片一级片免费看久久久久| 高清视频免费观看一区二区 | 日韩欧美精品免费久久| 欧美人与善性xxx| 免费电影在线观看免费观看| 男人狂女人下面高潮的视频| 久久热精品热| 卡戴珊不雅视频在线播放| 亚洲丝袜综合中文字幕| 婷婷色综合大香蕉| 一级黄片播放器| 网址你懂的国产日韩在线| 久久精品综合一区二区三区| 国产在线一区二区三区精 | 国产精品嫩草影院av在线观看| 欧美日韩国产亚洲二区| АⅤ资源中文在线天堂| 国产乱人偷精品视频| 亚洲精品日韩在线中文字幕| 久久午夜福利片| 国产伦精品一区二区三区四那| 久久精品夜夜夜夜夜久久蜜豆| 麻豆成人午夜福利视频| 小说图片视频综合网站| 青春草亚洲视频在线观看| 久久久久久久久久久丰满| 日日啪夜夜撸| 国产淫片久久久久久久久| 不卡视频在线观看欧美| 男女边吃奶边做爰视频| 国产综合懂色| 少妇丰满av| 在现免费观看毛片| 波多野结衣巨乳人妻| 男女啪啪激烈高潮av片| 91久久精品国产一区二区三区| 丝袜美腿在线中文| 欧美成人精品欧美一级黄| 大香蕉久久网| 国产av在哪里看| 激情 狠狠 欧美| 99热这里只有是精品在线观看| 高清视频免费观看一区二区 | 免费不卡的大黄色大毛片视频在线观看 | 看非洲黑人一级黄片| 国产免费视频播放在线视频 | 日韩视频在线欧美| 边亲边吃奶的免费视频| 成人特级av手机在线观看| 久久人人爽人人片av| 日本熟妇午夜| АⅤ资源中文在线天堂| 国产色爽女视频免费观看| 国产 一区精品| 成人av在线播放网站| 黄色日韩在线| 日日摸夜夜添夜夜爱| 久久久久久久久中文| 最近最新中文字幕免费大全7| 一级av片app| 国产免费一级a男人的天堂| 岛国在线免费视频观看| 久久久欧美国产精品| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 麻豆国产97在线/欧美| 久久精品国产鲁丝片午夜精品| 蜜臀久久99精品久久宅男| 国产亚洲av嫩草精品影院| 成人高潮视频无遮挡免费网站| 男女边吃奶边做爰视频| 中文字幕久久专区| 99久久中文字幕三级久久日本| 国产一级毛片在线| 汤姆久久久久久久影院中文字幕 | 国产精品一二三区在线看| 国内精品美女久久久久久| 国产免费视频播放在线视频 | 久久99蜜桃精品久久| 午夜精品在线福利| av.在线天堂| 97在线视频观看| 国产精品蜜桃在线观看| 精品99又大又爽又粗少妇毛片| 国产精品蜜桃在线观看| 亚洲成色77777| 男女国产视频网站| 一本一本综合久久| 天堂av国产一区二区熟女人妻| 亚洲激情五月婷婷啪啪| 午夜亚洲福利在线播放| 天堂影院成人在线观看| 欧美又色又爽又黄视频| 日韩高清综合在线| 中国美白少妇内射xxxbb| 色综合亚洲欧美另类图片| 国产麻豆成人av免费视频| 亚洲成色77777| 亚洲欧洲日产国产| 国产中年淑女户外野战色| 久久久久久久久久久免费av| 亚洲av.av天堂| 久久久色成人| 免费人成在线观看视频色| 日本午夜av视频| 一本久久精品| 国产午夜精品一二区理论片| 日韩精品青青久久久久久| 免费大片18禁| 直男gayav资源| 亚洲中文字幕一区二区三区有码在线看| 中文天堂在线官网| 欧美性猛交╳xxx乱大交人| 九九爱精品视频在线观看| 国产免费福利视频在线观看| 最近视频中文字幕2019在线8| 中文字幕亚洲精品专区| 国产熟女欧美一区二区| 一级毛片电影观看 | 亚洲av中文字字幕乱码综合| www日本黄色视频网| 插阴视频在线观看视频| 日韩在线高清观看一区二区三区| 国产成人91sexporn| 特大巨黑吊av在线直播| 国产一区二区三区av在线| 久久精品国产亚洲网站| 亚洲真实伦在线观看| 国产美女午夜福利| 在线观看av片永久免费下载| 国产成人aa在线观看| 国产成人精品久久久久久| 插阴视频在线观看视频| av.在线天堂| 一边摸一边抽搐一进一小说| 小说图片视频综合网站| 国产黄色小视频在线观看| 日日摸夜夜添夜夜添av毛片| 青春草亚洲视频在线观看| 波多野结衣巨乳人妻| 少妇的逼好多水| 男女那种视频在线观看| 国产成年人精品一区二区| 国产大屁股一区二区在线视频| 草草在线视频免费看| 三级国产精品片| 久久欧美精品欧美久久欧美| 一级黄片播放器| 日韩av不卡免费在线播放| 成人性生交大片免费视频hd| 日韩av不卡免费在线播放| 午夜爱爱视频在线播放| 日韩av在线免费看完整版不卡| 亚洲国产精品成人久久小说| 国产乱人视频| 午夜激情欧美在线| 亚洲精品,欧美精品| 国产精品一区www在线观看| 免费观看在线日韩| 亚洲欧美日韩高清专用| 永久免费av网站大全| 18禁在线无遮挡免费观看视频| 国产黄片美女视频| 成人无遮挡网站| 国产精品国产高清国产av| 观看美女的网站| 国产精品不卡视频一区二区| 观看美女的网站| 国产午夜精品论理片| 狂野欧美激情性xxxx在线观看| 免费一级毛片在线播放高清视频| 国产亚洲最大av| 亚洲av中文字字幕乱码综合| 亚洲av电影在线观看一区二区三区 | 成人漫画全彩无遮挡| a级一级毛片免费在线观看| 男女边吃奶边做爰视频| 日韩在线高清观看一区二区三区| 一级二级三级毛片免费看| 尾随美女入室| 久久99热这里只频精品6学生 | 亚洲国产欧美在线一区| av专区在线播放| 床上黄色一级片| 日韩一区二区视频免费看| 亚洲精品色激情综合| 嘟嘟电影网在线观看| 大香蕉久久网| 国产午夜精品一二区理论片| eeuss影院久久| 男插女下体视频免费在线播放| 97超碰精品成人国产| 国产精品国产三级国产av玫瑰| 亚洲一区高清亚洲精品| 国产女主播在线喷水免费视频网站 | 黄色配什么色好看| 高清av免费在线| 别揉我奶头 嗯啊视频| 亚洲精品亚洲一区二区| 卡戴珊不雅视频在线播放| 日韩在线高清观看一区二区三区| 亚洲怡红院男人天堂| 欧美3d第一页| 国产精品一区二区性色av| 天天躁日日操中文字幕| 麻豆成人av视频| 高清av免费在线| 99久久人妻综合| 亚洲自拍偷在线| 熟女人妻精品中文字幕| 2022亚洲国产成人精品| 免费看日本二区| 国产一区二区在线观看日韩| 又爽又黄a免费视频| 菩萨蛮人人尽说江南好唐韦庄 | 国产毛片a区久久久久| h日本视频在线播放| 国产精品综合久久久久久久免费| 一区二区三区乱码不卡18| 精品少妇黑人巨大在线播放 | 亚洲国产精品国产精品| 亚洲欧洲国产日韩| 亚洲欧美精品综合久久99| 国产精品无大码| 亚洲成人久久爱视频| 精品一区二区三区人妻视频| 少妇猛男粗大的猛烈进出视频 | 久久精品国产99精品国产亚洲性色| 午夜福利在线在线| 精品人妻熟女av久视频| 三级国产精品欧美在线观看| 一级毛片我不卡| 亚洲经典国产精华液单| 一级毛片久久久久久久久女| 成人午夜精彩视频在线观看| 99久久中文字幕三级久久日本| 国产真实乱freesex| 国产精品.久久久| 亚洲精品国产成人久久av| 中文字幕免费在线视频6| 亚洲av中文av极速乱| 麻豆成人午夜福利视频| 成人毛片60女人毛片免费| 亚洲四区av| 日本-黄色视频高清免费观看| 两个人视频免费观看高清| 日韩,欧美,国产一区二区三区 | 国产白丝娇喘喷水9色精品| 观看美女的网站| 国产精品人妻久久久久久| 国产亚洲午夜精品一区二区久久 | 亚洲av熟女| 特级一级黄色大片| 综合色丁香网| 欧美激情国产日韩精品一区| 一级二级三级毛片免费看| 建设人人有责人人尽责人人享有的 | 国产精品无大码| 18禁动态无遮挡网站| 亚洲av男天堂| 午夜视频国产福利| 综合色丁香网| 欧美区成人在线视频| 亚洲国产色片| 亚洲成人精品中文字幕电影| 精品久久久久久成人av| 亚洲熟妇中文字幕五十中出| 久久久久免费精品人妻一区二区| 六月丁香七月| av在线蜜桃| 久久久久性生活片| 婷婷六月久久综合丁香| 一级毛片aaaaaa免费看小| 激情 狠狠 欧美| 日本熟妇午夜| 男人的好看免费观看在线视频| 边亲边吃奶的免费视频| 亚洲精品456在线播放app| 精品一区二区免费观看| 免费一级毛片在线播放高清视频| 精品免费久久久久久久清纯| 在线天堂最新版资源| 最近2019中文字幕mv第一页| eeuss影院久久| 丝袜喷水一区| 我要看日韩黄色一级片| 国产精品1区2区在线观看.| 大香蕉97超碰在线| 日本与韩国留学比较| 女的被弄到高潮叫床怎么办| 美女黄网站色视频| 欧美日本视频| 午夜福利在线观看吧| 三级国产精品欧美在线观看| 女人被狂操c到高潮| 久久久久精品久久久久真实原创| 亚洲三级黄色毛片| 久久99热这里只频精品6学生 | 69av精品久久久久久| 99视频精品全部免费 在线| 国产美女午夜福利| 最新中文字幕久久久久| 欧美潮喷喷水| 嫩草影院新地址| 亚洲国产高清在线一区二区三| 99热6这里只有精品| 国内揄拍国产精品人妻在线| 国产色婷婷99| 日本爱情动作片www.在线观看| 欧美日韩一区二区视频在线观看视频在线 | 午夜福利网站1000一区二区三区| 欧美3d第一页| av免费在线看不卡| 97超碰精品成人国产| 美女内射精品一级片tv| 亚洲精品乱码久久久久久按摩| 国产精品蜜桃在线观看| 国产亚洲91精品色在线| 天天一区二区日本电影三级| 老司机影院成人| 欧美日本亚洲视频在线播放| 老女人水多毛片| 青春草亚洲视频在线观看| 一级黄色大片毛片| 久久精品国产亚洲网站| 99热网站在线观看| 欧美高清成人免费视频www| 亚洲五月天丁香| 男的添女的下面高潮视频| 老司机福利观看| 高清av免费在线| 久久鲁丝午夜福利片| 亚洲欧美日韩高清专用| 亚洲精品成人久久久久久| 成人综合一区亚洲| 久久6这里有精品| 国产高潮美女av| 成人亚洲精品av一区二区| 国产高清视频在线观看网站| 国产成人freesex在线| 国产亚洲av嫩草精品影院| 亚洲自拍偷在线| 国产一区有黄有色的免费视频 | 亚洲精品国产av成人精品| av在线观看视频网站免费| 乱码一卡2卡4卡精品| av在线天堂中文字幕| 在线天堂最新版资源| 日本黄色片子视频| 美女黄网站色视频| 国产探花在线观看一区二区| 成年女人永久免费观看视频| 欧美另类亚洲清纯唯美| 搞女人的毛片| 一级毛片久久久久久久久女| 国产亚洲精品av在线| 欧美性感艳星| 不卡视频在线观看欧美| 亚洲aⅴ乱码一区二区在线播放| 99久久九九国产精品国产免费| 好男人视频免费观看在线| 亚洲国产高清在线一区二区三| 老司机影院成人| 日韩一本色道免费dvd| 国产视频内射| 国产不卡一卡二| 日韩欧美国产在线观看| 中文精品一卡2卡3卡4更新| 99热精品在线国产| 我要看日韩黄色一级片| 国产探花极品一区二区| 性插视频无遮挡在线免费观看| 日本免费在线观看一区| 国产片特级美女逼逼视频| 亚洲精品自拍成人| 久久久国产成人精品二区| 亚洲av成人精品一二三区| 国产又黄又爽又无遮挡在线| АⅤ资源中文在线天堂| 国产精品不卡视频一区二区| 久久6这里有精品| 在线天堂最新版资源| 国产精品一区www在线观看| 2022亚洲国产成人精品| 极品教师在线视频| 高清日韩中文字幕在线| 午夜视频国产福利| 观看美女的网站| h日本视频在线播放| 中文字幕精品亚洲无线码一区| 日日撸夜夜添| 免费看光身美女| 欧美xxxx性猛交bbbb| 我的女老师完整版在线观看| 搞女人的毛片| 国产精品美女特级片免费视频播放器| 国产精品久久久久久精品电影| 伦精品一区二区三区| 亚洲欧美中文字幕日韩二区| 亚洲国产精品专区欧美| 内射极品少妇av片p| 国产精品嫩草影院av在线观看| 久久精品夜夜夜夜夜久久蜜豆| 91精品一卡2卡3卡4卡| 久久久色成人| 搡女人真爽免费视频火全软件| 日韩成人伦理影院| 国内精品宾馆在线| 搡女人真爽免费视频火全软件| 色5月婷婷丁香| 简卡轻食公司| 99久久中文字幕三级久久日本| 成人高潮视频无遮挡免费网站| 中文字幕精品亚洲无线码一区| 亚洲av免费在线观看| 午夜免费男女啪啪视频观看| 精品免费久久久久久久清纯| 人体艺术视频欧美日本| 国产在视频线精品| 亚洲天堂国产精品一区在线| 久久这里只有精品中国| 中文字幕久久专区| 久久综合国产亚洲精品| 亚洲色图av天堂| av女优亚洲男人天堂|