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

    Transient study of droplet oscillation characteristics driven by an electric field

    2023-12-15 11:48:20YanFeiGao高燕飛WeiFengHe何緯峰AdamAbdalazeemQiLeShi施其樂JiRongZhang張繼榮PengFeiSu蘇鵬飛SiYongYu俞思涌ZhaoHuiYao姚照輝andDongHan韓東
    Chinese Physics B 2023年12期
    關(guān)鍵詞:韓東鵬飛

    Yan-Fei Gao(高燕飛),Wei-Feng He(何緯峰),?,Adam Abdalazeem,Qi-Le Shi(施其樂), Ji-Rong Zhang(張繼榮),Peng-Fei Su(蘇鵬飛), Si-Yong Yu(俞思涌), Zhao-Hui Yao(姚照輝), and Dong Han(韓東)

    1Advanced Energy Conservation Research Group(AECRG),College of Energy and Power Engineering,Nanjing University of Aeronautics and Astronautics,Nanjing 210016,China

    2Dongfang Turbine Co.,Ltd,Deyang 618000,China

    Keywords: electrowetting,dynamic contact angle,level set model,droplet oscillation

    1.Introduction

    As widespread phenomena in various industrial fields(such as electrostatic spraying,[1]electric demulsification[2]and microfluidics[3]), droplet deformation, breakage and movement in an electric field have always been important areas of research.In particular, the wettability of droplets on solid surfaces has attracted much attention in the fields of microfluidics and heat transfer.

    The deformation, movement and wettability of droplets in an electric field are affected by many factors such as droplet location, quantity, diameter, electric field intensity,electric field frequency, fluid physical properties, surfactant and fluid type.[4]The wettability of droplets on a solid surface changes under the electric field; this is called the electrowetting phenomenon.[5]After years of development, electrowetting has been able to realize the transport of droplets between the control electrodes,[6]the directional transport of droplets in the oil phase[7]and the directional transport of underwater oil droplets.[8]Through electrowetting, droplet drive can be realized on a changing two-dimensional array electrode or surface, for example, separation, merging, mixing and transportation.Due to the advantages of in situ control, fast response, flexible manipulation and low energy consumption,electrowetting technology has been widely used in fields of microfluidics such as optics,displays,chip laboratories,printing and separation.[9]

    Some scholars have conducted many experiments and theoretical studies on electrowetting to understand how droplets interact with a solid surface under the influence of electric fields.Moreover, in recent years the study of electrowetting phenomena has begun to shift from surface characteristics to the internal mechanism.Zhanget al.[10]used a DC piezoelectric voltage to drive an electrowetting system to reduce the sliding resistance at a solid-liquid interface.The friction and slip behavior of the solid-liquid interface on a rough surface under the applied voltage was investigated.The results showed that the unstable region of rough slip behavior of droplets on the surface cannot be explained by the electrowetting equation.Heet al.[11]studied the wetting and dewetting behavior of nanodroplets with different molecular numbers on the surface of a nano-column array.The scale effect was tested by molecular dynamics (MD) simulation with or without an applied electric field.Ahmadet al.[12]reported numerical simulations and experiments on electrowetting-induced oscillatory droplets under different hydrophobic conditions,a phasefield based code with a dynamic contact angle model was used to study the droplet dynamics.Hydrophobicity of the substrate intensified the deformation and internal flow structures of the oscillating droplet.It also influenced the phase difference between the contact angle and the contact radius,producing three different phase regimes.The amplitude and frequency of the droplet oscillation change differently with the surface wettability in different phase regimes.Koet al.[13]observed hydrodynamic flow in droplets using an AC voltage from a configuration with air as the ambient phase.The flow was observed using a laser sheet and fluorescent tracers,and a toroidal vortex flow was reported.Furthermore,droplet oscillations were simultaneously observed and it was explained that the toroidal flow originates from such oscillations.Zhanget al.[14]used a DC piezoswing voltage to drive the electrowetting system to reduce the slip resistance of the solid-liquid interface.The results of experimental and theoretical analysis showed that with increase in the piezoswing rate,the slip amplitude of the contact line under electrowetting greatly increased while the corresponding slip frequency decreased.Kanget al.[15]found that under probe-type electrowetting there is vortex flow in the droplet, and this vortex flow changes direction with the frequency of applied voltage.Zhanget al.[16]conducted an MD simulation to study the static and dynamic wetting of water nanodroplets on the surface of nanostructures in the presence of a vertical electric field.Andreaset al.[17]studied the influence of parameters such as density, viscosity and mass on the electrowetting oscillation process and gave the simulation results; it turned out that only electrical grounding from below leads to utilizable oscillations.Markodimitrakiset al.[18]used experiments to prove and calculate the influence of the elasticity and thickness of the medium on the beginning of antenna saturation.The results show that the elastic effect is very important,especially when the thickness of the dielectric layer is less than 10 μm; for thicknesses greater than 20 μm the elastic effect is negligible.They attributed this finding to the effect of dielectric thickness on the electric field as well as the distribution of induced electrical stress near the threephase contact line(TCL).Malket al.[19]dealt with a characteristic hydrodynamic flow appearing in droplets actuated by electrowetting-on-dielectric(EWOD)with an AC voltage.In the coplanar electrode configuration,two pairs of vortex flows were observed to form in a droplet centered on the electrode gap.At the same time, droplet oscillations induced by AC EWOD were also revealed under stroboscopic lighting.These experiments show that vortex location can be controlled by frequency actuation with a fair degree of reproducibility.

    Based on the droplet dynamics, we can conclude that revealing the dynamic contact behavior and internal flow mechanism of the electrowetting process is essential for the performance prediction and design of advanced electrowetting devices.The effects of droplet physical characteristics,[17]electric field magnitude,[17,18]polarity[10,14]and wall hydrophobicity[12,13]on the electrowetting phenomenon under DC or AC electric fields have been extensively studied.However,the internal flow field mechanism of droplets under alternating current electrowetting(ACEW),as well as the effects of voltage frequency and initial phase angle on the dynamic contact behavior of droplets,have not been thoroughly and completely investigated.Meanwhile,the connection between physical parameters such as pressure,kinetic energy and force on droplets under ACEW and the oscillation process has not been looked at so far.To better simulate the variation of droplet contact angle,an axisymmetric model and dynamic contact angle theory are introduced in this paper.In addition, a transient numerical model of the force on the droplet is also proposed, and the variation of the droplet surface pressure is obtained.The dynamic contact behavior of droplets and the mechanism of variation of the internal flow field are first looked at from a new perspective by analyzing the force and pressure variation on droplets.Furthermore,the intrinsic connection between change in droplet kinetic energy and droplet oscillation is investigated,and the relationship between droplet velocity,energy and kinetic energy is analyzed.The proposed dynamic contact behavior and internal flow field changes of droplets under different environments provide theoretical support for the subsequent electrowetting research.

    2.Numerical models and solution methods

    2.1.Dynamic contact angle model

    The present work explores the electrowetting effect of a sessile droplet on a solid surface coated with a hydrophobic dielectric material.The basic principle of electrowetting on a medium is shown in Fig.1.In the electrowetting process the surface of the plate is usually coated with a dielectric layer to prevent electrolysis of the droplet under high voltage.In the same way a hydrophobic layer is applied to the dielectric layer to increase the hydrophobicity of the wall surface,[20]in order to increase the ability to change the contact angle:θs,OFFis the contact angle when no electric field is applied andθs,ONis the contact angle when an electric field is applied.The surface tensions between solid/liquid,liquid/gas and solid/gas are denoted byσSL,σLGandσSG,respectively.

    The surface tension between solid/liquid can be expressed by the Young equation

    whereε0,εdanddare the vacuum dielectric constant,dielectric constant and dielectric layer thickness, respectively.Formula (3) represents the contact angle of the droplet under an electric field with applied voltage(θs,ON).

    In electrowetting, an electric field force (Fe=ε0εdV2/(2d))drives the droplet along the TCL,friction resistance(Ff)and packing medium resistance(FD)as forward resistance.Considering a droplet as an unbounded sphere moving in air, its drag force is calculated asFD=0.5CDADρfU2,whereCD,AD,ρfandUdenote the drag coefficient,projected area,air density and droplet velocity,respectively.[13]According to the observation of Oprinset al.,[22]it is found that when air is the filling medium the resistance can be ignored.

    Fig.1.Schematic diagram of the principle of electrowetting on a medium.

    The imbalance of electric field force and friction force will cause acceleration on the contact line.Since the mass near the contact line is very small, the force generated by acceleration is not considered in deriving the dynamic contact angle.Thus,the resultant force(Fnet)of the oscillating droplet along the TCL direction can be written as

    Thus, the final expression for the instantaneous dynamic contact angle under electrowetting with voltageVcan be written as[14]

    Due to the singularity of the droplet at the contact line,obtaining an accurate contact linear velocityUcl(t) is a challenge.If the wall is set as a no-slip boundary, the contact line cannot move.In reality, however, the contact line does move,resulting in infinite viscous shear forces and divergent drag on the wall.To avoid this,the Navier sliding boundary condition(uslip=β?u/?y) is introduced to represent the contact line motion as a slip, whereβis the slip length, i.e., the distance from the boundary, where the velocity of linear extrapolation will reach zero.In addition, since the model is an axisymmetric regular geometry model, the contact linear velocity in Eq.(6) can be estimated as the rate of change of the droplet contact radius,which is given by[23]

    whereRis the radius of the droplet-wetting area.

    Moreover,when the friction law on the contact line is linear, the friction coefficientfcldepends on the tension interaction on the contact line, which has the same units as the viscosity and can be measured from experiments or molecular dynamics simulations.Meanwhile, Chenget al.[24]showed that the larger the slip length,the smaller the obstruction effect of the wall on the fluid.The wall resistance characteristics can be changed by adjusting the slip length.

    2.2.Numerical model

    2.2.1.Computational domain and boundary conditions

    In this paper, the numerical research is carried out in a two-dimensional axisymmetric computational domain.As shown in Fig.2(a),the droplet is assumed to be deionized water,the ambient medium is air and both the droplet and air are incompressible fluids.Layer 1 is the dielectric layer with a hydrophobic layer;layer 2 is the negative electrode and the probe is the positive electrode.Since the probe is very thin,it has little influence on the droplet flow field and is ignored in the calculation domain.After a voltage is applied, negative charges accumulate on the surface of layer 2 and positive charges accumulate on the surface of layer 1 inside the droplet; the closer to the TCL, the greater the density of positive charges accumulated on the surface of layer 2.An electric field is formed between the positive and negative charges,which disturbs the droplet.At the initial time, the droplet velocity is 0, and the air condition outside the droplet is atmospheric pressure.The droplet interface is tracked by the level set variable(φ),where the droplet has a value ofφequal to 1 and the surrounding air has a value ofφequal to 0.In addition,the numerical calculations are based on the following assumptions:

    (a)The shape of the droplet is assumed to be a spherical cap.

    (b) The voltage drop along the working fluid is ignored,thus,the droplet is assumed to be a conductor.

    (c)The effect of surfactants or any other foreign particles on deionized water is not considered.

    The mathematical model assumes that the fluid is laminar and incompressible.The laminar Navier-Stokes equations are established for droplet and air, respectively, and the effect of surface tension of the droplet is considered.[25]

    In the simulation of axisymmetric droplet diffusion,several boundary conditions are set to restore the actual scenario.S1 is the symmetric boundary,S2 is set as the wetted wall and S3 is set as the open boundary,so that the model converges by changing the system pressure release when the droplet moves.To accurately capture the influence of the stress tensor on the droplet, the mesh and wetting wall boundary near the droplet surface are specially refined,as shown in Fig.2(b).The final calculation results pass the grid independence numerical test,and the total number of grids is 14412.

    Fig.2.Geometric model for electrowetting simulation: (a) computational domains and boundary conditions;(b)finite element meshing of the model geometry.

    2.2.2.Governing equation

    The governing equations in this paper include continuity and momentum equations of level set variables.The details are as follows:

    whereg,pandFstare gravity, pressure and surface tension,respectively.

    The level set variable (φ) is used to track two phases,

    whereφ= 1 represents the secondary phase (i.e., droplet)andφ=0 represents the primary phase(i.e., air).Therefore,φ=0.5 denotes the droplet interface,where the level set variable(φ)can be calculated using the following equation:

    whereεis the interface thickness parameter andγis the reinitialization parameter, and for numerical calculations the appropriate choice isε=h,wherehis the grid cell size.

    The surface tension in Eq.(9)acting on the interface can be calculated by the normal vector of the interface and the curvature of the interface,as shown below:[23]

    Theδfunction is approximately a smoothing function

    whereσrepresents the surface tension coefficient,nrepresents the unit normal vector at the interface,krepresents the curvature of the fluid interface andδrepresents the Dirac function.

    In addition,the fluid domain is divided into two materials by the interface, and in each subdomain the material properties are constant.However, the density and viscosity are discontinuous throughout the solution domain.To simplify the calculation so that the values of density,viscosity and surface tension can continuously transition from the air to the droplet,the following formula is used:[25]

    where the subscripts g and l represent the gas and liquid,μlandμgindicate the dynamic viscosity of droplet and air, andρlandρgindicate the density of droplet and air,respectively.

    2.3.Model validation

    The friction coefficientfclcan be measured by fitting molecular dynamics simulation data to experimental data.In general, when contact angle hysteresis is not considered, it is found that the larger values (fcl= 0.2-0.4 N·s·m-2) fit the experimental data well,[26-28]otherwise a smaller value(fcl=0.18 N·s·m-2) is more suitable.[29,30]In order to verify the correctness and accuracy of the model constructed in this paper,other materials have the same characteristics under the same voltage(32 V DC)(see Table 1).In the simulation in this article,fcl=0.1 N·s·m-2andfcl=0.15 N·s·m-2are used,as shown in Fig.3.Whenβ=0.3μm andfcl=0.1 N·s·m-2,molecular dynamics simulation results data more accurately predict the experimental data of Chenget al.[24]The data obtained from the numerical simulation are basically consistent with the experimental data, and the error of the data is very small, which indicates the correctness of the model.Therefore this article choosesβ=0.3μm andfcl=0.1 N·s·m-2for subsequent simulation.

    Table 1.Material properties of liquids/substrates used in direct current electrowetting experiments and simulations.

    Fig.3.Comparison of experimental and numerical results at different times when the static contact angle is 118?.

    2.4.Force analysis model

    By analyzing the force balance of the droplet under electrowetting,the relationship between the force change and other parameters of the droplet can be obtained.Chenget al.[24]observed that viscous dissipation is not obvious in droplet oscillation dynamics with an AC electric field.In addition,Ahmadet al.[14]observed that the viscous force is negligible compared with the applied electricity and friction.Therefore, the influence of viscous force is ignored in this model.In the process of droplet oscillation,the friction forceFf=fclUclon the contact line during the oscillation of the droplet is the main part of the resistance.Therefore,Ffacts as a resistance to the electric field force in the mechanical analysis.In the process of droplet oscillation under electrowetting, droplet motion starts from the TCL.Therefore, in the force balance analysis, the dominant force along the TCL ignores other forces in the vertical direction.

    The applied sinusoidal waveform voltage(V)is given by

    whereK=ε0εd/(2d).

    3.Results and analysis

    In the simulation calculation,the droplet volume is 5μl.The initial shape of the droplet is a spherical crown, and the initial contact angle is 110?.The total thickness of the dielectric and hydrophobic layerd=2.5 μm and the total effective relative dielectric constantεd= 2.8, The droplet is deionized water,withρ=997 kg·m-3and kinematic viscosityμ=1.01×10-3Pa·s.The ambient gas is air,and the droplet is in an equilibrium state before the voltage is applied(t=0 ms).

    3.1.Dynamic contact characteristics of droplets under ACEW

    Figure 4 shows the variation of droplet contact radius with time,whereRwis the wetting radius of the droplet.The amplitude of the droplet is the displacement of the droplet contact line moving forward along the wall during diffusion,and is represented by the symbolAm.After completion of the first oscillation,the droplet maintains a new contact radius higher than the initial contact radius and continues to oscillate with the maximum possible diffusion; then, the time it takes to return to the original contact radius during a certain period of oscillation is called the stable oscillation periodT.As can be seen from Fig.4, the variation of the displacementRwof the droplet contact line tends to be a sinusoidal waveform.The AC amplitude is kept constant,and its magnitude is 75 V.When the AC frequency increases from 50 Hz to 500 Hz,the oscillating wave amplitudes after the droplet stabilization are 0.036 mm,0.016 mm,0.013 mm and 0.002 mm,respectively,and the oscillation periodsTof the droplet wetting radius are 11 ms, 4 ms, 2 ms and 1 ms, respectively.It can be seen that with the increase in AC frequency, the stable oscillation period of the droplet decreases stepwise,and the amplitude of the stable waveform after oscillation also decreases stepwise.It can be seen from the analysis that the surface of the droplet under stimulation with different frequencies presents different surface waves,which causes the movement of the contact line.In fact, there is a delay between the movement of the droplet contact line and the change in the electric field force.The higher the AC frequency is, the faster the rate of change the electric field force will be: if the electric field force changes too fast, it will hinder the process of droplet contact wire movement, which shows that the larger the AC frequency is,the smaller the oscillation amplitude of the droplet is.It can be seen from Fig.4(a)that whenf=500 Hz,the droplet oscillation amplitude is very small, which is similar to the variation trend of droplet wetting radius under direct current electrowetting(the DC voltage is the effective value corresponding to the AC amplitude).

    Fig.4.Dynamics of droplets under ACEW with time.(a)Variation of droplet wetting radius with time at different frequencies.(b)Oscillation period of droplet wetting radius at different frequencies.

    Figure 5 shows the time variation of different forces(electric field forceFe,contact line frictionFfand net forceFnet)on the droplet atV0=75 V,f=50 Hz and?=0?.As shown in Fig.5, the line of contact frictionFf=fclUcl, the line of contact friction force and the direction depend on the value of the linear velocity.The speed changes along the sine wave,so the line of contact friction also presents a sinusoidal waveform change, with the positive and negative friction corresponding to expanding and contracting behavior of the contact wires.As shown in Fig.4(a), withint=0-6 ms, the droplet amplitude rapidly increases to a peak,and aftert=6 ms,the droplet amplitude presents a waveform oscillation, corresponding to Fig.5.The droplet velocity reaches its maximum at the initial time, the contact line friction reaches its maximum and then the droplet size presents waveform fluctuation.As Ahmadet al.[12]mentioned, the magnitude and direction of friction depend on the speed of the contact line.The positive and negative values of friction represent the diffusion stage and the contraction stage, respectively.Friction increases with time,and reaches a maximum when the electric field force reaches a peak.At the same time,it is accompanied by the maximum change in diffusion radius.In addition,the staticFnetsize near the droplet contact line in Fig.5 also changes with a waveform oscillation, which further explains that the variation curve of stable droplet oscillation amplitude under the action of electrowetting is a sinusoidal waveform.

    Fig.5.Change of force on the contact line with time (V0 = 75 V,f =50 Hz,? =0?).

    Figure 6 shows the change in droplet contact behavior over time when the AC voltage is the same and the initial phase angle is different.As can be seen from Fig.6, when the initial phase angle is 0?, 45?or 90?there is no significant difference in the droplet wetting radius,indicating that the initial phase angle does not affect the change in the droplet wetting radius.However,the initial phase angle is related to the initial phase angle of the droplet oscillation curve.For?=45?, 0?and 90?,the droplet reaches stable oscillation successively and then continues to move with the same stable oscillation period.

    Fig.6.Change of droplet contact radius under different phase angles.

    3.2.Droplet internal flow under ACEW

    Fig.7.Change of average pressure at the gas-liquid interface with time.

    Fig.8.The change of droplet morphology with different frequency.(a)Internal flow field variation.(b)Vortex and wave crest relationship diagram.

    Under the action of an electric field force, the contact angle of stationary droplets changes, as shown in Fig.4(a).The droplet presents a sinusoidal oscillation motion along the wall.It can be seen that under the effect of electrowetting,the static internal flow of the droplet is changed into an internal eddy current under the stimulation of an electric field.Figure 7 shows the changes in the average pressure at the gas-liquid surface with time at different frequencies.It can be seen from Fig.7 that under AC stimulation, the pressure at the interface rapidly increases from close to atmospheric pressure to the maximum pressure value, and then the interface pressure changes to show waveforms.The reason for the pressure change can be obtained by referring to the change in the droplet contact radius over time in Fig.4(a).In the initial deformation stage of the droplet, the force balance of the droplet is broken under the stimulation of the electric field force,and a certain kinetic energy is accumulated,which eventually leads to the accumulation of pressure at the edge of the droplet.Under the action of pressure,the droplet reciprocates,resulting in the accelerated retreat of the contact line and oscillation of the wetting radius.At the same time,since the voltage applied to the droplet constantly changes with time,the pressure near the droplet contact line will also constantly change with the trend of voltage variation.In addition,the larger the frequencyfis, the smaller the accumulated pressure oscillation amplitude at the interface will be.The reason for this is that the pressure change has a delay compared with the voltage change.The larger the voltage frequency is,the faster the voltage change will be, so the change of electric field force will also be accelerated.Therefore,an electric field force that changes too quickly will interfere with the pressure accumulation process at the interface.The larger the frequencyfis,the stronger the interference will be,and the smaller the eventual change of the interface pressure will be.

    Figure 8 shows the change of droplet morphology over time with different frequencies.Figure 8(a)shows the velocity vector field and pressure distribution inside the droplet at different frequencies.In this figure,half an oscillation period is selected for modesP2andP4and one oscillation period is selected for modesP8.In the process of droplet oscillation,each peak value of droplet oscillation corresponds to a resonance modePn(n=2,4,...), and the value ofndepends on the number of wave peaks at the droplet interface.As can be seen from Fig.4, the higher the resonance mode is, the smaller the droplet oscillation amplitude is.In addition, the recirculating eddies can be identified from the streamline profile and the number of eddies matches their resonance modes.For example, two vortices can be observed inP2mode, four inP4mode and eight inP8mode,and most of the vortex centers are located on or near the free surface.It is also noted that the oscillation amplitude of the droplet decreases continuously as the frequencyfincreases from 50 Hz to 250 Hz.As can be seen from Fig.8(a), the oscillation process of the droplet decreases continuously as the oscillation mode of the droplet changes from modeP2to modeP8, and the size of the vortex at the interface also decreases continuously.As mentioned by Yiet al.,[23]when a droplet is actuated under ACEW,the droplet exhibits time-harmonic shape oscillations.The oscillation amplitude is found to decrease with increasing actuation frequency due to the dominant effect of inertia over other forces.The surface wave propagates with a similar phase velocity to capillary waves.

    The connection between the droplet surface pressure and the resonance mode can be investigated through Figs.7 and 8.Whenf=50 Hz,the pressure changes slowly,so the droplet surface forces change more slowly and can be seen as a more stable oscillation.Due to the slow change of droplet surface perturbation, there are fewer wave peaks and vortices on the surface.At the same time,whenf=250 Hz the rate of pressure change is accelerated,so the force balance on the droplet surface changes dramatically.Due to the strong change in perturbation on the droplet surface,the wave and vortex currents on the surface also increase.In addition, the pressure amplitude becomes smaller fromf=50 Hz to 250 Hz and the intensity of droplet surface perturbation decreases,resulting in a reduction in the size of wave crests and vortices.

    Figure 8(b) shows the quantitative correlation between the number of vorticesnon the droplet surface and the number of wave peakswat the interface during the oscillation.Resonance modes areP2corresponding to one wave,P4to two waves,P8to four waves andP12to six waves.It can be concluded from Fig.8(b) thatw=z/2, namely, the number of vortices on the droplet surfacezis always twice the number of wave peaks, and the interface wave peaks caused by droplet oscillation are always symmetric along the center line.The largerfis, the more frequently the interface is disturbed, but the intensity of the disturbance of the electric field force is decreasing.Therefore, under voltage stimulation, the larger frequency of the droplet the more surface crests and vortices there are and the smaller the size of ripples generated by the droplet surface is.

    4.Conclusions

    Based on our numerical simulation, the dynamic contact process and internal flow field of droplets under electrowetting are studied in this paper.The following conclusions are drawn.

    1.The dynamic contact angle model is used to simulate and predict the electrowetting process more reasonably, and the change in droplet contact radius under different frequencies and initial phase angles is studied.When the AC frequency increases from 50 Hz to 500 Hz, the oscillating wave amplitudes of droplet stabilization are 0.036 mm, 0.016 mm,0.013 mm and 0.002 mm,respectively,and the oscillation periodsTof the droplet wetting radius are 11 ms, 4 ms, 2 ms and 1 ms, respectively.With increase in frequency, the stable oscillation period of the droplet decreases stepwise, and the amplitude of the stable waveform after oscillation also decreases stepwise.In addition, when the initial phase angle is 0?,45?or 90?,there is no significant difference in the droplet wetting radius,indicating that the initial phase angle does not affect the change in the droplet wetting radius.

    2.A model of droplet force variation in the electrowetting process is proposed, and it is concluded that the friction of the contact line presents sinusoidal waveform changes.The positive and negative friction forces, respectively, correspond to the expansion and contraction behaviors of the contact line.The staticFnetnear the droplet contact line also changes as a waveform oscillation, which further explains why the fluctuation curve of stable droplet oscillation amplitude under the action of electrowetting is a sinusoidal waveform.

    3.In the case of ACEW, the droplet interface will show velocity vortices and oscillating ripples corresponding to the frequency.As the frequencyfincreases from 50 Hz to 250 Hz the oscillation amplitude of the droplet decreases continuously, and the oscillation mode of the droplet changes fromP2toP8mode.The oscillation process of the droplet decreases continuously,and the vortex size at the interface also decreases continuously.

    Acknowledgements

    The authors expressed their sincere gratitude to the Natural Science Foundation of Jiangsu Province (Grant No.BK2020194), the Basic Research Fund of Central University (Grant No.NS2022026), and the Graduate Research and Practice Innovation Program(Grant No.xcxjh20220215).

    猜你喜歡
    韓東鵬飛
    一匹馬
    詩林(2022年2期)2022-04-13 08:38:49
    可愛的敘事者——讀韓東短篇小說《素素與李蕓》
    都市(2022年1期)2022-03-08 02:23:34
    Quality Control for Traditional Medicines - Chinese Crude Drugs
    為了避嫌
    雜文月刊(2019年18期)2019-12-04 08:30:40
    懲“前”毖“后”
    21世紀(jì)(2019年10期)2019-11-02 03:17:02
    執(zhí)“迷”不悟
    21世紀(jì)(2019年10期)2019-11-02 03:17:02
    舉賢
    21世紀(jì)(2019年9期)2019-10-12 06:33:44
    重新做人
    詩潮(2017年4期)2017-12-05 10:16:18
    落后的家鄉(xiāng)
    上海故事(2017年4期)2017-04-18 16:30:59
    迷信的老爸
    三月三(2016年7期)2016-08-23 10:18:59
    婷婷精品国产亚洲av在线| 校园春色视频在线观看| 国内少妇人妻偷人精品xxx网站| 老司机在亚洲福利影院| 国产aⅴ精品一区二区三区波| 午夜老司机福利剧场| 久久久久久久久久黄片| 神马国产精品三级电影在线观看| 国产精品乱码一区二三区的特点| 级片在线观看| 男人的好看免费观看在线视频| 露出奶头的视频| 欧美在线黄色| 每晚都被弄得嗷嗷叫到高潮| 美女被艹到高潮喷水动态| 很黄的视频免费| 中亚洲国语对白在线视频| 99热只有精品国产| 午夜福利在线在线| 波野结衣二区三区在线 | 国产单亲对白刺激| 国模一区二区三区四区视频| 国产成人福利小说| 亚洲精品日韩av片在线观看 | 国产一区在线观看成人免费| 999久久久精品免费观看国产| 香蕉久久夜色| 亚洲天堂国产精品一区在线| 日韩有码中文字幕| 69人妻影院| www.999成人在线观看| 国产高清视频在线播放一区| 搡老熟女国产l中国老女人| 在线播放国产精品三级| 热99re8久久精品国产| 久久久久久久精品吃奶| 五月玫瑰六月丁香| 国产午夜精品久久久久久一区二区三区 | 午夜视频国产福利| 性欧美人与动物交配| 国产亚洲精品久久久com| 性色av乱码一区二区三区2| 亚洲av免费在线观看| 欧美又色又爽又黄视频| 精品不卡国产一区二区三区| 精品久久久久久久人妻蜜臀av| 午夜福利成人在线免费观看| 一个人免费在线观看电影| 波多野结衣高清作品| 制服人妻中文乱码| 亚洲人与动物交配视频| 中亚洲国语对白在线视频| 最近视频中文字幕2019在线8| 欧美bdsm另类| 91字幕亚洲| 久久香蕉精品热| 亚洲一区二区三区不卡视频| 一本精品99久久精品77| 国产成人欧美在线观看| 国产成人av教育| 午夜精品在线福利| 可以在线观看的亚洲视频| 中文字幕精品亚洲无线码一区| 成年女人毛片免费观看观看9| 国内久久婷婷六月综合欲色啪| 两个人看的免费小视频| 国产午夜精品论理片| 亚洲色图av天堂| 成人国产一区最新在线观看| 老熟妇仑乱视频hdxx| 成人三级黄色视频| 99精品欧美一区二区三区四区| 亚洲av一区综合| 午夜精品在线福利| 首页视频小说图片口味搜索| 内射极品少妇av片p| 在线视频色国产色| 久久久久久久精品吃奶| 亚洲天堂国产精品一区在线| 大型黄色视频在线免费观看| 51午夜福利影视在线观看| 欧美一级毛片孕妇| 一区二区三区激情视频| 亚洲熟妇中文字幕五十中出| 午夜激情福利司机影院| 99国产精品一区二区三区| 午夜福利欧美成人| 国产亚洲精品久久久com| 高潮久久久久久久久久久不卡| 人妻夜夜爽99麻豆av| 长腿黑丝高跟| www国产在线视频色| 在线十欧美十亚洲十日本专区| 亚洲国产精品999在线| 亚洲av第一区精品v没综合| 欧美另类亚洲清纯唯美| 麻豆成人av在线观看| 国产精品免费一区二区三区在线| 国产亚洲av嫩草精品影院| 观看免费一级毛片| 成人国产一区最新在线观看| 亚洲aⅴ乱码一区二区在线播放| 亚洲avbb在线观看| 亚洲av二区三区四区| 免费av毛片视频| 禁无遮挡网站| 一级毛片高清免费大全| 搡老岳熟女国产| av视频在线观看入口| 一卡2卡三卡四卡精品乱码亚洲| 老熟妇乱子伦视频在线观看| 高清日韩中文字幕在线| 亚洲人与动物交配视频| 免费无遮挡裸体视频| 午夜免费观看网址| 国产亚洲av嫩草精品影院| 国产精品电影一区二区三区| 在线视频色国产色| 国产成人影院久久av| 看免费av毛片| 两个人的视频大全免费| 欧美在线黄色| tocl精华| 少妇的丰满在线观看| 国产亚洲精品一区二区www| 欧美在线黄色| 1024手机看黄色片| 少妇熟女aⅴ在线视频| 精品国产三级普通话版| 别揉我奶头~嗯~啊~动态视频| 午夜日韩欧美国产| 午夜福利欧美成人| 午夜精品久久久久久毛片777| 18+在线观看网站| 美女大奶头视频| 欧美成人a在线观看| 国产又黄又爽又无遮挡在线| 女人高潮潮喷娇喘18禁视频| 午夜免费激情av| 啦啦啦观看免费观看视频高清| 国产爱豆传媒在线观看| 最近视频中文字幕2019在线8| 亚洲av不卡在线观看| 国产视频一区二区在线看| 亚洲精品色激情综合| 在线观看免费午夜福利视频| 欧美日韩精品网址| 深夜精品福利| 亚洲人成伊人成综合网2020| 婷婷精品国产亚洲av| 国产精品,欧美在线| 亚洲成人久久爱视频| 亚洲无线观看免费| 国产精品,欧美在线| bbb黄色大片| 成人永久免费在线观看视频| 欧美黑人欧美精品刺激| 国产69精品久久久久777片| 一个人看的www免费观看视频| 国产精品,欧美在线| 黄色片一级片一级黄色片| 欧美+日韩+精品| 国产精品久久久人人做人人爽| tocl精华| 久久精品亚洲精品国产色婷小说| 亚洲激情在线av| 午夜精品久久久久久毛片777| 国产伦一二天堂av在线观看| 精品乱码久久久久久99久播| 久久久国产成人精品二区| 免费av不卡在线播放| 亚洲18禁久久av| 美女 人体艺术 gogo| 日本熟妇午夜| 国产av不卡久久| www国产在线视频色| 成人特级黄色片久久久久久久| 国产三级黄色录像| 国产高清videossex| 大型黄色视频在线免费观看| 一边摸一边抽搐一进一小说| 一本综合久久免费| 久久久久久久久大av| 午夜激情欧美在线| 精品不卡国产一区二区三区| 深爱激情五月婷婷| 人妻丰满熟妇av一区二区三区| 亚洲乱码一区二区免费版| 成人无遮挡网站| 日日摸夜夜添夜夜添小说| www日本黄色视频网| 国产亚洲精品久久久com| 国产午夜精品论理片| 91麻豆精品激情在线观看国产| 在线观看免费午夜福利视频| 老熟妇乱子伦视频在线观看| 最近视频中文字幕2019在线8| 欧美国产日韩亚洲一区| 日本与韩国留学比较| 亚洲精华国产精华精| 亚洲av中文字字幕乱码综合| 久久精品夜夜夜夜夜久久蜜豆| 亚洲精品影视一区二区三区av| 国产色爽女视频免费观看| 精品人妻1区二区| 嫩草影院入口| 少妇的逼水好多| 草草在线视频免费看| 18禁黄网站禁片午夜丰满| 少妇熟女aⅴ在线视频| 国产精品一区二区三区四区久久| 欧美黄色淫秽网站| 国产精品综合久久久久久久免费| 午夜老司机福利剧场| 久久精品国产综合久久久| 亚洲激情在线av| av片东京热男人的天堂| 两个人的视频大全免费| 久久久久国内视频| 中文资源天堂在线| 国产一区二区在线av高清观看| 狂野欧美激情性xxxx| 久久精品91蜜桃| 制服丝袜大香蕉在线| 观看免费一级毛片| 午夜a级毛片| 99精品久久久久人妻精品| 国产精品国产高清国产av| 十八禁网站免费在线| 狂野欧美激情性xxxx| 悠悠久久av| 亚洲av熟女| 99热精品在线国产| 国产精品一及| 老司机福利观看| 91字幕亚洲| 成人无遮挡网站| 亚洲国产日韩欧美精品在线观看 | 午夜精品一区二区三区免费看| 一级a爱片免费观看的视频| 一本一本综合久久| 亚洲av五月六月丁香网| 免费看美女性在线毛片视频| 精品一区二区三区av网在线观看| 国内揄拍国产精品人妻在线| 成年免费大片在线观看| 成人av一区二区三区在线看| 在线观看午夜福利视频| 激情在线观看视频在线高清| 免费看美女性在线毛片视频| 看免费av毛片| 美女高潮喷水抽搐中文字幕| 99久久成人亚洲精品观看| 日韩成人在线观看一区二区三区| 国产美女午夜福利| 麻豆久久精品国产亚洲av| 日日干狠狠操夜夜爽| 成人精品一区二区免费| 欧美午夜高清在线| 国产精品久久电影中文字幕| 丁香欧美五月| e午夜精品久久久久久久| 身体一侧抽搐| www.www免费av| 伊人久久精品亚洲午夜| 国产精品女同一区二区软件 | 亚洲精品乱码久久久v下载方式 | 精品一区二区三区av网在线观看| 国产精品香港三级国产av潘金莲| 国产69精品久久久久777片| 99久久久亚洲精品蜜臀av| 免费人成视频x8x8入口观看| 无限看片的www在线观看| 亚洲内射少妇av| 51国产日韩欧美| 免费av观看视频| 国产精品久久久人人做人人爽| 嫩草影院精品99| 久久久久九九精品影院| 亚洲国产欧洲综合997久久,| 午夜福利在线观看吧| 一卡2卡三卡四卡精品乱码亚洲| 三级国产精品欧美在线观看| 国产精品一及| av女优亚洲男人天堂| 成人高潮视频无遮挡免费网站| 亚洲美女视频黄频| 久久天躁狠狠躁夜夜2o2o| 欧美一级毛片孕妇| 午夜两性在线视频| 麻豆国产97在线/欧美| 狂野欧美白嫩少妇大欣赏| 麻豆一二三区av精品| 亚洲国产精品sss在线观看| 色精品久久人妻99蜜桃| 中文在线观看免费www的网站| 亚洲 欧美 日韩 在线 免费| 哪里可以看免费的av片| 在线免费观看不下载黄p国产 | 69av精品久久久久久| 亚洲天堂国产精品一区在线| 国产一区二区三区视频了| 色综合站精品国产| 久久久国产成人精品二区| 变态另类成人亚洲欧美熟女| 日本精品一区二区三区蜜桃| 好男人在线观看高清免费视频| netflix在线观看网站| 丁香欧美五月| 欧美成人a在线观看| 亚洲在线观看片| 免费在线观看影片大全网站| 丰满乱子伦码专区| www国产在线视频色| av天堂中文字幕网| 一二三四社区在线视频社区8| 亚洲五月天丁香| 成年免费大片在线观看| 亚洲国产精品合色在线| 99热只有精品国产| 成人性生交大片免费视频hd| 亚洲无线观看免费| 草草在线视频免费看| 欧美乱妇无乱码| 9191精品国产免费久久| 精华霜和精华液先用哪个| 性色avwww在线观看| av专区在线播放| 亚洲精品在线观看二区| 狠狠狠狠99中文字幕| 法律面前人人平等表现在哪些方面| 国产毛片a区久久久久| 国产免费av片在线观看野外av| 午夜福利在线观看吧| 制服人妻中文乱码| 嫩草影视91久久| 91在线精品国自产拍蜜月 | 中文字幕av成人在线电影| 亚洲性夜色夜夜综合| 国产麻豆成人av免费视频| eeuss影院久久| 久久久久久国产a免费观看| 可以在线观看的亚洲视频| 欧美日本视频| 国产蜜桃级精品一区二区三区| 一本一本综合久久| 亚洲成人中文字幕在线播放| 69av精品久久久久久| av黄色大香蕉| 亚洲一区二区三区色噜噜| 18禁在线播放成人免费| 婷婷丁香在线五月| 欧美三级亚洲精品| 色综合站精品国产| 首页视频小说图片口味搜索| 亚洲av五月六月丁香网| 国产国拍精品亚洲av在线观看 | 夜夜躁狠狠躁天天躁| 日本免费一区二区三区高清不卡| av在线天堂中文字幕| 狂野欧美激情性xxxx| 日韩精品中文字幕看吧| 97超视频在线观看视频| 精品无人区乱码1区二区| 国产免费男女视频| 他把我摸到了高潮在线观看| 在线观看日韩欧美| 国产午夜精品论理片| 蜜桃亚洲精品一区二区三区| 在线a可以看的网站| 国产精品一及| 国产成人av教育| 国产午夜精品论理片| 久久久久久久久大av| 在线观看免费午夜福利视频| 天天一区二区日本电影三级| 一a级毛片在线观看| 亚洲成av人片免费观看| 在线播放无遮挡| av欧美777| 国产精品美女特级片免费视频播放器| 国产精品乱码一区二三区的特点| 午夜影院日韩av| 乱人视频在线观看| 国产精品久久久久久久久免 | 国产视频内射| 脱女人内裤的视频| 好男人电影高清在线观看| 亚洲精品美女久久久久99蜜臀| 日韩欧美精品免费久久 | 美女高潮的动态| 欧美日韩国产亚洲二区| 国产精品久久久久久亚洲av鲁大| 欧美日韩中文字幕国产精品一区二区三区| av片东京热男人的天堂| 中文字幕av在线有码专区| 又粗又爽又猛毛片免费看| av在线天堂中文字幕| 人妻久久中文字幕网| 又黄又爽又免费观看的视频| 97人妻精品一区二区三区麻豆| 久久亚洲真实| 国产高清videossex| 日韩中文字幕欧美一区二区| 伊人久久大香线蕉亚洲五| 欧美黄色片欧美黄色片| 欧美成人一区二区免费高清观看| 色哟哟哟哟哟哟| 草草在线视频免费看| 男插女下体视频免费在线播放| 国产亚洲精品综合一区在线观看| 国产欧美日韩一区二区三| 国产亚洲精品久久久com| 少妇熟女aⅴ在线视频| 男女床上黄色一级片免费看| 男人和女人高潮做爰伦理| 成人欧美大片| 搡老岳熟女国产| 欧美成人一区二区免费高清观看| 观看免费一级毛片| 日韩精品中文字幕看吧| 他把我摸到了高潮在线观看| 在线观看免费午夜福利视频| 少妇的逼水好多| 精华霜和精华液先用哪个| 亚洲最大成人手机在线| 成人精品一区二区免费| 国产成人系列免费观看| 毛片女人毛片| 日日夜夜操网爽| 日韩精品青青久久久久久| 色综合站精品国产| 成年女人毛片免费观看观看9| 亚洲精品一卡2卡三卡4卡5卡| 一二三四社区在线视频社区8| 亚洲av不卡在线观看| 两个人的视频大全免费| 日韩欧美免费精品| 久久久久久久亚洲中文字幕 | 高清毛片免费观看视频网站| 少妇的丰满在线观看| 亚洲电影在线观看av| 午夜亚洲福利在线播放| 国产色婷婷99| 国产真人三级小视频在线观看| 午夜精品久久久久久毛片777| 日韩精品中文字幕看吧| 淫秽高清视频在线观看| 国产精品永久免费网站| 国产精品久久久久久亚洲av鲁大| 在线观看午夜福利视频| 韩国av一区二区三区四区| 国产成年人精品一区二区| 国产精品久久电影中文字幕| 97碰自拍视频| 国产亚洲av嫩草精品影院| 久久久久国产精品人妻aⅴ院| 久久久久久国产a免费观看| 特大巨黑吊av在线直播| 精品人妻偷拍中文字幕| 免费人成视频x8x8入口观看| 51国产日韩欧美| 国产熟女xx| 亚洲精品美女久久久久99蜜臀| 色av中文字幕| 国产av在哪里看| 麻豆成人午夜福利视频| 久久精品亚洲精品国产色婷小说| 国产av不卡久久| 欧美大码av| 丁香欧美五月| 午夜激情福利司机影院| 婷婷精品国产亚洲av| 欧美一区二区亚洲| 国产精品野战在线观看| 黄片小视频在线播放| 丰满人妻熟妇乱又伦精品不卡| 亚洲欧美日韩无卡精品| 可以在线观看的亚洲视频| 麻豆成人av在线观看| 欧美日韩乱码在线| 淫妇啪啪啪对白视频| eeuss影院久久| e午夜精品久久久久久久| 人妻夜夜爽99麻豆av| 色精品久久人妻99蜜桃| 久久久久免费精品人妻一区二区| 亚洲内射少妇av| 国产在视频线在精品| av国产免费在线观看| 欧美色欧美亚洲另类二区| 三级国产精品欧美在线观看| 久久精品91无色码中文字幕| 脱女人内裤的视频| 2021天堂中文幕一二区在线观| 制服丝袜大香蕉在线| 熟女少妇亚洲综合色aaa.| 九九在线视频观看精品| 色综合亚洲欧美另类图片| 两个人视频免费观看高清| x7x7x7水蜜桃| 国产精品一区二区免费欧美| 亚洲国产高清在线一区二区三| 亚洲精品色激情综合| 人妻久久中文字幕网| 午夜福利欧美成人| 性欧美人与动物交配| 国产高潮美女av| 国产精品精品国产色婷婷| 99热精品在线国产| eeuss影院久久| 最近在线观看免费完整版| 日本熟妇午夜| 欧美性猛交黑人性爽| 黄色丝袜av网址大全| 搡老岳熟女国产| xxxwww97欧美| 欧美区成人在线视频| 日本免费a在线| 美女被艹到高潮喷水动态| 特级一级黄色大片| 亚洲av美国av| 高清在线国产一区| 美女免费视频网站| 午夜精品在线福利| 女同久久另类99精品国产91| av欧美777| 91麻豆av在线| 日本黄大片高清| 免费搜索国产男女视频| 中文资源天堂在线| 国产国拍精品亚洲av在线观看 | 小蜜桃在线观看免费完整版高清| 国产精品久久久久久精品电影| 热99re8久久精品国产| 不卡一级毛片| 一本久久中文字幕| 亚洲av二区三区四区| 国产精品女同一区二区软件 | 美女黄网站色视频| 亚洲精品乱码久久久v下载方式 | 久久久久久久亚洲中文字幕 | 91在线精品国自产拍蜜月 | 久久香蕉国产精品| 观看免费一级毛片| 亚洲片人在线观看| 久久久国产成人精品二区| 久久久久久人人人人人| 男人的好看免费观看在线视频| 一个人看的www免费观看视频| 日日干狠狠操夜夜爽| 国产成+人综合+亚洲专区| 美女高潮喷水抽搐中文字幕| 尤物成人国产欧美一区二区三区| av欧美777| 老司机午夜福利在线观看视频| 老熟妇仑乱视频hdxx| 免费观看的影片在线观看| 内射极品少妇av片p| 亚洲av日韩精品久久久久久密| 99热这里只有是精品50| 精华霜和精华液先用哪个| av片东京热男人的天堂| 亚洲专区中文字幕在线| 欧美成狂野欧美在线观看| 女同久久另类99精品国产91| 欧美激情在线99| 国内精品一区二区在线观看| 51午夜福利影视在线观看| 久久久久免费精品人妻一区二区| 淫秽高清视频在线观看| 精品久久久久久久人妻蜜臀av| 亚洲精品粉嫩美女一区| 窝窝影院91人妻| 亚洲精品粉嫩美女一区| 日日夜夜操网爽| 变态另类丝袜制服| 99热精品在线国产| 九色成人免费人妻av| 日韩大尺度精品在线看网址| 亚洲精品色激情综合| 国产成人av激情在线播放| 天堂动漫精品| 亚洲欧美日韩高清在线视频| 亚洲av免费在线观看| 亚洲在线自拍视频| 亚洲精品美女久久久久99蜜臀| 国产成人av激情在线播放| www.999成人在线观看| 一级毛片女人18水好多| 国产成人福利小说| 精品久久久久久久人妻蜜臀av| АⅤ资源中文在线天堂| 高潮久久久久久久久久久不卡| 中文字幕人妻丝袜一区二区| 丰满人妻熟妇乱又伦精品不卡| 精品国产三级普通话版| 国产久久久一区二区三区| 搡女人真爽免费视频火全软件 | 18美女黄网站色大片免费观看| 99视频精品全部免费 在线| 久久天躁狠狠躁夜夜2o2o| 国产精品亚洲一级av第二区| 黄片大片在线免费观看| 999久久久精品免费观看国产| 麻豆国产av国片精品| 国产乱人视频| 在线观看日韩欧美| 97碰自拍视频| 两人在一起打扑克的视频| 美女免费视频网站| 国产久久久一区二区三区| 久久精品国产清高在天天线| 淫妇啪啪啪对白视频| 国产一区二区激情短视频|