适用于混合网格的并行GMRES+LU-SGS方法

给出了一种适用于混合网格的并行无矩阵GMRES+LU-SGS隐式时间格式。首先采用LU-SGS方法迭代若干步以获得一个合适的初场,然后切换到GMRES方法在每一时间步内近似求解,并将LU-SGS方法作为其预处理器。为加速收敛,将CFL数随着残差的降低逐步放大;为减少存储量和计算量,通量Jacobian采用无矩阵处理。在保证与串行执行一致的前提下,采用基于共享内存的OpenMP方法实现了并行计算,并通过对网格的分组避免了内存争夺。算例验证表明,方法极大地提高了计算收敛效率,并行结果与串行结果完全一致,计算结果与实验结果吻合较好。     

超声速进气道可压及不可压流动数值模拟

采用Chio-Merkle预处理矩阵对可压NS(Navier-Stokes)方程组进行时间预处理,分析了预处理方法的物理和数学背景.用有限体积方法,结合LU-SGS(Lower-Upper-Symmetric-Gauss-Seidel)隐式时间积分和AUSM+(P)(Precondition Advection Upstream Splitting Method)格式、中心差分2种空间离散格式,求解该预处理NS方程组.通过圆弧凸包流动和二维方腔顶盖驱动流动的数值试验表明,该方法克服了传统时间迭代方法模拟低速流动时的刚性问题,加速了收敛过程,可以同时有效地模拟可压及不可压流动.用此方法与块结构化网格技术相结合,进行了在不同飞行马赫数、攻角和出口反压条件下,混压式轴对称超声速进气道三维流场数值模拟.计算结果表明:该方法能准确地捕获复杂激波系,清楚地揭示进气道三维流动现象;进气道性能随工作条件的变化规律,与理论分析一致.      

ADI-FDTD在激励源和完纯导体边界处的精确格式

证明了当采用交替方向隐式时域有限差分法(AD I-FDTD)计算时,完纯导体和激励源处的算法格式与无源空间的算法格式不同,并推导了AD I-FDTD算法在源和完纯导体所在网格处的精确格式。利用常规的FDTD、无源空间的AD I-FDTD以及在源和导体处采用精确格式的AD I-FDTD,分别计算了点源激励的在自由空间传播的二维TE波、TM波垂直入射于金属条带FSS的散射场。计算结果表明,如果计算域中包含导体或激励源,直接采用通常的AD I-FDTD格式会带来很大的计算误差,而采用该精确格式后,计算结果与FDTD的计算结果吻合,从而证明了本文导出的AD I-FDTD在导体和激励源处的精确格式的正确性和使用的必要性。     

三维Euler方程组隐式LU分解的高性能并行计算

许多非定常无粘流体力学问题的数值模拟都需要利用Euler方程组来进行计算,而由于在隐格式下,所选取的时间步长可以比在显格式下时大得多,所以隐格式越来越受到重视,其中隐式LU分解是最常用的方法之一。对三维Euler方程组,采用隐式LU分解进行计算时,网格点所在的各个对角阵面之间存在数据依赖关系,本文分析了采用区域分解且边界上用显格式代替隐格式进行计算的高效性,在长方体建筑物内的爆炸模拟表明,在有112个CPU的某MPP巨型机上,并行计算效率超过60%。本文还分析了计算结果与串行计算时的差异,以及利用区域重叠减小这种差异的方法,同时考虑了对处理器进行合理的逻辑组织,将计算网格映射到处理器网格,以最大限度减少通信开销的方法。文中最后以一个爆炸毁伤的例子实际说明了所述方法的可行性与高效性。     

几种典型迎风格式的分析与比较

重点分析了AUSM 格式的构造方法,同时将其与Van Leer, Roe格式等构造方法进行比较和分析,通过求解一维激波管问题,以及二维和三维跨音速流场的计算比较了AUSM 格式,VanLeer格式,Roe格式和中心格式的数值模拟结果和网格适应性,以及在经典四步Runge-Kutta显式时间推进和LU-SGS隐式时间推进下上述四种格式的计算效率.获取了对于基础研究和工程应用很有意义的结论.     

隐马尔可夫模型在对象定位中的应用与实现

以医学图像中脊柱的定位为应用背景,对经典Viterbi算法进行改进,根据全局最优搜索算法提出了可控的最优搜索算法,把隐马尔可夫模型(Hidden Markov model,HMM)应用于二维图像中线性结构对象的定位。以真实的人体脊柱X光照片为实验样本进行的定位实验,结果表明:HMM模型不需要对图像施加较多的限制,就可获得较成功的定位结果,同时引入了可控的最优搜索算法后,定位效果有了进一步的提高。还通过实验分析了可控制的最优搜索的两个控制参数对定位性能的影响,并提出了定位性能具有区间性的结论。     

一种改进的基于LDPC码的信源信道联合译码方法

提出一种改进的将LDPC迭代译码和隐马尔可夫估计相融合的迭代译码算法。一方面通过隐马尔可夫信源估计获取信源残留冗余,作为外信息提供给LDPC迭代译码;另一方面,通过LDPC迭代译码获得更精确的信源估计模型参数,加快迭代收敛速度,获得更好的迭代译码性能。其优点是无需知道任何关于信源的先验信息,计算复杂度低,并行操作性强,尤其适用于宽带多媒体传输。     

固体火箭发动机前封头人工脱粘层烧蚀研究

从理论分析的角度出发．通过对人工脱粘层和绝热层所构成的槽缝燃气温度场与烧蚀的耦合计算，得出了人工脱粘层的烧蚀规律。在计算中，采用了隐式差分格式及控制容积法等计算技术，并采用了先进的烧蚀模型。为模拟发动机静止与飞行两种不同状态，本文对不同宽度的槽缝人工脱粘层的烧蚀进行了计算，对其进行了比较。本文所采用的方法可供今后人工脱粘层的设计参考。     

基于双重隐马尔科夫模型的ISAR图像序列中散射中心关联算法

散射中心关联是雷达目标三维重构的关键技术,由于目标散射中心在连续观测视角下在序列图像间的运动轨迹及散射强度具有强相关特性,因此根据相关性建立的双重隐马尔科夫模型能够鲁棒地表征序列图像间散射中心关联状态的变化。大量的实验结果表明,该算法能够实现散射中心的鲁棒关联。