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

    Parametric study of high-frequency characteristics of plasma synthetic jet actuator

    2022-01-10 14:51:42GuozhengSONG宋國正HaohuaZONG宗豪華HuaLIANG梁華ZhiSU蘇志LikeXIE謝理科andXingZHENG鄭猩
    Plasma Science and Technology 2021年12期
    關(guān)鍵詞:蘇志理科

    Guozheng SONG (宋國正), Haohua ZONG (宗豪華), Hua LIANG (梁華),Zhi SU (蘇志), Like XIE (謝理科) and Xing ZHENG (鄭猩)

    1Science and Technology on Plasma Dynamics Laboratory, Air Force Engineering University, Xi’an 710038, People’s Republic of China

    2 School of Mechanical Engineering, Xi’an Jiaotong University, Xi’an 710049, People’s Republic of China

    Abstract A major issue of plasma synthetic jet actuator (PSJA) is the severe performance deterioration at high working frequency.In this study, experiments and numerical simulation are combined together to investigate the influence of thermal conductivity, throat length (Lth) and discharge duration (Td) on the high-frequency characteristics of PSJA.Results show that the variation of the actuator thermal conductivity and discharge duration will not alter the saturation frequency of the actuator, whereas decreasing the throat length results in an increase of the saturation frequency.For a short-duration capacitive discharge of 1.7 μs, a clear shock wave is issued from the orifice, followed by a weak jet.As a comparison, when the discharge duration is increased up to 202.6 μs, a strong jet column is formed and no obvious shock wave can be visualized.Based on numerical simulation results, it becomes clear that the long-duration pulse-DC discharge is able to heat the cavity gas to a much higher temperature (3141 K) than capacitive discharge, greatly improving the conversion efficiency of the arc discharge energy to the internal energy of the cavity gas.In addition, high-speed Schlieren imaging is deployed to study the performance degradation mechanism of PSJA at high working frequency.Monitor of the exit jet grayscale indicates that as long as the saturation frequency is exceeded, the actuator becomes unstable due to insufficient refresh time.The higher the discharge frequency, the more frequently the phenomenon of ‘misfires’ will occur, which explains well the decaying jet total pressure at above saturation frequency.

    Keywords: thermal conductivity, throat length, discharge duration, plasma synthetic jet,saturation frequency

    1.Introduction

    Plasma flow control(PFC)is a novel and promising active flow control technology, having the advantages of simple structure,fast response, wide working bandwidth, no mechanical moving parts, etc [1–3].Surface dielectric barrier discharge actuator(SDBDA) [4–6], arc discharge plasma actuator (ADPA) [7],plasma synthetic jet actuator (PSJA) [8, 9] are the three main actuation methods of PFC.Among them, plasma synthetic jet can generate high-energy jets with speeds up to hundreds of meters per second,and shock waves with approximate speed of sound.At present, plasma synthetic jet has demonstrated sufficient control capabilities in lift augmentation [10, 11], jet noise control [12], and shock wave/boundary layer interaction[13, 14].To realize better flow control effect at high-Reynolds number flow, it is necessary to conduct meticulous research on frequency characteristics of plasma synthetic jet actuator.

    Since Grossmanet al[15] proposed the first plasma synthetic jet actuator in 2003, PSJA has received growing attention from the active flow control community.According to previous researches, frequency characteristics of plasma synthetic jet are mainly affected by the geometric parameters of the actuator,which consist of the electrode distance,orifice diameter, throat length, and cavity volume.Cybyket al[16]used numerical simulations to study the performance characteristics (single-pulse) of PSJA with different orifice size,cavity volume, and energy deposition.They found that the average orifice exit velocity increases with the growth of energy deposition and orifice diameter, decreasing with increasing cavity volume.Popkinet al[17]focused on novel actuator designs that could potentially improve the reliability,efficiency and lifetime of individual actuators.Based on their result, larger electrode distances lead to higher peak cavity pressure efficiency up to 75%,and by increasing the electrode diameter, both electrode endurance and discharge reliability can be enhanced.

    Previous researches mainly focused on the influence of geometrical parameters and atmospheric parameters on the performance of PSJA, and the impact of actuator material’s thermal conductivity and discharge duration on frequency characteristics of PSJA remains unknown.What’s more, the performance degradation mechanism of PSJA at high frequency needs to be revealed for purpose of actuator optimization.This study fills these knowledge gaps.Three actuators made of different materials and three lids with different thickness are manufactured to study the influence of throat length and thermal conductivity on PSJA.The total pressure of PSJA is measured by a pressure transducer, and the working process of PSJA is studied in detail with a highspeed Schlieren system.Furthermore, numerical simulations are carried out to assist the analysis.The experimental setup and computational setup are introduced in section 2.Results and analyses are presented in section 3.

    2.Experimental setup

    2.1.Actuator and power supply

    As shown in figure 1(a),the plasma synthetic jet actuator is twoelectrode,consisting of ceramic lid,ceramic cavity,and tungsten needles.The lids and cavities of the actuators are made of three different ceramic materials,which are ZrO2(3 W/(m·K)),Al2O3(32 W/(m·K)), AlN (150 W/(m·K)).The cavity volume of the actuator is approximate226 mm ,3 the in-close electrode distance is 2.5 mm, the orifice is in the center of the ceramic lid,with a diameter of 1 mm,the range of throat length is 3 mm, 6 mm, and 9 mm.

    A sequential discharge (capacitive discharge) power supply is adopted to feed the plasma synthetic jet actuator,as shown in figure 2(a).A DC power supply (output voltage:U0=2 kV;power:8 kW),a resistor(R1,resistance:1.5 kΩ;power: 400 W), and an energy-storing capacitor (C1, capacitance:C1=0.1 μF)constitute the charging circuit.A highvoltage power supply (output voltage: 0–10 kV; width: 2 μs;output frequency: 1–4000 Hz) is used for breaking the electrode gap and controlling the discharge frequency.Immediately after breakdown, the energy stored in C1 is rapidly released into the discharge channel in the form of arc heating.Two high voltage diodes (D1 and D2) isolate the capacitivedischarge circuit and the trigger-discharge circuit from each other to ensure protection from back-current overloading.

    The DC power supply generates a relatively large current with high discharge frequency, producing a DC arc in the actuator and imposing threat of burning the actuator.In order to avoid DC arc and ensure higher working frequency, we connect a high voltage transistor switch in series into the charging circuit, using a DG535 to control the high voltage transistor switch and the high-voltage pulsed power supply with TTL signal 1 and signal 2 respectively.The timing diagram of signal 1 and signal 2 is displayed in figure 2(c).The low-level time of signal 1(15 μs)completely contains the high-level time of signal 2(10 μs)in single discharge period,meaning that the high voltage transistor switch is off and DC power supply cannot generate current into the discharge circuit during the whole process of arc discharge.With this circuit, the safe working frequency range of the plasma synthesis jet actuator is totally consistent with the output frequency of the high-voltage pulsed power supply without DC arc.In order to study the effect of discharge duration,we developed another sequential discharge (pulse-DC discharge)power supply, removing the energy-storing capacitor and connecting an inductor (L1, inductance:L1=40 mH) in series to the actuator discharge circuit.Because the inductor has the function of maintaining the current constant, after breakdown,the DC source forms a stable DC arc between the electrodes of the actuator, and the discharge duration can be controlled by the high voltage transistor switch.The timing diagram of signal 1 and signal 2 in this power supply is displayed in figure 2(d).

    2.2.Discharge characteristic measurement

    The electrical parameters acquisition system consists of a high voltage probe(Tektronix P6015A),a current probe(P6021A),and an oscilloscope (DPO4104B).The high voltage probe and current probe are set in P1 and P2 respectively to measure the discharge voltage and current of actuator, as shown in figure 2(a).The voltage and current waveforms pertaining to the ZrO2actuator obtained at a DC output voltage of 1.6 kV,a pulse voltage of 8 kV, and a repetition rate of 1000 Hz are shown in figure 3.Before breakdown, the voltage on the actuator increases with the output voltage of high-voltage pulsed power supply.When the voltage on the actuator reaches 7.7 kV, the air between the electrodes is ignited, and the voltage quickly drops.Meanwhile, the current goes up rapidly and reaches a maximum of 189 A at 0.8 μs.After breakdown, both the charging voltage and discharge voltage have obvious oscillations, resulting from impedance mismatch, which are second order transient responses to a fast voltage variation.The pulse width of the current is 1.7 μs,and the single pulse energy is integrated to be 32.8 mJ for figure 3,much less than the energy stored in capacitor C1 (136.4 mJ).

    The voltage and current waveforms of theAl2O3actuator with DC power supply outputting 1 kV and high-voltage power supply pulsed at 8 kV, and 1400 Hz are shown in figure 4.The average discharge current is 1.7 A, the pulse width of the current is 202.6 μs,and the single pulse energy is 20.4 mJ.In order to keep the single pulse energy on the actuator consistent, the capacitor charging voltage is set to 1.2 kV during the study of the influence of the discharge duration,and the single pulse energy is 19.8 mJ in capacitive discharge under 1400 Hz.

    2.3.Time-averaged total pressure measurement and Schlieren system

    Total pressure of plasma synthetic jet indicates the mechanical energy incorporated in the jet, and thus can reflect its authority to control external flow.A digital pressure transducer (DPT 6000, range: ?5 to 5 kPa, accuracy: 0.02% F.S)together with a thin total pressure probe was used to measure the total pressure at a distance of 1 mm from the exit orifice,as shown in figure 5.The pressure probe is connected to the pressure transducer with a 0.5 m long PVC tube.Coping with EMI generated by the arc discharge, the measures we have adopted include using two uninterruptible power supplies to drive PSJA and pressure transducer respectively; the highvoltage wire in the discharge circuit is reduced to the shortest;all signal wires and the high-voltage power supply are completely wrapped by shielding net;the high-voltage output wire and ground wire of the high-voltage power supply are twisted together to form a shielded twisted pair.After taking these measures, the influence of EMI on the pressure measurement is well eliminated.In each test,the sampling frequency of the pressure transducer is set as 100 Hz,and the sampling time is 2 s.Three sets of data are collected at each discharge frequency.From the variation of time-averaged total pressure with discharge frequency, the saturation frequency of the actuator can be directly picked up.

    Looking for a job in a foreign country is not easy, but I will never give up. After going through the past year, I become more confident in overcoming difficulties. The road is chosen by me, I have to keep going and never hesitate in front of cross of life. Having a goal of making an ideal career, I strongly believe that all of my efforts will bring me great success.

    A high-speed Schlieren system with Z-type light path is used for studying the performance of the actuator as shown in figure 6.A high-power xenon lamp is used as a continuous point light source.The vertical knife edge is chosen to change the amount of light entering the camera to improve the contrast of the captured image.The frame rate of CCD camera(Phantom v2512) is set as 100,000 fps, the exposure time is set as 2 μs,and the resolution 128×256 pixels.One pixel in the Schlieren image represents a physical distance of 0.22 mm.According to the Schlieren images, the working process of PSJA can be directly visualized.

    2.4.Numerical model and computational grids

    In order to understand the physical evolution process of PSJ in depth, the commercial software FLUENT was used to conduct a transient simulation of the energy deposition and jet ejection process in two-dimensional domain.The turbulence model is SSTk-ω model, and the calculation time step increases gradually from 10 ns at discharge energy deposition to 50 μs at jet expulsion.Based on a preliminary grid independence check, the number of grids is chosen as 210,000,above which further refinement would not bring any accuracy improvement.An adiabatic non-slip wall boundary condition and pressure-far-field boundary condition are set at the outer boundary, as shown in figure 7.The red part of the grids represents the energy deposition area,which is a circle with a diameter of 1 mm in the center of the cavity.In order to simulate the energy deposition process of the actuator, the power densityShof the energy deposition aera is added to the heat source term of the energy equation,and the definition ofShis shown in equation (1), whereηdenotes the energy conversion rate from electrical energy to heat (40%);Q0represents the single pulse energy(20 mJ);τpis the discharge duration;Vhis the total volume of heating region.

    3.Results and analyses

    3.1.The effect of thermal conductivity

    With the capacitive discharge power supply displayed in section 2, the variations of time-averaged total pressure with discharge frequency under different thermal conductivities are shown in figure 8.In order to better compare the variations,data presented are normalized by dividing the largest value of all these three curves.

    From figure 8(a), the time-averaged total pressure first increases and then decreases.The saturation frequency of PSJA is the discharge frequency corresponding to the total pressure peak.The saturation frequencies of AlN,Al2O3,and ZrO2actuators are 1800 Hz, 1700 Hz and 1800 Hz, respectively.In addition, the peak time-averaged total pressures of Al2O3and ZrO2actuators are 95.0% and 99.8% of that of AlN actuator.Theoretically, the greater the thermal conductivity of the actuator material, the greater the power of the actuator to dissipate heat to the outside atmosphere,while the electrical power of the actuators is basically the same, so the time-averaged total pressure of the actuator should decrease as the thermal conductivity increases.In our experiment, the time-averaged total pressure ofAlN actuator is larger than that ofAl2O3actuator.The reason may be due to the different positions of the total pressure probe relative to the orifice during the measurements.

    To compensate for the limitations of the experiment, we used the analytic model established by Zonget al[24] to study the influences of thermal conductivity and throat length on the frequency characteristics of the plasma synthetic jet actuator.Key thermal parameters set in the analytic model are shown in table 1.The single pulse energy is set as 33 mJ,the discharge duration is 2 μs, and the ratio of input energy to internal energy of the cavity air is set as 40% according to Popkinet al[25].Through analytic calculation,the variations of time-averaged total pressure with discharge frequency under different thermal conductivities are shown in figure 8(b).The saturation frequencies and total pressure of actuators with different thermal conductivities are basically the same.The effect of improving the thermal conductivity of the actuator material is similar to slightly reducing the input energy of the actuator,and the frequency characteristic of the actuator will not be changed.The peak time-averaged total pressures of ALN actuator and Al2O3actuator are 99.1%and 99.2% of that of ZrO2actuator.Higher thermal conductivity results in smaller time-averaged total pressure, but the difference is negligible.

    3.2.The effect of throat length

    With the capacitive discharge power supply,the variations of time-averaged total pressure with discharge frequency under different throat length are shown in figure 9(a).For a throat length of 3 mm, 6 mm and 9 mm, the value of the saturation frequency is determined to be 1700 Hz, 1200 Hz and 1100 Hz, respectively.The peak time-averaged total pressure of 6 mm actuator and 9 mm actuator is 81.2% and 59.5% of that of 3 mm actuator.

    Through the analytic model established by Zonget al[24], the variations of time-averaged total pressure with discharge frequency under different throat lengths are shown in figure 9(b).The saturation frequencies of 3 mm, 6 mm and 9 mm actuators are 1700 Hz, 1100 Hz and 900 Hz, respectively.The peak time-averaged total pressure of 6 mm actuator and9 mm actuator is 84.3% and 77.7% of3 mm actuator.

    Compared with experimental result, the calculated timeaveraged total pressure is larger,and the saturation frequency decreases more quickly with the throat length.The reason for the deviation between the analytic results and the experimental results is that the influence of gas viscosity in the throat length is ignored in the analytic model.The lowvelocity fluid in the throat will hinder the ejection of the highpressure gas in the cavity.Ignoring the effect of gas viscosity in the throat length will increase the jet velocity,make the jet velocity attenuate slower, and the actuator will enter the refresh stage later.The jet duration and time-averaged total pressure are correspondingly increased, making the actuator more sensitive to the change of discharge frequency and improving the speed of the saturation frequency changing with throat length.

    Table 1.Key thermal parameters set in the solver.

    3.3.The effect of discharge duration

    The pulse width of the current is 1.7 μs in capacitive discharge,and 202.6 μs in pulse-DC discharge.Variation of the time-average total pressure of with different discharge durations at different discharge frequencies is shown in figure 10.Although the time-averaged total pressure increases with the increase of the discharge duration, the saturation frequency remains approximately unchanged, both of which are 1800 Hz.The peak total pressure corresponding to capacitive discharge is only 26.8% of that in long-duration pulse-DC discharge.These observations can be explained as follows.PSJA is similar to a spring-mass-damper system, and its saturation frequency is dominant by the Helmholtz natural oscillation frequency (fh), which is related to both atmospheric parameters and geometrical parameters as[26, 27].Therefore, increasing the discharge duration does not change the saturation frequency of the actuator.The large total pressure of pulse-DC discharge at high frequency may be attributed to the increasing energy conversion efficiency of the actuator, and the corresponding mechanisms deserve further investigation.

    The working process of the actuator recorded by high-speed Schlieren system under different discharge duration is shown in figure 11.From the Schlieren images,it can be seen that the jet produced by pulse-DC discharge is stronger than that by capacitive discharge.In capacitive discharge, the actuator can produce obvious shock waves.As a comparison, no shock waves are observed in pulse-DC discharge.These Schlieren images are consistent with the pressure measurements.Under capacitive discharge,the gas in the cavity is heated rapidly and the spatial distribution of pressure and temperature is largely nonuniform,producing obvious shock wave and weak jet.However, as for pulse-DC discharge, the pressure and temperature of the gas increase more slowly and distribute more evenly, inducing stronger jet with negligible shock waves.

    Numerical simulations are further performed to understand the impact of discharge duration on inter-cavity heating.Figure 12 compares the average temperature in the cavity under different discharge duration, and the temperature distribution at different moments is shown in figure 13.At capacitive discharge,the gas in the cavity is rapidly heated to produce shock waves.The shock waves propagate back and forth in the cavity to dissipate heat and continuously take away a large amount of heat through the orifice to the outside,causing the temperature in the cavity to drop quickly.So obvious shock waves and a weak jet are formed one by one at the orifice.Compared to capacitive discharge, there is no obvious shock wave formed at pulse-DC discharge, and the gas in the cavity is gradually heated to a pretty high temperature (3141 K), greatly improving the efficiency of converting the heat generated by the arc into the internal energy of the gas in the cavity.Therefore, stronger jet is generated with large discharge duration.

    Figure 1.The configuration of plasma synthetic jet actuator.(a) Schematic diagram, (b) physical map.

    Figure 2.The power supply system and trigger signal.(a)Capacitive discharge, (b) pulse-DC discharge, (c) signal timing diagram in (a),(d) signal timing diagram in (b).

    Figure 3.Representative voltage and current waveforms obtained in capacitive discharge.

    Figure 4.Representative voltage and current waveforms obtained in pulse-DC discharge.

    Figure 5.The total pressure measurement system.

    Figure 6.The high-speed Schlieren system.

    Figure 7.Computational mesh.

    Figure 8.The variations of time-averaged total pressure with different thermal conductivities.(a) Experimental result, (b) analytic result.

    Figure 9.The variations of time-averaged total pressure with different throat lengths.(a) Experimental result, (b) analytic result.

    Figure 10.The variations of time-averaged total pressure with different discharge durations.

    Figure 11.The working process of the actuator at 1800 Hz under different discharge durations.(a)Capacitive discharge,(b)pulse-DC discharge.

    Figure 12.The average temperature in the cavity under different discharge durations.

    Figure 13.The temperature distributions in the cavity at different moments.

    Figure 14.The Schlieren images processing.(a) Monitoring area, (b) fd=1800 Hz, (c) fd=2200 Hz, (d) fd=2600 Hz.

    Figure 15.Variations of the mean grayscale in one period at increasing frequencies.(a) Capacitive discharge, (b) pulse-DC discharge.

    Figure 16.The misfire ratio at different working frequencies.

    3.4.The mechanism of the performance degradation at high frequencies

    To study the mechanism of the performance degradation at high frequencies,the mean grayscale of a 3 mm×3 mm area near the orifice in the Schlieren image is monitored, and the monitoring area is shown in figure 14(a).When the 3 mm actuator works at different discharge frequencies with pulse-DC discharge, the mean grayscale varies with time as shown in figures 14(b)–(d).Before reaching the saturation frequency,the mean grayscale curves have obvious periodicity, and the peak-to-peak values of the mean grayscale curves increase as the frequency increases.After the saturation frequency is exceeded, the periodicity of the mean grayscale curves decreases,and the peak-to-peak value decreases rapidly as the frequency increases.Under pulse-DC discharge and the discharge frequency of 2200 Hz,the mean grayscale curve has a significant decrease in amplitude after 0.58 s from the start of discharge, while it has a significant decrease in amplitude after 0.13 s after the beginning of discharge at a discharge frequency of 2600 Hz.In addition, when the discharge frequency is 2600 Hz, due to insufficient discharge stability,the gas between the electrodes is not ionized during some discharge cycle, which results in enough refresh time for PSJA,causing obvious fluctuations in the mean grayscale curve.These results show that when the actuator is working at high frequencies, due to the reduction of the refresh time, the higher the working frequency,the sooner the jet generated by the actuator will be weakened, and the performance of the actuator degrades rapidly.

    The instantaneous grayscale variation in figure 14 is averaged over multiple cycles, resulting in a phase-averaged variation of the jet grayscale in one period, as shown in figure 15.The phase-averaged result shows the process of PSJA’s working performance changing with the discharge frequency, which is in good agreement with the pressure measurements.As the discharge frequency exceeds the saturation frequency, the peak values of mean grayscale curves decrease rapidly, and PSJA exhibits obvious degradation in its authority.

    Previous analyses of the Schlieren images indicate that as long as the peak-to-peak value of grayscale variation in one cycle is less than 1.5, no obvious jet can be observed at the orifice.As such, we use this threshold grayscale variation(1.5) to evaluate whether the pulsed jet is fired normally in one cycle or abnormally as‘misfire’.Further,the misfire ratio at different working frequencies is calculated, as shown in figure 16.When the working frequency does not exceed the saturation frequency, the ratio of ‘misfires’ cycle to the total cycle is 0,and the jets appear periodically at the orifice of the actuator.As the working frequency reaches 2200 Hz,the ratio increases to 51.5%under capacitive discharge,and remains 0 under pulse-DC discharge.When the working frequency is improved to 2600 Hz, the ratio further increases to 67.4%under capacitive discharge, and 65.3% under pulse-DC discharge.So the performance of PSJA degrades rapidly as discharge frequency exceeds the saturation frequency,and the phenomenon of ‘misfires’ becomes more frequent.

    4.Conclusion

    In the present paper, experimental investigation and numerical simulation were conducted to study the effects of the thermal conductivity,throat length and discharge duration on PSJA.Furthermore, the mechanism of the performance degradation of PSJA at high working frequencies is analyzed in detail.

    It is found that the saturation frequency of PSJA is not affected by the thermal conductivity and discharge duration,but decreases with the increase of the throat length.This verifies that the saturation frequency of PSJA is determined by its own Helmholtz natural oscillation frequency(fh),which is related to both atmospheric parameters and geometrical parameters.With different thermal conductivities, the total pressure near the orifice remains basically unchanged.While the discharge duration is larger (202.6 μs), the total pressure increases significantly,and its peak reaches 3.7 times that under the discharge duration of 1.7 μs.The larger the discharge time, the longer the arc heating the air in the cavity, the more uniform the gas temperature and pressure distribution in the cavity, and no shock wave will be generated.But the temperature in the cavity will gradually rise to a higher level, and a stronger jet will be produced.In addition, as the working frequency exceeds the saturation frequency, due to the reduction of the refresh time,the higher the working frequency, the sooner the jet generated by the actuator will be weakened, the more frequent the phenomenon of ‘misfires’, and the performance of the actuator degrades rapidly.

    In future study, the effect of the discharge duration on PSJA should be more in-depth researched.As the discharge duration further increases,it is likely that the reduction of the refresh time will reduce the saturation frequency and the intensity of the jet.In addition, the saturation frequency is determined by the geometrical parameters of PSJA.The configuration of the actuator deserves to be further optimized to improve its performance.

    Acknowledgments

    This research work is supported by National Natural Science Foundation of China (No.12002384), the National Key Laboratory Foundation of China (No.614220210200112), and the Foundation Strengthening Project (No.2019-JCJQ-JJ-077).

    ORCID iDs

    猜你喜歡
    蘇志理科
    Airfoil friction drag reduction based on grid-type and super-dense array plasma actuators
    Experimental investigation on de-icing by an array of impact rod-type plasma synthetic jets
    Control of flow separation over a wing model with plasma synthetic jets
    和理科男談戀愛也太“有趣”啦
    意林(2021年21期)2021-11-26 20:27:37
    《高師理科學(xué)刊》征稿簡則
    文科不懂理科的傷悲
    2017年天津卷理科第19題的多種解法
    文科十里柔情,理科一朝驚夢
    創(chuàng)業(yè)家(2015年4期)2015-02-27 07:53:09
    奇葩“租婚”生二胎,假丈夫辨白婚史十萬火急
    精品久久久久久久末码| 欧美性感艳星| 男女之事视频高清在线观看| 精品久久久久久久久久免费视频| 久久婷婷人人爽人人干人人爱| 欧美日韩瑟瑟在线播放| 午夜福利成人在线免费观看| 午夜亚洲福利在线播放| 亚洲乱码一区二区免费版| 国产亚洲精品综合一区在线观看| 很黄的视频免费| 最近视频中文字幕2019在线8| 欧美日本亚洲视频在线播放| 国产精品1区2区在线观看.| 国产精品美女特级片免费视频播放器| 大型黄色视频在线免费观看| 又爽又黄a免费视频| 国产高清有码在线观看视频| 精华霜和精华液先用哪个| 一级黄色大片毛片| 国产伦精品一区二区三区四那| 毛片女人毛片| 狠狠狠狠99中文字幕| 色吧在线观看| 成人国产综合亚洲| 有码 亚洲区| 亚洲中文字幕一区二区三区有码在线看| 在线播放国产精品三级| 狂野欧美激情性xxxx在线观看| 国内精品一区二区在线观看| 成人亚洲精品av一区二区| 最近在线观看免费完整版| 2021天堂中文幕一二区在线观| 色5月婷婷丁香| 日本撒尿小便嘘嘘汇集6| 亚洲五月天丁香| 99热只有精品国产| 22中文网久久字幕| 国产aⅴ精品一区二区三区波| 天天躁日日操中文字幕| 精品国产三级普通话版| 看免费成人av毛片| 国产亚洲精品综合一区在线观看| 大又大粗又爽又黄少妇毛片口| 色精品久久人妻99蜜桃| 国产熟女欧美一区二区| 亚洲天堂国产精品一区在线| 午夜爱爱视频在线播放| 国产主播在线观看一区二区| 91在线观看av| 免费看av在线观看网站| 日韩 亚洲 欧美在线| 午夜福利在线观看吧| 99热网站在线观看| 日韩人妻高清精品专区| 一级黄色大片毛片| 看黄色毛片网站| 俺也久久电影网| 国内揄拍国产精品人妻在线| 日本熟妇午夜| 91久久精品国产一区二区成人| 神马国产精品三级电影在线观看| 国产精品电影一区二区三区| 99在线视频只有这里精品首页| 搡老妇女老女人老熟妇| 成人三级黄色视频| 亚洲色图av天堂| 成人综合一区亚洲| 中出人妻视频一区二区| 久9热在线精品视频| 日韩精品有码人妻一区| 看免费成人av毛片| 国产激情偷乱视频一区二区| 久久午夜亚洲精品久久| 国产男人的电影天堂91| 亚洲成人久久性| 亚洲av二区三区四区| 久久久久国产精品人妻aⅴ院| 国产v大片淫在线免费观看| 久久午夜福利片| 一夜夜www| 国国产精品蜜臀av免费| 麻豆成人av在线观看| 免费观看的影片在线观看| 亚洲熟妇中文字幕五十中出| 亚洲精品日韩av片在线观看| 亚洲avbb在线观看| 成人毛片a级毛片在线播放| 99视频精品全部免费 在线| 一个人看视频在线观看www免费| 国产成人a区在线观看| 一级黄片播放器| 国产高清视频在线观看网站| av中文乱码字幕在线| 高清在线国产一区| 日韩欧美一区二区三区在线观看| 国产乱人伦免费视频| 窝窝影院91人妻| 国产精品美女特级片免费视频播放器| 亚洲av一区综合| 88av欧美| 欧美黑人欧美精品刺激| 欧美一区二区精品小视频在线| 可以在线观看的亚洲视频| 亚洲av第一区精品v没综合| 亚洲精品粉嫩美女一区| 久久国产乱子免费精品| netflix在线观看网站| 国产美女午夜福利| 亚洲va在线va天堂va国产| 色精品久久人妻99蜜桃| 欧美xxxx黑人xx丫x性爽| 午夜福利欧美成人| 日本五十路高清| 久久热精品热| 午夜免费成人在线视频| 免费无遮挡裸体视频| 国产真实乱freesex| 淫秽高清视频在线观看| 国产精品日韩av在线免费观看| 日本免费一区二区三区高清不卡| 国产亚洲精品综合一区在线观看| 欧洲精品卡2卡3卡4卡5卡区| 婷婷丁香在线五月| 琪琪午夜伦伦电影理论片6080| 草草在线视频免费看| 亚洲自拍偷在线| 国产高清视频在线播放一区| 婷婷精品国产亚洲av| 久久久久免费精品人妻一区二区| 国产综合懂色| 成人亚洲精品av一区二区| 亚洲人成伊人成综合网2020| 久久国内精品自在自线图片| 精品久久久久久久末码| 成人特级av手机在线观看| 精品久久久久久,| 欧美日本亚洲视频在线播放| 亚洲精华国产精华液的使用体验 | 色尼玛亚洲综合影院| 免费黄网站久久成人精品| 男人狂女人下面高潮的视频| 1024手机看黄色片| 女生性感内裤真人,穿戴方法视频| 又粗又爽又猛毛片免费看| 中文亚洲av片在线观看爽| av在线蜜桃| 伦理电影大哥的女人| 亚洲av成人精品一区久久| 97超级碰碰碰精品色视频在线观看| 嫁个100分男人电影在线观看| 99久久中文字幕三级久久日本| 日韩一区二区视频免费看| 亚洲人成伊人成综合网2020| 成人性生交大片免费视频hd| 国产中年淑女户外野战色| 在线免费十八禁| 性插视频无遮挡在线免费观看| 99久久无色码亚洲精品果冻| 国产伦人伦偷精品视频| 国产精品三级大全| 99热只有精品国产| 免费看a级黄色片| 国产视频一区二区在线看| 不卡视频在线观看欧美| 色综合婷婷激情| 久久久久久久久久久丰满 | 国产精品久久久久久久久免| 欧美中文日本在线观看视频| 久久天躁狠狠躁夜夜2o2o| 亚洲av成人av| 成人国产麻豆网| 国产精品久久久久久精品电影| 国产熟女欧美一区二区| 国产精品乱码一区二三区的特点| 久久国产精品人妻蜜桃| 高清在线国产一区| 欧美一级a爱片免费观看看| 久久人妻av系列| 国内久久婷婷六月综合欲色啪| 校园春色视频在线观看| 又紧又爽又黄一区二区| 老师上课跳d突然被开到最大视频| 午夜激情欧美在线| 日韩欧美精品免费久久| 一级毛片久久久久久久久女| 国产av不卡久久| 综合色av麻豆| 搡女人真爽免费视频火全软件 | 欧美3d第一页| 最新在线观看一区二区三区| 亚洲国产精品成人综合色| 极品教师在线免费播放| 国产一区二区亚洲精品在线观看| 国产91精品成人一区二区三区| av国产免费在线观看| 3wmmmm亚洲av在线观看| 日韩中字成人| 他把我摸到了高潮在线观看| 欧美高清成人免费视频www| 精品久久久噜噜| 大又大粗又爽又黄少妇毛片口| 夜夜夜夜夜久久久久| 欧美色视频一区免费| 深爱激情五月婷婷| 亚洲黑人精品在线| 久9热在线精品视频| 特大巨黑吊av在线直播| 热99re8久久精品国产| 成年人黄色毛片网站| 国产成人影院久久av| 久久国内精品自在自线图片| 成年人黄色毛片网站| 欧美性猛交黑人性爽| 精华霜和精华液先用哪个| 日本色播在线视频| 午夜a级毛片| 亚洲 国产 在线| 久久久久免费精品人妻一区二区| 免费观看精品视频网站| 尾随美女入室| 99视频精品全部免费 在线| 亚洲国产精品合色在线| 两性午夜刺激爽爽歪歪视频在线观看| 看免费成人av毛片| 久久久久久久久久成人| 国产精品一区二区免费欧美| 色av中文字幕| 一卡2卡三卡四卡精品乱码亚洲| 一区二区三区四区激情视频 | 麻豆av噜噜一区二区三区| 在线免费观看不下载黄p国产 | x7x7x7水蜜桃| 露出奶头的视频| 在线观看午夜福利视频| 色精品久久人妻99蜜桃| 69人妻影院| 舔av片在线| 亚洲av一区综合| 91在线精品国自产拍蜜月| av视频在线观看入口| 亚洲av第一区精品v没综合| 免费人成视频x8x8入口观看| 18禁裸乳无遮挡免费网站照片| 国产中年淑女户外野战色| 国产在视频线在精品| 国产午夜精品论理片| 99久久精品一区二区三区| 国产在线精品亚洲第一网站| 99久久精品一区二区三区| 亚洲经典国产精华液单| 亚洲专区中文字幕在线| 国产乱人视频| 国产视频一区二区在线看| .国产精品久久| 99国产极品粉嫩在线观看| 国产免费av片在线观看野外av| 久久精品久久久久久噜噜老黄 | 在线观看舔阴道视频| 国产亚洲精品av在线| 免费av观看视频| 不卡视频在线观看欧美| 黄色日韩在线| 不卡一级毛片| 此物有八面人人有两片| 亚洲精品日韩av片在线观看| 国产主播在线观看一区二区| 在线免费观看的www视频| a级毛片免费高清观看在线播放| 俄罗斯特黄特色一大片| 欧美黑人巨大hd| 午夜老司机福利剧场| 久久午夜福利片| 亚洲美女黄片视频| 精品久久久久久久久亚洲 | a在线观看视频网站| 麻豆国产av国片精品| xxxwww97欧美| av专区在线播放| 国语自产精品视频在线第100页| 99国产极品粉嫩在线观看| 日韩欧美一区二区三区在线观看| 深夜精品福利| 狠狠狠狠99中文字幕| 麻豆精品久久久久久蜜桃| 国产在线精品亚洲第一网站| 不卡视频在线观看欧美| 又黄又爽又免费观看的视频| 国产精品人妻久久久久久| 亚洲va日本ⅴa欧美va伊人久久| 欧美不卡视频在线免费观看| 男人舔奶头视频| 国内少妇人妻偷人精品xxx网站| 国内精品美女久久久久久| 久久热精品热| 波多野结衣高清作品| 国产精品久久久久久av不卡| 中出人妻视频一区二区| 午夜精品一区二区三区免费看| 搞女人的毛片| 午夜福利高清视频| 男女做爰动态图高潮gif福利片| 久久精品影院6| 欧美黑人巨大hd| 麻豆成人午夜福利视频| 久久精品国产亚洲网站| 亚洲欧美日韩高清在线视频| 搡老熟女国产l中国老女人| 亚洲经典国产精华液单| 深夜a级毛片| 欧美精品国产亚洲| 制服丝袜大香蕉在线| 能在线免费观看的黄片| 成人综合一区亚洲| 成人特级黄色片久久久久久久| 婷婷精品国产亚洲av在线| 欧美日韩黄片免| av在线蜜桃| www.色视频.com| 成人午夜高清在线视频| 性欧美人与动物交配| 欧洲精品卡2卡3卡4卡5卡区| 成人综合一区亚洲| 看十八女毛片水多多多| 亚洲精品亚洲一区二区| 琪琪午夜伦伦电影理论片6080| 国内精品宾馆在线| 亚洲人成网站高清观看| 亚洲最大成人中文| 亚洲人成网站在线播放欧美日韩| 国产精品福利在线免费观看| 久久久久久大精品| 51国产日韩欧美| 国产色婷婷99| 波野结衣二区三区在线| 午夜福利在线观看免费完整高清在 | 蜜桃亚洲精品一区二区三区| 女的被弄到高潮叫床怎么办 | 日本与韩国留学比较| 91麻豆精品激情在线观看国产| 成人特级av手机在线观看| 桃红色精品国产亚洲av| 日本在线视频免费播放| 国产成人一区二区在线| 亚洲一区二区三区色噜噜| 91av网一区二区| 成人特级av手机在线观看| 国产精华一区二区三区| 午夜福利成人在线免费观看| 波多野结衣巨乳人妻| 久久香蕉精品热| 亚洲精品日韩av片在线观看| 欧美日韩中文字幕国产精品一区二区三区| 日韩欧美精品免费久久| 成人二区视频| 一本精品99久久精品77| 可以在线观看毛片的网站| .国产精品久久| 午夜a级毛片| 少妇丰满av| 波多野结衣高清作品| 淫秽高清视频在线观看| 99在线视频只有这里精品首页| 成人综合一区亚洲| 麻豆久久精品国产亚洲av| 精品久久久久久久久av| 国内精品宾馆在线| 又爽又黄a免费视频| 麻豆成人av在线观看| 日韩欧美三级三区| 一本一本综合久久| 啦啦啦观看免费观看视频高清| 欧美日韩黄片免| 午夜影院日韩av| 亚洲av免费在线观看| 久久久久久久精品吃奶| x7x7x7水蜜桃| 中文字幕av在线有码专区| 亚洲最大成人中文| 欧美日韩瑟瑟在线播放| 日韩人妻高清精品专区| 我的女老师完整版在线观看| 桃红色精品国产亚洲av| 十八禁国产超污无遮挡网站| 日本黄大片高清| 欧美精品国产亚洲| 久9热在线精品视频| 国产视频内射| 成人综合一区亚洲| 亚洲经典国产精华液单| 99热网站在线观看| 一区二区三区四区激情视频 | 婷婷亚洲欧美| 国产精品久久久久久亚洲av鲁大| 亚洲精品色激情综合| 亚洲人与动物交配视频| 国产大屁股一区二区在线视频| 欧美激情国产日韩精品一区| av国产免费在线观看| 岛国在线免费视频观看| 精品99又大又爽又粗少妇毛片 | 亚洲18禁久久av| 国国产精品蜜臀av免费| 能在线免费观看的黄片| 免费在线观看日本一区| 男女啪啪激烈高潮av片| or卡值多少钱| 久久久久九九精品影院| 国产真实乱freesex| 久久久久精品国产欧美久久久| 亚洲va日本ⅴa欧美va伊人久久| 午夜免费成人在线视频| 亚洲av第一区精品v没综合| 99热这里只有精品一区| 久久精品91蜜桃| 精品久久久久久久久av| 久9热在线精品视频| 国产激情偷乱视频一区二区| 男女啪啪激烈高潮av片| 嫩草影院新地址| 又爽又黄a免费视频| 亚洲成人久久性| 日韩强制内射视频| 国产真实乱freesex| 久久国产精品人妻蜜桃| 国产午夜福利久久久久久| 日韩欧美三级三区| 亚洲欧美日韩高清专用| 麻豆一二三区av精品| 色哟哟哟哟哟哟| 日韩精品有码人妻一区| 91在线精品国自产拍蜜月| 中出人妻视频一区二区| 欧美激情久久久久久爽电影| 亚洲av.av天堂| 国模一区二区三区四区视频| 美女高潮的动态| 欧美人与善性xxx| 国产 一区精品| 亚洲av熟女| 国产久久久一区二区三区| av福利片在线观看| 日韩欧美三级三区| 久久久色成人| 国产激情偷乱视频一区二区| 日韩精品青青久久久久久| 亚洲狠狠婷婷综合久久图片| 很黄的视频免费| 亚洲性久久影院| 精华霜和精华液先用哪个| 欧美在线一区亚洲| 内地一区二区视频在线| 女同久久另类99精品国产91| 亚洲av美国av| 国产精品一区www在线观看 | 久久亚洲精品不卡| 国产精品久久视频播放| 国产精品一区二区免费欧美| 成人av一区二区三区在线看| 91久久精品国产一区二区成人| 精品久久久久久久久久久久久| 一个人观看的视频www高清免费观看| 亚洲av成人av| 久久久国产成人精品二区| av福利片在线观看| 亚洲精品国产成人久久av| 欧洲精品卡2卡3卡4卡5卡区| 国产精品一区二区性色av| 国产精品爽爽va在线观看网站| 91麻豆av在线| 亚洲欧美日韩东京热| 草草在线视频免费看| 午夜精品久久久久久毛片777| 在线国产一区二区在线| 国产高清激情床上av| 毛片女人毛片| 看黄色毛片网站| 日本色播在线视频| 97热精品久久久久久| 大又大粗又爽又黄少妇毛片口| 人人妻人人澡欧美一区二区| 97热精品久久久久久| 午夜日韩欧美国产| 丰满的人妻完整版| 在线免费十八禁| 欧美绝顶高潮抽搐喷水| 97超视频在线观看视频| 亚洲一区高清亚洲精品| 国产高清不卡午夜福利| 床上黄色一级片| 成人亚洲精品av一区二区| 12—13女人毛片做爰片一| 观看美女的网站| 午夜久久久久精精品| av天堂在线播放| 一边摸一边抽搐一进一小说| 国产激情偷乱视频一区二区| 国产蜜桃级精品一区二区三区| www.www免费av| 韩国av一区二区三区四区| 香蕉av资源在线| 亚洲成人久久爱视频| 国语自产精品视频在线第100页| 国产精品日韩av在线免费观看| 亚洲国产精品久久男人天堂| 精品人妻偷拍中文字幕| 最近最新免费中文字幕在线| 国产亚洲精品综合一区在线观看| 色噜噜av男人的天堂激情| 久久久国产成人免费| 老司机深夜福利视频在线观看| 五月玫瑰六月丁香| 婷婷丁香在线五月| 亚洲专区中文字幕在线| 天天躁日日操中文字幕| 中文字幕av在线有码专区| 亚洲成人中文字幕在线播放| 午夜福利视频1000在线观看| 亚洲,欧美,日韩| 欧美黑人欧美精品刺激| 亚洲av不卡在线观看| 久久热精品热| 亚洲va在线va天堂va国产| 亚洲欧美日韩东京热| 国产一区二区激情短视频| 色尼玛亚洲综合影院| 久久人人精品亚洲av| 日本一二三区视频观看| 国产成人一区二区在线| 长腿黑丝高跟| 亚洲图色成人| 亚洲欧美日韩无卡精品| 简卡轻食公司| 小说图片视频综合网站| 亚洲av第一区精品v没综合| 啦啦啦韩国在线观看视频| 男人舔女人下体高潮全视频| 午夜爱爱视频在线播放| 99久久九九国产精品国产免费| 国产精品日韩av在线免费观看| 亚洲天堂国产精品一区在线| 精品人妻视频免费看| 国产成人aa在线观看| 日韩中文字幕欧美一区二区| 嫩草影院精品99| 精品人妻一区二区三区麻豆 | 美女 人体艺术 gogo| 在线天堂最新版资源| 校园人妻丝袜中文字幕| 成人一区二区视频在线观看| 一a级毛片在线观看| 可以在线观看的亚洲视频| 国产精品98久久久久久宅男小说| 琪琪午夜伦伦电影理论片6080| 别揉我奶头 嗯啊视频| 国产av在哪里看| 国内精品宾馆在线| 黄色一级大片看看| 九九热线精品视视频播放| 中文字幕人妻熟人妻熟丝袜美| 国产精品福利在线免费观看| 日韩亚洲欧美综合| 网址你懂的国产日韩在线| 国产精品自产拍在线观看55亚洲| 男女那种视频在线观看| 黄片wwwwww| 99热6这里只有精品| 欧美三级亚洲精品| 亚洲国产精品成人综合色| 亚洲午夜理论影院| 国产亚洲91精品色在线| 亚洲不卡免费看| 深夜a级毛片| 午夜激情欧美在线| 国产精品久久久久久久电影| 一区二区三区四区激情视频 | 天堂动漫精品| 亚洲国产精品久久男人天堂| 久久人人精品亚洲av| 国产精品一区www在线观看 | av视频在线观看入口| 色综合婷婷激情| 国产男人的电影天堂91| 精品人妻一区二区三区麻豆 | 日韩,欧美,国产一区二区三区 | 我要看日韩黄色一级片| 亚洲在线观看片| 成人国产麻豆网| 色视频www国产| 久久精品人妻少妇| 男人和女人高潮做爰伦理| 在线观看免费视频日本深夜| 99国产精品一区二区蜜桃av| 日本熟妇午夜| 观看美女的网站| 91狼人影院| 亚洲最大成人av| 小蜜桃在线观看免费完整版高清| 欧美精品国产亚洲| 国产综合懂色| 日韩欧美一区二区三区在线观看| 婷婷丁香在线五月| 91av网一区二区| 亚洲一区高清亚洲精品| av在线观看视频网站免费| 很黄的视频免费| 最好的美女福利视频网| 国产日本99.免费观看| 在线观看av片永久免费下载| 久9热在线精品视频| 日本-黄色视频高清免费观看| 嫩草影院入口| 久久久久国产精品人妻aⅴ院| or卡值多少钱| 18+在线观看网站|