黏度对流固界面滑移影响的试验研究

为了研究微纳米间隙下固液界面间流体的流动及输运特性,本文改进了商用的原子力显微镜,并利用其对微纳米间隙下固液界面的边界滑移现象进行了试验研究,重点考察了流体黏度对边界滑移的影响.固体壁面样品采用Si(100)表面,试验液体采用不同黏度的去离子水和蔗糖溶液.结果表明,Si(100)表面与去离子水和蔗糖溶液作用会发生边界滑移,而且随着溶液黏度的升高,滑移长度也随之升高,表现了边界滑移与流体黏度的相关性.所得结果对于微流体输运与控制有重要的理论意义与实际价值.     

生物芯片压电微流体泵液-固耦合系统模态分析

对压电微流体泵粘性流体周期流动进行厚度积分平均近似,得到包含粘性的,非线性浅水波动方程,并采用有限元法得到微泵液体压强矩阵方程．液体压强矩阵方程和压电硅片振动有限元方程耦合,得到一个包含微泵进出口扩散管的液-固耦合系统振动方程．液-固耦合系统的模态分析结果表明,做泵液-固耦合系统的自然频率比不耦合的硅片振动自然频率低很多．随着微泵厚度的减少,液体附加质量和粘性阻尼对耦合系统自然频率的影响更加明显．同时发现,对应的压电片振型函数在液-固耦合前后没有明显变化,还给出硅片-阶模态的振幅-频率特征曲线,对薄型无阀压电微流体泵,浅水波模型合理地表达了微泵液体流动和压电硅片振动的相互作用,以及液体附加质量和粘性阻尼对微泵液-固耦合系统动力特征的影响。     

Z切石英在氟化氢铵溶液中的腐蚀特性

对Z切石英在氟化氢铵溶液中的腐蚀特性进行了研究。首先研究了Z向腐蚀速率随腐蚀温度和腐蚀液浓度的变化关系,然后研究了Z向表面腐蚀粗糙度随腐蚀深度、腐蚀液体浓度、腐蚀温度的变化规律,最后将石英在氟化氢铵溶液中的腐蚀特性跟在BHF溶液中的腐蚀特性进行了简单对比。试验结果表明:腐蚀速率随腐蚀温度和腐蚀液浓度增加而增加;Z向腐蚀表面粗糙度随腐蚀时间的增加而增加;提高腐蚀液体浓度有利于减小腐蚀表面粗糙度;提高腐蚀液体温度有利于减小表面粗糙度;石英在氟化氢铵溶液腐蚀具有更高的腐蚀效率和更小的腐蚀表面粗糙度。本研究结果能够为石英MEMS器件的设计和工艺提供有益帮助。     

离散型湍流多相流动的研究进展和需求

离散型多相流动,指气体-颗粒(气-固)、液体-颗粒(液-固)、液体-气泡、气体-液雾以及气泡-液体-颗粒等两相或三相流动.这种类型的多相流动广泛存在于能源, 航天和航空, 化工和冶金,交通运输, 水利, 核能等领域.本文阐述了离散型多相流动的国内外基础研究,包括颗粒/液滴/气泡在流场中受流体动力作用力的研究, 颗粒-颗粒,液滴-液滴,气泡-气泡之间以及颗粒/液滴和壁面之间碰撞和聚集规律的研究,颗粒-气体和气泡-液体湍流相互作用的研究, 和数值模拟的研究,包括多相流动的雷诺平均模拟、大涡模拟和直接数值模拟的研究进展.最后, 归纳了目前尚待研究的需求.      

磨合过程中液体润滑实现超低摩擦的混合模型研究

针对球-盘滑动试验,在磨合过程中获得超低摩擦的液体润滑状态,建立耦合流体润滑、粗糙接触力学、Archard磨损方程和相关物理参数(液体黏度、表面粗糙度和磨损系数)时变函数的混合模型,研究磨合过程中液体润滑的摩擦系数演化.通过数值模拟结果可知：在磨合过程中,润滑介质等效黏度增大,形成流体动压润滑薄膜,有效隔开粗糙表面;其次在磨合过程中,新生成的表面粗糙度降低,减少粗糙峰承载比,实现超低摩擦润滑状态;最后在适当的液体黏度和提高表界面效应减少边界摩擦系数,可进一步实现液体超低摩擦润滑状态.为磨合过程宏观液体润滑性能演化所建立的混合数值模型对提高液体润滑超低摩擦设计效率具有重要价值意义.     

液体界面现象及表面张力的分子动力学模拟

本文利用分子动力学计算机模拟的方法,研究了氪的固-液、液-气和气-固的界面现象。显示了液体的零级径向分布是有序的,与固体相邻的液体部分有类晶结构,固体表面出现了气体吸附层。表面张力计算结果与实验值很好的吻合,证明Lennead-Jones势对球形原子的相互作用作了很好的描述,在表面张力计算中可以忽略去截断效应。     

NUMERICAL SIMULATION OF DRAG REDUCTION CAUSED BY PATTERNED WALL

用加入近邻相互作用势的格子玻尔兹曼方法模拟了表面修饰了去润湿性颗粒的二维管道内流体的流动. 结果显示,管壁的粗糙度、修饰颗粒的润湿性、管内压强以及管内流体黏滞性等都是影响管内流体流速的因素. 在一定压强下,管壁修饰疏水性的颗粒可以明显减小阻力,模拟结果也解释了开采地下油藏时存在启动压力梯度的原因.     

SIMULATION STUDIES OF VISCOSITIES OF Cu-H2O NANOFLUIDS BASED ON COARSE GRAINING WATER MOLECULES

分子动力学模拟是研究纳米流体的黏度特性的重要手段,但计算量庞大. 文章通过对基液水分子粗粒化,使得计算量大幅度减小,且计算精度与全原子模拟相当. 基于平衡态分子动力学,模拟研究了Cu-H₂O 纳米体系的微观运动特性,通过格林- 库博(Green-Kubo) 公式对Cu-H₂O 纳米流体的黏度进行了模拟计算,并考察了温度、体积分数、粒径和颗粒形状对于Cu-H₂O 纳米流体黏度的影响,对已有的悬浮液黏度经验公式进行了修正.     

湍流横掠固定颗粒的全尺度直接数值模拟

气固两相流模拟中,当固相尺度接近或大于Kolmogorov尺度时,普通的点源模型将不再适用,固体相的体积效应和表面效应将对流体相产生显著的影响。通过采用直接数值模拟方法,结合内嵌边界方法对湍流中不同湍流强度流体横掠大于Kolmogorov尺度的固相颗粒进行了全尺度模拟,讨论分析了在两种湍流度下方形颗粒对湍流的调制影响以及颗粒的受力情况。     

分子动力学模拟纳米镍单晶的表面效应

对单晶镍纳米丝、纳米薄膜零温准静态拉伸破坏过程进行了分子动力学模拟．模拟表明表面效应对单晶纳米材料的原子运动及整体力学行为有显著影响．自由表面增加纳米材料的塑性、降低其强度,影响纳米材料的变形机制．受表面效应的作用,纳米镍丝强度与弹性模量均低于纳米镍薄膜．纳米薄膜的断裂接近脆性断裂,断裂强度符合Griffith理想晶体脆断理论;纳米镍丝在断裂过程中表现出微弱塑性．     

Spontaneous Inertial Imbibition in Porous Media Using a Fractal Representation of Pore Wall Rugosity

Considering the separable phenomena of imbibition in complex fine porous media as a function of timescale, it is noted that there are two discrete imbibition rate regimes when expressed in the Lucas–Washburn (L–W) equation. Commonly, to account for this deviation from the single equivalent hydraulic capillary, experimentalists propose an effective contact angle change. In this work, we consider rather the general term of the Wilhelmy wetting force regarding the wetting line length, and apply a proposed increase in the liquid–solid contact line and wetting force provided by the introduction of surface meso/nanoscale structure to the pore wall roughness. An experimental surface pore wall feature size regarding the rugosity area is determined by means of capillary condensation during nitrogen gas sorption in a ground calcium carbonate tablet compact. On this nano size scale, a fractal structure of pore wall is proposed to characterize for the internal rugosity of the porous medium. Comparative models based on the Lucas–Washburn and Bosanquet inertial absorption equations, respectively, for the short timescale imbibition are constructed by applying the extended wetting line length and wetting force to the equivalent hydraulic capillary observed at the long timescale imbibition. The results comparing the models adopting the fractal structure with experimental imbibition rate suggest that the L–W equation at the short timescale cannot match experiment, but that the inertial plug flow in the Bosanquet equation matches the experimental results very well. If the fractal structure can be supported in nature, then this stresses the role of the inertial term in the initial stage of imbibition. Relaxation to a smooth-walled capillary then takes place over the longer timescale as the surface rugosity wetting is overwhelmed by the pore condensation and film flow of the liquid ahead of the bulk wetting front, and thus to a smooth walled capillary undergoing permeation viscosity-controlled flow.     

Experimental simulation of capillary effect on rough flat surfaces

The movement of the liquid column inside the slit was utilized to experimentally simulate the characteristics of the capillary force per unit length for different rough flat surfaces. The movement of the liquid column was achieved by continuously changing the slit interval. The maximum climb height and contact angle of the liquid column were observed during this process to study the relationship between capillary force and contact surface roughness. Based on the assumption that the microstructures on the rough surfaces are of the same form and continuously and uniformly distributed, it is shown that the capillary force per unit length under homogeneous wetting is independent of the roughness. For heterogeneous wetting, the capillary force per unit length is positively correlated with the roughness. The results also indicate that the appearance of “contact line pinning” is caused by insufficient capillary force along the direction of liquid column movement.     

Numerical Analysis of Coupled Heat and Mass Transfer Phenomena in Concrete at Elevated Temperatures

In this work, we present a novel methodology for incorporating the effect of fibre surface morphology on liquid water transport in polymer electrolyte membrane fuel cell gas diffusion layers (GDLs). Roughness features presented on the surface of the fibre are analysed using atomic force microscopy and are found to significantly impact the capillary pressure of liquid water pathways propagating through the GDL. A threshold capillary pressure was defined as the largest capillary pressure exhibited by the liquid water phase during the invasion of the throat. The threshold capillary pressures observed in the presence of roughness features are significantly greater than those in the absence of roughness features. Two-dimensional circumferential roughness models in cylindrical and converging-diverging throats are established, and an interfacial meniscus advancing algorithm is presented to determine the resulting threshold capillary pressures required for liquid water penetration. Revised Young–Laplace equations, which are particularly useful for pore network modeling, are suggested for calculating threshold capillary pressures that account for the effect of the roughness of throats with intrinsic contact angles greater than \(90^{\circ }\).     

A Variational Model of Disjoining Pressure: Liquid Film on a Nonplanar Surface

Variational methods have been successfully used in modelling thin liquid films in numerous theoretical studies of wettability. In this article, the variational model of the disjoining pressure is extended to the general case of a two-dimensional solid surface. The Helmholtz free energy functional depends both on the disjoining pressure isotherm and on the shape of the solid surface. The augmented Young–Laplace equation (AYLE) is a nonlinear second-order partial differential equation. A number of solutions describing wetting films on spherical grains have been obtained. In the case of cylindrical films, the phase portrait technique describes the entire variety of mathematically feasible solutions. It turns out that a periodic solution, which would describe wave-like wetting films, does not satisfy Jacobi’s condition of the classical calculus of variations. Therefore, such a solution is nonphysical. The roughness of the solid surface significantly affects liquid film stability. AYLE solutions suggest that film rupture is more likely at a location where the pore-wall surface is most exposed into the pore space, and the curvature is positive.     

Flow and heat transfer in the wake of a surface-mounted rib with a slit

An experimental study of flow and heat transfer downstream of a surface-mounted rib with a slit is reported. The open area ratios of the slit rib considered are 10, 20, 30, 40 and 50% with respect to the total projected rib area. Experiments were conducted in a wind tunnel, mostly at a hydraulic diameter based Reynolds number of 32,100. The surface Nusselt number distribution was determined by liquid crystal thermography. Results show that the slit inside the rib enhances heat transfer and reduces pressure penalty, with an optimum performance seen at an open area ratio of 20%. To explain this result, a qualitative picture of the flow field behind the rib was obtained by smoke visualization. Time averages and turbulent statistics of the velocity and temperature fluctuations were measured in detail, using hotwire anemometry and cold wire anemometry respectively. For open area ratios less than 30%, measurements show that the flow through the slit modifies the reattaching shear layer from the top of the rib. The resulting reattachment length is smaller, the peak in Nusselt number is higher, and the average heat transfer from the heated surface is enhanced. For the rib with an open area ratio greater than 40%, the lower portion behaves as an independent small rib with its own reattachment region. Simultaneously, the flow downstream of the upper rectangular part shows characteristics of vortex shedding. Thus, the size of the slit is seen to be an additional parameter that can be used to control heat transfer from the solid surface, in comparison to the solid rib.     

Three-Dimensional Lattice Boltzmann Simulations of Single-Phase Permeability in Random Fractal Porous Media with Rough Pore–Solid Interface

Single-phase permeability k has intensively been investigated over the past several decades by means of experiments, theories and simulations. Although the effect of surface roughness on fluid flow and permeability in single pores and fractures as well as networks of fractures was studied previously, its influence on permeability in a random mass fractal porous medium constructed of pores of different sizes remained as an open question. In this study, we, therefore, address the effect of pore–solid interface roughness on single-phase flow in random fractal porous media. For this purpose, we apply a mass fractal model to construct porous media with a priori known mass fractal dimensions \(2.579 \le D_{\mathrm{m}} \le 2.893\) characterizing both solid matrix and pore space. The pore–solid interface of the media is accordingly roughened using the Weierstrass–Mandelbrot approach and two parameters, i.e., surface fractal dimension \(D_{\mathrm{s}}\) and root-mean-square (rms) roughness height. A single-relaxation-time lattice Boltzmann method is applied to simulate single-phase permeability in the corresponding porous media. Results indicate that permeability decreases sharply with increasing \(D_{\mathrm{s}}\) from 1 to 1.1 regardless of \(D_{\mathrm{m}}\) value, while k may slightly increase or decrease, depending on \(D_{\mathrm{m}}\), as \(D_{\mathrm{s}}\) increases from 1.1 to 1.6.     

Estimation of adsorption-induced pore pressure and confinement in a nanoscopic slit pore by a density functional theory

This study aims at characterising the adsorption-induced pore pressure and confinement in nanoscopic pores by molecular non-local density functional theory (DFT). Considering its important potential industrial applications, the adsorption of methane in graphitic slit pores has been selected as the test case. While retaining the accuracy of molecular simulations at pore scale, DFT has a very low computational cost that allows obtaining highly resolved pore pressure maps as a function of both pore width and thermodynamic conditions. The dependency of pore pressure on these parameters (pore width, pressure and temperature) is carefully analysed in order to highlight the effect of each parameter on the confined fluid properties that impact the solid matrix.     

An elastoplasticity model in porous media theories

In this article, porous media theories are referred to as mixture theories extended by the well-known concept of volume fractions. This approach implies the diverse field functions of both the porous solid matrix and the pore fluid to be represented by average functions of the macroscale.The present investigations are based on a binary model of incompressible constituents, solid skeleton, and pore liquid, where, in the constitutive range, use is made of the second-grade character of general heterogeneous media. Within the framework of geometrically finite theories, the paper offers a set of constitutive equations for the solid matrix, the viscous pore liquid and the different interactions between the constituents. The constitutive model applies to saturated as well as to empty solid materials, taking into account the physical nonlinearities based on elasto-plastic solid deformations. In particular, the constitutive model concentrates on granular materials like soil or concrete, where the elastic deformations are usually small and the plastic range is governed by kinematically hardening properties.     

具有分形结构界面的固体膜/粘滞层/基底层状结构的变形

考虑固体膜/粘滞层/基底结构中粘滞层/基底界面不平整对结构的稳定性,特别是固体薄膜稳定时的褶皱变形产生的影响.考虑自仿射和峰状的粗糙界面.自仿射的粗糙界面由振幅和分形维表征,振幅和分形维越大对结构平衡的影响越大;峰状界面由振幅和平均峰间距表征,振幅越大对结构平衡影响越大,平均峰间距越小对结构平衡影响越大.进而,用有限元方法模拟分析该结构,得到界面平整和不平整两种情况下固体膜面内失稳力.结果表明,粘滞层/基底不平整情况下结构的平衡状态与假设该界面为理想平整面所得的结果有很大不同,在结构尺寸较小的情况下不可假设该界面为理想平整来考察该结构的稳定性.     

Investigation on effect of drag models on flow behavior of power-law fluid–solid two-phase flow in fluidized bed

In this study, a Eulerian-Eulerian two-fluid model combined with the kinetic theory of granular flow is adopted to simulate power-law fluid–solid two-phase flow in the fluidized bed. Two new power-law liquid–solid drag models are proposed based on the rheological equation of power-law fluid and pressure drop. One called model A is a modified drag model considering tortuosity of flow channel and ratio of the throat to pore, and the other called model B is a blending drag model combining drag coefficients of high and low particle concentrations. Predictions are compared with experimental data measured by Lali et al., where the computed porosities from model B are closer to the measured data than other models. Furthermore, the predicted pressure drop rises as liquid velocity increases, while it decreases with the increase of particle size. Simulation results indicate that the increases of consistency coefficient and flow behavior index lead to the decrease of drag coefficient, and particle concentration, granular temperature, granular pressure, and granular viscosity go down accordingly.