Transient flow around an impulsively started cylinder using a dynamic mesh method

The finite volume method (FVM) with a dynamic mesh method (DMM) to deal with the moving boundary was applied to the simulation of two-dimensional incompressible viscous flow past a circular cylinder that is impulsively started into rotation and translation. The non-dimensional rotating to translating speed ratio α is varied from 0.28 to 2.07, with the Reynolds number being 200 for the range of α. The computation covers a period, during which the cylinder translates seven times its diameter. The current scheme handles the impulsively moving boundary directly by DMM, which is implemented using both mesh deforming and local remeshing. The instantaneous asymmetrical flow configurations for various α are presented and compared with the experimental visualizations. Quantitatively, the velocity distributions with drag and lift coefficients are also compared with the experimental and numerical results. Results show that the flow is strongly influenced by the rotation. Comparisons are found to be satisfactory.     

A finite element approximation of steady flow through a rotating non-aligned straight tube

A finite element solution is developed for a penalty function formulation of the equations which govern the steady motion of a Newtonian fluid through a pipe that rotates about an axis not parallel to its own. The motion in this system is driven by the Coriolis acceleration, which has components in the axial direction as well as in the transverse plane of the pipe. The relative magnitudes of these components significantly affect the qualitative and quantitative nature of the primary and secondary flow field. The present results compare favourably with those of previously reported experimental and theoretical studies over a wide range of flow regimes.     

Comparison between the traction and pressure-imposed boundary conditions in finite element flow analysis

A formulation is developed to impose pressure-prescribed boundary conditions in the penalty finite element method. Some numerical experiments for the Poiseuille flow problem are performed to compare it with the conventional traction-prescribed boundary condition. Also the incorrectness of the traction-free outlet boundary condition for contained-flows is studied with explanatory numerical examples. Discussion is focused on the inlet and outlet boundary conditions to simulate fully developed flows. Finally, the three-dimensional flow in a bifurcated pipe is analysed with the proposed formulation.     

有限长大间隙环流中同心转子动特性系数研究

基于作者建立的大间隙环流中转子运动理论模型,用摄动法推导了有限长大间隙环流流场非线性控制方程的零阶和一阶摄动方程,研究了摄动方程的数值求解方法,并用该数值方法深入研究了有限长大间隙环流中同心转子的动特性系数以及壁面粗糙度、入口压力、长径比和入口预旋等参数的影响,研究结果表明,系统参数对有限长大间隙环流中同心转子动特性系数的影响是流体惯性效应、旋流效应、摩擦耗散效应和Lomakin效应综合影响的因素。     

Transient mixed convection flow instabilities in a vertical pipe

The present paper investigates a numerical study of flow instabilities in transient mixed convection in a vertical pipe. At the entrance of the pipe, the flow is suddenly submitted to a temperature step. The convection heat transfer on the outer surface of the pipe is taken into account. The governing equations are solved using a finite difference explicit scheme. The numerical results show that the time development of streamlines and isotherms is strongly dependent on the inlet temperature steps. For positive temperature steps, the unsteady vortex is significant in the vicinity of the wall and the reversal flow appears below the wave instability. In the case of negative temperature steps and especially for the low Reynolds number, the reversal flow appears on top of the wave instability. During the transient, the apparition of the vortical structures along the wall leads to the wall boundary layer instability. This phenomenon is due to the transient mixed convection flows. The temperature step effects on the heat transfer of the flow are presented in our paper.     

The energy transfer mechanism of a perturbed solid-body rotation flow in a rotating pipe

Three-dimensional direct numerical simulations of a solid-body rotation superposed on a uniform axial flow entering a rotating constant-area pipe of finite length are presented. Steady in time profiles of the radial, axial, and circumferential velocities are imposed at the pipe inlet. Convective boundary conditions are imposed at the pipe outlet.The Wang and Rusak(Phys. Fluids 8:1007–1016, 1996.doi:10.1063/1.86882) axisymmetric instability mechanism is retrieved at certain operational conditions in terms of incoming flow swirl levels and the Reynolds number. However, at other operational conditions there exists a dominant,three-dimensional spiral type of instability mode that is consistent with the linear stability theory of Wang et al.(J. Fluid Mech. 797: 284–321, 2016). The growth of this mode leads to a spiral type of flow roll-up that subsequently nonlinearly saturates on a large amplitude rotating spiral wave. The energy transfer mechanism between the bulk of the flow and the perturbations is studied by the Reynolds-Orr equation. The production or loss of the perturbation kinetic energy is combined of three components: the viscous loss, the convective loss at the pipe outlet, and the gain of energy at the outlet through the work done by the pressure perturbation. Theenergy transfer in the nonlinear stage is shown to be a natural extension of the linear stage with a nonlinear saturated process.     

旋转圆管湍流的大涡模拟数值研究

采用大涡模拟（LES）方法,并结合动力学亚格子尺度应力（SGS）模型,通过数值求解柱坐标系下的滤波Navier-Stokes方程,研究了绕管轴旋转圆管内的湍流流动特性.为验证计算的可靠性,以及动力学SGS模型对于旋转湍流的适用性,将大涡模拟计算所得的结果,与相应的直接模拟（DNS）结果和实验数据进行了对比验证,吻合良好.进一步对旋转圆管湍流的物理机理进行了探讨,研究了湍流特性随旋转速率的变化规律.当旋转速率增加时,湍流流动有层流化的发展趋势.基于湍动能变化的关系,分析了旋转效应对湍流脉动生成的抑制作用.     

A three-dimensional transient two-phase flow analysis with a density pertubation finite volume method

In this paper, we present the analysis of a full three dimensional transient two-phase flow with strong disequilibria. This analysis uses a finite volume scheme with a recently developed density perturbation method. The application test case is a three dimensional extension of a one-dimensional experiment on blowdown pipe. The results are stable with mesh refinements and show realistic three-dimensional effects.     

单柱单锥型液—液旋流分离管内流场的LDV诊断

应用二维激光多普勒仪（ＬＤＶ）对一种单柱单锥型液－液旋流分离管内流场进行了测量,考察了流量、溢流比、压力比和气芯等参数对流场的影响。测量结果表明：切向速度分布呈典型的Ｒａｎｋｉｎｅ涡结构,沿轴向衰减很少,表明所用锥角是合适的;因该旋流管的水力直径较大,切向速度的总体水平较低,由于对了离特性带来了不利影响。此外,没有观察到切向速度分布的的双峰分布现象。轴向速度的总体水平较低,尤其是在锥形管的上游更为     

Computation of fluid flow with a parallel multigrid solver

A finite volume numerical method for the prediction of fluid flow and heat transfer in simple geometries was parallelized using a domain decomposition approach. The method is implicit, uses a colocated arrangement of variables and is based on the SIMPLE algorithm for pressure-velocity coupling. Discretization is based on second-order central difference approximations. The algebraic equation systems are solved by the ILU method of Stone.¹ To accelerate the convergence, a multigrid technique was used. The efficiency was examined on three different parallel computers for laminar flow in a pipe with an orifice and natural convection in a closed cavity. It is shown that the total efficiency is made up of three major factors: numerical efficiency, parallel efficiency and load-balancing efficiency. The first two factors were thoroughly investigated, and a model for predicting the parallel efficiency on various computers is presented. Test calculations indicate reasonable total efficiency and favourable dependence on grid size and the number of processors.     

Rotation effects on a fully-developed turbulent pipe flow

A fully-developed turbulent pipe flow is allowed to pass through a rotating pipe section, whose axis of rotation coincides with the pipe axis. At the exit end of the rotating section, the flow passes into a stationary pipe. As a result of the relaxation of surface rotation, the turbulent flow near the pipe wall is affected by extra turbulence production created by the large circumferential shear strain set up by the rapid decrease of the rotational velocity to zero at the wall. However, the flow in the most part of the pipe is absent of this extra turbulence production because the circumferential strain is zero as a result of the solid-body rotation imparted to the flow by the rotating pipe section. The combined effect of these two phenomena on the flow is investigated in detail using hot-wire anemometry techniques. Both mean and turbulence fields are measured, together with the wall shear and the turbulent burst behavior at the wall. A number of experiments at different rotational speeds are carried out. Therefore, the effects of rotation on the behavior of wall shear, turbulent burst at the wall, turbulence production and the near-wall flow can be documented and analysed in detail.     

Vortex simulation of the intake flow in a planar piston-chamber device

A numerical scheme based on the application of the vortex method to update the vorticity field and the implementation of the finite element method to satisfy the normal velocity boundary condition inside a complex time-dependent geometry is applied to simulate the flow produced by a piston sliding out of a chamber equipped with single or multiple intakes. This unsteady confined vortex flow is of interest in many applications. We use the idealization that the flow is incompressible, two-dimensional and planar and we analyse the results to study the flow during the intake process inside a model of an engine cylinder. The chamber top is fitted with an inlet channel, an inlet port or an inlet valve. In all cases when the intake channel axis coincides with that of the chamber, the flow in each side of the chamber consists essentially of two large counter-rotating eddies of almost the same size. The computed structures of these flows resemble qualitatively those which have been observed experimentally. The fluid motion is also computed for the case of a chamber equipped with an intake whose axis is not aligned with the chamber axis. In this case the flow at the end of the stroke is dominated by a single large eddy produced by the merging of the two eddies forming on the sides of the port.     

The influence of inlet flow condition on the frequency of self-excited jet precession

A precessing jet flow can be generated naturally by a fluidic nozzle comprising a cylindrical nozzle-chamber with a large sudden expansion at its inlet and a small lip at its outlet. Such a precessing jet flow is offset with respect to the chamber axis, about which it rotates. The aim of the present study is to investigate the influence of the chamber-inlet configuration on the frequency of such precession. Three different inlet configurations, classified as long pipe, smooth contraction, and sharp-edged orifice plate, are tested. It is found that the frequency of precession from the orifice is highest, whereas that of the pipe jet is lowest. These differences appear to result partly from the distinct differences in their respective initial boundary layers.     

Flow field modelling near a well with a conductive fracture

In the technology of oil recovery the oil production rate can be increased by generation of a vertical sand-filld conductive fracture on the wall of the well. Oil diffuses through the conductive fracture to the well. In this paper the seepage flow and isothermal deformation fields in both the formation and fracture and the oil production rate at the well are studied by modelling the formation as an infinite poroelastic medium saturated with a one-phase compressible fluid. The fracture is treated as a one-dimensional poroelastic medium. Darcy flows are considered in both the formation and fracture. The plane strain condition is imposed. Our solution is obtained numerically by a finite element method based on a variational principle. The accuracy of the analysis is studied by comparison of the numerical solutions of some problems with their analytical solutions. Since we are dealing with the transient flow problem of an infinite region, an extrapolation technique is employed to find the finite element solution. The production rate of a well with the conductive fracture is compared with that of a well without the conductive fracture.     

Simulation of 2D external viscous flows by means of a domain decomposition method using an influence matrix technique

Two-dimensional external viscous flows are numerically approximated by means of a domain decomposition method which combines a vortex method and a finite difference method. The vortex method is used in the flow region which is dominated by convective effects, whereas the finite difference method is used in the flow region where viscous diffusion effects are dominant. An influence matrix technique combined with the uniformity condition of the pressure is used to enforce the tangential velocity boundary condition. Comparisons between numerical and experimental data show that the method is well adapted for simulating two-dimensional flows.     

Finite element computation of a turbulent flow over a two-dimensional backward-facing step

This paper is devoted to the computation of turbulent flows by a Galerkin finite element method. Effects of turbulence on the mean field are taken into account by means of a k-? turbulence model. The wall region is treated through wall laws and, more specifically, Reichardt's law. An inlet profile for ? is proposed as a numerical treatment for physically meaningless values of k and ?. Results obtained for a recirculating flow in a two-dimensional channel with a sudden expansion in width are presented and compared with experimental values.     

Finite element computation of a turbulent flow over a two-dimensional backward-facing step

The present paper is devoted to the computation of turbulent flows by a Galerkin finite element method. Effects of turbulence on the mean field are taken into account by means of a (k-ε) turbulence model. The wall region is treated through wall laws and, more specifically, Reichardt's law. An inlet profile for ε is proposed as a numerical treatment for physically meaningless values of k and ε. Results obtained for a recirculating flow in a two-dimensional channel with a sudden expansion in width are presented and compared with experimental values.     

On the implementation of symmetric and antisymmetric periodic boundary conditions for incompressible flow

In this paper we consider symmetric and antisymmetric periodic boundary conditions for flows governed by the incompressible Navier-Stokes equations. Classical periodic boundary conditions are studied as well as symmetric and antisymmetric periodic boundary conditions in which there is a pressure difference between inlet and outlet. The implementation of this type of boundary conditions in a finite element code using the penalty function formulation is treated and also the implementation in a finite volume code based on pressure correction. The methods are demonstrated by computation of a flow through a staggered tube bundle.     

Effect of weak rotation on natural convection in a horizontal porous cylinder

Natural convection in a fluid saturated porous medium confined in a horizontal circular cylinder and rotating about its axis, with isothermal boundary conditions and uniform internal heat sink, is studied by both numerical and perturbation methods. No symmetry with respect to the vertical diameter is expected for the flow and temperature fields and the whole region must be involved in the computation. Only the weak rotation regime, for which the centrifugal force is negligible compared to gravity, is considered. Governing equations for the two-dimensional flow field are solved in both rotating and non-rotating coordinate systems. Results indicate that rotation significantly decreases the radial amplitude of fluid particle trajectories in the radial direction and thus reduces the overall heat transfer.