不可压流动的二阶迭代投影法

投影方法的计算,需要引入一个中间速度场及相应的人工边界条件,因而带来了分裂数值误差和数值边界层.为了解决这些不足,提出了一种迭代投影方法.在每个时间步,采用投影方法作为该方法的子迭代过程,当迭代收敛时便构成了完全耦合的数值方法.既然中间速度场是对真实速度场的逐步逼近,因此,就无需人工边界条件,或者说人工边界条件即为物理边界条件.数值试验表明: 迭代投影方法可以显著地减小数值边界层; 经过1～4次迭代后,速度和压力在时间方向上都可以达到二阶精度.     

求解不可压Navier-Stokes方程的三阶精度投影方法

在数值模拟不可压缩流动中,投影方法得到了越来越广泛的应用,一般认为,目前的投影方法在时间方向上仅仅停留在二阶精度.该文通过分析发现,压力更新公式在高阶投影方法的构造中具有非常重要的作用,它不仅影响压力的精度,同时还影响格式的稳定性.在此基础上结合相容的压力更新公式,提出了三阶精度的投影方法.正则模态分析的结果表明,该文提出的投影方法速度具有三阶精度.当采用相容的压力更新公式时,压力也具有三阶精度,而且格式是稳定的; 而采用传统压力更新公式,压力在边界附近只有二阶精度,且格式不稳定.数值结果验证了上述结论.     

基于反欠松弛方法的不可压流动压力修正算法

对交错网格上不可压流动的压力修正算法进行了研究，利用离散的动量方程和连续方程建立了压力方程，通过检查迭代中速度场的散度来建立压力修正方程和速度修正公式。提出了加速收敛和稳定数值计算的反欠松弛方程，给出了决定反欠松弛函数的准则。数值计算获得了较好的收敛特性。     

不可压粘性流动N-S耦合方程的有限元解法—I应用

应用有限元方法求解不可压粘性流动Ｎ－Ｓ耦合方程，将此方法用于驱动方腔和突扩管道内的层流分析，结果表明此方法具有良好的椭圆型性质，能预测分离流动．用该方法结合两层湍流代数模型，对压气机叶栅内不可压流动进行模拟分析，与实验结果的比较证明方法是有效的．     

求解不可压三维湍流的隐式SMAC方法

该文基于SM AC(s im p lified m arker and ce ll)方法,发展了一种在任意曲线坐标系中求解三维粘性不可压湍流R eyno lds时均方程的全隐式数值方法。基本方程是以逆变速度为变量的R eyno lds时均动量方程和椭圆型压力Po isson方程,并采用标准k-ε湍流模型封闭方程组。压力Po isson方程用T schebyscheff SLOR方法交替方向迭代求解。R eyno lds时均动量方程、k方程和ε方程对流项均采用Chakravarthy-O sher TVD格式离散,该格式不但有助于提高数值稳定性,而且能有效消除网格扭曲较大的地方产生的非物理振荡误差。用自编程序对后台阶方腔流场进行了计算,计算结果和实验结果吻合较好。     

不可压Navier-Stokes方程的扩散松弛近似(英文)

通过扩散松弛机制,研究了二维环T2上不可压Navier-Stokes方程一个双曲奇异摄动问题.运用可对称化双曲系统的能量方法,严格证明了Navier-Stokes方程的扩散松弛极限.     

管内四孔通道后不可压湍流流动的LDV测量

运用偏振分离型ＬＤＶ系统对流经管内四孔通道的不可压湍流流动过程进行了测量研究。结果表明，通孔下游附近，流动具有湍动射流特征；管壁附近存在明显的回流再附区、在轴向时均速度分布趋于稳定前后，轴向湍流度由大于切向湍流度变化至两者近似相等。     

针对不可展曲面的近似展开方法浅析

随着市场发展和人们生活质量的提升,现代工业产品造型越来越复杂多样。复杂曲面大多通过平面材料加工而成,对不可展曲面的近似展开已成为工业设计的重要课题。该文对近20年来常用的不可展曲面的近似展开方法进行了分析、总结和归纳。     

基于投影收缩的SA方法求解随机变分不等式问题

求解变分不等式的各种算法中,投影收缩算法易于执行、稳健、而且可以处理大规模问题,因此发展迅速.何炳生教授根据变分不等式及投影算子的性质确定的三个不等式,提出了求解变分不等式的投影收缩算法,此方法简单易行,且便于实现.用随机近似方法来求解随机变分不等式和随机优化问题已经被广泛的研究,其中函数值和一阶导数不可求,但可以用近似的方法得到.将投影收缩算法应用到求解随机变分不等式当中,在一些适当的条件下,可得到全局收敛的结果.     

用改进的强隐式方法求解三维不可压粘性流动问题

提出了一种用改进的强隐式方法求解三维不可压缩粘性流动的交错网格有限体积法．该方法在ＳＩＭＰＥＲ算法的基础上，对Ｎ－Ｓ方程的对流项采用ＱＵＩＣＫ格式离散，利用ＭＳＩ方法求解代数方程．数值计算表明本方法具有计算稳定、收敛速度快等特点．     

求解不可压缩流动的分步有限元格式

提出求解不可压缩 Navier- Stokes方程的分步有限元格式 ,该格式没有高阶微分项产生、程序编制简单 ,适用于非线性的多维复杂流动。应用该方法实际模拟了二维圆柱绕流的旋涡形成与脱落过程 ,得出了不同 Re情况下圆柱绕流的流速分布。计算得到的不同 Re下的旋涡脱落频率(Strouhal数 )与前人已有的经典解答符合良好     

Bessel函数近似解法比较

利用差分方法求Bessel方程的数值解;并与Bessel函数的2种级数表达式的计算结果进行比较,通过数值计算及相应的分析得出,Bessel方程的数值解与级数解各有利弊;利用计算机,由差分方法计算,无需考虑自变量的大小,但其解是离散的,不可进行解析运算;而由2种级数表达式计算时,必须关注自变量的大小,但可以利用它们进行解析运算.     

A projection method for LES of incompressible turbulent combustion

~~A projection method for LES of incompressible turbulent combustion~~     

用分步有限元法求解三维不可压缩流动

将计算二维不可压缩流动的分步有限元格式扩展到三维情况 ,由于该格式没有引入新的高阶空间导数项 ,适用于多维空间的非线性问题。实际求解了三维空腔流 ,在 Re数较低的情况下得到的流态不随时间而变化 ;在高 Re数(Re=40 0 0 )的情况下 ,流态十分复杂并且是非定常的 ,在水平断面具有多个非定常旋涡 ,在垂直于流动方向的立面断面上可以模拟到 TGL涡。本格式的迎风效应是 Taylor展开式的高阶精度项 ,没有人工粘性引入 ,得到的结果可靠性好     

粘性不可压流场数值模拟的无网格方法

采用无网格伽辽金法(elementfreeGalerkinmethod,EFG)对二维不可压粘性流场进行了数值模拟,主要计算了几个经典的流体力学问题:管道流、绕台阶流以及方腔驱动流。计算中采用标准的Galerkin离散,使用规则分割产生的背景网格来形成积分单元,并与有限元、有限差分法计算所得的速度剖面图进行了比较,结果相差很小。更重要的是,与有限元法相比,无网格虽然计算量比较大,但是它能够处理结点分布畸形甚至重合的情况,从而为解决多相流动直接数值模拟中网格重构的问题奠定基础。     

求解凸集约束优化问题的共轭梯度的GLP投影算法

利用GLP投影技术 ,对凸约束的非线性规划问题构造了一个共轭梯度的GLP投影算法 ,在一维精确步长搜索下 ,给出了算法较强的全局收敛性结果 ,由于算法需要较小的存储量 ,特别适合于计算大规模的约束优化问题。该算法提高了梯度投影法的收敛速度。