It is known that contact lines keep relatively still on solids until static contact angles exceed an interval of hysteresis of static contact angle (HSCA), and contact angles keep changing as contact lines relatively ...It is known that contact lines keep relatively still on solids until static contact angles exceed an interval of hysteresis of static contact angle (HSCA), and contact angles keep changing as contact lines relatively slide on the solid. Here, the effects of HSCA and boundary slip were first distinguished on the micro-curvature force (MCF) on the seta. Hence, the total MCF is partitioned into static and dynamic MCFs correspondingly. The static MCF was found proportional to the HSCA and related with the asymmetry of the micro-meniscus near the seta. The dynamic MCF, exerting on the relatively sliding contact line, is aroused by the boundary slip. Based on the Blake-Haynes mechanism, the dynamic MCF was proved important for water walking insects with legs slower than the minimum wave speed . As insects brush the water by laterally swinging legs backwards, setae on the front side of the leg are pulled and the ones on the back side are pushed to cooperatively propel bodies forward. If they pierce the water surface by vertically swinging legs downwards, setae on the upside of the legs are pulled, and the ones on the downside are pushed to cooperatively obtain a jumping force. Based on the dependency between the slip length and shear rate, the dynamic MCF was found correlated with the leg speed U, as , where and are determined by the dimple depth. Discrete points on this curve could give fitted relations as (Suter et al., J. Exp. Biol. 200, 2523-2538, 1997). Finally, the axial torque on the inclined and partially submerged seta was found determined by the surface tension, contact angle, HSCA, seta width, and tilt angle. The torque direction coincides with the orientation of the spiral grooves of the seta, which encourages us to surmise it is a mechanical incentive for the formation of the spiral morphology of the setae of water striders.展开更多
Wetting phenomena are widespread in nature and industrial applications. In general, systems concerning wetting phenomena are typical multicomponent/multiphase complex fluid systems. Simulating the behavior of such sys...Wetting phenomena are widespread in nature and industrial applications. In general, systems concerning wetting phenomena are typical multicomponent/multiphase complex fluid systems. Simulating the behavior of such systems is important to both scientific research and practical applications. It is challenging due to the complexity of the phenomena and difficulties in choosing an appropriate numerical method. To provide some detailed guidelines for selecting a suitable multiphase lattice Boltzmann model, two kinds of lattice Boltzmann multiphase models, the modified S-C model and the H-C-Z model, are used in this paper to investigate the static contact angle on solid surfaces with different wettability combined with the geometric formulation(Ding, H. and Spelt, P.D. M. Wetting condition in diffuse interface simulations of contact line motion. Physical Review E, 75(4), 046708(2007)). The specific characteristics and computational performance of these two lattice Boltzmann method(LBM) multiphase models are analyzed including relationship between surface tension and the control parameters, the achievable range of the static contact angle, the maximum magnitude of the spurious currents(MMSC), and most importantly, the convergence rate of the two models on simulating the static contact angle. The results show that a wide range of static contact angles from wetting to non-wetting can be realized for both models. MMSC mainly depends on the surface tension. With the numerical parameters used in this work, the maximum magnitudes of the spurious currents of the two models are on the same order of magnitude. MMSC of the S-C model is universally larger than that of the H-C-Z model. The convergence rate of the S-C model is much faster than that of the H-C-Z model. The major foci in this work are the frequently-omitted important details in simulating wetting phenomena. Thus, the major findings in this work can provide suggestions for simulating wetting phenomena with LBM multiphase models along with the geometric formulation.展开更多
There are five most widely used contact angle schemes in the pseudopotential lattice Boltzmann(LB)model for simulating the wetting phenomenon:The pseudopotential-based scheme(PB scheme),the improved virtualdensity sch...There are five most widely used contact angle schemes in the pseudopotential lattice Boltzmann(LB)model for simulating the wetting phenomenon:The pseudopotential-based scheme(PB scheme),the improved virtualdensity scheme(IVD scheme),the modified pseudopotential-based scheme with a ghost fluid layer constructed by using the fluid layer density above the wall(MPB-C scheme),the modified pseudopotential-based scheme with a ghost fluid layer constructed by using the weighted average density of surrounding fluid nodes(MPB-W scheme)and the geometric formulation scheme(GF scheme).But the numerical stability and accuracy of the schemes for wetting simulation remain unclear in the past.In this paper,the numerical stability and accuracy of these schemes are clarified for the first time,by applying the five widely used contact angle schemes to simulate a two-dimensional(2D)sessile droplet on wall and capillary imbibition in a 2D channel as the examples of static wetting and dynamic wetting simulations respectively.(i)It is shown that the simulated contact angles by the GF scheme are consistent at different density ratios for the same prescribed contact angle,but the simulated contact angles by the PB scheme,IVD scheme,MPB-C scheme and MPB-W scheme change with density ratios for the same fluid-solid interaction strength.The PB scheme is found to be the most unstable scheme for simulating static wetting at increased density ratios.(ii)Although the spurious velocity increases with the increased liquid/vapor density ratio for all the contact angle schemes,the magnitude of the spurious velocity in the PB scheme,IVD scheme and GF scheme are smaller than that in the MPB-C scheme and MPB-W scheme.(iii)The fluid density variation near the wall in the PB scheme is the most significant,and the variation can be diminished in the IVD scheme,MPB-C scheme andMPBWscheme.The variation totally disappeared in the GF scheme.(iv)For the simulation of capillary imbibition,the MPB-C scheme,MPB-Wscheme and GF scheme simulate the dynamics of the liquid-vapor interface well,with the GF scheme being the most accurate.The accuracy of the IVD scheme is low at a small contact angle(44 degrees)but gets high at a large contact angle(60 degrees).However,the PB scheme is the most inaccurate in simulating the dynamics of the liquid-vapor interface.As a whole,it is most suggested to apply the GF scheme to simulate static wetting or dynamic wetting,while it is the least suggested to use the PB scheme to simulate static wetting or dynamic wetting.展开更多
为设计和优化面粉输送设备,应用离散元法对面粉进行准确地工程建模和分析,需要对其接触参数进行必要的标定。该研究依据颗粒缩放理论,用“Hertz-Mindlin with Johnson-Kendall-Roberts”接触模型表征面粉颗粒间黏性的影响,提出了一种基...为设计和优化面粉输送设备,应用离散元法对面粉进行准确地工程建模和分析,需要对其接触参数进行必要的标定。该研究依据颗粒缩放理论,用“Hertz-Mindlin with Johnson-Kendall-Roberts”接触模型表征面粉颗粒间黏性的影响,提出了一种基于静/动态休止角的接触参数标定方法。运用正交试验方法,对接触参数的敏感性和方差分析,表明面粉颗粒间的滚动摩擦系数、面粉颗粒与不锈钢表面间的静摩擦系数、表面能对静态休止角的影响极显著(P<0.01),并且多组接触参数都可以模拟出与试验相同的静态休止角。进一步研究表明,面粉颗粒与不锈钢表面间的静摩擦系数的合理取值范围为0.2~0.4。通过2种填充率、4种转速下基于动态休止角的参数标定,将其中与试验最为吻合的一组参数作为标定结果,其值如下:面粉颗粒之间恢复系数为0.6、面粉颗粒之间静摩擦系数为0.2、面粉颗粒之间滚动摩擦系数为0.1、面粉颗粒与不锈钢容器表面之间恢复系数为0.6、面粉颗粒与不锈钢容器表面之间静摩擦系数为0.6、面粉颗粒与不锈钢容器表面之间滚动摩擦系数为0.5、表面能为0.12 J/m^(2)。使用该组参数对矩形容器中物料自由坍塌试验进行仿真,其结果与试验结果相符,验证了该标定方法的有效性。该研究提出的标定方法简单、易执行,对粉料输送设备的设计及优化具有一定的工程应用价值。展开更多
为了提高接触角计算的准确性和自动化程度,基于MATLAB的图形用户界面(graphical user interface,GUI)实现了一套动态和静态接触角计算软件系统,实现了手动方式和自动方式获得液滴接触角的方法,自动方式中基于水平集方法自动获得液滴边缘...为了提高接触角计算的准确性和自动化程度,基于MATLAB的图形用户界面(graphical user interface,GUI)实现了一套动态和静态接触角计算软件系统,实现了手动方式和自动方式获得液滴接触角的方法,自动方式中基于水平集方法自动获得液滴边缘,该方法基于偏微分方程,能获得全局最优的液滴边缘,抗干扰能力强,后续处理方便,另外还实现了基于阈值的自动边缘提取算法。接触角计算中实现了切线法、圆拟合法、椭圆拟合法、量角法、量高法、多项式拟合法、Young-Laplace拟合法等算法。分析了典型算法的原理、特点和适用范围。拟合和计算结果能以数字、图像方式显示和存储,同时自动存入数据库。将该软件应用于硅橡胶憎水性检测,结果表明在不同情况下该系统均能自动、准确获得动静态接触角,可在很大程度上节省接触角计算的工作量。展开更多
文摘It is known that contact lines keep relatively still on solids until static contact angles exceed an interval of hysteresis of static contact angle (HSCA), and contact angles keep changing as contact lines relatively slide on the solid. Here, the effects of HSCA and boundary slip were first distinguished on the micro-curvature force (MCF) on the seta. Hence, the total MCF is partitioned into static and dynamic MCFs correspondingly. The static MCF was found proportional to the HSCA and related with the asymmetry of the micro-meniscus near the seta. The dynamic MCF, exerting on the relatively sliding contact line, is aroused by the boundary slip. Based on the Blake-Haynes mechanism, the dynamic MCF was proved important for water walking insects with legs slower than the minimum wave speed . As insects brush the water by laterally swinging legs backwards, setae on the front side of the leg are pulled and the ones on the back side are pushed to cooperatively propel bodies forward. If they pierce the water surface by vertically swinging legs downwards, setae on the upside of the legs are pulled, and the ones on the downside are pushed to cooperatively obtain a jumping force. Based on the dependency between the slip length and shear rate, the dynamic MCF was found correlated with the leg speed U, as , where and are determined by the dimple depth. Discrete points on this curve could give fitted relations as (Suter et al., J. Exp. Biol. 200, 2523-2538, 1997). Finally, the axial torque on the inclined and partially submerged seta was found determined by the surface tension, contact angle, HSCA, seta width, and tilt angle. The torque direction coincides with the orientation of the spiral grooves of the seta, which encourages us to surmise it is a mechanical incentive for the formation of the spiral morphology of the setae of water striders.
基金Project supported by the National Natural Science Foundation of China(Nos.50874071 and51704191)the Shanghai Leading Academic Discipline Project(No.S30106)+1 种基金the Key Program of Science and Technology Commission of Shanghai Municipality(No.12160500200)the PetroChina Innovation Foundation(No.2017D-5007-0209)
文摘Wetting phenomena are widespread in nature and industrial applications. In general, systems concerning wetting phenomena are typical multicomponent/multiphase complex fluid systems. Simulating the behavior of such systems is important to both scientific research and practical applications. It is challenging due to the complexity of the phenomena and difficulties in choosing an appropriate numerical method. To provide some detailed guidelines for selecting a suitable multiphase lattice Boltzmann model, two kinds of lattice Boltzmann multiphase models, the modified S-C model and the H-C-Z model, are used in this paper to investigate the static contact angle on solid surfaces with different wettability combined with the geometric formulation(Ding, H. and Spelt, P.D. M. Wetting condition in diffuse interface simulations of contact line motion. Physical Review E, 75(4), 046708(2007)). The specific characteristics and computational performance of these two lattice Boltzmann method(LBM) multiphase models are analyzed including relationship between surface tension and the control parameters, the achievable range of the static contact angle, the maximum magnitude of the spurious currents(MMSC), and most importantly, the convergence rate of the two models on simulating the static contact angle. The results show that a wide range of static contact angles from wetting to non-wetting can be realized for both models. MMSC mainly depends on the surface tension. With the numerical parameters used in this work, the maximum magnitudes of the spurious currents of the two models are on the same order of magnitude. MMSC of the S-C model is universally larger than that of the H-C-Z model. The convergence rate of the S-C model is much faster than that of the H-C-Z model. The major foci in this work are the frequently-omitted important details in simulating wetting phenomena. Thus, the major findings in this work can provide suggestions for simulating wetting phenomena with LBM multiphase models along with the geometric formulation.
基金sponsored by the National Natural Science Foundation of China under Grant No.52206101Shanghai Sailing Program under Grant No.20YF1431200the Experiments for Space Exploration Program and the Qian Xuesen Laboratory,China Academy of Space Technology under Grant No.TKTSPY-2020-01-01.
文摘There are five most widely used contact angle schemes in the pseudopotential lattice Boltzmann(LB)model for simulating the wetting phenomenon:The pseudopotential-based scheme(PB scheme),the improved virtualdensity scheme(IVD scheme),the modified pseudopotential-based scheme with a ghost fluid layer constructed by using the fluid layer density above the wall(MPB-C scheme),the modified pseudopotential-based scheme with a ghost fluid layer constructed by using the weighted average density of surrounding fluid nodes(MPB-W scheme)and the geometric formulation scheme(GF scheme).But the numerical stability and accuracy of the schemes for wetting simulation remain unclear in the past.In this paper,the numerical stability and accuracy of these schemes are clarified for the first time,by applying the five widely used contact angle schemes to simulate a two-dimensional(2D)sessile droplet on wall and capillary imbibition in a 2D channel as the examples of static wetting and dynamic wetting simulations respectively.(i)It is shown that the simulated contact angles by the GF scheme are consistent at different density ratios for the same prescribed contact angle,but the simulated contact angles by the PB scheme,IVD scheme,MPB-C scheme and MPB-W scheme change with density ratios for the same fluid-solid interaction strength.The PB scheme is found to be the most unstable scheme for simulating static wetting at increased density ratios.(ii)Although the spurious velocity increases with the increased liquid/vapor density ratio for all the contact angle schemes,the magnitude of the spurious velocity in the PB scheme,IVD scheme and GF scheme are smaller than that in the MPB-C scheme and MPB-W scheme.(iii)The fluid density variation near the wall in the PB scheme is the most significant,and the variation can be diminished in the IVD scheme,MPB-C scheme andMPBWscheme.The variation totally disappeared in the GF scheme.(iv)For the simulation of capillary imbibition,the MPB-C scheme,MPB-Wscheme and GF scheme simulate the dynamics of the liquid-vapor interface well,with the GF scheme being the most accurate.The accuracy of the IVD scheme is low at a small contact angle(44 degrees)but gets high at a large contact angle(60 degrees).However,the PB scheme is the most inaccurate in simulating the dynamics of the liquid-vapor interface.As a whole,it is most suggested to apply the GF scheme to simulate static wetting or dynamic wetting,while it is the least suggested to use the PB scheme to simulate static wetting or dynamic wetting.
文摘为设计和优化面粉输送设备,应用离散元法对面粉进行准确地工程建模和分析,需要对其接触参数进行必要的标定。该研究依据颗粒缩放理论,用“Hertz-Mindlin with Johnson-Kendall-Roberts”接触模型表征面粉颗粒间黏性的影响,提出了一种基于静/动态休止角的接触参数标定方法。运用正交试验方法,对接触参数的敏感性和方差分析,表明面粉颗粒间的滚动摩擦系数、面粉颗粒与不锈钢表面间的静摩擦系数、表面能对静态休止角的影响极显著(P<0.01),并且多组接触参数都可以模拟出与试验相同的静态休止角。进一步研究表明,面粉颗粒与不锈钢表面间的静摩擦系数的合理取值范围为0.2~0.4。通过2种填充率、4种转速下基于动态休止角的参数标定,将其中与试验最为吻合的一组参数作为标定结果,其值如下:面粉颗粒之间恢复系数为0.6、面粉颗粒之间静摩擦系数为0.2、面粉颗粒之间滚动摩擦系数为0.1、面粉颗粒与不锈钢容器表面之间恢复系数为0.6、面粉颗粒与不锈钢容器表面之间静摩擦系数为0.6、面粉颗粒与不锈钢容器表面之间滚动摩擦系数为0.5、表面能为0.12 J/m^(2)。使用该组参数对矩形容器中物料自由坍塌试验进行仿真,其结果与试验结果相符,验证了该标定方法的有效性。该研究提出的标定方法简单、易执行,对粉料输送设备的设计及优化具有一定的工程应用价值。
文摘为了提高接触角计算的准确性和自动化程度,基于MATLAB的图形用户界面(graphical user interface,GUI)实现了一套动态和静态接触角计算软件系统,实现了手动方式和自动方式获得液滴接触角的方法,自动方式中基于水平集方法自动获得液滴边缘,该方法基于偏微分方程,能获得全局最优的液滴边缘,抗干扰能力强,后续处理方便,另外还实现了基于阈值的自动边缘提取算法。接触角计算中实现了切线法、圆拟合法、椭圆拟合法、量角法、量高法、多项式拟合法、Young-Laplace拟合法等算法。分析了典型算法的原理、特点和适用范围。拟合和计算结果能以数字、图像方式显示和存储,同时自动存入数据库。将该软件应用于硅橡胶憎水性检测,结果表明在不同情况下该系统均能自动、准确获得动静态接触角,可在很大程度上节省接触角计算的工作量。
