COMPUTATION OF DILUTE PARTICULATE LAMINAR FLOW OVER A BACKWARD-FACING STEP

Particle-laden flows are calculated for a classical laminar backward- facing step problem. The particle tracks are calculated using a recently developed exponential Lagrangian tracking scheme. The behaviour of the particle-laden flow is considered for various inlet for Reynolds number, Stokes numbers and void fractions. Doping the flow with low-Stokes-number particles has the effect of increasing the inlet inertia of the flow and this increases the strength of the recirculation behind the step. High-Stokes- number particles are dominated by gravitational effects which affect the flow accordingly. Differences between the single-phase flow and the particle-laden flows are therefore dependent on the Stokes number and increase linearly with void fraction.     

紧致插值曲线方法在流向强迫振荡圆柱绕流中的应用

利用紧致插值曲线(constrained interpolation profile method in Zhejiang University, CIP-ZJU) 数学模型, 对低科勒冈-卡朋特(Keulegan–Carpenter) 数KC 静止流体中振荡圆柱以及雷诺数Re = 200 时流向强迫振荡圆柱绕流进行了数值模拟. 模型在直角坐标系统下建立, 采用紧致插值曲线方法作为流场的基本求解器离散了纳维-斯托克斯方程, 基于多相流的理论实现流固耦合同步求解, 利用浸入边界方法处理固体边界. 模拟结果与现有文献结果进行比较, 二者吻合情况较好, 验证了此方法对于计算复杂流动问题的可靠性.      

A NUMERICAL STUDY OF FLOW PAST A ROTATIONALLY OSCILLATING CIRCULAR CYLINDER

Vortex shedding from a rotationally oscillating circular cylinder in a uniform flow is studied by numerical solutions of the two-dimensional incompressible Navier-Stokes equations using primitive variables. To demonstrate the viability and accuracy of the method, we calculate the approach flow past a rotating cylinder with constant angular velocity. For approach flow past a rotationally oscillating cylinder, the object of the study is to examine the effect of oscillating rotation on the flow structure. In the present study, some basic patterns of vortex shedding can be identified according to our calculated results and are in good agreement with available experiments. In addition, the influence of the oscillating frequency and amplitude on the forces acting on the cylinder is also investigated.     

APPLICATION OF THE LES METHOD TO THE OSCILLATING FLOW PAST A CIRCULAR CYLINDER

The oscillating flow of an incompressible fluid over a circular cylinder is investigated by solving the vorticity/stream-function version of the two-dimensional Navier Stokes equations using a finite-difference Large Eddy Simulation (LES) method. Two different subgrid scale (SGS) models are tested. They are the classical Smagorinsky model and Yoshizawa's Two-Scale Direct-Interaction Approximation (TSDIA) model, both of which require the input of constants in the specification of the subgrid scale viscosity. In addition, a solution was obtained with no subgrid scale model for comparison purposes. Yoshizawa's model is based on the combination of Kraichman's Direct Interaction Analysis (1964) with the separation of scales of mean (grid scale) and fluctuating (subgrid scale) fields. This TSDIA model has the length scale as a function of space and time so that it incorporates the effects of local time rate of change and advection of the resolvable field. Both models are applied to the classical problem of an oscillating two-dimensional flow over a circular cylinder at various KC values atβ =1035.Calculated lift and in-line force coefficients compare favorably with available experimental data for KC≤2 with the best solutions developed from the Yoshizawa model. The solution for KC>2 demonstrates that the two-dimensional model is inadequate to describe the flow accurately.     

弹性支撑圆柱绕流稳定性分析

基于CFD 技术,采用系统辨识方法,建立了亚临界雷诺数(Re < 47) 下绕圆柱流动的非定常气动力模型(reduced order model, ROM). 耦合结构运动方程和降阶气动力模型,建立了弹性支撑圆柱绕流的稳定性分析模型. 算例分析了亚临界雷诺数下,结构固有频率、质量比等参数以及支撑方式对弹性系统稳定性的影响. 对于单自由度横向支撑圆柱,当结构固有频率趋近流动最不稳定模态频率时,弹性系统会在一定频率范围内失稳,这种现象最低可在Re~20 时出现. 旋转自由度的释放能够进一步降低系统的稳定性,可将临界雷诺数进一步降低至18 左右. ROM 方法不仅具有很高的效率,而且清晰地指出了弹性系统失稳的根本原因:流动模态和结构模态耦合作用导致结构模态失稳所致. 因此,失稳状态下系统振荡频率锁定于结构固有频率. 基于ROM 技术预测的失稳边界与直接CFD/CSD 仿真结果吻合,证明了该方法的正确性和精度.      

UNSTEADY CONFINED VISCOUS FLOWS WITH OSCILLATING WALLS AND MULTIPLE SEPARATION REGIONS OVER A DOWNSTREAM-FACING STEP

This paper presents computational solutions for unsteady viscous flows in channels with a downstream-facing step, followed by an oscillating floor. These solutions of the unsteady Navier–Stokes equations are obtained with a time-integration method using artificial compressibility in a fixed computational domain, which is obtained via a time-dependent coordinate transformation from the fluid domain with moving boundaries. The computational method is first validated for steady viscous flows past a downstream-facing step by comparison with previous numerical solutions and experimental results. This method is then used to obtain solutions for unsteady viscous flows with multiple separation regions over a downstream-facing step with oscillating walls, for which there are no previously known solutions. Thus, the present results may be used as benchmark solutions for the unsteady viscous flows with multiple separation regions between fixed and oscillating walls.     

VORTICES INCIDENT UPON AN OSCILLATING CYLINDER: FLOW STRUCTURE AND LOADING

Interaction of an incident vortex street with an oscillating cylinder is addressed using high-image-density particle image velocimetry and simultaneous force measurements. This approach reveals that the timing of the incident vortices relative to the cylinder motion controls the large-scale vortex formation in the near-wake, and thereby the phase shift between the loading on the cylinder and its motion. As a consequence, it is possible to change the sign of the fluid-dynamic work done by the fluid on the cylinder. The incident vortices dramatically shorten the formation length of vortices in the near-wake and yield values of lift coefficient up to a factor of five larger than that for an isolated cylinder subjected to controlled oscillations in the absence of incident vortices. These alterations of the wake structure and loading occur in conjunction with globally locked-on patterns of incident and shed vortices with respect to the cylinder oscillation. Different states of global lock-on are attainable for different values of timing of the incident vortices.     

纳米尺度圆柱绕流尾迹区流动形式模拟研究

采用非平衡分子动力学模拟方法,对微尺度低{Re}数下的圆柱绕流问题进行了研究,模拟结果表明:当{Re}<12时,圆柱下游形成对称、无分离的定常流;当{Re}>20时,圆柱下游形成周期性交替出现的对称涡;当12相似文献   

基于波尔兹曼方法研究幂律流体的椭圆柱绕流

基于格子波尔兹曼方法(lattice Boltzmann method) 和幂律流体本构方程, 建立二维流动模型, 将充分发展的速度分布与理论解进行对比, 吻合良好. 对幂律流体的圆柱绕流进行模拟, 采用了反弹格式的无滑移圆柱边界, 并使用应力积分法计算阻力系数, 分析了稳态圆柱绕流时, 阻力系数随幂律指数n 以及雷诺数Re 的变化规律. 分析了椭圆横轴/纵轴长度比和幂律指数n, 对压力系数Cp 和黏度系数Cv 的影响. 得到的变化规律与有限元方法规律一致, 验证了格子波尔兹曼模型的可行性.     

流向强迫振荡圆柱绕流的涡脱落模态分析

通过求解采用ALE方法描述的运动坐标系Navier-Stokes方程组,分析均匀来流下雷诺 数为150的静止和流向振荡的圆柱绕流. 主要研究了强迫振荡频率和较大振幅比 (A/D=0.3-1.2)对圆柱升力、阻力变化特性以及涡脱落模态的影响. 研究表 明,流向振荡圆柱绕流存在多种涡脱落模态,如对称S以及反对称A-I, A-III, A-IV等多种形式;比较研究结果,拓展了各模态下对应的锁定区域,并将其分为5个 子区;A-I模态中圆柱受力较以前所知更复杂;通过分析计算结果,发现最大加速度 比Af_{c}^{2}/Df_{s0}^{2}可能是涡脱落模态(尤其是对称S模态)最有效的控制参数.     

振动翼型和襟翼的绕流特性研究

用数值模拟手段详细地研究了振动翼型和襟翼的绕流问题,数值模拟的出发方程为Euler和N-S方程,格式为Bcam-Warming格式的改进型。数值实验主要针对流场的二大特性进行的,即振动对激波的影响和振动对分离的抑制作用,结果表明:(1)随翼型或襟翼的振动激波强度和位置也相应地变化但这一变化滞后于攻角的变化;(2)振幅加大激波强度的变化和激波运动范围也加大;(3)振动频率越高对激波的影响反而较低频时要小;(4)流动条件的不同可使升力回线的走向发生变化;(5)振动对分离有明显的抑制作用。     

应用共轴型二维激光测速系统测量孔板管流的湍流特性

本文发表了一种共轴型二维激光测速系统,可同时测量由三束入射光组成的平面内的二维速度分量。讨论了主要的测量误差并提出了一种修正共轴分量角度偏差的方法。应用该系统详细测量了单孔板和双孔板管流的轴向和径向平均速度。湍流度和雷诺切应力分布,表明来流条件对孔板下游的湍流特性有强烈影响。     

后向台阶湍流拟序结构的实验研究

在水槽中用氢气泡法观测了一个后向台阶的回流区、再附区和边界层再生区的大涡和湍流团的发展,以及近壁区的流动结构,观察到大涡结构对湍流团产生的作用。     

EXPERIMENTAL FLOW STUDY WITHIN A SELF OSCILLATING COLLAPSIBLE TUBE

The flow field within a self-excited flexible tube was studied by employing flow visualization, velocity and pressure measurements. Under low positive tansmural pressures at the tube inlet (of the order of 50 mm H₂O) the tube was set to an oscillatory motion, the initiation of which was due to a flow asymmetry. Namely, although initially the flow was separated from both sides of the formed divergent nozzle close to the tube outlet, at an arbitrary instant this became attached to one side, but stayed separated in the remaining part of the nozzle. When this flow asymmetry occurred, the walls approached each other and the tube neck formed there by started oscillating streamwise, setting the tube to an oscillatory motion. During the downstream motion of the tube neck, it was shrunk, causing a twofold increase of the local velocities compared to the inlet ones, which remained almost constant. On the contrary, when moving upstream, the tube neck expanded, causing a flow reversal in this area and a flow deceleration in the remaining part of the tube. The pressure signals upstream and downstream of the tube exhibited a phase difference, the latter leading, taking an order of magnitude higher values than the first one.     

蜿蜒边界下平面流的线性流动稳定性

为便于数值分析,蜿蜒河流水动力和演变模型中一般隐性假设二次时均流-二次涡的关系与明渠流时均流-明渠湍流的关系相同,但由于高雷诺数下的DNS算力限制和实验尺度限制,这种隐含假设是否成立目前尚无相关湍流研究来支撑.文章试图通过分析明渠湍流和二次湍流发展初期的研究,侧面揭示其湍流结构的异同.通过对曲线正交坐标系下的平面二维NS方程使用双参数摄动的方法,建立了一种求解蜿蜒边界弱非线性层流的摄动解法,并推导得出一个适用于蜿蜒边界的EOS方程以及其特征值问题的解法.蜿蜒边界下弱非线性层流解为一系列蜿蜒谐波分量的叠加,其中线性部分使得两壁产生流速差,非线性部分随着雷诺数增大呈指数增长.水流的扰动增长率特征谱的第一模态与直道流相似,由3条曲线、4个波段合成,但其长波段和短波段的扰动流场与直道流不同,所有短波段的扰动流速近似于KH涡.蜿蜒边界对内部水流扰动有一定的选择性.偏角幅值越大扰动增长越快;蜿蜒波数的影响则为先增后减,有一个使扰动增长最快的蜿蜒波数.扰动流场由一个典型的TS波和一对波包形式的二次涡叠加而成,波包只有纵向流速分量,包络线由蜿蜒波数控制,波包内是与直道扰动波参数相同的TS波.     

基于比拟理论的翼型扰流声场数值模拟

从二维模型方程的全离散形式出发,重点分析了差分格式的色散特性和各向异性效应,证实迎风紧致格式比对称格式有更好的色散和各向同性特性,故有利于声场的数值模拟,并采用三阶迎风紧致格式（ＵＣＤ３）和四阶对称紧致格式（ＳＣＤ４）计算了绕ＮＡＣＡ００１２翼型的可压缩非定常流场,并将此流场作为近场声源,运用声学比拟理论对气动声进行模拟。     

壁面传质扰动有纵向速度的流动问题

研究壁面有周期性定向抽吸－引射且在壁面上形成纵向速度的二维渠道流动,壁面上的纵向平均速度为〈uω〉．数值结果表明,抽吸－引射的传质倾角θ对流场的性质和壁面上的切应力等有重要的影响．与速度为〈uω〉的运动壁无扰动流动相比,它的阻力偏低,这种阻力偏低与流场的扰动特性有关,即与扰动速度分量的二重相关积分I成正比．阻力在θ＞-24°的范围内有减阻效果;能量在θ＞8°时有净能量减小的效果．     

NUMERICAL SIMULATION OF STEADY FLOW IN A COMPLIANT TUBE OR CHANNEL WITH TAPERED WALL THICKNESS

Flow through compliant tubes with linear taper in wall thickness is numerically simulated by finite element analysis. Two models are examined: a compliant channel and an axisymmetric tube. For verification of the numerical method, flow through a compliant stenotic vessel is simulated and compared to existing experimental data. Steady two-dimensional flow in a collapsible channel with initial tension is also simulated and the results are compared with numerical solutions from the literature. Computational results for an axisymmetric tube show that as cross-sectional area falls with a reduction in downstream pressure, flow rate increases and reaches a maximum when the speed index (mean velocity divided by wave speed) is near unity at the point of minimum cross-sectional area, indicative of wave-speed flow limitation or “choking” (flow speed equals wave speed) in previous one-dimensional studies. For further reductions in downstream pressure, the flow rate decreases. Cross-sectional narrowing is significant but localized. For the particular wall and fluid properties used in these simulations, the area throat is located near the downstream end when the ratio of downstream-to-upstream wall thickness is 2; as wall taper is increased to 3, the constriction moves to the upstream end of the tube. In the planar two-dimensional channel, area reduction and flow limitation are also observed when outlet pressure is decreased. In contrast to the axisymmetric case, however, the elastic wall in the two-dimensional channel forms a smooth concave surface with the area throat located near the mid-point of the elastic wall. Though flow rate reaches a maximum and then falls, the flow does not appear to be choked.     

NUMERICAL STUDY OF HIGH-REYNOLDS-NUMBER FLOW PAST A BLUFF OBJECT

A semi-explicit finite difference scheme is proposed to study unsteady two-dimensional, incompressible flow past a bluff object at high Reynolds number. The bluff object comes from a class of elliptical cylinders in which the aspect ratio and the angle of attack are two controlled parameters. Associated with the streamfunction–vorticity formulation, the interior vorticity, streamfunction and wall vorticity are updated in turn for each time step. The streamfunction and wall vorticity are solved by means of a multigrid method and a projection method respectively. In regard to the vorticity transport equation, implicitness is merely associated with the diffusion operator, which can be made semi-explicit via approximate factorization. Low-diffusive upwinding is devised to handle the convection part. Numerical results are reported for Reynolds numbers up to 40,000. Comparisons with other numerical or physical experiments are included.