Stokes问题的集中质量非协调有限元法

针对其平面有界凸区域上的一类非定常不可压Stokes方程，提出了一种新的有限元方法，借助于所谓的速度一压力混合型公式，讨论了非定常不可压Stokes方程的质量集中非协调有限元逼近格式(全离散情形)。首先，给出了非定常不可压Stokes方程的质量集中非协调Galerkin有限元逼近的全离散格式，其次，对所讨论问题的解与其所给出的离散问题的解之间的误差进行了分析研究；最后，利用Stokes投影算子的性质和离散的LBB条件，得到了非定常不可压Stokes方程关于速度L2模和能量模及压力L2模方面的最优阶误差估计。     

一种求解背包问题的混合遗传算法

将启发式贪婪算法与简单遗传法相结合构成一种混合遗传算法,用该混合遗传算法求解大规模背包问题,含有５０个物品的背包问题的试算结果表明,用这种方法求解大规模背包问题,其解的质量和求解性能较简单遗传算法和贪婪算法都有很大的改善。     

求解广义特征值问题的并行保域多分法

近年来，随着并行机的发展，提出了代数特征值问题的并行多分法，但国内外的研究工作迄今仅限于对称三对角矩阵的标准特征值问题。在科学与工程众多领域内有着重要应用的广义特征值问题的多分法，因难度大等方面原因尚无人研究。本文提出广义特征值问题的并行保域多分法，该算法适用于大型稀疏实对称矩阵广义特征值问题的求解，它克服了传统的广义特征值问题的对分法（行列式查找法）出现的漏根或迭代不收敛等缺点，并保持其优点。作者在ＹＨ－１向量机上对这一算法进行了数值实验，并与并行保域行列式查找法作了比较。数值结果表明，该算法具有较高的加速比，当系统自由度为２１１４、求解特征对个数为３时，加速比可达７．７；且当问题规模较大时，并行保域多分法优于并行保域行列式查找法。     

求解广义特征值反问题的数值方法

讨论一类广义特征值反问题的数值解法,这类问题包括加法、乘法和经典特征值反问题作为其特殊情况。基于行列式和最小奇异值的计算,文中给出了求解这类问题的两个二次收敛的数值方法,描述了在出现重特征值的情况下如何改进其中的一个方法以保持二次收敛性,并且给出了两个数值例子以解释收敛性结果     

多层介质涂覆腔体电磁散射特性的IBC/FEM分析

将Leontovich阻抗边界条件(IBC)与矢量有限元法(EB-FEM)相结合对内壁涂覆多层介质开口腔体的电磁散射特性进行分析。将IBC和腔体口面边界积分方程分别转化为第三类边界条件形式，从而导出边值问题的相应泛函。表面输入阻抗通过广义反射系数得到，计入了介质内部电磁波的多次反射，提高了IBC／FEM混合法的精度。数值结果表明这种IBC／FEM方法在分析涂覆多层介质开口腔体电磁散射问题中的便捷性与精确性。     

基于修正偶应力理论的正弦剪切变形圆板的弯曲屈曲研究

利用修正偶应力理论(Modified couple stress theory,MCST)和正弦剪切变形理论,分析了圆形微孔板的轴对称弯曲和屈曲行为.通过最小总势能原理推导得到微分形式的控制方程和边界条件,并用引入的标称变量表示.应用广义微分求积法,将微分形式的控制方程和边界条件离散化,得到一组线性方程.通过求解线性方...     

亚声速大迎角模型试验洞壁干扰修正方法研究

通过测量洞壁附近的压力分布来模拟透气壁试验段的洞壁边界条件,采用数值求解Euler方程的方法模拟模型在风洞中的绕流场,然后将洞壁边界条件用远场边界条件替换,用同样方法计算出模型在自由流中的绕流场,从而计算出洞壁干扰对模型的绕流场和气动力的影响。针对大迎角模型的洞壁干扰问题的特殊性,采用重叠网格技术,对飞机大迎角标模和SB-03模型进行了亚声速大迎角的洞壁干扰计算与修正。     

基于MPP编程环境的一种并行子空间迭代法

子空间迭代法是科学与工程计算中求解广义特征值问题的有效方法 ,针对向量机和共享内存的多处理机 ,前人已成功地作了并行处理。文中给出了适合 MPP大规模并行计算机的并行子空间迭代法。该算法将广义特征值问题转换为一般特征值问题 ,其计算工作量主要体现在矩阵乘法 ,通过对该方法作并行处理 ,使矩阵求逆及一部分乘法运算转换为各结点机上三角形方程组的并行求解。在大规模并行计算机 PA R95上结合 J8- II机翼的动力特性问题对该算法作了数值试验 ,结果说明所给算法是非常有效的     

非对称广义特征值问题的拟-Eberlein算法及其并行化

非对称广义特征值问题的并行计算,目前在国内外研究得很少, Ｇ． Ｗ ． Ｓｔｅｗ ａｒｔ 和 Ｐ． Ｊ． Ｅｂｅｒｌｅｉｎ 曾分别研究非 Ｈｅｒｍ ｉｔｅ 矩阵标准特征值的并行拟 Ｊａｃｏｂｉ算法,１９８９ 年 Ｊ． Ｐ． Ｃｈａｒｌｉｅｒ 和 Ｐ． Ｖａｎ Ｄｏｏｒｅｎ 在 Ｇ． Ｗ ． Ｓｔｅｗ ａｒｔ 的工作基础上提出了求解非对称广义特征值问题的拟 Ｊａｃｏｂｉ算法（简称 Ｃ Ｖ 算法）与并行拟 Ｊａｃｏｂｉ算法。文中以 Ｊ． Ｐ． Ｃｈａｒｌｉｅｒ 等人的工作为基础,提出求解大型非对称广义特征值问题的拟 Ｅｂｅｒｌｅｉｎ 算法与并行拟 Ｅｂｅｒｌｅｉｎ 算法, Ｃｈａｌｌｅｎｇｅ Ｌ 并行系统上的数值试验表明,不仅并行效率很高,且敛速远优于 Ｃ Ｖ 算法     

集值广义强向量拟均衡问题系统解的存在性

研究集值广义强向量拟均衡问题系统,利用KFG不动点定理得到了具有控制结构的集值广义强向量拟均衡问题系统解的存在性.当集值广义强向量拟均衡问题系统解退化为集值广义强向量拟均衡问题解时,即为特例.     

A FINITE ELEMENT ADAPTIVE METHOD FOR SOLVING GENERALIZED STOKES PROSLEM

ＡＦＩＮＩＴＥＥＬＥＭＥＮＴＡＤＡＰＴＩＶＥＭＥＴＨＯＤＦＯＲＳＯＬＶＩＮＧＧＥＮＥＲＡＬＩＺＥＤＳＴＯＫＥＳＰＲＯＳＬＥＭＤ．Ａ．ＦｏｋｉｎＹａｎｇＺｕｏｓｈｅｎｇ（ＤｅｐａｒｔｍｅｎｔｏｆＡｅｒｏｄｙｎａｍｉｃｓ，ＮＵＡＡ２９ＹｕｄａｏＳｔｒｅｅ...     

二维多介质流动问题的界面处理方法

就二维可压缩多介质流动问题的数值模拟,给出了一种新的界面处理方法。通过在界面处构造Riemann问题,利用Riemann问题的解分别定义界面两边流体的边界条件,由于Riemann问题的解准确地描述了界面处流体的流动状态,因此得到了更加准确的界面边界条件。本文将由Riemann问题的解得到的界面速度外推到整个流场,重新定义速度场,避免了由于速度的大梯度变化而导致的Level—Set等值线相互交错,因而得到了更加精确的界面位置。利用该方法对水下激波与柱型气泡相互作用问题进行数值模拟,结果表明该方法能准确地捕捉各种物理现象。     

双重互易边界元法在求解N—S方程中的应用

提出了用双互易边界元法求解Ｎ－Ｓ方程的新思路。用Ｐｅａｃｅｍａｎ－Ｒａｃｈｆｏｒｄ算子分理解将时间相依的Ｎ－Ｓ方程分解为线性和非线性子问题，线性问题用共轭梯度法解除压强－速度的耦合；非线性问题进行局部线化。对所得的Ｐａｓｓｉｏｎ方程及相关类型方程采用双重互易的边界元法求解，消除了传统Ｎ－Ｓ方程边界元解法的区域积分问题。     

外壁热流响应计算的导热反问题方法及其验证

研究组合结构外壁热响应特性,根据Beck’S序列函数法的思想,优化算法各个环节,建立应用多个温度测点、多个未来时间步的计算模型。通过假设热流呈三角波形变化开展数值验证,计算结构内部温度场,并取其中数点温度信息作为反问题输入条件反演得到外壁热流响应,得到热流误差不超过O．5％、温度误差不超过0．01％的理想结果,并得到测点数量及位置等因素对计算误差的影响。该计算方法可应用于包装结构在模拟火灾等特殊环境中的外壁热响应研究。     

基于PIV技术的压力场重构算法实现与研究

介绍了有限容积法、直接积分法和Poisson方程法3种基于PIV瞬时速度场重构压力场的基本原理以及相应的计算方法，选取管流突扩流场和偏置方块绕流流场两个不可压缩流场的瞬时速度场数据，采用上述3种压力场重构算法，分别研究了图像噪声、速度场精度、插值算法以及边界条件的类型与精度对重构压力场的影响。最后针对管流突扩过程第20ms的流场，给出了3种重构算法下的压力场云图以及对应的CFD模拟结果。研究表明，有限容积法和直接积分法容易受到噪声的影响而产生剧烈震荡，但是可以在较大的速度场误差范围内保持较高的精度，通过采用双线性插值可以获得更高精度的重构压力场；Poisson方程法不易受到噪声的影响而产生震荡，同时在高精度PIV速度场下的优势较为突出，通过采用双三次差值可以获得更高精度的重构压力场；混合边界条件仅仅测定边界上有限个点的压力值，就获得了接近狄利克雷边界条件下重构压力场的精度，远高于诺依曼边界条件；边界条件的误差严重降低重构压力场的精度，其影响程度比速度场误差还要大。     

Banach空间中椭圆方程Dirichlet问题的迭代解

以半序理论为工具，在Ｂａｎａｃｈ空间中用上下解方法研究了椭圆方程Ｄｉｒｉｃｈｌｅｔ问题的迭代解，并证明了存在两个迭代序列分别一致收敛于该问题的最小解和最大解。     

ON SPECTRAL PROBLEM FOR VISCOUS SHEAR FLOWS

INTRODUCTIONLarge- scale complex eigenvalue problemsarising in continuum physics and engineering ap-plications are important not only because theyyield fundamental knowledge on the systems be-ing studied,they also provide a challenging con-textforthe developmentof numerical and analyt-ical techniques in their solution.In this paper,we re- examine the well- known Orr- Sommerfeldproblem which governs the linear stability of par-allel shear flows.More specifically,we are con-cerned with the …