Simulation of single‐ and two‐phase flows on sliding unstructured meshes using finite volume method

The paper presents an efficient finite volume method for unstructured grids with rotating sliding parts composed of arbitrary polyhedral elements for both single‐ and two‐phase flows. Mathematical model used in computations is based on the ensemble averaged conservation equations. These equations are solved for each phase and in case of single‐phase flow reduce to the transient Reynolds‐averaged Navier–Stokes (TRANS) equations. Transient flow induced by rotating impellers is thus resolved in time. The use of unstructured grids allows an easy and flexible meshing for the entire flow domain. Polyhedral cell volumes are created on the arbitrary mesh interface placed between rotating and static parts. Cells within the rotating parts move each time step and the new faces are created on the arbitrary interfaces only, while the rest of the domain remain ‘topologically’ unchanged. Implicit discretization scheme allows a wide range of time‐step sizes, which further reduce the computational effort. Special attention is given to the interpolation practices used for the reconstruction of the face quantities. Mass fluxes are recalculated at the beginning of each time step by using an interpolation scheme, which enhances the coupling between the pressure and velocity fields. The model has been implemented into the commercially available CFD code AVL SWIFT (AVL AST, SWIFT Manual 3.1, AVL List GmbH, Graz, Austria, 2002). Single‐phase flow in a mixing vessel stirred by a six‐bladed Rushton‐type turbine and two‐phase flow in aerated stirred vessel with the four‐blade Rushton impeller are simulated. The results are compared with the available experimental data, and good agreement is observed. The proposed algorithm is proved to be both stable and accurate for single‐phase as well as for the two‐phase flows calculations. Copyright 2004 John Wiley & Sons, Ltd.     

Grid‐free surface vorticity method applied to flow induced vibration of flexible cylinders

In order to study cross flow induced vibration of heat exchanger tube bundles, a new fluid–structure interaction model based on surface vorticity method is proposed. With this model, the vibration of a flexible cylinder is simulated at Re=2.67 × 10⁴, the computational results of the cylinder response, the fluid force, the vibration frequency, and the vorticity map are presented. The numerical results reproduce the amplitude‐limiting and non‐linear (lock‐in) characteristics of flow‐induced vibration. The maximum vibration amplitude as well as its corresponding lock‐in frequency is in good agreement with experimental results. The amplitude of vibration can be as high as 0.88D for the case investigated. As vibration amplitude increases, the amplitude of the lift force also increases. With enhancement of vibration amplitude, the vortex pattern in the near wake changes significantly. This fluid–structure interaction model is further applied to simulate flow‐induced vibration of two tandem cylinders and two side‐by‐side cylinders at similar Reynolds number. Promising and reasonable results and predictions are obtained. It is hopeful that with this relatively simple and computer time saving method, flow induced vibration of a large number of flexible tube bundles can be successfully simulated. Copyright © 2004 John Wiley & Sons, Ltd.     

Transient free‐surface flow of a viscoelastic fluid in a narrow channel

The interplay between inertia and elasticity is examined for transient free‐surface flow inside a narrow channel. The lubrication theory is extended for the flow of viscoelastic fluids of the Oldroyd‐B type (consisting of a Newtonian solvent and a polymeric solute). While the general formulation accounts for non‐linearities stemming from inertia effects in the momentum conservation equation, and the upper‐convected terms in the constitutive equation, only the front movement contributes to non‐linear coupling for a flow inside a straight channel. In this case, it is possible to implement a spectral representation in the depthwise direction for the velocity and stress. The evolution of the flow field is obtained locally, but the front movement is captured only in the mean sense. The influence of inertia, elasticity and viscosity ratio is examined for pressure‐induced flow. The front appears to progress monotonically with time. However, the velocity and stress exhibit typically a strong overshoot upon inception, accompanied by a plug‐flow behaviour in the channel core. The flow intensity eventually diminishes with time, tending asymptotically to Poiseuille conditions. For highly elastic liquids the front movement becomes oscillatory, experiencing strong deceleration periodically. A multiple‐scale solution is obtained for fluids with no inertia and small elasticity. Comparison with the exact (numerical) solution indicates a wide range of validity for the analytical result. Copyright © 2004 John Wiley & Sons, Ltd.     

Inviscid,laminar and turbulent opposed flows

This paper attempts to reproduce numerically previous experimental findings with opposed flows and extends their range to quantify the effects of upstream pipes and nozzles with inviscid, laminar and turbulent flows. The choice of conservation equations, boundary conditions, algorithms for their solution, the degree of grid dependence, numerical diffusion and the validity of numerical approximations are justified with supporting calculations where necessary. The results of all calculations on the stagnation plane show maximum strain rates close to the annular exit from the nozzles and pipes for lower separations and it can be expected that corresponding reacting flows will tend to extinguish in this region with the extinction moving towards the axis. With laminar flows, the maximum strain rate increased with Reynolds number and the maximum values were generally greater than with inviscid flows and smaller than with turbulent flows. With large separations, the strain rates varied less and this explains some results with reacting flows where the extinction appeared to begin on the axis. The turbulent‐flow calculations allowed comparison of three common variants of a two‐equation first‐moment closure. They provided reasonable and useful indications of strain rates but none correctly represented the rms of velocity fluctuations on the axis and close to the stagnation plane. As expected, those designed to deal with this problem produced results in better agreement with experiment but were still imperfect. Copyright © 2004 John Wiley & Sons, Ltd.     

低速二维流动的粒子仿真研究

本文运用信息保存法对低速二维的流动现象进行模拟，考察了低速条件下的有限平板绕流以及微槽道气体流动问题。研究表明：在对低速流动的模拟过程中，运用IP法在能够获得较好的结果的同时，具有比DSMC方法更高的计算效率。     

Temperature factor effect on the structure of the separated flow within a supersonic gas stream

The effect of the temperature factor, that is, the ratio of the body temperature to the freestream stagnation temperature, on the structure of the separated flow formed in the presence of a concave corner in a supersonic stream is studied. The strong influence of the temperature factor on the separation zone length and the flow-generated aerodynamic characteristics is established. It is shown that for fairly large deflection angles this flow cannot be described by free interaction, or triple deck, theory.     

Analysis and finite element simulations of a second‐order fluid model in a bounded domain

This article is concerned with the equations governing the steady motion of a viscoelastic incompressible second‐order fluid in a bounded domain. A new proof of existence and uniqueness of strong solutions is given. In addition, using appropriate finite element methods to approximate a coupled equivalent problem, sharp error estimates are obtained using a fixed point argument. The method is applied to the two‐dimensional lid‐driven cavity problem, at low Reynolds number and in a certain range of values of the viscoelastic parameters, to analyze the combined effects of inertia and viscoelasticity on the flow. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

乙醇增敏光泽精－H_2O_2－Pb（Ⅱ）化学发光体系的研究及应用

本文研究了Ｈ＿２Ｏ＿２氧化光泽精发光反应及乙醇的敏化作用，发现加入乙醇可使铅催化Ｈ＿２Ｏ＿２，氧化光泽精的发光强度有较大幅度的增强；铅量在３×１０￣（－７）～８×１０￣（－５）ｇ／ｍＬ。范围内与发光强度成正比．方法检出限为０．１０μｇ／ｍＬ；对含４μｇ／ｍＬＰｂ（Ⅱ）１１次重复测定，其ＲＳＤ为２．３％．对３２种常见的阴阳离子进行了干扰实验，表明方法具有较好的选择性。本法应用于水样及沉积物标样中痕量铅的测定，结果良好。     

Use of dynamic trees in a network simplex algorithm for the maximum flow problem

Goldfarb and Hao (1990) have proposed a pivot rule for the primal network simplex algorithm that will solve a maximum flow problem on ann-vertex,m-arc network in at mostnm pivots and O(n ² m) time. In this paper we describe how to extend the dynamic tree data structure of Sleator and Tarjan (1983, 1985) to reduce the running time of this algorithm to O(nm logn). This bound is less than a logarithmic factor larger than those of the fastest known algorithms for the problem. Our extension of dynamic trees is interesting in its own right and may well have additional applications.Research partially supported by a Presidential Young Investigator Award from the National Science Foundation, Grant No. CCR-8858097, an IBM Faculty Development Award, and AT&T Bell Laboratories.Research partially supported by the Office of Naval Research, Contract No. N00014-87-K-0467.Research partially supported by the National Science Foundation, Grant No. DCR-8605961, and the Office of Naval Research, Contract No. N00014-87-K-0467.     

Preparation of Papers

We motivate the study of a vector variational inequality by a practical flow equilibrium problem on a network, namely a generalization of the well-known Wardrop equilibrium principle. Both weak and strong forms of the vector variational inequality are discussed and their relationships to a vector optimization problem are established under various convexity assumptions.