数值解法在燃烧室出口燃气温度计算中的应用

为了求解燃气温度,利用双变量迭代法和赫夫法求解根据物质守恒、化学平衡、压力和能量守恒方程建立的非线性方程组,通过算例验证了2种数值解法求解化学平衡成分和温度的计算结果基本一致,在燃烧室出口温度场试验中,2种数值解法都能很好地满足测试需求,但双变量迭代法的计算收敛更快。试验结果表明:燃气分析的温度计算结果与热电偶的测温结果沿着燃烧室周向变化的趋势基本一致,燃气分析计算的温度基本大于热电偶测试结果。     

迭代动力缩聚法的收敛性分析

 利用Lyapunov矩阵方程和Riccati矩阵方程解的理论,对迭代动力缩聚法的收敛性进行了分析证明,并给出了迭代收敛的充分条件。揭示了动力缩聚法与经典的子空间迭代法的内在关系,阐明了各自的优缺点。迭代动力缩聚法实质上是子空间迭代法的变形,它需要人为选择主辅自由度,而子空间迭代法需要人为选定初始迭代向量。从理论上讲,只有主辅自由度选择满足收敛的充分条件要求,才能保证迭代结果收敛到理论上的精确解。给出了一个数值算例,对几种算法进行了对比,并验证了本文的论点。     

非对称纤维增强板室温形状计算的一种数值解

在研究非对称层合薄板室温形状时,考虑大挠度变形,应用最小势能原理,最后得到一组非线性代数方程组。Hyer在求解时采用的是拟牛顿法,迭代运算常常不收敛。为了解决这一问题,在对Hyer的正交层合板问题解加以推广的同时,对计算方法也作了一些探讨。 我们所需求解的问题有其特殊性,因此,仅单独使用增量法或迭代法进行数值计算往往难以达到计算的目的。为此,我们把增量法和迭代法联合使用,以达到迭代收敛,求解准确的目的。计算结果表明;这种方法能满足求解上述非线性方程组的要求,得到满意的数值计算方法。     

线性规划问题的一个数值解法

通过用一组非线性方程组来逼近线性规划的Kuhn-mcker条件方程，从而将求解线性规划问题转化为一个用Newton迭代法来求解非线性方程组的问题。经数值计算表明，该算法在一般情况下能实现快速收敛。     

Newton迭代法的P．C．格式

基于Newton迭代法 ,提出了一种有效的预估校正 (P .C .)迭代格式。本方法把牛顿法中的微分以差商替代 ,因而 ,既可求解具离散根的方程 ,也适用于有重根的方程求解 ,大量的试验结果表明 ,本文算法不仅可有效解决重根问题 ,而且有较高的收敛速度。     

电弧喷射推力器流动区域电磁场数值解法

为了有效地进行电弧喷射推力器流动区域全场(流场、电场、磁场和化学反应)耦合数值模拟并揭示其电磁场特征，建立了适当的电磁场模型，对电磁场控制方程的数值解法进行了详细研究。模型基于推力器稳态工作特征，控制方程由麦克斯韦方程组简化得到，采用有限控制容积积分方法离散椭圆型控制方程，采用9种不同的迭代方法求解离散方程。给出了不同迭代方法的收敛速度、数值稳定性和最终精度。研究表明，Gauss—Seidel逐线超松弛迭代法是推力器流动区域电磁场离散方程的一种快速有效的数值解法。     

一种不需要迭代的发动机辅助变量建模方法

现代航空发动机建模计算时常用牛顿迭代法求解非线性方程组,但由于发动机模型非线性太强且十分复杂而不易在全包线收敛。本文针对这一不足之处提出了一种不需要迭代就可以求解此非线性方程组的方法:引入航空发动机容积动力学中的变量及方程作为辅助变量及辅助方程,从而使整个非线性方程组闭合而可以直接求解。本文将此方法成功应用于某涡扇发动机,取得了良好的效果,验证了本方法的可行性。      

全变差正则化方法对一类热传导反问题的应用

热传导终端反问题是典型的不适定问题,应用基于全变差的正则化方法来获得它的稳定数值解,将近似解的求解范围由连续函数空间扩大到了有界变差函数空间,并将固定点迭代法用于求解非线性的Euler方程,并进行了理论分析和数值实验。结果表明:该方法与算法具有高效稳定全局收敛等优点。     

基于物面弱边界条件的并行全隐式流场求解方法

提出了一种基于物面弱边界条件的Navier-Stokes方程并行全隐式求解方法。首先,实现了非结构网格格点格式有限体积法的粘性物面弱边界条件。以此为基础,推导了相应的粘性Jacobian矩阵。相比传统的物面强边界条件,弱边界条件的应用简化了该推导过程。接着,为了解决Spalart-Allmaras(S-A)模型的收敛和保正性问题,提出了基于物面弱边界条件的S-A模型无条件保正性收敛方法,通过构造特别的隐式矩阵同时保证S-A模型的收敛和湍流工作变量的正值。流场方程全隐式时间推进和无条件保正性收敛方法所得到的线性方程组采用多色高斯-塞德尔迭代法进行求解,并由MPI实现并行化。通过若干二维和三维算例对本文所提出的方法进行了测试,其结果表明:1)随着罚系数的增大,基于弱边界条件的结果趋向于基于强边界条件的结果;2)对于力的计算,弱边界条件自动满足通量一致关系而无需特别处理。     

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

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

二维不可压分离流计算的约束迭代法

本文提出约束迭代法以改进不可压分离流计算中的迭代收敛性。在每个时间步的迭代中,对涡量传输方程的边界涡量,预先进行近似估计。因此它相当于一种隐式的边值算法。 两个翼型的算例表明,约束迭代法具有高的收敛速度。由于无需寻找最优松弛参数,因而节省人力,便于使用。     

超声速化学反应流动的LU-SGS伪时间迭代空间推进求解

介绍了求解超声速多组分有限速率化学反应流动的伪时间迭代(lower-upper symmetric-Gauss-Seidel,LU-SGS)方法.空间推进求解多组分抛物化(parabolized Navier-Stokes,PNS)方程时,在一个推进面上采用修正的LU-SGS方法迭代至收敛,把得到的结果作为初值赋给下一个推进面.沿推进面依次迭代直至求解完整个流场.采用伪时间迭代LU-SGS方法求解化学反应PNS方程,计算结果的准确性和时间迭代求解完全Navier-Stokes(N-S)方程相当,求解效率提高一个数量级.      

轴对称跨音速进气道外罩绕流的高效差分算法

 本文考虑轴对称跨音速进气道外流场的有限差分计算方法。采用守恒型位流方程、贴体坐标网格和精确边界条件,根据最佳收敛准则对轴对称情形设计出新的近似因式分解迭代格式,并将所设计的格式应用到计算皮托式进气口的轴对称跨音速流场。对几种典型进气口的计算表明,本文格式收敛快,计算结果与实验符合很好。     

ITERATION SOLUTION OF THE MIXED FORMULATION IN NONLINEAR FEM

ＩＴＥＲＡＴＩＯＮ　ＳＯＬＵＴＩＯＮ　ＯＦ　ＴＨＥ　ＭＩＸＥＤ　ＦＯＲＭＵＬＡＴＩＯＮ　ＩＮ　ＮＯＮＬＩＮＥＡＲ　ＦＥＭＷａｎｇＳｈｏｕｍｅｉ（ＢｅｉｊｉｎｇＵｎｉｖｅｒｓｉｔｙｏｆＡｅｒｏｎａｕｔｉｃｓａｎｄＡｓｔｒｏｎａｕｔｉｃｓ，Ｂｅｉｊｉｎｇ...     

流体网络节点残量修正算法的改进

为了提高流体网络求解的收敛性与稳定性,对网络节点残量修正算法进行改进。基于理论分析,提出使残量修正算法大范围收敛的初值赋值规律,并通过简单模型的数值分析加以验证;在节点残量修正显式迭代的基础上,提出隐式梯形迭代格式,扩大最大迭代步长。基于改进方法,以某典型航空发动机空气系统为例进行计算,获得空气系统各节点的压力、温度以及各支路的质量流量。该计算过程显示,相对于残量修正显式迭代,改进计算振荡大幅度减弱,收敛速度提高了80%以上。结果表明:改进算法有效地提高了流体网络求解的稳定性和收敛性。     

NEW TIME-EXPLICIT ASYMPTOTIC METHOD FOR THE SOLUTION OF AN EXACT CONTROLLABILITY PROBLEM OF SCATTERING WAVES

In this paper a new time-explicit asymptotic method is presented through introducing an auxiliary parameter for the solution of an exact controllability problem of scattering waves. Based on an exact control function, the asymptotic iteration is controlled by the auxiliary parameter and the optimal auxiliary parameter is updated during the iteration based on the existing or current iterated solutions of the wave equation. The numerical results show that the new method presented has a significant advantage over the purely asymptotic method in the history of convergence and has the ability to solve the scattering by the multi-bodies.     

解大型结构特征值问题的同时迭代的一个新算法

本文提出了同时迭代法的一种变型算法,用来求解大型结构动力分析中的广义特征值问题。使用了E_k子空间和特征方向的概念,证明了该算法的收敛率,并讨论了如何在计算机上实现的某些细节问题。通过统计每步迭代所需运算量和大量的实际计算,表明该算法与当前几种常用的同时迭代法或子空间迭代法相比,收敛率相同,但减少了每步迭代所需的计算量,提高了计算效率。     

高阶DG格式多重网格方法构造中的不同隐式策略影响研究

DG方法是一种非常具有潜力的高精度方法,但其在对复杂外形的数值模拟方面仍存在内存需求量大、计算量巨大等不足.为了进一步提高DG方法求解Euler方程的效率,在传统p型多重网格的基础上,结合LU-SGS和GMRES两种隐式迭代方法,研究其整体加速性能.p型多重网格方法通过对不同阶次多项式近似解进行递归迭代求解,来达到加速收敛的目的.高阶近似(p＞0)使用显式龙格库塔格式,最低阶近似(p=0)使用隐式格式.对NACA0012翼型和ONERA M6机翼跨音速无粘流动进行数值模拟,结果表明:与显式TVD-RKDG时间格式相比,DG(p0)层上采用LU-SGS和GMRES的p型多重网格方法收敛速度均得到明显提高,且GMRES迭代法性能最佳,LU-SGS迭代法次之.     

三维跨音速流有限差分计算中的某些问题

一、基本方法简介 1.三维跨音速小扰动方程 柱坐标下的三维跨音速小扰动方程为