纳米尺度下流体剪切流动的分子动力学模拟

利用分子动力学方法模拟流体在两无限大平板间剪切流动过程。研究通道内加入不同体积分数纳米颗粒、体积分数相同纳米颗粒数目不同以及剪切速度对流体密度、速度以及界面滑移的影响。结果表明:近壁区流体数密度呈衰减振荡分布,由近壁区到主流区振幅逐渐减小,颗粒和流体的整体数密度在中心主流区呈抛物线分布。流体剪切应变率随颗粒体积分数的增加逐渐减小,同时剪切应变率和滑移速度之间呈近似线性分布。体积分数相同颗粒数目不同,颗粒在运动过程中呈线,性排列时,剪切应变率最大。随剪切速度的增加,流体滑移速度和滑移长度随之增大,但滑移长度增加量相对较小。     

矩形微通道中气–非牛顿流体两相流实验研究

本文通过实验研究了空气与非牛顿流体在小尺寸矩形微通道中的两相流动特性,给出了在不同气液两相流速下的流型图。同时比较了气液两相流速,液相黏度,表面张力,微通道尺寸对各流型分布区域以及Taylor气泡/液柱长度的影响。发现气液两相流流速变化对流型区域的影响极为明显,而黏性力和表面张力只是在局部范围改变了流型分布。Taylor气泡长度随气液两相流速比增大而增大,随液相黏度及表面张力增大而减小。液柱长度随液相黏度增大而增大,随两相流速比及表面张力增大而减小。最后,我们给出了基于无量纲数JG/JL,Re和Ca的Taylor气泡/液柱长度预测公式,公式预测结果与实验结果基本一致。     

圆管等壁温振荡流动对流换热理论研究

对等壁温边界条件下圆管内层流振荡流动对流换热问题进行理论推导。在随时间做正弦变化的速度入口条件下,经求解圆管内流体振荡时的能量方程,得到振荡流动圆管内温度场分布的解析解表达式,并利用Matlab对表达式进行直观表示。计算结果表明:等壁温边界条件下速度入口脉冲流动可引起温度随时间波动,波动随无量纲振荡频率和速度振荡幅度增加,换热效果加强,且速度振荡幅度和无量纲振荡频率存在协同作用。     

纳米通道内气体剪切流动的分子动力学模拟

采用分子动力学模拟方法研究了表面力场对纳米通道内气体剪切流动的影响规律.结果显示通道内的气体流动分为两个区域:受壁面力场影响的近壁区域和不受壁面力场影响的主流区域.近壁区域内,气体流动特性和气体动力学理论预测差别很大,密度和速度急剧增大并出现峰值,正应力变化剧烈且各向异性,剪切应力在距壁面一个分子直径处出现突变.主流区域的气体流动特性与气体动力学理论预测相符合,该区域内的密度、正应力与剪切应力均为恒定值,速度分布亦符合应力-应变的线性响应关系.不同通道高度及密度下,近壁区域的归一化密度、速度及应力分布一致,表明近壁区域的气体流动特性仅由壁面力场所决定.随着壁面对气体分子势能作用的增强,气体分子在近壁区域的密度和速度随之增大,直至形成吸附层,导致速度滑移消失.通过剪切应力与切向动量适应系数(TMAC)的关系,得到不同壁面势能作用下的TMAC值,结果表明壁面对气体分子的势能作用越强,气体分子越容易在壁面发生漫反射.     

管内变密度流体流动与换热的充分发展及解析解

本文提出了圆管内变密度流体流动与换热充分发展的概念和相关假设及其推论.在所提出的概念和假设基础上对理想气体在圆管内的流动与换热进行了解析求解,得到了充分发展时的无量纲速度抛物分布及无量纲温度分布.对圆管内密度随温度变化的氩气的流动与换热数值模拟所获得的无量纲速度分布及无量纲温度分布在入口段后与上述解析解非常吻合,从而验证了在本文的计算条件范围内,圆管内变密度流体充分发展流动与换热的概念及其解析解是合理的.     

基于拟沸腾理论的超临界CO2管内传热恶化量纲分析

超临界流体广泛应用于工程技术领域,其流动传热特性对工程设计具有重要意义,但是,由于超临界流体的物理微观和宏观行为的机理尚不清晰,所以其异常的流动传热特性并未得到很好的解决.普遍认为超临界流体在分子尺度上可分为类气和类液两种不同的特性,直到最近通过实验在宏观上监测到超临界水类液和类气之间的转变,且这一过程与拟沸腾理论一致,使得问题逐渐变得清晰.本文基于拟沸腾理论对超临界CO2异常流动传热行为进行了研究,在假设类液和类气转换过程不均匀的情况下,从经典的量纲分析和亚临界过冷沸腾理论模型出发,提出了一个适用于超临界流体拟沸腾换热过程的分析方法.通过引入表征类气膜生长速度与流体主流平均流速之比π=(qw·ρ1)/(G·Δi·ρg)和表征近壁区类气膜温度梯度π13=(qw·βpc·di)/λg两个无量纲数,来表征拟沸腾如何导致传热恶化,解释了超临界CO2竖直向上加热流动过程中的异常换热特性,即较大的类气膜生长速度使近壁区快速聚集了较多的高温流体,而较大的类气膜温度梯度使类气膜覆盖在壁面.当核心的冷类液不能充分润湿热壁面时,传热恶化.新无量纲数较好的诠释了超临界流体拟沸腾诱导传热恶化机制,为超临界拟沸腾传热研究提供了理论依据.     

大曲率弯道中拟塑性流体的湍流数值计算研究

本文采用一种耦合流体流变性的修正的低雷诺数k-ε模型,结合前期建立的适用于拟塑性非牛顿流体的高精度算法,数值模拟研究了90°大曲率方截面弯道中拟塑性流体的湍流流动,并与实验数据进行比较。结果表明,主流速度的预测值与实验结果呈现出较为一致的分布趋势,表明了方法与模型的有效性。在弯道及其下游平直恢复段,拟塑性流体湍流的脉动动能最大值均出现在弯道的侧壁附近。而在弯道的下游恢复段,轴向和径向湍流脉动强度在从截面中心到弯管外壁的方向上分布呈现出先减小后增大的趋势。     

非稳态流动对管内湍流换热的增强机理

通过对非稳态管内湍流流动及换热的数值模拟,发现决定非稳态信号强弱的因素对换热效果的影响不同。频率因素的影响非常有限,而振幅因素的影响在无量纲幅值小于1的情况下也不明显,但是当无量纲幅值大于1时,管内流动的换热强化非常显著,同时管内流动的阻力系数也有显著增加。由于幅值大于1带来的流动反向导致管道内的流动剧烈振荡,靠近壁面处的速度梯度在短期内发生大幅度改变并降低壁温,从而导致换热效果的显著增强。     

����Һ̬������ṹ������ģ��

应用SolidWorks和ANSYS软件设计了自由液态锡表面流体结构模型,计算了初始流速和热通量不同时液态锡的速度变化和温度变化,得到了流动液态锡的速度分布和温度分布。结果表明,垂直流动方向液态锡流速较为均匀,沿流动方向液态锡流速逐渐增大、液态锡液面厚度逐渐变薄。初始温度为600K的条件下,热通量为1MW·m?2时,液态锡出口温度为623.38K;热通量为5MW·m?2时,液态锡出口温度为720.18K。在相同条件下使用液态锂作为计算流体,结果表明出口处液态锂的温度低于液态锡的温度。     

自由液态锡表面结构的流体模型

应用SolidWorks 和ANSYS 软件设计了自由液态锡表面流体结构模型,计算了初始流速和热通量不同时液态锡的速度变化和温度变化,得到了流动液态锡的速度分布和温度分布。结果表明,垂直流动方向液态锡流速较为均匀,沿流动方向液态锡流速逐渐增大、液态锡液面厚度逐渐变薄。初始温度为600K 的条件下,热通量为1MW·m⁻² 时,液态锡出口温度为623.38K;热通量为5MW·m⁻² 时,液态锡出口温度为720.18K。在相同条件下使用液态锂作为计算流体,结果表明出口处液态锂的温度低于液态锡的温度。     

Direct numerical simulation of the turbulent flows of power-law fluids in a circular pipe

Fully developed turbulent pipe flows of power-law fluids are studied by means of direct numerical simulation. Two series of calculations at generalised Reynolds numbers of approximately 10000 and 20000 were carried out. Five different power law indexes n from 0.4 to 1 were considered. The distributions of components of Reynolds stress tensor, averaged viscosity, viscosity fluctuations, and measures of turbulent anisotropy are presented. The friction coefficient predicted by the simulations is in a good agreement with the correlation obtained from experiment. Flows of power-law fluids exhibit stronger anisotropy of the Reynolds stress tensor compared with the flow of Newtonian fluid. The turbulence anisotropy becomes more significant with the decreasing flow index n. An increase in apparent viscosity away from the wall leads to the damping of the wall-normal velocity pulsations. The suppression of the turbulent energy redistribution between the Reynolds stress tensor components observed in the simulations leads to a strong domination of the axial velocity pulsations. The damping of wall-normal velocity pulsations leads to a reduction of the fluctuating transport of momentum from the core toward the wall, which explains the effect of drag reduction.     

Tracer dispersion in power law fluids flow through porous media: Evidence of a cross-over from a logarithmic to a power law behaviour

An analytical model is presented to describe the dispersion of tracers in a power-law fluid flowing through a statistically homogeneous and isotropic porous medium. The model is an extension of Saffman's approach to the case of non-Newtonian fluids. It is shown that an effective viscosity depending on the pressure gradient and on the characteristics of the fluid, must be introduced to satisfy Darcy's law. An analytical expression of the longitudinal dispersivity is given as a function of the Peclet number Pe and of the power-law index n that characterizes the dependence of the viscosity on the shear rate . As the flow velocity increases the dispersivity obeys an asymptotic power law: . This asymptotic regime is achieved at moderate Peclet numbers with strongly non-Newtonian fluids and on the contrary at very large values when n goes to 1 ( for n=0.9). This reflects the cross-over from a scaling behaviour for towards a logarithmic behaviour predicted for Newtonian fluids (n=1). Received: 22 July 1997 / Revised and Accepted: 2 July 1998     

Shear-thinning-induced chaos in taylor-couette flow

The effect of weak shear thinning on the stability of the Taylor-Couette flow is explored for a Carreau-Bird fluid in the narrow-gap limit. The Galerkin projection method is used to derive a low-order dynamical system from the conservation of mass and momentum equations. In comparison with the Newtonian system, the present equations include additional nonlinear coupling in the velocity components through the viscosity. It is found that the critical Taylor number, corresponding to the loss of stability of the base (Couette) flow, becomes lower as the shear-thinning effect increases. That is, shear thinning tends to precipitate the onset of Taylor vortex flow. Similar to Newtonian fluids, there is an exchange of stability between the Couette and Taylor vortex flows, which coincides with the onset of a supercritical bifurcation. However, unlike the Newtonian model, the Taylor vortex cellular structure loses its stability in turn as the Taylor number reaches a critical value. At this point, a Hopf bifurcation emerges, which exists only for shear-thinning fluids.     

Shear viscosity for inelastic hard spheres

The hydrodynamic equations of the Enskog theory for inelastic hard spheres is considered as a model for rapid flow granular fluids at finite densities. A detailed analysis of the shear viscosity of the granular fluid has been done using homogenous cooling state (HCS) and uniform shear flow (USF) models. It is found that shear viscosity is sensitive to the coefficient of restitution α and pair correlation function at contact. The collisional part of the Newtonian shear viscosity is found to be dominant than its kinetic part.     

分形油藏非牛顿幂律流体低速非达西不稳定渗流的组合数学模型

结合分形理论与渗流理论,对分形油藏非牛顿幂律流体低速非达西不稳定渗流的试井分析问题的数学模型进行了推导.该分形油藏模型由内域为非牛顿幂律流体低速非达西渗流,外域为非牛顿幂律流体达西渗流的同心圆域组成.在考虑井筒储存、表皮效应影响下,建立了该油藏的不稳定渗流有效井径组合数学模型,在3种外边界条件下求出了两个区域内压力在Laplace空间的解析解,应用Stehfest数值反演方法求得井底的无因次压力,分析了井底压力动态特征和参数影响.非牛顿幂律流体的幂律指数、分形参数均对典型曲线产生较大的影响,呈现出与牛顿流体和均质油藏明显不同的特征.这对非均质油藏非牛顿流体的不稳定试井分析及研究其非线性渗流特征均十分重要.     

剪切变稀液滴撞击不同浸润性壁面的数值模拟研究

基于有限元法,采用水平集方法捕捉相界面的移动,构建了液滴撞击固体壁面的数值模型.通过修正的幂律模型描述流体的非牛顿剪切变稀特性,探讨了剪切变稀特性对液滴撞击固体壁面后铺展行为的影响,分析了撞击不同浸润性壁面时剪切变稀特性对液滴撞击壁面行为的影响差异.研究结果表明:随着幂律指数m的减小,液滴撞击过程中的黏性耗散减小,液滴的形貌变化及无量纲参数变化更为显著.接触角为55°的情况下:当m降低至0.85时,液滴铺展过程中开始出现显著区别于牛顿流体液滴的振荡现象;当m降低至0.80时,液滴在回缩过程中会出现中心液膜断裂的情况.接触角为100°时,剪切变稀液滴均会出现振荡行为,振荡幅度随着m的减小而增大.接触角为160°时,牛顿流体液滴与剪切变稀液滴均会在回缩过程中弹起,但剪切变稀液滴的弹起速度更快.此外,基于数值计算结果,本文提出了接触角为55°情况下剪切变稀液滴撞击壁面后的最大无量纲铺展直径预测模型.     

Drag and Lift Forces Acting on a Sphere in Shear Flow of Power-Law Fluid

Laminar flow of a power-law fluid over a sphere is considered for unbounded shear flow. The Navier–Stokes equations with power-law viscosity are solved numerically using an in-house developed CFD package. Vorticities structures downstream of particle are suppressed for powerlaw fluid. The shear rate influence on drag force is negligible for power index close to unit, and the drag force appreciably decreases with falling power index. For small Reynolds numbers, the lift force coefficient monotonically decreases against the power index and exhibits an opposite behavior for moderate values of Reynolds numbers. The results of the parametric studies are used to derive correlations for the drag force and to detect the hydrodynamic differences from uniform flow. The investigation parameters varied within the following ranges: power-law index 0.3 ≤ n ≤ 1, Reynolds number 0 < Re ≤ 150, and dimensionless shear rate 0.05 ≤ s ≤ 0.4.     

刮膜蒸发器内非牛顿流体流场特性数值模拟

刮膜蒸发器是通过旋转刮板强制成膜,可实现高黏度非牛顿流体类物料平稳蒸发的新型高效蒸发器.蒸发器内流体的流动、分布与传输机制直接决定了蒸发器的蒸发效率与功耗.不同于现有研究主要基于牛顿流体开展,本文针对不同黏度的非牛顿流体,建立蒸发器三维计算流体动力学模型,系统研究了蒸发器内的流场分布特性和成膜机理.结果表明:低黏非牛顿流体的流场分布特性和牛顿流体类似,物料可在壁面形成均匀且连续的液膜;随着黏度的增加,液膜的均匀性和连续性逐渐变差.通过对流场分布与传输形式的研究,结合液膜分布、速度分布、剪应变率分布,以及黏度分布进行对比分析发现,蒸发器内部结构与运行状态形成的剪切场与黏度分布是蒸发器良好成膜的关键.此外,提出对刮板前缘进行弯折可辅助高黏流体液膜铺展,并对最佳弯折角度进行探索.本研究为刮膜蒸发器的设计和应用提供了理论指导与依据.     

Measurement of the viscosity-density product using multiple reflections of ultrasonic shear horizontal waves

We have developed an on-line computer-controlled sensor, based on ultrasound reflection measurements, to determine the product of the viscosity and density of a liquid or slurry for Newtonian fluids and the shear impedance of the liquid for non-Newtonian fluids. A 14 MHz shear wave transducer is bonded to one side of a 45-90 degrees fused silica wedge and the base is in contract with the liquid. Twenty-eight echoes were observed due to the multiple reflections of an ultrasonic shear horizontal (SH) wave within the wedge. The fast Fourier transform of each echo was obtained for a liquid and for water, which serves as the calibration fluid, and the reflection coefficient at the solid-liquid interface was obtained. Data were obtained for 11 sugar water solutions ranging in concentration from 10% to 66% by weight. The viscosity values are shown to be in good agreement with those obtained independently using a laboratory viscometer. The data acquisition time is 14s and this can be reduced by judicious selection of the echoes for determining the reflection coefficient. The measurement of the density results in a determination of the viscosity for Newtonian fluids or the shear wave velocity for non-Newtonian fluids. The sensor can be deployed for process control in a pipeline, with the base of the wedge as part of the pipeline wall, or immersed in a tank.     

Shear Stress in a Couette Flow of Liquid-Particle Suspensions

The mechanisms of momentum transfer and shear stress of liquid-particle suspensions in two-dimensional Couette flow are studied using direct numerical simulation by lattice-Boltzmann techniques. The results obtained display complex flow phenomena that arise from the two-phase nature of the fluid including a nonlinear velocity profile, layering of particles, and apparent slip near the solid walls. The general rheological behaviour of the suspension is dilatant. A detailed study of the various momentum transfer mechanisms that contribute to the total shear stress indicates that the observed shear thickening is related to enhanced relative solid phase stress for increasing shear rates.