粒子输运方程的线性间断有限元方法

将空间线性间断有限元方法应用于动态粒子输运方程的求解.数值算例表明,空间线性间断有限元方法在网格边界的数值精度方面明显高于指数格式和菱形格式,并且通量在时间上的微分曲线相对光滑,避免了指数格式、菱形格式数值解的非物理振荡现象.     

耦合FE/WB法在声分析中的应用

简要描叙FE法(finite element method)和WB法(wave based method)的理论背景以及耦合FE/WB法的数学基础.耦合FE/WB法利用两者的优势——FE法的广泛应用和WB法的高收敛特性,将FE模型中较大且几何简单的部分采用WB法代替.耦合模型具有相对较少的自由度.对于较高的频率还可以进行细分得到更高的计算精度,并利用模态缩减法进一步减少自由度数.数值算例结果表明,该耦合方法有潜力覆盖中频段的声分析.     

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.     

Concentration of Laplace Eigenfunctions and Stabilization of Weakly Damped Wave Equation

In this article, we prove some universal bounds on the speed of concentration on small (frequency-dependent) neighbourhoods of sub-manifolds of L²-norms of quasi modes for Laplace operators on compact manifolds. We deduce new results on the rate of decay of weakly damped wave equations.     

Finite Difference/Element Method for a Two-Dimensional Modified Fractional Diffusion Equation

We present the finite difference/element method for a two-dimensional modified fractional diffusion equation. The analysis is carried out first for the time semi-discrete scheme, and then for the full discrete scheme. The time discretization is based on the $L1$-approximation for the fractional derivative terms and the second-order backward differentiation formula for the classical first order derivative term. We use finite element method for the spatial approximation in full discrete scheme. We show that both the semi-discrete and full discrete schemes are unconditionally stable and convergent. Moreover, the optimal convergence rate is obtained. Finally, some numerical examples are tested in the case of one and two space dimensions and the numerical results confirm our theoretical analysis.     

一维多介质可压缩Euler方程的高精度RKDG有限元方法

采用RKDG有限元目的、Level Set目的和改进的带"Isentropic"修正的Ghost Fluid目的模拟了一维多介质可压缩Euler方程,其中Euler方程、Level Set方程和重新初始化方程都采用了三阶精度的RKDG有限元目的进行离散,并对一维两种介质可压缩流体进行了数值实验,得到了较高分辨率的计算结果.     

A Mixed Finite Element Method for Stationary Magneto-Heat Coupling System with Variable Coefficients

In this article, a mixed finite element method for thermally coupled, stationary incompressible MHD problems with physical parameters dependent on temperature in the Lipschitz domain is considered. Due to the variable coefficients of the MHD model, the nonlinearity of the system is increased. A stationary discrete scheme based on the coefficients dependent temperature is proposed, in which the magnetic equation is approximated by Nédélec edge elements, and the thermal and Navier–Stokes equations are approximated by the mixed finite elements. We rigorously establish the optimal error estimates for velocity, pressure, temperature, magnetic induction and Lagrange multiplier with the hypothesis of a low regularity for the exact solution. Finally, a numerical experiment is provided to illustrate the performance and convergence rates of our numerical scheme.     

辐射传输方程的三阶球谐展开(P3)近似的有限元法求解

从辐射传输方程出发,应用球谐函数方法对辐射率进行多项式展开,并利用球谐函数的正交性和递推性,在二维笛卡儿坐标下,导出了三阶展开的球谐函数微分方程组(P3近似),改进了以往成像文献中忽略各项异性因子的P3近似,并用有限元方法对二维圆域均匀和非均匀两种情况做了数值模拟.与漫射近似模型相比较,P3近似能更准确地描述光源附近及...     

辐射传输方程的三阶球谐展开(P3)近似的有限元法求解

从辐射传输方程出发,应用球谐函数方法对辐射率进行多项式展开,并利用球谐函数的正交性和递推性,在二维笛卡儿坐标下,导出了三阶展开的球谐函数微分方程组（P3近似）,改进了以往成像文献中忽略各项异性因子的P3近似,并用有限元方法对二维圆域均匀和非均匀两种情况做了数值模拟.与漫射近似模型相比较,P3近似能更准确地描述光源附近及吸收较强情况下边界的光辐射分布情况.     

求解洞穴散射问题的带有吸收边界层的有限元方法

给出二维洞穴散射问题的带有吸收边界层的有限元方法.通过吸收边界层将无界域问题截断,得到一个有界域问题,求解此有界域问题来代替原问题.进一步分析收敛性并进行数值模拟.结果表明了此方法的有效性.     

间断有限元方法求解一维非平衡辐射扩散方程

研究一维非平衡辐射扩散方程的数值方法.通过求解间断系数热传导方程的广义黎曼问题,得到一种带加权数值流量,基于该数值流量构造了一类新型的间断有限元方法.在时间离散上采用向后Euler方法,形成的非线性方程组采用Picard迭代求解.数值试验表明该方法具有捕捉大梯度的能力,而且能适应扩散系数间断的情形.     

A Comparative Study of Finite Element and Finite Difference Methods for Two-Dimensional Space-Fractional Advection-Dispersion Equation

The paper makes a comparative study of the finite element method (FEM) and the finite difference method (FDM) for two-dimensional fractional advection-dispersion equation (FADE) which has recently been considered a promising tool in modeling non-Fickian solute transport in groundwater. Due to the non-local property of integro-differential operator of the space-fractional derivative, numerical solution of FADE is very challenging and little has been reported in literature, especially for high-dimensional case. In order to effectively apply the FEM and the FDM to the FADE on a rectangular domain, a backward-distance algorithm is presented to extend the triangular elements to generic polygon elements in the finite element analysis, and a variable-step vector Grünwald formula is proposed to improve the solution accuracy of the conventional finite difference scheme. Numerical investigation shows that the FEM compares favorably with the FDM in terms of accuracy and convergence rate whereas the latter enjoys less computational effort.     

Simulation of Electromagnetic Wave Logging Response in Deviated Wells Based on Vector Finite Element Method

The vector finite element method of tetrahedral elements is used to model 3D electromagnetic wave logging response. The tangential component of the vector field at the mesh edges is used as a degree of freedom to overcome the shortcomings of node-based finite element methods. The algorithm can simulate inhomogeneous media with arbitrary distribution of conductivity and magnetic permeability. The electromagnetic response of well logging tools are studied in dipping bed layers with the borehole and invasion included. In order to simulate realistic logging tools, we take the transmitter antennas consisting of circular wire loops instead of magnetic dipoles. We also investigate the apparent resistivity of inhomogeneous formation for different dip angles.     

微波元件与加速结构的有限元分析

研究了有限元方法在二维均匀结构和轴对称加速结构中的应用,采用了带宽优化技术和子空间求解特征值方法,并给出了部分例子,结果表明该方法精度高、速度快.这些方法可直接用于三维程序、并行计算.     

Parameters Extraction of Millimeter Wave Circuits by Hybrid Vector Finite Element and Moment Method Combined with SOC Technique

This paper presents a hybrid method, which couples the vector finite element method (FEM) and method of moment (MOM) for analyzing the field and current distribution of the millimeter wave circuits. The FEM is applied to handle the interior region of dielectric bodies and MOM is used to solve surface integral equations. Then, These integral expressions are coupled into the FEM equations through the continuity of the tangential fields across the connection boundaries. Simultaneously, the short-open calibration (SOC) technique is used for predicting accurately the scattering parameters of the circuits. Numerical results are well compared with those published in the previous literatures.     

An Inverse Finite Element Method for Pure and Binary Solidification Problems

A 2D axisymmetric formulation for the solution of a directional solidification problem using an inverse finite-element method (IFEM) is presented. An algorithm developed by A. N. Alexandrou (Int. J. Numer. Methods Eng.28, 2383, 1989) has been modified and extended to include more general boundary conditions. The latter includes the explicit presence of an ampoule (with a complex shape) that contains the solid and the melt from which it is growing. Heat transfer between the ampoule and the external environment, time-dependent thermal boundary conditions, nonmonotonic temperature distributions, and species diffusion in the melt and crystal are also taken into account. Thus, our extended formulation encompasses a wider class of solidification problems than previous IFEM methods. Numerical experiments that illustrate the suitability of the extended IFEM are presented. In particular, we present a simulation of the directional solidification of zinc cadmium telluride using boundary conditions corresponding to an actual experiment scenario.     

Jacobian Elliptic Function Method and Solitary Wave Solutions for Hybrid Lattice Equation

In this paper, we have successfully extended the Jacobian elliptic function expansion approach to nonlinear differential-difference equations. The Hybrid lattice equation is chosen to illustrate this approach. As a consequence, twelve families of Jacobian elliptic function solutions with different parameters of the Hybrid lattice equation are obtained. When the modulus m→1 or 0, doubly-periodic solutions degenerate to solitonic solutions and trigonometric function solutions, respectively.     

A Posteriori Error Estimates of a Combined Mixed Finite Element and Discontinuous Galerkin Method for a Kind of Compressible Miscible Displacement Problems

A kind of compressible miscible displacement problems which include molecular diffusion and dispersion in porous media are investigated. The mixed finite element method is applied to the flow equation, and the transport one is solved by the symmetric interior penalty discontinuous Galerkin method. Based on a duality argument, employing projection estimates and approximation properties, a posteriori residual-type $hp$ error estimates for the coupled system are presented, which is often used for guiding adaptivity. Comparing with the error analysis carried out by Yang (Int. J. Numer. Meth. Fluids, 65(7) (2011), pp. 781-797), the current work is more complicated and challenging.     

Analysis of Two-Grid Methods for Nonlinear Parabolic Equations by Expanded Mixed Finite Element Methods

In this paper, we present an efficient method of two-grid scheme for the approximation of two-dimensional nonlinear parabolic equations using an expanded mixed finite element method. We use two Newton iterations on the fine grid in our methods. Firstly, we solve an original nonlinear problem on the coarse nonlinear grid, then we use Newton iterations on the fine grid twice. The two-grid idea is from Xu$'$s work [SIAM J. Numer. Anal., 33 (1996), pp. 1759-1777] on standard finite method. We also obtain the error estimates for the algorithms of the two-grid method. It is shown that the algorithm achieves asymptotically optimal approximation rate with the two-grid methods as long as the mesh sizes satisfy $h=\mathcal{O}(H^{(4k+1)/(k+1)})$.