方柱绕流的大涡模拟

采用大涡模拟(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,即相位转换现象.该文还将计算结果和文献中的计算结果进行了比较和评价.     

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

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

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

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

二维浅水流动的离散涡方法

根据二维浅水流动中旋涡既可产生于流动边界层,也可在主流场中自生自灭的特点,建立起Eulerian-Lagrangian框架下的二维浅水流动离散涡方法;利用该方法数值仿真了高雷诺数下的二维浅水方柱绕流,得到了流动的速度场、涡谱,计算结果中明显可见流动分离、旋涡以及旋涡随时间的发展演化,即模拟出流动所具有的非定常不稳定等非线性特征。     

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

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

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

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

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方法对方柱绕流的三种情况进行详细数值模拟，在第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.     

用双硫腙分光光度法同时测定Zn~(2 )、Cd~(2 )和Pb~(2 )

介绍 CPA法测定多种组分的原理 ,并应用该原理 ,对用双硫腙分光光度法同时测定混合样品中 Zn2 、Cd2 和Pb2 的质量浓度进行了研究 ,取得了较为有意义的结果     

笛卡尔切割单元法在湖泊二维流场模拟中的应用

针对湖泊几何边界复杂的特点,为了高效并准确地计算流场,采用了基于有限体积法的笛卡尔切割单元法.该方法在流场主要区域采用笛卡尔网格,而对与固体边界相交处的单元采用特殊的处理.这样既保留了笛卡尔网格求解速度快的优点,同时在边界处又保证了精度,避免了传统差分法处理该类问题需采用阶梯形网格逼近而导致的失真问题.通过与一个标准算例的贴体网格计算结果对比,显示计算结果一致.最后运用笛卡尔切割单元法,求解了杭州西湖的二维流场,结果表明该方法适用于具有复杂几何边界的湖泊流场求解,且精度较高.     

用两个正交拉丁幻方构造2n+1阶完美幻方的一种简便方法

先构造两个2n 1阶正交拉丁幻方,再经一系列列变换得到另外两个正交拉丁幻方,进而构造出2n 1阶完美幻方．