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

    適用于HCCI燃燒研究的甲苯參比燃料化學(xué)動(dòng)力學(xué)簡(jiǎn)化模型

    2011-12-12 02:43:08張慶峰鄭朝蕾何祖威
    物理化學(xué)學(xué)報(bào) 2011年3期
    關(guān)鍵詞:重慶大學(xué)甲苯品位

    張慶峰 鄭朝蕾 何祖威 王 迎

    (重慶大學(xué)低品位能源利用技術(shù)及系統(tǒng)教育部重點(diǎn)實(shí)驗(yàn)室,重慶400030)

    1 Introduction

    Homogeneous charge compression ignition(HCCI)combustion process holds the promise of both good fuel economy and very low emissions of nitrogen oxides and soot.Therefore, HCCI is considered as a high-efficiency alternative to sparkignited(SI)gasoline operation and as a low-emission alternative to traditional diesel compression ignition(CI)combustion. In a HCCI engine,a homogeneous mixture of fuel and air is compressed sufficiently to raise the mixture temperature above its ignition point.Ignition is distributed throughout the cylinder volume,starting at the hottest regions.As the mixture burns, the expanding gases cause further compression of the unburned mixture and subsequent ignition without flame propagation.More detailed information about progress and recent trends in HCCI engines,as dicussed by Yao et al.1

    The combustion in HCCI engine is controlled by the chemical kinetics.2Hence,the chemical kinetics modeling is one of an important tool to the research of HCCI combustion.The combustion process is initiated by autoignition in HCCI engines.Autoignition depends on the pressure and temperature history of the fuel-air mixture.For the characterization of fuels the autoignition delay is a significant observable and is a critical parameter for HCCI engines.Tanaka et al.3performed experiments on HCCI combustion in a rapid compression machine with complex fuels,such as cyclic paraffins,olefins,and aromatics,which exist in petroleum-based fuels.They concluded that ignition delay time and burn rate could be independently controlled using various fuel mixtures and additives for HCCI combustion.

    Gasoline is one of the most important fuels in internal combustion engine,and it is a complex mixture of hundreds of hydrocarbons,such as paraffins,olefins,aromatics,and cycloalkanes.Unfortunately,it is not currently possible to represent the complex chemistry of full blend of gasoline in a detailed chemical kinetics model.In general,the term surrogate denotes a simpler representation of a fully-blended fuel.The use of iso-octane is the simplest gaoline surrogate.Binary blends of iso-octane and n-heptane,pimary reference fuel(PRF)is applied widespread as gasoline surrogates for variable octane number fuel.With rapid growth of the HCCI combustion research,the study on chemical kinetics mechanism of iso-octane and PRF become very active.However,it is clear that the values of octane number are insufficient in characterzing autoignition behavior under HCCI combustion conditions.Due to differing physical and chemical properties,the behavior of PRF mixtures in homogenous charge and spark ignition engines is significantly different than that of gasoline itself.Consequently,more complex mixtures of individual species representing different molecular groups have been suggested as surrogate fuels for gasoline.As explained above,gasoline contains significant quantities of aromatic hydrocarbons.Toluene is the largest single aromatic component of the majority of fuels and the mole fraction is as high as 35%in premium fuels, and toluene is formed during the oxidation of other hydrocarbons and commercial fuels.4Therefore,toluene should be concerned.The relative importance of fuel components for gasoline,diesel and jet fuel surrogates were reported by Pitz et al.,5the result shows that a ternary mixture of iso-octane,n-heptane and toluene is suitable as a gasoline surrogate fuel for the HCCI engine simulation.This ternary blend is usually called toluene reference fuel(TRF).

    The prior published mechanisms for iso-octane and n-heptane developed by Curran et al.6,7were minimized,optimized and combined to reproduce PRF results.Chaos et al.8constructed a detailed chemical kinetic mechanism for the TRF oxidation consisting of 469 species and 1221 reactions,and they performed the validated experiment in a variable pressure flow reactor at 1.27 MPa.Andrae et al.9added the toluene sub-mechanism to the PRF detailed chemical kinetic mechanism and developed a detailed chemical kinetics model(1083 species and 4635 reactions)for the autoignition of TRF oxidation,which was validated using shock tube autoignition delay time data. Afterwards,they presented a semidetailed mechanism10(137 species and 633 reactions)in a HCCI engine on the autoignition of TRF.Sakai et al.11mechanisms for iso-octane and n-heptane were added to a detailed toluene submechansim.The model showed generally good agreement with ignition delay times measured in a shock tube and a rapid compression machine. Sakai et al.proposed a toluene sub-mechanism started from the toluene mechanism developed by Pitz et al.,12which was merged to the previous PRF model.The TRF oxidation mechanism(783 species and 2883 reactions)was validated by existing shock tube and flow tube data.Anderlohr et al.13proposed a TRF oxidation mechanism(536 species and 3000 reactions) which was based on the PRF autoignition model coupled with the model for the oxidation of toluene.Flow rate and sensitivity analysis were performed in order to explain the low temperature chemical kinetics,especially the impact of NOxon hydrocarbon oxidation.

    At present,most of TRF mechanisms are the detailed chemical kinetic model,and the development of TRF oxidation mechanisms is still in primary stage.The detailed chemical kinetic model is large and difficult to use in engine models.It spends too much money and time to apply the 3D-CFD combustion simulation with the detailed mechanisms.In sum,the research on the TRF oxidation mechanism is not enough,specially the reduced mechansim.The necessity of developing the reduced TRF mechanism is brought out.In this study,a new reduced mechanism is suitable for generalized numerical simulation which developed and validated by different experiments related to HCCI combustion.

    2 Construction of TRF kinetic model

    Because this mechanism is mainly for HCCI engine application,it is used over the range of low equivalence ratio and high pressure.Since the ignition of HCCI engine is very sensitive to temperature,the mechanism should include both low temperature and high temperature kinetics.

    The reaction mechanism steps of iso-octane and n-heptane oxidation are similar.However,the rate coefficients for many of the individual reactions are different reflecting the very different structure of these two molecules.In following reactions RH denotes paraffins,and Q denotes olefins structure.A chemical kinetic mechanism proposed by Tanaka et al.3is used as a starting model,here the low temperature mechanism includes the following reactions.

    The mechanism in the case of iso-octane is as follows.Fuel oxidation is initiated to form the alkyl radical R and HO2when H atoms are abstracted by oxygen molecules RH+O2?R+HO2(Reaction(R1,R14)).At low temperature,in the sequence of reactions:an alkyl radical is oxidized by an oxygen molecule R+O2?RO2(the first oxygen addition process,Reaction(R2, R15)),an internal hydrogen is rearranged RO2?QOOH(Reaction(R3,R16)),another oxygen molecule is added to a hydroperoxy alkyl radical QOOH+OB2?OOQOOH(the second oxygen addition process,Reaction(R4,R17)),and the product decomposes irreversibly OOQOOH?OH+OROOH(Reaction (R5,R18)),as a result,an OH radical and an ketohydroperoxide are formed.After the reaction of an OH radical with a fuel molecule RH+OH?R+H2O(Reaction(R6,R19)),highly exothermic chain-cycle reactions from reaction(R2)to(R6)are subsequently completed.Then the branching reaction started, lots of heat is released and the temperature rises rapidly,until the following competing reactions R+O2?Q+HOR2(R7,R20) becomes faster than reaction(R2).

    Tanaka et al.3included three global reactions and directly decomposed OC8H15O and OC7H13O into CO and H2O.Because the presence of the important intermediate product is ignored, it is difficult for Tanaka model to predict HC and CO emissions.At the same time,their sensitivity analysis shows that these reaction paths are relatively unimportant to ignition delay times and heat release ratio using the species OC8H15O and OC7H13O.Tsurushima14pointed out that CH2O and HO2were both key intermediate species,and the fuel consumption rate during the low-temperature chain reaction was too small with Tanaka model.In order to raise prediction level,the following several reactions are included and the large molecules formed during the low temperature stage accordingly break up into smaller molecules.In these reactions,R′denotes paraffins with C5or C6.The decomposition process of ketohydroperoxide,the reaction OROOH?R′CO+CH2O+OH(R8,R21)and R′CO+O2?small molecule+CO+HO2(R9,R22)were modified.The reaction Q+O2?R′+CH2+HCO(R10,R23)was added to condider CH2O formation.Then the reaction RH+HO2?R+H2O2(R11,R24)was added to maintain the balance of production and consumption of alkyl radical.In this model,betascission of alkyl radicals is considered in the reaction R?R′+ small molecules(R12,R25).Afterwards,smaller alkyl radical R′decomposes to small molecules through the reaction R′?small molecules(R13,R26).

    Fig.1 Major reaction branches of toluene oxidation

    From the existing mechanisms of the toluene,similarity in reaction pathways can be noted.The overall reaction proceeds primarily through C6H5CH3→C6H5CH2→C6H5CHO→C6H5CO→C6H5→C6H5O→C6H5OH,followed by the ring-breaking reactions.Fig.1 presents the chemical reaction paths of toluene. The initial mechanism of the toluene oxidation is from Sivara-makrishnan et al.4Since benzene is a key intermediate in the oxidation of toluene,the benzene submodel should be considered.Selected benzene reactions were mainly based on Andrae et al.10from Alzueta model,15which was validated by comparison with the experimental data as well as flow reactor and mixed reactor data from literature.The validated experiments were also performed under plug-flow conditions,at excess air ratios ranging from close to stoichiometric to very lean,which are relevant to HCCI combustion.Some additional reactions including the species CH2OH,CH2CO were added,based on the work by UCSD.16The toluene sub-mechanism is shown in Appendix(available free of charge via the internet at http://www. whxb.pku.edu.cn).

    Some rate constantsare updated mainly fortoluene (C6H5CH3)and benzyl radical(C6H5CH2)reactions taken from recent studies.The rate coefficients of the toluene reactions are derived from published studies of Oehlschlaeger et al.17-19and Seta et al.20Oehlschlaeger et al.investigated some reactions of toluene using ultraviolet laser absorption of benzyl radicals in shock tube experiments.The rate coefficient results were in excellent agreement with the previous rate determinations.They thought the reaction of toluene with H-atoms producing benzyl radicals and molecular hydrogen is an important elementary reaction in the toluene pyrolysis and oxidation reaction systems. The reaction retards ignition during toluene oxidation by removing an H-atom that otherwise would lead to chain branching via reaction with O2.The decomposition of toluene is an important initiation step in the oxidation of toluene at moderateto high-temperatures(>1200 K),and the toluene decomposition takes place via two channels:C6H5CH3?C6H5CH2+H (R27)and C6H5CH3?C6H5+CH3(R28).The reaction of toluene with molecular oxygen is also the primary initiation reaction during the oxidation of toluene,which has been studied.Seta et al.measured the overall rate constants of the toluene reaction with OH radical and estimated channel specific rate coefficients of this reaction on the basis of RRKM-master equation analysis.This reaction is seen to account for the major consumption of toluene during the induction period.

    Benzyl is a resonantly stabilized radical that commonly occurs as an intermediate in the oxidation and pyrolysis of toluene.Benzyl is generated as an intermediate at significant concentrations during the oxidation of toluene via the reaction of toluene with molecular oxygen,the reaction of toluene with radical species,and the toluene thermal decomposition channel leading to benzyl radicals and H-atoms.Silva et al.21pointed out that the bimolecular reactions of benzyl with HO2were important in the oxidation of toluene,especially at low to moderate temperatures.We also updated rate constants of the benzyl radical reactions according to the studies of Andrae et al.10In simulation for the experiment of Davidson et al.,22the reaction C6H5CH2OO?C6H5CHO+OH(R48)was found to be very sensitive to the overall reactivity for neat toluene-air mixtures and the rate constant was tuned for best fit of model to the shock tube data that proposed by Andrae et al.,10and Klotz et al.23pointed out that considering the formation of 1,3-butadiene from toluene was necessary and the addition of iC4H5reactions significantly improved predictions of 1,3-butadiene and acetylene.Seven reactions about iC4H5were added into the toluene submechanism.Andrae et al.10performed simulations with and without the reactions to checked the importance of cross reactions.The introduction of cross reactions into a detailed model for TRF was found to improve the results in those particular HCCI simulations.In this work,five cross reactions (R147-R151)were merged to the TRF mechanism based on previous work of Andrae.Cross reactions between phenyl radicals and the primary reference fuels are also of importance, which were considered in the new model.

    Unlike previous TRF mechanisms incorporating detailedmechanism,in order to keep the mechanism size compact,thereactions are taken from a reduced mechanism for HCCI engines developed by Patel et al.24.The main branching reaction is H2O2+M?OH+OH+M(R173)at the temperatures from 900 to 1100 K,and the reaction H+O2+M?O+OH+ M(R174)is the significant branching reaction at higher temperatures,when H radical becomes an important species.It may be noted,the role of the competition between this branching reaction(R174)and the recombination process reaction H+ O2+M?HO2+M(R175)cannot be disregarded.Reaction (R174)is added to the reactions of C1-C3oxidation developed by Patel.Then,the whole toluene reference fuels oxidation mechanism is completed.The final version of the TRF mechanism consists of 70 species and 196 reactions.

    3 Model validation

    3.1 Validation of ignition delay

    A parameter often used in autoignition studies is the ignition delay,which is the time required for autoignition to happen once the fuel-air mixture is raised to a given pressure and temperature and held at that condition in a shock tube or a rapid compression machine.We have compared the model predictions with some representative experimental results.All calculations were performed by using CHEMKIN 4.1 program package,25and an adiabatic homogeneous reactor is assumed for simulations of ignition delay times.All the figures of ignition delay are shown in semi-log plot.In terms of detailed studies, the whole validation can be divided into three categories:primary reference fuels,toluene,and toluene reference fuels.

    3.1.1 Primary reference fuels

    Mechanism validation of primary reference fuels were performed for the ignition delay times from the shock tube experiments in previous studies,26,27The ignition delay time is defined in the experiment as the time interval between the arrival of the reflected shockwave and the onset of ignition in combustion chamber.And computation was done under completely homogeneous,constant volume,and adiabatic conditions.

    Fig.2(a)is the comparison of ignition delays between the cal-culation and the shock tube experiments of Fieweger et al.26at initial pressures of 1.3,3.4,4.5 MPa with stoichiometric iso-octane-air mixture.At the higher pressures,the ignition delay time becomes shorter,but the slope of the curves is essentially unchanged.Fig.2(b)is computed ignition delays of iso-octane for the different equivalence ratios at a pressure of 4 MPa also compared to experimental data.The ignition delay times increase with decreasing equivalence ratio.The values of ignition delay show a pronounced negative temperature coefficient (NTC)behaviour for stoichiometric mixtures and rich mixtures.The ignition delay with the new model seems good agreements with the experiments.In Fig.3 a comparison of ignition delay times for stoichiometric fuel/air mixture with different octane numbers at an initial pressure of 4 MPa with comparisons to the shock tube data of Fieweger.For temperatures above 1000 K,there is no significant difference among the ignition delay times of the mixtures.Overall,however,the ignition timing advances with the decrease of the octane number.In this graph,the sensitivity of ignition delay to the octane number shows excellent agreement to the shock tube data.

    Fig.2 Comparison of experimental26and computed ignition delay of iso-octane-air mixture at different pressures(a)and equivalence ratios(b)

    Fig.3 Comparison of experimental26and calculated ignition delay of PRF-air mixture

    Fig.4 Comparison of experimental27and calculated ignition delay of n-heptane-air mixture at different pressures(a)and equivalence ratios(b)

    Fig.4(a)shows computed ignition delays of stoichiometric n-heptane/air mixture at initial pressures of 0.65,1.35,4.2 MPa,with comparison to shock tube data of Ciezki et al.27Fig.4 (b)gives computed ignition delays of n-heptane at a pressure of 4 MPa and equivalence ratios of 0.5,1.0,and 2.0.In these figures,it can be seen that the ignition delays of PRF with new model shows very good agreement with the experimental data. The calculation is sensitive to changes in temperature and mixture composition and captures the NTC region.

    3.1.2 Toluene

    Fig.5(a)shows the comparison between the calculated results and the experimental results measured by Davidson et al.22in shock tube for toluene.The ignition delay time is defined in the experiment as the time interval between the arrival of the reflected shockwave and the onset of ignition at the sidewall observation location.The arrival of the reflected shockwave was determined by the step rise in pressure,and the onset of ignition was determined by monitoring either the pressure history or the emitted light corresponding to an intermediate species.The onset of ignition from the pressure history as well as both CH*and OH*emission were defined by locating the time of steepest rise and linearly extrapolating back in time to the pre-ignition baseline.As can be seen from Fig.5(a),good agreement between experimental and calculated data is obtained at a pressure of 5 MPa,especially for the conditions in relative high temperature regime.But the calculation tends to predict longer ignition delays than the experimental data in relative low temperature regimes.The calculation over-predicts the ignition delay times below 1250 K for φ=0.5 and below 1080 K for φ=1.0.

    Fig.5 Comparison of experimental[Davidson et al.22(a),Pitz et al.12(b)]and calculated ignition delay of toluene-air mixture

    In Pitz's shock tube experiment,1210%of the maximum OH emission was used as an indication of ignition.The toluene concentration in the reactants kept constant at 1.25%,and the values of equivalence ratio were 0.5,1.0,and 1.5.The reflected shock pressures and temperatures ranged from 8.0 to 9.4 MPa and 1300 to 1900 K,respectively.Using the new reduced TRF mechanism,simulation was carried out in accordance with Pitz's experimental conditions.Fig.5(b)presents the comparison of the computed ignition delays and the measuresed ones for the different equivalence ratios.As a whole,the predictions are satisfactory,and the trends of the calculated curves are in agreement with the experiment.

    3.1.3 Toluene reference fuels

    In Gauthier's shock tube study,28the two ternary surrogates used in the ignition delay time experiments consisted of iso-octane,toluene,and n-heptane in the following proportions:surrogateAand surrogate B.Fuels and their properties are presented in Table 1.Ignition delay times of toluene reference fuels were measured at temperatures 850 to 1200 K.The definition of ignition delay time is the same as Davidson et al.22(see Section 3.1.2).Fig.6 shows a comparison of ignition delay times for stoichiometric fuel-air mixture with comparisons to the shock tube data.At this temperature range and higher pressure conditions,the variation range of TRF ignition delay time is narrower than that of PRF in Fig.3.The computed ignition delay with new model shows good agreement with the experimental data.

    3.2 Validation in HCCI engines

    A single zone model is used in the computation assuming that temperature,pressure and species concentrations are uniform in the chamber.The wall heat transfer,residual gasesfrom previous cycle,and the temperature inhomogeneity in the boundary and crevice zones are ignored in the model.Despite its simplicity,the single zone model is a useful tool for investigating certain fundamental aspects of HCCI combustion.The operating conditions of each engine experiment in the simulation are given in Table 2.

    Table 1 Compostion(liquid volume fraction(%))of TRF investigated

    Fig.6 Comparison of experimental28and calculated ignition delay of TRF-air mixture(a)surrogateA,(b)surrogate B;er=equivalence ratio

    Table 2 Engine operating conditions

    Dec et al.29have performed a series of experiments using iso-octane as fuel in HCCI engine.Fig.7 shows the cylinder pressure profiles of modeling and experiment for various equivalence ratios at the compression ratio of 18,engine speed of 1200 r·min-1,and initial pressure of 0.12 MPa.Other operating parameters were held constant.The highest pressure increases with the increases of fuel load.Because of the assumption of adiabatic condition in the single zone model,the maximum pressures are over predicted,but the trends of the calculation results are in accordance with experiment results for three loads.

    HCCI combustion data30for PRF fuel-air mixtures tested in a single-cylinder CFR research engine were used to validate the predictions of PRF autoignition under engine conditions. The compression ratio and the engine speed were set at 16.55, and 900 r·min-1respectively for all the cases.The intake pressure was maintained at 0.1 MPa.The results of the fuel/air mixture(surrogate C and surrogate D)HCCI experiments are shown in Fig.8,and fuels and their properties are presented in Table 1.The cylinder pressure histories for surrogate C at the initial temperature of 320 K and the equivalence ratio of 0.25 are shown in Fig.8(a).Surrogate D at the initial temperature of 340 K and the equivalence ratio of 0.29 shown in Fig.8(b).Curran et al.6,7developed detailed reaction mechanism for iso-octane and n-heptane,which were widespreadly applied.The reduced mechanism and detailed mechanism are compared by computation.To be noted,the pressures predicted by the present reduced mechanism are almost the same as that by the PRF detailed mechanism.However,under the computing environment of 2.51 GHz CPU and 2 GB DDR,it only took 2 s on average with the reduced mechanism,which is far less than the time 420 seconds with the detailed mechanism.

    Fig.7 Cylinder pressure profiles for the modeling and the experiment29with iso-octane at different aquivalence ratios

    HCCI experiments have been conducted with TRF in the HCCI engine by Andrae et al.10Fig.9 shows both the experimental results and the simulation results with surrogate B(see Table 1)at the engine speed of 1200 r·min-1and the compression ratio of 14.04.The start of calculations is at 99 degree before top dead center.Two different initial conditions are set: the first condition is at the initial temperature of 489 K and initial pressure of 0.15 MPa(Fig.9(a));the second condition is at the initial temperature of 435 K and initial pressure of 0.33 MPa(Fig.9(b)).The results show good agreements both at the two conditions.Fig.9 also shows heat production predicted by the model.Maximum heat production of two operating conditions is at the crank angle of-1.8 and 1.8 degree after top dead center,respectively.The ignition delay time is more advanced at condition 1(high temperature and low pressure)compared with condition 2(low temperature and high pressure).However,maximum heat production of condition 2 is higher.

    Fig.8 Cylinder pressure profiles for the modeling and the experiment with PRF(a)surrogate C,(b)surrogate D

    Fig.9 Cylinder pressure profiles for the simulation and the experiment10with TRFcondition:(a)high temperature and low pressure;(b)low temperature and high pressure

    Fig.10 Normalized sensitivity coefficients for TRF at different conditionscondition:(a)high temperature and low pressure;(b)low temperature and high pressure

    Fig.10 presents a bar plot of normalized sensitivity coefficients for surrogate B.The sensitivity calculation was done for the 10 most sensitive reactions at the crank angle for maximum heat production.The conditions for the calculations are conditions 1 and 2 as in Fig.9.For condition 1(see Fig.9(a)),the reaction with two of the highest positive sensitivity are alkyl radical isomerization in C7H15OO?C7H14OOH(R16)and the toluene reaction with O radical in C6H5CH3+O=C6H5CH2+OH (R32).

    For the case of condition 2(see Fig.9(b)),on the other hand, reactions(R16)and(R32)are not among the 10 most sensitive reactions.Instead,the reaction(R80)of phenol radicals with O2has the highest positive sensitivity.The sensitivity results for higher pressure are the same as shock tube experiments performed by Sivaramakrishnan et al.4

    For both operating conditions,the following five reactions show higher sensitivity.

    Although there are some small variations,the trend in sensitivity to reactions is the same with changing operating conditions.Reaction(R173)initializes the thermal explosion,and the reaction of fuel oxidation with the product OH radicals happens immediately,which causes the temperature rise rapidly. Fig.10 shows that the reactions(R4)and(R8)have high positive sensitivity,especially in the condition of low temperature and high pressure.The reaction(R178)of formaldehyde (CH2O)with OH radicals shows a higher negative sensitivity, and the reaction(R179)of CH2O with HO2radicals also shows a higher positive sensitivity.Formaldehyde is very important intermediate,which should not be neglected in constructing mechanism.

    4 Conclusions

    (1)A reduced toluene reference fuel(TRF)oxidation mechanism for HCCI combustion comprised of 70 species and 196 reactions was developed.The mechanism has been validated by comparing computed results with experimental data within the HCCI combustion mode.

    (2)In shock tube conditions,computed ignition delay time with the new TRF mechanism was validated with the experimental data in literature with iso-octane,n-heptane,and toluene as fuels at different conditions.Good agreements of ignition delay times are obtained.

    (3)In HCCI engine conditions,the comparison of the computed and measured results with PRF/TRF as fuels indicates that the combustion process can be well predicted with this mechanism.A sensitivity analysis from low intake pressureor high intake temperature to high intake pressure or low intake temperature shows that reaction(R80)of phenol radicals with O2has the higher sensitivity as the pressure rising,and formaldehyde(CH2O)is very important intermediate,which should not be neglected.

    (4)As a conclusion,the new TRF model is believed to have enough accuracy for the HCCI combustion simulation.It is suitable for toluene/PRF/TRF as fuels.And the present reduced model can be a useful tool to study HCCI engine operation.

    (1)Yao,M.F.;Zheng,Z.L.;Liu,H.F.Prog.Energ.Combust. 2009,35,398.

    (2)Westbrook,C.K.Proc.Combust.Inst.2000,28,1563.

    (3)Tanaka,S.;Ayala,F.;Keck,J.C.Combust.Flame 2003,134, 219.

    (4) Sivaramakrishnan,R.;Tranter,R.S.;Brezinsky,K.Proc. Combust.Inst.2005,30,1165.

    (5) Pitz,W.J.;Cernansky,N.P.;Dryer,F.L.;Egolfopoulos,F.N.; Farrell,J.T.;Friend,D.G.;Pitsch,H.SAE Tech.Pap.Ser.2007, 2007-01-0175.

    (6) Curran,H.J.;Gaffuri,P.;Pitz,W.J.;Westbrook,C.K.Combust. Flame 1998,114,149.

    (7) Curran,H.J.;Gaffuri,P.;Pitz,W.J.;Westbrook,C.K.Combust. Flame 2002,129,253.

    (8) Chaos,M.;Zhao,Z.;Kazakov,A.;Gokulakrishnan,P.; Angioletti,M.;Dryer,F.L.APRF+Toluene Surrogate Fuel Model for Simulating Gasoline Kinetics.In 5th U.S. Combustion Meeting,March 25-28,2007;University of California,San Diego,California,Paper#E26.

    (9)Andrae,J.C.G.;Bj?rnbom,P.;Cracknell,R.F.;Kalghatgi,G. T.Combust.Flame 2007,149,2.

    (10)Andrae,J.C.G.;Brinck,T.;Kalghatgi,G.T.Combust.Flame 2008,155,696.

    (11) Sakai,Y.;Miyoshi,A.;Koshi,M.;Pitz,W.J.Proc.Combust. Inst.2009,32,411.

    (12) Pitz,W.J.;Seiser,R.;Bozzelli,J.W.;Seshadri,K.;Chen,C.J.; Costa,D.;Fournet,R.;Billaud,F.;Battin-Leclerc,F.; Weatbrook,C.K.Chemical Kinetic Study of Toluene Oxidation.In 29th International Symposium on Combustion, Hokkaido University,Sapporo,Japan,July 21-26,2002; Elsevier:New York,2002,UCRL-JC-125890.

    (13)Anderlohr,J.M.;Bounaceur,R.;Pires Da Cruz,A.; Battin-Leclerc,F.Combust.Flame 2009,156,505.

    (14) Tsurushima,T.Proc.Combust.Inst.2009,32,2835.

    (15)Alzueta,M.U.;Glarborg,P.;Dam-Johansen,K.Int.J.Chem. Kinet.2000,32,498.

    (16)Welcome to Chemical-Kinetic Mechanisms for Combustion Applications.http://maeweb.ucsd.edu/~combustion/cermech/ index.html(accessedAug 27,2010).

    (17) Oehlschlaeger,M.A.;Davidson,D.F.;Hanson,R.K.Proc. Combust.Inst.2007,31,211.

    (18) Oehlschlaeger,M.A.;Davidson,D.F.;Hanson,R.K.J.Phys. Chem.A 2006,110,9867.

    (19)Oehlschlaeger,M.A.;Davidson,D.F.;Hanson,R.K.Combust. Flame 2006,147,195.

    (20) Seta,T.;Nakajima,M.;Miyoshi,J.Phys.Chem.A 2006,110, 5081.

    (21) Silva,G.;Bozzelli,J.W.Proc.Combust.Inst.2009,32,287.

    (22)Davidson,D.F.;Gauthier,B.M.;Hanson,R.K.Proc.Combust. Inst.2005,30,1175.

    (23) Klotz,S.D.;Brezinsky,K.;Glassman,I.Modeling the Combustion of Toluene-Butane Blends.In 27th Symposium (International)on Combustion,University of Colorado at boulder,USA;The Combustion Institute:Pittsburgh,PA,1998; 337.

    (24) Patel,A.;Kong,S.C.;Reitz,R.D.SAE Tech.Pap.Ser.2004, 2004-01-0558.

    (25)Kee,R.J.;Rupley,F.M.;Miller,J.A.;et al.CHEMKIN Release 4.1,Reaction Design:San Diego,CA,2006.

    (26)Fieweger,K.;Blumenthal,K.R.;Adomeit,G.Combust.Flame 1997,109,599.

    (27)Ciezki,H.K.;Adomeit,G.Combust.Flame 1993,93,421.

    (28) Gauthier,B.M.;Davidson,D.F.;Hanson,R.K.Combust. Flame 2004,139,300.

    (29) Dec,J.E.;Sj?berg,M.SAE Tech.Pap.Ser.2003,2003-01-0752.

    (30)Aroonsrisopon,T.;Sohm,V.;Werner,P.;Foster,D.E.; Morikawa,T.;Lida,M.SAE Tech.Pap.Ser.2002, 2002-01-2830.

    猜你喜歡
    重慶大學(xué)甲苯品位
    品位(外一首)
    重慶大學(xué)學(xué)報(bào)征稿簡(jiǎn)則
    鐘情山花爛漫 品位幸福時(shí)光
    高效液相色譜法測(cè)定降糖藥甲苯磺丁脲片中甲苯磺丁脲的含量
    金堆城鉬礦床硫元素分布規(guī)律研究
    1-(對(duì)甲苯基)-2-(三對(duì)甲苯基-5-亞磷?;?乙醛的汞(Ⅱ)配合物的X射線晶體學(xué)、光譜表征和理論計(jì)算研究
    Who Is The Master?
    大東方(2018年9期)2018-10-21 15:29:02
    “精益管理五原則”在高校圖書館社區(qū)服務(wù)中的應(yīng)用——以重慶大學(xué)城為例
    萊鋼3200 m3高爐低品位冶煉實(shí)踐
    山東冶金(2015年5期)2015-12-10 03:27:38
    甲苯-4-磺酸催化高效合成尼泊金正丁酯防腐劑
    日韩欧美精品免费久久| 日韩欧美三级三区| 国产 一区精品| 波多野结衣高清无吗| 免费观看在线日韩| 九九在线视频观看精品| 不卡视频在线观看欧美| 国产亚洲精品av在线| kizo精华| 看片在线看免费视频| 久久人妻av系列| 成人美女网站在线观看视频| 成人美女网站在线观看视频| 级片在线观看| 亚洲天堂国产精品一区在线| 亚洲精品乱码久久久v下载方式| 亚洲成a人片在线一区二区| 久久人人爽人人片av| 成人二区视频| 精品99又大又爽又粗少妇毛片| 国产精品伦人一区二区| 国产中年淑女户外野战色| 国产麻豆成人av免费视频| 三级毛片av免费| 中文字幕精品亚洲无线码一区| 狂野欧美白嫩少妇大欣赏| 国产视频内射| 99热这里只有是精品在线观看| 男人狂女人下面高潮的视频| 校园春色视频在线观看| 麻豆av噜噜一区二区三区| 韩国av在线不卡| 女人被狂操c到高潮| 成人鲁丝片一二三区免费| 国产精品伦人一区二区| 久久久国产成人精品二区| 国产视频内射| 午夜激情欧美在线| 中出人妻视频一区二区| 欧美成人免费av一区二区三区| 久久精品国产亚洲av天美| 亚洲不卡免费看| 99久久成人亚洲精品观看| 老师上课跳d突然被开到最大视频| 国产 一区 欧美 日韩| 级片在线观看| 久久这里有精品视频免费| 赤兔流量卡办理| 岛国毛片在线播放| 人人妻人人澡人人爽人人夜夜 | 在线a可以看的网站| 夜夜看夜夜爽夜夜摸| 国产三级中文精品| 欧美日韩精品成人综合77777| 26uuu在线亚洲综合色| 亚州av有码| 国产一区二区三区在线臀色熟女| 最近中文字幕高清免费大全6| 一个人观看的视频www高清免费观看| 18禁在线无遮挡免费观看视频| 亚洲最大成人手机在线| 亚洲四区av| 亚洲人与动物交配视频| 干丝袜人妻中文字幕| 欧美成人精品欧美一级黄| 欧美极品一区二区三区四区| 老师上课跳d突然被开到最大视频| 日韩一区二区视频免费看| 中文字幕av在线有码专区| 国产精品野战在线观看| 亚洲欧美成人精品一区二区| 搞女人的毛片| 亚洲,欧美,日韩| 岛国毛片在线播放| 在线a可以看的网站| 国产高潮美女av| 久久精品国产亚洲av涩爱 | 欧美性猛交╳xxx乱大交人| 日本爱情动作片www.在线观看| 精品久久久久久久末码| 精品欧美国产一区二区三| 在线免费十八禁| 麻豆av噜噜一区二区三区| 久久久久久久久久久丰满| 日韩大尺度精品在线看网址| kizo精华| 网址你懂的国产日韩在线| 永久网站在线| 男人舔奶头视频| 看非洲黑人一级黄片| 久久精品国产清高在天天线| 九色成人免费人妻av| 狂野欧美激情性xxxx在线观看| 免费大片18禁| 一边摸一边抽搐一进一小说| 免费在线观看成人毛片| 嘟嘟电影网在线观看| 久久草成人影院| 国产av一区在线观看免费| 又爽又黄a免费视频| 看片在线看免费视频| 欧美不卡视频在线免费观看| av.在线天堂| 欧美+亚洲+日韩+国产| 成人国产麻豆网| 国产高清视频在线观看网站| 听说在线观看完整版免费高清| 欧美成人a在线观看| 久99久视频精品免费| 成年免费大片在线观看| 欧洲精品卡2卡3卡4卡5卡区| 国产综合懂色| 日韩在线高清观看一区二区三区| 看十八女毛片水多多多| 人妻夜夜爽99麻豆av| 国产一区二区在线观看日韩| 亚州av有码| 我要搜黄色片| 男插女下体视频免费在线播放| 秋霞在线观看毛片| 亚洲精品乱码久久久v下载方式| 美女被艹到高潮喷水动态| 欧美人与善性xxx| 亚洲中文字幕日韩| 天堂中文最新版在线下载 | 老熟妇乱子伦视频在线观看| 人妻夜夜爽99麻豆av| 97人妻精品一区二区三区麻豆| av免费在线看不卡| 久久综合国产亚洲精品| 久久午夜亚洲精品久久| 国产真实乱freesex| 亚洲婷婷狠狠爱综合网| 久久人人爽人人片av| 99久久成人亚洲精品观看| 又爽又黄无遮挡网站| 真实男女啪啪啪动态图| 免费黄网站久久成人精品| 欧美高清性xxxxhd video| a级一级毛片免费在线观看| 黄色日韩在线| 亚洲激情五月婷婷啪啪| 久久精品夜夜夜夜夜久久蜜豆| 国产黄a三级三级三级人| 久久6这里有精品| 精品人妻熟女av久视频| 99热这里只有精品一区| 中文字幕av在线有码专区| 欧美日韩综合久久久久久| 日本与韩国留学比较| 亚洲av成人av| 亚洲欧美日韩卡通动漫| 久99久视频精品免费| 欧美成人精品欧美一级黄| 国产欧美日韩精品一区二区| 白带黄色成豆腐渣| 大型黄色视频在线免费观看| 久久久成人免费电影| 激情 狠狠 欧美| 亚洲婷婷狠狠爱综合网| 成人性生交大片免费视频hd| 午夜精品国产一区二区电影 | 真实男女啪啪啪动态图| 欧美激情国产日韩精品一区| 三级国产精品欧美在线观看| 国产蜜桃级精品一区二区三区| 久久久久久久午夜电影| 少妇猛男粗大的猛烈进出视频 | 青春草视频在线免费观看| 亚洲色图av天堂| 天天躁夜夜躁狠狠久久av| 深夜精品福利| 夜夜爽天天搞| 女人被狂操c到高潮| 天堂av国产一区二区熟女人妻| 欧美区成人在线视频| 91aial.com中文字幕在线观看| 长腿黑丝高跟| 欧美成人一区二区免费高清观看| 嘟嘟电影网在线观看| 91精品一卡2卡3卡4卡| av在线蜜桃| 丰满人妻一区二区三区视频av| 69av精品久久久久久| 国产黄色视频一区二区在线观看 | 精品久久久久久久久久免费视频| 国产免费男女视频| 亚洲成人久久性| 久久亚洲精品不卡| 成人漫画全彩无遮挡| 99久国产av精品国产电影| 特大巨黑吊av在线直播| 国产欧美日韩精品一区二区| 全区人妻精品视频| eeuss影院久久| av在线天堂中文字幕| 99在线视频只有这里精品首页| 欧美成人免费av一区二区三区| 国产激情偷乱视频一区二区| 国语自产精品视频在线第100页| 亚洲欧美精品专区久久| 99在线人妻在线中文字幕| 欧美成人一区二区免费高清观看| 亚洲不卡免费看| 99久久人妻综合| 亚洲成av人片在线播放无| 悠悠久久av| 日韩欧美三级三区| 日日摸夜夜添夜夜爱| 亚洲一区高清亚洲精品| 国产精品野战在线观看| av.在线天堂| 亚洲欧洲国产日韩| 天美传媒精品一区二区| 久久精品久久久久久久性| 亚洲成av人片在线播放无| 天天躁夜夜躁狠狠久久av| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 神马国产精品三级电影在线观看| 久久中文看片网| 能在线免费看毛片的网站| 国产精品久久久久久久久免| 能在线免费观看的黄片| 国产色婷婷99| 亚洲中文字幕一区二区三区有码在线看| 日日撸夜夜添| 18禁裸乳无遮挡免费网站照片| 国产精品永久免费网站| 黄色配什么色好看| 深夜a级毛片| 精品免费久久久久久久清纯| 欧美3d第一页| 国产亚洲精品av在线| 91久久精品电影网| 少妇的逼好多水| 麻豆国产av国片精品| 男女那种视频在线观看| 一级黄色大片毛片| 久久久久久久久大av| 欧美xxxx黑人xx丫x性爽| 别揉我奶头 嗯啊视频| 伊人久久精品亚洲午夜| 2021天堂中文幕一二区在线观| a级毛片免费高清观看在线播放| 久久久色成人| 精品久久久久久久久久久久久| 国产黄片美女视频| 婷婷色综合大香蕉| 成人午夜高清在线视频| 秋霞在线观看毛片| 99精品在免费线老司机午夜| 高清日韩中文字幕在线| 国产精品久久久久久av不卡| 蜜桃亚洲精品一区二区三区| 国产精品99久久久久久久久| 欧美又色又爽又黄视频| 欧美日本亚洲视频在线播放| 99久久九九国产精品国产免费| 国产成人91sexporn| 高清毛片免费看| 美女脱内裤让男人舔精品视频 | 国产精品久久电影中文字幕| 国产伦精品一区二区三区四那| 国产乱人视频| 色综合站精品国产| 毛片女人毛片| 日韩视频在线欧美| 国产免费一级a男人的天堂| 亚洲精品色激情综合| 国产亚洲av嫩草精品影院| 日韩欧美国产在线观看| 九九在线视频观看精品| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲成人中文字幕在线播放| 亚洲av一区综合| 免费看av在线观看网站| 国产伦在线观看视频一区| 欧美日本亚洲视频在线播放| 观看美女的网站| 成人鲁丝片一二三区免费| 免费看光身美女| 五月伊人婷婷丁香| 亚洲aⅴ乱码一区二区在线播放| 亚洲欧美精品综合久久99| 九草在线视频观看| 精品久久久久久久久亚洲| 91久久精品国产一区二区成人| 蜜桃亚洲精品一区二区三区| www.av在线官网国产| 久久精品国产亚洲av香蕉五月| 免费搜索国产男女视频| 亚洲精品乱码久久久v下载方式| 亚洲成人av在线免费| 亚州av有码| 久久久色成人| 国产激情偷乱视频一区二区| 国产三级在线视频| 在线观看美女被高潮喷水网站| 婷婷色av中文字幕| 搞女人的毛片| 在线观看一区二区三区| 中国美女看黄片| 联通29元200g的流量卡| 91久久精品电影网| 亚洲欧洲日产国产| 少妇人妻精品综合一区二区 | 蜜臀久久99精品久久宅男| 一夜夜www| 精品人妻一区二区三区麻豆| 亚洲第一区二区三区不卡| 亚洲熟妇中文字幕五十中出| 亚洲欧美日韩高清专用| 久久这里有精品视频免费| 蜜臀久久99精品久久宅男| 国产精品女同一区二区软件| 人妻久久中文字幕网| 亚洲最大成人中文| 少妇被粗大猛烈的视频| 国产高清有码在线观看视频| 亚洲图色成人| 精品熟女少妇av免费看| 日韩欧美三级三区| 国产精品嫩草影院av在线观看| av专区在线播放| 亚洲精品乱码久久久久久按摩| 午夜亚洲福利在线播放| 中文字幕制服av| 成人亚洲欧美一区二区av| 国国产精品蜜臀av免费| 麻豆av噜噜一区二区三区| 在线a可以看的网站| 国产精品国产三级国产av玫瑰| 国产淫片久久久久久久久| 91精品一卡2卡3卡4卡| 亚洲乱码一区二区免费版| 国产高清有码在线观看视频| 国产在视频线在精品| 中文亚洲av片在线观看爽| 国产三级在线视频| 久久综合国产亚洲精品| 日韩 亚洲 欧美在线| 国产探花在线观看一区二区| 亚洲精品亚洲一区二区| 日韩av在线大香蕉| 日韩欧美国产在线观看| 一区二区三区高清视频在线| 超碰av人人做人人爽久久| 亚洲久久久久久中文字幕| 国产午夜福利久久久久久| 国产视频首页在线观看| 国产精品久久电影中文字幕| 精品一区二区三区人妻视频| 深爱激情五月婷婷| 成熟少妇高潮喷水视频| 亚洲电影在线观看av| 少妇的逼好多水| 特大巨黑吊av在线直播| 日本-黄色视频高清免费观看| 色哟哟·www| 免费黄网站久久成人精品| 亚洲成人久久性| 噜噜噜噜噜久久久久久91| 亚洲精品乱码久久久v下载方式| 久久热精品热| 亚洲av第一区精品v没综合| 老师上课跳d突然被开到最大视频| 亚洲精品亚洲一区二区| 亚洲aⅴ乱码一区二区在线播放| 欧美一级a爱片免费观看看| av在线观看视频网站免费| 嫩草影院入口| 看黄色毛片网站| 色吧在线观看| 99热6这里只有精品| 久久精品91蜜桃| 欧美日韩乱码在线| 亚洲av成人精品一区久久| 国产三级在线视频| 午夜免费激情av| 精品少妇黑人巨大在线播放 | 少妇高潮的动态图| 少妇的逼好多水| 免费黄网站久久成人精品| 国产精品久久久久久亚洲av鲁大| 夜夜爽天天搞| 中国美白少妇内射xxxbb| 长腿黑丝高跟| 国产精品永久免费网站| 2022亚洲国产成人精品| 国产成人影院久久av| 亚洲性久久影院| 国产精品一区二区三区四区久久| 男人狂女人下面高潮的视频| 国产免费一级a男人的天堂| 狠狠狠狠99中文字幕| 国产欧美日韩精品一区二区| 久久99热这里只有精品18| a级毛色黄片| 天堂中文最新版在线下载 | 丰满乱子伦码专区| 国产精品永久免费网站| 欧美+亚洲+日韩+国产| 男女视频在线观看网站免费| 欧美高清成人免费视频www| 丰满的人妻完整版| 亚洲精品456在线播放app| 永久网站在线| 亚洲欧美日韩卡通动漫| 看片在线看免费视频| 免费电影在线观看免费观看| 嘟嘟电影网在线观看| 日韩一本色道免费dvd| 一本久久精品| 伦精品一区二区三区| АⅤ资源中文在线天堂| 亚洲不卡免费看| 国产亚洲5aaaaa淫片| 久久久久网色| 女同久久另类99精品国产91| 黄色视频,在线免费观看| 国产色婷婷99| 日韩中字成人| 最近的中文字幕免费完整| 在线观看66精品国产| 搡老妇女老女人老熟妇| 免费观看精品视频网站| 久久婷婷人人爽人人干人人爱| 欧美色欧美亚洲另类二区| 国内精品美女久久久久久| 成人午夜高清在线视频| 日韩欧美在线乱码| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 99riav亚洲国产免费| 亚洲国产日韩欧美精品在线观看| 精品午夜福利在线看| 成年av动漫网址| 人人妻人人澡欧美一区二区| 亚洲av免费高清在线观看| 白带黄色成豆腐渣| 我要看日韩黄色一级片| 联通29元200g的流量卡| 色综合亚洲欧美另类图片| 六月丁香七月| 亚洲国产日韩欧美精品在线观看| 国产精品人妻久久久影院| 如何舔出高潮| 日产精品乱码卡一卡2卡三| 国产老妇伦熟女老妇高清| 免费av不卡在线播放| 午夜亚洲福利在线播放| 亚洲人成网站在线观看播放| 久久久久久久久久黄片| 国产久久久一区二区三区| 亚洲欧美日韩无卡精品| 国产成人aa在线观看| 大香蕉久久网| 国产精品三级大全| 亚洲国产精品国产精品| 成人三级黄色视频| 免费在线观看成人毛片| 人妻夜夜爽99麻豆av| 成人无遮挡网站| 成人高潮视频无遮挡免费网站| 九九热线精品视视频播放| 国产成人精品一,二区 | 欧美成人免费av一区二区三区| 波多野结衣高清无吗| av女优亚洲男人天堂| 插逼视频在线观看| 成人午夜高清在线视频| 亚洲精华国产精华液的使用体验 | 久久久精品大字幕| 日韩制服骚丝袜av| 欧美极品一区二区三区四区| 亚洲成人精品中文字幕电影| 美女 人体艺术 gogo| 国内精品美女久久久久久| 日本一二三区视频观看| 日本熟妇午夜| 18禁在线播放成人免费| 亚洲人成网站在线观看播放| 亚洲自偷自拍三级| av在线老鸭窝| 亚洲综合色惰| 国产精品久久电影中文字幕| 久久国产乱子免费精品| 国产精品一二三区在线看| 悠悠久久av| 男人和女人高潮做爰伦理| 亚洲欧美精品自产自拍| 男人的好看免费观看在线视频| 免费av不卡在线播放| 卡戴珊不雅视频在线播放| 久久久久国产网址| 免费大片18禁| 色5月婷婷丁香| 国产一区二区三区在线臀色熟女| 成人欧美大片| 国产av不卡久久| 国产色婷婷99| 亚洲精品色激情综合| 最近最新中文字幕大全电影3| 最新中文字幕久久久久| 成熟少妇高潮喷水视频| 久久国产乱子免费精品| 久久99热6这里只有精品| 久久精品人妻少妇| 熟女人妻精品中文字幕| 99久久久亚洲精品蜜臀av| 国产老妇伦熟女老妇高清| 日本-黄色视频高清免费观看| 最近手机中文字幕大全| av.在线天堂| 亚洲欧美清纯卡通| 亚洲电影在线观看av| 春色校园在线视频观看| 人体艺术视频欧美日本| 亚洲婷婷狠狠爱综合网| 亚洲av二区三区四区| 国产精品久久视频播放| 婷婷精品国产亚洲av| 国产成人a∨麻豆精品| 黄色一级大片看看| 一区二区三区四区激情视频 | 午夜亚洲福利在线播放| 国产伦理片在线播放av一区 | 中出人妻视频一区二区| 国产真实伦视频高清在线观看| 欧美丝袜亚洲另类| 成人毛片a级毛片在线播放| 麻豆国产av国片精品| 国产精品综合久久久久久久免费| 亚洲最大成人手机在线| 亚洲欧美精品自产自拍| 亚洲欧美日韩高清专用| 日韩欧美 国产精品| 日韩 亚洲 欧美在线| 中文字幕av成人在线电影| 国产真实乱freesex| ponron亚洲| 熟女人妻精品中文字幕| 日韩av不卡免费在线播放| 少妇人妻精品综合一区二区 | 国产精品爽爽va在线观看网站| 欧美高清成人免费视频www| or卡值多少钱| 国产午夜精品论理片| 国产黄片视频在线免费观看| 久久精品国产亚洲av涩爱 | 免费看光身美女| 免费搜索国产男女视频| h日本视频在线播放| 99久久无色码亚洲精品果冻| 伦理电影大哥的女人| 性插视频无遮挡在线免费观看| 黄色日韩在线| 一级毛片我不卡| 成人特级黄色片久久久久久久| 欧美成人一区二区免费高清观看| 熟女电影av网| 国产片特级美女逼逼视频| 大型黄色视频在线免费观看| 麻豆精品久久久久久蜜桃| 国产精品无大码| 91精品国产九色| 久久亚洲精品不卡| 中文在线观看免费www的网站| 国产成年人精品一区二区| 欧美最黄视频在线播放免费| 久久久久久久午夜电影| 国产精品不卡视频一区二区| 日本免费a在线| 大香蕉久久网| 日韩欧美精品免费久久| 1000部很黄的大片| 在线免费观看不下载黄p国产| 亚洲欧美精品专区久久| 如何舔出高潮| 久久亚洲精品不卡| 午夜久久久久精精品| 亚洲av成人av| 国产午夜精品一二区理论片| 黄色视频,在线免费观看| 九草在线视频观看| 美女内射精品一级片tv| 日本-黄色视频高清免费观看| 熟女人妻精品中文字幕| 日产精品乱码卡一卡2卡三| 嘟嘟电影网在线观看| 亚洲美女视频黄频| 国产色爽女视频免费观看| АⅤ资源中文在线天堂| 亚洲成av人片在线播放无| 天美传媒精品一区二区| 婷婷精品国产亚洲av| 午夜精品在线福利| 国产三级中文精品| 欧美成人精品欧美一级黄| 日日撸夜夜添| 大又大粗又爽又黄少妇毛片口| 一级毛片电影观看 | 成熟少妇高潮喷水视频| 精品一区二区免费观看| 亚洲在线自拍视频| 综合色av麻豆| 婷婷亚洲欧美| 久久国内精品自在自线图片| 男女啪啪激烈高潮av片| 免费看日本二区| 91麻豆精品激情在线观看国产| 啦啦啦观看免费观看视频高清| 国产一区二区三区在线臀色熟女| 黄色一级大片看看| 午夜福利视频1000在线观看| 久久精品人妻少妇|