多级轴流压气机叶栅内三维紊流流场的数值模拟

对多级轴流压气机叶栅内三维紊流流场进行数值模拟,采用高精度高分辨率的三阶ENN格式以保证对激波的捕捉和对紊流特征的正确模拟,利用LU-SGS隐式解法提高了计算速度,从而构成了一种既准确又高效的多级跨声速轴流压气机紊流流动数值求解系统。重点研究了动静网格的交接方法及相应的动静交接面处理模型。计算了某带进口导叶的三级轴流压气机,计算结果与实验数据吻合较好。     

高阶Schr(o)dinger型方程的两层高精度恒稳差分格式

众所周知, 高阶Schrodinger方程在量子力学、非线性光学及流体力学中都有广泛的应用.本文对高阶Schrodinger型方程(eu/et)＝I(-1)m(e2mu)/(ex2m)(其中I＝-1,m为正整数),利用待定系数法,构造出一个两层高精度的隐式差分格式.其截断误差阶为O((Δt)2+(Δx)6),比同类格式精度高2～4阶,并用Fourier分析法证明了它是绝对稳定的.最后,数值例子表明本文格式比著名的Crank-Nicolson格式精度高10-2～10-7,这说明我们的格式是有效的,理论分析与实际计算相吻合.     

推进－迭代算法计算效率的研究

本文对求解三维定常超音速动性流场的一次空间推进，在每一个推进站沿伪时间层局部迭代的推进－迭代算法作了进一步的研究．在每一推进站（侧向平面）沿伪时间层局部迭代时，给出了四种不同的隐式迭代方法，即沿侧面两个方向（法向和周向）全用隐式；法向隐式而周向采用Ｇａｕｓｓ－Ｓｉｌｄｌｅ来回扫描迭代；法向隐式而周向显式及以系数矩阵谱半径代替系数矩阵的简化标量隐式算法．用这四种算法模拟了三维球锥黏性绕流，给出了四种不同算法的计算效率和收敛特性比较．     

混凝土本构模型及其动态有限元算法研究

混凝土的力学性质及其本构研究是当前力学领域的一个重要课题。混凝土的实验数据表明其力学性质表现出很强的各向异性及并且在高加载速率的情况下，混凝土的性质物美价廉 时时明显不同。这就要求我们采用应变率相关的弹粘弹性本构模型来描述混凝土的力学行为。本语文从Ｏｔｔｏｓｅｎ的四参数混凝土模型出发，引进损伤和应变率的影响，参考关联塑性的理论，改进了混凝土的本构模型。为了把混凝土的数值模拟推向实用，本文概括混凝土     

分块隐式有限差分法计算弯管紊流

本文利用贴体坐标系中分块隐式有限差分法计算矩形截面90°弯管中不可压恒定紊流.在计算中,雷诺方程的数值离散采用混合差分格式,局部联立求解雷诺方程和连续方程而得到速度压力解.在全流场的迭求解过程中采用对称联立Gauss—Seidel法.利用标准K-ε紊流模型模拟紊流.计算结果与有关试验进行了对比.     

一种改进的分块隐式数值方法

提出了一种改进的分块隐式数值方法，在贴体坐标和交错网格下以逆变速度分量和压力Ｕ，Ｖ，Ｗ，ｐ为基本求解变量，由此克服了原分块隐式数值方法求解复杂边界流动时的困难．９０°弯管流动数值计算初步表明，本文提出的方法合理、可行     

Three-stage Stiffly Accurate Runge-Kutta Methods for Stiff Stochastic Differential Equations

In this paper we discuss diagonally implicit and semi-implicit methods based on the three-stage stiffly accurate Runge-Kutta methods for solving Stratonovich stochastic differential equations（SDEs）.Two methods,a three-stage stiffly accurate semi-implicit（SASI3） method and a three-stage stiffly accurate diagonally implicit （SADI3） method,are constructed in this paper.In particular,the truncated random variable is used in the implicit method.The stability properties and numerical results show the effectiveness of these methods in the pathwise approximation of stiff SDEs.     

基于混合网格Navier-Stokes方程的并行隐式计算方法研究

针对结构网格很难处理复杂外形和非结构网格无法计算具有边界层的粘性流动的缺点,发展了基于混合网格格点的隐式算法,成功地解决了在工程应用中难于处理的复杂外形粘性流场计算和效率问题。同时针对大规模的工程问题,发展了基于MPI通信技术的染色分层通讯并行计算方法。其中空间离散采用基于Roe格式发展的三阶迎风HLLEW(Harten-Lax-Van Leer-Einfeldt-Wada)或AUSM格式,湍流模型采用k??两方程湍流模型,时间推进考虑到LU-SGS并行等效较困难则采用基于DP-LUR(Data-Parallel Lower-Upper Relaxation)格式的隐式算法,计算CFL数可取到105量级,从2个到128个CPU的并行加速效率都保持在90%以上,大大提高了计算效率。算例对标模M6机翼模型流场进行计算,验证了方法的可靠性;然后对标模DLR-F6翼身组合体进行混合网格粘性与无粘计算结果进行比较,进一步验证混合网格方法;最后计算了DLR-WBNP外挂发动机翼身组合体模型,准确模拟了外挂和超临界机翼的相互干扰流动问题,采用4 CPU 16 CORE到24 CPU 96 CORE,2000步计算时间都不超...     

基于隐马尔可夫模型的地磁匹配算法

为了在地磁特征微弱区域内实现地磁辅助导航以及提高惯性导航系统在地磁特征明显区域内的定位精度,提出了基于隐马尔可夫模型的地磁匹配算法.首先以惯性导航系统定位误差为隐状态,以实时测量的地磁强度为观测量,建立了地磁匹配的隐马尔可夫模型;其次,针对该模型,使用Viterbi算法来确定最优状态序列,给出了惯性导航系统的当前定位误差.仿真结果表明,该算法可以实现地磁辅助导航,导航误差优于EKF算法,组合导航系统的定位误差在50 m左右.     

基于移动粒子半隐式法的表面张力模拟

采用移动粒子半隐式法(MPS)模拟了受表面张力影响的自由面流动。表面张力的计算采取了一种较适合于MPS方法的表面自由能模型。方形液滴振荡和射流断裂的模拟结果分别与理论分析和试验结果一致,同时进行了三维射流注水模拟,从而验证了MPS方法结合该表面张力模型可以有效、方便地进行自由面流动中表面张力现象的模拟。