用Chebyshev多项式加速的子空间迭代法

研究计算大型稀疏对称矩阵的若干个最大或最小特征值的问题，首先引入了求解大型对称特征值问题的子空间迭代法和Chebyshev迭代法，并对后者作了理论分析。为了加速子空间迭代法的收敛速度，作者用Chebyshev多项式来改进原始的子空间迭代法，即讨论Chebyshev迭代法对子空间迭代法的应用，从而给出了Chebyshev-子空间迭代法。最后把原始的方法和改进的方法计算数值例子的结果进行了比较，其结果表明Chebyshev-子空间迭代法比子空间迭代法优越，不仅收敛速度快，并且减少了计算量和计算时间。     

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

子空间迭代法是科学与工程计算中求解广义特征值问题的有效方法，针对向量机和共享内存的多处理机，前人已成功地作了并行处理。文中给出了适合ＭＰＰ大规模并行计算机的并行子空间迭代法。该算法将广义特征值问题转换为一般特征值问题，其计算工作量主要体现在矩阵乘法，通过对该方法作并行处理，使矩阵求逆及一部分乘法运算转换为各结点机上三角形方程组的并行求解。在大规模并行计算机ＰＡＲ９５上结合Ｊ８－ＩＩ机翼的动力特性问     

CFD问题显式有限差分方法的并行化

基于ＭＩＭＤ并行计算机模型，本文讨论了ＣＦＤ问题应用显式有限差分方法的并行化问题。利用区域分解法将计算问题分解为多个子问题，每个子问题由不同的处理器分别处理。针对ＳＣＢ计算格式的特点，尽量降低各处理器间的数据通讯，提高了并行计算效率。最后给出二维Ｅｕｌｅｒ方程组计算实例，计算结果令人满意，且算法有较好的可扩放性。     

瞬态热传导问题的隐式EBE并行算法

结合有限元法，研究了ＥＢＥ策略在热传导问题并行计算方面的应用，给出了瞬态热传导问题的隐式ＥＢＥ并行算法。根据所用ＣＰＵ个数并行计算单元刚度矩阵，而总刚度矩阵Ｋ不需要组装，仅在“单元级”上进行各种计算。就模型实验规模１０２４阶而言，机器存存储量比传统算法降低了４０％左右，在拥有共享内存的并行计算机Ｃｈａｌｌｅｎｇｅ上的数值实验表明，问题的规模越大，并行效率越高，所给算法非常有效。     

自动电位滴定法测定啤酒中氯离子含量

本文讨论了自动电位滴定法测定啤酒中的离子含量，且将原用二阶微商法求终点电位，用计算机ＢＡＳＥＩＣ语文程序代替，简化了复杂的计算。用此法测得市购南昌啤酒中氯离子含量为１４１．８－１４２．５ｍｇ．Ｌ＾－１。回收率的范围为：９８．６－１０２．４％。方法精密度：连续测定同一样品１０次，标准偏差０．２９，变异系数０．２％。     

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

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

非对称线性方程组的块广义最小向后误差算法

许多科学领域都需要求多个右边值的大型非对称线性方程组，使用块方法同时计算所有的方程组比分别计算每一个方程要有效得多。因此，能同时计算所有方程的块迭代方法比单独计算每一个方程的迭代法要有效得多。本提出了一个块GMBACK方法求解有多个右边值的大型非对称线性方程组，该方法利用块Arnoldi过程构造Krylov子空间来求解Xm∈X0 Km(A,R0)使得矩阵A的扰动范数最小。     

求解大型对称线性方程组的循环收缩Lanczos算法

向Krylov子空间中加入一些模接近于零的特征值对应的特征向量能够加快收敛速度，事实上，对于这些模接近于零的特征值对应的特征向量，可以用Krylov子空间方法得到，并且在新的Krylov子空间形成的过程中，近似特征向量的近似度会不断提高，特别在标准Krylov子空间方法中，如果因为这些特征向量而减缓了收敛速度，则随着这些特征向量的近似度的提高，用增广Krylov子空间方法解线性方程组的收敛速度会明显加快。Lanczos算法是求解大型对称不定线性方程组的有效方法之一。但在计算过程中由于Lanczos向量失去正交性减慢了收敛速度。本文根据增广Krylov子空间方法提出循环收缩Lanczos算法，新算法充分利用Lanczos过程所得到的谱信息，确定预处理，从而加速Lanczos算法的收敛速度。     

进口速度分布对短突扩压器性能的影响

在任意曲线坐标系下，对环型燃烧室的三维流场进行了数值计算，研究了四种不同的进口速度分布对扩压器性能的影响。计算结果表明，进口速度大小和方向都对扩压器性能有重要影响，进口速度分布均匀以及方向贴近扩压器壁面，有利于提高扩压器的性能。计算中采用标准ｋ－ε双方程紊流模型，采用控制容积法进行离散，在非交错网格体系下用ＳＩＭＰＬＥ法求解     

管式电除尘器中粉法颗粒运动轨迹的数值模拟

对管式电除尘器中的气固两相流动流场进行数值模拟，实质上是研究粉尘颗粒在有电场力作用下的运动轨迹。计算中，将气相作为连续介质，采用Ｋ－ε双方程湍流模型，并用ＳＩＭＰＬＥ算法对流场进行数值模拟；将固相作为离散体系，采用颗粒轨迹法计算其运动轨迹。计算时分别选取４种电压、４种气流速度和４种颗粒直径为计算工况，计算结果显示出颗粒运动特性和除尘效率，并指出其与电压、气流速度和颗粒粒径等３个主要因素相关。计算结     

雷达散射截面值的多机并行计算

为了提高计算效率,本文讨论了飞行器的RCS值的并行计算方法。基于自带存储器的多处理机系统,对多任务模式和农场主-雇员模式都作了探讨。发现这两种模式都适用于RCS值的并行计算,因此开发了计算RCS值的并行处理程序PARCS。实例表明其并行效率是很高的。     

Weak GPS L1 Signal Acquisition Based on BPDC

It is difficult to achieve accurate acquisition of weak global positioning system（GPS） signals with traditional methods. A weak signal acquisition strategy based on block processing and differentially coherent （BPDC） is put forward after analyzing the advantages and disadvantages of coherent and non-coherent integration algorithms. Code phase parallel search of the pre-coherent integration is conducted by using fast Fourier transform（FFT）, and the results are then differential coherent processed and block processed. BPDC method reduces computation cost compared with coherent and non-coherent（CNC） algorithm. The performance of the two algorithms is also compared based on simulated signals. The result shows that the noise suppression effect of BPDC algorithms is superior to that of traditional CNC algorithm, and the superiority of BPDC is more apparent with the reduction of carrier to noise ratio （CNR）. In the case that the pre-coherent integration length is 4 ms and CNR is reduced to 28 dB-Hz, CNC algorithm cannot yet acquire signal correctly while BPDC has well acquisition performance. Therefore, for weak GPS signal acquisition, BPDC algorithm can acquire the signal with lower CNR and has better acquisition property.     

有限元结构分析并行计算的若干研究进展

有限元结构分析并行处理在大规模科学与工程计算中占有重要地位，本文简要介绍作者在这一领域内的研究工作及其成果。涉及单元刚度矩阵的并行计算、有限元方程组的并行直接解法、有限元方程组的并行迭代解法、结构动力分析并行直接积分法、广义特征值问题的并行算法以及EBE（Element-By-Element）技术在有限元结构分析并行处理中的应用等研究领域。特别提出了“伪单元向量”与“单元分组技术”等概念与新方法；首先考虑了矩阵向量积的多自由度问题；系统、深入而又全面地研究了EBE策略在结构分析并行处理中的应用，填补了国内空白并突破了国际现有成果。     

基于学习的并行免疫量子进化算法

提出了基于学习的多宇宙并行免疫量子进化算法,算法中将种群分成若干个独立的子群体,称为宇宙。并给出了多宇宙的并行拓扑结构,其中各宇宙独立演化．宇宙内采用免疫量子进化算法,宇宙间采用基于学习机制的移民、模拟量子纠缠的种群交叉等信息交互方式．使得进化算法具有更好的种群多样性,更快的收敛速度和全局寻优能力。仿真实验结果表明该算法比串行的免疫量子进化算法运算效率更高。     

湍流燃烧模拟中化学反应的加速算法研究进展

为研究湍流燃烧数值模拟中化学反应机理计算的加速方法，讨论了动态自适应化学（Dynamic Adaptive Chemistry，DAC）方法和Krylov子空间近似的指数格式的应用情况。在湍流火焰大涡模拟中，使用DAC简化可以加速化学反应计算。然而，在并行燃烧数值模拟中，处理器核心的负载极度不平衡，加速效果有限。而Krylov子空间近似的指数格式的加速效果可以作用于每个处理器核心，更有利于整体计算效率的提高。在同等精度下，相比于隐式格式耦合DAC和MTS加速方法，Krylov子空间近似的指数积分格式对化学反应计算的加速效果更为显著。     

超长大直径钻孔灌注桩桩端承载力研究

阐述了超长大直径钻孔灌注桩桩端承载力的影响因素,包括施工工艺、深度效应和尺寸效应等。分析和对比了目前国内常用的各规范所采用的计算方法。文中算例的计算结果表明,对于同等条件的桩,各规范的计算结果相差较大。采用自平衡试桩法获得了几个重大工程的超长大直径钻孔灌注桩的桩端承载力的实测值。将工程实例的测试结果与各规范的计算结果进行对比,证实了目前规范的计算方法不尽合理。目前各规范中桩端承载力的计算方法是在针对小直径的中短桩的研究基础上建立的。如今,大型桥梁工程中普遍采用超长大直径钻孔灌注桩的情况下,对其桩端承载力的计算应当进行相应的调整,使设计更为合理。     

Numerical Simulation of ATPS Parachute Transient Dynamics Using Fluid-Structure Interaction Method

In order to simulate and analyze the dynamic characteristics of the parachute from advanced tactical parachute system(ATPS),a nonlinear finite element algorithm and a preconditioning finite volume method are employed and developed to construct three dimensional parachute fluid-structure interaction(FSI)model.Parachute fabric material is represented by membrane-cable elements,and geometrical nonlinear algorithm is employed with wrinkling technique embedded to simulate the large deformations of parachute structure by applying the NewtonRaphson iteration method.On the other hand,the time-dependent flow surrounding parachute canopy is simulated using preconditioned lower-upper symmetric Gauss-Seidel(LU-SGS)method.The pseudo solid dynamic mesh algorithm is employed to update the flow-field mesh based on the complex and arbitrary motion of parachute canopy.Due to the large amount of computation during the FSI simulation,massage passing interface(MPI)parallel computation technique is used for all those three modules to improve the performance of the FSI code.The FSI method is tested to simulate one kind of ATPS parachutes to predict the parachute configuration and anticipate the parachute descent speeds.The comparison of results between the proposed method and those in literatures demonstrates the method to be a useful tool for parachute designers.     

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

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