消除粘性项高阶离散数值振荡的半结点-结点交错方法

通常的粘性项处理方法是连续两次采用结点型中心差分格式求得.但是,通过两次采用高阶结点型中心格式求得的粘性项易在流场中的"间断"附近产生数值振荡.本文采用多种格式,通过求解一维Burgers方程来充分展示了这种振荡现象.为了消除这种数值振荡,一种半结点型的高阶紧致格式被用来求解粘性项.Fourier频谱分析表明,这种格式具有非常优越的保频谱性能.一维和二维的数值计算结果表明,通过"半结点-结点"交错采用这种半结点型格式的方法可以非常有效地避免粘性项高阶离散可能导致的数值振荡.     

双时间步长加权ENO-强紧致高分辨率格式及在叶轮机械非定常流动中的应用

在非结构网格下,针对非定常三维N-S方程组发展了一种双时间步长高精度快速迭代格式。该格式在时间上具有二阶精度,在空间上将r=3的加权ENO格式与强紧致格式相结合去处理N-S方程中的对流项以及离散方程的右端项,并用四阶精度的紧致格式去计算N-S方程中的粘性项。典型的3个算例从不同侧面对格式进行了考核。计算表明:该算法具有高效率与高分辨率的特征,所得的计算结果与相关实验数据比较吻合,初步表明了该算法可以在非结构网格下模拟非定常流动的物理过程。      

烧蚀外形方程差分计算方法研究

分析了烧蚀外形方程的性质,研究了人工粘性项对烧蚀外形计算结果的影响,发现人工粘性在抹平非物理波动的同时,会严重影响计算精度.借鉴于CFD中NND格式的思想,构造了求解烧蚀外形方程的无波动、无自由参数的耗散差分格式,得到了满意的计算结果.     

应用TVD格式分离求解不可压N-S流动

提出一种新的将连续方程与动量方程分离求解的TVD引用方法。对动量方程中的对流项采用三阶精度、迎风TVD格式离散, 粘性项和压力项采用中心差分格式离散。压力采用压力替代法求解。计算所用的网格是具有交错网格主要特点的修正非交错网格, 并采用一种压力修正法和G-S中心对称迭代法加速收敛。给出了二维管流和二维空穴流动在不同状态下的计算结果, 并与其它解法的计算结果或精确解进行比较。结果表明, TVD格式不仅能用于不可压流数值计算, 而且还能得到比非TVD差分法精度更高的结果。      

基于Van Leer＋AUSM混合格式求解三维N-S方程

基于多块粘性结构网格,开展了三维N-S方程数值算法的研究。控制方程的空间离散采用有限体积法,在前人工作的基础上,发展了Van Leer＋AUSM混合格式并应用到对流通量的离散中,粘性项采用中心格式离散并利用格林定理计算粘性通量中的导数项,时间推进采用五步R-K法,湍流模型为S-A一方程模型。最后,以M6机翼和某超声速弹丸的粘性流场作为数值算例,计算表明：发展的数值算法对跨声速、超声速流场均具有较高的分辨率,适用于跨、超声速流场的数值模拟。     

带副翼三维机翼绕流的N-S方程解

提出了一种求解带副翼偏转三维机翼绕流的N-S方程计算方法。采用对接分区网格与分区求解算法的结合,有效地求解绕此外形的复杂流动。提出了一种满足通量守恒的内边界耦合条件。数值方法中选用vanLeer分裂格式离散无粘通量项,并构造了一种Limiter函数以保证TVD性质,采用中心差分格式来离散粘性通量项。数值算例表明该方法是求解带操纵面偏转的机翼绕流的有效方法。     

高阶加权紧致非线性格式（WCNS）在二维流动计算中的加速收敛研究

以二维圆柱超声速无粘和粘性绕流的数值模拟为例,空间项采用高阶精度格式WCNS离散,对比研究了LU-SGS、高斯-赛德尔点松弛、线松弛以及GMRES等隐式求解方法的收敛性,并对右端项中的无粘通矢量、GMRES方法中的预处理和子迭代等影响作了对比计算.结果表明,右端项采用Steger-Warming无粘通矢量的收敛性优于其他通失量方法,采用了准确的解析雅克比矩阵的点、线松弛的收敛速度优于LU-SGS,以线松弛为预处理的GMRES算法具有良好的收敛特性.     

跨声速轴流压气机特性全三维粘性流动分析

应用Steger-Ｗarming通量分裂技术,将守恒型方程中的流通向量分裂成两部分,应用 NND格式对无粘项进行数值离散,粘性项采用中心差分格式,采用LU-SGS隐式推进迭代建立了跨声速轴流压气机性能及内流场三维N-S方程的高分辨率和高效率的数值方法。数值模拟给出NASA ROTOR37和NASA ROTOR67性能及内流场流动结构。并进一步基于混合平面法建立了多级轴流压气机性能计算模型,详细分析了某一级跨声速轴流压气机性能和内流特性,计算与实验结果比较表明了本模型与方法是可行的和可靠的。     

求解守恒律及对流扩散方程的中心迎风方法

提出了一种新的求解双曲守恒律方程(组)的四阶半离散中心迎风差分方法。空间导数项的离散采用四阶CWENO的构造方法,使所得到的新方法在提高精度的同时,具有更高的分辨率。使用该方法产生的数值粘性同阶要比交错的中心格式小,而且由于数值粘性与时间步长无关,从而时间步长可根据稳定性需要尽可能的小。     

扩压器内二维湍流流动的时间推进数值解

ＡＦ有限差分理论。指出薄层Ｎ－Ｓ方程抛物化后右端粘性项是主元素变量矢的一个简单线性函数,这对实际计算很有用。应用隐式ＡＦ方法计算了为发动机进气道设计的典型扩压器内的复杂湍流流场。研究表明,计算结果与实验结果基本一致,这为三维非定常湍流流动进一步研究提供了基础。     

TVD格式计算一维Euler方程间断解存在的几个理论问题

本文讨论了TVD格式按本文第二节A段意义推广到拟线性方程组的几个理论问题。指出:1.对于Riemann间断分解,Euler方程的解和Riemann不变量都没有TVD性质,因此不能将TVD格式自然推广;2.TVD格式是提高分辨率的一种途径,它不能保证得到物理解,能得到物理解的TVD格式一定是耗散守恒格式,因此将TVD格式推广时,必须保证格式是耗散的,否则仍会得到非物理解。并给出了一阶迎风格式,TVD2,UNO2得到的若干非物理解算例;3.非物理解现象只是发生在具有剧烈膨胀的单波区的计算上,而某些TVD格式若会得到非物理解,则一定发生在特征值变号的单波区的计算上。     

用虚拟黏性法构造Navier-Stokes方程黏性项的隐式格式

提出的虚拟黏性法是一种关于Navier-Stokes方程黏性项(黏性应力项和传热项)的隐式新方法。通过引入虚拟时间和虚拟黏性项,将隐式格式的构造大大简化,从而避免了大型复杂隐式差分方程组的常规求解。在虚拟时间推进过程中所需求解的方程组的系数矩阵是一个三对角矩阵,它具有计算简单且计算量小的优点。用模型方程和Navier-Stokes方程进行了数值仿真,研究各种参数对计算的影响,并在精度和效率上与显式方法进行比较,证实了算法的正确性和优势。最后分析了本文格式的适用范围。     

钝锥超声速粘性绕流的隐式推进求解方法研究

本文提出了一种新的推进求解方法。通过对模型方程的稳定性分析和钝锥超声速绕流的数值模拟,表明该方法既克服了文献中推进格式步长有下界和稳定性差的不足,也克服了亚声速区压力不好处理的困难。文中建立的一阶、二阶隐式推进格式,稳定性好,推进步长无下界,计算结果的精度较高。     

不安抗辩权与预期违约之比较

《合同法》在第68条和第69条规定了不安抗辩权制度．在第108条规定了明示预期违约。但《合同法》对英美法系的预期违约制度和大陆法系的不安抗辩权制度融合的不够彻底。使预期违约和不安抗辩权在适用上产生了混乱。提出了通过对两大法系相关制度的比较来完善我国合同立法的思路。     

高精度的对流差分格式

对于广义对流差分格式,本文应用两种方法(截断误差分析法及事后估计误差法)导出了使其误差最小的最优α值。和以前的方法不同,它们是从所研究的对流-扩散方程本身出发导得的,而无需事先知道方程的解析解。大量算例结果表明,截断误差分析法有很高的精度。事后估计误差法的误差,除个别点外,也都比一般格式的误差要小或者小得多。但是,对于复杂的问题,前者的应用可能会有困难,而后者却很易实现。     

乘方和 QUICK 格式在非对称扩压器内流动计算的比较

采用Ｓｐｅｚｉａｌｅ的非线性ｋ－ε模型计算了非对称扩压器内的流动,着重用数值方法考察了两种典型的对流差分格式（乘方格式和ＱＵＩＣＫ格式）对模型预测性能的影响,得出两种紊流模型对对流差分格式并不敏感的结论。     

旋转体跨音速零升波阻和压力的差分计算

应用Murman-Cole的有限差分法,求解具有纵向大扰动而横向小扰动的跨音速轴对称速势方程,由此计算旋转体跨音速零升力时的压力分布和波阻力,以及激波位置。物面边界条件被转移到物体轴上。远场边界条件由无穷远处的条件近似代替。计算物面压力系数时,用细长体理论进行物面速势插值。 速势的差分方程用沿半径方向线超松弛改进迭代求解。网格取62×16,迭代初场取零,达到收敛的迭代次数对M_∞<1,M_∞>1以及M_∞≈1分别大约为150,40和300次。松弛因子取为:M_∞<1时,0.9≤ω_b≤1.7,0.9≤ω_p<1.0;M_∞≥1时,0.8≤ω_b≤0.9,0.8≤ω_p≤0.9,这里ω_b,ω_p分别为局部亚音速点和超音速点的松弛因子。 算例为七种不同外形的细长体,计算结果与实验符合尚好。 文中对网格、初场、迭代方法、松弛因子等有关收敛性、收敛速度问题进行了探讨。在局部线化条件下,对定常小扰动轴对称势流的差分方程,进行了线超松弛迭代的稳定性和收敛性分析。数值计算经验与理论分析所得结论相符。     

COMPACT FINITE VOLUME SCHEMES AND THEIR APPLICATIONS

ＣＯＭＰＡＣＴＦＩＮＩＴＥＶＯＬＵＭＥＳＣＨＥＭＥＳＡＮＤＴＨＥＩＲＡＰＰＬＩＣＡＴＩＯＮＳＪｕＨｏｎｇｂｉｎ（ＳｃｈｏｏｌｏｆＰｏｗｅｒａｎｄＥｎｅｒｇｙＥｎｇｉｎｅｅｒｉｎｇ，ＳｈａｎｇｈａｉＪｉａｏｔｏｎｇＵｎｉｖｅｒｓｉｔｙ，Ｓｈａｎｇｈａｉ...     

串件拼修对策下两级备件维修供应系统动态管理模型

 在备件维修供应系统日常运作中,通过对库存和维修资源分配进行动态管理,能够在现有库存配置方案下提高保障效能。针对串件拼修对策下两级维修供应系统,分别建立了有限维修资源约束下多维修渠道故障件的送修和库存分配动态管理模型。根据备件库存动态分析模型(Dyna-METRIC)原理,建立了保障系统的蒙特卡罗仿真模型,并将动态管理模型嵌入其中。使用传统先到先服务(FCFS)原则和指定可用度目标确定不同初始库存配置方案,将这些方案作为仿真模型输入进行仿真试验。结果表明,采用本文模型后保障系统效能相比传统FCFS原则有了明显提高,证明了模型的正确性。     

槽道流动中湍斑的数值模拟

采用Navier-Stokes方程直接数值模拟槽道流动中湍斑的形成和发展.建立高精度、高分辨率的三维耦合中心差分格式,并和时间分裂法、高阶迎风差分格式一起组成高效的高阶紧致差分方法;克服一般高阶中心差分格式不适用于边界邻域的困难,可用于包括边界邻域在内的整个计算区域;在边界上满足Navier-Stokes方程的无反射和压力条件,保证了格式的精度和一致性.通过模拟槽道流动中从壁面脉冲到湍斑的演化过程,显示出湍斑内雷诺应力和速度、压力脉动等湍流的基本特征,表明了湍斑的生成与壁面脉冲的强度和Re数大小直接相关.