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

    Experimental characteristics of oblique shock train upstream propagation

    2017-11-20 12:07:17WangChengpengXueLongshengTianXuang
    CHINESE JOURNAL OF AERONAUTICS 2017年2期

    Wang Chengpeng,Xue Longsheng,Tian Xuang

    College of Aerospace Engineering,Nanjing University of Aeronautics and Astronautics,Nanjing 210016,China

    Experimental characteristics of oblique shock train upstream propagation

    Wang Chengpeng*,Xue Longsheng,Tian Xuang

    College of Aerospace Engineering,Nanjing University of Aeronautics and Astronautics,Nanjing 210016,China

    Boundary layer separation;Compressible flow;Shock train;Oscillating flow;Wind tunnel testing

    The structure and dynamics of an oblique shock train in a duct model are investigated experimentally in a hypersonic wind tunnel.Measurements of the pressure distribution in front of and across the oblique shock train have been taken and the dynamics of upstream propagation of the oblique shock train have been analyzed from the synchronized schlieren imaging with the dynamic pressure measurements.The formation and propagation of the oblique shock train are initiated by the throttling device at the downstream end of the duct model.Multiple reflected shocks,expansion fans and separated flow bubbles exist in the unthrottled flow,causing three adversepressure-gradient phases and three favorable-pressure-gradient phases upstream the oblique shock train.The leading edge of the oblique shock train propagates upstream,and translates to be asymmetric with the increase of backpressure.The upstream propagation rate of the oblique shock train increases rapidly when the leading edge of the oblique shock train encounters the separation bubble near the shock reflection point and the adverse-pressure-gradient phase,while the oblique shock train slow movement when the leading edge of the oblique shock train is in the favorablepressure-gradient phase for unthrottled flow.The asymmetric flow pattern and oscillatory nature of the oblique shock train are observed throughout the whole upstream propagation process.

    1.Introduction

    The oblique shock train(OST)1flow diffusion phenomenon,which involves an interaction between the duct’s peripheral boundary layer and central oblique-shock-wave field(Fig.1(d)),usually appears in constant or nearly constant cross-sectional area supersonic/hypersonic duct flows.The study of such a complex flow structure in a confined duct under afiniteadversepressuregradienthasimportantimplicationsfor thedesign andoperationofa varietyofdevices includinghypersonic vehicle inlet/isolator,wind tunnel diffusers,supersonic ejectors,etc.To develop the design methods and control strategies of the flow devices,it is necessary to fully understand the mechanismoftheOST.Differentflowconditionswillleadtodifferent kinds of shock train structure:normal shock train occurs for lower incoming Mach numbers of about 1.2–2.2(Fig.1(c)),whereas oblique shock train(Fig.1(d))occurs for higher incoming Mach numbers.Yet the mechanisms of compressionof supersonic gas flow in normal shock train and oblique shock trainflowdiffusionmaybedifferent.Bothofthemcoincidewith different topology of shock boundary layer interaction.Many studies have been conducted on the normal shock train fluid structuretobetterunderstandtheinfluence ofdifferent pressure levels and Mach numbers.2–12Nevertheless,studies on the OST are not very abundant in the literature.13,14

    According to the previous studies,the flow pattern of the shock wave/boundary layer interaction in a duct can be classified into different con figurations(Fig.1).While the freestream Mach number just upstream of the shock is less than about 1.27(Fig.1(a)),15the shock is close to a normal shock and no flow separation occurs at the foot of the shock.For the Mach number between about 1.3 and 1.5(Fig.1(b)),16a nearly normal shock with bifurcated foots is observed and a small boundary layer separation region exists immediately downstream of the bifurcated shock.As the Mach number increases further,the interaction between shock waves and boundary layer is significant,a Mach reflection(MR)is established at the center region of the duct flow(Fig.1(c)),and more than one shock appear in central shock wave field of the duct flow.A series of shocks as illustrated in Fig.1(c)has been called‘normal shock train” in the paper1to indicate such a sequence of shocks.The initial normal shocks bifurcate for the normal shock train case and interact with the separating wall boundary layer at the foot of the bifurcated shocks.The pattern shown in Fig.1(d),named after an ‘oblique shock train”1,tends to appear at higher freestream Mach number than that of a normal shock train and regular reflection(RR)does,and it starts with initial oblique shock waves that separate the boundary layer.The shock-induced separation flow shapes are also different besides the shock pattern between both the shock trains.The restricted shock separation(RSS)17is present in the normal shock train case,and RSS is characterized by a restricted length scale of the separation region in which the mean flow tilts away from the wall before the flow reattaches and continues downstream as an attached boundary layer.Depending on the strength of the secondary shock,more than one separation bubble may be present.For an oblique shock train case,the flow pattern is usually asymmetric and both the RSS and free shock separation(FSS)17exist simultaneously(Fig.1(d)).In FSS,the oblique-shock-induced separation region on the random upper or lower wall fails to reattach on the downstream length of the duct.The separated shear layer forms as it is convected downstream.

    In recent years,the dual-mode scramjet has received more and more research interests and an isolator,also called a supersonic diffuser,is regarded as an essential component of a high-speed dual-mode ramjet engine.The isolator portion is actually a constant-area duct located between the engine inlet and combustor section,and its purpose is to contain the shock train caused by the downstream rise in pressure from combustion.In the scramjet mode of operation,an OST often forms in the isolator for higher freestream Mach number.13,18Wagner et al.studied the dynamics of the unstart shock system in an isolator model in Mach 5 flow by controlling unstart through the deflection of a flap located at the exit plane of the isolator.Unstart was seen to progress upstream in the form of an ‘unstart shock system”.This ‘unstart shock system”took the form of OST but not the form of the normal shock trains as reported by O’Byrne et al.18Wagner et al.used schlieren and high frequency pressure sensors to investigate the OST feature.The data were used to calculate the unstart shock system speed moving upstream and characterize the flow structure.The particle image velocimetry(PIV)measurement was also used to show the velocity field of the started and unstarted OST flow.19The quasi-steady state properties of an OST in a Mach 2.75 ducted flow were studied by Klomparens et al.14including the shock train location,the amplitude of shock displacement and the shape of time-averaged pressure distribution.The fine flow structures of the shock train in an isolator flow were revealed by Chen et al.20,21who performed Nano-tracer Planar Laser Scattering(NPLS),which makes it propitious for studying flow mechanism.Tan et al.22and Li et al.23investigated the transient flow structures of the unstarted flow by downstream choking in generic rectangular hypersonic inlets.For the purpose of preventing combustorinlet interaction,the detection of the location of the isolator shock train is of utmost importance.Le et al.24found that monitoring the standard deviation of wall pressure appears to be the best method.Chang et al.25observed one of two novel oscillatory patterns of inlet buzz:a non-oscillatory violent pattern in which the OST stays at its most upstream position and oscillates just slightly.

    Fig.1 Schematic of shock wave/boundary layer interaction in duct.

    Although substantial work has been performed,the OST dynamics is an area which is still not very well understood.For example,one of key issues is the in fluence of background waves on the OST.Tan et al.26in 2012 thought that the actual flow condition in a scramjet isolator is much more complex since background waves exist,including the impingement shock,separation shock and expansion waves.The interaction of the background waves and the leading shock of the OST would make the pointing direction of the core flow switches upwards and downwards repeatedly during the upstream propagation of the OST.Experimental results of Tian27in 2013 show that steady and rapid forward motion of the OST exists with the increased backpressure.Chang et al.28in 2016 investigated the influence of background waves on the forward motion of the OST using unsteady Reynolds-averaged Navier-Stokes(RANS),and they described the mechanism of motion state of the OST changes from slow forward movement to sharp forward movement and the prediction method of the occurrence of sharp forward movement.

    Another important issue is the asymmetric flow structure of the OST.The unstart OST system is seen to be asymmetric with respect to the duct horizontal centerline in the paper13and Wagner et al.thought that it may be owing to the asymmetric thickness distributions of the inflow boundary layer on the floor and the ceiling.Meanwhile,they found that the propagation of the unstart OST system is related to boundary layer separation in the isolator.However,a full symmetric test model with respect to the horizontal centerline is used in current work,and an asymmetric OST structure is still observed.

    There is no affirmative answer as to how the symmetric level of background waves influences the dynamic characteristics of the OST.In other word,will the OST structure be symmetric and will the OST flip be between the upper wall and lower wall during the upstream propagation process if the background waves are symmetric to the duct horizontal centerline?To answer these questions,a relatively large-scale experimental model(the height and width of the duct are 27.5 mm and 60 mm,respectively)has been designed to research the interaction phenomenon between the OST and the symmetric background waves.Especially,what is the dynamic characteristic of the OST upstream propagation when the leading edge of the OST meets the separation bubbles in the background waves.

    The current paper describes the results of an experimental study of the OST in a special model that was mounted in a hypersonic wind tunnel.Optical visualization of the flow in the duct was performed,and distributions of static pressure in typical longitudinal sections were obtained.Patterns of the OST in the duct were identified.The forward-propagation process of the OST has been analyzed based on the above data.The other purpose of this study is to deepen the understanding of the dynamic characteristics of the OST in the duct and provide a fundamental research for the engineering applications.

    2.Experimental setup

    Fig.2 Hypersonic Wind Tunnel of Nanjing University of Aeronautics and Astronautics(NHW).

    A schematic of the measurement system setup is shown in Fig.4.The time-accurate pressure measurement system consisted of 10 fast-response transducers(CYG 1053AT,Double Bridge Inc.)with respective amplifiers and data acquisition cards(National Instruments PXIe 6358).During the wind tests,the data sampling frequency was 10 kHz per channel and the amplified signals of ten transducers were digitized at a rate of 20 kHz with two A/D cards mentioned above mounted in a personal computer in which there is a controlling program by LabVIEW code.The pressure sampling process lasts for more than 9 s to cover the entire operating period of the NHW wind tunnel.

    To investigate the unsteady behavior of the flow and the process of the OST generation and propagation,the high speed schlieren imaging system was used with a typical Z-style optical apparatus to visualize flow phenomena in an effort to further the understanding of the transient shock/boundary layer interaction.A NAC(NAC Image Technology)Hotshot High Speed Camera with a 1/500000 s electronic shutter was used to capture the schlieren imagery.The sampling rate of the schlieren visualization can reach as high as 5000 frames per second.A horizontal knife edge and a Tungsten Halogen lamp were used for the schlieren system.

    For the flow choking system,two flaps at the exit of the duct section were driven by a servo motor to open and close at the same time.In the experimental process,the controller of the servo motor is selected as the main control hub.The pulse/delay signal emitted by the main controller was used to synchronize the qualitative visual schlieren,quantitative pressure data measurements and motor drive control systems.

    3.Results and analysis

    3.1.Reflected shock structures for undisturbed flow

    Before we discuss the OST structure in the test model under downstream throttled conditions,it is necessary to understand carefully the flow pattern and the corresponding parameter distributions of the reflected shock structures downstream spurred by the multiple shock/boundary layer interactions.Fig.5(a)shows the schlieren image for the reflected shock structures for undisturbed flow.Owing to the size of the schlieren mirrors and the effect of the bracket of the mirrors,the upstream flow field of the test model cannot be visualized fully.Therefore,the special schlieren experiment of the separated flow region 1 induced by the re flected shocks was again conducted after all the runs of the wind tunnel have been finished(Fig.5(b)).Although these two fields of view are not acquiredsimultaneously,on different runs of the wind tunnel,the incoming freestream parameters and the shock structure are almost coincident.

    Table 1 Test conditions.

    Fig.3 Schematic of test model.

    Fig.5(c)gives the corresponding numerical image to understand the flow characteristics further.Here,the full Navier-Stokes equations for two-dimensional flow are numerically solved by the finite-volume method.A second order total variation diminishing method was used to discretize the convection terms.The time terms are discretized with a second order implicit scheme and the SST(Shear Stress Transport)k-x model was used to model the turbulent flow with a compressibility correction for high speed flow.The structured mesh of 45000 elements was used for the numerical simulation andy+at the wall surface is less than 1.

    Fig.4 Schematic of test model in NHW with schlieren imaging and data acquisition(DAQ)system setup.

    Fig.5 Reflected shock structures for undisturbed flow.

    Fig.6 Static pressure distributions of lower wall of duct section and propagation rate of OST.

    Fig.7 Pressure time histories of transducers T1,T6 and T10 under varying backpressure conditions.

    As can be seen in Fig.5(c),two initial oblique shocks from two opposite 7℃ompression ramps intersect and reflect visually in the upstream convergent section and the pair of reflected shocks is captured near the shoulder of the test model in Fig.5(b).The reflected shocks impinge on the lower and upper wall near the shoulder(the joint connection between the convergent section and duct section),which results in an increase in boundary layer thickness and the formation of a separation bubble.The view(Fig.5(a)and(b))was obtained using the horizontal knife edge to reveal density gradients including those due to shocks and boundary layers,so two distinct separation bubbles ‘1” induced by two reflected shocks are seen in Fig.5(b).Two separated shocks ‘A” and ‘B” form and two reattached shocks ‘C” and ‘D” are seen to form downstream of the separated shocksAandB.The expansion fans ‘a(chǎn)” and ‘b” in Fig.5(c)and d originating at the shoulder can also be seen upstream of the separated shocks ‘A” and‘B”.Another pair of expansion fan ‘c” and ‘d” in Fig.5(c)induced by the separation bubble reattachment can be seen between the separated shocks ‘A”, ‘B” and reattached shocks‘C”,‘D”.In Fig.5(b),the expansion fan ‘c” can also be seen as a white-fan region between‘B” and ‘D”.Meanwhile,the shear layer surrounding the separation bubble ‘1”is seen as a white curve like a bridge in Fig.5(b).

    The subsequent reflected shocks are seen as the signs ‘E”and ‘F”in Fig.5(a).Another separation bubble ‘3”is located near the reflection point and the dimension of the separation bubble is smaller than the separation bubble ‘1” due to the decrease of the strength of the impinged shock.The expansion fan ‘e” and ‘f” originating at the peak of the separation bubble ‘3” can also be seen in Fig.5(c).At last,the reflected shocks ‘E”and ‘F” impinged just on the exit of the duct section and the third impingement point is formed.

    Fig.8 Power spectral density contours of pressure signals of T1.

    Fig.6 is the normalized static pressure distributions on the lower wall of the flow path.Good consistency can be seen between the calculated values and the experimental measured values.As a result of the shock/boundary layer/expansion fan interactions,the static pressure distributions along the lower wall of the duct section exhibit three adverse-pressuregradient phases(apg 1,apg 2,apg 3)and three favorablepressure-gradient phases(fpg 1,fpg 2,fpg 3).In the range of separation bubble ‘1”,the values of the pressure of T1 and T2 nearly keep constant.The formation of the ‘a(chǎn)pg 1” and‘a(chǎn)pg 2” phases is mainly due to the reattached shock ‘D”and ‘F”.For the ‘a(chǎn)pg 3” phase,the adverse pressure gradient should be owing to the reflected shock because the CFD data indicate that no separated flow exists in this region.Both phases of ‘fpg 1” and ‘fpg 2” locate in the range between the reattached shocks.The ‘fpg 3” phase is just at the entrance of the downstream expansion flow path(Fig.3).In addition,the rms pressures of T3,T8 and T10 have the higher values than those of the other transducers.T3 and T8 are just downstream of separation bubbles,whereas T10 is located quite near the shock reflection point.The shock wave/turbulent boundary-layer interactions should be responsible for the increased pressure fluctuations.This phenomenon has been found in the experimental results of Wagner et al.13As the OST crossed over T1–T10,a pressure rise was perceived,which made it possible to calculate the propagation rates of the OST.The upstream propagation rates of the OST are also given out in Fig.6(pis pressure,p1is freestream pressure,the right ordinate represents the average velocity of the OST propagating between transducers).The rates are plotted at these positions corresponding to the midpoints between consecutive transducers.The analysis of the propagation process of the OST is presented in Section 3.2.

    Fig.9 Power spectral density contours of pressure signals of T6.

    Fig.10 Power spectral density contours of pressure signals of T10.

    Fig.11 Pressure time histories obtained at T1-T10 for whole process of wind tunnel test.

    These results in Figs.5 and 6 indicate that the flow field on the symmetrical plane of the test model is not influenced by three-dimensional effects and the flow pattern is entirely symmetric about the horizontal centerline of the test model.As a research baseline,it is important for the following section to acquire such information mentioned above.

    Fig.12 Power spectral density contours and phase differences of all pressure signals.

    3.2.Formation and propagation of OST

    As is mentioned above,two flaps linked with gear system can be closed synchronously.The control strategy of the flaps is as follows:during this sequence,the inner surface of two flaps is closed to 9.7?to the horizontal plane in the initial stage of the experiment.At 1.5 s,the servo-motor begins to drive two flaps to close gradually and the flaps stop moving at about 6 s till the run is terminated,where the inner surface of the flap is leaning at a 17.6?angle to the horizontal plane.The drive rod under the flap rises at a constant speed of 2 mm/s during the period of the flaps closing.When the flaps are fully closed,the angle between the inner surface of the lower flap and the horizontal plane is 18.2?27.

    Fig.7 shows the time histories of the wall pressure measured by transducers number T1,T6 and T10 in Fig.5 during a sequence of flaps closing.During this sequence,events such as wind tunnel start,steady state operation,OST propagation,and OST oscillation are clearly observed through static pressure data.The starting shock of the wind tunnel is observed at about 0.35 s.The flaps begin in the open position and are then raised at 1.5 s to initiate downstream flow choking.At about 3.5 s,the leading edge of the OST firstly arrives at the position of the transducer T10 located most downstream of the duct section,and the static pressure distributions at T10 exhibit a sharp increase trend at this moment.At about 6 s,two flaps at the duct section exit stop moving and then keep immovable till 8 s.From 6 to 8 s,the stable pressure oscillations are seen at three transducers T1,T6 and T10.Finally,the flow in the wind tunnel is cut off at 8 s.

    The Power Spectrum Density(PSD)contours of T1,T6 and T10 pressure signals are shown in Figs.8–10,and the left figures are calculated by Wavelet-Transform(WT)during 0–8 s to illustrate the time-frequency characteristics of OST and the right ones are calculated by Fast Fourier-Transform(FFT)aroundt=7 s.Att=5.6 s,a sudden sharp pressure increase occurs when the leading edge of the OST sweeps the transducer T1.As seen in Figs.8–10,a dominant frequency is about 578 Hz att=6–8 s for three transducers T1,T6,and T10.During the upstream propagation of the OST from 4.5–6 s,the trend of the pressure distributions at T6 is monotonic increase,and the centralized frequency cannot be observed during this period(Fig.9).

    Furthermore,more attention is paid to the static pressure gradient at T1–T10.The static pressure-time histories(filtered by cutoff frequency 1 000 Hz)obtained at T1–T10 for the whole test process during the flaps closing are plotted in Fig.11 to illustrate the characteristics of the OST propagation.The right plot of Fig.11 shows the pressure time histories with a smaller time window fort=3–6 s than that used in the left plot,which emphasizes the propagation of the OST.Fig.12(a)-(g)show the PSD contours of the pressure signals for T2–T5 and T7–T9 during 0–8 s to illustrate the time-frequency characteristics.

    Fig.13 Typical schlieren images and pressure-time histories at t=3.3–4.3 s.

    During the propagation process,an increase in throttling degree is expected to drive the OST leading edge upstream through the duct section.This is reflected by the results of Figs.11 and 12.As the throttling degree increases,each station of high-frequency transducers T1-T10 registers the arrival of the OST as an increase in pressure when the leading edge shock moves upstream.Initially,the steady values ofp/p1are measured by ten transducers T1-T10 fort<3.5 s.A sharp pressure rise is seen to occur at T10 as mentioned above fort=3.5 s,which indicates that the leading edge of the OST has influenced the most downstream transducer in the duct section.Meanwhile,a small peak value in pressure at T9 has been found,but the substantial pressure rise of transducer T9 occurs att=3.65 s.Maybe it is not very suitable for monitoring the standard deviation of instantaneous pressure to determine the location of the OST leading edge.Att=4.2 s,two transducers at station T8 and T7 register the arrival of the OST as the pressure increases sharply at approximately the same time and the PSD in Fig.12(e)and(f)has almost the same distribution.As throttling degree increases further,a gradual rise in pressure is sensed by T6 att=4.2–5.6 s.This indicates that the OST moved upstream further and has increased in length.The trend is shown further as throttling degree increases and the same gradual rise in pressure is seen by T5 att=5–5.6 s.At aboutt=5.7 s,a sharp increase in the pressure rise behavior is simultaneously observed at T1-T4.Att=5.7–6 s,the pressure at T1-10 rises gradually.At last,the pressure distribution at T1-T10 switches to the oscillating state with a stabilized averaged value at 6–8 s(Figs.8–10,and 12).

    From Figs.11 and 12,it is shown that the pressure rise generated by the OST exhibits complicated distribution while the OST shifts upstream in the duct section,and even the pressure at T9 decreases at 4.2–5 s.The oscillations with different levels in amplitude and frequency can be observed on the pressure distributions and PSD contours at all the transducers and all the period oft=3.5–8 s.This indicates the complicated process of the OST propagation in the duct section.Therefore,the schlieren image associated with the transient pressure data is presented and analyzed carefully as stated below.

    Fig.14 Typical schlieren images and pressure-time histories at t=4.1–5.8 s.

    (1)Att=3.3–4.3 s,the onset of the OST was defined by using both the schlieren imaging and the pressure time history of T10 to aid in the characterization of the OST dynamics at the rear of the duct section.The initiation of the OST under rising throttling conditions is often considered to be associated with boundary layer separation near the duct exit.13The most upstream pair of separation oblique shocks(GandHin Fig.13(a))crosses to form the leading ‘X” structure of the OST.Att=3.636 s,these oblique shocksGandHappear in the downstream view of the duct section.Both shocksGandHinteract at their impingement location on the lower and upper wall,and the schlieren image and pressure data in Fig.13(c)indicate that the impingement point of the shockHis located in the middle of T9 and T10 att=3.636 s.Att=3.676 s,the shock H just sweeps upstream the T9 transducer and the pressure of T9 takes on a sharply increasing trend from 3.645 s to 3.657 s.Att=4.176 s,the separation shockGand the reflection shockEimpinge and interact together near the upper wall.The OST gradually moves upstream,a pair of separation shocksGandHare slightly asymmetric during the period oft=3.636–4.200 s,and the OST forward propagation process is a stable motion.The observed propagation phenomenon of the OST is related to the pressure distribution of ‘fpg 2” phase(Fig.6)in this region for undisturbed flow.

    (2)The OST shifts upstream continuously and the leading separated shockHarrives at T8 att=4.19 s(Fig.14(c)).T8 is close to the impingement point of the reflection shockCwhere a separation flow region 3 shown in Fig.5 forms.Then,the schlieren results show that the leading edge of the OST moves sharply across the T7 and T8 att=4.190–4.242 s.As a result,the shock H sweeps quickly over the measuring points T8 and T7 in turn and the pressure of T8 and T7 almost rises sharply at one time(Fig.14(c)).Meanwhile,the pressure at the measuring point of T9 also increases sharply att=4.198 s,but a gorge att=4.198–4.270 s is found in the pressure distribution of T10.As one can see,the core flow of the OST turns toward the lower wall of the duct section and an asymmetric flow structure of the OST appears as the schlieren images in Fig.14(a).

    During the period oft=4.242–5.504 s,the OST propagation reenters a stable motion process like that int=3.636–4.190 s.The OST moves short distances in this period and the shocks in the OST shrink to the lower wall.As a possible result,the small separation region exists near the lower wall,while the larger separation region forms near the upper wall.The value of T9 almost keeps constant in this period because the downstream flow structure isolates the effects of downstream backpressure.The propagation of downstream back-pressure increase occurs mainly in the flow region near the upper wall,so the increasing backpressure squeezes the shocks in the OST onto the lower wall.

    Fig.15 Typical schlieren images and pressure-time histories at t=5.5–6.0 s.

    During the period oft=4.40–5.65 s,the pressure of T6 and T5 rises gradually and this trend of pressure signals is distinctly different from that of T8 and T7.The corresponding movement of the OST is relatively slow.As shown in Fig.6,the pressure distribution in front of the leading separated shock H is located in the ‘fpg 1”phase for undisturbed flow.During the periods of the forward movements of the OST from T6 to T4,the pressure in front of the leading separated shockHincreases,so it is hard to separate the boundary layer flow near the wall and the OST moves forward slowly.

    The above experimental results show that the OST sharp movement occurs when the leading edge of the OST is coming to the shock reflection point where a separation bubble exists,while the OST slow movement occurs when the leading edge of the OST is in the ‘fpg” phase for undisturbed flow.It appears that both the separated bubble near the shock reflection point and the ‘a(chǎn)pg” phase for the flow undisturbed by the OST will be jointly responsible for the sharp forward movement of the OST.

    (3)As shown in Fig.15,the lower branchHof the leading shock of the OST comes to the measuring point of T4 att=5.662 s.Then,the leading edge of the OST enters the region of flow separation 1(Fig.5(b))att=5.676 s.As a result,an obvious sharp movement of the OST is obtained,and the leading edge of the OST moves rapidly across the T1–T4 where the local detached boundary layer and the impingement point exist.The corresponding pressure of T1-T4 rises sharply att=5.662–5.686 s.Meanwhile,the pressure at the measure point of T5 rises,whereas the pressure of T6 and T7 decreases at that time.

    Att=5.694 s,the OST structure suddenly retreats to the downstream station in the middle of T3 and T4(Fig.15(a)),and the core flow of the OST is deflected upwards abruptly.The corresponding pressure at T1-T3 decreases dramatically to the value of pressure where the flow is not disturbed by the OST(Fig.15(c)).Thereafter,the OST moves upstream quickly again att=5.7 s.Att=5.9 s,the downstream throttling system stops working.There are only two shear layers observed in optical window(Fig.15(a)).

    This sharp forward movement is similar to the above sharp forward movement near the transducers T7 and T8(Fig.14),but the leading edge shockHof the OST sweeps over four transducers T1–T4 at this time.It is believed that the combined effect of the large scale separation bubble 1(Fig.5(b))and the steeper ‘a(chǎn)pg 1” phase(Fig.6)results in the transient movement of the OST.

    Fig.16 Typical schlieren images and pressure-time histories at t=6–8 s.

    (4)Att=6–8 s,the pressure distribution at T1-T10 keeps stable(Fig.16(c))and the schlieren images are similar to the image att=5.952 s(Fig.16(a)).In these schlieren images in Fig.16(a),there exists no distinct shock,and the pressure increases gradually(Fig.16(b)).This region behind the shock train is called ‘mixing region”in Ref.16where the mixing of non-uniform flow gives rise to pressure rise.

    4.Conclusions

    The dynamic characteristics of the oblique shock trains in a duct model are investigated experimentally in a hypersonic wind tunnel.The flow patterns and transient pressure distributions in the undisturbed and oblique shock train flow are obtained under varying throttling conditions with the synchronized schlieren imaging and dynamic pressure measurements.The main results are summarized as follows:

    (1)When no backpressure is generated from the throttling device,an unthrottled flow without the oblique shock train forms in the ducted model.It is found that multiple reflected shocks,expansion fans and separated flow bubbles exist in the unthrottled flow,causing three adversepressure-gradient phases and three favorable-pressuregradient phases upstream the oblique shock train.The experimental and numerical results indicate that the unthrottled flow on the symmetrical plane of the test model is a two-dimensional flow and the flow pattern is entirely symmetric about the horizontal centerline of the test model.

    (2)As the throttling degree increases,oblique shock train is observed to propagate upstream.The upstream propagation rate of the oblique shock train increases rapidly when the leading edge of the oblique shock train encounters the separation bubble near the shock reflection point and the adverse-pressure-gradient phase mentioned above,while the oblique shock train slow movement when the leading edge of the oblique shock train is in the favorable-pressure-gradient phase for undisturbed flow.The asymmetric flow pattern and oscillatory back and forth nature of the oblique shock train are observed throughout the whole upstream propagation process.

    Acknowledgements

    This study was supported by the National Natural Science Foundation of China(Nos.51476076 and 10702029).

    1.Heiser WH, Pratt DT.Hypersonicairbreathingpropulsion.Reston:AIAA Inc.;1994,p.251–256.

    2.Tamaki T,Tomita Y,Yamane R.A study of pseudo-shock,1st Report,k-type pseudo-shock.Bull JSME1970;13(55):51–8.

    3.Carroll BF,Dutton JC.Turbulence phenomena in a multiple normal shock wave/turbulent boundary-layer interaction.AIAA J1992;30(1):43–8.

    4.Carroll BF,Dutton JC.Multiple normal shock wave/turbulent boundary-layer interactions.J Propul Power1992;8(2):441–8.

    5.Wang CP,Tian XA,Cheng KM,Wu YZ.Numerical analysis of pseudo-shock flow diffusion phenomenon in variable cross-section duct.Proc Inst Mech Eng Part G:J Aerosp Eng2008;222(8):1109–21.

    6.Weiss A,Grzona A,Oliver H.Behavior of shock trains in a diverging duct.Exp Fluids2010;49(2):355–65.

    7.Weiss A,Oliver H.Behaviour of a shock train under the influence of boundary-layer suction by a normal slot.Exp Fluids2012;52(2):273–87.

    8.Gawehn T,Gu¨lhan A,Al-Hasan NS,Schnerr GH.Experimental and numerical analysis of the structure of pseudo-shock systems in laval nozzles with parallel side walls.Shock Waves2010;20(4):297–306.

    9.Lee HJ,Lee BJ,Kim SD,Jeung IS.Investigation of the pseudoshock wave in a two-dimensional supersonic inlet.J Visual2010;13(1):25–32.

    10.Sun LQ,Sugiyama H,Mizobata K,Fukuda K.Numerical and experimental investigations on the Mach 2 Pseudo-shock wave in a square duct.J Visual2003;6(4):363–70.

    11.Koike S,Osawa J,Nakakita K,Kato H,Kameda M.Simultaneous pressure and velocity field measurement of pseudo-shockwave using PSP and PIV.Reston:AIAA;2012.Report No.:AIAA-2012-1190.

    12.Oka T,Ono D,Miyazato Y.Study of shock trains and pseudoshock waves in constant area ducts.Reston:AIAA;2014.Report No.:AIAA-2014-0949.

    13.Wagner JL,Yuceil KB,Valdivia A,Clemens NT,Dolling DS.Experimental investigation of unstart in an inlet/isolator model in Mach 5 flow.AIAA J2009;47(6):1528–42.

    14.Klomparens RL,Driscoll JF,Gamba M.Unsteadiness characteristics and pressure distribution of an oblique shock train.Reston:AIAA;2015.Report No.:AIAA-2015-1519.

    15.Bogar TJ,Sajben M,Kroutil JC.Characteristic frequencies of transonic diffuser flow oscillations.AIAA J1983;21(9):1232–40.

    16.Matsuo K,Miyazato Y,Kim HD.Shock train and pseudo-shock phenomena in internal gas flows.Prog Aerosp Sci1999;35(1):33–100.

    17.Olson BJ,Lele SK.A mechanism for unsteady separation in overexpanded nozzle flow.Phys Fluids2013;25(11):080901–80924.

    18.O’Byrne S,Doolan M,Olsen SR,Houwing AFP.Analysis of transient thermal choking processes in a model scramjet engine.J Propul Power2000;16(5):808–14.

    19.Wagner JL,Yuceil KB,Clemens NT.PIV measurements of unstart of an inlet-isolator model in a Mach 5 flow.Reston:AIAA;2009.Report No.:AIAA-2009-4209.

    20.Chen Z,Yi SH,Zhu YZ,Wu Y,Zhang QH,Quan PC.Investigation on flows in a supersonic isolator with an adjustable cowl convergence angle.Exp Thermal Fluid Sci2014;52(1):182–90.

    21.Chen Z,Yi SH,Ahmed NA,Kong XP,Zhang QH,Quan PC.Transient behavior of shock train with and without controlling.Exp Thermal Fluid Sci2015;66:79–96.

    22.Tan HJ,Li LG,Wen YF,Zhang QF.Experimental investigation of the unstart process of a generic hypersonic inlet.AIAA J2011;49(2):279–88.

    23.Li ZF,Gao WZ,Jiang HL,Yang JM.Unsteady behaviors of a hypersonic inlet caused by throttling in shock tunnel.AIAA J2013;51(10):2485–92.

    24.Le DB,Goyne CP,Krauss RH.Shock train leading-edge detection in a dualmode scramjet.J Propul Power2008;24(5):1035–41.

    25.Chang JT,Wang L,Bao W,Qin J,Niu J,Xue W.Novel oscillatory patterns of hypersonic inlet buzz.J Propul Power2012;28(6):1214–21.

    26.Tan HJ,Sun S,Huang HX.Behavior of shock trains in a hypersonic inlet/isolator model with complex background waves.Exp Fluids2012;53(6):1647–61.

    27.Tian XA.Inlet design and flowfield analysis of integrated airbreathing hypersonic vehicle[dissertation].Nanjing:Nanjing University of Aeronautics and Astronautics;2013[Chinese].

    28.Xu KJ,Chang JT,Zhou WX,Yu DR.Mechanism and prediction for occurrence of shock-train sharp forward movement.AIAA J2016;54(4):1403–12.

    29.Wang CP,Cheng KM,Wang JF,Wu YZ.Flow field calibrations and force tests of an AGARD model in a Hypersonic Wind tunnel of Nanjing University of Aeronautics and Astronautics(NHW).Proceedings of China 1st modern aerodynamics&aerothermodynamics conference;2006 August 23–26;Mianyang,China.Beijing:National defense industry press;2006.p.356–61[Chinese].

    30.Wang CP,Tian XA,Yan LF,Xue LS,Cheng KM.Preliminary integrated design of hypersonic vehicle configurations including inward-turning inlets.J Aerosp Eng2015;28(6),040141431-11.

    3 May 2016;revised 27 October 2016;accepted 1 December 2016

    Available online 14 February 2017

    *Corresponding author.

    E-mail address:wangcp@nuaa.edu.cn(C.Wang).

    Peer review under responsibility of Editorial Committee of CJA.

    性高湖久久久久久久久免费观看| 日产精品乱码卡一卡2卡三| 免费播放大片免费观看视频在线观看| 综合色丁香网| 婷婷色av中文字幕| 亚洲男人天堂网一区| 精品人妻偷拍中文字幕| 最近最新中文字幕免费大全7| 午夜免费鲁丝| 超碰97精品在线观看| 少妇被粗大的猛进出69影院| 最新的欧美精品一区二区| 国产男女内射视频| 曰老女人黄片| 人妻 亚洲 视频| 最近中文字幕2019免费版| 国产精品无大码| 久久久精品94久久精品| 午夜激情av网站| 国产色婷婷99| 边亲边吃奶的免费视频| 女人精品久久久久毛片| 欧美 亚洲 国产 日韩一| 久久精品国产a三级三级三级| 老汉色av国产亚洲站长工具| av福利片在线| 在线观看免费高清a一片| 丝袜在线中文字幕| 国产免费现黄频在线看| 久久久久久伊人网av| 在线观看免费高清a一片| 在线亚洲精品国产二区图片欧美| 成人亚洲精品一区在线观看| 一级毛片 在线播放| 亚洲成人手机| 日本91视频免费播放| 亚洲人成77777在线视频| 26uuu在线亚洲综合色| 观看美女的网站| 国产乱来视频区| 久久毛片免费看一区二区三区| 欧美精品一区二区大全| 色哟哟·www| 男人爽女人下面视频在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 高清黄色对白视频在线免费看| 爱豆传媒免费全集在线观看| a 毛片基地| 男人添女人高潮全过程视频| 亚洲精品视频女| 久久久精品94久久精品| 成人二区视频| 色94色欧美一区二区| 又大又黄又爽视频免费| 最近手机中文字幕大全| 精品视频人人做人人爽| 欧美亚洲 丝袜 人妻 在线| av一本久久久久| 亚洲五月色婷婷综合| 亚洲欧美成人精品一区二区| 久久99蜜桃精品久久| 国产高清国产精品国产三级| 精品国产超薄肉色丝袜足j| 亚洲国产精品999| 亚洲精品久久久久久婷婷小说| 欧美国产精品一级二级三级| 三级国产精品片| 亚洲精品一二三| 麻豆av在线久日| 有码 亚洲区| 亚洲精品日本国产第一区| 国产精品麻豆人妻色哟哟久久| 欧美日韩综合久久久久久| 日本wwww免费看| 久久 成人 亚洲| 中文字幕亚洲精品专区| 国产成人精品在线电影| 免费不卡的大黄色大毛片视频在线观看| av卡一久久| 999久久久国产精品视频| 亚洲精品美女久久av网站| 亚洲熟女精品中文字幕| 色播在线永久视频| 欧美av亚洲av综合av国产av | 日韩av不卡免费在线播放| 亚洲伊人久久精品综合| 国产精品亚洲av一区麻豆 | 欧美人与性动交α欧美软件| 亚洲av免费高清在线观看| videossex国产| 午夜91福利影院| 国产在视频线精品| 高清欧美精品videossex| 久久久久精品性色| 久久久久网色| 狂野欧美激情性bbbbbb| 国产无遮挡羞羞视频在线观看| 欧美精品高潮呻吟av久久| 午夜影院在线不卡| 午夜久久久在线观看| 国产男人的电影天堂91| 777久久人妻少妇嫩草av网站| 久久久精品免费免费高清| 永久免费av网站大全| 国产激情久久老熟女| 一级a爱视频在线免费观看| 欧美变态另类bdsm刘玥| 18在线观看网站| 久久久久久久久免费视频了| 日本91视频免费播放| 免费观看av网站的网址| 国产激情久久老熟女| 黄色视频在线播放观看不卡| 99国产精品免费福利视频| 女性被躁到高潮视频| 男女无遮挡免费网站观看| 婷婷色麻豆天堂久久| 亚洲欧洲精品一区二区精品久久久 | 免费播放大片免费观看视频在线观看| 啦啦啦视频在线资源免费观看| 晚上一个人看的免费电影| 久久女婷五月综合色啪小说| 午夜久久久在线观看| 韩国高清视频一区二区三区| 免费高清在线观看日韩| 丰满乱子伦码专区| 老汉色∧v一级毛片| 美女高潮到喷水免费观看| 少妇人妻久久综合中文| 国产精品久久久久成人av| 日韩精品免费视频一区二区三区| 国产熟女午夜一区二区三区| 2021少妇久久久久久久久久久| 精品视频人人做人人爽| 精品国产乱码久久久久久小说| 性色avwww在线观看| 天天影视国产精品| 久久99精品国语久久久| 婷婷色综合www| 午夜免费观看性视频| 大香蕉久久成人网| 亚洲欧美清纯卡通| 国产毛片在线视频| 少妇的丰满在线观看| 久久久久国产精品人妻一区二区| 亚洲精品美女久久av网站| 精品国产一区二区久久| 亚洲一区中文字幕在线| 国产精品不卡视频一区二区| 久久久精品国产亚洲av高清涩受| 女性被躁到高潮视频| 宅男免费午夜| 色吧在线观看| 在线观看www视频免费| 国产精品.久久久| 国产成人免费观看mmmm| 亚洲精品第二区| 欧美日本中文国产一区发布| 18禁裸乳无遮挡动漫免费视频| 少妇 在线观看| 亚洲在久久综合| 七月丁香在线播放| 下体分泌物呈黄色| 亚洲av电影在线进入| 国产亚洲欧美精品永久| 狠狠精品人妻久久久久久综合| 在线天堂最新版资源| 国产片内射在线| 中文字幕制服av| av.在线天堂| 制服人妻中文乱码| 午夜福利视频在线观看免费| 日产精品乱码卡一卡2卡三| 国产成人免费无遮挡视频| 欧美av亚洲av综合av国产av | 日本猛色少妇xxxxx猛交久久| 少妇熟女欧美另类| 免费在线观看完整版高清| 少妇被粗大猛烈的视频| 在现免费观看毛片| 国产白丝娇喘喷水9色精品| 亚洲综合色网址| 黄网站色视频无遮挡免费观看| videosex国产| 一级毛片电影观看| 亚洲精品一区蜜桃| 丝袜人妻中文字幕| 欧美人与性动交α欧美软件| 熟女av电影| 亚洲国产av影院在线观看| 又黄又粗又硬又大视频| 青青草视频在线视频观看| 中文字幕亚洲精品专区| 婷婷色麻豆天堂久久| 丝袜美足系列| 亚洲欧美成人精品一区二区| av免费在线看不卡| 亚洲精品在线美女| 边亲边吃奶的免费视频| 精品少妇内射三级| 黄片小视频在线播放| 男人操女人黄网站| 不卡av一区二区三区| 久久av网站| 亚洲综合色惰| 久久久久久免费高清国产稀缺| 成人免费观看视频高清| 亚洲av电影在线观看一区二区三区| 18禁观看日本| 最新的欧美精品一区二区| av国产久精品久网站免费入址| 一区二区日韩欧美中文字幕| 精品福利永久在线观看| 国产日韩一区二区三区精品不卡| 欧美成人午夜免费资源| 成人亚洲精品一区在线观看| 爱豆传媒免费全集在线观看| 婷婷成人精品国产| 午夜福利网站1000一区二区三区| 99九九在线精品视频| 热99久久久久精品小说推荐| 日韩三级伦理在线观看| 亚洲欧美清纯卡通| 97人妻天天添夜夜摸| 欧美另类一区| 啦啦啦在线免费观看视频4| 女人高潮潮喷娇喘18禁视频| 街头女战士在线观看网站| 久久久久人妻精品一区果冻| 中文字幕色久视频| 亚洲一码二码三码区别大吗| 欧美人与性动交α欧美精品济南到 | 亚洲,一卡二卡三卡| 美女午夜性视频免费| av免费观看日本| 9热在线视频观看99| 国产精品.久久久| 国产亚洲最大av| 97在线人人人人妻| 久久久久精品久久久久真实原创| 啦啦啦啦在线视频资源| 91午夜精品亚洲一区二区三区| 99久久综合免费| 国产精品国产av在线观看| 99久久中文字幕三级久久日本| 亚洲精品在线美女| 丝袜人妻中文字幕| 国产黄色视频一区二区在线观看| 久久99一区二区三区| 久久国产亚洲av麻豆专区| 国产又色又爽无遮挡免| 自拍欧美九色日韩亚洲蝌蚪91| 欧美成人精品欧美一级黄| 伊人亚洲综合成人网| 精品卡一卡二卡四卡免费| 宅男免费午夜| 男女国产视频网站| 美女中出高潮动态图| 精品久久久久久电影网| 女的被弄到高潮叫床怎么办| 午夜久久久在线观看| 国产乱人偷精品视频| 国产一区二区三区av在线| 亚洲欧美中文字幕日韩二区| 亚洲成人av在线免费| 日本爱情动作片www.在线观看| 日韩精品免费视频一区二区三区| av又黄又爽大尺度在线免费看| 中文字幕制服av| 国产极品天堂在线| 亚洲av免费高清在线观看| 亚洲,欧美,日韩| 男女边吃奶边做爰视频| 亚洲国产看品久久| 国产男人的电影天堂91| 亚洲精品一区蜜桃| 国产熟女午夜一区二区三区| 秋霞在线观看毛片| 亚洲av中文av极速乱| 又黄又粗又硬又大视频| 国产有黄有色有爽视频| 亚洲成人手机| 久久久国产精品麻豆| 亚洲国产精品一区三区| 一区二区日韩欧美中文字幕| 一本色道久久久久久精品综合| 女人被躁到高潮嗷嗷叫费观| 一区在线观看完整版| 精品少妇一区二区三区视频日本电影 | 成年女人在线观看亚洲视频| 日韩三级伦理在线观看| 精品人妻一区二区三区麻豆| 国产一区二区在线观看av| 国产成人欧美| 香蕉国产在线看| 最近的中文字幕免费完整| 日本-黄色视频高清免费观看| 国产1区2区3区精品| 亚洲欧美成人精品一区二区| 搡老乐熟女国产| 欧美在线黄色| 成年女人毛片免费观看观看9 | 2021少妇久久久久久久久久久| 看十八女毛片水多多多| 观看美女的网站| 午夜免费男女啪啪视频观看| 自拍欧美九色日韩亚洲蝌蚪91| 新久久久久国产一级毛片| 欧美bdsm另类| 久久久久久人人人人人| 日本午夜av视频| 少妇人妻精品综合一区二区| 91午夜精品亚洲一区二区三区| 在线 av 中文字幕| 国产免费福利视频在线观看| 免费在线观看视频国产中文字幕亚洲 | 我的亚洲天堂| 91精品三级在线观看| 亚洲欧洲精品一区二区精品久久久 | av线在线观看网站| 国产国语露脸激情在线看| 超碰97精品在线观看| 汤姆久久久久久久影院中文字幕| 少妇人妻精品综合一区二区| 这个男人来自地球电影免费观看 | 91成人精品电影| 日韩精品有码人妻一区| 亚洲精品中文字幕在线视频| 精品人妻偷拍中文字幕| 亚洲欧美中文字幕日韩二区| 精品人妻偷拍中文字幕| 在现免费观看毛片| 国产女主播在线喷水免费视频网站| 久久韩国三级中文字幕| 国产精品二区激情视频| 国产日韩一区二区三区精品不卡| 免费在线观看完整版高清| 国产人伦9x9x在线观看 | 日本av手机在线免费观看| 久久久欧美国产精品| 黑人猛操日本美女一级片| 精品一品国产午夜福利视频| 男女边吃奶边做爰视频| 成年女人在线观看亚洲视频| 一级毛片电影观看| 日韩在线高清观看一区二区三区| 国产精品欧美亚洲77777| 精品一区二区三卡| 一边亲一边摸免费视频| h视频一区二区三区| 久久久久国产一级毛片高清牌| 免费黄频网站在线观看国产| 自拍欧美九色日韩亚洲蝌蚪91| 国产精品久久久av美女十八| 中文精品一卡2卡3卡4更新| 天天影视国产精品| 又粗又硬又长又爽又黄的视频| 丝袜喷水一区| 大片电影免费在线观看免费| 一区福利在线观看| 国产人伦9x9x在线观看 | 亚洲成人手机| 叶爱在线成人免费视频播放| 久久国内精品自在自线图片| 大话2 男鬼变身卡| www日本在线高清视频| 免费观看在线日韩| 女性被躁到高潮视频| 大话2 男鬼变身卡| 九色亚洲精品在线播放| 汤姆久久久久久久影院中文字幕| 侵犯人妻中文字幕一二三四区| 一级片免费观看大全| av女优亚洲男人天堂| 亚洲精品aⅴ在线观看| 亚洲欧美中文字幕日韩二区| 国产精品 国内视频| 国精品久久久久久国模美| 一区福利在线观看| 久久久精品国产亚洲av高清涩受| 激情视频va一区二区三区| 美女大奶头黄色视频| 亚洲国产最新在线播放| 精品国产一区二区三区久久久樱花| 午夜日韩欧美国产| 深夜精品福利| 中文精品一卡2卡3卡4更新| 中文字幕亚洲精品专区| 亚洲av国产av综合av卡| 亚洲,一卡二卡三卡| 国产av一区二区精品久久| 丝袜人妻中文字幕| 亚洲第一区二区三区不卡| 国产亚洲最大av| 999久久久国产精品视频| 看十八女毛片水多多多| 精品久久蜜臀av无| 如何舔出高潮| 亚洲av欧美aⅴ国产| 99九九在线精品视频| 99re6热这里在线精品视频| 夫妻午夜视频| 日韩制服丝袜自拍偷拍| 久久精品亚洲av国产电影网| 亚洲一区中文字幕在线| 90打野战视频偷拍视频| 美女xxoo啪啪120秒动态图| 欧美97在线视频| 久久这里只有精品19| 亚洲成人一二三区av| 成人毛片a级毛片在线播放| 精品少妇内射三级| 久久精品国产亚洲av涩爱| 丰满乱子伦码专区| 欧美在线黄色| 国产成人91sexporn| 18禁动态无遮挡网站| 亚洲av国产av综合av卡| 精品少妇内射三级| 我要看黄色一级片免费的| 秋霞伦理黄片| 爱豆传媒免费全集在线观看| 欧美 日韩 精品 国产| 亚洲天堂av无毛| 一二三四在线观看免费中文在| 午夜福利在线观看免费完整高清在| 亚洲精品av麻豆狂野| 宅男免费午夜| 欧美精品高潮呻吟av久久| 亚洲精品日韩在线中文字幕| 99久久精品国产国产毛片| 欧美老熟妇乱子伦牲交| 欧美精品国产亚洲| 少妇的逼水好多| 午夜免费鲁丝| 美女脱内裤让男人舔精品视频| 色婷婷久久久亚洲欧美| 婷婷色麻豆天堂久久| 美女视频免费永久观看网站| 亚洲av成人精品一二三区| av在线老鸭窝| 欧美精品av麻豆av| 五月开心婷婷网| 美女福利国产在线| 国产有黄有色有爽视频| 一区二区三区乱码不卡18| 亚洲人成77777在线视频| 免费日韩欧美在线观看| 亚洲精品,欧美精品| 999久久久国产精品视频| 最近最新中文字幕免费大全7| 午夜福利乱码中文字幕| 久久99热这里只频精品6学生| 久久久国产精品麻豆| 伦理电影免费视频| 国产又色又爽无遮挡免| 久久 成人 亚洲| av.在线天堂| 亚洲欧美精品综合一区二区三区 | av在线播放精品| 美女脱内裤让男人舔精品视频| 色婷婷久久久亚洲欧美| 国产黄频视频在线观看| 精品亚洲乱码少妇综合久久| 夫妻午夜视频| 黄色一级大片看看| 一区二区三区乱码不卡18| 国产熟女欧美一区二区| 亚洲欧美日韩另类电影网站| 中文字幕人妻熟女乱码| 日本欧美视频一区| 国产综合精华液| 欧美激情极品国产一区二区三区| 一级毛片我不卡| 日韩人妻精品一区2区三区| 精品国产超薄肉色丝袜足j| 日韩制服骚丝袜av| 久久精品熟女亚洲av麻豆精品| 久久人人爽av亚洲精品天堂| 美女脱内裤让男人舔精品视频| 九色亚洲精品在线播放| 国产欧美日韩综合在线一区二区| 日本欧美视频一区| 人人妻人人添人人爽欧美一区卜| 中文字幕另类日韩欧美亚洲嫩草| 99热全是精品| 性少妇av在线| 精品人妻在线不人妻| 最近中文字幕高清免费大全6| 精品国产一区二区三区久久久樱花| 国产精品99久久99久久久不卡 | 啦啦啦中文免费视频观看日本| 国产精品av久久久久免费| 欧美日韩精品网址| 26uuu在线亚洲综合色| 巨乳人妻的诱惑在线观看| 精品第一国产精品| 精品国产露脸久久av麻豆| 国产精品99久久99久久久不卡 | 亚洲成人手机| 久久久亚洲精品成人影院| 如何舔出高潮| 大片电影免费在线观看免费| 久久鲁丝午夜福利片| 一本大道久久a久久精品| 超色免费av| 女人高潮潮喷娇喘18禁视频| 搡老乐熟女国产| 老司机影院毛片| 久久久久国产精品人妻一区二区| 国产精品秋霞免费鲁丝片| 精品人妻在线不人妻| 国产成人免费无遮挡视频| 交换朋友夫妻互换小说| 午夜精品国产一区二区电影| 99久久中文字幕三级久久日本| 亚洲精品一二三| 久久精品国产综合久久久| 1024视频免费在线观看| 精品亚洲乱码少妇综合久久| av天堂久久9| 99精国产麻豆久久婷婷| 青春草国产在线视频| 精品一区二区免费观看| kizo精华| 亚洲精品国产av成人精品| 免费高清在线观看视频在线观看| 女人高潮潮喷娇喘18禁视频| 国产视频首页在线观看| 亚洲成色77777| 午夜激情av网站| 免费播放大片免费观看视频在线观看| 午夜免费观看性视频| 国产成人a∨麻豆精品| 男女啪啪激烈高潮av片| 日韩制服丝袜自拍偷拍| 国产人伦9x9x在线观看 | 不卡av一区二区三区| 久久99蜜桃精品久久| 男人舔女人的私密视频| 一边亲一边摸免费视频| 国产欧美亚洲国产| 777久久人妻少妇嫩草av网站| 亚洲欧洲国产日韩| 国产又色又爽无遮挡免| av电影中文网址| av在线观看视频网站免费| 女的被弄到高潮叫床怎么办| 亚洲欧洲日产国产| 看免费成人av毛片| 高清欧美精品videossex| 男人操女人黄网站| 日韩精品免费视频一区二区三区| 91精品伊人久久大香线蕉| 视频在线观看一区二区三区| 国产野战对白在线观看| 观看美女的网站| 国产女主播在线喷水免费视频网站| 狠狠精品人妻久久久久久综合| 免费av中文字幕在线| 久久久久久久久久人人人人人人| 一级黄片播放器| 少妇被粗大猛烈的视频| 午夜福利视频在线观看免费| 色哟哟·www| 亚洲av电影在线进入| 麻豆精品久久久久久蜜桃| 国产 一区精品| 伊人久久国产一区二区| 亚洲内射少妇av| 亚洲av电影在线进入| 丰满饥渴人妻一区二区三| 久久久久视频综合| 国产在线免费精品| 亚洲五月色婷婷综合| 免费不卡的大黄色大毛片视频在线观看| 十八禁网站网址无遮挡| 欧美精品亚洲一区二区| 一级毛片 在线播放| 精品少妇黑人巨大在线播放| 欧美日韩精品网址| 国产成人91sexporn| 在线精品无人区一区二区三| 国产精品久久久久久精品电影小说| 精品一区二区三区四区五区乱码 | 国产精品久久久久久久久免| 夫妻午夜视频| 精品国产乱码久久久久久男人| 飞空精品影院首页| 丰满乱子伦码专区| 中文字幕人妻丝袜制服| 亚洲精品成人av观看孕妇| 美女主播在线视频| av在线app专区| av网站免费在线观看视频| 9色porny在线观看| 国产探花极品一区二区| 啦啦啦视频在线资源免费观看| 最近中文字幕2019免费版| 在现免费观看毛片| freevideosex欧美| 国产精品蜜桃在线观看| 一级毛片 在线播放| 亚洲人成网站在线观看播放| 91在线精品国自产拍蜜月| 波野结衣二区三区在线| 色婷婷久久久亚洲欧美| 只有这里有精品99| 日韩制服丝袜自拍偷拍| 最近最新中文字幕免费大全7| 久久久精品免费免费高清| 91在线精品国自产拍蜜月| 国产有黄有色有爽视频| 亚洲精品第二区| 精品国产乱码久久久久久小说|