• <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.

    一区在线观看完整版| 桃红色精品国产亚洲av| 亚洲精品乱久久久久久| 国产成+人综合+亚洲专区| 老司机靠b影院| 另类亚洲欧美激情| 99国产综合亚洲精品| 99国产精品免费福利视频| cao死你这个sao货| 国产精品欧美亚洲77777| 天天躁日日躁夜夜躁夜夜| 999久久久国产精品视频| 亚洲精品一二三| 老鸭窝网址在线观看| 狂野欧美激情性xxxx| 在线观看免费高清a一片| 老司机靠b影院| 中文字幕人妻丝袜一区二区| 成人手机av| 国产激情久久老熟女| 两性午夜刺激爽爽歪歪视频在线观看 | 一本大道久久a久久精品| 伊人久久大香线蕉亚洲五| 一二三四社区在线视频社区8| 在线观看午夜福利视频| av一本久久久久| 欧美日韩av久久| 在线观看一区二区三区激情| 母亲3免费完整高清在线观看| 久久久久国产一级毛片高清牌| 天天躁日日躁夜夜躁夜夜| 亚洲伊人色综图| 久久久久久久午夜电影 | 免费观看精品视频网站| 国产欧美日韩一区二区三区在线| 精品久久久久久久毛片微露脸| 在线观看免费午夜福利视频| 中文字幕av电影在线播放| 亚洲精品国产区一区二| 大香蕉久久网| 两人在一起打扑克的视频| 国产成人精品久久二区二区免费| 极品少妇高潮喷水抽搐| 亚洲精品中文字幕一二三四区| 在线av久久热| 日韩免费高清中文字幕av| 午夜免费成人在线视频| 久久久国产精品麻豆| 精品熟女少妇八av免费久了| 久久精品国产综合久久久| 国产亚洲精品久久久久久毛片 | 母亲3免费完整高清在线观看| 精品一品国产午夜福利视频| 久久国产乱子伦精品免费另类| 免费少妇av软件| 无限看片的www在线观看| 老司机靠b影院| 国产精品免费一区二区三区在线 | 99久久国产精品久久久| 国产aⅴ精品一区二区三区波| 欧美 亚洲 国产 日韩一| 中文字幕最新亚洲高清| 中文字幕最新亚洲高清| 亚洲情色 制服丝袜| 日日摸夜夜添夜夜添小说| 国产1区2区3区精品| 女性被躁到高潮视频| 国产视频一区二区在线看| 免费在线观看影片大全网站| 精品亚洲成国产av| 精品国产一区二区三区四区第35| 欧美老熟妇乱子伦牲交| 在线免费观看的www视频| 久久国产精品大桥未久av| av网站免费在线观看视频| 久久九九热精品免费| 婷婷丁香在线五月| 精品少妇一区二区三区视频日本电影| 一本一本久久a久久精品综合妖精| 一区二区三区精品91| 亚洲情色 制服丝袜| 老司机靠b影院| 黄色视频不卡| 亚洲成国产人片在线观看| 欧美 亚洲 国产 日韩一| 国产真人三级小视频在线观看| 男女下面插进去视频免费观看| √禁漫天堂资源中文www| 一级黄色大片毛片| 亚洲精品粉嫩美女一区| 91麻豆av在线| 国产视频一区二区在线看| 99精品欧美一区二区三区四区| 大香蕉久久成人网| 午夜日韩欧美国产| 亚洲国产看品久久| 久久九九热精品免费| 亚洲 欧美一区二区三区| 天天躁夜夜躁狠狠躁躁| 一个人免费在线观看的高清视频| 欧美中文综合在线视频| 亚洲国产欧美日韩在线播放| 免费女性裸体啪啪无遮挡网站| 成人三级做爰电影| 国产99白浆流出| 国产欧美亚洲国产| 老司机午夜福利在线观看视频| 美女扒开内裤让男人捅视频| 国产无遮挡羞羞视频在线观看| 国产一区二区三区综合在线观看| 成人国语在线视频| 国产亚洲av高清不卡| 91精品三级在线观看| 国产精品美女特级片免费视频播放器 | 亚洲色图 男人天堂 中文字幕| 日日夜夜操网爽| 国产精品综合久久久久久久免费 | 涩涩av久久男人的天堂| 免费观看人在逋| 人妻久久中文字幕网| 又紧又爽又黄一区二区| 国产一区二区三区在线臀色熟女 | 美女国产高潮福利片在线看| 亚洲精品美女久久av网站| 精品久久久久久电影网| 国产乱人伦免费视频| 宅男免费午夜| 成人手机av| 一区二区三区国产精品乱码| 日韩精品免费视频一区二区三区| tocl精华| 亚洲va日本ⅴa欧美va伊人久久| 久久久精品免费免费高清| 大香蕉久久成人网| 在线观看免费高清a一片| 亚洲国产毛片av蜜桃av| tocl精华| 免费日韩欧美在线观看| 99re6热这里在线精品视频| 成年人黄色毛片网站| 一区二区三区国产精品乱码| 国产又爽黄色视频| 午夜精品久久久久久毛片777| 91在线观看av| 亚洲视频免费观看视频| 视频区图区小说| 又黄又粗又硬又大视频| 国产精品国产av在线观看| 久久性视频一级片| 国产成+人综合+亚洲专区| 国产免费av片在线观看野外av| 制服诱惑二区| 成年版毛片免费区| av有码第一页| 男女高潮啪啪啪动态图| 国产在线一区二区三区精| 欧美在线一区亚洲| svipshipincom国产片| 免费在线观看影片大全网站| 人成视频在线观看免费观看| 久久 成人 亚洲| tube8黄色片| 在线观看午夜福利视频| 99久久99久久久精品蜜桃| 欧美在线黄色| a在线观看视频网站| xxxhd国产人妻xxx| 国产淫语在线视频| 久久久国产成人免费| x7x7x7水蜜桃| 亚洲精品成人av观看孕妇| 岛国毛片在线播放| 欧美一级毛片孕妇| 18禁观看日本| 国产黄色免费在线视频| 在线观看舔阴道视频| 国产蜜桃级精品一区二区三区 | 久久国产乱子伦精品免费另类| 麻豆av在线久日| 黄色毛片三级朝国网站| 国产淫语在线视频| aaaaa片日本免费| 一级黄色大片毛片| 国产有黄有色有爽视频| 搡老乐熟女国产| 午夜免费成人在线视频| 国产av又大| 国产精品久久久久成人av| 久久天堂一区二区三区四区| 老司机午夜福利在线观看视频| 丰满饥渴人妻一区二区三| 国产一区二区三区在线臀色熟女 | 亚洲欧美激情在线| 日本a在线网址| 午夜免费鲁丝| 精品视频人人做人人爽| 亚洲人成伊人成综合网2020| 久久久久久免费高清国产稀缺| 大陆偷拍与自拍| 成人影院久久| 国产日韩欧美亚洲二区| 99香蕉大伊视频| 女性被躁到高潮视频| 捣出白浆h1v1| 国产高清videossex| 好男人电影高清在线观看| 欧美在线一区亚洲| 一本综合久久免费| 久久久久精品人妻al黑| 熟女少妇亚洲综合色aaa.| 三上悠亚av全集在线观看| 99国产精品一区二区蜜桃av | 黄频高清免费视频| 真人做人爱边吃奶动态| 宅男免费午夜| 欧洲精品卡2卡3卡4卡5卡区| 精品亚洲成国产av| 视频区图区小说| 一级黄色大片毛片| 久久精品国产亚洲av高清一级| 看免费av毛片| 亚洲精品一二三| 女人被狂操c到高潮| 999精品在线视频| 大型av网站在线播放| 国产精品电影一区二区三区 | 人妻丰满熟妇av一区二区三区 | 国产麻豆69| 侵犯人妻中文字幕一二三四区| 91成年电影在线观看| 黄色丝袜av网址大全| 极品少妇高潮喷水抽搐| bbb黄色大片| 在线观看免费视频日本深夜| 男女下面插进去视频免费观看| 免费一级毛片在线播放高清视频 | 国产真人三级小视频在线观看| 久久性视频一级片| 亚洲欧美精品综合一区二区三区| 男女之事视频高清在线观看| 日日摸夜夜添夜夜添小说| 国产蜜桃级精品一区二区三区 | 国产不卡av网站在线观看| 国产成人免费观看mmmm| 精品电影一区二区在线| 两性夫妻黄色片| 在线十欧美十亚洲十日本专区| 成年人免费黄色播放视频| 免费av中文字幕在线| 99久久人妻综合| 又大又爽又粗| 亚洲色图av天堂| 国产成人精品在线电影| 亚洲第一av免费看| 精品一区二区三区四区五区乱码| 飞空精品影院首页| 成人手机av| 国产精品国产高清国产av | 又黄又爽又免费观看的视频| 91av网站免费观看| 天天躁日日躁夜夜躁夜夜| 亚洲熟妇中文字幕五十中出 | 精品少妇一区二区三区视频日本电影| 欧美在线一区亚洲| 亚洲欧美一区二区三区久久| 国产欧美日韩精品亚洲av| 国产精品欧美亚洲77777| 国产男靠女视频免费网站| 一区在线观看完整版| 一级a爱视频在线免费观看| 极品教师在线免费播放| 精品高清国产在线一区| av有码第一页| 日日爽夜夜爽网站| 亚洲久久久国产精品| 欧美国产精品一级二级三级| 国产xxxxx性猛交| 国产成人一区二区三区免费视频网站| 老司机靠b影院| 又大又爽又粗| 亚洲五月色婷婷综合| 高清黄色对白视频在线免费看| 99国产精品99久久久久| 亚洲精品国产精品久久久不卡| 国产精品永久免费网站| av欧美777| 日本vs欧美在线观看视频| 久久婷婷成人综合色麻豆| 日韩人妻精品一区2区三区| 色综合婷婷激情| 天堂√8在线中文| 亚洲欧洲精品一区二区精品久久久| 一级a爱片免费观看的视频| bbb黄色大片| 超碰97精品在线观看| 中亚洲国语对白在线视频| 欧美国产精品va在线观看不卡| 老司机影院毛片| 一级毛片女人18水好多| 高清视频免费观看一区二区| 国产日韩一区二区三区精品不卡| av天堂在线播放| 久99久视频精品免费| 亚洲精品在线美女| 成年人黄色毛片网站| 亚洲欧美日韩另类电影网站| 日韩人妻精品一区2区三区| 99国产精品一区二区三区| 久久国产精品大桥未久av| 一级毛片精品| 久久久水蜜桃国产精品网| 成人18禁高潮啪啪吃奶动态图| 国产主播在线观看一区二区| 午夜成年电影在线免费观看| 久久精品aⅴ一区二区三区四区| 久久草成人影院| av线在线观看网站| 老司机午夜十八禁免费视频| 精品电影一区二区在线| 热99re8久久精品国产| 亚洲精品国产区一区二| 麻豆国产av国片精品| 麻豆乱淫一区二区| 高清毛片免费观看视频网站 | 飞空精品影院首页| 99香蕉大伊视频| 999久久久精品免费观看国产| 中亚洲国语对白在线视频| 成年人午夜在线观看视频| 久久人妻熟女aⅴ| 欧美乱码精品一区二区三区| 精品国产乱码久久久久久男人| 成人精品一区二区免费| 国产男女超爽视频在线观看| 精品卡一卡二卡四卡免费| 亚洲熟妇熟女久久| 一级毛片精品| 国产成人精品在线电影| 91九色精品人成在线观看| 成人免费观看视频高清| 99国产精品99久久久久| 国产亚洲欧美在线一区二区| 亚洲成人手机| 亚洲九九香蕉| 在线十欧美十亚洲十日本专区| 国产伦人伦偷精品视频| 757午夜福利合集在线观看| 精品免费久久久久久久清纯 | 黄片大片在线免费观看| 999久久久精品免费观看国产| 男人的好看免费观看在线视频 | 精品一区二区三卡| 国产成人免费观看mmmm| 精品一区二区三区视频在线观看免费 | 日韩欧美国产一区二区入口| av视频免费观看在线观看| 桃红色精品国产亚洲av| 亚洲欧洲精品一区二区精品久久久| 久久婷婷成人综合色麻豆| 久久人妻熟女aⅴ| 国产人伦9x9x在线观看| 国产精品影院久久| 又大又爽又粗| 美国免费a级毛片| 日韩欧美一区二区三区在线观看 | 一区二区三区国产精品乱码| 欧美黑人精品巨大| 亚洲人成电影免费在线| 亚洲国产精品合色在线| 国产一区二区三区综合在线观看| 高清av免费在线| 国产精品综合久久久久久久免费 | 亚洲一区二区三区不卡视频| 久久精品人人爽人人爽视色| 亚洲午夜理论影院| 久久中文字幕一级| 欧美激情久久久久久爽电影 | 亚洲九九香蕉| 欧美日韩瑟瑟在线播放| 欧美精品人与动牲交sv欧美| 中文亚洲av片在线观看爽 | 日韩欧美国产一区二区入口| 女人被狂操c到高潮| 亚洲国产精品sss在线观看 | 精品国产乱子伦一区二区三区| 99精品欧美一区二区三区四区| 高潮久久久久久久久久久不卡| 亚洲精品美女久久av网站| 亚洲精品av麻豆狂野| 欧美日韩中文字幕国产精品一区二区三区 | 欧美日韩亚洲综合一区二区三区_| 搡老熟女国产l中国老女人| 国产有黄有色有爽视频| 激情在线观看视频在线高清 | 18禁观看日本| 欧美黄色淫秽网站| 成年动漫av网址| 亚洲国产中文字幕在线视频| 久久久精品免费免费高清| 国产人伦9x9x在线观看| 少妇被粗大的猛进出69影院| 中出人妻视频一区二区| 搡老岳熟女国产| 欧美一级毛片孕妇| 国产亚洲精品第一综合不卡| 99香蕉大伊视频| 国产高清视频在线播放一区| 精品无人区乱码1区二区| 极品少妇高潮喷水抽搐| 一进一出抽搐动态| 亚洲精品一卡2卡三卡4卡5卡| 久久狼人影院| tocl精华| 在线观看舔阴道视频| 天天躁夜夜躁狠狠躁躁| 69av精品久久久久久| 免费在线观看日本一区| 日韩精品免费视频一区二区三区| 免费在线观看黄色视频的| 精品少妇一区二区三区视频日本电影| 成在线人永久免费视频| 黄频高清免费视频| 在线视频色国产色| 亚洲人成77777在线视频| 青草久久国产| 免费日韩欧美在线观看| 五月开心婷婷网| bbb黄色大片| 免费看a级黄色片| 老司机福利观看| 国产淫语在线视频| 久久人人爽av亚洲精品天堂| 老汉色av国产亚洲站长工具| 亚洲中文日韩欧美视频| 一夜夜www| 手机成人av网站| 欧美日韩成人在线一区二区| 18禁观看日本| 自拍欧美九色日韩亚洲蝌蚪91| 精品国产一区二区三区四区第35| 久久人妻av系列| 国产一区二区三区综合在线观看| 好看av亚洲va欧美ⅴa在| 在线看a的网站| 黑人巨大精品欧美一区二区mp4| 国产色视频综合| 欧美+亚洲+日韩+国产| 女人精品久久久久毛片| 国产精品香港三级国产av潘金莲| 精品一区二区三区av网在线观看| 久久国产精品大桥未久av| 国产精品乱码一区二三区的特点 | 精品国产乱码久久久久久男人| 色老头精品视频在线观看| 18禁裸乳无遮挡免费网站照片 | 欧美乱妇无乱码| 91字幕亚洲| 亚洲av熟女| 国产精品国产高清国产av | av线在线观看网站| 熟女少妇亚洲综合色aaa.| 欧美激情 高清一区二区三区| 一级毛片精品| 亚洲七黄色美女视频| 亚洲精品一卡2卡三卡4卡5卡| 如日韩欧美国产精品一区二区三区| www.自偷自拍.com| 大片电影免费在线观看免费| 日本精品一区二区三区蜜桃| 丁香六月欧美| 最近最新中文字幕大全免费视频| 在线观看www视频免费| 视频区图区小说| 巨乳人妻的诱惑在线观看| 亚洲精品国产色婷婷电影| 精品人妻熟女毛片av久久网站| 在线国产一区二区在线| 亚洲一码二码三码区别大吗| 超碰97精品在线观看| 日韩一卡2卡3卡4卡2021年| 日韩 欧美 亚洲 中文字幕| 欧美久久黑人一区二区| 黄色视频,在线免费观看| 欧美精品亚洲一区二区| 免费观看a级毛片全部| 久久久国产欧美日韩av| 黄色丝袜av网址大全| 国产精品九九99| 欧美国产精品va在线观看不卡| 欧美日韩成人在线一区二区| 午夜免费成人在线视频| 国产片内射在线| 国产欧美亚洲国产| 女性被躁到高潮视频| 久久九九热精品免费| 女人高潮潮喷娇喘18禁视频| 欧美国产精品一级二级三级| 国产精品.久久久| 欧美激情久久久久久爽电影 | 久久久久久久久久久久大奶| 女人被躁到高潮嗷嗷叫费观| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲av美国av| 欧美乱码精品一区二区三区| 日韩大码丰满熟妇| av中文乱码字幕在线| svipshipincom国产片| 国产精品九九99| 亚洲一区中文字幕在线| 亚洲熟女精品中文字幕| 免费黄频网站在线观看国产| 国产精品一区二区免费欧美| 91av网站免费观看| 亚洲精品在线美女| 久久精品国产99精品国产亚洲性色 | 久久香蕉激情| 99re在线观看精品视频| 黄色怎么调成土黄色| 18禁美女被吸乳视频| 中文亚洲av片在线观看爽 | 日韩欧美一区二区三区在线观看 | www.999成人在线观看| av电影中文网址| 嫁个100分男人电影在线观看| av天堂在线播放| 正在播放国产对白刺激| 黑丝袜美女国产一区| 亚洲精品成人av观看孕妇| 精品无人区乱码1区二区| av视频免费观看在线观看| 亚洲情色 制服丝袜| 啦啦啦 在线观看视频| 女性被躁到高潮视频| 麻豆国产av国片精品| 国产xxxxx性猛交| 免费在线观看日本一区| 69av精品久久久久久| 国产欧美日韩一区二区三区在线| 亚洲国产中文字幕在线视频| 久久精品人人爽人人爽视色| 两个人免费观看高清视频| 亚洲中文日韩欧美视频| 这个男人来自地球电影免费观看| 国产日韩欧美亚洲二区| 黄片播放在线免费| 亚洲国产精品一区二区三区在线| 午夜视频精品福利| 午夜成年电影在线免费观看| 露出奶头的视频| 一二三四在线观看免费中文在| 日本wwww免费看| 91麻豆精品激情在线观看国产 | 新久久久久国产一级毛片| 久久国产精品人妻蜜桃| 黄色视频,在线免费观看| 建设人人有责人人尽责人人享有的| 色婷婷av一区二区三区视频| 久久精品国产清高在天天线| 天天躁日日躁夜夜躁夜夜| 久久性视频一级片| 日本欧美视频一区| 精品国内亚洲2022精品成人 | 精品久久久精品久久久| a在线观看视频网站| 久久国产乱子伦精品免费另类| 人人妻人人添人人爽欧美一区卜| 国产成人精品久久二区二区91| 精品电影一区二区在线| 国产片内射在线| 国内毛片毛片毛片毛片毛片| 亚洲欧美日韩另类电影网站| 首页视频小说图片口味搜索| 天堂√8在线中文| 97人妻天天添夜夜摸| 99久久国产精品久久久| 国产一卡二卡三卡精品| 久久精品亚洲精品国产色婷小说| 麻豆国产av国片精品| 久久精品aⅴ一区二区三区四区| 亚洲av片天天在线观看| 青草久久国产| 午夜福利欧美成人| 午夜日韩欧美国产| 久久性视频一级片| 一a级毛片在线观看| 在线国产一区二区在线| 青草久久国产| 精品视频人人做人人爽| 亚洲午夜理论影院| 制服诱惑二区| 久热这里只有精品99| 黄色毛片三级朝国网站| 国产成人影院久久av| 亚洲美女黄片视频| av国产精品久久久久影院| 亚洲精华国产精华精| 久久久国产欧美日韩av| 日日爽夜夜爽网站| 最新的欧美精品一区二区| 天天影视国产精品| 日日爽夜夜爽网站| 99riav亚洲国产免费| 巨乳人妻的诱惑在线观看| 亚洲中文日韩欧美视频| 青草久久国产| 午夜福利乱码中文字幕| 亚洲一区二区三区不卡视频| 中文字幕人妻熟女乱码| 亚洲精品一二三| 女性被躁到高潮视频| 久久久精品国产亚洲av高清涩受| 人人妻人人爽人人添夜夜欢视频| 欧美日韩国产mv在线观看视频| 久久九九热精品免费| 91成人精品电影|