圆管绕流旋涡脱落诱导振动的数值模拟

利用计算流体力学软件Ansys／Flotran CFD,首先对粘性不可压缩流体的固定圆管绕流进行了数值模拟,然后结合逐步积分法完成了同时考虑纵横两向弹性支撑圆管绕流旋涡脱落诱导振动的数值模拟,并通过快速傅立叶变换,得到了弹性支撑圆管和固定圆管的升力及弹性支承圆管横向位移响应的功率谱．通过计算结果分析,得出了一些有价值的结论,可供从事具有圆管绕流构件设备设计的工程技术人员参考．     

Slow viscous flow around two particles in a cylinder

This article describes the motion of two arbitrarily located free moving particles in a cylindrical tube with background Poiseuille flow at low Reynolds number. We employ the Lamb’s general solution based on spherical harmonics and construct a framework based on cylindrical harmonics to solve the flow field around the particles and the flow within the tube, respectively. The two solutions are performed in an iterated framework using the method of reflections. We compute the drag force and torque coefficients of the particles which are dependent on the distances among the cylinder wall and the two particles. In addition, we provide detailed flow field in the vicinity of the two particles including streamlines and velocity contour. Our analysis reveals that the particle–particle interaction can be neglected when the separation distance is three times larger than the sum of particles radii when the two particles are identical. Furthermore, the direction of Poiseuille flow, the particle position relative to the axis and the particle size can make the two particles attract or repel. Unlike the single particle case, the two particles can move laterally due to the hydrodynamic interaction. Such analysis can give insights to understand the mechanisms of collision and aggregation of particles in microchannels.     

Laminar viscous flow fields in moving curved channels

The laminar steady flow of air through moving narrow deep circumferential grooves on a rotor is studied. Fluid enters the moving grooves from a stationary inlet nozzle with uniform velocity and exits into a stationary diffuser after approximately 150 degrees of rotation. Three dimensional developing velocity profiles and pressure distributions are presented for a variety of operating conditions. Computed overall head-flow characteristics are found to be in relatively close agreement with a simple Poiseuille model where both walls are taken to be moving at the same rate and in the same direction at a velocity different from the entrance velocity.     

Calculating viscous incompressible flow past an elliptical cylinder on a parallel computer

An algorithm developed for a parallel computer is described. The algorithm has been simulated on an ES-1060 machine, and the results are compared with those obtained on BÉSM-6 computer.Translated from Kibernetika, No. 1, pp. 64–68, January–February, 1990.     

Numerical simulation of viscous flow over non-smooth surfaces

The incompressible viscous flow over several non-smooth surfaces is simulated numerically by using the lattice Boltzmann method. A numerical strategy for dealing with a curved boundary with second-order accuracy for velocity field is presented. The drag evaluation is performed by the momentum-exchange method. The streamline contours are obtained over surfaces with different shapes, including circular concave, circular convex, triangular concave, triangular convex, regular sinusoidal wavy and irregular sinusoidal wavy, are obtained. The flow patterns are discussed in detail. The velocity profiles over different kinds of non-smooth surfaces are investigated. The numerical results show that the lattice Boltzmann method is reliable, accurate, easy to implement, and can provide valuable help for some engineering practices.     

Numerical study of the viscous flow in oscillatory spherical annuli

The method of fractional steps is applied for investigating the viscous flow field in the oscillatory spherical annuli. The numerical solution obtained enlarges the intervals of the frequency parameter M and amplitude of the torsional oscillations ε in comparison with known solutions. Typical flow fields are shown graphically.     

小型轴流风扇内外部流场的数值研究

小型轴流风扇的性能及降噪研究建立在流场研究的基础之上。通过对风扇VD8025的Fluent数值模拟,展开了其内外部流场的分析研究,结论如下：在叶片吸力面（非工作面）侧存在主要由边界层分离引起的分离涡,在叶片压力面（工作面）侧顶部间隙存在二次流引起的叶顶间隙涡,沿轴向存在有限叶高诱导双旋涡。对风扇内部流场的速度分布进行了定量分析,分析的结果说明了上述旋涡等流动情况的存在。由风扇出口管道速度分布的定量分析可得,对于所研究的风扇,距离风扇出口120 mm处为排风最佳位置。     

Steady flow around and through a permeable circular cylinder

The steady flow around and through a porous circular cylinder was studied numerically. The effects of the two important parameters, the Reynolds and Darcy numbers, on the flow were investigated in details. The recirculating wake existing downstream of the cylinder is found to either penetrate into or be completely detached from the cylinder. It is also found that, contrary to that of the solid cylinder, the recirculating wake develops downstream of or within the porous cylinder, but not from the surface of it. These new findings provide additional evidence to Leal’s conclusion (Leal LG. Vorticity transport and wake structure for bluff bodies at finite Reynolds number. Phys Fluids A 1989;1:124) that the appearance of recirculating wakes at finite Reynolds number is due to vorticity accumulation, but not a result of the same physical phenomena associated with separation in boundary layers in adverse pressure gradients. Also presented in the current study are the variation of the critical Reynolds number for the onset of a recirculating wake as a function of Darcy number and the variation of a newly defined parameter, the penetration depth, as a function of the Reynolds number and Darcy number.     

Simulation of a viscous compressible flow past a circular cylinder with high-order discontinuous Galerkin methods

This paper is devoted to the simulation of viscous compressible flows with high-order accurate discontinuous Galerkin methods on bidimensional unstructured meshes. The formulation for the solution of the Navier–Stokes equations is due to Oden et al. [An hp-adaptive discontinuous finite element method for computational fluid dynamics. PhD thesis, The University of Texas at Austin, 1997; J Comput Phys 1998;146:491–519]. It involves a weak imposition of continuity conditions on the state variables and on fluxes across interelement boundaries. It does not make use of any auxiliary variables and then the cost for the implementation is reasonable. The method is coupled with a limiting procedure recently developed by the authors to suppress oscillations near large gradients. The limiter is totally free of problem dependence and maintains the convergence order for errors measured in the L¹-norm. This paper presents detailed numerical results of a viscous compressible flow past a circular cylinder at a Reynolds number of 100 for the cases of subsonic and supersonic regimes. The proposed simulations suggest that the method is very robust and is able to produce very accurate results on unstructured meshes.     

Numerical study on two-phase flow through fractured porous media

Based on the flux equivalent principle of a single fracture,the discrete fracture concept was developed,in which the macroscopic fractures are explicitly described as(n-1)dimensional geometry element.On the fundamental of this simplification,the discrete-fractured model was developed which is suitable for all types of fractured porous media.The principle of discrete-fractured model was introduced in detail,and the general mathematical model was expressed subsequently.The fully coupling discrete-fractured ma...     

基于遗传算法的圆柱阵列稀疏方法

由于共形阵列所具有的特性,使其正得到日益广泛的应用,但其耗用较多的阵元,方向图具有相对主瓣较高的旁瓣电平。为此,针对基本的共形阵列——圆柱阵列的天线阵元,应用经典遗传算法,以阵元的工作状态为优化参量,对其进行稀疏。减少了阵元数量,仿真结果表明,该方法能有效地降低圆柱阵列的旁瓣电平。     

Computer analysis of the aerodynamics and heat exchange of a cylinder in a viscous oscillating flow

The paper considers the results of analyzing the aerodynamics and heat exchange of a fixed cylinder under nonstationary airflow with average-velocity alternating pulsations numerically. It is shown that the pulsations significantly affect the vortex separation frequency, the Karman vortex street behind the cylinder, the pulsation frequency of the aerodynamic and heat exchange coefficients, and their average and amplitude values. The authors determined the engineering applications.     

Numerical studies of steady,viscous, incompressible flow between two rotating spheres

A new numerical method, developed for the study of secondary flow in a curved tube, is adapted and extended to the study of viscous, incompressible, steady flow between two rotating spheres. The Navier-Stokes equations are approximated by a triple sequence of linear problems, each of which has a diagonally dominant coefficient matrix. Computer examples are described and discussed.     

Numerical solution of incompressible flow through branched channels

The work deals with numerical solution of 3D turbulent flow in straight channel and branched channels with two outlets. The mathematical model of the flow is based on Reynolds-averaged Navier–Stokes equations for incompressible flow in 3D with explicit algebraic Reynolds stress turbulence model (EARSM). The mathematical model is solved by artificial compressibility method with implicit finite volume discretization. The channels have constant square or circular cross-section, where the hydraulic diameter is same in order to enable comparison between these numerical simulations. First, developed flow in a straight channel of square cross-section is presented in order to show the ability of the used EARSM turbulence model to capture secondary corner vortices, which are not predicted by eddy viscosity models. Next the flow through channels with perpendicular branch is simulated. Methods of setting the flow rate are discussed. The numerical results are presented for two flow rates in the branch.     

Time-dependent solutions of the viscous incompressible flow past a circular cylinder by the method of series truncation

Semi-analytic solutions of the Navier-Stokes equations are calculated for two-dimensional, symmetrical, viscous incompressible flow past a circular cylinder. The stream and vorticity functions are expanded in the finite Fourier series and then substituted in the Navier-Stokes equations. This led to a system of coupled parabolic partial differential equations which are solved numerically. More terms of the series are required as Reynolds number increases and the present calculations were terminated at Reynolds number 600 with 60 terms of Fourier series. The results are compared with similar calculations and experimental data for Reynolds numbers 60, 100, 200, 500, 550 and 600. At the termination of the calculations for Reynolds numbers 60 and 100, the separation angle, the wake length, the drag coefficient, and the vorticity distributions around the surface were very close to their steady-state values. A secondary vortex appeared on the surface of the cylinder in the case of Reynolds numbers 500, 550 and 600. The wake length, the drag coefficient and the separation angle differ slightly at a given instant in the case of Reynolds numbers 500, 550 and 600.     

轻汽油醚化反应器内流动状况的数值模拟

为了解决某石化公司轻汽油醚化反应器的放大问题,本研究采用CFD方法对轻汽油醚化反应器内的流动状况进行模拟研究。通过模拟计算得到滴流床反应器内气液两相流动状况,为反应器的选型与设计提供依据。研究表明,该醚化反应器中的流动处于滴流床流型,气液两相均以连续相的形式存在。     

Numerical parameter identification in multiphase flow through porous media

Multiphase flow is of high interest for the investigation of the behavior of waste in groundwater. The high nonlinearity of the model equations pose special problems. Here, a new parameter identification technique in this context is proposed which takes advantage of recently developed highly efficient numerical simulation techniques. It is based on a reduced Gauss-Newton technique in combination with an efficient gradient computation. Numerical experiments are performed for the McWhorter model problem.     

Numerical simulation of flow through microchannels with designed roughness

A three-dimensional numerical simulation of flow through serpentine microchannels with designed roughness in form of obstructions placed along the channels walls is conducted here. CFD-ACE+ is used for the numerical simulations. The effect of the roughness height (surface roughness), geometry, Reynolds number on the friction factor is investigated. It is found that the friction factor increases in a nonlinear fashion with the increase in obstruction height. The friction factor is more for rectangular and triangular obstructions and it decreases as the obstruction geometry is changed to trapezoidal. It is observed that the obstruction geometry, i.e., aspect ratio plays an important role in prediction of friction factor in rough channels. It is also found that the pressure drop decreases with the increase in the roughness pitch. Hence, the roughness pitch is an important design parameter for microchannels.     

Numerical simulation of oscillating flow through idealized sclerotic arteries

Numerical and approximate computational techniques are presented for studying the oscillating flow in a two-dimensional channel with an irregular surface. A finite element algorithm is formed in order to evaluate the influence of the surface roughness on the basic flow characteristics, while emphasis is given to the development and variation of the shear stresses on the walls. In addition, an approximate solution is obtained when the wall roughness is small compared with the mean height of the channel. Comparison of the two methods in an appropriate range of wall roughness values gives satisfactory results.     

Distributed finite element calculations on transputer arrays and the DAP

An equation solver based on the preconditioned conjugate gradient method for the sparse system arising from finite element analysis is presented. The preconditioning matrix has been designed to take advantage of the domain decomposition approach used on local memory multiprocessor computers. The method has been implemented on a transputer array and on the DAP; results are given for these computers. A simple domain decomposition algorithm is also presented. This method is suitable for the decomposition of finite element meshes for the transputer based analysis program.