基于概率加权共轭梯度算法的混凝土超声波层析成像

在混凝土层析成像中,为了提高反演的准确性和计算效率,针对共轭梯度算法提出一种加权算法——概率加权共轭梯度算法.新算法不同于常规的加权算法,权重是加在成像单元上而不是方程上.为取得较好的权重因子和较好的迭代初始值,采用IART算法的权重和迭代初始值的选取方法.模拟算例和混凝土试验均表明这种加权算法的可行性.     

Cell Conservative Flux Recovery and a Posteriori Error Estimate of Vertex-Centered Finite Volume Methods

A cell conservative flux recovery technique is developed here for vertex-centered finite volume methods of second order elliptic equations. It is based on solving a local Neumann problem on each control volume using mixed finite element methods. The recovered flux is used to construct a constant free a posteriori error estimator which is proven to be reliable and efficient. Some numerical tests are presented to confirm the theoretical results. Our method works for general order finite volume methods and the recovery-based and residual-based a posteriori error estimators are the first result on a posteriori error estimators for high order finite volume methods.     

Application of Robust Incomplete Cholesky Factorizations Preconditioned CG Algorithms for FEM Analysis of Millimeter Wave Filters

The Incomplete Cholesky factorizations preconditioning scheme is applied to the conjugate gradient (CG) method for solving a large system of linear equations resulting from finite element method (FEM) analysis of millimeter wave filters. As is well known, the convergence of CG method deteriorates with increasing EM wave number and in millimeter wave band the eigen-values of A are more and more scattered between both the right and the left half-plane. The efficient implementation of this preconditioned CG (PCG) algorithm is described in details for Complex coefficient matrix. With incomplete factorization preconditioning scheme in the conjugate gradient algorithm, this PCG approach can reach convergence in 20 times CPU time shorter than CG for several typical millimeter wave structures.     

Interleaving Gradient Magnetic Field Method for Diffusion Weighted Spectroscopy

Diffusion-weighted magnetic resonance spectroscopy （DWS） has considerable potential in clinical and research applications. However, it is seldom implemented in conventional magnetic resonance imaging （MRI） scanners due to the strict hardware requirements. We propose an interleaving gradient magnetic field （IGMF） method based on point resolved spectroscopy （PRESS）. Four interlaced powerful diffusion sensitive gradient magnetic fields are positioned around the two π refocusing rf pulses in the PRESS sequence. This method utilizes the interval time in the PRESS pulse sequence, doubles the duration time of the diffusion sensitive gradient magnetic field and decreases the detrimental effect of the induced eddy current. The results of theoretical analysis and experimental observation demonstrate that the IGMF method is suitable for conventional MRI scanners.     

Some Properties of Weighted Tsallis and Kaniadakis Divergences

We are concerned with the weighted Tsallis and Kaniadakis divergences between two measures. More precisely, we find inequalities between these divergences and Tsallis and Kaniadakis logarithms, prove that they are limited by similar bounds with those that limit Kullback–Leibler divergence and show that are pseudo-additive.     

共轭梯度法求解多宗量瞬态热传导反问题

提出瞬态热传导反问题中热物性参数和边界条件多宗量反演的一般数值求解模式,导出了敏度计算公式,并用共轭梯度法求解,探讨了测点数目、测量误差和变量初值对反演结果的影响,并进行了数值验证.     

非均匀采样干涉数据光谱反演技术研究

干涉光谱成像仪获取的干涉数据是一种中间数据,需要进行光谱反演才能够得到目标光谱数据,傅里叶变换方法是常规的光谱反演方法。由于干涉数据中存在非均匀采样问题,若忽略光谱混叠,直接采用快速傅里叶变换会导致反演光谱的失真,难以满足实时处理需求。针对非均匀采样干涉数据的光谱反演需求,将插值及非均匀快速傅里叶变换(NUFFT)方法应用到光谱反演处理中,对过采样及部分欠采样情况下的非均匀采样干涉数据,提出了相应的光谱反演方法,并分析了方法的适用性。最后对过采样及部分欠采样情况下的光谱反演方法进行计算机仿真,过采样情况下采用NUFFT方法反演光谱的精度要明显高于插值方法,而部分欠采样情况下插值方法反演光谱的精度要明显高于NUFFT方法,并对欠采样造成的光谱混叠有一定的修正,验证了方法的有效性。     

Moving Finite Element Methods for a System of Semi-Linear Fractional Diffusion Equations

This paper studies a system of semi-linear fractional diffusion equations which arise in competitive predator-prey models by replacing the second-order derivatives in the spatial variables with fractional derivatives of order less than two. Moving finite element methods are proposed to solve the system of fractional diffusion equations and the convergence rates of the methods are proved. Numerical examples are carried out to confirm the theoretical findings. Some applications in anomalous diffusive Lotka-Volterra and Michaelis-Menten-Holling predator-prey models are studied.     

Numerical-experimental hybrid method for stress separation in digital gradient sensing method

A numerical-experimental hybrid method for the stress separation in the digital gradient sensing (DGS) method is proposed in this study. In the proposed hybrid method, boundary conditions for a local finite element model, that is, nodal force along boundaries are inversely determined from experimental values obtained by the digital gradient sensing method. The hybrid method follows two stages. In stage 1, the DGS method measures the Cartesian stress gradient components directly and, subsequently, the sum in Cartesian stresses at all interesting points on the surface; stress sum are used to compute the unknown boundary conditions for the local model. In stage 2, the individual stress components are calculated by the direct finite element method using the computed boundary conditions from stage 1. The effectiveness is demonstrated by applying the proposed method to a stress concentration problem involving concentrated load acting on an edge of a large planar sheet. The individual stress components thus determined are summed and compared with analytical stress sum, confirming the effectiveness and accuracy of the proposed technique.     

积分梯度法的两种格式

在二维柱坐标系下Lagrange流体力学的计算中,积分梯度法是动量方程的一种有效离散方法.积分梯度法中,IGT(Integral Gradient Total)格式不能保持柱几何下一维球对称性;IGA(Integral Gradient Average)格式可以保持一维球对称性,但当相邻网格质量相差比较大时,会得到远远脱离真实物理现象的加速度.深入研究IGA和IGT格式发现,当相邻网格边界压力取为质量加权时,即使相邻网格质量相差较大,对于一维平面和一维柱问题,IGT与IGA等价;在二维情形下,可以缩小IGT和IGA之间的差异.理论证明,IGA格式不能保持系统的动量守恒,IGT格式能保持系统的动量守恒性.数值模拟结果进一步显示了这两个格式的优缺点.     

Predictions of pulsed field gradient NMR echo-decays for molecules diffusing in various restrictive geometries. Simulations of diffusion propagators based on a finite element method

Pulsed field gradient NMR diffusometry is a promising tool for investigating structures of porous material through determinations of dynamic displacements of molecules in porous systems. A problem with this approach is the lack of closed analytical expressions for echo-decays in anything but idealized pore geometries. We present here an approach based on calculating the appropriate diffusion propagator by means of finite element calculations. The suggested method is quite general, and can be applied to arbitrary porous systems. The protocol for the calculations is outlined and we show results from some different cases: diffusion in confined geometries and in systems that are spatially inhomogeneous with respect to concentration.     

Efficient Proximal Gradient Algorithms for Joint Graphical Lasso

We consider learning as an undirected graphical model from sparse data. While several efficient algorithms have been proposed for graphical lasso (GL), the alternating direction method of multipliers (ADMM) is the main approach taken concerning joint graphical lasso (JGL). We propose proximal gradient procedures with and without a backtracking option for the JGL. These procedures are first-order methods and relatively simple, and the subproblems are solved efficiently in closed form. We further show the boundedness for the solution of the JGL problem and the iterates in the algorithms. The numerical results indicate that the proposed algorithms can achieve high accuracy and precision, and their efficiency is competitive with state-of-the-art algorithms.     

负离子源引出系统束光学的数值模拟

根据在表面-等离子体型负离子源中,负离子主要由转换电极表面产生这一特点,参考正离子源引出系统数值计算程序,建立了负离子源引出系统的数值计算模型和计算程序;对表面-等离子体型桶式负离子源引出系统束光学的性质进行了数值模拟。对计算结果的检验和分析表明,这个计算模型和程序反映了负离子源引出系统束光学的基本特性。     

含噪声包裹相位图的加权最小二乘相位展开算法研究

二维相位展开广泛应用在精密光学测量、自适应光学、合成孔径雷达、图像处理等领域中。为处理含噪声包裹相位图,以预条件共轭斜量法求解权重最小二乘相位展开方程。引入非加权二维离散余弦变换求解泊松方程得到的最小二乘相位解作为共轭斜量法的初始解,从而加快了收敛速度,同时提出一种新质量图确定算法求解过程中的权重项。计算机模拟和试验表明算法计算速度快,能有效地消除传统路径积分法在处理信噪比低包裹相位图时的"拉线"现象,是一种有效的相位展开方法。     

孔径平均对自由空间光通信误码率的影响

研究了接收端孔径平均效应对自由空间光通信误码率的影响。对数辐射方差是反映接收孔径上某一点闪烁起伏的参量。如果自由空间光通信的误码率仅仅是和对数辐射方差有关,那么得到的就是点接收器的结果。对数辐射的协方差是与接收平面上某两点间距离有关的闪烁起伏参量。孔径平均效应意味着要将误码率和对数辐射的协方差联系起来。文中共考虑了两种模型:球面波和平面波,并假设接收孔径直径等于第一菲涅耳区尺度。分析的结果显示,孔径平均后的误码率有明显的减小。     

Finite element method for the two-dimensional atmospheric radiative transfer

The finite element method is applied to the solution of the two-dimensional atmospheric radiative transfer. The analysis is mainly focussed on the derivation of the cell or element equation. The Galerkin method and several hybrid methods using the integral and finite difference form of the radiative transfer equation are employed to obtain the cell equation. The assembled system of equations relating the radiances at the lower and upper boundary of the domain is solved by a direct method.     

高频腔本征值问题的三维数值模拟

 运用有限积分法软件MAFIA和有限元软件ANSYS对两类高频腔的本征值问题进行数值模拟，将模拟结果与部分实验数据进行验证，得到了比较一致的结果，其平均计算精度在5%以内，证实了两种电磁软件都可用于高频腔本征值问题的数值模拟，在计算精度上为满足工程设计需要做了一些研究。另外将两种模拟软件在计算时间上作了对比。     

求解中子动力学方程的加权蒙特卡罗方法

为了实现基于蒙特卡罗方法的中子动力学计算,在传统的直接蒙特卡罗动力学方法的基础上,提出了一种加权蒙特卡罗动力学方法。该方法通过引入粒子权重的概念,隐式考虑中子俘获反应和裂变反应过程中中子数目的变化,避免了模拟粒子的数目随时间的变化,降低了统计偏差,消除了程序计算过程中粒子的存库操作,提高了计算精度。基于单能点堆模型,开发了中子动力学计算程序NECP-Dandi,进行了大量数值验证与分析,包括无缓发中子、单组缓发中子、六组缓发中子、正阶跃反应性引入、负阶跃反应性引入、正脉冲反应性、负脉冲反应性和正线性反应性引入等情况。数值结果表明,相比于直接蒙特卡罗动力学方法,加权蒙特卡罗动力学方法在计算结果的精度和计算效率上有较为明显的改进,程序结构更为简洁。     

Two-Grid Discretization Scheme for Nonlinear Reaction Diffusion Equation by Mixed Finite Element Methods

In this paper, we study an efficient scheme for nonlinear reaction-diffusion equations discretized by mixed finite element methods. We mainly concern the case when pressure coefficients and source terms are nonlinear. To linearize the nonlinear mixed equations, we use the two-grid algorithm. We first solve the nonlinear equations on the coarse grid, then, on the fine mesh, we solve a linearized problem using Newton iteration once. It is shown that the algorithm can achieve asymptotically optimal approximation as long as the mesh sizes satisfy $H =\mathcal{O}(h^{\frac{1}{2}})$. As a result, solving such a large class of nonlinear equations will not be much more difficult than getting solutions of one linearized system.     

共轭梯度算法在画幅式层析成像光谱仪中的应用

介绍了画幅式层析成像光谱仪的基本原理及共轭梯度算法,解决了共轭梯度法在实际应用中的几个问题,如图像分割,先验算子的引入.通过对实际数据的计算,得到了95×95的像素,28个波段的数据立方体.实际证明本算法及程序在层析成像光谱仪上的应用是成功的.