A Multistep Finite Volume Method with Penalty for Thermal Convection Problems in Two or Three Dimensions

The author considers a thermal convection problem with infinite Prandtl number in two or three dimensions. The mathematical model of such problem is described as an initial boundary value problem made up of three partial differential equations. One equation of the convection-dominated diffusion type for the temperature, and another two of the Stokes type for the normalized velocity and pressure. The approximate solution is obtained by a penalty finite volume method for the Stokes equation and a multistep upwind finite volume method for the convection-diffusion equation. Under suitable smoothness of the exact solution, error estimates in some discrete norms are derived.     

L ∞-estimates of mixed finite element methods for general nonlinear optimal control problems

This paper investigates L~∞-estimates for the general optimal control problems governed by two-dimensional nonlinear elliptic equations with pointwise control constraints using mixed finite element methods.The state and the co-state are approximated by the lowest order Raviart-Thomas mixed finite element spaces and the control is approximated by piecewise constant functions.The authors derive L~∞-estimates for the mixed finite element approximation of nonlinear optimal control problems.Finally,the numerical examples are given.     

Finite element approximation for a class of parameter estimation problems

This paper investigates the finite element approximation of a class of parameter estimation problems which is the form of performance as the optimal control problems governed by bilinear parabolic equations,where the state and co-state are discretized by piecewise linear functions and control is approximated by piecewise constant functions.The authors derive some a priori error estimates for both the control and state approximations.Finally,the numerical experiments verify the theoretical results.     

Efficient analysis of dielectric radomes using multilevel fast multipole algorithm with CRWG basis

A full-wave analysis of the electromagnetic problem of a three-dimensional （3-D） antenna radiating through a 3-D dielectric radome is preserued. The problem is formulated using the Poggio-Miller-Chang-Harrington- Wu（PMCHW） approach for homogeneous dielectric objects and the electric field integral equation for conducting objects. The integral equations are discretized by the method of moment （MoM）, in which the conducting and dielectric surface/interfaces are represented by curvilinear triangular patches and the unknown equivalent electric and magnetic currents are expanded using curvilinear RWG basis functions. The resultant matrix equation is then solved by the multilevel fast multipole algorithm （MLFMA） and fast far-field approximation （FAFFA） is used to further accelerate the computation. The radiation patterns of dipole arrays in the presence of radomes are presented. The numerical results demonstrate the accuracy and versatility of this method.     

Variational discretization for optimal control problems governed by parabolic equations

This paper considers the variational discretization for the constrained optimal control problem governed by linear parabolic equations. The state and co-state are approximated by Raviart- Thomas mixed finite element spaces, and the authors do not discretize the space of admissible control but implicitly utilize the relation between co-state and control for the discretization of the control. A priori error estimates are derived for the state, the co-state, and the control. Some numerical examples are presented to confirm thetheoretical investigations.     

RESIDUAL A POSTERIORI ERROR ESTIMATE OF A NEW TWO-LEVEL METHOD FOR STEADY NAVIER-STOKES EQUATIONS

Residual-based a posteriori error estimate for conforming finite element solutions of incom-pressible Navier-Stokes equations,which is computed with a new two-level method that is differentfrom Volker John,is derived.A posteriori error estimate contains additional terms in comparison tothe estimate for the solution obtained by the standard finite element method.The importance of theadditional terms in the error estimates is investigated by studying their asymptotic behavior.For opti-mal scaled meshes,these bounds are not of higher order than the convergence of discrete solution.Thetwo-level method aims to solve the nonlinear problem on a coarse grid with less computational work,then to solve the linear problem on a fine grid,which is superior to the usual finite element methodsolving a similar nonlinear problem on the fine grid.     

A discontinuous Galerkin method combined with mixed finite element for seawater intrusion problem

<正> Seawater intrusion problem is considered in this paper.Its mathematical model is anonlinear coupled system of partial differential equations with initial boundary problem.It consistsof the water head equation and the salt concentration equation.A combined method is developedto approximate the water head equation by mixed finite element method and concentration equationby discontinuous Galerkin method.The scheme is continuous in time and optimal order estimates inH~1-norm and L~2-norm are derived for the errors.     

FINITE ELEMENT ERROR EXPANSIONS FOR THE EIGENVALUE APPROXIMATION TO THE MULTIGROUP DIFFUSION EQUATIONS

Asymptotic expansions for the finite element approximation to the eigenvalues of the mult-igroup diffusion equations on Ω(?)R~n(n=2,3)of reactor theory are given,firstly for a piecewiseuniform triangulation,then for nonuniform quadrilateral meshes,and finally for nonuniformhexahedral meshes.The effect of certain classes of numerical integrations is studied.As applicationsof the expansions,several extrapolation formulas and a posteriori error estimates are obtained.     

The Necessary and Sufficient Conditions of Exponential Stability for Heat and Wave Equations with Memory

This paper is addressed to a study of the stability of heat and wave equations with memory.The necessary and sufficient conditions of the exponential stability are investigated by the theory of Laplace transform. The results show that the stability depends on the decay rate and the coefficient of the kernel functions of the memory. Besides, the feedback stabilization of the heat equation is obtained by constructing finite dimensional controller according to unstable eigenvalues. This stabilizing...     

求解非均匀介质问题的双共轭梯度方法

分析非均匀介质条件下的电成像问题 ,在说明电成像仪的测量环境和测量原理之后 ,对这种复杂条件的电磁场问题 ,采用三维有限元方法进行分析。为了保证计算精度 ,在分析过程中 ,需要划分较多的空间网格 ,从而生成大型的有限元矩阵方程。对此大型矩阵方程 ,采用计算效率较高的双共轭梯度方法求解 ,给出了双共轭梯度方法的算法和利用该算法求解有限元矩阵方程时的收敛速度曲线 ,并对水平分层和倾斜分层两种典型情况下的非均匀介质成像问题进行分析 ,给出了模拟测量成像结果。     

INTERPOLATED FINITE ELEMENT METHODS FOR SECOND ORDER HYPERBOLIC EQUATIONS AND THEIR GLOBAL SUPERCONVERGENCE

We develop the interpolated finite element method to solve second-order hy-perbolic equations. The standard linear finite element solution is used to generate a newsolution by quadratic interpolation over adjacent elements. We prove that this interpo-lated finite element solution has superconvergence. This method can easily be applied togenerating more accurate gradient either locally or globally, depending on the applications.This method is also completely vectorizable and parallelizable to take the advantages ofmodern computer structures. Several numerical examples are presented to confirm ourtheoretical analysis.     

AN EXTRAPOLATION METHODFOR FINITE ELEMENT APPROXIMATION OF INTEGRODIFFERENTIAL EQUATIONS WITH PARAMETERS

In this paper,we discuss the accelerating convergence method for finite elementapproximation of integro-differential equations with parameters.As applications,we give theerror estimates of finite element for the first kind of Fredholm integral equation,particularlyfor the Volterra integral equation with kernel condition k(x,x)=0.     

Vacuum Fluctuations, Isolated Systems and Hierarchical Cosmology

1.THEAETHER,UNITS,ANDARBrrRARINESSOFGEOMETaY1.1TheAetherNullresultsofaetherdriftexperimentsdonotprovethenon--existenceofanaether.TheMichelson--Moneyexperimentcomparesround--triptransittimesovertwopathsoflengthsLllandLIwhicharerespectivelyparallelandperpendiculartoaethervelocityvg.AlongLlltheexpectedtimewasLll/(c v,) Lll/(c--v,)andalongLIitwasZL./c(l--vg/.')"'.Introducingacv./candy.=(l--p:)"',theexpectedtimesarerespectivelyZyUlj/cand.Zy/II/c.ThenullresultwasexplainedbyLor…     

高次有限元方程的一种并行预条件子

偏微分方程(PDEs)是大规模科学工程计算与数值仿真中的基本数学模型,有限元方法是数值求解偏微分方程的一类重要的离散化方法,高次有限元又是其中的一类常用有限元。针对一类基本的PDE模型(Poission方程)的高次有限元方程,设计了一种基于辅助变分问题的并行预条件子,并从理论上严格证明了该预条件子的条件数的一致有界性,数值实验验证了理论结果的正确性及相应预条件共轭梯度(PCG)法的高效性和鲁棒性。     

矢量波动方程的直接差分解法

对电磁场波动方程进行了分析,研究了一种基于电场波动方程的时域数值方法———矢量波动方程的直接差分解法,给出了电场和磁场的计算公式。采用线性插值方法获得了导体表面法向电场的计算公式,给出了算法的稳定性条件。与时域有限差分法相比,由于该方法中电场分量的迭代更新无需磁场参与,从而使算法的复杂性大大降低,是一种具有良好应用前景的电磁场时域数值分析方法,最后通过一个数值计算实例证明了该算法的有效性。     

基于有限单元的贝叶斯滤波估计算法

现代工程系统具有较强的非线性特性,针对这类非线性系统的状态估计问题,提出基于有限单元的贝叶斯原理估计的非线性滤波方法。采用有限单元法逼近系统状态的先验概率解,即前向Kolmogorov方程的解,通过贝叶斯估计得到状态的后验信息。将其方法应用到惯性/地形组合导航系统中,仿真结果表明该方法的可行性。     

解析与数值的意义及所描述的世界

基本科学规律一般指动力学规律即状态时间变化的微观规律,呈现为“边界条件＋动力学方程”形式。动力学方程是潜无限意义下的微分方程,少数简单情况下有解析解,初等连续函数形式的解析解描述的是一个确定性、连续性的世界;多数情况下没有解析解,那是一个不确定性和非线性放大机制并存的世界;潜无限动力学方程不能完全描述的是实无限的连续性世界。在运动学层面,无论有无解析解也无论存不存在动力学方程,都可以用数值方法研究动力学现象的运动学行为,数值解描述的是有限性、确定性、离散性、历史性的世界。