一种新的红外焦平面阵列非均匀性代数校正算法

 针对传统红外图像非均匀性代数校正算法收敛速度慢、运动估计精度不高的缺点,提出一种多尺度光流帧间运动估计的非均匀代数校正算法。通过时域低通滤波,采用多尺度光流估计下一帧图像,将所得图像对进行代数校正。该方法性能在自主研发的热像仪中得到验证,和传统的红外图像非均匀性代数校正算法相比,算法收敛速度从60帧提高到25帧,提高一倍左右,运动估计精度从校正后图像方差76.539减少到32.482。     

用代数多重网格法求解一维分裂格式的Euler方程

提出了代数多重网格法（ＡＭＧ）的一种新算法。新算法改进了插值公式和粗网格方程，并把它应用到求解一维的分裂格式Ｅｕｌｅｒ方程。数值结果表明，对于具有高ＣＦＬ条件数的Ｅｕｌｅｒ方程，代数多重网格法可以求解；对于Ｇａｕｓｓ－Ｓｅｉｄｅｌ方法求解不能收敛的代数方程组，代数多重网格法求解可以收敛。新算法改进了代数多重网格法的收敛性和扩展和它的应用范围，数值结果表明了它的有效性和强壮性。     

基于LU-SGS的非结构弹簧网格迭代算法

针对运动间断拟合中需频繁更新网格点位置的特点,提出一种基于LU-SGS(lower-upper symmetricGauss-Seidel)迭代方法的非结构弹簧网格运动算法.根据弹簧网格原理构建与网格拓扑关系相对应的稀疏系数矩阵,将LU-SGS思想成功引入动网格迭代算法,并辅以合理的网格运动管理策略,实现动网格的快速迭代.研究表明,在非结构网格下,LU-SGS算法可以满足运动间断拟合的需求,在流场隐式时间推进时,仍能保证获得稳定解;与传统的SOR方法相比,计算时耗减少20%以上.     

一种可扩展的三维半导体器件并行数值模拟算法

在基于漂移-扩散模型的三维半导体器件数值模拟中,通过有限体积法进行数值离散,采用完全耦合的牛顿迭代求解非线性代数方程组,并使用基于代数多重网格预条件子的GMRES方法求解牛顿迭代中的线性方程组,构造一种稳健且高度可扩展的非结构四面体网格上求解半导体方程的并行算法.基于PHG平台实现该算法的并行计算程序,并对PN结和MOS场效应晶体管等问题进行了最大网格规模达到5亿单元、最大并行规模达到1 024进程的大规模数值模拟实验,结果表明,该算法计算效率高,可扩展性好.     

求解二维三温辐射扩散方程组的一种代数两层迭代方法

在二维三温辐射扩散方程离散代数方程组的求解中,由于光子、电子和离子温度之间存在耦合关系,而且三个温度在同种介质中有不同的扩散性质,使得经典的代数多重网格(AMC)方法难以直接应用.基于特殊粗化策略,在粗网格层解除了这种耦合关系,得到一种代数两层网格方法,而粗网格方程由经典AMG方法求解.将这一算法具体应用于JFNK(Jacobian自由的Newton-Krylov)框架中预处理方程的求解,并基于该框架求解二维三温辐射扩散方程组.数值结果显示了算法的可扩展性和健壮性.     

基于局部信息优先的代数多重网格粗化策略

网格粗化是实施代数多重网格方法的重要步骤之一。在网格粗化过程中,当出现多个最大集合势相同的网格点时,经典网格粗化策略中的集合势并不能精确表征网格点的整体强弱连接强度,导致网格粗化路径不唯一并影响后续求解阶段的计算效率。为克服这一不足,本文提出了一种基于局部信息优先原则的网格粗化策略,并发展了基于局部信息的整体强弱连接系数判别法和基于局部信息的二次粗化阈值判别法两种实施方法。数值计算表明:与经典网格粗化策略相比,本文提出的网格粗化策略可精确表征网格点之间的整体强弱连接强度,有效优化网格粗化路径,并在一定程度上提高了求解阶段的计算效率。     

光谱法分析代数迭代重建算法及改进

采用光谱分析技术,系列研究了已有代数迭代重建算法,并进一步研制了改进算法一简单自相关代数迭代重建算法(SSART).采用数值模拟的方法模拟重建一个含遮挡物场,研究了SSART模拟效果和重建精度.作为对照,分析了已有的传统代数迭代算法同条件下的重建结果,包括基本代数迭代重建算法(ART)、联合代数迭代重建法(SART)和改进的联合代数迭代重建法(MSART).重建精度用均方误差(MSE),绝对平均误差(AVE)和峰值相对误差(PE)三种误差指数标定.结果发现,SSART重建结果的三种误差指数都明显减小,其中MSE在10-4数量级上比ART的降低了26.6%,PE在10-2数量级降低了77.6%.因此,SSART重建含遮挡物模型场时,重建精度大幅度提高,重建效果明显改善,对传统代数迭代重建算法显著改进.     

非结构子网格压力方法

纯拉氏流体力学计算失败的根源是砂漏运动和虚假涡旋,运用子网格拉氏质量与其相关的子网格密度得到子网格压力,子网格压力产生的扰动隅角力能消除这些虚假运动,以便消除困扰拉氏流体力学算法的网格畸变带来的数值计算困难.将包括子网格压力的相容流体算法推广到非结构网格,研制包括非结构子网格压力的相容流体程序,对Saltzman活塞问题进行数值模拟,从数值结果分析可知,非结构子网格压力方法能抑制虚假涡旋,消除砂漏畸变.     

人造地球卫星运动方程的代数动力学算法数值解

代数动力学算法首次被用于求人造地球卫星运动方程的数值解，在四阶算法下，与Runge-Kutta算法和辛算法的计算结果作了比较.结果表明，代数动力学算法对于人造地球卫星长期轨道的计算有较高的精度.并讨论了地球四极和八极带谐项对卫星轨道的影响. 关键词： 人造地球卫星运动方程 代数动力学算法数值解 地球四极和八极带谐项对卫星轨道的影响     

基于最小刚性图代数特性的无线网络拓扑优化算法

对于能量受限的无线传感器网络,拓扑优化能够降低能耗,优化通信链路结构.本文基于最小刚性图原理提出了一种新的拓扑优化算法,算法综合考虑了生成拓扑链路图中通信链路的权值与生成刚性图的代数特性问题,既保证了通信链路较短,有利于延长网络的生命周期,同时使生成的通信链路图结构更加稳定,网络具有较好的鲁棒性.仿真实验表明,与相关算法比较,提出的算法中通信链路较短,具有较好的网络连通性与结构稳定性,同时生成刚性图矩阵的迹较大,具有较好的刚度代数性能.     

Efficient algebraic multigrid for migration–diffusion–convection–reaction systems arising in electrochemical simulations

The article discusses components and performance of an algebraic multigrid (AMG) preconditioner for the fully coupled multi-ion transport and reaction model (MITReM) with nonlinear boundary conditions, important for electrochemical modeling. The governing partial differential equations (PDEs) are discretized in space by a combined finite element and residual distribution method. Solution of the discrete system is obtained by means of a Newton-based nonlinear solver, and an AMG-preconditioned BICGSTAB Krylov linear solver. The presented AMG preconditioner is based on so-called point-based classical AMG. The linear solver is compared to a standard direct and several one-level iterative solvers for a range of geometries and chemical systems with scientific and industrial relevance. The results indicate that point-based AMG methods, carefully designed, are an attractive alternative to more commonly employed numerical methods for the simulation of complex electrochemical processes.     

Algebraic Method of Studying Excited Vibrational Transitions

The algebraic method of perturbation theory is used to solve problems in molecular spectroscopy. Formulas for the shifts in the band centers, changes of the form factors of excited vibrational transitions, and changes of the components of the kinematic interaction matrix are derived by this method for the nth order of perturbation theory.     

由扩散张量导出的各向异性扩散模型的隐式数值模拟

研究了由扩散张量导出的各向异性扩散的图像处理模型,并构造了隐式差分格式,形成了有13条对角线的大型稀疏矩阵.利用代数多重网格法求解了这个线性代数方程组.并进行了数值试验.     

An algebraic multigrid solver for transonic flow problems

This article presents the latest developments of an algebraic multigrid (AMG) based on full potential equation (FPE) solver for transonic flow problems with emphasis on advanced applications. The mathematical difficulties of the problem are associated with the fact that the governing equation changes its type from elliptic (subsonic flow) to hyperbolic (supersonic flow). The flow solver is capable of dealing with flows from subsonic to transonic and supersonic conditions and is based on structured body-fitted grids approach for treating complex geometries. The computational method was demonstrated on a variety of problems to be capable of predicting the shock formation and achieving residual reduction of roughly an order of magnitude per cycle both for elliptic and hyperbolic problems, through the entire range of flow regimes, independent of the problem size (resolution).     

Algebraic isometric embeddings of charged spherically symmetric space-times

Recent results concerning isometric embeddings of charged spherically symmetric space-times as algebraic 4-surfaces of a pseudo-Euclidean ⁿ are given. Several problems arising in algebraic isometric embeddings of space-times with horizons are analyzed by studying the Reissner-Nordström solution. The connection between the images of radial null geodesics and the possibility of extending the algebraic isometric embedding through a horizon is studied in detail.     

基于代数滑移模型的光催化反应器内催化剂颗粒分布特性研究

针对光催化制氢反应器内液固两相流的特点,运用代数滑移模型,对反应器内的催化剂颗粒-水两相流动进行了数值模拟。采用分块的结构化网格,并运用多重网格方法进行求解。通过与前人实验中压力梯度结果的比较,对模型进行了验证,并最终得到了反应器内催化剂颗粒的典型分布及相应的类型图。     

Algebraic Approach for Shape Invariant Potentials in?Klein-Gordon Equation

The Shape invariant method has the algebraic structure and its algebras are infinite-dimensional. These algebras are converted into finite-dimensional under conditions. Based on the property of this method we obtain the algebraic structure of some physical potentials in the s-wave Klein-Gordon equation with scalar and vector potential which satisfy the Shape invariant condition.     

Algebraic treatment of quantum-mechanical models with modified Coulomb potentials

Nonrelativistic models with modified Coulomb potentials are solved by an algebraic method based on SO(4,2) dynamical group. The nonrelativistic model with Yukawa long-range potential and the Stark effect in the nonrelativistic hydrogen atom in the homogeneous and inhomogeneous external electrostatic fields are studied in details. The algebraic method for the Dirac and Klein-Gordon relativistic hydrogen atom as well as relativistic models with the rotationally symmetric modified Coulomb potentials are also discussed.     

Algebraic computing in general relativity

The purpose of this paper is to bring to the attention of potential users the existence of algebraic computing systems, and to illustrate their use by reviewing a number of problems for which such a system has been successfully used in General Relativity. In addition, some remarks are included which may be of help in the future design of these systems.     

Algebraic Rossby Solitary Waves Excited by Non-Stationary External Source

The paper deals with the effects of non-stationary external source forcing and dissipation on algebraic Rossby solitary waves.From quasi-geostrophic potential vorticity equation,basing on the multiple-scale method,an inhomogeneous Korteweg-de Vries-Benjamin-Ono-Burgers(KdV-B-O-Burgers) equation is obtained.This equation has not been previously derived for Rossby waves.By analysis and calculation,four conservation laws associated with the above equation are first obtained.With the help of pseudo-spectral method,the waterfall plots are obtained and the evolutional characters of algebraic Rossby solitary waves are studied.The results show that non-stationary external source and dissipation have great effect on the generation and evolution of algebraic solitary Rossby waves.