999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

Determination of interfacial tension and viscosity under dripping flow in a step T-junction microdevice

2022-04-01 05:47:04LiMaYongjinCuiLinShengChencanDuJianDengGuangshengLuo

Li Ma,Yongjin Cui,Lin Sheng,Chencan Du,Jian Deng,Guangsheng Luo

State Key Laboratory of Chemical Engineering,Department of Chemical Engineering,Tsinghua University,Beijing 100084,China

Keywords:Step T-junction Interfacial tension Viscosity Force balance

ABSTRACT Microfluidic approaches for the determination of interfacial tension and viscosity of liquid-liquid systems still face some challenges.One of them is liquid-liquid systems with low interfacial and high viscosity,because dripping flow in normal microdevices can’t be easily realized for the systems.In this work,we designed a capillary embedded step T-junction microdevice to develop a modified microfluidic approach to determine the interfacial tension of several systems,specially,for the systems with low interfacial tension and high viscosity.This method combines a classical T-junction geometry with a step to strengthen the shear force further to form monodispersed water/oil(w/o)or aqueous two-phase(ATP)droplet under dripping flow.For systems with low interfacial tension and high viscosity,the operating range for dripping flow is relative narrow whereas a wider dripping flow operating range can be realized in this step Tjunction microdevice when the capillary number of the continuous phase is in the range of 0.01 to 0.7.Additionally,the viscosity of the continuous phase was also measured in the same microdevice.Several different systems with an interfacial tension from 1.0 to 8.0 mN·m-1 and a viscosity from 0.9 to 10 mPa·s were measured accurately.The experimental results are in good agreement with the data obtained from a commercial interfacial tensiometer and a spinning digital viscometer.This work could extend the application of microfluidic flows.

1.Introduction

The property of interfacial tension is one of the most important physical properties for liquid-liquid systems.This property could greatly determine the phase separation and droplet coalescence[1],and therefore influence mass/heat transfer [2-4],multiphase chemical reaction[5,6],functional materials[7,8]and pharmaceutical materials preparation [9,10],biomedicine and biochemical analysis [11-13],protein crystallization [14],and many other applications.Over the past two decades,various methods have been developed to measure the interfacial tension of liquid-liquid systems,such as the direct measurement technique using a microbalance and some typical techniques based on capillary pressure,capillary gravity forces or gravity-distorted drops[15].Unfortunately,all the equipment are relatively complicated and hard to operate.Especially,it is not very reliable for the systems with low or ultralow interfacial tension.Therefore,more facile and reliable methods for the measurement of interfacial tension are required.With the growth of microfluidic technologies,many microdevices and approaches have been developed for the preparation of highly monodispersed microdroplets [16,17],including T-junction microchannels [18],flow-focusing microchannels [18],coaxial microchannels,and so on.Droplets generation methods mainly include cross-flowing shearing [19,20],perpendicular flowing shearing [21],hydrodynamic flow focusing [22,23],and geometry-dominated breakup [24].It has been proven that the droplet size can be controlled accurately in microdevices.Therefore,measuring the interfacial tension based on the droplet size[13],the droplet shape [25,26],and the pressure drop [25,27] has attracted widespread attention owning to the advantages of facile[28],reliable,high efficiency,and accuracy.

To measure the interfacial tension,the first microfluidic method was proposed and implemented by Hudson and co-workers[29,30],which is mainly based on the deformation and retraction dynamics of drops under an extensional flow in a microchannel.In this approach,the interfacial tension values can be drawn rapidly by calculating the rate change in the deformation of the droplet.Xuet al.[13]applied a coaxial microfluidic device to measure the interfacial tension based on the balance of the viscous shear force and interfacial tension during the droplet formation process.Nguyenet al.[31]proposed a microfluidic technique using a T-junction microdevice to measure some physical properties such as viscosity,surface tension,and interfacial tension.Additionally,some studies on test of interfacial tension and dynamic interfacial tension using microfluidic method by Luo’s group have been reported [6,32-35].All these microfluidic methods are considered to facile and convenient.But there are some limitations regarding the measurement of systems with low interfacial tension or high viscosity.In these situations,all the techniques are hardly applicable or is unsuitable very much.Additionally,compared with normal oil-water systems,the interfacial tension of aqueous twophase systems is extremely low,sometimes as low as 1-100 μN·m-1[36].The low interfacial tension with w/w systems tends to form all-aqueous jets rather than droplets due to the slow growth of the Rayleigh-Plateau instability [36].Thus,it is a challenge to generate aqueous two-phase (ATP) droplets in microdevices.Therefore,microfluidic measurement techniques with more flexible and much wider applicability must be further developed.

At present,the microfluidic approaches used to measure the interfacial tension are mostly operated under liquid-liquid dripping flow,because compared with jetting flow,the generated droplets under dripping flow can be more stable with a highly uniform size.In this way,the measurement accuracy can be guaranteed.To the best of our knowledge,microdroplets with uniformity and monodispersity can be generated under a dripping flow for normal liquid-liquid systems in T-junction,coaxial and flow-focusing microchannels [18].However,the operating range for generating microdroplets under a dripping flow is relatively narrow [37,17],and especially,liquid-liquid dripping flow is even hard to realize for systems with a low interfacial tension and high viscosity,such as ATP system.Therefore,it can be known that the limitations of all the reported methods are mainly come from the structures of microdevices.Therefore,specially designed microdevices must be developed and applied.

In this work,we attempt to measure the interfacial tension and viscosity of liquid-liquid systems by using a novel capillary embedded step T-junction microchannel,which was specially designed by Luo’s group [37].In this new microdevice an asymmetric flow field can be formed,and it can strengthen the shear force in the microdispersion process.Thus,liquid-liquid microdispersion process under dripping flow can be easily realized.In this way we may developed a more reliable microfluidic measurement technique,which could have more applicability for different systems in a much wider operating range.More than 8 liquid-liquid systems including water/oil (w/o) systems and an ATP system were selected as the working systems.A liquid-liquid microflow pattern and a mathematical model under dripping flow to calculate the interfacial tension and viscosity were provided.The experimental results were compared with the values measured by commercial instruments to valid the reliability of the microfluidic technique.In addition,the advantages of the new microfluidic technique were summarized.

2.Experimental

2.1.Materials and reagents

1-Butanol (≥99%) andn-hexane (≥97.0%) were obtained from Greagent.Chitosan with degree of deacetylation 80%-95%was purchased from Sinopharm Chemical Reagent Co,Ltd,n-octanol(99%),and Caprolactam (99.5%) were purchased from Shanghai Meryer.H2SO4(AR)was purchased from Sinopharm Chemical Reagent Beijing Co,Ltd.NH3·H2O was purchased from Shanghai Titan Scientific Co,Ltd.Span 80 used as the surfactant dissolved inn-octanol purchased from Shanghai Meryer.Deionized water was obtained from an ultrapure water machine (Center120FV-S,ZEBULA Technology,China).The composition and physical properties of liquid-liquid systems in this work are listed in Table 1.The viscosity and density of different liquid-liquid systems were measured by a commercial spinning digital viscometer (DV-IIa+P,Brookfield,USA) and an electronic densimeter (LC-MDJ-600G,LICHEN Technology,China),respectively.The interfacial tension of the liquid-liquid systems was measured by a commercial interfacial tensiometer(OCAH200,DataPhysics Instruments,Germany) based on pendant drop technology.

Table 1Physical properties of the liquid-liquid system (25 °C)

Table 2The experimental values of the interfacial tension σ for different systems

2.2.Fabrication of the microdevice

In this work,A capillary embedded step T-junction microchannels was specially designed and used to generate microdroplets,which was fabricated on a polymethyl methacrylate (PMMA) chip with a size of 20 mm×20 mm×3 mm.The main upstream channel was 400 μm in width(W)and 365 μm in height(h).In addition,the size of the main downstream channel and the side channel both in width and in height were 365 μm.A quart capillary with an outer diameter of 365 μm and an inner diameter of 75 μm was embedded in the side channel and extended into the main channel.After the microdevice was fabricated,it was sealed using another PMMA chip at 75°C for 3 min using a high-pressure thermal sealing technique.The structures of the microdevice are shown in Fig.1.

2.3.Experimental schematic

The experimental schematic diagram is shown in Fig.2.The dispersed phase was injected into the side channel by a syringe pump(LSP01-1BH,Longer,China).The continuous phase was injected into the main channel by a syringe pump (LSP01-1BH,Longer,China) simultaneously.The surfactants of Span 80 and Tween 20 were added into the continuous phase and dispersed phase to reduce the interfacial tension and stabilize the droplets.The surfactant concentrations were all kept at 2.0% (mass) in the experiment,which is much higher than the critical micelle concentration to eliminate the mass transfer influence.To minimize the measurement error,the flow rates of the two pumps were calibrated using deionized water.The droplet generation process was carried out with a microscope connected to a high-speed CCD video camera(i-SPEED TR,Olympus,Japan).The images were recorded with different frame rates adjusted to match the droplet generation frequency.The droplet size was measured by i-speed software.The entire process of droplet formation requires 30 s of equilibration time after changing the flow parameters.To ensure the measurement accuracy,100 microdroplets were measured to determine the average size.All experiments were conducted at room temperature and atmospheric pressure.For the measurement of the viscosity of the continuous phase,different mass concentrations of chitosan were added to an aqueous solution of 2%(mass) acetic acid.

Fig.1.Structures of the capillary embedded step T-junction microchannel.

Fig.2.The experimental schematic diagram.

2.4.The measurement fundamentals and mathematical model

During the liquid-liquid microdispersion process under dripping flow,the interfacial tension and viscous shear force are the dominant forces [38-40].Other forces,such as inertial force and buoyancy force can be neglected.In this capillary embedded step T-junction microchannel,there are two typical flow regimes:dripping and jetting.In the process of droplet formation,compared with jetting flow,the generated droplets under dripping flow are more stable,more uniform and highly monodispersed.In addition,droplet generation is mainly affected by the flow rates,the viscosity of the two phases,and the interfacial tension.

Under dripping flow,the dynamics of droplet breakup are dominated by the force balance between the interfacial tension and the viscous force [38-40],as shown in Fig.3.Therefore,the equation during droplet generation can be written as Eq.(1):

whereFDrepresents the viscous shear force during the process droplet formation induced by continuous phase flow,andFσ represents the interfacial tension between the two phases.

2.4.1.Viscous shear force

According to the literature,the values of the viscous shear force can be approximately expressed by Eq.(2)at low Reynolds numberRe[38,39,41,42]:

whereuc-udis the relative local velocity of the continuous phase flowucwith respect to the velocity of the dropud,ddrepresents the average diameter of the droplets,and μcis the velocity of the continuous phase.It is worth note that,the Reynolds number of continuous phase is in the range of 0.23 <Rec<58.44 in experiment.For simplicity,we can use the average velocity of the continuous and dispersed phase flowto replace the local velocity,then,Eq.(2) can be replaced by Eq.(3):

wherekDis a constant for a fixed geometry microdevice.The values ofcan be calculated by Eqs.(4) and (5):

whereQdandQcin Eqs.(4)and(5)are the dispersed and continuous phase volume flow rates,respectively.WandHare the width and height of the main microchannel,respectively.

Fig.3.Force analysis of droplet generation.

2.4.2.Interfacial tension

The interfacial tension can be calculated by Eq.(6)according to the literature reported [38,41,42]:

where σ anddnare the interfacial tension and the inner diameter of the capillary,respectively.Eq.(6)based on the assumption that the Laplace pressure of the liquid droplet is at its maximum value when the diameter of the liquid drop curvature is equal to the capillary diameter,whereas in the actual situation,there will exist a departure of the interfacial from the rim of the capillary during the detachment process.Therefore,the interfacial tension is expected to be greatly reduced,and Eq.(6) is no longer invalid to predict the interfacial tension.Then,Eq.(6) was modified based on the equilibrium droplet diameter instead of the dispersed phase channel width by Husny and Cooper-White [43] and can be further described by Eq.(7):it should be noted that the applicable range for measuring the interfacial tension and viscosity based on Eq.(3) and Eq.(7) is under dripping flow.In our experiment,the measurement of interfacial tension and viscosity were carried out under dripping flow that is when the capillary number of the continuous phase is in the range of 0.01 to 0.7 as shown in Fig.4.

Combining Eqs.(1),(3) and (7),Eq.(8) can be obtained:

Furthermore,we can obtain a mathematical model to predict the interfacial tension and viscosity:

wherekis a parameter of the fixed microchannel.In Eq.(9),the value ofkmainly depends on the geometry of the microchannel,the degree of confinement,the surfactant concentration distributions and the Marangoni effects in droplet flow [6,13,44].Since the surfactant concentration far exceeds the critical micelle concentration,thus,the adsorption of surfactant to the interfacial surface could be finished in a short time,and surfactant distributions and Marangoni effects could be neglected in this work[6,33].Moreover,the wetting condition can also be neglected in this step T-junction microchannel.Thus,kcan be regarded as a constant for a fixed microchannel.Then,according to Eq.(9),kcan be easily obtained by a standard system if we know the value μcof the continuous phase and σ of the system.

Moreover,the capillary numberCacof the continuous phase and Weber number of the dispersed phase (Wed) can be described by Eq.(10) and Eq.(11),respectively.

where μcis the viscosity of continuous phase,ρ andlin Eq.(11)are density and characteristic size of the dispersed phase,ucandudcan be calculated by Eq.(12) and Eq.(13).

Combing Eqs.(9)and(10),the droplet size can be approximated by a power law that is described as Eq.(14).

Then,according to Eq.(9),the interfacial tension can be gained if the viscosity of the continuous phase and the droplet size are known.If the interfacial tension and the droplet size are known,then the viscosity of the continuous phase can be calculated.

3.Results and Discussion

3.1.Flow regime in the step T-junction microchannel

As illustrated in Fig.4,the liquid-liquid microflow regime can be mainly divided into squeezing,dripping and jetting flow.Compared with jetting flow,dripping flow has a much stronger stability and the size of the generated droplets has a very narrow distribution.Therefore,we carried out the measurements under the dripping regime in this work.Dripping flow can be observed when 0.01 ≤Cac≤0.70 in this step T-junction microdevice.Because the asymmetric flow field with this microdevice strengthen the shear force during the droplet formation,stable dripping flow can be easily realized even with low interfacial tension systems as shown in Fig.4.The experimental operating range is much larger than that of other microdevices,such as flow-focusing and Tjunction microdevices [43,45].The wider operating range facilitates the measurement of the interfacial tension in a wider range,even for systems with low interfacial tension.

Fig.4.The flow pattern in this microdevice.

Fig.5.Process of droplet generation (continuous phase:deionized water-saturated butanol,dispersed phase:deionized water; Qc=100 μl·min-1, Qd=10 μl·min-1).

3.2.Generation of the monodisperse droplets

Fig.6.Distribution of droplet size (continuous phase:deionized water-saturated butanol,dispersed phase:deionized water; Qc=100 μl·min-1, Qd=10 μl·min-1).

Fig.7.Effect of the continuous phase flow rate on the droplet size(Qd=5 μl·min-1).

Before the measurement,it is very important to determine the flow pattern between the two phases,which determines whether Eq.(9) can be used to measure the interfacial tension within the operating range.To ensure an accurate relationship between the average size of droplet and the interfacial tension,n-butanol/deionized water was served as the controlled experimental system because of constant interfacial tension.The flow rate of the two phases is a significant parameter in liquid-liquid microdispersion.By changing the two-phase flow rate,monodispersed microdroplets with high uniformity under a dripping flow can be generated as depicted in Fig.5.Dripping flow can be realized when thedispersed phase flow rate was 5-50 μl·min-1and the corresponding continuous phase flow rate was 50-1400 μl·min-1.The coefficient of variation (CV) of the droplet size is defined as CV=δ/dav×100%,where δ is the standard deviation anddavis the average droplet diameter.As shown in Fig.6,through the statistics of 100 droplets diameter in the main channel near the step structure,the CV value is less than 5%,thus,the droplet generation processes can repeat quite well.The Stable droplets with highly monodispersity,uniformity generation in a wider operating range under dripping flow can provide the fundamentals for the measurements of interfacial tension and viscosity.

Fig.8.Determination of the value of k, (d d/dn)2 is a function ofthe slope is equal to the parameter k.

Fig.9.Effect of the continuous phase flow rate on the droplet size for the different systems.(Qd=5 μl·min-1).

3.3.Determination of the value of parameter k

The standard system,deionized water/1-butanol,was used to determine the value of parameterkof the capillary embedded step T-junction microchannel.It is worth noted that 1-butanol and water are mutually saturated which excludes mass transfer during the droplet formation process[46].Fig.7 shows that under a given dispersed phase flow rate of 5 μl·min-1,the droplet size decreases with the increase of the continuous phase flow rates.It is because an increase in the continuous phase flow rate means an increase in the shear force,which leads to a continuous decrease in the droplet size.

Combined with Eq.(9),according to the linear fit ofto(dd/dn)2,the slope is the value of para meterk,as shown in Fig.8,and the value of parameterkis 0.483 for the capillary embedded step T-junction microchannel.

3.4.Measurement of interfacial tension σ

Several liquid-liquid systems and an ATP system with different interfacial tension were selected as the working systems.Firstly,four w/o systems(2%(mass)Span 80 inn-octanol/water,2%(mass)Span80n-hexane/water,n-octanol/water,andn-hexane/2%(mass)Tween 20 in water) with constant interfacial tension were tested.The organic phase was used as the continuous phase,and the aqueous phase was used as the dispersed phase in these working systems.Fig.9 represents the effects of the continuous phase flow rates on the droplet size at a given dispersed flow rate of 5 μl·min-1.The trends of the four systems are similar when the continuous phase flow rates range from 50 to 800 μl·min-1.Then,Fig.10 shows(dd/dn)2as a function ofk/,and the slopes of the fitting lines are the interfacial tensions for different systems.The measured interfacial tensions of the four systems are (4.43 ± 0.04),(2.94 ± 0.10),(8.20± 0.30),and (1.68 ± 0.05) mN·m-1,respectively.

Fig.10.Determination of the values of σ: (d d/dn)2 as a function ofThe slope is equal to the interfacial tension σ(mN·m-1).((a):black square,2%(mass)Span 80 in n-octanol/water;red circle, n-octanol/water;(b):black square,2% (mass) Span80 n-hexane/water;red circle, n-hexane/2% (mass) Tween 20 in water).

Fig.11.(a) Effect of the continuous phase flow rate on the droplet size (Qd=5 μl·min-1,ε-caprolactam aqueous solution as the continuous phase and (NH4)2SO4 aqueoussolution as the dispersed phase);(b) Determination of the values of σ: (d d/dn)2 as a function ofThe slope is equal to the interfacial tension σ (mN·m-1).

Then,Caprolactam-(NH4)2SO4-water with mutually saturated could form liquid-liquid two phase system,and was selected as an ATP system to prove the effectiveness and feasibility of this microfluidic method.The microdroplets also can be generated under a dripping flow in the step T-junction microdevice,and then the interfacial tension can be obtained.Fig.11(a)shows the droplet size of the aqueous two-phase system.According to Eq.(9),(dd/dn)2was fit to a linear function ofand the value of the slope of (2.48 ± 0.05) mN·m-1is the value of the interfacial tension of the system (Fig.11(b)).

In Table 2,the experimental data were compared with the data acquired from commercial interfacial tensiometer.It can be seen they are in good consistency.Additionally,the experimental values are in good agreement with the literature reported[13,35].Therefore,the effectiveness and accuracy of this equipment to measure interfacial tension have been further proved.

Fig.12.Effect of the continuous phase flow rate on the relative droplet size for the different systems.

Subscript mea represents the values of the interfacial tension σ were measured by the commercial pendant-droplet interfacial tensiometer

3.5.Measurement of the continuous phase viscosity

To determine the viscosity of the continuous phase,several systems were selected for testing in this work.First,2%(mass)span 80 inn-hexane/deionized water was initially used as the standard experimental system to ensure an accuracy and fundamentals of other system’s measurements.Then,we chose chitosan aqueous solution with different mass concentrations used as the dispersed phase andn-octanol,with 2% (mass) Span 80 added to stabilize the droplets,as the continuous phase to prove the repeatability of the method.Fig.12 illustrates that the microdroplets size varies with the continuous phase flow rates at a given dispersed flow rate of 5 μl·min-1.Similarly,with increasing continuous phase flow rates,the size of the microdroplets decreases,according to Eq.(9),is a function of(dd/dn)2,and the slope is equal to μc(mPa·s).The measurement results are shown in Fig.13(a).Eventually,we chose 2% (mass) span80 inn-octanol/deionized water andn-octanol/deionized water to confirm the applicability and suitability of this microfluidic method and the measurements results are shown in Fig.13(b).

Eventually,we compared the experimental results with the data obtained from the commercial spinning digital viscometer and the results are given in Table 3.From Table 3 we can conclude that both are in good agreement.

Subscript mea represents the values of the viscosity μ were measured by the spinning digital viscometer.

Through this microfluidic equipment,the interfacial tension and the viscosity of different systems were tested with high accuracy,good repeatability and a small sample volume.The proper operating condition range for interfacial tension in this microdevice have been provided and have compared the operating range with other microdevices,the comparison results are shown in Table 4.

Fig.13.Determination of the viscosity μ:kσ/(-) as a function of (d d/dn)2.The slope is equal to μc(mPa·s).((a) blue triangle,red circle,black square:represents three systems with same continuous phase that 2%(mass)Span80 in n-octanol,1.2,0.8,0.4%(mass)chitosan aqueous solution used as dispersed phase;green triangle represents 2%(mass)Span80 in n-octanol/deionized water system;(b)black square,n-octanol/deionized water;red circle,2%(mass)span80 in n-hexane/deionized water;blue triangle,n-hexane/2% (mass) Tween 20 in deionized water).

Table 3The experimental values of the continuous phase viscosity for the different systems

Table 4The proper operating range of different microdevices

4.Conclusions

In this work,we designed a novel microfluidic device that is a capillary embedded step T-junction microchannel to measure the interfacial tension of two immiscible liquids including several different w/o systems and an ATP system.Additionally,the viscosity of the continuous phases was also tested using the same microdevice.All the measurements were based on the force balance between the interfacial tension and the viscous shear force during the droplet generation under the dripping regime.The fabrication of the microdevice and the mathematical derivation process of the measurement fundamentals are provided in this work.Several systems with interfacial tension from 1.0 to 8.0 mN·m-1and the viscosity of the continuous phase from 0.9 to 10 mPa·s were accurately measured with this microdevice.Then we compared the experimental data with data acquired from the interfacial tensiometer and the spinning digital viscometer,the comparison results show that they are in good agreement.Finally,the accuracy,effectiveness,applicability and suitability of this microfluidic approach used to test low interfacial tension and viscosity have been proved.

CDediT Authorship Contribution Statement

Li Ma:Writing -original draft,Investigation,Methodology.Yongjin Cui:Investigation,Validation.Lin Sheng:Validation,Formal analysis.Chencan Du:Writing -review &editing,Conceptualization.Jian Deng:Writing - review &editing,Conceptualization.Guangsheng Luo:Writing -review &editing,Conceptualization,Supervision.

Declaration of Competing Interest

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.

Acknowledgements

This work was financially supported by the National Natural Science Foundation of China (21991104).

Nomenclature

Caccapillary number of the continuous phase

dddiameter of droplet,μm

dninner diameter of capillary

FPcinertia force of the continuous force,N

FPdinertia force of the dispersed phase,N

Hheight of the microchannel,μm

Qcvolumetric flow rate of the continuous phase,μl·min-1

Qdvolumetric flow rate of the dispersed phase,μl·min-1

ucsuperficial velocity of the continuous phase,m·s-1

udsuperficial velocity of the dispersed phase,m·s-1

Wwidth of the microchannel,μm

Wedweber number of the dispersed phase

wwidth of the orifice,μm

μ viscosity of different system mPa·s

σ interfacial tension between two phases,mN·m-1

主站蜘蛛池模板: 天堂av高清一区二区三区| a在线亚洲男人的天堂试看| 免费观看成人久久网免费观看| 麻豆精品在线| 东京热一区二区三区无码视频| 国产中文一区a级毛片视频| 国产在线一二三区| 日韩精品视频久久| 国产精品思思热在线| 18黑白丝水手服自慰喷水网站| 久久这里只有精品2| 午夜性刺激在线观看免费| 一区二区三区四区日韩| 亚洲高清在线天堂精品| av一区二区无码在线| 国产午夜精品一区二区三| 国产在线拍偷自揄观看视频网站| 亚洲欧美在线综合图区| 激情综合五月网| 麻豆AV网站免费进入| 亚洲91精品视频| 欧美日本在线一区二区三区| 欧美一区日韩一区中文字幕页| 欧美精品伊人久久| 99这里只有精品在线| 国产精品性| 国产毛片片精品天天看视频| 国产91无码福利在线| 国产精品视频白浆免费视频| 亚洲AV无码久久精品色欲| a级毛片免费网站| 成人精品视频一区二区在线| 国产在线98福利播放视频免费| 91网在线| 成年人免费国产视频| 无套av在线| av无码久久精品| 成人在线第一页| 国产女人喷水视频| 成人一区在线| 色婷婷国产精品视频| 色视频国产| 午夜福利视频一区| 日韩国产综合精选| 日韩专区欧美| 国产高清在线丝袜精品一区| 国产女人在线观看| av一区二区无码在线| 欧美日韩精品一区二区在线线| 在线国产毛片手机小视频| 日本三级精品| 精品国产自| 九色视频线上播放| 超碰aⅴ人人做人人爽欧美| 成人午夜天| 99国产精品一区二区| 色婷婷亚洲十月十月色天| 成人久久精品一区二区三区| 亚洲无码视频一区二区三区| 午夜人性色福利无码视频在线观看| 国产自视频| 国产激情第一页| 国产精品浪潮Av| 成人福利在线视频免费观看| 91九色国产porny| 免费人成网站在线观看欧美| av在线5g无码天天| 97se亚洲综合不卡| 色综合久久88色综合天天提莫| 手机精品视频在线观看免费| 无码高潮喷水专区久久| 55夜色66夜色国产精品视频| 久久黄色视频影| 日日噜噜夜夜狠狠视频| 蜜臀av性久久久久蜜臀aⅴ麻豆| 日韩一区二区三免费高清| 小说区 亚洲 自拍 另类| 波多野结衣一区二区三视频| 国产日韩欧美在线视频免费观看| 日韩不卡高清视频| 一级高清毛片免费a级高清毛片| 亚洲自拍另类|