微圆管内液体流动与换热特性的实验研究

以去离子水作为工质,在当量直径300-1570μm的微圆管(1cr18Ni9Ti)内进行层流流动与换热特性的实验研究。实验结果表明当壁面粗糙度超过1％时摩擦系数随着壁面相对粗糙度的增加而增加。在给定外壁面热流密度的条件下,当雷诺数小于100时,管壁轴向导热对换热特性的影响随着壁厚与内径比的增加而增大;当雷诺数大于100时,轴向导热的影响随着雷诺数的增加逐渐减弱,充分发展区的Nu数接近常规值4．36。     

热对流条件下颗粒沉降的格子Boltzmann方法模拟

采用格子Boltzmann方法模拟了在热对流条件下的颗粒沉降问题, 在研究单颗粒在等温流体、热流体和冷流体中运动的基础上, 进一步模拟了两个不同温度的颗粒在流体中的沉降.结果表明:两等温颗粒的沉降方式与雷诺数Re以及格拉晓夫数Gr密切相关, 而两不同温度的颗粒与两等温颗粒的沉降规律有显著不同.无论初始位置如何, 冷颗粒最终总位于热颗粒下方运动, Re较大时, 发生连续的拖曳、接触现象, 而Re较小时, 冷颗粒会以较大的沉降速度远离热颗粒. 关键词： 格子Boltzmann方法 颗粒沉降 热对流     

通道中串列旋转圆柱的格子Boltzmann-虚拟区域方法的模拟

通过结合格子Boltzmann方法(LBM)和虚拟区域(Fictitiou sDomain)思想,建立格子Boltzmann-虚拟区域(LB-DF/FD)方法.采用两套网格系统,欧拉网格用于流体,拉格朗日网格用于固体.原有的LBM在计算运动固体的受力方面存在数据振荡,LB-DF/FD方法改进了此缺陷.为验证该方法,模拟圆柱绕流、圆形颗粒在无限长通道中平动及在无限大流场中转动三种情况,结果与其他数值解及理论解符合得很好.利用该方法模拟低雷诺数下通道中串列旋转圆柱周围的流场,分析圆柱间距(g)及雷诺数(Re)对流场结构的影响.给出Re=0.001,0.1和10下,0.2≤g≤8.0的流线结构、圆柱升力、阻力以及力矩等数值结果.结果表明,g对流场的结构及圆柱的受力有显著影响,Re对圆柱阻力及Stokes单元数目的影响较大.     

方形颗粒在黏弹性流体中的沉降特性

采用有限元任意拉格朗日-欧拉(ALE)法对方形颗粒在黏弹性流体中的沉降特性进行研究。通过直接数值模拟得到了不同弹性数下方形颗粒的稳定取向角的变化情况,并讨论了颗粒长宽比和通道宽度对其沉降特性的影响。结果表明,当方形颗粒在黏弹性流体中沉降时,弹性数存在一个临界值。当弹性数小于临界值时,颗粒的稳定取向为长轴方向垂直于重力方向;当弹性数大于临界值时,颗粒的稳定取向为长轴方向平行于重力方向。颗粒长宽比和通道宽度对其沉降特性都有一定的影响。长宽比大的颗粒在沉降过程中的取向角和横向漂移的振幅更大。弹性数的临界值随着长宽比的增大而减小,随着阻塞比的增大而增大。     

颗粒滚动摩擦系数对堆积特性的影响

为探讨颗粒摩擦系数对堆积特性的影响,利用离散单元法,模拟静摩擦系数固定时,变化滚动摩擦系数对椭球形颗粒堆积角及堆积体的影响.利用自制斜面仪测定了颗粒静摩擦系数,并对滚动摩擦系数与堆积角建立线性方程.结果表明,滚动摩擦系数对颗粒堆积特性有显著影响,颗粒堆积角随滚动摩擦系数增大而增大,边界圆与连续圆间的间隙随滚动摩擦系数增大而减小.依据颗粒堆积过程中旋转动能的变化可以阐述上述结果,建立的滚动摩擦系数与堆积角的线性方程可为具体颗粒物料滚动摩擦系数的获取提供数值测量思路.模拟堆积的过程可为散体物料一些不易测量的参数进行虚拟实验标定.     

微细多孔介质中流动及换热实验研究

本文对空气流过烧结微细多孔介质的流动和对流换热进行了实验研究。分析了不同颗粒直径条件下摩擦阻力系数与等效雷诺数的关系,并通过实验数据拟和得到了渗透率K和惯性常数F。结果表明:当颗粒直径比较大时,摩擦阻力系数的实验结果与计算关联式基本吻合;随着颗粒直径的减小,二者的差别增大。通过实验得到了微细烧结多孔介质内部体积平均对流换热系数,并与已有关联式的计算结果基本吻合。     

颗粒沉降运动的虚拟区域法直接数值模拟

文采用基于四边形网格的分布式拉格朗日乘子／虚拟区域方法(DLM／FD method)对二维方槽内775个圆形颗粒在流体中的沉降过程进行了直接数值模拟。得到了颗粒流沉降过程中流体和颗粒速度和涡量分布、流场压力分布等流动细节,展示了颗粒在沉降过程中由于相间的相互作用以及颗粒间的作用,使得颗粒流在流场内形成大小不一的旋流区,颗粒回旋着沉降,同时颗粒的尾涡影响附近颗粒的运动．本文的结果说明分布式拉格朗日乘子／虚拟区域方法对模拟存在很多颗粒的悬浮体流动是可行的。     

双分散颗粒体系在临界堵塞态的结构特征

堵塞行为是颗粒体系中一种常见的现象,其力学性质与堆积结构的关联非常复杂.本文采用离散元法研究了由两种不同半径颗粒组成的二维双分散无摩擦球形颗粒体系在临界堵塞态所呈现的结构特征,讨论了大小颗粒粒径比与大颗粒百分比对临界堵塞态的影响.数值模拟结果表明,当粒径比小于1.4时,临界平均接触数与大颗粒百分比关系不大,当粒径比大于1.4时随着大颗粒百分比的增大临界平均接触数先减小再增大.而临界体积分数在粒径比小于1.8时随着大颗粒百分比的增加先减小后增大,大于1.8时又基本不随大颗粒百分比而变化.大颗粒百分比在接近0或1时,系统近似为单分散体系,临界平均接触数与体积分数基本不随半径比的增大而变化;在接近0.5时,临界平均接触数随着半径比的增大逐渐减小,而临界体积分数则是先减小后增大.文中对大-小颗粒这一接触类型的百分比也进行了探讨,其值随着大颗粒百分比的增大呈二次函数的变化趋势,粒径比对这一变化趋势只有较小的影响.     

低阶格式在大涡模拟计算中的适用性

采用三维Taylor-Green涡作为研究对象,利用工程中常用的低阶数值格式,研究格式本身的数值误差对大涡模拟计算的影响.结果表明:三种数值格式的数值耗散行为都与亚格子模型行为类似,即在小雷诺数下,流场比较光滑时,耗散很小,当雷诺数增加,流动转捩为湍流,流场梯度增大,耗散显著增大.对于MUSCL格式和二阶有界中心格式,在高雷诺数下,亚格子尺度模型没有明显改善计算结果,但也没有使计算结果恶化.中心格式相比其它两种格式,数值耗散最小,但是在高雷诺数湍流情况下,中心格式的数值耗散仍然主导了能量的耗散,再添加亚格子模型,计算结果反而变得稍差.对于工程中的低阶格式而言,采用中心格式计算大涡模拟是比较好的选择,而且在计算不存在稳定性问题时,采用不添加亚格子模型的隐式大涡模拟效果更好.     

颗粒趋壁沉积的直接数值模拟

气相采用直接数值模拟方法、颗粒相采用拉格朗日轨道模型方法对摩擦雷诺数为180的充分发展槽道内Stokes数为1～104的颗粒趋壁沉积现象进行了研究.研究表明,惯性颗粒在湍流涡的作用下趋向壁面运动,并以一定速率沉积.在"吸收壁面"条件下沉积速率首先随着颗粒弛豫时间增大而增大,在弛豫时间大于100时,沉积速率保持不变.     

高超临界雷诺数区间内二维圆柱绕流的实测研究

实测强风工况下高度167 m的徐州彭城电厂冷却塔的表面风荷载,并归纳历史上其他研究人员给出的实测结果,以丰富高超临界雷诺数(Re)区间二维圆柱绕流的试验成果.在低湍流度均匀流场和高湍流度大气边界层流场中分别开展4种风速8类粗糙度条件下的冷却塔刚性模型测压风洞试验,通过对比低雷诺数(Re=2.1×10~5—4.19×10~5)条件下的风洞试验结果和高雷诺数(Re=5.4×10~7—1×10~8)条件下的现场实测结果研究各种静动态绕流特征随雷诺数的变化规律,重点考察雷诺数无关现象的产生条件.研究结果表明,对于物表相对粗糙度在0.01以上的圆柱绕流,雷诺数不相关现象存在于很宽的雷诺数范围(2×10~5Re1×10~8)内;增大来流湍流度亦能引起的雷诺数无关现象,但此时该现象可能仅存在于一个较窄的低雷诺数范围内.     

基于间断有限元方法的并列圆柱层流流动特性

基于高阶的间断有限元方法, 数值模拟低马赫数下并列圆柱的可压缩层流流动, 捕捉并列圆柱流场中的漩涡结构, 以便分析并列圆柱尾流的流动特性. 针对二维圆柱的边界形式, 采用曲边三角形单元构造二维圆柱的曲面边界, 以适应高阶离散格式的精度. 在验证方法合理性的基础上, 分析圆柱间距及雷诺数对漩涡脱落及受力特性的影响规律. 研究结果表明: 并列圆柱的间距是影响流场流动特性的一个主要因素, 它会改变圆柱漩涡脱落的形式. 随着圆柱间距的增加, 上下圆柱的平均阻力系数及平均升力系数的绝对值随之显著下降. 雷诺数对于平均阻力系数的影响相对较小. 但随着雷诺数的增加, 上下圆柱的平均升力系数会随之降低, 而漩涡的脱落频率会随之增大.     

The Oseen drag on a circular cylinder revisited

A scheme is presented to evaluate the drag on an infinite circular cylinder, at rest in a perpendicular uniform stationary flow. This scheme leads already in second approximation to very satisfactory values for the drag coefficient up to Reynolds number R = 10 and probably even up to R = 10².     

Predetermined control of turbulent boundary layer with a piezoelectric oscillator

With a piezoelectric(PZT) oscillator, the predetermined controls of the turbulent boundary layer(TBL) are effective in reducing the drag force. The stream-wise velocities in the TBL are accurately measured downstream of the oscillator driven by an adjustable power source. The mean velocity profiles in the inner and outer scales are reported and the skin friction stresses with different voltage parameters are compared. Reduction of integral spatial scales in the inner region below y+of 30 suggests that the oscillator at work breaks up the near-wall stream-wise vortices responsible for high skin friction. For the TBL at Reθof 2183, the controls with a frequency of 160 Hz are superior among our experiments and a relative drag reduction rate of 26.83% is exciting. Wavelet analyses provide a reason why the controls with this special frequency perform best.     

基于浸入式边界方法的串联双矩形柱绕流数值模拟

基于浸入式边界算法(Virtual Boundary Method)中力源反馈边界的思想,改进其原有内部流体处理方法以减少计算耗费,并结合非等间距网格以便工程应用计算,模拟雷诺数范围内(Re=200–10³)串联双矩形柱绕流,研究表明:Re=200–300时,前柱尾流涡脱处于双剪切层控制阶段;柱间涡街为Karman类涡街,在小间距条件下被抑制,形成涡环;前柱对后柱屏蔽效应体现为后柱阻力系数远小于前柱;临界间距时柱间涡街充分发展,后柱阻力系数等气动参数亦在此发生跃升,但仍小于前柱值;随雷诺数升高,尾流涡街尺寸缩小,临界间距及跃升幅度变小. Re=400时,前柱尾流涡脱进入冲击剪切层控制阶段,阻力系数不再呈现规律性振荡;此后随雷诺数升高,冲击剪切层逐步完善,前柱流动分离使其表面产生更多附着涡,导致尾流旋涡尺寸进一步减小,屏蔽效应消失,涡脱更为剧烈,进而对后柱产生脉动冲击效应;适当间距比条件下此类脉动冲击效应使得后柱阻力系数发生跃升,并略高于前柱. 关键词： 浸入式边界算法 串联双矩形柱 屏蔽效应 临界间距     

同位旋对用裂变截面探测鞍点前摩擦的影响

耦合了轻粒子发射的随机Langevin模型被用于计算三个Bi裂变系统因核耗散效应引起的相对于标准统计模型预计的裂变截面值的变化,σ_fdrop作为激发能的函数。发现209Bi的σ_fdrop要大于202Bi和195Bi的σ_fdrop,这表明高同位旋系统209Bi的裂变截面对摩擦强度有更大的敏感性。该结果建议,当实验上使用裂变激发函数来更精确地提取鞍点前摩擦强度时,应产生具有高同位旋的复合核。进而,通过Langevin计算拟合p+206Pb和p+209Bi反应的裂变激发函数数据,提取的鞍点前摩擦强度为（3~5）×1021 s-1。Langevin models that are coupled with particle emission are applied to calculate the drop of fission cross section caused by friction over its standard statistical-model value, σ_fdrop, as a function of excitation energy for three Bi nuclei with different isospins. We find that σ_fdrop become larger for 209Bi than for 202Bi and 195Bi, indicating a greater sensitivity of fission cross section to friction for high-isospin 209Bi. These results suggest that, to more accurately determine the presaddle dissipation strength by measuring fission excitation functions, it is best to populate those compound systems with a high isospin. Furthermore, we analyze the data of fission excitation functions measured in p + 209Bi and p + 208Pb reactions, which yield high-isospin 210Po and 209Bi systems, and extract a presaddle friction strength of (3-5)×1021 s-1 by comparing Langevin simulations with these experimental data.     

Response of force behaviors of a spherical particle to an oscillating flow

Combined multi-direct forcing technique and the immersed boundary method is applied to investigate the response of force behaviors of a fixed spherical particle to an oscillating flow. The influences of the inlet oscillating flow at a mean Reynolds number of 300 with six different oscillating frequencies and three oscillating amplitudes are investigated. Three different zones with different behaviors are identified and the specific behaviors of the drag and the lift coefficients are analyzed. The averaged drag coefficient and the maximum lift force exerted on the particle increase with the increment of the oscillating Reynolds number or the oscillating amplitude. When the frequency of the inlet flow is equal to the natural vortex shedding frequency, the average drag coefficient reaches the maximum value. A linear relation between the gradient of the maximum lift coefficient and the frequency as well as the amplitude of the inlet oscillating flow is observed in certain regime.     

A study on drag reduction of a rotationally oscillating circular cylinder at low Reynolds number

The flow field around a rotationally oscillating circular cylinder in a uniform flow is studied by using a particle image velocimetry to understand the mechanism of drag reduction and the corresponding suppression of vortex shedding in the cylinder wake at low Reynolds number. Experiments are conducted on the flow around the circular cylinder under rotational oscillation at forcing Strouhal number 1, rotational amplitude 2 and Reynolds number 2,000. It is found from the flow measurement by PIV that the width of the wake is narrowed and the velocity fluctuations are reduced by the rotational oscillation of the cylinder, which results in the drag reduction rate of 30%. The mechanism of drag reduction is studied by phase-averaged PIV measurement, which indicates the formation of periodic small-scale vortices from both sides of the cylinder. It is found from the cross-correlation measurement between the velocity fluctuations that the large-scale structure of vortex shedding is almost removed in the cylinder wake, when the small-scale vortices are generated at the unstable frequency of shear layer by the influence of rotational oscillation.     

Numerical simulation of flow around a circular cylinder with curved sectional grooves

In a circular cylinder with uniform flow, a sudden decrease in the drag force occurs at a high Reynolds numbers; however, it is known that the same phenomenon occurs at a lower Reynolds number in the case where there exist grooves or roughness on the circular cylinder surface. To clarify the flow characteristics around a circular cylinder in the case of changing the shape of grooves we analyzed the drag coefficient, lift coefficient, turbulent kinetic energy, vorticity and pressure by applying the RNGk-ε turbulent model. The shapes of the grooves were arced, triangulated and curved. The results showed that the separation point for a circular cylinder with curved sectional grooves shifts to the most downstream side and the drag coefficient becomes the smallest among circular cylinders with grooves.