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

    Effects of Mg/PTFE pyrotechnic compositions on reignition characteristics of base bleed propellants and heating mechanism

    2022-02-15 03:47:54XiochunXueLongzeYonggngYu
    Defence Technology 2022年1期

    Xio-chun Xue ,Long-ze M ,Yong-gng Yu ,*

    a School of Energy and Power Engineering,Nanjing University of Science and Technology,Nanjing,Jiangsu,210094,China

    b BYD Company Limited,Shenzhen,Guangdong,518118,China

    Keywords:Thermal effect High-temperature combustion gas Mg/PTFE pyrotechnic Compositions Reignition

    ABSTRACT The effects of magnesium/polytetrafluoroethylene (Mg/PTFE) pyrotechnic compositions on the coupled flow field and reignition mechanism are important aspects governing the performance and range of base bleed projectiles (BBPs).Owing to a decrease in pressure and temperature when the BBP leaves the muzzle,rapid depressurization occurs,which extinguishes the base bleed propellant.The Mg/PTFE pyrotechnic composition pressed in the igniter of the base bleed unit (BBU) provides additional energy to the BBU via a chemical reaction.Thus,the extinguished base bleed propellant is reignited under the effect of high-temperature combustion gas jets from the igniter.In this study,a numerical analysis is conducted to evaluate the effects of PTFE and Mg granularity as well as Mg/PTFE pyrotechnic compositions.Owing to the rapid depressurization,the temperature and pressure was found to decrease for different Mg/PTFE pyrotechnic compositions.However,the depressurization time increased as the PTFE granularity increased,the Mg granularity decreased,and the Mg content increased.When the pressure in the combustion chamber of the BBU decreased to the atmospheric pressure,the combustion gas jets from the igniter expand upstream(rather than downstream).However,these combustion gas jets exhibit different axial and radial expansion characteristics depending on the pyrotechnic compositions used.The results show that the reignition delay time,td,of the base bleed propellant was 377.608,94.27,387.243,523.966,and 221.094 ms for cases A-E,respectively.Therefore,it was concluded that the Mg/PTFE pyrotechnic composition of case B was the most beneficial for the reignition of the base bleed propellant,with the earliest addition of energy and mass to the BBP.

    1.Introduction

    Aerodynamic projectiles such as artillery shells experience significant drag owing to the formation of a low-pressure and lowvelocity recirculation region near their bases [1,2].This drag,also termed as base drag,can account for up to two-thirds of the total drag acting on the projectile when it is launched from the muzzle.As a solution,base bleed units(BBUs)are conventionally used to reduce base drag,which,in turn,helps increase the range of the projectile[3-5].Generally,projectiles equipped with a BBU are termed as base bleed projectiles (BBPs).The base bleed jets in such projectiles are mainly composed of the combustion gases of the ammonium perchlorate/hydroxyterminated polybutadiene(AP/HTPB)propellant in the BBU.However,in case of an unsteady flow field in the BBU,rapid depressurization occurs,and the AP/HTPB propellant is extinguished as soon as the BBP exits the muzzle [6].However,this extinguished AP/HTPB propellant can be reignited,depending on the coupled flow field of the high-temperature combustion gas from the igniter in the BBU and the propellant gas in the gun.In an actual BBU,magnesium/polytetrafluoroethylene(Mg/PTFE) pyrotechnic compositions[7]are pressed as the igniter.Once the AP/HTPB propellant is reignited by combustion gas jets from the igniter,the base drag continues to decrease with increasing base bleed.However,uncertainties in the reignition mechanism of the AP/HTPB propellant can reduce the accuracy of the BBP and cause significant variations in its range.Consequently,the effects of different Mg/PTFE pyrotechnic compositions on the complex flow field characteristics of a BBU and the thermal mechanisms of the AP/HTPB propellant during rapid depressurization are essential parameters that need to be considered for improving the accuracy of BBPs.

    As mentioned above,the AP/HTPB propellant in the BBU also undergoes depressurization when the BBP exits the muzzle.Focusing on the combustion of the AP/HTPB propellant under unsteady conditions,Steinz [8] carried out depressurization extinguishment experiments;using a high-speed video camera,it was observed that the gaseous flame was extinguished immediately owing to adiabatic expansion.Suhas [9] presented the unsteady combustion analysis of composite propellants and proposed a mechanism for the unsteady combustion of composite propellants subjected to rapid depressurization.Surzhikov [10] analyzed combustion behaviors under periodic high-frequency variations in pressure and found that the distribution of the average mass burning rate was non-uniform throughout the surface of the solid propellant.Wang [11] presented a laser technique for measuring the transient burning rate of a solid propellant;the aim of that study was to develop a depressurization system,where the starting time of depressurization could be set in advance by simply altering the initial rate of depressurization.Furthermore,Jackon [12]established an experiment to measure the burning rate of a solid propellant under a decreasing pressurization rate.Ye[6]conducted unsteady combustion experiments on the AP/HTPB propellant subjected to depressurization by using a semi-hermetic explosive device;three types of unsteady combustion characteristics-continuous combustion,flame extinguishment,and flame reignition-were observed.

    Despite such studies on the unsteady combustion mechanism and the thermal effect of solid propellants,research pertaining to the combustion of Mg/PTFE-based pyrotechnic compositions has been significantly limited.Kuo [13] summarized the ignition and combustion characteristics of pyrotechnic compositions for a solid rocket engine.Deyong [14,15] reported theoretical studies on simulations of the combustion of pyrotechnic compositions.The effects of changing the formulations on the distribution of reaction temperatures and the heat effect were determined by providing a new chemical kinetic mechanism of 18 reactions and 17 species for the combustion of pyrotechnic compositions.Based on this,the detailed equilibrium chemistry of Mg/PTFE-based pyrotechnic compositions was examined for over 47 species and various initial conditions[16].Our previous studies[17,18]also modeled the flame evolution characteristics of an igniter by adopting Mg/PTFE-based pyrotechnic compositions under atmospheric conditions by simplifying Deyong’s detailed kinetics mechanism,and the obtained simulation results were in agreement with the experimental results.However,the abovementioned studies only focused on the burning characteristics of Mg/PTFE-based pyrotechnic compositions under a steady environment.Details regarding the coupled flow field characteristics for various combustion gas interactions under rapid pressurization processes are still scarce.Lu[19]studied the transient flame characteristics and performances of different types of igniters under rapid depressurization through semi-closed bomb experiments and found that the flame structure of the Mg/PTFE-based pyrotechnic igniter was the most stable for the longest duration.However,it is difficult to determine the thermal effect and reignition mechanism of the propellant in the Mg/PTFEbased pyrotechnic igniter of high-temperature and high-pressure BBUs under rapid depressurization by using nonvisual experiments alone;numerical investigations are also necessary to better understand such a complex system.

    The remainder of this paper is organized as follows.First,the governing equations are presented.Here,we assume that the unsteady compressible flow equations account for the combustion of Mg/PTFE pyrotechnic compositions.Subsequently,the flow simulations are discussed in detail.The results highlight the complex coupled flow field characteristics for various combustion gases in the BBU with different Mg/PTFE pyrotechnic compositions when rapid depressurization occurs.Finally,the mechanism of the thermal effects on the extinguished AP/HTPB propellant is investigated by determining the optimal Mg/PTFE composition that is beneficial for reigniting the AP/HTPB propellant and facilitating the base bleed in projectiles.

    2.Physical model for computation

    2.1.Problem description

    Details of the BBU in the BBP are presented in Fig.1;the BBU comprises a shell,base bleed propellant,combustion chamber,nozzle,and pyrotechnic igniter.There are two types of combustibles in the BBU:the AP/HTPB propellant,which is used as the base bleed propellant,and the ignition powder in the igniter.Generally,Mg/PTFE is adopted as the ignition powder and pressed in the igniter.Initially,when the BBP is in the muzzle,the hightemperature and high-pressure gas generated by the combustion of the gun propellant from the gun breech can enter the BBU through the nozzle.Here,it is assumed that the gun propellant is completely combusted in the gun breech;hence,only the combustion of Mg/PTFE in the igniter is considered,and the chemical reaction rate is calculated using the Arrhenius law.In the present study,we focus on the effects of the granularity and mass fraction of Mg/PTFE components on the coupled flow field of hightemperature combustion gases from the gun propellant and the Mg/PTFE pyrotechnic composition in the combustion chamber.With regard to an actual BBU,when the BBP exits the muzzle,rapid depressurization occurs in the combustion chamber,accompanied by a decrease in pressure and temperature,owing to the pressure difference with respect to the atmosphere.According to the experimental results in Ref.[6],when the BBP leaves the muzzle,AP/HTPB is extinguished within a significantly short period of time due to the flame separation effect.Thus,the combustion of the base bleed propellant is neglected in this study.However,the base bleed propellant is reignited under the effect of the high-temperature combustion gas from the igniter at the later stages of depressurization.Accordingly,the thermal effects of different igniter combustion gases on the base bleed propellant are considered in this study.

    2.2.Governing equations

    An ideal compressible and viscous gas in two dimensions is described by the continuity and transient Naiver-Stokes equations,as follows:

    Fig.1. Schematic of a base bleed unit.

    where W,F,and Gdenote the vectors of the conservative variables and vectors of the fluxes,respectively.

    where E is the total energy;τ,τ,τ,and τare the shear stress tensors;and qand qare the vectors of heat flux.

    where the subscripts x and r represent the axial and radial directions of combustion,respectively;u and v are the axial and radial velocities,respectively;ρ is the density of the combustion gas mixture;ρ,h,and ωrepresent the density,static enthalpy,and net production rate of species i produced via chemical reactions,respectively;and p is the pressure,which can be expressed by the ideal gas state equation according to the abovementioned assumptions:

    where R and T are the universal gas constant and temperature,respectively;Yis the mass fraction,and it can be expressed as Y=ρ/ρ.Mis the molar mass of species i.

    In a mixture of gases,the mass diffusion coefficient,D,of species i can be expressed as follows:

    where μand Scare the turbulent dynamic viscosity and the turbulent Schmidt number,respectively.For species i,Drepresents the mean mass diffusion coefficient as

    Here,Xis the mole fraction of species i,and Xis the mole fraction of species j.Drepresents the binary mass diffusion coefficient of species i in species j.The thermal conductivity coefficient λ is calculated as follows:

    Pris the turbulence Prandtl number,and λis the mean mixture thermal conductivity,which is expressed as

    The specific heat capacity of the mixture is expressed as

    cdenotes the specific heat capacity of species i,and it can be calculated as

    where μ is the dynamic viscosity.

    The mean mixture dynamic viscosity μ,is expressed as

    Here,μis the viscosity coefficient of species i.

    2.3.Turbulence model

    In this study,the compressible flow in the BBU was numerically investigated.Computations were conducted using axisymmetric mass averaged Navier-Stokes equations with the realizable k-ε turbulence model [20].The results for this turbulent model were validated using the experiments described in the following section.

    2.4.Chemical kinetic model

    A careful review of existing experimental results[21]suggested that the combustion gases are strongly pressed on the surface of the igniter by the high-pressure combustion gas formed by the hightemperature and high-pressure gun propellant during the depressurization process.Consequently,fluctuations in the turbulence are weak.For a premixed reactant,combustion mainly occurs near the exit of the igniter,and the chemical reactions are key factors influencing this combustion process.Accordingly,the chemical reaction source term,ω,is calculated using a finite-rate kinetic model,and it can be expressed as

    where γrepresents the third-body efficiency of the jth species in the rth reaction.

    The forward rate constant is calculated using the Arrhenius law for reaction r:

    Here,A,βand Erepresent the pre-exponential factor,temperature exponent,and activation energy,respectively.

    In the present study,we assume that the products of the solid evaporation and decomposition of Mg/PTFE pyrotechnic compounds are only Mg vapor and CFgas injected into the combustion chamber of the BBU and as the initial gaseous reactants in the combustion chamber.According to Ref.[17],the chemical reactions are reversible;they are presented in Table 1.For reaction r,the backward rate constant is calculated using the forward rate constant:

    Table 1 Chemical reaction mechanisms.

    where the equilibrium constant,K,is expressed as

    Here,pis 101325 Pa Sis the entropy of species i,and hrepresents the enthalpy of the ith species.

    3.Computational solution

    3.1.Computational methodology

    The finite volume scheme [22],which involves integrating the governing equations for each control volume and obtaining discrete equations that conserve each quantity on a control volume basis,is applied to simulate the coupled flow field characteristics of the high-temperature combustion gases generated by the gun propellant and the igniter pyrotechnics when subjected to rapid depressurization.The governing equations are discretized using the following methods with the AUSM+scheme to achieve the best accuracy [23].

    3.2.Computational domain and grid arrangement

    Fig.2 depicts the computational domain with a twodimensional axisymmetric model and the results of grid independence.Combined with the actual BBU,the height and length of the igniter surface are ma=15 mm and ab=20 mm,respectively.The length of the base bleed propellant was cd=60 mm and ef=9 mm.The length of the axis is lm=1500 mm,and the radial distance of the computational domain is kl=800 mm.The remaining sizes are shown in Fig.2.For convenience,we divided the entire domain into two zones:Zone 1 represents the high-pressure inner flow zone and Zone 2 represents the outer flow zone.Computations were conducted on the grids in the range of 130000,420000,and 950000 nodes by using a block-structured grid.After grid independence verification,as shown in Fig.2(c),grid sizes with 420000 nodes or more afford less changes in the profiles of pressure distribution in the flow field of the BBU.Consequently,it is considered that 420000 nodes are sufficient to obtain grid-independent solutions for the present study.

    3.3.Initial condition

    Fig.2. Computational model.

    At the time when the BBU has just left the muzzle,the hightemperature and high-pressure gun propellant combustion gases occupy Zone 1 completely,with the total pressure p,temperature T,and specific heat ratio γ.Zone 2 comprises the atmospheric environment,where the absolute pressure and temperature are pand T,respectively.These initial parameters are listed in Table 2.Moreover,ρis the density of the pyrotechnics,and Ais the end face area of the igniter for an actual BBU.

    3.4.Boundary conditions

    As shown in Fig.2,the central axis (ml) is the axisymmetric boundary,and the shell wall and surface of the base bleed propellant are adiabatic,without sliding walls.The boundary of the flow field for the computational domain is the pressure outlet,and the end face of the igniter is the inlet of the mass flow rate defined by ˙m=ρ˙rA.Here,the burning rate of the Mg/PTFE pyrotechnic composition,defined according to the experiments conducted by Kubota [24],is given by ˙r=ap,where a is the burning rate coefficient,p is the absolute pressure in Zone 1,and n is the pressure exponent.

    To define the effect of the particle sizes of Mg and PTFE and the content of Mg on the burning rate of the Mg/PTFE pyrotechnic grain,Kubota [24] studied the burning rate characteristics of various Mg/PTFE pyrotechnic compositions.The test results of five such compositions are shown in Fig.3.Thus,the burning rate formulas for these Mg/PTFE compositions can be obtained by fitting the experimental results in this figure,as presented in Table 3.In this study,we focus on the effect of these compositions of Mg/PTFE on the flow field characteristics of the BBU as well as the thermal feedback mechanism affecting the base bleed propellant under depressurization.

    3.5.Validation of the models

    A simulation considering the experimental conditions was performed in advance in order to validate the present models and computational solutions for the prediction of the thermal effect of high-temperature Mg/PTFE compositions with different granularities and mass ratios in the BBU when subjected to rapid depressurization.The experiment conducted by Ma [21],where the BBU was evaluated under such extreme conditions,was chosen as a reference;further details and information regarding the experimental device and conditions can be found in Ref.[21].Fig.4(a)shows the pressures and depressurization rates in the combustion chamber of the BBU,obtained using the present models and the experiment,with respect to time.Evolutions of the flame pattern near the nozzle are also obtained,as depicted in Fig.4(b),using the abovementioned experimental results and the corresponding numerical results below.From the figures,overall,it can be considered that the present numerical models were validated successfully based on the previous experimental study.In particular,the evolution of the flame pattern was in good agreement with the experimental observations.

    Fig.3. Burning rate of Mg/PTFE.

    Table 3 Parameters of Mg/PTFE compositions and the burning rate.

    Table 2 Initial parameters.

    Fig.4. Comparison of computational and experimental data.

    4.Results and discussion

    In actual BBUs,the flow field structure,particularly the thermal effect of the high-temperature combustion gas from the igniter of the AP/HTPB base bleed propellant,plays a significant role in the performance and accuracy of BBPs.The choice of pyrotechnic composition pressed in the igniter is an important factor determining this flow field structure and the thermal feedback affecting the base bleed propellant.Consequently,through numerical simulations,we study the detailed effects of the granularity and mass fraction of the pyrotechnic composition in the BBU when the BBP exits the muzzle.The results are presented in the following sections.

    4.1.Depressurization characteristics in BBU

    Fig.5 depicts the pressure distributions along the central axis of the combustion chamber at 0.9 ms and 1.1 ms for different Mg/PTFE compositions.Once the BBP leaves the muzzle,the pressure in the combustion chamber decreases significantly.The pressure variations under all these cases demonstrate identical trends;the pressure in the combustion chamber continues to decrease,falling below the atmospheric pressure,due to the continuous expansion of the combustion gas during the later stages of depressurization.Consequently,the igniter combustion gas jets move upstream due to the negative pressure difference,thereby improving radial expansion;this is beneficial for the heating of the base bleed propellant(located on both sides of the BBU)by the high-temperature igniter combustion.However,the depressurization times for these five cases are different.Table 4 presents the time required for the combustion chamber pressure to decrease to the atmospheric pressure when using different granularities and mass fractions of the pyrotechnic compositions;these results were obtained by analyzing the pressure distributions at different times on the basis of the simulation results.The depressurization times are in the range of 1.0-1.1 ms.The depressurization time was noted to increase with an increase in the PTFE granularity,Mg granularity,and mass fraction,as presented in Table 5.According to Kubota[24],the burning rate of the Mg/PTFE igniter increases with a decrease in the Mg particle size and an increase in the PTFE particle size and Mg content.The rate of combustion gas production increases with an increase in the burning rate of Mg/PTFE.Consequently,the depressurization time as well as the increasing rates of the depressurization time increase with combustion gas production,due to an increase in the Mg/PTFE reaction rate.In other words,the depressurization time is determined by the burning rate of Mg/PTFE.This conclusion in combination with Fig.3 indicates that the burning rate for case E is the highest;accordingly,the depressurization time and the increasing rate of the depressurization time in this case are the highest,as presented in Tables 4 and 5.

    Table 4 Depressurization time.

    Fig.5. Pressure distributions along the central axis of the BBU.

    Table 5 Increasing rates of depressurization time.

    4.2.Flow field characteristics of BBU under depressurization

    The observed thermal feedback affecting the base bleed propellant in the combustion chamber of the BBU is attributed to the production of the igniter combustion gas and the coupled interaction effect between the gun propellant combustion gas and the igniter combustion gas.When the BBP has just exited the muzzle,the combustion chamber is filled with high-temperature gun propellant combustion gases.With the production of the hightemperature combustion gas from the igniter,the temperature of the combustion chamber continues to increase.However,rapid depressurization occurs due to the difference in the combustion chamber pressure and the atmospheric pressure.Thus,the pressure and temperature of the gun propellant combustion gas in the combustion chamber rapidly decrease with the expansion of the combustion gas (and its subsequent release into the atmosphere).In addition,the axial expansion of the igniter combustion gas is greater than its radial expansion due to the axial power provided by the pressure difference.Consequently,the complex coupled flow field in the BBU determines whether the continuous heat from the combustion in the igniter is received by the base bleed propellant.With regard to the effects of different Mg/PTFE pyrotechnic compositions on the flow field characteristics of the BBU,Fig.6 depicts the density and temperature distributions in the BBU under depressurization,while Fig.7 depicts the mass fraction distribution of MgF,the product of the high-temperature igniter combustion,for different cases.

    When depressurization occurs,the temperature in the combustion chamber decreases with the expansion of the gun propellant gas.However,the high-temperature combustion gases from the igniter are pressed on the end face of the igniter initially,and gradually expand with the expansion of the gun propellant combustion gas at a later time.The axial expansion of the igniter combustion gas is considerably faster than its radial expansions;therefore,the gun propellant is gradually surrounded by lowtemperature gun propellant gases.For cases A,C,and D,the characteristics of the combustion gas jets from the igniter remain the same even under depressurization,and small-and large-sized vortices exist when the pressure in the combustion chamber decreases to the atmospheric pressure.This induces stronger Kelvin-Helmholtz instability due to the velocity difference between the gun propellant combustion gas and the igniter combustion gas.However,for cases B and E,the small-sized vortices disappear and only large-sized vortices remain,resulting in weaker Kelvin-Helmholtz instability.The radial expansion of the combustion gas is significantly faster in cases B and E,as compared to that that for the other cases,which implies that the igniter heats the base bleed propellant earlier.This is because,as shown in Fig.3 and Table 4,the burning rates of the igniter for cases B and E occupy the absolute advantage gradually under the depressurization process.Subsequently,the quantities of combustion gases produced by the igniter also increase with the depressurization time.In summary,with an increase in PTFE granularity,a decrease in Mg granularity,and an increase in the Mg content of the Mg/PTFE pyrotechnic composition,the Kelvin-Helmholtz instability is weakened,accompanied by significantly quicker expansions along the axial and radial directions of the combustion chamber.

    Fig.8 and Fig.9 depict the distributions of the combustion chamber temperature and the mass fraction of the main product from the igniter combustion,i.e.,MgF,at high temperatures along the central axis of the combustion chamber at different times,respectively.As depicted in Figs.6 and 7,because the hightemperature combustion gases from the igniter are pressed onto the end surface of the igniter and cannot completely expand along the combustion chamber at 0.1 ms,the mass fractions of MgF,which is the main product of the igniter combustion gas,are the highest near the end surface of the igniter for all five cases.The temperatures increase sharply to 2180 K near the igniter and then decrease significantly along the central axis of the combustion chamber for cases A,C,and D at 0.1 ms.In contrast,the temperatures for cases B and E increase sharply to 2050 K and then decrease gradually at 0.1 ms.This is because the burning rates for cases B and E increase with an increase in the production rate of the combustion gas,and the combustion gas acquires greater pressure,allowing it to expand downstream.This results in the hightemperature zones expanding downstream along the axial direction,instead of these zones accumulating near the igniter at the initial stages.At 0.6 ms,the mass fraction peaks of MgFfor case E do not indicate any change and maintain 0.5,indicating that the igniter combustion gases are in a state of chemical equilibrium.The corresponding temperature peak of the igniter combustion gas for case E is also the highest under the highest burning rate and the chemical equilibrium state.From 0.9 ms onwards,except for case C,which exhibits the lowest burning rate,the igniter combustion gases are in a state of chemical equilibrium,and the corresponding mass fraction peaks of MgFdo not exhibit any changes with respect to time.At 1.1 ms,case C reaches chemical equilibrium,with a constant mass fraction peak of MgF.In summary,the depressurization time increases with an increase in the PTFE particle size and a decrease in the Mg particle size(as presented in Table 4),and the pressure of the flow field decreases to the atmospheric pressure the earliest,accompanied by a decrease in the PTFE particle size and an increase in the Mg particle size.Thus,the temperature of the igniter combustion gas decreases with a decrease in the pressure and the depressurization time in the flow field;however,the slowest decrease in temperature was observed for the slowest decrease in pressure.Consequently,the temperature of the igniter combustion gas increases with an increase in the PTFE particle size and a decrease in the Mg particle size under depressurization.

    Fig.6. Density and temperature distributions in the BBU.

    Fig.7. Mass fraction distributions of MgF2 in the BBU.

    Fig.8. Temperature distributions along the central axis of the BBU.

    Fig.9. Mass fraction distributions of MgF2 along the central axis of the BBU.

    With regard to the secondary ignition phenomenon,the temperatures of the combustion gas acting on the base bleed propellant surface along the r=0.022 m axial direction are calculated at different times under depressurization for the five cases;the results are depicted in Fig.10.At 0.1 ms,the pyrotechnic compositions on this temperature have a negligible effect.At 0.6 ms,the temperature of the combustion gas acting on the base bleed propellant surface is the highest for case E and the lowest for case C.The difference between the temperatures acting on the base bleed surface in cases E and C is approximately 4 K.At 0.9 ms,the temperature for case E is the highest,while the temperature difference between cases E and C increases to approximately 12 K.At 1.1 ms,the combustion gas temperature acting on the base bleed propellant surface for cases B and E are significantly higher than those for the other cases.Furthermore,the temperature difference between cases E and C increases to approximately 23 K.In other words,as the PTFE granularity decreases and the Mg granularity and content increase,the temperature of the combustion gas acting on the base bleed propellant surface decreases.As a result,the heating of the base bleed propellant afforded by the combustion gas is decreased,and the temperature of the base bleed propellant surface decreases with the temperature of the combustion gas.As reported in Ref.[25,26],the critical ignition temperature of the base bleed propellant is 650 K;thus,the base bleed propellant is reignited when the local temperature increases to 650 K;thereafter,this temperature remains above 650 K under continuous heating.

    Fig.10. Temperature distributions on the surface of the base bleed propellant.

    4.3.Reignition delay under the heating effect of igniter combustion gas

    According to the abovementioned analysis,the hightemperature combustion gas jets from the igniter continuously heat and even ignite the base bleed propellant as the depressurization process progresses.To reveal the reignition characteristics of the base bleed propellant under this continuous heating afforded by the high-temperature igniter combustion gas under depressurization for different cases,Fig.11 depicts the distributions of the mass fraction of MgF,the product of the high-temperature combustion in the igniter,and the temperature in the combustion chamber after its pressure decreases to the atmospheric pressure.The combustion gas jets from the igniter move upstream (rather than downstream) under the negative pressure difference due to the inertia effect of the combustion gas.The axial and radial expansions of the igniter combustion gas jets for different pyrotechnic compositions exhibit two significantly different trends.One trend refers to cases A,C,and D,where small,high-temperature vortices are formed on the end face of the igniter and then expand to the base bleed surface under the heating effect along the radial direction.It is clear that the combustion gas jets move back into the BBU along the axial direction with the upstream expansion of the igniter combustion gas jets,which is caused by the negative pressure difference.Thus,radial expansion is relatively enhanced.The second trend refers to cases B and E.The igniter combustion gas jets directly expand toward the base bleed propellant along the radial direction of the BBU and undergo heating with the downstream expansion of the igniter combustion gas jets.Thereafter,the igniter combustion gas jets along the axial direction of the BBU shrink,making them extremely weak.However,the igniter combustion gas jets can still be released into the atmosphere from the BBU.Even under different characteristics of the igniter combustion jet,the igniter combustion gas jets all reach the base bleed surface and heat it from 1.35 ms to 1.40 ms for all the five cases studied.

    To better understand the heating effects under different pyrotechnic compositions,the reignition process is divided into two stages depending on the moment when the combustion gas jets from the igniter contact the surface of the base bleed propellant.The two stages are as follows:the delay stage of the combustion gas jet from the igniter,t,and the heating time of the base bleed propellant,t.Here,tis the period between the time when the BBU exits the muzzle and the time when the high-temperature igniter combustion gases come in contact with the surface of the base bleed propellant.The heating time for the base bleed propellant,t,refers to the period between the time when the igniter combustion gases come in contact with the surface of the base bleed propellant and the time when the base bleed propellant is reignited under the continuous heating.Thus,the reignition delay time,t,can be expressed as t=t+t.Table 6 lists the delay stages of the combustion gas jet from the igniter,t,for all the five cases considered.Fig.12 depicts the corresponding highest temperature distributions of the combustion gas on the surface of the base bleed propellant for different times.As shown in Fig.12,the equations for the highest temperature on the surface of the base bleed propellant with respect to time can be fitted for different pyrotechnic compositions,as shown in Eqs.(35)-(39).

    To determine the heating time of the base bleed propellant,t,we assume that the surface temperature of the base bleed propellant is T.The heat flux from the igniter combustion gas jets acting on the base bleed surface can then be expressed as

    Fig.11. Distributions of the mass fractions of MgF2 and temperatures after rapid depressurization.

    Table 6 Delay time of combustion gas jets from the igniter.

    5.Conclusion

    Fig.12. Highest temperature distributions for r=0.022 m after rapid depressurization.

    Fig.13. Highest heat flux variations in igniter combustion gas jets with respect to time.

    A clear understanding of the effect of different Mg/PTFE pyrotechnic compositions on the coupled flow field characteristics of the gun propellant combustion gas and the igniter combustion gas is essential.In particular,the heating of the base bleed propellant by the high-temperature igniter combustion gas is crucial for achieving suitable performance of the BBPs.In this work,numerical models were developed,and their accuracy was validated using a previously reported experiment.Based on the numerical results,depressurization occurs as soon as the BBP exits the muzzle.When this depressurization occurs,the temperature and pressure were found to decrease for the different Mg/PTFE pyrotechnic compositions studied;moreover,the depressurization time increased when the PTFE granularity increased,the Mg granularity decreased,and the Mg content increased.Meanwhile,the lower the PTFE granularity and the higher the Mg granularity and content,the lower the temperature of the combustion gas acting on the base bleed propellant surface.Thus,the heating of the base bleed propellant by the combustion gas is reduced;the temperature of the base bleed propellant surface decreases with that of the combustion gas.

    In particular,after the pressure in the combustion chamber decreases to the atmospheric pressure,the combustion gas jets from the igniter tend to move upstream(rather than downstream)under the negative pressure difference resulting from the inertia effect of the combustion gas.The axial and radial expansions of the igniter combustion gas jets are significantly different for the various pyrotechnic compositions used.Finally,the reignition delay times tof the base bleed propellants were calculated to be 377.608,94.27,387.243,523.966,and 221.094 ms for cases A-E,respectively.These results indicate that the reignition delay time decreases as the PTFE granularity increases and the Mg granularity decreases.Consequently,case B was determined to be the most beneficial for the reignition of the base bleed propellant in the BBU,and for adding mass and energy to the BBP.

    Future research should focus on identifying the coupled flow field characteristics and the effect of adding energy and mass after the ignition of the base bleed propellant AP/HTPB under the heating afforded by the high-temperature igniter combustion gas.Moreover,studies should also attempt to enhance the performance and range of BBPs using measures for drag reduction.

    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.

    This work was supported by the Fundamental Research Funds for the Central Universities (No.30918011324).

    日本 av在线| 在线播放国产精品三级| 国产精品久久视频播放| 夜夜爽天天搞| 免费看十八禁软件| 午夜福利视频1000在线观看| 国产精品嫩草影院av在线观看 | 精品99又大又爽又粗少妇毛片 | 性色av乱码一区二区三区2| 欧美成人a在线观看| 久久草成人影院| netflix在线观看网站| 一进一出抽搐gif免费好疼| 久久这里只有精品中国| 99久久无色码亚洲精品果冻| 亚洲国产中文字幕在线视频| 国内精品一区二区在线观看| 3wmmmm亚洲av在线观看| 性色av乱码一区二区三区2| 日韩有码中文字幕| 国产av麻豆久久久久久久| 精品国产亚洲在线| 色综合亚洲欧美另类图片| 久久亚洲真实| 亚洲第一电影网av| 黄色日韩在线| 啦啦啦免费观看视频1| 99国产综合亚洲精品| 操出白浆在线播放| 99久久精品一区二区三区| 国产一级毛片七仙女欲春2| 久久精品国产清高在天天线| 十八禁人妻一区二区| 久久伊人香网站| 国产精品一区二区三区四区免费观看 | 99在线人妻在线中文字幕| 亚洲最大成人手机在线| 一级毛片女人18水好多| 国产高清视频在线播放一区| 波多野结衣高清作品| 亚洲av免费高清在线观看| 九色国产91popny在线| 麻豆成人av在线观看| 婷婷六月久久综合丁香| www.www免费av| svipshipincom国产片| 国产精品一区二区三区四区免费观看 | 国产麻豆成人av免费视频| 窝窝影院91人妻| 精品国产超薄肉色丝袜足j| 亚洲,欧美精品.| 18禁裸乳无遮挡免费网站照片| 成人18禁在线播放| 亚洲成人免费电影在线观看| 十八禁网站免费在线| 成年人黄色毛片网站| 国产高清三级在线| 精品一区二区三区视频在线观看免费| 一二三四社区在线视频社区8| 国产精品嫩草影院av在线观看 | 日韩有码中文字幕| 麻豆一二三区av精品| 国产乱人视频| 国产精品,欧美在线| www.色视频.com| 午夜视频国产福利| 美女被艹到高潮喷水动态| 欧美日韩国产亚洲二区| 在线观看免费视频日本深夜| 亚洲欧美日韩高清在线视频| а√天堂www在线а√下载| 最新中文字幕久久久久| 亚洲精品粉嫩美女一区| 在线观看午夜福利视频| 老司机午夜十八禁免费视频| 国产视频内射| 欧美色欧美亚洲另类二区| 久久九九热精品免费| 久久久久久久精品吃奶| 午夜亚洲福利在线播放| 一边摸一边抽搐一进一小说| 毛片女人毛片| 亚洲五月天丁香| 三级毛片av免费| 最近最新中文字幕大全电影3| 在线观看舔阴道视频| 国产老妇女一区| av福利片在线观看| 亚洲av成人不卡在线观看播放网| 在线观看一区二区三区| 久久精品综合一区二区三区| 嫩草影院入口| 国产视频内射| 日韩欧美在线乱码| 99视频精品全部免费 在线| 国产真实伦视频高清在线观看 | 国产久久久一区二区三区| 国产又黄又爽又无遮挡在线| 国产精品99久久99久久久不卡| 国产乱人视频| 搡老岳熟女国产| 久久亚洲精品不卡| 舔av片在线| 桃红色精品国产亚洲av| xxxwww97欧美| 精品久久久久久久久久久久久| 精品99又大又爽又粗少妇毛片 | 久久亚洲真实| 欧美成人免费av一区二区三区| 欧美日本视频| 婷婷精品国产亚洲av| 色哟哟哟哟哟哟| 午夜影院日韩av| 日本 av在线| 90打野战视频偷拍视频| 狠狠狠狠99中文字幕| 成人无遮挡网站| 啪啪无遮挡十八禁网站| 99在线视频只有这里精品首页| 亚洲成人久久爱视频| 狂野欧美激情性xxxx| 国内久久婷婷六月综合欲色啪| 成人一区二区视频在线观看| 老司机深夜福利视频在线观看| 亚洲人成网站在线播放欧美日韩| 又粗又爽又猛毛片免费看| 欧美在线黄色| 高清日韩中文字幕在线| 噜噜噜噜噜久久久久久91| 十八禁人妻一区二区| 欧美色欧美亚洲另类二区| 久久伊人香网站| 亚洲av中文字字幕乱码综合| av天堂中文字幕网| 一区二区三区高清视频在线| 宅男免费午夜| 岛国在线免费视频观看| 在线观看免费午夜福利视频| 欧美极品一区二区三区四区| 真实男女啪啪啪动态图| 日韩欧美在线乱码| 日本黄色视频三级网站网址| av中文乱码字幕在线| 久久国产乱子伦精品免费另类| 亚洲av电影在线进入| 观看美女的网站| 久久久久九九精品影院| 老司机午夜十八禁免费视频| 高潮久久久久久久久久久不卡| 成人精品一区二区免费| 两性午夜刺激爽爽歪歪视频在线观看| 麻豆成人av在线观看| 女人高潮潮喷娇喘18禁视频| 久久久久久久久久黄片| 51午夜福利影视在线观看| 久久午夜亚洲精品久久| 久久久精品大字幕| 成人特级黄色片久久久久久久| 18禁黄网站禁片免费观看直播| 在线观看美女被高潮喷水网站 | 变态另类丝袜制服| 91在线精品国自产拍蜜月 | 男女午夜视频在线观看| 男女做爰动态图高潮gif福利片| 日韩欧美免费精品| 少妇人妻精品综合一区二区 | 久久精品影院6| 久久亚洲真实| 亚洲天堂国产精品一区在线| 国产精品久久久久久久电影 | 国产麻豆成人av免费视频| 国产亚洲精品综合一区在线观看| 激情在线观看视频在线高清| 日韩免费av在线播放| 草草在线视频免费看| 极品教师在线免费播放| 非洲黑人性xxxx精品又粗又长| 国产熟女xx| 老熟妇仑乱视频hdxx| 国产一区二区在线观看日韩 | 久久香蕉国产精品| 中文字幕熟女人妻在线| 一区二区三区免费毛片| 一区二区三区免费毛片| 悠悠久久av| 亚洲精品在线观看二区| 麻豆国产97在线/欧美| 国产淫片久久久久久久久 | 精品日产1卡2卡| av专区在线播放| 精品电影一区二区在线| 88av欧美| 一边摸一边抽搐一进一小说| 一级黄片播放器| 一级毛片女人18水好多| 日韩免费av在线播放| 一二三四社区在线视频社区8| 69av精品久久久久久| 首页视频小说图片口味搜索| 亚洲av中文字字幕乱码综合| 精品一区二区三区人妻视频| 国产高清激情床上av| 制服丝袜大香蕉在线| www.www免费av| 欧美一级毛片孕妇| 亚洲国产日韩欧美精品在线观看 | 欧美在线一区亚洲| 91麻豆av在线| 亚洲成av人片在线播放无| 亚洲国产高清在线一区二区三| 免费电影在线观看免费观看| 制服丝袜大香蕉在线| 免费在线观看亚洲国产| 亚洲中文日韩欧美视频| 村上凉子中文字幕在线| 狠狠狠狠99中文字幕| 久久伊人香网站| 免费av不卡在线播放| 少妇人妻一区二区三区视频| 国产伦人伦偷精品视频| 欧美另类亚洲清纯唯美| 老司机在亚洲福利影院| 午夜福利在线在线| 精品国产三级普通话版| 少妇裸体淫交视频免费看高清| 午夜老司机福利剧场| 国产97色在线日韩免费| 内射极品少妇av片p| 91av网一区二区| 18禁美女被吸乳视频| 成年免费大片在线观看| 国产激情欧美一区二区| 我的老师免费观看完整版| 哪里可以看免费的av片| 亚洲五月婷婷丁香| www.熟女人妻精品国产| 老熟妇乱子伦视频在线观看| 亚洲乱码一区二区免费版| 国产一区二区亚洲精品在线观看| 国产日本99.免费观看| 精品国内亚洲2022精品成人| 欧美三级亚洲精品| 亚洲av美国av| 男人和女人高潮做爰伦理| www.熟女人妻精品国产| 又黄又粗又硬又大视频| 国产蜜桃级精品一区二区三区| 好男人在线观看高清免费视频| 最新美女视频免费是黄的| 国产探花极品一区二区| 18禁黄网站禁片免费观看直播| 岛国在线免费视频观看| 在线播放国产精品三级| 亚洲一区二区三区色噜噜| 国产三级黄色录像| 久久久久性生活片| 日本熟妇午夜| 久久婷婷人人爽人人干人人爱| 久久天躁狠狠躁夜夜2o2o| 国产在视频线在精品| 亚洲一区二区三区色噜噜| 有码 亚洲区| 国产97色在线日韩免费| 日韩精品中文字幕看吧| 国产免费av片在线观看野外av| 一本一本综合久久| 欧美日本视频| www.999成人在线观看| 欧美另类亚洲清纯唯美| 精华霜和精华液先用哪个| bbb黄色大片| 女同久久另类99精品国产91| 99久久无色码亚洲精品果冻| 欧美日韩瑟瑟在线播放| 深爱激情五月婷婷| 99久久99久久久精品蜜桃| 中文字幕人妻丝袜一区二区| 久久人妻av系列| 啦啦啦免费观看视频1| 在线免费观看不下载黄p国产 | 国产熟女xx| 亚洲内射少妇av| 在线播放国产精品三级| 国产伦精品一区二区三区视频9 | 日韩欧美国产一区二区入口| 亚洲成人久久爱视频| 久久精品91蜜桃| 欧美zozozo另类| 搡老熟女国产l中国老女人| 国产黄a三级三级三级人| 高清在线国产一区| 成人国产一区最新在线观看| 岛国在线免费视频观看| 亚洲片人在线观看| or卡值多少钱| 女人被狂操c到高潮| 成人国产综合亚洲| 免费av毛片视频| 麻豆国产av国片精品| 久久精品国产综合久久久| avwww免费| 国产亚洲精品久久久com| 老司机深夜福利视频在线观看| 男女床上黄色一级片免费看| 很黄的视频免费| 手机成人av网站| 欧美黄色片欧美黄色片| 在线免费观看不下载黄p国产 | 一个人观看的视频www高清免费观看| 国产在线精品亚洲第一网站| 亚洲av成人不卡在线观看播放网| av视频在线观看入口| 亚洲国产精品久久男人天堂| 国产主播在线观看一区二区| 欧美中文综合在线视频| 草草在线视频免费看| 99久国产av精品| 国产男靠女视频免费网站| 久99久视频精品免费| 色哟哟哟哟哟哟| 嫁个100分男人电影在线观看| 欧美日韩精品网址| 亚洲黑人精品在线| 日韩国内少妇激情av| 叶爱在线成人免费视频播放| 亚洲久久久久久中文字幕| 亚洲真实伦在线观看| 熟女人妻精品中文字幕| 中文字幕av成人在线电影| 精品欧美国产一区二区三| 日本成人三级电影网站| 18美女黄网站色大片免费观看| 手机成人av网站| 亚洲七黄色美女视频| 欧美日本视频| 香蕉丝袜av| 在线观看舔阴道视频| 90打野战视频偷拍视频| 在线播放无遮挡| 九九热线精品视视频播放| 中文字幕久久专区| 久久香蕉国产精品| 久久久久久久久中文| 久久香蕉国产精品| 午夜免费男女啪啪视频观看 | 在线十欧美十亚洲十日本专区| 叶爱在线成人免费视频播放| 国产精品99久久久久久久久| 一夜夜www| 内射极品少妇av片p| 香蕉丝袜av| 欧美日韩国产亚洲二区| 18禁美女被吸乳视频| 欧美绝顶高潮抽搐喷水| 亚洲av熟女| 老司机深夜福利视频在线观看| 日韩欧美免费精品| 欧美日韩瑟瑟在线播放| 欧美成人a在线观看| 激情在线观看视频在线高清| 色精品久久人妻99蜜桃| 亚洲中文字幕日韩| 一本一本综合久久| 精品电影一区二区在线| 欧美色欧美亚洲另类二区| 久久久久亚洲av毛片大全| 精品无人区乱码1区二区| 国产精品久久久久久亚洲av鲁大| 在线看三级毛片| 午夜福利在线在线| aaaaa片日本免费| 欧美又色又爽又黄视频| 女人被狂操c到高潮| 亚洲国产精品999在线| 成年女人毛片免费观看观看9| 精品久久久久久,| 免费看日本二区| 97人妻精品一区二区三区麻豆| 国产精品亚洲av一区麻豆| 2021天堂中文幕一二区在线观| 国产亚洲欧美98| 亚洲自拍偷在线| 国产精品亚洲av一区麻豆| 日韩成人在线观看一区二区三区| 欧美日韩亚洲国产一区二区在线观看| 国产蜜桃级精品一区二区三区| 尤物成人国产欧美一区二区三区| 深夜精品福利| av专区在线播放| 成人特级av手机在线观看| 国产亚洲精品综合一区在线观看| 午夜亚洲福利在线播放| 国产乱人伦免费视频| 欧美成人一区二区免费高清观看| 亚洲精品在线美女| 女警被强在线播放| 2021天堂中文幕一二区在线观| 亚洲一区高清亚洲精品| av欧美777| 色尼玛亚洲综合影院| 国产蜜桃级精品一区二区三区| 天堂影院成人在线观看| 一二三四社区在线视频社区8| 国产激情欧美一区二区| av在线天堂中文字幕| 两性午夜刺激爽爽歪歪视频在线观看| 母亲3免费完整高清在线观看| 麻豆成人av在线观看| 男人的好看免费观看在线视频| 日本黄大片高清| 在线观看免费午夜福利视频| 日本与韩国留学比较| 亚洲第一欧美日韩一区二区三区| 黄色女人牲交| 欧美一级a爱片免费观看看| www.999成人在线观看| 国产激情偷乱视频一区二区| 俄罗斯特黄特色一大片| 婷婷亚洲欧美| 国产亚洲av嫩草精品影院| 国产高潮美女av| 一进一出抽搐动态| 亚洲国产精品sss在线观看| 99热这里只有精品一区| 亚洲av免费在线观看| 国产激情欧美一区二区| 日韩欧美免费精品| av福利片在线观看| 国产淫片久久久久久久久 | 婷婷六月久久综合丁香| 久久久精品大字幕| 天天添夜夜摸| 国产精品国产高清国产av| 精品久久久久久久末码| 国产极品精品免费视频能看的| 亚洲国产色片| av中文乱码字幕在线| 欧美成人免费av一区二区三区| 久久久成人免费电影| 99热这里只有是精品50| 丰满乱子伦码专区| 欧美绝顶高潮抽搐喷水| 午夜老司机福利剧场| 人妻夜夜爽99麻豆av| 非洲黑人性xxxx精品又粗又长| 国产真实伦视频高清在线观看 | 国产成人欧美在线观看| 久久久成人免费电影| 天堂网av新在线| a级毛片a级免费在线| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲精品久久国产高清桃花| 悠悠久久av| av片东京热男人的天堂| 在线观看一区二区三区| 夜夜躁狠狠躁天天躁| 国产精品野战在线观看| 国产高清视频在线观看网站| 亚洲国产日韩欧美精品在线观看 | 亚洲中文字幕日韩| 中文字幕精品亚洲无线码一区| 97人妻精品一区二区三区麻豆| 91麻豆av在线| 亚洲国产精品999在线| 桃色一区二区三区在线观看| 亚洲 国产 在线| 天天躁日日操中文字幕| 欧美乱码精品一区二区三区| 国产黄a三级三级三级人| 精品人妻一区二区三区麻豆 | 91av网一区二区| 两人在一起打扑克的视频| 99国产极品粉嫩在线观看| 黄色片一级片一级黄色片| 国产免费一级a男人的天堂| 禁无遮挡网站| 久久香蕉精品热| 成人无遮挡网站| 黄色日韩在线| 亚洲黑人精品在线| 国产精品一区二区三区四区免费观看 | 亚洲无线在线观看| www日本在线高清视频| 男人的好看免费观看在线视频| 大型黄色视频在线免费观看| 91麻豆av在线| 国产高清视频在线观看网站| 国产伦一二天堂av在线观看| www.www免费av| 好男人电影高清在线观看| 精品日产1卡2卡| 亚洲精品色激情综合| 欧美日本视频| 午夜老司机福利剧场| 在线观看av片永久免费下载| 色综合欧美亚洲国产小说| 国产精品免费一区二区三区在线| 亚洲成av人片在线播放无| 男女视频在线观看网站免费| 悠悠久久av| 成年女人毛片免费观看观看9| 国产精品免费一区二区三区在线| 国产精品99久久久久久久久| 国产高清有码在线观看视频| 动漫黄色视频在线观看| 最新在线观看一区二区三区| 黄片小视频在线播放| 中文字幕高清在线视频| 久久亚洲真实| 国产乱人视频| 国产成人啪精品午夜网站| 19禁男女啪啪无遮挡网站| 在线十欧美十亚洲十日本专区| 欧美乱妇无乱码| 免费看a级黄色片| 级片在线观看| 听说在线观看完整版免费高清| av福利片在线观看| 精品无人区乱码1区二区| 国产伦精品一区二区三区视频9 | 日韩欧美三级三区| 国产亚洲精品av在线| 中国美女看黄片| 亚洲av成人不卡在线观看播放网| 亚洲国产精品久久男人天堂| 成人特级av手机在线观看| 网址你懂的国产日韩在线| 亚洲成人久久爱视频| av欧美777| 国产伦精品一区二区三区四那| 丝袜美腿在线中文| 精品电影一区二区在线| 日本三级黄在线观看| 久久香蕉国产精品| 黄色日韩在线| 国内精品久久久久久久电影| 他把我摸到了高潮在线观看| 亚洲国产日韩欧美精品在线观看 | 97碰自拍视频| 亚洲第一电影网av| 亚洲性夜色夜夜综合| 熟女电影av网| av片东京热男人的天堂| 久久久久久久久久黄片| 一个人看视频在线观看www免费 | 97人妻精品一区二区三区麻豆| 天堂网av新在线| 狠狠狠狠99中文字幕| 一进一出抽搐gif免费好疼| 亚洲人成网站在线播| 有码 亚洲区| 国产乱人伦免费视频| 午夜激情福利司机影院| 一进一出好大好爽视频| 一本综合久久免费| 一本一本综合久久| 久久精品国产自在天天线| 美女被艹到高潮喷水动态| 欧美成人a在线观看| 神马国产精品三级电影在线观看| 国产av一区在线观看免费| 国产精品爽爽va在线观看网站| 欧美区成人在线视频| 欧美黄色淫秽网站| 亚洲熟妇中文字幕五十中出| 国产精品一区二区三区四区免费观看 | 婷婷精品国产亚洲av| 一级黄片播放器| 免费电影在线观看免费观看| 老司机深夜福利视频在线观看| 久久久国产成人免费| 日本 av在线| 9191精品国产免费久久| 国产亚洲av嫩草精品影院| 18禁黄网站禁片免费观看直播| 国产激情欧美一区二区| 欧美xxxx黑人xx丫x性爽| 欧美日本亚洲视频在线播放| 精品一区二区三区av网在线观看| 亚洲aⅴ乱码一区二区在线播放| 又粗又爽又猛毛片免费看| 中文在线观看免费www的网站| 超碰av人人做人人爽久久 | 中国美女看黄片| 亚洲自拍偷在线| 免费人成在线观看视频色| 日本 欧美在线| 精品人妻偷拍中文字幕| 一本精品99久久精品77| 国产黄色小视频在线观看| 国产爱豆传媒在线观看| 午夜久久久久精精品| 51午夜福利影视在线观看| 亚洲专区中文字幕在线| 国产综合懂色| av视频在线观看入口| 亚洲av成人不卡在线观看播放网| 两个人看的免费小视频| 欧美成人性av电影在线观看| 97碰自拍视频| 精品久久久久久久久久免费视频| 免费人成视频x8x8入口观看| 神马国产精品三级电影在线观看| 精品久久久久久久久久免费视频| 亚洲黑人精品在线| 人人妻人人澡欧美一区二区| 怎么达到女性高潮| 欧美黑人巨大hd| 午夜影院日韩av| av专区在线播放| 桃红色精品国产亚洲av| 一a级毛片在线观看| 国产熟女xx| 精品熟女少妇八av免费久了| 国产精品久久久人人做人人爽|