更全的杂志信息网

Study of droplet flow in a T-shape microchannel with bottom wall fluctuation

更新时间:2016-07-05

1 Introduction

Droplet microfluidics, in which nanoliter to picoliter droplets act as individual compartments,has enabled a wide range of applications.These types of systems can save the reagent and improve the mixing or testing effect[1],prevent important fluids from pollution or chemical reaction with other materials[2],fabricate nanoparticles for photochemical synthesis due to the high surface-to-volume-ratios[3]and so on.As a result,there is a growing interest in the microfluidic techniques,especially for the formation and control of droplets or bubbles.

作者罗琳大量运用夸张手法,深刻、生动地揭示事物的本质,增强语言的感染力,给人以深刻的印象。从而使读者对所表达事物产生鲜明的印象,产生强烈的感情,引起共鸣。同时结合读者丰富的想象力,更好地突出事物的特征,从人物内心深处感悟角色和故事情节。

Droplet size and velocity are important input parameters in the microfluidic systems[4,5].Generally,the droplet size is primarily determined by the dimensions of the fluid flowing part of the chip[6].In addition, flow parameters,such as the flow rate[7,8],viscosity[5],interfacial tension[9],surface wettability[10],surface roughness[9],affect the droplet generation process.The work of Garstecki et al.[7]firstly documented the relation of the droplet length(L)and width of the main channel(w)in a T-shape microchannel as L/w=1+ αQ d/Q c,where α is a constant depending on the geometry of the channel,Q d is the volumetric flow rate of the dispersed phase,and Q c is the volumetric flow rate of the continuous phase.After that,the droplet size is scaled in some other channels,such as the flow-focusing channel[11],and many modified and specific scaling laws[12–14]are reported to describe the droplet size basing on the emerging process. Many important properties of droplets,including stability,chemical reactivity,and physical characteristics are highly dependent on the parameters of the physical size and the size distribution of droplets[15,16].

Since controlling the size of droplets to raise the specific performance is in great need,there is growing interests in such investigations.In order to generate dynamically droplets with tunable size,the tunable membrane is adopted and controlled by the pressure of the chamber beside them[15,17,18].The membrane changes the width of the droplet emerging channel and finally alter the droplet size.The study by Li et al.[19]involved the perturbation into the driving of the dispersed phase in the flow-focusing channel.They found that the higher perturbation frequency reduced the droplet size,and the response to the change in flow conditions toward the hydrodynamic perturbation was significantly faster than that triggered by changes in the flow rate using typical syringe pumps.Then Huang et al.[20]also adopted this type of perturbation in the high throughput size-tunable generation system of cell containing collagen microbeads.For tunable and unequal droplets,Link et al.[21]proposed two methods,using an obstacle inside a straight channel whose drawback was that the breakup droplets flow together and an additional carefully designed separating system was needed,and a T-junction with unequal length narrow branches whose pressure drop could be increased.Similarly,Bedram and Moosavi[22]proposed a T-junction with unequal width branches.In addition,a thermally mediated breakup method in which aqueous droplets were manipulated by thermally induced surface tension gradients[23]and a valve-based method that relies on the breakup of droplets in a T-junction with a valve in one of the minor branches[24]were adopted.

All droplet microfluidic techniques rely on transport of droplets in channels[5].Typically,the hydrodynamic resistance of the flow,which is the final expression for the flow field of the continuous phase,depends on the speed of the confinement droplet in a microchannel[25].The droplet speed can also incorporate variation to the inside velocity field.For example,Harshe et al.[26]investigated the dependence on droplet mixing,which is highly controlled by the droplet velocity field,the droplet traveling velocity in a microfluidic serpentine channel,and then assessed principles of the droplet velocity influencing the mixing time.In addition,the droplet velocity was shown to be significant to the droplet fusion and breakup[21,27–29].

With the investigations of droplet size and velocity,the basic affecting factors,as well as the manipulating methods of the dropletsize and speed are achieved in the microfluidics.Based on those,the size tunable droplets emerging system,which can provide a wide droplet sizes range to match the application could be further studied.In this work,the effects of the forced deformation of the bottom wall during the droplet generation in a microchannel with amain width of 50 μm have been depicted.A new protocol to tune the droplet size and velocity is proposed,where the droplet size distribution law is scaled.Several typical shapes of the channel during the bottom wall fluctuation process are investigated to figure out the effects of the shape without the fluctuation deformation.To indicate the effect of the movable bottom wall,different fluctuation models are involved in the numerical simulation.

2 Calculation model and parameters

In this study,three dimentional(3D)geometry is used to investigate the effects of bottom wall deformation on droplet generation in the T-junction microchannel.The immiscible liquids are pumped into the microchannel through the water inlet and the oil inlet,respectively,and then merge at the right angle to generate droplets(Fig.1a),where water is set as the dispersed phase and perfluoropolyether oil as the continuous phase.The channel length is 960 μm in the x-direction,the lengths of both inlets are L o=L w=200 μm,and the channel height is 50μm.Widths of the side channel and main channel are ws=25 μm and w=50 μm respectively.

In order to figure out the influence of the bottom wall deformation, fluctuation is applied at a specific area on the bottom wall of main channel’s downstream.The fluctuation’s period and amplitude are varied to investigate their effects on droplet generation.The basic shape of the fluctuation can be described by Eq.(1),which is the deformation equation of the thin plate with four boundaries fixed under a uniform load.

where wz is the z-direction deformation,x,y are the x-direction and y-direction positions of the bottom wall respectively as shown in Fig.1b,D is the bending stiffness of the deformable bottom wall,l is the length of the fluctuant area on the wall(l=100 μm),q0 is the uniform load which varies with time t and is controlled by the period T and load factor q f with the following equation.

(4)水行政主管部门。加强立法和制度措施的制订。加强水资源论证管理工作的立法工作,完善水资源论证管理法律体系,加强制度措施的制订,严肃水资源论证程序。以法律为主,制度为辅,保障使水资源论证工作实现有法可依、执法必行的局面,确保水资源论证工作的顺利施行。

The deformation wz varies with the location(x and y)and reaches its maximum shown as A at the center of the deformable wall in Fig.1b.Here,q f=81700 is adopted in the basic case to achieve the maximum deformation of the wall A=8 μm.It can be described in short that the deformable part of the bottom wall moves in the shape of Eq.(1)with time correlation of Eq.(2).

3 Validation of numerical simulation

Before the deformable wall channel is investigated,droplet generations in the normal solid channel with the same general geometry are analyzed and compared with the experimental and numerical methods,to validated the numerical method in this study.

More details of the droplet generation results are shown in Fig.4,including the average droplet length,the length difference and droplet generation time(frequency).As shown in Fig.4a,the average length is distributed around 76.1 μm within an error less than±0.3%.Equally,the droplet generation time brings out a similar trend as shown in Fig.4b since all cases are driven with the same liquid flow rate.Furthermore,for all the droplet emerging processes in this part,the droplet emerging frequencies are relatively stable at about 104 Hz, although the droplet size occupies a wide range about 73–82 μm, which indicates that the fluctuation of the bottom wall has no effect on the droplet generation time.

Fig.1 Models of the T-junction microchannel.a Schematic of liquids flowing to generate droplets in the T-junction microchannel;b details of the bottom wall deformation are shown in dashed black ellipse of a.A is the maximum height of the dome,which is also the amplitude of the fluctuation in the deformable wall cases.When A is negative,zero and positive,the wall of the channel is accordingly expanded,straight and reduced;c schematic of oil–water contact angle θow;d details of meshes in case 4 where the boundary layer near walls of the main channel is incorporated;e details of meshes in case 6 where the local refinement of 1.25 μm in the x-direction of the junction in the main channel is incorporated with a smooth transition between two grid sizes

Table1 Properties of the immiscible liquids

Liquids Density(kg·m−3) Viscosity(Pa·s) Interfacial tension(N·m−1)Water 998.2 0.001003 0.04 Oil 1880 0.11468

3.1 Experimental and numerical methods

In this model,liquid properties are shown in Table 1.The experimental results with the polydimethylsiloxane(PDMS)T-shape channel are adopted to validate the simulation results.The general geometry,liquids and flow conditions are the same as the simulation.Liquids are pumped by syringe pumps(KD Scientific,KDS-100)at a flow rate of 0.04 m L/h for water and 0.08 m L/h for per fluoropolyether oil.The shadow processes of all the flows are recorded by a highs peed camera(Phantom v7.3,Vision Research Inc.).

The droplet sizes show some patterns along the droplet sequence as shown in Fig.3,which seems like a periodic process.Motivated by these observations,this pattern is documented to infer the scaling law of the droplet size with the deformable bottom wall.

3.2 Mesh independence validation

For droplet generation in the T-junction microchannel,the droplet length is controlled by the liquid flow rates and the geometry parameters,especially the channel’s width[6],in the way that the droplet length grows with increasing width of the channel.Based on the fluctuation curve of Eq.(1),nine types of wall shapes,which are varied by setting the A(Fig.1b)to different values separately called A f without time-dependent deformation,are involved to analyze the effects of the shapes on the emerging droplet.In addition,a rationality of the 3D simulation is provided to get trends of the droplet generation with different corresponding width before investigation of the bottom wall fluctuation.

As mentioned in Eq.(3),the droplet length has a corelationship with the height of the channel.In the varied fluctuation period cases,the amplitude A,which corresponds to A f in the fixed shape channel in the former section,varies by the following equation in the basic case:

When the grid number remains,the mesh type is verified by adopting the boundary layer near the wall(Fig.1d)as shown in case 3 and case 4.About 6.1%variation of length in case4 is induced as compared with case3,which is derived from the mesh type(the boundary layer)because grid numbers are the same in every corresponding orientation.As the grid number of the cross-section produces a significant variation,case 5 is involved and finally yields a similar result with case 4,even though the mesh accuracy in y-direction and z-direction is improved by 50%and the cell amount is doubled.It is indicated that the meshes of case 4 and case 5 are proper in this type of simulation and the mesh independence is valid.To confirm result,case 6 is also simulated because the junction part has relation with the two liquids flow and the interface evolution.Taking computing performance into consideration,case 5 is chosen in this study.In addition,the experimental result(as shown in case 7)is also used to give some proof that case 5 is closer to real flow and simulations are in good agreement.The continuous phase film near the wall,which can be solid-like[32]and result in a longer droplet,is obviously formed in the experiment,but hard to obtain in the simulation as shown in Table 2.

4 Drop let generation under different fixed shapes of the bottom wall

Mesh accuracy and independence should be validated first to analyze the simulation reasonably.Results of meshes with different meshing types and grid numbers are compared in Table 2.Mesh sizes in the x-direction are documented in case 1 and case 2.Even though the size is reduced from 5 to 1.25 μm with a fourfold improvement in the accuracy and increase in the cell number,corresponding variations,which are about 7.1%in droplet length,1.6%in frequency and a slight difference of the droplet shape near the the simulations.The accuracy of y-direction and z-direction in case 3 is improved from 5 to 2.5 μm,and the cell amount is also increased to 90000 equally to case 2,but the droplet length has a 12.0%extension and the frequency has an 8.7%reduction as compared to case 1,which varies more than case 2.Therefore,the simulation is more sensitive to the grid size in the cross-section of the channel(y-direction and z-direction)since the larger variation is involved under this mesh accuracy improvement.

该合作社建立“粮食银行”的主要目的是想让农民把粮食卖给合作社,合作社能够获得稳定的优质粮源,加工优质大米出售。这个“粮食银行”主要有以下几个特点:

During the droplet flow,the monitored barycenter location can be transferred to the barycenter velocity,which is also the velocity of the droplet.Under different periods of the forced bottom wall fluctuation,the droplet velocity evaluates the curve with the same trend as the movement of the wall,which can also obtain some evidence from Fig.4d that the flow rate of the downstream cross-section(x=0.1 mm)follows the movement of the wall fluctuation.As shown in Fig.6,the velocity changes in the form of sine whose period is derived from fitted velocity curves is completely consistent with the moving bottom wall.The velocity plot is relatively flat about 17000–20000 μm/s when the fluctuation period is T=0.040s as shown in Fig.6a.As the period is reduced,the droplet velocity grows gradually within a wider varying range and when T is reduced to 0.0096 s,the velocity has a lower limit about 14000 μm/s and an upper limit about 23000 μm/s.

where A f is the maximum height of the dome in the fixed shape channel corresponds to A shown in the Fig.1b in the deformable wall cases and can take values in the range of−16 to 16 μm in this simulation and correlation.Although the maximum deformation near the junction wall varies from−16 to 16 μm,which is from −32%to 32%of the main channel height,the droplet length is only changed within 7%.

Since the flow of the dispersed and continuous phases are uniform separately, the droplet emerging frequency grows to cover the reduction of the droplet size.Besides,the shape affects the droplet size chiefly,but not the monodispersity of the droplet.

5 Drop let length under different bottom wall fluctuations

To determine the influence of fluctuation,the period and amplitude,which are the key parameters of the wall shape during the droplet generation in the channel,and the droplet generation results,including the droplet length and droplet emerging frequency,are documented under these different boundary conditions.

Fig.2 The droplet length and droplet emerging frequency of various bottom wall shapes.a Droplet length;b droplet average emerging frequency

5.1 Droplet length distribution under different time periods

To document the influence of time periods on the droplet generation, fluctuations with periods in the range of 0.0096–0.0800 s are simulated.Droplet length distributions during the droplet emerging sequence are shown in Fig.3.Lengths of cases with the larger fluctuation period occupy a wider range of the droplet length,which means a worse monodispersity.The droplet sizes of the fixed shape of A f= −8,8μm,which are simulated in last section,become the two stable limits for lengths in all fluctuating cases.A ll the droplets of different fluctuation periods in this study are seated between the size of the A f=8 μm case and that of the A f= −8 μm case.The reason could be that for these two cases,these are the limited shapes of the bottom wall in the fluctuation,that is, the movable bottom wall is deformed between them to construct the fluctuation. Although the droplet size can fill the range between the two limited lengths, the droplet emerging frequency does not perform discriminately, but stays at 103.8 Hz,which means that droplets of different sizes emerge in a similar time interval.

The liquid flow process is simulated with the commercial software FLUENT(Ansys Inc.USA)using the volume of fluid(VOF)method to track the interface.The hydrophobic boundary condition is applied to all solid walls with an oil–water contact angle of 140°.The non-slip boundary is set for all walls,which means the velocity normal to the cross-section is zero for the continuous phase(the oil flow)at the wall.The fluctuation of the bottom wall,which is described by Eqs.(1)and(2),is applied by the dynamic mesh module with the user defined file at the designed location.The liquid flow in the microchannel is laminar with the velocity inlet condition of 0.00889 m/s(0.04 m L/h)for the water inlet,0.00889 m/s(0.08 m L/h)for the oil inlet,the pressure outlet condition of 0 Pa and the reference pressure of 101325Pa for the outlet.In the numerical simulation,the flow is modeled using the Navier–Stokes equation and the continuity equation.Simulation methods,including the pressure-implicit with splitting of operators(PISO)method,the geo-reconstruct scheme(piecewise-linear interface construction,PLIC),and the control-volume-based technique are adopted to discretize the governing equations and algebraic equations to perform the surface tension calculations,which are mentioned in the literature[30,31].

P∈Rn×d表示用户特征矩阵,第u个行向量pu对应用户u的一个d维特征描述;Q∈Rd×m表示项目特征矩阵,第i个列向量qi对应项目i的一个d维特征描述,d通常取低维.Bu=[b1,b2,…,bn]′是用户评分的偏置向量,其分量bu对应用户u的评分偏好相对于所有用户对所有项目给出的评分的平均分数μ的偏置;Bi=[b1,b2,…,bm]′是项目得分的偏置向量,其分量bi对应项目i的得分相对于μ的偏置.Iu记录u给出了评分的项目的集合,|Iu|表示u给出评分的项目数.对每个j∈Iu,d维向量yj对应了基于u对j的评分行为而得到的隐式特征反馈,Y=I1∪I2∪…∪In是所有隐式反馈向量集合之并.

When we suppose the Eq.(3)can also be used to predict the droplet length in the fluctuating cases,the length can be derived from Eqs.(3)–(5).

where tx is the shift to achieve that the time t of the first recorded droplet emerged is set to be 0.The droplet generation time t could not be continuous but an arithmetic progression with the droplet emerging period as the common difference.So droplet lengths can be a sequence of discontinuous values seating in the predicted equation and varying with the generation time.

Droplet length sequences under various fluctuation are fitted with the form of Eq.(6),which results from the assumption of Eq.(5).

where L ave is the average length of the droplet sequence,A L is the length distribution amplitude,and T L is the length distribution period.As shown in Table 3,the fitted period completely matches the wall moving period in most cases,the average length of droplets has a narrow error with the supposed Eq.(6),which means the droplet length in the fluctuation channel is mainly controlled by the apparent heightof the channel at the emerging time.From the basic relation of the droplet length[7],it is mainly affected by the geometry and flow rates,so the fluctuation should have little influence on flow rates of liquids during the droplet emerging process,which is also proved in Fig.4d that the flow rate of the upstream(x=0 mm)is not changed by the fluctuation.

Table3 Fitted factors of the droplet length distribution

T(s) L ave(μm) T L(s) A L 0.0160 76.200 0.0242 1.061 0.0192 76.027 0.0193 1.018 0.0200 76.126 0.0200 2.028 0.0288 76.039 0.0287 2.559 0.0300 76.094 0.0300 2.471 0.0384 76.027 0.0384 2.478 0.0400 76.310 0.0400 2.521 0.0500 76.086 0.0500 2.444 0.0600 76.090 0.0600 2.441 0.0700 76.089 0.0700 2.240 0.0800 76.180 0.0800 2.366 Eq.(6)76.889–2.176

Predicted droplet lengths of a recorded droplet sequence are plotted in Fig.3(the common difference is about 0.0096 for this investigation),with which the numerical lengths are compared.When the fluctuation period is relatively long(which means the moving speed of the deformable wall is slow),such as the cases with the moving period T above 0.0192 s,the numerical droplet length meets Eq.(6)with low standard errors.For the short period fluctuation,the deformable wall affects the droplet emerging process in a way that the droplet length does not follow the channel width relation Eq.(6),as well as the long period cases.In addition,the distribution range of the numerical droplet length will be reduced as shown in Fig.3e and f where the period is under about 0.020 s.

“今儿个呢,王爷?”晚饭毕,眼看着王爷已跨出饭馆门,戏班子几个人赶紧起身,用手胡乱抹几下嘴,追了出去。

Fig.3 Curves of droplet length via relative flow time under different fluctuation periods.The solid circle stands for the numerical result,the circle stands for the predicted result from the fitted equation and the line is the plot of Eq.(6)under different periods.a A max = 8 μm,T = 0.060 s.b A max=8 μm,T=0.04 s.c A max=8 μm,T=0.0288 s.d A max=8 μm,T=0.020 s.e A max=8 μm,T=0.0192 s.f A max=8 μm,T=0.016 s

Fig.4 Droplet generation performance under the fluctuation wall.a Average droplet length,the solid line is the average length of all cases.b Average droplet generation frequency,the solid line is the average frequency of all cases.c Droplet length difference,which is the amplitude of the fitted relation of the droplet length under different periods and also is half of the difference between the minimum and maximum values.The dash line is the average length difference in the larger period cases,which is supposed to be the stable magnitude of length difference.The solid line is the limitation from the fixed shape cases.d The x=0 mm and x=0.1 mm lines stand for the flow rate of cross-section of the main channel at the position of x=0 mm and x=0.1 mm,which can be called Q0 and Q0.1,respectively.The net flow rate is Q0−Q0.1.The fluctuation of A is the variation of A with time during the droplet process

Fig.5 Droplet sizes under various amplitude of the bottom wall fluctuation

The documented length difference is the amplitude of the fitted plots in Fig. 3,which is half of the real length distributed range in that case. The length difference shows variation with the changing fluctuation period(Fig.4c),which is also indicated in Fig.3.That is,the length difference grows with the fluctuation period when T is below about 0.028 s,and then is stable at a value less than the half difference of droplet lengths between two limited shapes in the fixed wall cases.

仓库可以定期地进行打扫,但是那么多的地贴码却无法每次运作时都去打扫一遍。针对灰尘较多、场地较大、打扫比较困难的问题,可以在机器人设计上进行改进。在机器人下方靠近地贴码探头附近安装一个小型的清扫装置,每次经过一个地贴码时就会进行一遍快速的清扫。这样即使周边灰尘再多,也不会影响机器人的运作了。

5.2 Droplet length under various amplitude of bottom wall fluctuation

Three models with a droplet generation time period of 0.0096 s are simulated under the amplitudes of A=8,12,16 μm,respectively.The droplet length distribution and average droplet sizes of these models show the differences in the way that the larger fluctuation amplitude produces the larger droplet(Fig.5).For the fluctuating cases,droplets generate unstably at the beginning of the flow and become a stable generation gradually(Fig.5)with a long time average length,as shown by the dash line in Fig.5,and the emerging time is still kept at about 0.0096 s.In other words,the droplets are stable at some length after a developing process of the droplet generation when the bottom wall movement is applied on the channel. Otherwise, the droplet will keep the size since the beginning of the droplet flow as shown by the non-fluctuation case. The difference of the droplet size under various fluctuation amplitudes is slight within a relative difference of 4% of the droplet length, where the A = 8 μm case is used as the basic one mentioned in Sect. 2.

这里,我们为大家介绍四种适合用来处理人像照片的色彩效果——当然了,用来处理任意一种类型的照片也没有问题。这些色彩风格最棒的一点就是它们虽然最适合于人像,但并不局限于此。

For the droplet length, the bottom wall fluctuation has limited effects and the length variation(the difference between the maximum and minimum value)is less than 7%of the average length,which is almost the variation range of the two limits cases.That is,the fixed shape cases of the same wall moving amplitude where the local height difference is 32%of the channel’s width.Fluctuations involved in this study have a gentle influence on the distribution of droplet sizes,even though the lengths fall on a predicted equation.In addition,when the fluctuation has the same frequency with the droplet generation,the monodispersity of the length is not affected,but the average length varies slightly.

6 Drop let velocity under different bottom wall fluctuation

When the bottom wall is deforming,the volume of the channel is varied,and the obvious effect is to pump the liquid flow inside the channel and finally change the liquid velocity.The droplet velocity is a typical parameter in the flow,which is indicated by the velocity of the droplet’s barycenter.In the constricted droplet flow,the droplet is considered to be symmetric,so the barycenter of the droplet should lie in the midsection of the channel.The droplet curve in the midsection of the z-direction is used to access and monitor the barycenter’s movement during the droplet flow,and finally the droplet velocity is documented.

石油大数据技术已经在石油工业的上中下游各个领域中体现出体现非常重要的支撑力和推动力,能够有效解决数据管理部门多、数据存储系统分散、数据应用领域广等诸多问题,石油大数据技术的广泛应用必将大幅度提高石油企业各个部门之间的数据共享度和业务协同度,同时还能够逐步提升我国智慧油田建设水平。因此,尽快着手构建我国石油大数据技术应用体系进而在数据密集型计算一体化研究方法指导下建立石油大数据云平台是十分必要的。[课题与项目:(1)陕西省社科界2017年度重大理论与现实问题研究项目“陕西‘智慧油田’建设中的油气大数据平台架构研究”

6.1 Effects of time periods on the droplet velocity

Accordingly,the shape of the channel’s bottom wall near the T-junction,which is shown in Fig.1b and supposed to be deformable in the following sections,is fixed to several shapes including that of wall expanded(A f<0),straight(A f=0),and reduced(A f>0).As a result,the height of the channel near the junction is changed and the average droplet length is tuned within a range about 74–82 μm orders in this analysis as shown in Fig.2.This relation shows the consistency with the result of the tunable junction width channel in Ref.[17]and also the scaling law of Garstecki’s that the droplet length is in direct proportion to the channel width[6].For this study,the droplet length L is fitted with A f in Eq.(3).

Fig.6 Evolutions of the average droplet velocity in the main channel from x=150 μm to x=800 μm.a A max=8 μm,T=0.04 s.b A max=8 μm,T=0.0288 s.c A max=8 μm,T=0.020 s.d A max=8 μm,T=0.0192 s.e A max=8 μm,T=0.016 s.f A max=8 μm,T=0.0096 s

Fig.7 Droplet velocities under different fluctuation periods.a Average droplet length;b droplet velocity difference(half of the velocity distribution range)and relative variation

The average droplet velocity and the velocity difference for every fluctuation are also plotted as shown in Fig. 7. The average velocity is relatively uniform and all fall with in a narrow extent of 1% around 18000 μm/s, while amplitudes of the droplet velocity difference increase significantly with the reducing period of the forced movement of the bottom wall.Since the liquid driving flow rate is uniform and the average droplet size is almost stable,as mentioned in Fig.4,the average droplet velocity is similar under different fluctuations to match the uniform dispersed phase flow rate.However,the movement of the bottom wall practices the variation in the instantaneous velocity,that is,the velocity evolves with time.For A max=8 μm cases,when the period varies from about 0.038–0.0096 s,the velocity amplitude which is half of the overall velocity difference increases from about 6%–22%.That means the velocity variation range could be doubled to 44%of the average droplet velocity.

6.2 Effects of amplitude on the droplet velocity

Besides the fluctuation period,the amplitude of the moving bottom wall is also proved to be a significant factor on the droplet velocity as shown in Fig.8.The larger droplet velocity range is obtained for the larger amplitude of the forced fluctuation.When A max=16μm,the velocity amplitude can reach above 40%,where the velocity range during the flow is shown to cover about 10000–26000μm/s.The droplet velocity amplitude increases linearly with the A max of the forced bottom wall movement as shown in Table 4.

Fig.8 Droplet velocity under different fluctuation amplitudes

Table 4 Droplet velocity amplitude under different fluctuation amplitudes

A max(μm) 8 12 16 Velocity amplitude(μm/s) 3881.9 5780.8 7715.9

The droplet flow can be affected by both the period and the amplitude obviously in the way that the higher speed of the moving wall causes the larger velocity variation during the flow.The relevant velocity variation(the difference of the maximum and minimum droplet velocity in the droplet flow)could be up to 80%of the average velocity in this study,whose wall moving amplitude is 32%(A max=16 μm)of the channel’s width and the average droplet length is stable at 78.21μm with a standard error of0.2%.In addition,A max=12 μm and A max=8 μm cases can also achieve a large droplet velocity variation(about 60% and 40%, respectively)with a relatively uniform droplet size.In other words,the fluctuation of the bottom wall can force an obvious velocity variation on the droplet with little change of the generated droplet,when its period is the same as the droplet generation time.

随着我国经济的不断发展,市政工程施工技术也在逐渐提升,采用绿色施工技术不仅可以达到绿色环保的目的,还能提升市政工程的施工质量。市政工程施工对人们的生活有着极大影响,绿色施工技术体现绿色环保理念的同时还能够满足人们对于工程的根本需求,最大程度上避免对周边环境造成影响,促使市政工程施工更加安全和稳定,使我国建筑行业发展的更加迅速,为我国经济发展奠定基础。

第三种最为方便,通过远程办税端。2019年1月1日,个人所得税手机APP和电子税务局PC端(https://its.sdn-tax.gov.cn/)将正式上线,届时,纳税人可以通过手机或电脑填写相关信息,直接可以将信息传递到税务部门,由税务部门网络推送到任职单位的客户端中。

7 Conclusions

In this article,a protocol based on bottom wall fluctuations for the droplet generation and flow control is presented.The bottom wall shapes,the period and the amplitude of wall fluctuations are investigated to document the effects of the deformation and the fluctuation in droplets generation separately.

1.The fluctuations are demonstrated to expand the droplet size range within about 7%but not change the average size and droplet emerging frequency under most of our investigating conditions,while the fixed shape of the bottom wall changes the average droplet size significantly.Droplet sizes could be predicted by the droplet-sequence length formula along the relative time when the fluctuation is forced on the bottom wall near the T-junction.

2.The wall moving fluctuation,whose frequency is the same as that of the droplet generation,can remain the droplet size stable at a specific value.A little variation(within4%)of this longtime steady magnitude is incorporated by the double-fold amplitude of the wall fluctuation after the flow development stage at the beginning.

3.The forced moving wall changes the droplet velocity in a relatively large range in the way that the velocity evolution has the same sine period with the moving wall..When the wall moving frequency is equal to the droplet emerging frequency,the droplet size has a superior stability.Therefore,the large variation range of droplet velocity(up to 80%)can be imposed by the bottom wall fluctuation without altering the droplet size and emerging interval.

Acknowledgements The authors are grateful for the support of the National Natural Science Foundation of China(11572013 and 11702007)and the China Postdoctoral Science Foundation(2017M 610725).

在进行森林抚育操作过程中,对树木进行修剪枝条非常重要,不仅能够有效的改善目标树木的生长情况,还能够提升其发育的状态,更好的让其得到生长。一般来说,在修去树冠下部的1~2轮活枝时,为了不对整个树木的枝干产生一定的伤害,需要对树干的枝桩尽量进行平整的修剪。并且,在对不同生长阶段的林木进行修剪过程中,科学的保留对应的冠长非常重要。

References

1.Niu,X.,Gielen,F.,Edel,J.B.,et al.:A microdroplet dilutor for high-throughput screening.Nat.Chem.3,437–442(2011)

2.Edd,J.F.,Di Carlo,D.,Humphry,K.J.,et al.:Controlled encapsulation of single-cells into monodisperse picolitre drops.Lab Chip 8,1262–1264(2008)

3.Hafermann,L.,Michael Köhler,J.:Small gold nanoparticles formed by rapid photochemical flow-through synthesis using micro fluid segment technique.J.Nanopart.Res.17,99(2015)

4.Rosenfeld,L.,Lin,T.,Derda,R.,et al.:Review and analysis of performance metrics of droplet microfluidics systems.Micro fluid.Nanofluidics 16,921–939(2014)

5.Jakiela,S.,Makulska,S.,Korczyk,P.M.,et al.:Speed of flow of individual droplets in microfluidic channels as a function of the capillary number,volume of droplets and contrast of viscosities.Lab Chip 11,3603–3608(2011)

6.Maan,A.A.,Nazir,A.,Khan,M.K.I.,et al.:Microfluidic emulsification in food processing.J.Food Eng.147,1–7(2015)

7.Garstecki,P.,Fuerstman,M.J.,Stone,H.A.,et al.:Formation of droplets and bubbles in a microfluidic T-junction-scaling and mechanism of break-up.Lab Chip 6,437–446(2006)

8.Li,Z.,Leshansky,A.M.,Metais,S.,et al.:Step-emulsification in a microfluidic device.Lab Chip 15,1023–1031(2015)

9.Maan,A.A.,Schroën,K.,Boom,R.:Spontaneous droplet formation techniques for monodisperse emulsions preparation—Perspectives for food applications.J.Food Eng.107,334–346(2011)

10.Gu,H.,Malloggi,F.,Vanapalli,S.A.,et al.:Electrowetting enhanced microfluidic device for drop generation.Appl.Phys.Lett. 93,183507(2008)

11.Lee,W.,Walker,L.M.,Anna,S.L.:Role of geometry and fluid properties in droplet and thread formation processes in planar flow focusing.Phys.Fluids 21,032103(2009)

12.Steegmans,M.L.J.,Schroën,C.G.P.H.,Boom,R.M.:Generalised insights in droplet formation at T-junctions through statistical analysis.Chem.Eng.Sci.64,3042–3050(2009)

13.Christopher,G.F.,Noharuddin,N.N.,Taylor,J.A.,et al.:Experimental observations of the squeezing-to-dripping transition in T-shaped microfluidic junctions.Phys.Rev.E 78,036317(2008)

14.Xu,J.H.,Li,S.W.,Tan,J.,et al.:Correlations of droplet formation in T-junction microfluidic devices:from squeezing to dripping.Microfluid.Nanofluidics 5,711–717(2008)

15.Lee,C.H.,Hsiung,S.K.,Lee,G.B.:A tunable micro flow focusing device utilizing controllable moving walls and its applications for formation of m icro-droplets in liquids.J.M icromech.M icroeng.17,1121–1129(2007)

16.Pang,Y.,Kim,H.,Liu,Z.,et al.:A soft microchannel decreases polydispersity of droplet generation.Lab Chip 14,4029–4034(2014)

17.Wu,H.W.,Huang,Y.C.,Wu,C.L.,et al.:Exploitation of a microfluidic device capable of generating size-tunable droplets for gene delivery.Micro fluid.Nanofluidics 7,45–56(2008)

18.Wang,J.H.,Lee,G.B.:Formation of tunable,emulsion microdroplets utilizing flow-focusing channels and a normally-closed micro-valve.Micromach.4,306–320(2013)

19.Li,J.,Mittal,N.,Mak,S.Y.,et al.:Perturbation-induced droplets for manipulating droplet structure and configuration in microfluidics.J.Micromech.and Microeng.25,084009(2015)

20.Huang,S.B.,Chang,Y.H.,Lee,H.C.,et al.:A pneumatically-driven microfluidic system for size-tunable generation of uniform cellencapsulating collagen microbeads with the ultrastructure similar to native collagen.Biomed.Microdevices 16,345–354(2014)

21.Link,D.R.,Anna,S.L.,Weitz,D.A.,et al.:Geometrically mediated breakup of drops in microfluidic devices.Phys.Rev.Lett.92,054503(2004)

22.Bedram,A.,Moosavi,A.:Droplet breakup in an asymmetric microfluidic T junction.Eur.Phys.J.E 34,78(2011)

23.Ting,T.H.,Yap,Y.F.,Nguyen,N.T.,et al.:Thermally mediated breakup of drops in microchannels.Appl.Phys.Lett.89,234101(2006)

24.Bedram,A.,Darabi,A.E.,Moosavi,A.,et al.:Numerical investigation of an efficient method(T-junction with valve)for producing unequal-sized dropletsinmicro-andnano- fluidic systems.J.Fluids Eng.Trans.ASME 137,9(2015)

25.Labrot,V.,Schindler,M.,Guillot,P.,et al.:Extracting the hydrodynamic resistance of droplets from their behavior in microchannel networks.Biomicro fluidics 3,12804(2009)

26.Harshe,Y.M.,van Eijk,M.J.,Kleijn,C.R.,et al.:Scaling of mixing time for droplets of different sizes traveling through a serpentine microchannel.RSC Adv.6,98812–98815(2016)

27.Liu,Z.,Cao,R.,Pang,Y.,et al.:The influence of channel intersection angle on droplets coalescence process.Exp.Fluids 56,24(2015)

28.Liu,Z.,Wang,X.,Cao,R.,et al.:Dropletcoalescence at microchannel intersection chambers with different shapes.Soft Matter 12,5797–5807(2016)

29.Liu,Z.,Zhang,L.,Pang,Y.,et al.:Micro-PIV investigation of the internal flow transitions inside droplets traveling in a rectangular microchannel.Micro fluid.Nanofluidics 21,180(2017)

30.Pang,Y.,Liu,Z.,Zhao,F.:Downstream pressure and elastic wall reflection of droplet flow in a T-junction microchannel.Acta Mech.Sin.PRC 32,579–587(2016)

31.Sivasamy,J.,Wong,T.N.,Nguyen,N.T.,et al.:An investigation on the mechanism of droplet formation in a microfluidic T-junction.Micro fluid.Nanofluidics 11,1–10(2011)

32.Yuan,Q.,Zhao,Y.:Precursor film in dynamic wetting,electrowetting,and electro-elasto-capillarity.Phys.Rev. Lett.104,246101(2010)

YanPang,XiangWang,ZhaomiaoLiu
《Acta Mechanica Sinica》2018年第4期文献

服务严谨可靠 7×14小时在线支持 支持宝特邀商家 不满意退款

本站非杂志社官网,上千家国家级期刊、省级期刊、北大核心、南大核心、专业的职称论文发表网站。
职称论文发表、杂志论文发表、期刊征稿、期刊投稿,论文发表指导正规机构。是您首选最可靠,最快速的期刊论文发表网站。
免责声明:本网站部分资源、信息来源于网络,完全免费共享,仅供学习和研究使用,版权和著作权归原作者所有
如有不愿意被转载的情况,请通知我们删除已转载的信息 粤ICP备2023046998号