两并列方柱绕流相互干扰的数值研究

用数值方法研究了绕流两并列方柱的气动力及其相互干扰.使用了陈素琴等人提出的对流项用三阶逆风格式的带多重网格的改进的MAC方法,对两并列方柱的流场进行了模拟.结果表明:两并列方柱在不同的间距比时其干扰特征有很大的不同,小间距比时出现双稳态偏流,流动并不对称于双方柱间隙的中心轴线,而是偏向其中的一个方柱.     

含压力基Navier-Stokes方程最优动力系统建模和分析

研究了采用压力基函数和速度基函数的Navier-Stokes方程的最优截断低维动力系统建模理论.在黏性不可压缩流体中模拟了并排三方柱绕流流场,对此流场进行了含压力基函数和速度基函数的Navier-Stokes方程的最优动力系统建模,并以此为工具分析了三方柱绕流最优动力系统的动力学特性.该研究得到了如下结论:三方柱绕流的最优动力系统的动力学行为为混沌,它与双方柱绕流场的极限环动力学特性有着本质的区别,因此可以通过多柱绕流增进尾流的复杂性,从而促进流体混合.     

Navier-Stokes方程的脉动速度方程的最优动力系统建模和动力学分析

研究了基于Navier-Stokes方程的脉动速度方程的最优低维动力系统建模理论.最优目标泛函为脉动速度基函数的不可压缩性和正交性.数值计算了充分发展的并排双方柱绕流问题,并基于双尺度全局最优化方法,建立了它的脉动速度的最优动力系统模型.对其相空间轨道、Poincaré截面、分岔特性、功率谱和Lyapunov指数集等动力学特性进行了分析.随着Reynolds数的增加,双方柱绕流的脉动速度方程最优动力系统具有复杂的类倍周期分岔行为.     

近壁湍流脉动的概率分布函数

采用大涡模拟方法,模拟槽道湍流,获得了不同雷诺数情形下的槽道流大涡模拟数据库．在此基础上,获得了流向和垂向脉动速度的概率分布函数,并运用假设检验,分析了其与正态分布的定量差别．进一步计算了流向和垂向脉动速度的偏斜度、平坦度,讨论了二者在粘性子层、过渡区和对数律区的变化．同时,讨论了粘性子层、过渡区和对数律区流向和垂向脉动速度概率分布函数的特点及其与湍流猝发的高速流下扫和低速流喷发事件的关系．最后,分析了雷诺数对流向、垂向脉动速度分布的影响．     

三维弹性管的涡致振动特性分析

采用有限体积法联合大涡模拟方法求解三维湍流流场,采用有限元法离散弹性管结构,对Re=1.35×104的湍流流动作用下三维弹性管的涡致振动进行了数值模拟,结构的动力学响应用Newmark算法来求解,管的运动采用基于扩散光顺方法的动网格模型来实现.利用建立的数值模型,分析了升力系数、阻力系数、位移、涡脱频率、相位差随频率比的变化特征,成功扑捉到锁定、相位开关,并联合运动轨迹、相图及Poincaré截面映射,研究了升力系数与横向位移的极限环与分叉等非线性特性.研究结果表明,在阻力系数的最小值处,横向振幅达到最大值,同时,横向响应的"锁定"也始于阻力系数最小值处;在"锁定"范围内,横向振幅随着频率比的增大而逐渐减小;在升力系数的最小值处,升力系数与位移间的相位由反相变为同相;在均匀湍流流动作用下,三维弹性管的升力与横向位移并未出现周期解的分叉.     

黏性分层流中椭球体俯仰振荡研究

以黏性密度分层流下椭球体自由俯仰振荡衰减过程为研究内容,建立了密度连续分层流数值计算模型.通过对经典小球黏性绕流场的数值模拟和增阻系数的计算验证了数值模型的正确性.以初始45°攻角下的椭球体俯仰振荡过程为研究对象,采用基于Aitken亚松弛适应算法的双向流固耦合方法,数值模拟了不同内Froude(弗汝德)数Fri下椭球体俯仰衰减振荡的动态绕流场.数值研究结果表明,俯仰振荡将上下搅动周围流体,在椭球体上下两侧对称形成四个密度涡环,密度的垂向分层限制了涡环的垂向传播,也加速了涡环的消失,这种限制助长了水平运动的发展,远场尾涡流场将以水平波动的形式传播.在较高的内Froude数Fri和Reynolds(雷诺)数Re下,双向耦合抑制了数值震荡.研究还发现,随着来流速度的增加,阻力系数不增反降,这说明,对于自由俯仰振荡的椭球体,负阻尼现象仍然会出现.     

前后排列旋转圆柱体对流热交换的格子Boltzmann方法解

用格子Boltzmann方法,数值研究流过前后排列两旋转圆柱体的二维层流．用二阶精度的速度场和温度场,数值化涉及运动的曲线边界．在Reynolds数为100,Prandtl数为0.71时,研究旋转速度比的变化和不同间距的影响．在4种不同间距(3, 1.5, 0.7, 0.2)下,研究旋转速度比的不同范围．结果表明,当间距取大数值时,第1个圆柱体的升力和阻力系数,与单个圆柱体相类似;对所有间距(除间距3以外),第2个圆柱体的升力系数,随着角速度的增加而减小,而阻力系数反而增加．圆柱体表面平均周期Nusselt数的结果表明,当两圆柱体间距小且角速度又低时,热传导是主要的传热机理,而当间距大且角速度又高时,对流是主要的传热机理．     

三角形二次插值系数有限元法解半线性椭圆问题的超收敛性

基于均匀三角形的剖分求解一类二阶半线性椭圆问题,用插值系数有限元方法比经典有限元法更容易实现,与经典二次有限元一样,二次插值系数有限元方法在对称点处也有四阶超收敛精度,数值计算表明这些结论是正确的.     

速度-涡量法数值求解具有表面吹吸圆柱的绕流问题

用速度-涡量法数值求解了具有表面吹吸圆柱的绕流问题.所得高阶隐式差分方程,采用以修正的不完全LU分解作预处理器的共轭梯度法(MILU-CG),高效解出.研究了雷诺数Re=100时,各种吹吸位置、吹吸强度对圆柱尾流涡旋结构和阻力、升力系数的影响规律.指出,在圆柱肩部的吸气和在圆柱尾部的吹气,可有效地抑制尾流涡旋结构在垂直来流方向上的非对称性,达到减小升力的目的.对在圆柱肩部吸气的情形,合适选择吸气强度,还可有效减小圆柱在来流方向上所受的阻力.     

蠕动流中圆球阻力系数的入口效应

本文研究了圆球在半无穷长圆管入口处的蠕动流。得到了速度分布,压力分布和流函数的无穷级数形式的分析解.采用配置法将无穷级数截断并确定出级数中各项系数.求出了均匀入口流绕静止圆球以及圆球以瞬时速度在管内静止流体中运动这两种情形下圆球的阻力系数以及圆球表面上的应力分布.结果表明,当圆球在入口处运动时会遭受到较无穷圆管内为大的阻力.本文还对配置法的收敛性进行了数值实验.试验证明,该法具有好的收敛性.     

Nonlinear characteristics of vortex-induced vibration at low Reynolds number

Numerical simulations of vortex-induced vibration of a two-dimensional elastic circular cylinder under the uniform flow are calculated when Reynolds number is 200. In order to achieve the vortex-induced vibration, two-dimensional incompressible Navier–Stokes equations are solved with the space–time finite element method, the equations of the cylinder motion are solved with the new explicit integral method and the emeshing is achieved by the spring analogy technology. Considering vortex-induced vibration with the low reduced damping parameters, the variety trends of the lift coefficient, the drag coefficient, the displacement of cylinder are analyzed under different oscillating frequencies of cylinder. The nonlinear phenomena of locked-in, beat and phaseswith are captured successfully. The limit cycle and bifurcation of lift coefficient and displacement are analyzed. Besides, the Poincare sections of the lift coefficient are used for discussing the bifurcation of periodic solution. There are some differences in nonlinear characteristics between the results of the one degree of freedom cylinder model and those of the two degrees of freedom cylinder model. The streamwise vibration has a certain effect on the lateral vibration.     

MILU-CG方法和绕旋转圆柱流动的数值研究

本文采用以修正的不完全LU分解作预处理器的共轭梯度法(MILU-CG),结合高阶隐式差分格式,改进了作者(1992)提出的基于区域分解、有限差分法与涡法杂交的数值方法(HDV).系统地研究了雷诺数Re=1000,200,旋转速度比α∈(0.5,3.25)范围内,绕旋转圆柱从突然起动到充分发展,长时间内尾流旋涡结构和阻力、升力系数的变化规律.计算所得流线与实验流场显示相比,完全吻合.首次揭示了临界状态时的旋涡结构特性,并指出最佳升阻比就在该状态附近得到.     

一个新的边界层方程及其在形状控制中的应用

本文给出固壁边界上(即一个二维流形上) 的流体速度梯度和压力的二阶偏微分方程, 从而也给出边界上法向应力, 以及流体中运动物体所受的阻力和升力的计算公式. 本方法的创新在于边界上法向速度梯度不是通过在边界层内速度梯度的数值微分达到, 而是通过它与其他变量一起作为一组偏微分方程的解而得到, 证明边界层方程组的适定性问题, 并且给出解关于边界形状的Gâteaux 导数所满足的偏微分方程. 本文将本方法应用于飞机外形的形状最优控制, 给出阻力泛函关于形状第一变分的可计算形式. 数值例子表明, 用本方法得到的阻力精度比通用程序得到要高.     

A second order characteristics finite element scheme for natural convection problems

In this paper a second order characteristics finite element scheme is applied to the numerical solution of natural convection problems. Firstly, after recalling the mathematical model, a second order time discretization of the material time derivative is introduced. Next, fully discretized schemes are proposed by using finite element methods. Numerical results for the two-dimensional problem of buoyancy-driven flow in a square cavity with differentially heated side walls are given and compared with a reference solution.     

CFD simulation of gas–solid flow with a cluster structure-dependent drag coefficient model in circulating fluidized beds

To model the effect of clusters on hydrodynamics of gas and particles phases in risers, the interfacial drag coefficient is taken into account in computational fluid dynamic simulations by means of a two-fluid model. The momentum and energy balances that characterize the clusters in the dense phase and dispersed particles in the dilute phase are described by the multi-scale resolution approach. The model of cluster structure-dependent (CSD) drag coefficient is proposed on the basis of the minimization of energy dissipation by heterogeneous drag (MEDHD) in the full range of Reynolds number. The model of CSD drag coefficient is then incorporated into the two-fluid model to simulate flow behavior of gas and particles in a riser. The distributions of volume fraction and velocity of particles are predicted. Simulated results are in agreement with experimental data published in the literature.     

On the accuracy of the viscous form in simulations of incompressible flow problems

We present a numerical study of drag/lift and flux estimates using two forms of Navier‐Stokes equations (NSE) that are equivalent in the continuum formulation but not in the discrete finite element formulation. The two investigated forms of the NSE differ in the viscous term that is represented in one form by νΔ u with ν being the viscosity and 2ν?·?^S u in the other form where ?^S represents the deformation tensor. The study consists of numerical analysis of the two forms and computations of drag/lift, pressure drop on the cylinder problem and computations of flux for the Poiseuille flow. The main objective is to provide a clear comparison of the reference values for the maximal drag and lift coefficient at the cylinder and for the pressure difference between the front and the back of the cylinder at the final time for the two forms of NSEs. Our computational results of the reference values do not differ significantly between the two forms, but the differences are there. For the Poiseuille flow, the differences in the flux computations were much smaller, and this agreed with the computationally obtained results of the divergence of the velocity field. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 28: 523–541, 2012     

Longitudinal permeability of spatially periodic rectangular arrays of circular cylinders II. An arbitrary distribution of cylinders inside the unit cell

We study the longitudinal permeability of unidirectional disjoint circular cylinders, when a Newtonian fluid is flowing at low Reynolds number along these cylinders; the longitudinal velocity satisfies the Poisson equation. The cylinders are arranged according to a doubly periodic structure. The number of cylinders in each rectangle can be arbitrary as well as their positions and radii. The method of functional equations yields analytical formulae for permeability in terms of these quantities. These formulae are written also in continuous form to study the flow for large numbers of cylinders. Special attention is paid to the case of the square unit cell, equal radii and lognormal distribution of radii.     

Numerical modeling of the pressure coefficient of the circular cylinder

The article examines the possibilities of numerical solution of chimney load from the effect of wind. The shear-stress transport (SST) k-ω turbulence model in ANSYS Fluent software is used to evaluate the task of the flow around the circumference of the rough cylinder. Calculations are performed on two different meshes that lead to the solution using wall function and near wall modeling. These two solution approaches in terms of defining wall roughness are presented in the paper by evaluating of the time dependence of the mean pressure coefficient distribution at the circumference, drag coefficient, and lift coefficient. The accuracy of the calculations is verified with parameters determined according to valid standards.