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.     

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

干涉光谱成像仪获取的干涉数据是一种中间数据,需要进行光谱反演才能够得到目标光谱数据,傅里叶变换方法是常规的光谱反演方法。由于干涉数据中存在非均匀采样问题,若忽略光谱混叠,直接采用快速傅里叶变换会导致反演光谱的失真,难以满足实时处理需求。针对非均匀采样干涉数据的光谱反演需求,将插值及非均匀快速傅里叶变换(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.     

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.     

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

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

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

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

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.     

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

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

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.     

渐变折射率波导的BPM分析

利用BPM对MEMS波分复用系统中的渐变折射率波导的性能进行分析。详细推导了BPM的微分方程,并给出了BPM微分方程的有限差分解的形式。并根据FDBPM对已知折射率分布的径向渐变折射率波导进行数值计算分析,其结果表明,在波长为1.55μm的高斯光束垂直入射的情况下,中心折射率为1.5901,聚焦参数为0.052μm-1,并且半径为3μm的平方律渐变波导中具有快速稳定的、高传输率的高斯基模分布。     

A Fast Finite Element Analysis of Millimeter Wave Waveguide Filters

The multifrontal method is applied for solving a large system of linear equations resulting from the use of edge-based finite element method (FEM). The finite element method combined with the perfectly matched layers (PML) is given for simulation of microwave filters and the algorithm of multifrontal method is described in detail. The reflection and insertion losses of millimeter wave filter are analyzed as the examples and the obtained results are compared with those obtained from literature. In order to demonstrate the efficiency of the multifrontal method, the computational time is compared with that of both symmetric successive overrelaxation (SSOR) preconditioned conjugate gradient (PCG) and conjugate gradient methods (CG) for the thick-iris waveguide bandpass filer.     

权基本无振荡格式求解图像恢复问题

详细描述权基本无振荡格式解决基于全变分方法的图像恢复问题,考虑高斯核和不同的边界条件.采用权基本无振荡格式、维纳算法、权基本无振荡格式和维纳算法的联合三种方法对恢复问题进行数值试验,并对结果进行比较.试验结果显示权基本无振荡格式和维纳算法的联合方法的有效性.     

A three-dimensional finite element approach for predicting the transmission loss in mufflers and silencers with no mean flow

A three-dimensional finite element method has been implemented to predict the transmission loss of a packed muffler and a parallel baffle silencer for a given frequency range. Iso-parametric quadratic tetrahedral elements have been chosen due to their flexibility and accuracy in modeling geometries with curved surfaces. For accurate physical representation, perforated plates are modeled with complex acoustic impedance while absorption linings are modeled as a bulk media with a complex speed of sound and mean density. Domain decomposition and parallel processing techniques are applied to address the high computational and memory requirements. The comparison of the computationally predicted and the experimentally measured transmission loss shows a good agreement.     

基于梯度翅片换热器的热电发电系统性能研究

本文建立了热电发电系统(TEG)多物理场数值模型,并充分考虑换热器流体影响,综合研究了具有不同热侧换热器翅片结构的TEG系统性能。在雷诺数为1000~10000范围内,分析了流体沿程温度分布特征、泵功及热电发电模块的能量转换特性.所研究的三种翅片结构包括:全流道等高度直翅片(Fin-1)、下游强化梯度翅片(Fin-2)以及上游强化梯度翅片(Fin-3).研究表明,通道长高比及热电材料覆盖率一定,热电发电功率及转换效率随流量呈二次曲线变化关系,存在最匹配流量使得系统发电性能最佳。等高度直翅片对流量的变化敏感,随流量增大,则压损增大,导致系统净输出功率及发电效率无收益.而梯度翅片可以在更大范围内产生正收益;下游强化梯度翅片具有最佳的流体沿程温度均匀性,但沿程局部热阻却最大.综合考虑沿程局部热阻分布及泵功消耗,上游强化梯度翅片TEG系统净转换效率最高,因此局部热阻分布及泵功综合因素应为TEG内的换热器合理设计的关键。     

不同方法测量折射率结构常数的比较

对几种测量折射率结构常数（Ｃ＾２ｎ）的方法进行了分析和比较,指出这几种方法的测量结果有时不一致的原因。结果表明,和双点法测量相比,温度谱幂率的变化对闪烁法测量结果造成较大偏差,而对到达角起伏测量和声雷达测量结果差别较小;外悄度发迹对闪烁法测量几乎没有影响,但当外尺度变小时,利用双点法测量,到达角起伏测量和声雷达测量结果都有较大的影响。给出了双点测量Ｃ＾２ｎ和外尺度的关系的实验验证,并给出了两固定点     

The Gradient Condition for One-Dimensional Symmetric Exclusion Processes

For every Gibbs measure on the one dimensional lattice Z with translation-invariant potential of finite range, an exchange rate for one-dimensional lattice gas which satisfy both the detailed balance condition relative to the Gibbs measure and the gradient condition is constructed.     

The finite element method for micro-scale modeling of ultrasound propagation in cancellous bone

Quantitative ultrasound for bone assessment is based on the correlations between ultrasonic parameters and the properties (mechanical and physical) of cancellous bone. To elucidate the correlations, understanding the physics of ultrasound in cancellous bone is demanded. Micro-scale modeling of ultrasound propagation in cancellous bone using the finite-difference time-domain (FDTD) method has been so far utilized as one of the approaches in this regard. However, the FDTD method accompanies two disadvantages: staircase sampling of cancellous bone by finite difference grids leads to generation of wave artifacts at the solid–fluid interface inside the bone; additionally, this method cannot explicitly satisfy the needed perfect-slip conditions at the interface. To overcome these disadvantages, the finite element method (FEM) is proposed in this study. Three-dimensional finite element models of six water-saturated cancellous bone samples with different bone volume were created. The values of speed of sound (SOS) and broadband ultrasound attenuation (BUA) were calculated through the finite element simulations of ultrasound propagation in each sample. Comparing the results with other experimental and simulation studies demonstrated the capabilities of the FEM for micro-scale modeling of ultrasound in water-saturated cancellous bone.