Observation of fluid layering and reverse motion in double-walled carbon nanotubes

Most modelling-based research in the field of carbon nanotube-related nano-fluidics has been concerned with the fluid flow in single-walled carbon nanotubes (SWCNTs), showing that the dynamics of the channel affect the structure and behaviour of the fluid. We have extended this work by modelling the flow of Ar in a double-walled carbon nanotube, and have modelled the flow in both the inner shell and the outer annular region of such a nanotube. We have found that the flows in these channels are strongly correlated, such that the fluid moves in opposite directions in these two regions. This phenomenon can give rise to a circulatory motion which can be exploited in nano-fluidic devices. Fluid layering phenomenon, that is usually associated with the flow of fluids in nano-scale channels, is also observed. Furthermore, we have also found that the fluid velocity in dynamic channels is smaller than in static channels, in line with the findings reported for single-walled carbon nanotubes.     

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.     

Effect of side walls on the motion of a viscous fluid induced by an infinite plate that applies an oscillating shear stress to the fluid

The velocity field corresponding to the unsteady motion of a viscous fluid between two side walls perpendicular to a plate is determined by means of the Fourier transforms. The motion of the fluid is produced by the plate which after the time t = 0, applies an oscillating shear stress to the fluid. The solutions that have been obtained, presented as a sum of the steady-state and transient solutions satisfy the governing equation and all imposed initial and boundary conditions. In the absence of the side walls they are reduced to the similar solutions corresponding to the motion over an infinite plate. Finally, the influence of the side walls on the fluid motion, the required time to reach the steady-state, as well as the distance between the walls for which the velocity of the fluid in the middle of the channel is unaffected by their presence, are established by means of graphical illustrations.     

Viscous near-wall flow in a wake of circular cylinder at moderate Reynolds numbers

Here we present the results of experimental investigation of a cross flow around a circular cylinder mounted near the wall of a channel with rectangular cross section. The experiments were carried out in the range of Reynolds numbers corresponding to the transition to turbulence in a wake of the cylinder. Flow visualization and SIV-measurements of instantaneous velocity fields were carried out. Evolution of the flow pattern behind the cylinder and formation of the regular vortex structures were analyzed. It is shown that in the case of flow around the cylinder, there is no spiral motion of fluid from the side walls of the channel towards its symmetry plane, typical of the flow around a spanwise rib located on the channel wall. The laminar-turbulent transition in the wake of the cylinder is caused by the shear layer instability.     

Vortex generation in the cross-flow around a cylinder attached to an end-wall

A cylinder attached to an end-wall normal to its axis is a common feature of many practical flow systems, e.g. in turbo-machinery or when a bridge is supported by a pillar from the bed of a river. In this situation, the nominally two-dimensional boundary layer flow incident upon the cylinder develops strong three-dimensional features and a very pronounced vortex structure may arise in the upstream flow close to the wall. For the appropriate Reynolds number range, the upstream vortical structure is nominally steady and is commonly referred to as the “horseshoe vortex system”. In contrast, the flow downstream is unsteady and periodic over a wide range of Reynolds numbers and vortices aligned with the cylinder axis are shed at a regular frequency into the wake. The generation of both these vortex systems requires energy to be extracted from the incident flow with the result that the drag force on the cylinder is increased.This paper concentrates on the upstream region of the cylinder and discusses an investigation in which two-component Particle Image Velocimetry (PIV) has been used to visualise the flow behaviour for a circular cylinder on a plane end-wall. The use of PIV has enabled two orthogonal velocity components to be measured in planes defined by the upstream flow direction and the axis of the cylinder. The third (out-of-plane) velocity component was then calculated by integrating the continuity equation. Subsequently, the velocity field information has been manipulated and converted into time-averaged information.Discussion of the measured results confirms that colour displays are an invaluable aid to understanding this complex fluid flow situation since they reveal substantially more information than grey-scale plots of the same data. In particular, the source of the horseshoe vortex system can be identified when colour plots of the time-averaged velocity and vorticity distributions are obtained. A limited amount of information on the unsteady vortex structures appearing in the end-wall region upstream of the cylinder is also presented. Finally, the experimental findings are discussed in relation to the results of previous workers.     

从非定常效应看对旋风扇/压气机的工程实现

首先利用全隐格式的非稳态压力修正求解方法；计算对旋风扇级间流场非定常流动情况，并与实验结果进行了有益的对比；在动态实验研究中则通过单斜丝热线风速仪对对旋风扇级间速度场进行了详细分析，并揭示其级间非定常效应的独特规律：下游转子的非定常位势流作用在对旋级间流场的周期性非定常影响中占主导地位，级间径向速度的大小则主要与前级转子叶尖区和叶根区的二次流动旋涡强度有关．上述结果表明：对旋风扇／压气机的非定常效应对其性能有很大影响，如果很好地组织对旋叶轮的非定常流动气动布局，可进一步提高增压系统的压比、效率和稳定工作范围。     

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

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

Numerical Study of Unsteady MHD Flow and Entropy Generation in a Rotating Permeable Channel with Slip and Hall Effects

This article investigates an unbiased analysis for the unsteady two-dimensional laminar flow of an incompressible, electrically and thermally conducting fluid across the space separated by two infinite rotating permeable walls.The influence of entropy generation, Hall and slip effects are considered within the flow analysis. The problem is modeled based on valid physical arguments and the unsteady system of dimensionless PDEs (partial differential equations) are solved with the help of Finite Difference Scheme. In the presence of pertinent parameters, the precise movement of the flow in terms of velocity, temperature, entropy generation rate, and Bejan numbers are presented graphically, which are parabolic in nature. Streamline profiles are also presented, which exemplify the accurate movement of the flow. The current study is one of the infrequent contributions to the existing literature as previous studies have not attempted to solve the system of high order non-linear PDEs for the unsteady flow with entropy generation and Hall effects in a permeable rotating channel. It is expected that the current analysis would provide a platform for solving the system of nonlinear PDEs of the other unexplored models that are associated to the two-dimensional unsteady flow in a rotating channel.     

Exact Solutions for the Flow of Fractional Maxwell Fluid in Pipe-Like Domains

This paper presents an analysis of unsteady flow of incompressible fractional Maxwell fluid filled in the annular region between two infinite coaxial circular cylinders. The fluid motion is created by the inner cylinder that applies a longitudinal time-dependent shear stress and the outer cylinder that is moving at a constant velocity. The velocity field and shear stress are determined using the Laplace and finite Hankel transforms. Obtained solutions are presented in terms of the generalized G and R functions. We also obtain the solutions for ordinary Maxwell fluid and Newtonian fluid as special cases of generalized solutions. The influence of different parameters on the velocity field and shear stress is also presented using graphical illustration. Finally, a comparison is drawn between motions of fractional Maxwell fluid, ordinary Maxwell fluid and Newtonian fluid.     

The design secret of Kyokusui-no-En’s meandering channel

The purpose of this study is to investigate the characteristics of the flow through the Jonangu channel which is used for ceremonial game called as ‘Kyokusui-no-En’ in Japanese. The geometry of the channel is measured, a visualization technique is used to measure the actual flow characteristics, and then a numerical flow model is used to represent the flow including unsteady flow characteristics. Numerical model of drifting cup is introduced to investigate an interaction between flow and motion of the cup. Finally, the intention of the channel design is interpreted from the viewpoint of fluid mechanics using observed and calculated results.     

Transmission loss prediction of reactive silencers using 3-D time-domain CFD approach and plane wave decomposition technique

A predictive method is proposed to determine the transmission loss of reactive silencers using the three-dimensional (3-D) time-domain computational fluid dynamics (CFD) approach and the plane wave decomposition technique. Firstly, a steady flow computation is performed with a mass-flow-inlet boundary condition, which provides an initial condition for the following two unsteady flow computations. The first unsteady flow computation is conducted by imposing an impulse (acoustic excitation) superimposed on the constant mass flow at the inlet of the model and then adding the non-reflecting boundary condition (NRBC) when the impulse completely propagates into the silencer. The second unsteady flow computation is conducted for the case without acoustic excitation at the inlet. The time histories of pressure and velocity at the upstream monitoring point as well as history of pressure at the downstream monitoring point are recorded during the two transient computations. The differences between the two unsteady flow computational results are the corresponding acoustic quantities. Therefore, the incident sound pressure signal is obtained by using plane wave decomposition at upstream, while the transmitted sound pressure signal is just the sound pressure at downstream. Finally, those two sound pressure signals in the time-domain are transformed into the frequency-domain by Fast Fourier Transform (FFT) and then the transmission loss (TL) of silencer is determined. For the straight-through perforated tube silencers with and without flow, the numerical results agree well with the published measurements.     

Brownian motion and thermophoresis effects on unsteady stagnation point flow of Eyring-Powell nanofluid: a Galerkin approach

This article concerns the analysis of an unsteady stagnation point flow of Eyring-Powell nanofluid over a stretching sheet. The influence of thermophoresis and Brownian motion is also considered in transport equations. The nonlinear ODE set is obtained from the governing nonlinear equations via suitable transformations. The numerical experiments are performed using the Galerkin scheme. A tabular form comparison analysis of outcomes attained via the Galerkin approach and numerical scheme (RK-4) is available to show the credibility of the Galerkin method. The numerical exploration is carried out for various governing parameters, namely, Brownian motion, steadiness, thermophoresis, stretching ratio, velocity slip, concentration slip, thermal slip, and fluid parameters, and Hartmann, Prandtl and Schmidt numbers. The velocity of fluid enhances with an increase in fluid and magnetic parameters for the case of opposing, but the behavior is reversed for assisting cases. The Brownian motion and thermophoresis parameters cause an increase in temperature for both cases (assisting and opposing). The Brownian motion parameter provides a drop-in concentration while an increase is noticed for the thermophoresis parameter. All the outcomes and the behavior of emerging parameters are illustrated graphically. The comparison analysis and graphical plots endorse the appropriateness of the Galerkin method. It is concluded that said method could be extended to other problems of a complex nature.     

An improved penalty immersed boundary method for fluid–flexible body interaction

An improved penalty immersed boundary (pIB) method has been proposed for simulation of fluid–flexible body interaction problems. In the proposed method, the fluid motion is defined on the Eulerian domain, while the solid motion is described by the Lagrangian variables. To account for the interaction, the flexible body is assumed to be composed of two parts: massive material points and massless material points, which are assumed to be linked closely by a stiff spring with damping. The massive material points are subjected to the elastic force of solid deformation but do not interact with the fluid directly, while the massless material points interact with the fluid by moving with the local fluid velocity. The flow solver and the solid solver are coupled in this framework and are developed separately by different methods. The fractional step method is adopted to solve the incompressible fluid motion on a staggered Cartesian grid, while the finite element method is developed to simulate the solid motion using an unstructured triangular mesh. The interaction force is just the restoring force of the stiff spring with damping, and is spread from the Lagrangian coordinates to the Eulerian grids by a smoothed approximation of the Dirac delta function. In the numerical simulations, we first validate the solid solver by using a vibrating circular ring in vacuum, and a second-order spatial accuracy is observed. Then both two- and three-dimensional simulations of fluid–flexible body interaction are carried out, including a circular disk in a linear shear flow, an elastic circular disk moving through a constricted channel, a spherical capsule in a linear shear flow, and a windsock in a uniform flow. The spatial accuracy is shown to be between first-order and second-order for both the fluid velocities and the solid positions. Comparisons between the numerical results and the theoretical solutions are also presented.     

On the damping of unsteady flow by compliant boundaries

This paper concerns unsteady motion of an incompressible inviscid fluid near a flexible surface which, in responding to the surface pressure field, absorbs energy. The modification of the flow consequent on energy removal at the boundaries is examined. Energy absorption always occurs when the mechanical surface properties include an element of dissipation. But surface dissipation is not essential; surface waves have a similar property. Unsteady fluid induced forces excite surface waves which carry with them energy that must have originated in the flow. The question of how flow characteristics change as energy is gradually given up to the boundary is examined through a particular model problem from which it becomes evident that surface motion draws vorticity towards the surface. The model chosen is that of a rectilinear vortex adjacent to a weakly responding boundary. Surface motion induces a velocity perturbation which is shown to move the vortex towards the surface whenever the fluid gives energy to that surface.     

用谱方法数值模拟槽道内的气固两相流动

在数值模拟领域内,谱方法具有收敛快、分辨率高和精度高的优点．谱方法处理边界方便,随着数值方法的改善和计算机的发展,它在数值模拟中的作用愈加重要．这里采用谱方法数值求解三维Ｎ—Ｓ方程,用这一方法计算了直槽道内流体的流动．计算得到的层流和湍流结果与理论结果符合地较好．在此基础之上进一步模拟了几种不同槽道内的三维粘性流体层流流动,特别在弯曲槽道内的流动计算中,发展了源项处理方法,正确地反映了弯曲固壁对流体流动的影响．通过对湍流计算获得的脉动速度场的统计可以得到湍流运动的许多统计量,正确地反映了湍流运动的特征,说明可以用模拟得到的流场来代替真实的流场．进行了气固两相流动的研究,由直接模拟得到的流体瞬时速度场对固体粒子的作用进行了粒子运动的模拟计算,得到了颗粒在真实流场中运动的浓度,轨道等有用信息和运动特性,得到了令人鼓舞的结果．     

Multi-layer flows of immiscible fractional Maxwell fluids in a cylindrical domain

Laminar unsteady multilayer axial flows of fractional immiscible Maxwell fluids in a circular cylinder are investigated. The flow of fluids is generated by a time-dependent pressure gradient in the axial direction and by the translational motion of a cylinder along his axis. The considered mathematical model is based on the fractional constitutive equation of Maxwell fluids with Caputo time-fractional derivatives. Analytical solutions for the fractional differential equations of the velocity fields with boundary and interfaces conditions have been determined by using the Laplace transform coupled with the Hankel transform of order zero and the Weber transform of order zero. The influence of the memory effects on the motion of the fluid has been investigated for the particular case of three fractional Maxwell fluids. It is found that for increasing values of the fractional parameter the fluid velocity is decreasing. The memory effects have a stronger influence on the velocity of the second layer.     

微扩张管道内幂律流体非定常电渗流动

研究了电渗驱动下幂律流体在有限长微扩张管道内非稳态流动特性.基于Ostwald-de Wael幂律模型,采用高精度紧致差分离散二维Poisson-Nernst-Planck方程及修正的Cauchy动量方程,数值模拟了初始及稳态时刻微扩张管道内幂律流体电渗流流场分布情况,研究了管道截面改变对幂律流体无量纲剪切应变率及无量纲表观黏度的影响,以及无量纲表观黏度对拟塑性流体与胀流型流体流速分布的影响.数值模拟结果显示,当扩张角和无量纲电动宽度一定时,电场驱动下的幂律流体在近壁区域速度响应都很快;初始时刻,近壁处表观黏度的变化受到剪切应变率变化的影响,从而影响了三种幂律流体速度峰值的分布,出现拟塑性流体流速在扩张段上游及扩张段近壁处速度峰值均为幂律流体中最大、而在扩张段下游三种幂律流体速度峰值相近的现象;稳态时刻,幂律流体速度剖面呈现塞型分布,且满足连续性条件下,幂律流体流速随扩张管半径增大而减小,牛顿流体流动规律与宏观尺度下流动规律相同;初始时刻,在相同电动宽度、不同壁面电势作用下,幂律流体在扩张管近壁处剪切应变率分布的差异导致表观黏度分布的差异,并最终导致拟塑性流体与胀流型流体流速分布的差异.     

Analytical and semi-analytical solutions to flows of two immiscible Maxwell fluids between moving plates

Unsteady flows of two immiscible Maxwell fluids in a rectangular channel bounded by two moving parallel plates are studied. The fluid motion is generated by a time-dependent pressure gradient and by the translational motions of the channel walls in their planes. Analytical solutions for velocity and shear stress fields have been obtained by using the Laplace transform coupled with the finite sine-Fourier transform. These analytical solutions are new in the literature and the method developed in this paper can be generalized to unsteady flows of n-layers of immiscible fluids. By using the Laplace transform and classical method for ordinary differential equations, the second form of the Laplace transforms of velocity and shear stress are determined. For the numerical Laplace inversion, two accuracy numerical algorithms, namely the Talbot algorithm and the improved Talbot algorithm are used.     

微平行管道内Jeffrey流体的非定常电渗流动

研究了微平行管道内线性黏弹性流体的非定常电渗流动, 其中线性黏弹性流体的本构关系是由Jeffrey流体模型来描述的. 利用Laplace变换法, 求解了线性化的Poisson-Boltzmann方程、 非定常的柯西动量方程和Jeffrey流体本构方程, 给出了黏弹性Jeffrey流体电渗速度的解析表达式, 分析了无量纲弛豫时间λ₁和滞后时间λ₂对速度剖面的影响. 发现滞后时间为零时, 弛豫时间越小, 速度剖面图越接近牛顿流体的速度剖面图; 随着弛豫时间和滞后时间的增加, 速度振幅也变得越来越大, 随着时间的增加, 速度逐渐趋于恒定. 关键词： 双电层 微平行管道 Jeffrey流体 非定常电渗流动     

Effect of the ion slip on the MHD flow of a dusty fluid with heat transfer under exponential decaying pressure gradient

In the present study, the unsteady Hartmann flow with heat transfer of a dusty viscous incompressible electrically conducting fluid under the influence of an exponentially decreasing pressure gradient is studied without neglecting the ion slip. The parallel plates are assumed to be porous and subjected to a uniform suction from above and injection from below while the fluid is acted upon by an external uniform magnetic field applied perpendicular to the plates. The equations of motion are solved analytically to yield the velocity distributions for both the fluid and dust particles. The energy equations for both the fluid and dust particles including the viscous and Joule dissipation terms, are solved numerically using finite differences to get the temperature distributions.