方柱绕流的大涡模拟

采用大涡模拟(LES)方法对雷诺数为2.2×104的方柱绕流流场进行了数值模拟.使用非交错网格的有限差分法,分别对准三维物理模型和真实三维物理模型求解不可压N-S方程,将沿流向方向方柱水平中心线上的时均速度的计算结果与实验数据进行了比较,结果表明,三维模型的模拟结果优于准三维模型的模拟结果.比较升力系数和阻力系数发现,与二维模拟(RNG)方法相比, 三维模拟的结果更加接近实验测试数据.     

方柱绕流的大涡模拟

采用物理分步法对方柱绕流的流动分离现象进行了大涡模拟(LES).物理分步法将N-S方程剖分为对流、扩散及波动过程,对不同的物理过程分别采用适宜的方法进行求解.对流过程采用逆风——中心差分格式求解;扩散过程采用中心差分格式求解;压力方程采用SOR迭代法求解.计算的方柱四周的平均压力系数的分布与物理模型试验结果吻合良好;模拟的流动现象也得到了物理模型试验的证实.     

二维方柱绕流阻塞效应的大涡模拟

在均匀流条件下对阻塞比(R)分别为4.2%、10%、16.7%、20%和25%的二维方柱绕流进行了大涡模拟(LES)研究,对比了各工况下柱体的整体气动力和表面风压分布随阻塞比的变化规律,分析了阻塞比对柱体附近区域内流场特性和瞬时涡量的影响.研究表明,当阻塞比由4.2%增大到25%,柱体的平均阻力系数、脉动升力系数和Strouhal数分别增加了27.6%、43.6%和31%.随着阻塞比的增大,侧面及背风面负风压系数的均值绝对值以及脉动值均显著增加,柱体两侧附近流场特性也有明显的变化.在大阻塞比工况中,洞壁附近会生成二次旋涡,并与柱体脱落的旋涡相互掺混.     

二维翼型绕流计算中预处理和多重网格方法的应用

结合Turkel矩阵预处理方法和多重网格方法,发展一种适合低速粘性流动计算的高效数值方法,通过对Navier-Stokes方程的时间导数项实施Turkel矩阵预处理,使得可压缩Navier-Stokes方程在低速情况下的系统刚性得到改善,为进一步加速收敛,提高计算效率,采用多重网格的3层V循环方式,对RAE2822超临...     

高雷诺数条件下二维方柱涡激振动的数值模拟

利用作者编制的基于松耦合方法的求解气动弹性问题的数值模拟程序,计算了二维方柱在高雷诺数条件下,的横向涡激振动.结果表明:随着折减风速的增加,在特定折减风速条件下,方柱的涡激振动出现锁定现象;由于方柱的自然振动频率和涡激振动系统的旋涡脱落频率的相互作用,涡激振动出现拍的现象;随着折减风速的增加,涡激振动系统的旋涡脱落模式从异相2S (2-single)变化到同相2S,即相位转换现象.该文还将计算结果和文献中的计算结果进行了比较和评价.     

二维翼型绕流的无网格差分解法

 根据作者提出的基于全方向导数的无网格方法,运用方向全导数公式构造样点的所有偏导数,并用该公式对第二类边界条件进行离散处理.它不需要任何网格或单元,所以是非常简便彻底的无网格方法.数值算例表明,该方法具有很高的精度.     

二维可压翼型绕流的新型无网格解法

根据一种新型的基于方向全导数的无网格方法,运用方向全导数公式构造样点的相应偏导数,并用该公式对第二类边界条件进行离散处理,这样不需要构造任何网格或单元,就可以对多种流动问题进行求解,所以是非常简便和彻底的无网格方法.通过求解两个可压缩流场中二维翼型的数值算例表明,该方法对于解决可压流场中的问题具有较高的精度.     

二维非定常有涡流动的数值模拟方法

为有效地模拟二维有分离现象的粘性流动,发展了求解涡量流函数方程的数值方法。对涡量输运方程的时间项运用四阶Runge-Kutta法离散,从而将方程分解为4个计算步分别求解。在每一计算步对方程进行Steger-Warming近似因式分解,从而使涡量的计算在空间的两个方向上分别进行。对流项采用Chakravaythy-Oscher总变差减少(TVD)格式离散。流函数的Poisson方程运用Tschebyscheff SLOR方法交替方向迭代求解。运用该方法对突然起动圆柱和高雷诺数时弯曲薄翼的非定常有涡流动进行了数值模拟,并将计算结果与其它方法的计算结果和试验结果进行了对比,结果表明本方法具有精度高、数值稳定性好和计算效率高的优点。     

二维非定常对流扩散方程的对角占优隐格式及多重网格算法

提出了一种数值求解二维非定常变系数对流扩散方程的对角占优、空间为二阶精度的隐格式。利用Fourier分析方法证明了该格式是无条件稳定的．由于格式具有对角占优性,因此适用于大梯度（高雷诺数）问题的数值求解．另一方面,为了克服传统迭代法在求解隐格式时收敛速度慢的缺陷,采用了多重网格加速技术,大大加快了迭代收敛速度,提高了求解效率．数值实验结果证明了该格式的精确性、稳定性和对高网格雷诺数问题的强适应性．     

运动平板附近圆柱绕流的数值模拟

用格子Boltzmann方法模拟运动平板附近的圆柱绕流问题, 给出一种精确确定临界间隙率的综合判定方法, 并分析了流场特性的内在本质以及各种物理现象之间的联系. 将圆柱置于运动平板上方, 平板运动速度与入口处均匀来流的速度保持一致, 模拟的雷诺数为1 000. 定义间隙率为G/D, 其中G为圆柱边界到运动平板的最小距离, D为圆柱的直径. 结果表明: 当间隙率取值范围不同时, 流场特性有较大差异; 与孤立圆柱相比, 本文中的升力和阻力有明显增加, 并且旋涡脱落也受到平板抑制.     

Numerical Simulation of Uniform Flow past an In-Line Oscillating Circular Cylinder

IntroductionOscillation of cylinders in the in- line or cross- flowdirection relative to an incident flow is a commonoccurrence in industry. Understanding vortex-induced vibrations is of great importance in thedesign of a variety of offshore engineeringstructures and nuclear plant components,forexample,heat- exchangers. To understand thebehavior of this complex interactive vibrationsystem,a number of researchers have examinedthe somewhat simpler and more tractable problemof circular cylinder f…     

用格子Boltzmann方法计算剪切流的方柱绕流问题

用格子Boltzmann方法计算来流为水平剪切流的方柱绕流问题, 得到了在不同速度梯度条件下方柱绕流的流线和等涡线图. 发现在圆柱尾部形成两排涡, 当来流速度梯度较大时, 两排涡有很大的不同. 计算结果表明, 用格子Boltzmann方法计算剪切流的方柱绕流问题是可行的, 计算结果与理论分析相符.     

数值模拟放置附属圆柱的主圆柱绕流

利用格子Boltzmann方法, 对主圆柱尾流区域内放置附属圆柱的绕流进行数值模拟. 结果表明: 放置单个附属圆柱时, 主圆柱所受阻力减小, 阻力
系数发生周期性改变; 放置两个附属圆柱的减阻效果更好.     

隔水管附属管控制流动的离散涡模拟分析

采用离散涡方法对单隔水管和带附属管的隔水管的二维绕流问题分别进行了数值模拟分析.首先对4种不同附属管管数的配置情况进行了模拟,结果表明,除了管数为4、来流角度为45°的情况外,附属管能在一定程度上降低升阻力,并抑制漩涡脱落,在附属管管数较多的情况下,抑制效果更好;然后对附属管管数分别为4和10的2种布局,进行了来流方向的敏感性分析并与实验结果对比,发现多管数对来流方向的敏感性要低得多,相对于单隔水管,10根管能使时均阻力系数和脉动升力系数分别获得约26%和79%的降低.采用离散涡方法,不仅计算时间短,而且模拟结果比较可靠,具有较好的参考价值,丰富了控制隔水管绕流流动的研究手段.     

非均匀格子Boltzmann方法模拟方柱绕流

应用非均匀格子Boltzmann方法对方柱绕流的三种情况进行详细数值模拟，在第1种情况中，方柱位于流场中央，模拟了卡门涡街现象，给出了斯特鲁哈数随雷诺数变化曲线；在第2种情况中，方柱位于流场壁面，分析了雷诺数对方柱后回流区的影响；在第3种情况中，两方柱并列在流场中央，考察了方柱间距对流场的影响。     

改进镶边法构造任意奇阶同心幻方

首先引入了双关联等差数列的概念,借此提出了一个构造n阶幻方的充分条件,然后将奇阶幻方分为n=4m-1阶与n=4m+1(m=1,2,…,m∈N)阶两类,介绍了一种改进的镶边法,分别构造两类奇阶幻方,并给出了严格的证明.此构造法简单易行,灵活多变,所构造出的幻方具有独特的性质.     

计算广义逆A_(T,S)~((2))的初等变换法

利用矩阵秩的性质,讨论了一类加边矩阵M的逆,给出了逆矩阵M-1的一种表示．由计算逆矩阵的初等变换方法,给出了计算矩阵A的广逆AT,S的初等变换方法．     

Predication of Plastic Flow Characteristics in Ferrite/Pearlite Steel Using a Fern Unit Cell Method

The flow stress of ferrite/pearlite steel under uni-axial tension was simulated with finite element method (FEM) by applying commercial software MARC/MENTAT. Flow stress curves of ferrite/pearlite steels were calculated based on unit cell model. The effects of volume fraction, distribution and the aspect ratio of pearlite on tensile properties have been investigated.