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

    Application of Br?nsted acid ionic liquids as green catalyst in the synthesis of 2-propanol with reactive distillation☆

    2016-06-07 05:44:20JinbeiYangDongrenCaiTingZengLihuaZhouLingLiRuoyuHongTingQiu
    Chinese Journal of Chemical Engineering 2016年11期

    Jinbei Yang ,Dongren Cai,Ting Zeng ,Lihua Zhou ,Ling Li,Ruoyu Hong ,Ting Qiu ,*

    1 School of Chemical Engineering,Fuzhou University,Fuzhou,Fujian 350108,China

    2 School of Ocean Science and Biochemistry Engineering,Fuqing Branch of Fujian Normal University,Fuzhou,Fujian 350300,China

    1.Introduction

    2-Propanol(IPOH),which is an excellent organic raw material and solvent,is widely used as an intermediate in chemical synthesis,or as an important solvent in the ink,paints,dyes,cosmetics,pharmaceutical industry and food industry[1–6].Ink and paints are the mainly applications which account for about 50%of the total consumption of 2-propanol.Several methods are available for manufacture of 2-propanol.Among them,direct hydration and indirect hydration of propylene are used most widely[7].Propylene and water are used as raw materials of both the processes.Indirect hydration is based on a two step process,in which an ester is formed and then hydrolyzed to 2-propanol.Nevertheless,this route has been gradually phased out in China since the 1980s due to its complex process,low selectivity,high energy requirement and serious corrosion to equipment.Compared to indirect hydration,direct hydration of propylene avoids some corrosion and environment problems and it has been deemed as the preferred process for production of IPOH.However,the process remainssome disadvantages such as high consumption of energy and cost for separation of IPOH from the azeotropic mixture(IPOH and water)[8].Hydrogenation of acetone is another method for producing IPOH.Nevertheless,it is used rarely because of its high demand on raw materials and catalyst.Therefore,it is urgent to develop an efficient,stable,environmentally friendly catalyst and green process for production of IPOH.

    The preparation of IPOH via transesterification of is opropyl acetate(IPAc)with methanol(MeOH)as raw materials is a novel process as shown in Eq.(1),which has been investigated by our team using sodium methoxide as catalyst for the first time[3].The high conversion of IPAc(above 99%)was gotten under the optimal conditions.Hunan Zhongchuang Chemical Co.,Ltd.produce IPAc by the addition reaction of acetic acid with propylene on the 20000 t·a?1scale since 2006.This route has high conversion and selectivity,and greatly reduces manufacturing cost.Thus,IPAc is suitable as a feedstock for IPOH production via transesterification.However,sodium methoxide is not the best choice because it is easily inactivated by air or water and cannot be reused.

    Ionic liquids(ILs),kinds of environmentally friendly solvent and catalyst,have attracted significant attention of scholars from various fields owing to its adjustable physical and chemical properties[9–16].The main driving force to explore ILs is the fact that they offer distinctive set of attributes,such as non flammable thermally,negligible vapor pressure,excellent thermal stability,high catalytic activity,and excellent recyclability and designability[17–20].Reactions with acid-functionalized ionic liquids as catalysts have gained desired results,such as nitration[21],transesteri fication[22,23],hydrogenation[24],polymerization[25],and Beckmann rearrangement[26].Nevertheless,to the best of our knowledge,their application as catalyst for the synthesis of IPOH via transesterification of IPAc with MeOH hasrarely reported in the literature.

    In this work, five Br?nsted acidic ionic liquids were prepared and characterized.Their uses for the synthesis of IPOH via the transesterification of IPAc with MeOH were investigated.Then,the reaction kinetic behavior was studied using[Ps-mim]HSO4as catalyst,which performed best.Finally,a batch reactive distillation was proposed,and the reusability of catalyst was also studied.

    2.Experimental Section

    2.1.Materials

    Isopropyl acetate(>99.5 wt%analytical grade)was supplied by Hunan Zhongchuang Chemical Co.,Ltd.Other chemicals(AR grade)were commercial products and used without further purification.

    2.2.Preparation of Br?nsted acidic ILs

    In this paper, five Br?nsted acidic ILs were synthesized according to the procedures previously reported in the literature[27–30],which included 1-propylsulfonate-3-methyliminazole hydrogen sulfate([Ps-mim]HSO4),1-propylsulfonate pyridine hydrogensulfate([Ps-Py]HSO4),propylsulfonate triethylamine hydrogensulfate([Ps-N(Et)3]HSO4),N,N-dimethyl benzylamine-N-propylsulfonate hydrogensulfate([Ps-N-CH2C6H5Me2]HSO4)and N,N-dimethylcyclohexylamine-N-propylsulfonate hydrogensulfate ([Ps-NC6H11Me2]HSO).The structures of ILs were illustrated in Fig.1.

    2.3.Characterization of Br?nsted acidic ILs

    NMR spectra were obtained on a Bruker AV500 spectrometer(Bruker Co.,Ltd.,Switzerland)in DMSO-d6 and calibrated with TMS as the internal reference.IR measurements of ILs were recorded using a Spectrum 2000 FT-IR absorption spectrometer(Perkin Elmer Inc.,American)for KBr pellets in the frequency range of 4000–400 cm?1.

    The acidity of ILs was determined by using the theoretical model and UV–visible spectroscopic methods.UV–visible spectra were performed on an Elmer Lambda 900 spectrometer(Perkin Elmer Inc.,American)with a basic indicator(4-nitroanline)by following the literature reported previously[31,32],and the solvent was methanol.

    2.3.1.NMR analysis

    The spectraldata of1HNMRand13C NMRfor the five Br?nsted acidic ILs were shown as follows:

    2.3.2.FT-IR analysis

    The spectral data of FT-IR for the five Br?nsted acidic ILs were listed as follows:

    Fig.1.Structures of five Br?nsted acidic ILs.

    The structures of the five Br?nsted acidic ILs were investigated using1HNMR,13CNMR and FT-IR spectroscopy,which were coincided with the theoretical structure(Fig.1).Moreover,1HNMR spectra showed no impurities,which demonstrate that high-purity five ILs were synthesized.

    2.3.3.Hammett acidity evaluation

    The acidity strength of Br?nsted acidic ILs can be efficiently expressed by the Hammett acidity function,H0,which can be calculated by the following equation:

    According to the Lambert–Beer's Law,the[I]/[IH+](I represents indicator)could be determined and calculated from the changes of UV–visible absorbance of 4-nitroanline after the addition of Br?nsted acidic ILs.In this paper,the maximal absorbance of the unprotonated form of the indicator was observed at 370 nm under the same concentration of 4-nitroanline(54 μmol·L?1,p Ka=0.99)and ILs(80 mmol·L?1)in methanol.The absorption curves and the calculated of Hammett acidity values were shown in Fig.2 and Table 1,respectively.

    Fig.2.UV visible spectra of the five Br?nsted acidic ILs.(a)blank;(b)[Ps-NCH2C6H5Me2]HSO4;(c)[Ps-N(Et)3]HSO4;(d)[Ps-N-C6H11Me2]HSO4;(e)[Ps-Py]HSO4;(f)[Ps-mim]HSO4.

    Table 1 Calculation of Hammettacidity function(H0)values of different ILs in methanol(293.15 K)

    As shown in Table 1 and Fig.2,the results showed that the order of the H0values of the five ILs is as follows:[Ps-N-CH2C6H5Me2]HSO4>[Ps-N(Et)3]HSO4>[Ps-N-C6H11Me2]HSO4>[Ps-Py]HSO4>[Ps-mim]HSO4.On the contrary,the acidity of[Ps-mim]HSO4is strongest.It is clearly shown that the acidity of the ILs depended on the types of cations,although the influence was insignificant.

    2.4.Experimental apparatus and procedure for 2-propanol synthesis

    In this study,the boiling point of the reaction system was relatively low at atmospheric pressure(about 338.15 K).Therefore,the experiment of the transesterification was carried out at high pressure(higher than 0.2 MPa)in order to increase the temperature to enhance the reaction rate.

    2.4.1.Kinetic experiments

    The kinetic experiments of the transesterification were conducted in a stainless steel reactor(reaction volume:500 ml,material:316 L)equipped with an agitation and temperature-controlling device(±0.1 K)(Fig.3).A HPLC pump was used for injection of the raw material into the reactor.

    The desired amount mixture ofIPAc and MeOHwas charged into the reactor firstly.After the reactor was sealed up,nitrogen was introduced into the reactor to ensure the reactants keep liquid under the desired reaction temperature.Then,the reactants were heated to specified temperature with moderate agitation.Once the reaction temperature was reached,the preheated catalyst(dissolved in methanol)was fed into the reactor quickly by the HPLC pump,at the same time,the stirrer was set at an appointed rate and the time was regarded as the initial time.Samples were withdrawn at a fixed time interval,cooled down quickly and then analyzed by gas chromatography(GC2014,Shimadzu Corporation,Japan).

    The experiments at the reaction temperature range of 338.15–393.15 K,the initial molar ratio of MeOH to IPAc of 2.0–5.0 and the catalyst dosage of 1.0 wt%,1.5 wt%and 2.0 wt%were studied.The reaction could be considered as reaching chemical equilibrium when the composition of the reaction mixture was nearly constant.All samples were performed in triplicate and mean values were quoted as results.

    2.4.2.Batch reactive distillation experiments

    The batch reactive distillation experiment of the transesterification was performed in a stainless steel column(internal diameter:20 mm,height:3 m,bottom volume:500 ml)equipped with an electronic reflux splitter to control reflux ratio and a heater in the bottom to heat the reaction mixture.The column was packed with θ stainless steel packing(φ3×3).There was a condenser located at the top of column for cooling down the vapor.The equipment diagram was shown in Fig.4.

    Firstly,the reactants including catalyst ILs were proportionately added to the reboiler.The column pressure was increased to the specified value by introducing nitrogen after all correlative equipment was sealed up.Secondly,the reactants were heated and then total reflux operation was lasted for 10 min after the distillate appeared at the top.Subsequently,the distillates were continuously removed out.Samples were withdrawn at different time intervals from the reboiler,cooled down quickly and then analyzed in triplicate by gas chromatography until the pointed reaction time was reached.The used catalyst ILs could be easily separated from mixture of the bottom by removing the sample mixture with the rotary evaporator(353 K,2000 Pa),and reused for the next time.

    2.4.3.Sample analysis method

    All samples were quantitatively analyzed by gas chromatography(GC 2014,Shimadzu Corporation,Japan)equipped with a FID detector and AT-FFAP capillary column(50 m×0.32 mm×0.5μm),using 1,4-dioxane as internal standard and N,N-dimethyl for mamide as solvent.The temperature program of column was as follows:started at 328.15 K for 1 min,increased to 338.15 K at 2.5 K·min?1and maintained for 1 min,and then increased to 353.15 K at 30 K·min?1,the temperature was increased to 453.15 Kat40 K·min?1finally.Injector and detector temperatures were both maintained at 533.15 K.

    3.Results and Discussion

    3.1.Performance of five IL catalysts

    Fig.3.Equipment diagram of kinetic experimental.1—heating mantle,2—high pressure reactor,3—thermocouple,4—mechanical stirrer,5—pressure tap,6—HPLC pump,7—control cabinet,8—exhaust valve,9—condenser,10—nitrogen cylinder.

    Fig.4.Equipment diagram of high-pressure reactive distillation.1—reboiler,2—stainless steel packing,3—column,4—condenser,5—reflux splitter,6—re flux splitter controller,7—stainless steel receiving tanks,8—exhaust valve,9—nitrogen cylinder.

    The five IL catalysts were tested for their activity under the initial MeOH/IPAc molar ratio of 3.0,catalyst dosage of 1.5%(mass fraction,based on the total reactant mass),reaction temperature of 353.15 K,and reaction time of 240 min.The results were depicted in Table 2.As seen in Table 2,all ILs are active for the transesterification,and the catalytic activity(conversion of the IPAc)of these ILs was listed in the following order:[Ps-mim]HSO4>[Ps-Py]HSO4>[Ps-NC6H11Me2]HSO4>[Ps-N(Et)3]HSO4>[Ps-N-CH2C6H5Me2]HSO4,which is consistent with the order of acidity of these ILs(Table 1).

    Based on the Hammett acidity evaluation and the performance testing of five acidic IL catalysts,it can be concluded that the catalyticactivity of the ILs is influenced by the type of cation,although the effect was insignificant.The effect of cation on the catalytic activity was as follows:[Ps-mim]+>[Ps-Py]+>[Ps-N-C6H11Me2]+>[Ps-N(Et)3]+>[Ps-N-CH2C6H5Me2]+.Obviously,[Ps-mim]HSO4shows highest catalytic activity and the yield of[Ps-mim]HSO4is highest among the five acidic ILs.

    Table 2 Catalytic activity of different ILs on the conversion of IPAc①

    As reported in earlier literature[28,33],the thermal decomposition temperature of[Ps-mim]HSO4is above 320°C,which shows that[Ps-mim]HSO4has high thermal stability,wide liquid ranges and suitable to be used as catalyst under high temperature compared to the resin,solid acid catalyst etc.Therefore,[Ps-mim]HSO4was chosen as catalyst for reaction kinetics and batch reactive distillation experiment in follow-up experiments.

    3.2.Reaction kinetics

    3.2.1.Effect of reaction temperature

    The effect of reaction temperature is important to the kinetic experiments,from which the activation energy of the transesterification of IPAc with MeOHcan be obtained.During the course of the experiments,the reaction temperature was varied in the range of 338.15–393.15 K,and the results were illustrated in Fig.5.As can be seen in Fig.5,the increment of temperature is apparently favorable to accelerate the reaction rate of transesterification,while the equilibrium conversion decreases slightly with the increasing of reaction temperature.Hence,it can be concluded the transesterification reaction is a marginal exothermal reaction.

    Fig.5.Effect of temperature on the conversion of IPAc.Conditions:catalyst,[Psmim]HSO4;initial mole ratio of MeOH to IPAc,3:1;catalyst dosage,1.5 wt%.The dots represent experiment results;the lines represent model results.

    3.2.2.Effect of initial molar ratio

    The effect of initial molar ratio of MeOH to IPAc was studied experimentally.The mole ratio varied from 2:1 to 5:1,and the conversions of IPAc were obtained and shown in Fig.6.As shown in Fig.6,the reaction rate and equilibrium conversion are sensitive to the initial reactant molar ratio.The more the MeOH was added,the higher the conversion of IPAc was obtained in the same reaction time.Among them,the equilibrium conversion is 90.7%with MeOH to IPAc mole rate of 4:1 in 180 min.Further increase of this mole ratio makes weak influence to the conversion and reaction rate.Theoretically,the equilibrium constant should be the same at the same temperature.Equilibriumconstants calculated by Eq.(4)were listed in Table 3,from which we can see that they are all similar to each other.This is in a good agreement with the theory.

    Fig.6.Effect of molar ratio of MeOH to IPAc on the conversion of IPAc.Conditions:catalyst,[Ps-mim]HSO4;reaction temperature,363.15 K;catalyst dosage,1.5 wt%.The dots represent experiment results;the lines represent model results.

    Table 3 Equilibrium constants K eq at different molar ratio

    3.2.3.Effect of catalyst dosage

    A homogeneous catalyst,[Ps-mim]HSO4(soluble in methanol),was used as catalyst for the transesterification reaction.The catalyst dosage,defined as the mass percentage of the catalyst to the total feed mass,was varied over a range of 1.0 wt%–2.0 wt%at the temperature of 363.15 K and the initial molar ratio of MeOH to IPAc of 3.0.The results were exhibited in Fig.7.As illustrated in the figure,the reaction rate of IPAc increases with the increase of catalyst dosage from 1.0 wt%to 1.5 wt%,yetthe equilibrium conversion is constant.Furthermore,the reaction rate is improved slightly using 2.0 wt%catalyst.Therefore,a suitable catalyst dosage should be determined.

    3.2.4.Chemical equilibrium

    The equilibrium constants of the transesterification reaction catalyzed by[Ps-mim]HSO4in this work were defined by Eq.(3).The equilibrium constant Keqcan be deduced from the equilibrium composition obtained from the experiment:

    Fig.7.Effect of catalyst dosage on the conversion of IPAc.Conditions:catalyst,[Psmim]HSO4;reaction temperature,363.15 K;initial mole ratio of MeOH to IPAc,3:1.The dots represent experiment results;the lines represent model results.

    where ΔH is the activation energy,J·mol?1.

    In a small temperature interval,ΔH can be considered to be constant and therefore,

    where C is a constant.

    A plot of ln Keqversus 1/T was shown in Fig.8.The temperature dependence of the equilibrium constant can be expressed by Eq.(7).

    According to Eq.(7),the activation energy is?5.26 kJ·mol?1,which means that the reaction is weakly exothermal.The mean squared error between experimental and calculated equilibrium constants at different temperatures is 2.36×10?4.

    3.2.5.Reaction mechanism

    Fig.8.Temperature dependence of the chemical equilibrium constants:ln K eq versus 1/T(■)and the best fit of our data(—).

    The reaction mechanism of the transesterification catalyzed by acid-functionalized ionic liquids is similar to that of proton acid.Therefore,the possible reaction mechanism can be deduced as shown in Fig.9.Steps 1 and 3 are relatively fast reactions and they are always assumed to be in chemical equilibrium.Step 2 is considered as a rate-controlling step for the slowest reaction rate.All in all,the reaction mechanism can be predigested into two steps,as shown in Fig.10.

    Fig.9.Reaction mechanism of the transesterification of IPAc and MeOH.

    Fig.10.Simplified reaction mechanism of the transesterification of IPAc and MeOH.

    3.2.6.Kinetic modeling

    Based on the mechanism of the reaction as mentioned above(Fig.10),a kinetic model can be deduced:

    whereare the forward and backward reaction rate constants of Step 2,respectively.

    As nucleophilic substitution reaction are relatively fast,and considered to be in equilibrium,the rates of Step 1 can be expressed as:

    whereare the forward and backward reaction rate constants of step 1,respectively.Eqs.(8)and(9)lead to

    whereis the initial concentration ofIPAc.cH+is the concentration of catalyst.

    3.2.7.Estimation of the reaction rate constants

    The relationship between the rate constant k and temperature T can be expressed by the Arrhenius equation:

    where Eais the activation energy,J·mol?1,and A is the pre-exponential factor or apparent frequency factor,L·min?1mol?1.

    According to the nonlinear Eq.(11),a fourth-order Runge–Kutta method is used to calculate the conversion of IPAc under different time.An objective function(OF)is to minimize the squared differences between the calculated value Xcaland the experimental value Xexp,as shown in Eq.(13).The optimal parameters for the kinetics are estimated by using the Nelder–Mead simplex method,as shown in Table 4.

    The fitted parameters were selected to calculate the conversion of IPAc by the homogeneous reaction kinetic model,and the results were shown in Figs.5,6 and 7.There is close agreement between the calculated and experimental values.Therefore,the kinetic model can describe the kinetic behavior of the system reliably.

    3.3.Batch reactive distillation experiments

    According to the previous work,the transesterification of IPAc with MeOH is a reversible reaction.High conversion of IPAc can be obtained by removing one product out of the system in time.The transesterification was carried out in a stainless steel column as shown in Fig.4.During the experimental period,several operation conditions were studied:reaction pressure,catalyst dosage,initial molar ratio and reaction time.Moreover,the stability of the catalyst was also investigated.

    3.3.1.Effect of reaction pressure

    The reaction temperature and evaporation are both associated with operating pressure.The effect of reaction pressure on the conversion of IPAc was examined by fixing the catalyst dosage of 1.0 wt%,reaction molar ratio of MeOH to IPAc of 3:1 and reaction time of 3.0 h.As can be seen in Fig.12(a),the conversion of IPAc significantly increasesfrom 94.5%to the maximum value of 98.1%with the increase of the reaction pressure from 0.1 to 0.2 MPa.However,with further increase of the reaction pressure,the conversion of IPAc decreases.A possible reason is that the increase of reaction pressure leads to high reaction temperature,which accelerates the reaction rate and then the high conversion is obtained.Nevertheless,the amount of vaporization is decreasing with the increase of reaction pressure,which is not conducive to removing products out.Then,the forward reaction is inhibited.Therefore,the optimum reaction pressure is between 0.2 and 0.3 MPa,and the conversion of IPAc is above 97.5%.

    Table 4 Chemical reaction rate constant of the transesterification at various temperatures

    Fig.11.Arrhenius plotfor the reaction rate constants of IPAc with MeOH.Forward reaction k+(■)and backward reaction k? (▲).The lines represent the results of the linear regression.

    3.3.2.Effect of catalyst dosage

    The effect of catalyst dosage on the conversion of IPAc was studied and illustrated in Fig.12(b).The reaction conditions are:catalyst,[Psmim]HSO4;reaction pressure,0.3 MPa;initial mole ratio of MeOH to IPAc,3:1;and reaction time,3.0 h.As the catalyst dosage increased from 0.25 to 1.0 wt%,the conversion of IPAc was found to increase from 92.4%to 97.6%due to the increase of catalytic active sites.Nevertheless,with further increase of the catalyst dosage,the conversion of IPAc changes insignificantly.As a result,in consequence of the cost and activity of catalyst,the suitable catalyst dosage is around 1.0 wt%.

    Fig.12.Conversion of IPAc by transesterification over[Ps-mim]HSO4 as a function of(a)reaction pressure,(b)catalyst dosage,(c)reaction molar ratio,and(d)reaction time.

    3.3.3.Effect of reaction molar ratio

    Fig.12(c)shows the effect of molar ratio of reactants(MeOH/IPAc)on the conversion of IPAc.The reaction conditions are:catalyst,[Psmim]HSO4;catalyst dosage,1.0 wt%;reaction pressure,0.3 MPa;and reaction time,3.0 h.Overall,the reaction molar ratio has great effect on the conversion of IPAc as depicted in Fig.12(c).As the reaction molar ratio increased from 2.0 to 4.0,the conversion of IPAc raises obviously from 92.0%to 99.4%.However,when the reaction molar ratio is larger than 4,the effect on conversion is weak.Actually,the increase of initial molar ratio can lead to the increase of energy consumption of subsequent separation.Therefore,in the consideration of the conversion of IPAc and energy consumption of subsequent separation,a molar ratio(MeOH/IPAc)of 4.0 is the most suitable.

    3.3.4.Effect of reaction time

    The suitable reaction time was explored in a range of 1.0–3.5 h,as shown in Fig.12(d).The reaction conditions are:catalyst,[Psmim]HSO4;catalyst dosage,1.0 wt%;reaction pressure,0.3 MPa;and initial mole ratio of MeOH to IPAc,4:1.It can be seen from Fig.12(d)that the conversion of IPAc increases with the increase of reaction time.When the reaction time was longer than 3.0 h,the conversion remains invariant.This phenomenon can be explained by the fact that the transesterification of IPAc with MeOH is nearly completed with a time of around 3.0 h.Meanwhile,the temperature and mass fraction of MeAc and IPAc at the bottom were recorded and analyzed,respectively.The results were listed in Fig.13(a)and(b).It was observed that the temperature at the bottom of column is gradually increased with reaction time due to the continuous change of composition at the bottom of column(Fig.13(a)).As shown in Fig.13(b),the top temperature was close to the binary azeo tropic point(MeOH/MeAc).Hence,the preferred reactant time is 3.0 h.

    3.3.5.Repeatability experiment

    From the experimental study done in this work,the optimum batch reactive distillation operation conditions can be concluded as follows:reaction pressure,0.3 MPa;catalyst dosage,1.0 wt%;initial mole ratio of MeOH to IPAc,4:1;and reaction time,3.0 h.Under the optimum conditions,the conversion of IPAc is above 99%.

    For the purpose of verifying the reliability of the experimental data,a repeatability experiment was performed.The results were listed in Table 5.As depicted in Table 5,there is no significant difference in the three parallel experiments,which indicates that the experimental data are reliable and accurate.

    3.3.6.Reusability of catalyst

    In order to evaluate the potential reusability of catalyst in the transesterification of IPAc with MeOH,a series of recycle experimentswere conducted in stainless steel column under the optimum conditions.The reuse performances of the ILs[Ps-mim]HSO4were shown in Table 6.It can be seen thatthe conversion ofIPAc only slightly decreases from 99.4%to 99.0%after the catalyst was repeatedly used for four times.Furthermore,the FT-IR,1H NMR and13C NMR spectroscopy analyses of recycled ILs were carried out to verify that there are no new functional groups formed for the used IL catalyst.The results of recycle experiments and characterization both indicate that the ILs[Psmim]HSO4as the catalyst for the transesterification of IPAc with MeOH is recyclable and stable.

    Table 5 Repeatability experiment results

    Table 6 Effect of recycle times on the conversion of IPAc①

    3.4.Comparison of[Ps-mim]HSO4 and traditional catalyst

    According to the previous study of our team[3],sodium alkoxide was used as catalyst in the transesterification of IPAc with MeOH.The reaction reached equilibrium in 120 min at 333 K and the conversion of IPAc is about 88.0%with catalyst dosage of 0.3 wt%and initial mole ratio of MeOH to IPAc of 3:1.However,sodium methoxide is sensitive to water and air,and it cannot be reusable,which would produce large amounts of solid waste and cause environmental pollution.In this paper,[Ps-mim]HSO4was used as catalyst.The reaction reached equilibrium in 480 min at 353 K and the conversion of IPAc is 87.4%with catalyst dosage of 1.5 wt%and initial mole ratio of MeOH to IPAc of 3:1.Compared to sodium alkoxide,[Ps-mim]HSO4as catalyst can acquire almost the same equilibrium conversion,the reaction rate is relatively low.However,ionic liquids offer a new possibility for developing environmentally-friendly acidic catalysts owing its advantages,such as stability in air and water,excellent thermal stability,high catalytic activity,and easy reusability.In this work,the high conversion of IPAc was obtained under optimum conditions using[Ps-mim]HSO4as catalyst,and the catalyst was recycled easily by removing the reaction mixture with the rotary evaporator and reused without any further procedure.Ionic liquids are promising to replace conventional catalysts and therefore have a great potential for industrial application.

    Fig.13.Effect of reaction time on(a)the purity of MeAc and IPAc in bottom and(b)the bottom and top temperature.

    4.Conclusions

    In this work,we have prepared and characterized five Br?nsted acidic ionic liquids.The experimental results show that[Ps-mim]HSO4is more active than the other ILs in the synthesis of 2-propanol(IPOH)via the transesterification of isopropyl acetate(IPAc)with methanol(MeOH),which was chosen as catalyst for further study.The reaction kinetics of IPOH synthesis was studied.A homogeneous second order kinetic model for the transesterification was established,and it has been found that the kinetic model established can be well used in describing the transesterification.The high-pressure batch reactive distillation experiment of the transesterification was performed in a stainless steel column.As a result,a high conversion of IPAc of 99.4%was obtained under the optimal reaction conditions.Moreover,the catalyst[Psmim]HSO4can be recycled easily by the rotary evaporator and reused without any further treatment.The catalyst had been repeatedly used for four times and no obvious changes in the structure of catalyst could be observed.

    [1]E.Grousseau,J.Lu,N.Gorret,S.E.Guillouet,A.J.Sinskey,Isopropanolproduction with engineered Cupriavidus necator as bioproduction platform,Appl.Microbiol.Biotechnol.98(9)(2014)4277–4290.

    [2]K.Inokuma,J.C.Liao,M.Okamoto,T.Hanai,Improvement of isopropanol production by metabolically engineered Escherichia coli using gas stripping,J.Biosci.Bioeng.110(6)(2010)696–701.

    [3]T.Qiu,P.Zhang,J.B.Yang,L.Xiao,C.S.Ye,Novel procedure for production of isopropanol by transesterification of isopropyl acetate with reactive distillation,Ind.Eng.Chem.Res.53(36)(2014)13881–13891.

    [4]L.Xiao,Q.L.Wang,H.X.Wang,T.Qiu,Isobaric vapor-liquid equilibrium data for the binary system methyl acetate+isopropylacetate and the quaternary system methyl acetate+methanol+isopropanol+isopropyl acetate at 101.3 kPa,Fluid Phase Equilib.344(2013)79–83.

    [5]S.J.Wang,D.S.Wong,Control of reactive distillation production of high-purity isopropanol,J.Process Control 16(4)(2006)385–394.

    [6]A.Rahman,S.Salem,A review on reduction of acetone to isopropanol with Ni nano superactive,heterogeneous catalysts as an environmentally benevolent approach,Appl.Catal.A Gen.469(2014)517–523.

    [7]Y.Xu,K.T.Chuang,A.R.Sanger,Design of a process for production of isopropyl alcohol by hydration of propylene in a catalytic distillation column,Chem.Eng.Res.Des.80(6)(2002)686–694.

    [8]S.M.Mahajani,M.M.Sharma,T.Sridhar,Direct hydration of propylene in liquid phase and under supercritical conditions in the presence of solid acid catalysts,Chem.Eng.Sci.57(22)(2002)4877–4882.

    [9]M.I.Kim,S.J.Choi,D.W.Kim,D.W.Park,Catalytic performance of zinc containing ionic liquids immobilized on silica for the synthesis of cyclic carbonates,J.Ind.Eng.Chem.20(2014)3102–3107.

    [10]L.Zhang,M.Xian,Y.He,L.Li,J.Yang,S.Yu,X.Xu,A Br?nsted acidic ionic liquid as an efficient and environmentally benign catalyst for biodiesel synthesis from free fatty acids and alcohols,Bioresour.Technol.100(19)(2009)4368–4373.

    [11]N.Muhammad,Y.A.Elsheikh,M.I.A.Mutalib,A.A.Bazmi,R.A.Khan,H.Khan,Z.J.Man,An overview of the role of ionic liquids in biodiesel reactions,J.Ind.Eng.Chem.21(2015)1–10.

    [12]T.Qiu,W.L.Tang,C.G.Li,C.M.Wu,L.Li,Reaction kinetics for synthesis of sec-butyl alcohol catalyzed by acid-functionalized ionic liquid,Chin.J.Chem.Eng.23(1)(2015)106–111.

    [13]H.X.Wang,C.M.Wu,X.W.Bu,W.L.Tang,L.Li,T.Qiu,A benign preparation of secbutanol via transesterification from sec-butyl acetate using the acidic imidazolium ionic liquids as catalysts,Chem.Eng.J.246(2014)366–372.

    [14]A.Ghorbani-Choghamarani,M.Norouzi,Synthesis and characterization of ionic liquid immobilized on magnetic nanoparticles:A recyclable heterogeneous organocatalyst for the acetylation of alcohols,J.Magn.Magn.Mater.401(2016)832–840.

    [15]M.Olkiewicz,N.V.Plechkova,M.J.Earle,A.Fabregat,F.Stüber,A.Fortuny,J.Font,C.Bengoa,Biodieselproduction from sewage sludge lipids catalysed by Br?nsted acidic ionic liquids,Appl.Catal.B Environ.181(2016)738–746.

    [16]A.Chinnappan,A.H.Tamboli,W.J.Chung,H.Kim,Green synthesis,characterization and catalytic efficiency of hypercross-linked porous polymeric ionic liquid networks towards 4-nitrophenol reduction,Chem.Eng.J.285(2016)554–561.

    [17]K.Scott,N.Basov,R.J.J.Jachuck,N.Winterton,A.Cooper,C.Davies,Reactor studies of supported ionic liquids:Rhodium-catalysed hydrogenation of propene,Chem.Eng.Res.Des.83(10)(2005)1179–1185.

    [18]J.P.Hallett,T.Welton,Room-temperature ionic liquids:Solvents for synthesis and catalysis.2,Chem.Rev.111(5)(2011)3508–3576.

    [19]J.B.Yang,L.H.Zhou,X.T.Guo,L.Li,P.Zhang,R.Y.Hong,T.Qiu,Study on the esterification for ethylene glycol diacetate using supported ionic liquids as catalyst:Catalysts preparation,characterization,and reaction kinetics,process,Chem.Eng.J.280(2015)147–157.

    [20]M.Deetlefs,K.R.Seddon,M.Shara,Predicting physical properties of ionic liquids,Phys.Chem.Chem.Phys.8(5)(2006)642–649.

    [21]G.Cheng,X.Duan,X.Qi,C.Lu,Nitration of aromatic compounds with NO2/air catalyzed by sulfonic acid-functionalized ionic liquids,Catal.Commun.10(2)(2008)201–204.

    [22]K.X.Li,L.Chen,Z.C.Yan,H.L.Wang,Application of pyridinium ionic liquid as a recyclable catalyst for acid-catalyzed transesterification of Jatropha oil,Catal.Lett.139(3–4)(2010)151–156.

    [23]T.Qiu,X.T.Guo,J.B.Yang,L.H.Zhou,L.Li,H.X.Wang,Y.Niu,The synthesis of biodiesel from coconut oil using novel Br?nsted acidic ionic liquid as green catalyst,Chem.Eng.J.296(2016)71–78.

    [24]I.Cerna,P.Kluson,M.Bendova,T.Floris,H.Pelantova,T.Pekarek,Intensification of the use of ionic liquids as efficient reaction co-solvents in asymmetric hydrogenations,Chem.Eng.Process.50(3)(2011)264–272.

    [25]W.Och?dzan-Siod?ak,K.Dziubek,D.Siod?ak,Biphasic ethylene polymerisation using 1-n-alkyl-3-methylimidazolium tetrachloroaluminate ionic liquid as a medium of the Cp2TiCl2titanocene catalyst,Eur.Polym.J.44(11)(2008)3608–3614.

    [26]X.Liu,L.Xiao,H.Wu,Z.Li,J.Chen,C.Xia,Novel acidic ionic liquids mediated zinc chloride:Highly effective catalysts for the Beckmann rearrangement,Catal.Commun.10(5)(2009)424–427.

    [27]A.C.Cole,J.L.Jensen,I.Ntai,K.L.T.Tran,K.J.Weaver,D.C.Forbes,J.H.Davis,Novel Br?nsted acidic ionic liquids and their use as dual solvent-catalysts,J.Am.Chem.Soc.124(21)(2002)5962–5963.

    [28]H.Li,S.Yu,F.Liu,C.Xie,L.Li,Synthesis of dioctyl phthalate using acid functionalized ionic liquid as catalyst,Catal.Commun.8(11)(2007)1759–1762.

    [29]G.Schmidt-Naake,A.Schmalfu?,I.Woecht,Free radical polymerization in ionic liquids—In fluence of the IL-concentration and temperature,Chem.Eng.Res.Des.86(7)(2007)765–774.

    [30]J.Fraga-Dubreuil,K.Bourahla,M.Rahmouni,J.P.Bazureau,J.Hamelin,Catalysed esterifications in room temperature ionic liquids with acidic counteranion as recyclable reaction media,Catal.Commun.3(5)(2002)185–190.

    [31]C.Thomazeau,H.Olivier-Bourbigou,L.Magna,S.Luts,Determination of an acidic scale in room temperature ionic liquids,J.Am.Chem.Soc.125(18)(2003)5264–5265.

    [32]Y.Wang,D.Jiang,L.Dai,Novel Br?nsted acidic ionic liquids based on benzimidazolium cation:Synthesis and catalyzed acetalization of aromatic aldehydes with diols,Catal.Commun.9(15)(2008)2475–2480.

    [33]Q.Wu,H.Chen,M.Han,D.Wang,J.Wang,transesterification of cottonseed oil catalyzed by Br?nsted acidic ionic liquids,Ind.Eng.Chem.Res.46(24)(2007)7955–7960.

    9色porny在线观看| 高清在线视频一区二区三区| 香蕉精品网在线| 久久狼人影院| 色94色欧美一区二区| 免费高清在线观看视频在线观看| 大码成人一级视频| 一区二区三区四区激情视频| 久久国内精品自在自线图片| 久久影院123| 国产高清国产精品国产三级| av国产久精品久网站免费入址| 国产日韩欧美视频二区| 乱码一卡2卡4卡精品| 亚洲在久久综合| 免费女性裸体啪啪无遮挡网站| 91aial.com中文字幕在线观看| 中文字幕人妻丝袜制服| 黄色配什么色好看| 岛国毛片在线播放| 日本免费在线观看一区| 免费观看a级毛片全部| 丰满乱子伦码专区| 2021少妇久久久久久久久久久| 飞空精品影院首页| 18+在线观看网站| 亚洲第一区二区三区不卡| 免费播放大片免费观看视频在线观看| 亚洲国产欧美在线一区| 观看美女的网站| 亚洲 欧美一区二区三区| 在线免费观看不下载黄p国产| 午夜福利视频在线观看免费| 伦理电影大哥的女人| 日本与韩国留学比较| 高清av免费在线| 欧美日韩国产mv在线观看视频| 国产 一区精品| 国产精品麻豆人妻色哟哟久久| 女的被弄到高潮叫床怎么办| 亚洲精品av麻豆狂野| 国产精品99久久99久久久不卡 | 国产男人的电影天堂91| 精品国产国语对白av| 999精品在线视频| 久久久欧美国产精品| 在线看a的网站| 少妇人妻久久综合中文| 日日爽夜夜爽网站| 丰满饥渴人妻一区二区三| 午夜福利乱码中文字幕| 精品熟女少妇av免费看| 精品国产露脸久久av麻豆| 国产精品国产三级专区第一集| 成人手机av| 一级a做视频免费观看| 欧美国产精品va在线观看不卡| 99久久人妻综合| 国产在线免费精品| 国产精品无大码| 狂野欧美激情性xxxx在线观看| 国产片特级美女逼逼视频| 国产成人欧美| 日韩精品免费视频一区二区三区 | 男女国产视频网站| 国产老妇伦熟女老妇高清| 97在线人人人人妻| 黄色一级大片看看| 五月伊人婷婷丁香| 久久综合国产亚洲精品| 国产精品国产三级国产专区5o| 久久这里有精品视频免费| 国产成人精品久久久久久| 日日撸夜夜添| 啦啦啦视频在线资源免费观看| 九色亚洲精品在线播放| 天天躁夜夜躁狠狠躁躁| 亚洲国产看品久久| 你懂的网址亚洲精品在线观看| 精品国产一区二区久久| 黄色怎么调成土黄色| 建设人人有责人人尽责人人享有的| 久久久精品免费免费高清| 你懂的网址亚洲精品在线观看| 99热网站在线观看| 热re99久久国产66热| 中文字幕免费在线视频6| 国产欧美亚洲国产| 日本色播在线视频| 久久人人爽av亚洲精品天堂| 少妇熟女欧美另类| 少妇的逼水好多| 国产高清不卡午夜福利| 国产极品粉嫩免费观看在线| 亚洲精品久久午夜乱码| 蜜桃国产av成人99| 欧美+日韩+精品| 亚洲精品,欧美精品| 成人影院久久| 婷婷色综合www| 晚上一个人看的免费电影| 国产无遮挡羞羞视频在线观看| 日日啪夜夜爽| 久久亚洲国产成人精品v| 国产欧美亚洲国产| 中文天堂在线官网| 制服诱惑二区| 大码成人一级视频| 国产片特级美女逼逼视频| 国产淫语在线视频| 边亲边吃奶的免费视频| 色哟哟·www| 一本大道久久a久久精品| 丝袜喷水一区| 亚洲av电影在线进入| 草草在线视频免费看| 91精品伊人久久大香线蕉| 中文字幕免费在线视频6| videosex国产| 26uuu在线亚洲综合色| 中文字幕av电影在线播放| 国产精品蜜桃在线观看| 亚洲国产日韩一区二区| 少妇被粗大猛烈的视频| 一二三四中文在线观看免费高清| 中文字幕免费在线视频6| 久久精品人人爽人人爽视色| 欧美亚洲 丝袜 人妻 在线| 极品人妻少妇av视频| 亚洲成国产人片在线观看| 一级毛片电影观看| 午夜91福利影院| 99热这里只有是精品在线观看| 亚洲精品中文字幕在线视频| 国产伦理片在线播放av一区| 中文欧美无线码| 免费大片18禁| 亚洲av综合色区一区| 99国产精品免费福利视频| 18在线观看网站| 成人18禁高潮啪啪吃奶动态图| 秋霞在线观看毛片| 亚洲精华国产精华液的使用体验| 欧美亚洲 丝袜 人妻 在线| 视频中文字幕在线观看| 欧美3d第一页| 成年人免费黄色播放视频| 成人手机av| 精品久久久久久电影网| 久久综合国产亚洲精品| 亚洲国产欧美日韩在线播放| 免费大片黄手机在线观看| 国产国语露脸激情在线看| 丝袜美足系列| 在线观看人妻少妇| 国产成人精品无人区| 精品久久国产蜜桃| 国产日韩欧美视频二区| 精品人妻在线不人妻| 国产欧美亚洲国产| 亚洲美女黄色视频免费看| 日韩成人av中文字幕在线观看| 18禁动态无遮挡网站| 国产白丝娇喘喷水9色精品| www.色视频.com| 两性夫妻黄色片 | 国产精品蜜桃在线观看| 人人妻人人澡人人爽人人夜夜| 少妇人妻久久综合中文| 老司机影院成人| 亚洲精品久久成人aⅴ小说| 免费人妻精品一区二区三区视频| av国产精品久久久久影院| 内地一区二区视频在线| 九九在线视频观看精品| 免费不卡的大黄色大毛片视频在线观看| 午夜福利影视在线免费观看| 国产亚洲av片在线观看秒播厂| 国产有黄有色有爽视频| 国产亚洲欧美精品永久| 亚洲人与动物交配视频| 亚洲在久久综合| 色婷婷av一区二区三区视频| 大香蕉97超碰在线| 另类亚洲欧美激情| 国产成人精品无人区| 国产精品秋霞免费鲁丝片| 亚洲情色 制服丝袜| 亚洲国产欧美日韩在线播放| 老司机影院毛片| 国产精品人妻久久久影院| 黄色视频在线播放观看不卡| 伦精品一区二区三区| 国产极品粉嫩免费观看在线| 亚洲精品国产av成人精品| 国产淫语在线视频| 两性夫妻黄色片 | 精品人妻在线不人妻| 欧美精品一区二区免费开放| 国产1区2区3区精品| 久久久久久久亚洲中文字幕| 一级片免费观看大全| 亚洲精品久久午夜乱码| 极品少妇高潮喷水抽搐| 国产女主播在线喷水免费视频网站| 丰满迷人的少妇在线观看| 纯流量卡能插随身wifi吗| 欧美日韩一区二区视频在线观看视频在线| 九色成人免费人妻av| 久久人人爽人人片av| 国产在线视频一区二区| 少妇人妻久久综合中文| 精品少妇久久久久久888优播| 秋霞在线观看毛片| 在线免费观看不下载黄p国产| 久久亚洲国产成人精品v| 国产欧美日韩综合在线一区二区| 亚洲av电影在线观看一区二区三区| 欧美3d第一页| 99国产综合亚洲精品| 欧美少妇被猛烈插入视频| 国产精品久久久久久久久免| 最近中文字幕2019免费版| 大片免费播放器 马上看| 中文字幕亚洲精品专区| 巨乳人妻的诱惑在线观看| 日韩熟女老妇一区二区性免费视频| tube8黄色片| 91久久精品国产一区二区三区| 亚洲精品中文字幕在线视频| 91精品三级在线观看| 国产极品天堂在线| 黄色怎么调成土黄色| 午夜福利视频在线观看免费| 汤姆久久久久久久影院中文字幕| 一二三四中文在线观看免费高清| 国产色婷婷99| 两性夫妻黄色片 | av线在线观看网站| 国产永久视频网站| 国精品久久久久久国模美| 亚洲av电影在线观看一区二区三区| 日本av免费视频播放| 亚洲婷婷狠狠爱综合网| 80岁老熟妇乱子伦牲交| av福利片在线| 国产片特级美女逼逼视频| 不卡视频在线观看欧美| 色吧在线观看| 久久99热6这里只有精品| 国产成人精品在线电影| 免费久久久久久久精品成人欧美视频 | 午夜久久久在线观看| 国产一级毛片在线| 狠狠婷婷综合久久久久久88av| 亚洲欧洲精品一区二区精品久久久 | 国产有黄有色有爽视频| 制服人妻中文乱码| 丝瓜视频免费看黄片| 人成视频在线观看免费观看| 啦啦啦中文免费视频观看日本| 欧美精品高潮呻吟av久久| 欧美最新免费一区二区三区| 亚洲国产精品专区欧美| 人妻人人澡人人爽人人| 日本91视频免费播放| 国产免费现黄频在线看| 一级毛片 在线播放| 高清欧美精品videossex| 在线亚洲精品国产二区图片欧美| 久久鲁丝午夜福利片| 国产极品粉嫩免费观看在线| 在线 av 中文字幕| 99国产精品免费福利视频| 日韩三级伦理在线观看| 亚洲国产欧美日韩在线播放| 国产女主播在线喷水免费视频网站| 日日爽夜夜爽网站| 人妻 亚洲 视频| www.色视频.com| 国产亚洲一区二区精品| 日本-黄色视频高清免费观看| 综合色丁香网| 亚洲国产看品久久| 国产极品粉嫩免费观看在线| 狂野欧美激情性bbbbbb| 九草在线视频观看| 一级片'在线观看视频| 国产精品蜜桃在线观看| 久久久久久久久久久久大奶| 高清不卡的av网站| 午夜免费观看性视频| 成人国产麻豆网| 久久精品夜色国产| 岛国毛片在线播放| 亚洲欧洲精品一区二区精品久久久 | 一级片免费观看大全| 美女内射精品一级片tv| 午夜视频国产福利| 少妇人妻久久综合中文| 建设人人有责人人尽责人人享有的| 在线天堂最新版资源| 欧美性感艳星| 日本欧美视频一区| 欧美bdsm另类| 熟女人妻精品中文字幕| 90打野战视频偷拍视频| 国产精品一区二区在线不卡| 久久精品国产亚洲av天美| 欧美精品一区二区大全| 午夜av观看不卡| 国产高清国产精品国产三级| 免费看av在线观看网站| 精品人妻熟女毛片av久久网站| 最新中文字幕久久久久| 亚洲精品成人av观看孕妇| 亚洲成人av在线免费| 人人澡人人妻人| 午夜影院在线不卡| 99热国产这里只有精品6| 亚洲精品第二区| av.在线天堂| 亚洲欧美日韩卡通动漫| 日韩欧美一区视频在线观看| 成人黄色视频免费在线看| 久久久久网色| 五月开心婷婷网| www.熟女人妻精品国产 | 一本色道久久久久久精品综合| 另类精品久久| 交换朋友夫妻互换小说| 亚洲精品av麻豆狂野| 成人午夜精彩视频在线观看| 两性夫妻黄色片 | 黑人巨大精品欧美一区二区蜜桃 | 另类精品久久| 交换朋友夫妻互换小说| av国产久精品久网站免费入址| 国产成人一区二区在线| 一区二区av电影网| 国产成人aa在线观看| 插逼视频在线观看| 久久97久久精品| 少妇的逼水好多| 七月丁香在线播放| 久久久久久久大尺度免费视频| 亚洲一码二码三码区别大吗| 亚洲精品av麻豆狂野| 国产精品人妻久久久久久| 中国三级夫妇交换| 纯流量卡能插随身wifi吗| 汤姆久久久久久久影院中文字幕| 日韩免费高清中文字幕av| 日韩成人伦理影院| 日日撸夜夜添| 超碰97精品在线观看| 中文字幕另类日韩欧美亚洲嫩草| 欧美精品av麻豆av| 新久久久久国产一级毛片| 国产一区二区三区av在线| 国产黄色视频一区二区在线观看| 色94色欧美一区二区| 色视频在线一区二区三区| h视频一区二区三区| 国产精品国产三级国产av玫瑰| 婷婷色av中文字幕| 国产永久视频网站| 99久久精品国产国产毛片| 色94色欧美一区二区| 久久影院123| 男人添女人高潮全过程视频| 草草在线视频免费看| 美女视频免费永久观看网站| 天天躁夜夜躁狠狠久久av| 午夜91福利影院| av有码第一页| 中国国产av一级| 人妻 亚洲 视频| 中文精品一卡2卡3卡4更新| 精品人妻一区二区三区麻豆| 日韩伦理黄色片| 黄网站色视频无遮挡免费观看| 国产福利在线免费观看视频| 性色avwww在线观看| 看十八女毛片水多多多| 色视频在线一区二区三区| 国产精品人妻久久久久久| 一边亲一边摸免费视频| av在线app专区| 亚洲av成人精品一二三区| 久久精品久久久久久久性| 午夜日本视频在线| 中文乱码字字幕精品一区二区三区| 欧美最新免费一区二区三区| 桃花免费在线播放| 少妇熟女欧美另类| 黑丝袜美女国产一区| 精品一区二区免费观看| 三级国产精品片| 亚洲第一av免费看| 观看av在线不卡| 建设人人有责人人尽责人人享有的| 国产免费福利视频在线观看| 一区二区三区乱码不卡18| 在线天堂中文资源库| 青青草视频在线视频观看| 午夜免费男女啪啪视频观看| √禁漫天堂资源中文www| 一级爰片在线观看| 久久久久久伊人网av| av线在线观看网站| 亚洲欧美成人综合另类久久久| 欧美日韩亚洲高清精品| 成人黄色视频免费在线看| 性色av一级| 2022亚洲国产成人精品| 久久久久网色| 99国产综合亚洲精品| 有码 亚洲区| 久久精品国产鲁丝片午夜精品| 免费高清在线观看视频在线观看| 精品国产国语对白av| 免费看不卡的av| 一级a做视频免费观看| 男人添女人高潮全过程视频| av又黄又爽大尺度在线免费看| 午夜福利视频在线观看免费| 我的女老师完整版在线观看| 新久久久久国产一级毛片| 亚洲伊人久久精品综合| 内地一区二区视频在线| 免费播放大片免费观看视频在线观看| 18禁动态无遮挡网站| 婷婷色综合大香蕉| 国产精品国产三级国产专区5o| 亚洲精品日韩在线中文字幕| 日韩 亚洲 欧美在线| 高清黄色对白视频在线免费看| 欧美 亚洲 国产 日韩一| 岛国毛片在线播放| 免费观看性生交大片5| 亚洲精华国产精华液的使用体验| 日韩中文字幕视频在线看片| 国产深夜福利视频在线观看| 如日韩欧美国产精品一区二区三区| 国产 精品1| 亚洲第一区二区三区不卡| 亚洲精品国产av成人精品| 香蕉丝袜av| 婷婷色综合www| 免费女性裸体啪啪无遮挡网站| 又黄又粗又硬又大视频| 免费高清在线观看日韩| 在现免费观看毛片| 啦啦啦视频在线资源免费观看| 亚洲性久久影院| 男女高潮啪啪啪动态图| 国产免费又黄又爽又色| 爱豆传媒免费全集在线观看| 男女午夜视频在线观看 | 热re99久久国产66热| 日韩视频在线欧美| 搡女人真爽免费视频火全软件| 久久久久久久大尺度免费视频| 欧美bdsm另类| 国产在线一区二区三区精| av片东京热男人的天堂| 国产日韩欧美亚洲二区| 一级毛片 在线播放| 国产精品久久久av美女十八| 国产国拍精品亚洲av在线观看| 久久精品国产鲁丝片午夜精品| 亚洲美女视频黄频| 内地一区二区视频在线| 亚洲国产精品专区欧美| 久久久国产欧美日韩av| 亚洲av福利一区| 精品一品国产午夜福利视频| 久久精品久久久久久久性| 男女免费视频国产| 国产精品一区www在线观看| 亚洲av国产av综合av卡| 国产亚洲一区二区精品| av免费观看日本| 高清毛片免费看| 日韩大片免费观看网站| 男人添女人高潮全过程视频| 欧美国产精品一级二级三级| kizo精华| 看免费av毛片| 90打野战视频偷拍视频| 婷婷色综合www| 成人毛片a级毛片在线播放| 人妻少妇偷人精品九色| 蜜臀久久99精品久久宅男| 妹子高潮喷水视频| 日韩中字成人| 精品亚洲乱码少妇综合久久| 日韩中文字幕视频在线看片| 久久这里有精品视频免费| 成人18禁高潮啪啪吃奶动态图| 久久久久精品人妻al黑| 久久国产亚洲av麻豆专区| 99视频精品全部免费 在线| 99久国产av精品国产电影| 交换朋友夫妻互换小说| 久久久久国产精品人妻一区二区| 欧美精品人与动牲交sv欧美| 久久久久精品性色| 老司机影院毛片| 美国免费a级毛片| 插逼视频在线观看| 中文字幕精品免费在线观看视频 | 日本爱情动作片www.在线观看| 女人久久www免费人成看片| 性色avwww在线观看| 欧美精品一区二区大全| 国产精品国产三级国产专区5o| 美女主播在线视频| 久久精品国产鲁丝片午夜精品| 精品午夜福利在线看| 两个人免费观看高清视频| 国产亚洲欧美精品永久| 国产av一区二区精品久久| 免费久久久久久久精品成人欧美视频 | 超碰97精品在线观看| 成人黄色视频免费在线看| 自线自在国产av| 亚洲精品色激情综合| 日韩一区二区视频免费看| 黄片播放在线免费| 国产成人精品在线电影| www日本在线高清视频| 汤姆久久久久久久影院中文字幕| 欧美成人精品欧美一级黄| 女人精品久久久久毛片| h视频一区二区三区| 人人澡人人妻人| 夫妻性生交免费视频一级片| 国产老妇伦熟女老妇高清| 日韩av在线免费看完整版不卡| 婷婷色av中文字幕| 亚洲国产精品专区欧美| h视频一区二区三区| 综合色丁香网| 最近中文字幕2019免费版| 黄色视频在线播放观看不卡| 国产 一区精品| 精品人妻一区二区三区麻豆| 有码 亚洲区| 亚洲精品美女久久久久99蜜臀 | 精品人妻一区二区三区麻豆| 国产一区二区在线观看av| av在线app专区| 国产精品秋霞免费鲁丝片| 韩国精品一区二区三区 | 欧美日韩亚洲高清精品| 国产极品粉嫩免费观看在线| 十八禁高潮呻吟视频| 99香蕉大伊视频| 久久久久久人妻| 搡女人真爽免费视频火全软件| 97在线人人人人妻| 成人国语在线视频| 夜夜爽夜夜爽视频| 日韩中文字幕视频在线看片| 久热这里只有精品99| 中文字幕最新亚洲高清| 看十八女毛片水多多多| 久久久久网色| 超色免费av| 丝瓜视频免费看黄片| 国产男女内射视频| 在线精品无人区一区二区三| 美国免费a级毛片| 亚洲av国产av综合av卡| 自线自在国产av| 国产男女内射视频| 日韩,欧美,国产一区二区三区| 又黄又粗又硬又大视频| 欧美日韩成人在线一区二区| 99久久人妻综合| 一级毛片电影观看| 久久久国产精品麻豆| 久久久久精品性色| 中文字幕人妻熟女乱码| 日本欧美视频一区| 国产成人av激情在线播放| 成人亚洲欧美一区二区av| 最近的中文字幕免费完整| 久久精品国产自在天天线| 精品国产一区二区三区四区第35| 人人澡人人妻人| 一级毛片黄色毛片免费观看视频| 人妻系列 视频| 99热网站在线观看| kizo精华| 日韩欧美精品免费久久| 丝袜人妻中文字幕| 国产在线一区二区三区精| 男女高潮啪啪啪动态图| 色94色欧美一区二区| 亚洲欧美一区二区三区国产| 成人国产麻豆网| 久热这里只有精品99| 男男h啪啪无遮挡| 精品国产一区二区三区久久久樱花| 曰老女人黄片| 欧美精品人与动牲交sv欧美| 晚上一个人看的免费电影| 色婷婷av一区二区三区视频| 久久久欧美国产精品| 男男h啪啪无遮挡| 一级,二级,三级黄色视频| 国产深夜福利视频在线观看|