波纹通道中流动非稳态特性的数值研究

本文对两种不同间距的波纹通道中的流动进行了非定常的数值模拟,给出了间距和雷诺数对反映涡脱落频率的无量纲斯特劳哈尔数的影响情况,比较了二维和三维数值模拟的差别。增大雷诺数或者减小间距将会使斯特劳哈尔数增大,而斯特劳哈尔数随雷诺数增大而增加的趋势会趋于平缓。雷诺数相同的情况下,三维模拟会得到振幅更大且形式更复杂的速度随时间振荡的曲线,而且波纹通道间距越大,速度振荡曲线的振幅也越大且形式越复杂。     

Mixing in a T-shaped micromixer at moderate Reynolds numbers

In the present work, the regimes of the flow and mixing of fluids in a T-shaped micromixer in the range of the Reynolds numbers from 1 to 1000 are investigated systematically with the aid of numerical modeling. The flow and mixing regimes are shown to alter substantially with increasing Reynolds numbers. Five different flow regimes have been identified in the total. The dependencies of the friction coefficient and mixing efficiency on the Reynolds number are obtained. A sharp increase in the mixing efficiency at a flow transition from the symmetric to asymmetric steady regime is shown. On the other hand, the mixing efficiency slightly drops in the laminar-turbulent transition region. A substantial influence of the slip presence on walls on flow structure in the channel and mixing efficiency has been revealed.     

Orientation distribution of fibres in a channel flow of fibre suspension

The orientation and concentration distributions of fibres in laminar and turbulent channel flows are investigated numerically. The obtained results are in good agreement with the experimental data. In the laminar flow regime, more fibres orient to the flow direction as the Reynolds number increases. The shear rate of fluid around a fibre plays an important role in determining the orientation distribution of fibres, while the fibre density and the fibre aspect-ratio have marginal influence on the orientation distribution. In the turbulent regime, the orientation distribution of fibres becomes more homogeneous with the increase of Reynolds number, and the concentration profile is flatter than that in the laminar regime. The fluctuating intensity of fibre velocity in the downstream direction is larger than that in the lateral directions.     

Large-eddy simulations of temporally accelerating turbulent channel flow

The unsteady turbulent channel flow subject to the temporal acceleration is considered in this study. Large-eddy simulations were performed to study the response of the turbulent flow to the temporal acceleration. The simulations were started with the fully developed turbulent channel flow at an initial Reynolds number of Re₀ = 3500 (based on the channel half-height and the bulk-mean velocity), and then a constant temporal acceleration was applied. During the acceleration, the Reynolds number of the channel flow increased linearly from the initial Reynolds number to the final Reynolds number of Re₁ = 22,600. The effect of grid resolution, domain size, time step size on the simulation results was assessed in a preliminary study using simulations of the accelerating turbulent flow as well as simulations of the steady turbulent channel flow at various Reynolds numbers. Simulation parameters were carefully chosen from the preliminary study to ascertain the accuracy of the simulation. From the accelerating turbulent flow simulations, the delays in the response of various flow properties to the temporal acceleration were measured. The distinctive features of the delays responsible for turbulence production, energy redistribution, and radial propagation were identified. Detailed turbulence statistics including the wall shear stress response during the acceleration were examined. The results reveal the changes in the near-wall structures during the acceleration. A self-sustaining mechanism of turbulence is proposed to explain the response of the turbulent flow to the temporal acceleration. Although the overall flow characteristics are similar between the channel and pipe flows, some differences were observed between the two flows.     

Flow separation behind a rib in a channel with laminar flow

Visualization data and results of combined measurements of flow quantities in flow with separation past a rib at nominally laminar regime of channel flow are reported. In the separation region, the flow is found to be essentially three-dimensional and unsteady, exhibiting a distinct cellular structure and flow zones with transverse motion. It is shown that the rib-induced flow separation gives rise to low-frequency fluctuations of flow velocity and initiates the turbulence transition in the channel flow. The critical Reynolds number at which flow instability starts developing in the channel is estimated. It is shown that at Reynolds numbers higher than the critical Reynolds number the linear integral scale of flow velocity fluctuations in the channel is defined by the duct size.     

二维颗粒流通道宽度效应的分子动力学模拟

通过用分子动力学方法对颗粒物质流的计算机模拟，研究发现增大通道宽度可以使二维颗粒流从稀疏流转变为密集流状态.通过对不同通道宽度下，固定开口为9.5d的颗粒流和 漏斗口以上9.5d×8d区域记录的模拟结果分析，发现随通道宽度增大，密度变大、温 度降低.当“颗粒温度”T较低时(T/m<0.05 J/kg),颗粒流内部接触数开始超过1.2 ，同时出现较为牢固的横向链状颗粒团簇，是造成流量突变以及密集流的原因. 关键词： 颗粒物质 颗粒流 计算机模拟     

边界条件对二维斜面颗粒流颗粒分布的影响

实验研究了粗糙边界条件、不同通道宽度以及不同斜面倾斜角下二维颗粒流的颗粒分布 ，发现颗粒流在稀疏流状态下尽管在现象上表现出类似流体的性质，但由于颗粒流系统的能 量耗散,所以颗粒流的速度和密度分布和牛顿流体完全不同，颗粒流在通道中的分布既受通 道宽度的影响同时也受重力场的影响，颗粒在通道中以通道中轴线呈对称分布，通道两侧的 颗粒密度明显高于通道中间，由于通道中颗粒受剪切力的作用，导致颗粒在通道中的密度分 布的变化，通道边界粗糙度明显影响着颗粒流横向分布. 关键词： 二维颗粒流 颗粒物质 稀疏流     

微喷管流的连续介质模型及其适用性

基于无滑移和有滑移的连续介质模型,对微喷管内的超声速冷态气体流场进行了二维和三维数值模拟,利用DSMC方法验证微喷管流中的连续介质模型,并重点分析微喷管流的低雷诺数效应、三维端面效应及其推进性能.研究表明,局部流场的模拟对模型和边界条件的要求要高于推进性能的估算,在努森数小于0.03时,可以使用无滑移的N-S方程预测推进性能;雷诺数是表征低雷诺数效应和推进性能的特征参数,提高工作压力可以改善微喷管的粘性损失和推进性能;在雷诺数大于1000时,若蚀刻深度和喉部宽度的比值超过13,微喷管具备很好的二维特性.     

Hopf Bifurcations,Drops in the Lid-Driven Square Cavity Flow

The lid-driven square cavity flow is investigated by numerical experiments. It is found that from $ \mathrm{Re} $$=$ $5,000 $ to $ \mathrm{Re} $$=$$ 7,307.75 $ the solution is stationary, but at $ \mathrm{Re}$$=$$7,308 $ the solution is time periodic. So the critical Reynolds number for the first Hopf bifurcation localizes between $ \mathrm{Re} $$=$$ 7,307.75 $ and $ \mathrm{Re} $$=$$ 7,308 $. Time periodical behavior begins smoothly, imperceptibly at the bottom left corner at a tiny tertiary vortex; all other vortices stay still, and then it spreads to the three relevant corners of the square cavity so that all small vortices at all levels move periodically. The primary vortex stays still. At $ \mathrm{Re} $$=$$ 13,393.5 $ the solution is time periodic; the long-term integration carried out past $ t_{\infty} $$=$$ 126,562.5 $ and the fluctuations of the kinetic energy look periodic except slight defects. However, at $ \mathrm{Re} $$=$$ 13,393.75 $ the solution is not time periodic anymore: losing unambiguously, abruptly time periodicity, it becomes chaotic. So the critical Reynolds number for the second Hopf bifurcation localizes between $ \mathrm{Re} $$=$$ 13,393.5 $ and $ \mathrm{Re} $$=$$ 13,393.75 $. At high Reynolds numbers $ \mathrm{Re} $$=$$ 20,000 $ until $ \mathrm{Re} $$=$$ 30,000 $ the solution becomes chaotic. The long-term integration is carried out past the long time $ t_{\infty} $$=$$ 150,000 $, expecting the time asymptotic regime of the flow has been reached. The distinctive feature of the flow is then the appearance of drops: tiny portions of fluid produced by splitting of a secondary vortex, becoming loose and then fading away or being absorbed by another secondary vortex promptly. At $ \mathrm{Re} $$=$$ 30,000 $ another phenomenon arises—the abrupt appearance at the bottom left corner of a tiny secondary vortex, not produced by splitting of a secondary vortex.     

Crisis of hydrodynamic drag of drops in the two-phase turbulent flow of a spray produced by a mechanical nozzle at transition reynolds Numbers

When processing experimental data for the hydrodynamics of a two-phase flow in a spray produced by a mechanical nozzle, we revealed an anomaly in the behavior of the hydrodynamic drag of drops: the drag coefficient turns out to be four to seven times lower than the previously known values. Several hypotheses are put forward to explain the anomaly. It is found that, when the gas flows around drops under highly turbulent conditions, an “early” (i.e., observed even at transition Reynolds numbers, Re>50) crisis of drag resistance of drops takes place. This new physical phenomenon allows us to account for a number of features of the two-phase flow that are observed in the experiment. Among these features is, in particular, the fact that the momentum transferred to the gas is roughly half the initial momentum of the liquid jet.     

Surface roughness effects on the turbulence statistics in a low Reynolds number channel flow

A high-resolution particle image velocimetry was used to characterize a low Reynolds number turbulent flow in a channel. Experiments were conducted over a sand grain-coated surface of large relative roughness, and the results were compared with measurements over a smooth surface. The roughness perturbation significantly modified the outer layer. Even though the streamwise Reynolds stress shows less sensitivity in the outer layer to the boundary condition, significant enhancements were observed in the wall-normal Reynolds stress and the Reynolds shear stress. These modifications were considered as footprints of the larger-scale eddies transporting intense wall-normal motions away from the rough wall. A quadrant decomposition shows that strong and more frequent ejections are responsible for the larger values of the mean Reynolds shear stress over the rough wall. The results also indicate that spanwise vortex cores with mean vorticity of the same sign as the mean shear are the dominant smaller-scale vortical structures over the smooth and rough walls. A linear stochastic estimation-based analysis shows that the average larger-scale structure associated with these vortices is a shear layer that strongly connects the outer layer flow to the near-wall flow. A proper orthogonal decomposition of the flow suggests that the large-scale eddy is more energetic for the rough wall, and contributes more significantly to the resolved turbulent kinetic energy and the Reynolds shear stress than the smooth wall.     

The numerical investigation of heat transfer and pressure drop of turbulent flow in a triangular microchannel

In this presentation, the flow and heat transfer inside a microchannel with a triangular section, have been numerically simulated. In this three-dimensional simulation, the flow has been considered turbulent. In order to increase the heat transfer of the channel walls, the semi-truncated and semi-attached ribs have been placed inside the channel and the effect of forms and numbers of ribs has been studied. In this research, the base fluid is Water and the effect of volume fraction of Al₂O₃ nanoparticles on the amount of heat transfer and physics of flow have been investigated. The presented results are including of the distribution of Nusselt number in the channel, friction coefficient and Performance Evaluation Criterion of each different arrangement. The results indicate that, the ribs affect the physics of flow and their influence is absolutely related to Reynolds number of flow. Also, the investigation of the used semi-truncated and semi-attached ribs in Reynolds number indicates that, although heat transfer increases, but more pressure drop arises. Therefore, in this method, in order to improve the heat transfer from the walls of microchannel on the constant heat flux, using the pump is demanded.     

Numerical simulation of laminar forced convection of water-CuO nanofluid inside a triangular duct

In this article, distilled water and CuO particles with volume fraction of 1%, 2% and 4% are numerically studied. The steady state flow regime is considered laminar with Reynolds number of 100, and nano-particles diameters are assumed 20 nm and 80 nm. The hydraulic diameter and the length of equilateral triangular channel are 8 mm and 1000 mm, respectively. The problem is solved for two different boundary conditions; firstly, constant heat flux for all sides as a validation approach; and secondly, constant heat flux for two sides and constant temperature for one side (hot plate). Convective heat transfer coefficient, Nusselt number, pressure loss through the channel, velocity distribution in cross section and temperature distribution on walls are investigated in detail. The fluid flow is supposed to be one-phase flow. It can be observed that nano-fluid leads to a remarkable enhancement on heat transfer coefficient. Furthermore, CuO particles increase pressure loss through the channel and velocity distribution in fully developed cross section of channel, as well. The computations reveal that the size of nano-particles has no significant influence on heat transfer properties. Besides, the study shows a good agreement between provided outcomes and experimental data available in the literature.     

Electrical discharge characteristics of a dielectric liquid under external flow in a microchannel with planar electrode configuration

The effect of fluid flow on discharge characteristics in a dielectric liquid was determined. Tests were performed for flow in a 100 micron high, 5 mm wide, 50 mm long microchannel at Reynolds numbers up to 13 and applied DC potentials up to 1.5 kV through 100 micro planar electrode pairs along the floor of the channel. The current within the conduction regime increased as the flow rate increased for the case where the flow was in the same direction as the applied electric field. In the injection regime, there was an optimum flow rate for maximum current.     

Statistics of the energy dissipation rate and local enstrophy in turbulent channel flow

Using high-resolution direct numerical simulations, the height and Reynolds number dependence of high-order statistics of the energy dissipation rate and local enstrophy are examined in incompressible, fully developed turbulent channel flow. The statistics are studied over a range of wall distances, spanning the viscous sublayer to the channel flow centerline, for friction Reynolds numbers Re_τ=180 and Re_τ=381. The high resolution of the simulations allows dissipation and enstrophy moments up to fourth order to be calculated. These moments show a dependence on wall distance, and Reynolds number effects are observed at the edge of the logarithmic layer. Conditional analyses based on locations of intense rotation are also carried out in order to determine the contribution of vortical structures to the dissipation and enstrophy moments. Our analysis shows that, for the simulation at the larger Reynolds number, small-scale fluctuations of both dissipation and enstrophy show relatively small variations for z⁺?100.     

Kolmogorov scaling of turbulent flow in the vicinity of the wall

How to scale even the simplest of turbulent flows continues to be a cause for considerable controversy. In the present research, a data base compiling results from channel flow direct numerical simulations and turbulent boundary layer experiments is employed to investigate the properties of shear and normal Reynolds stresses very close to the wall. Two types of scaling based on Kolmogorov length and velocity scales are analyzed. It is shown that it is highly likely that large length scales of the order of the channel half-width or the boundary layer thickness play an important role even in the innermost regions of wall-bounded turbulent flows, which hints at the persistence of Reynolds number effects in even high Reynolds number flows.     

Dense granular flow around a penetrating object: experiment and hydrodynamic model

We present in this Letter experimental results on the bidimensional flow field around a cylinder penetrating into dense granular matter, together with drag force measurements. A hydrodynamic model based on extended kinetic theory for dense granular flow reproduces well the flow localization close to the cylinder and the corresponding scalings of the drag force, which is found to not depend on velocity, but linearly on the pressure and on the cylinder diameter and weakly on the grain size. Such a regime is found to be valid at a low enough "granular" Reynolds number.     

Numerical heat transfer analysis in turbulent channel flow over a side-by-side triangular prism pair

Heat transfer and flow behavior in a channel fitted with a transverse triangular prism pair is numerically investigated in the turbulent flow regime for the Reynolds number ranging from 10000 to 50000. The aspect ratio of channel height to the prism base is fixed at 4.0 throughout the study. The Navier-Stokes equation, along with the energy equation, is solved using a finite volume method with the SIMPLE technique and the QUICK numerical scheme for coupling the discretized equations while the standard k-ɛ turbulence model is used for closure of the problem. The numerical result reveals that heat transfer augmentation in the channel with the built-in prism pair can be obtained. It is observed that as compared to a channel, the heat transfer is enhanced by about 17% for a single triangular prism and by some 85% for a triangular prism pair mounted on the channel wall. Effects of the clearance between the prisms on the heat transfer augmentation are presented. The heat transfer enhancement is due to the vortex formation or recirculation zone downstream of the prism elements. However, the presence of the prisms also leads to higher values of friction loss over the channel.     

矩形微通道中流阻特性的实验研究

以工程上常用的66%的乙二醇水溶液作为工质,对几何特性相似而高宽比不同的4种纯铝矩形微通道内的流动特性进行了实验研究,得到了微通道冷板基础性的设计数据。实验测量了Reynolds数在50~500之间的流动阻力系数。实验结果表明:通道高度H与宽度W之比对微通道流阻特性有显著的影响;当Re数小于100时,在实验误差内,流动阻力系数的值近似等于经典理论计算值;随着Re数的增大和高宽比的变化,f的值远大于理论值,这可能是由微通道内部壁面粗糙度效应所导致的。     

A Stochastic Description of Transition Between Granular Flow States

Two-dimensional granular flow in a channel with small exit is studied by molecular dyhamics simulations. We firstly define a key area near the exit, which is considered to be the choke area of the system. Then we observe the time variation of the local packing fraction and flow rate in this area for several fixed inflow rate, and find that these quantities change abruptly when the transition from dilute flow state to dense flow state happens. A relationship between the local flow rate and the local packing fraction in the key area is also given. The relationship is a continuous function under the fixed particle number condition, and has the characteristic that the flow rate has a maximum at a moderate packing fraction and the packing fraction is terminated at a high value with negative slope. By use of the relationship, the properties of the flow states under the fixed inflow rate condition are discussed in detail, and the discontinuities and the complex time variation behavior observed＇in the preexisting works are naturally explained by a stochastic process.