FINITE ELEMENT METHOD FOR INCOMPRESSIBLE MISCIBLE DISPLACEMENT IN POROUS MEDIA

Under the assumptions of nonlinear finite element and △t=o(h),Ewing and Wheeler discussed a Galerkin method for the single phase incompressible miscible displacement of one fluid by another in porous media.In the present paper we give a finite element scheme which weakens the △t=o(h)-restriction to △t=o(h~ε),0<ε≤1/2.Furthermore,this scheme is suitable for both linear element and nonlinear element.We also derive the optimal approximation estimates for concentration c,its gradient ▽c and the gradient ▽p of the pressure p.     

THE CONVERGENCE OF A DOMAIN DECOMPOSITION ALGORITHM FOR PARABOLIC PROBLEMS

At recent, Hourgat et gave a domain decomposition algorithm for elliptic problems which can be implemented in parallel. Many numerical experiments have illustrated its efficiency. In the present paper, we apply this algorithm to solve the discrete parabolic problems, analyse its convergence and show that its convergence rale is about (1 - 2p + σp2 ) which is nearly optimal and independent of the parameter τ, where σ τ O((1 +H )(1 + ln(H / h))2 ). 0 < p < 1 / σ,τ,h,H are the time step size, finite element parameter and subdomain diameter, respectively.     

有约束多元回归分析中的C_p统计量

In this paper, Mallow's C_p statistic in multiple linear regression analysis is genera-lized to the case with linear equality constrained conditions.The C_(p,m_1) statistic is introduced asC_(p,m_1)=RS_p~2/σ~2-m 2(p-m_1),where m and m_1 are respectively the numbers of the observed data of the variables andof the constrained conditions, p>m_1 is the number of the remaining regression variablesafter q variables are rejected from the n regression variables (p q=n) RS_p~2 is the sumof squares of residuals of the regression for the p remaining variables, σ~2 is an unbiasedestimate of variance σ~2. The properties of the C_(p,m_1) statistic and a method for selectingand analyzing variables in constrained multiple linear regression analysis by using theC_(p,m_1) statistic are discussed in detail.     

High accuracy nonconforming finite elements for fourth order problems

The approach of nonconforming finite element method admits users to solve the partial differential equations with lower complexity,but the accuracy is usually low.In this paper,we present a family of highaccuracy nonconforming finite element methods for fourth order problems in arbitrary dimensions.The finite element methods are given in a unified way with respect to the dimension.This is an effort to reveal the balance between the accuracy and the complexity of finite element methods.     

THE OPTIMAL ORDER ERROR ESTIMATES FOR FINITE ELEMENT APPROXIMATIONS TO HYPERBOLIC PROBLEMS

In this paper, the linear finite element approximation to the positive and symmetric,linear hyperbolic systems is analyzed and an O(h^2) order error estimate is established under the conditions of strongly regular triangulation and the H^3-regularity for the exact solutions. The convergence analysis is based on some superclose estimates derived in this paper. Our method and result here are also applicable to general hyperbolic problems.Finally, we discuss the linearized shallow water system of equations.     

Exponential polynomials as solutions of certain nonlinear difference equations

Recently, C.-C. Yang and I. Laine have investigated finite order entire solutions f of nonlinear differential-difference equations of the form fn + L(z, f ) = h, where n ≥ 2 is an integer. In particular, it is known that the equation f(z)2 + q(z)f (z + 1) = p(z), where p(z), q(z) are polynomials, has no transcendental entire solutions of finite order. Assuming that Q(z) is also a polynomial and c ∈ C, equations of the form f(z)n + q(z)e Q(z) f(z + c) = p(z) do posses finite order entire solutions. A classification of these solutions in terms of growth and zero distribution will be given. In particular, it is shown that any exponential polynomial solution must reduce to a rather specific form. This reasoning relies on an earlier paper due to N. Steinmetz.     

SURFACE FINITE ELEMENTS FOR

In this article we define a surface finite element method （SFEM） for the numerical solution of parabolic partial differential equations on hypersurfaces F in R^n＋1. The key idea is based on the approximation of F by a polyhedral surface Гh consisting of a union of simplices （triangles for n = 2, intervals for n = 1） with vertices on Г. A finite element space of functions is then defined by taking the continuous functions on Гh which are linear affine on each simplex of the polygonal surface. We use surface gradients to define weak forms of elliptic operators and naturally generate weak formulations of elliptic and parabolic equations on Г. Our finite element method is applied to weak forms of the equations. The computation of the mass and element stiffness matrices are simple and straightforward. We give an example of error bounds in the case of semi-discretization in space for a fourth order linear problem. Numerical experiments are described for several linear and nonlinear partial differential equations. In particular the power of the method is demorrstrated by employing it to solve highly nonlinear second and fourth order problems such as surface Allen-Cahn and Cahn-Hilliard equations and surface level set equations for geodesic mean curvature flow.     

Three-Dimensional Finite Element Superconvergent Gradient Recovery on Par6 Patterns

<正>In this paper,we present a theoretical analysis for linear finite element superconvergent gradient recovery on Par6 mesh,the dual of which is centroidal Voronoi tessellations with the lowest energy per unit volume and is the congruent cell predicted by the three-dimensional Gersho's conjecture.We show that the linear finite element solution u_h and the linear interpolation u_1 have superclose gradient on Par6 meshes. Consequently,the gradient recovered from the finite element solution by using the superconvergence patch recovery method is superconvergent to▽u.A numerical example is presented to verify the theoretical result.     

ON THE FINITE VOLUME ELEMENT VERSION OF RITZ-VOLTERRA PROJECTION AND APPLICATIONS TO RELATED EQUATIONS

In this paper, we present a general error analysis framework for the finite volume element (FVE) approximation to the Ritz-Volterra projection, the Sobolev equations and parabolic integro-differential equations. The main idea in our paper is to consider the FVE methods as perturbations of standard finite element methods which enables us to derive the optimal L2 and H1 norm error estimates, and the L∞ and W∞1 norm error estimates by means of the time dependent Green functions. Our disc ussions also include elliptic and parabolic problems as the special cases.     

Vector subdivision schemes in (L_p(R~5))~r(1≤p≤∞) spaces

The purpose of this paper is to investigate the refinement equations of the formwhere the vector of functions = (1, … ,r)T is in (LP(R8))T,1 ≤ p ≤∞, α(α),α ∈ Z5, is a finitely supported sequence of r × r matrices called the refinement mask, and M is an s x a integer matrix such that limn→ ∞ M-n = 0. In order to solve the refinement equation mentioned above, we start with a vector of compactly supported functions (0 ∈ (LP(R8))r and use the iteration schemes fn := Qan0,n = 1,2,…, where Qa is the linear operator defined on (Lp(R8))r given byThis iteration scheme is called a subdivision scheme or cascade algorithm. In this paper, we characterize the Lp-convergence of subdivision schemes in terms of the p-norm joint spectral radius of a finite collection of some linear operators determined by the sequence a and the set B restricted to a certain invariant subspace, where the set B is a complete set of representatives of the distinct cosets of the quotient group Z8/MZ8 containing 0.     

Ancora sulla permanenza della struttura lacunare nel prolungamento analitico

Annali di Matematica Pura ed Applicata (1923 -) - Si considerano le serie di potenze f(z;0) lacunari pure con la successione di lacune {ph, qh}, (h=1, 2, 3, ...), qh/ph →1 + Λ, (0...     

A Finite—dimensional Integrable System with Time Parameter

ＡＦｉｎｉｔｅ－ｄｉｍｅｎｓｉｏｎａｌＩｎｔｅｇｒａｂｌｅＳｙｓｔｅｍｗｉｔｈＴｉｍｅＰａｒａｍｅｔｅｒＺｈａｎｇＪｉｎｓｈｕｎ（张金顺）；ＷａｎｇＨｏｎｇｙｅ（王鸿业）（ＤｅｐａｒｔｍｅｎｔｏｆＭａｔｈｅｍａｔｉｃｓ，ＺｈｅｎｇｚｈｏｕＵｎｉｖｅｒ...     

Navier-Stokes方程的最佳有限元非线性Galerkin算法

1．引言对于Navier－Stokes方程有限元数值求解方面的研究已有很多的文章和专著,多数是采用有限元Galerkin算法,例见文献[1-4]．然而,由于Navier-Stokes方程在大雷诺数时有其强的非线性性和对时间土的长期依赖性,用计算机求解Navier－Stokes方程在速度和容量方面是难以承受的．为了克服这些困难,最近人们提出了有限元非线性Galerkin算法,见文献卜8］,然而这种算法只是在某一有限时刻之后具有好的收敛速度,在初始时刻的某一区间不能达到好的收敛速度．本文应用Taylor展开技术导出了数值求解二维非定常Navier－Stokes方程的最佳…     

各向异性网格下抛物方程一个新的非协调混合元收敛性分析

本文将 Crouzeix-Raviart 型非协调线性三角形元应用到抛物方程,建立了一个新的混合元格式.在抛弃传统有限元分析的必要工具 Ritz 投影算子的前提下,直接利用单元的插值性质和导数转移技巧, 分别得到了各向异性剖分下关于原始变量u 的H^-1-模和积分意义下L²-模以及通量p=-▽u 在L²-模下的最优阶误差估计.数值结果与我们的理论分析是相吻合的.     

球面稳定同伦群中非平凡元素γ_sξ_n

当p≥5, n≥0时,(i_1i_0)_*(h_n)∈Ext_■~(1,p~nq)(H~*K,Z_p)在Adams谱序列中是永久循环,并且收敛到π_(p~nq-1)K中的非零元.本文在此基础上,考虑了涉及第三希腊字母类乘积元素的收敛性,并且扩大了球面稳定同伦群中非平凡元素滤子s+1的取值范围,即当p+1 s+1 2p时,■_sh_n∈Ext_■~(s+1,t)(Z_p,Z_p)在Adams谱序列中是永久循环,并且收敛到π_(t-s-1)S中的非零元γ_sξ_n,其中p≥7, n≥3, t=p~nq+sp~2q+(s-1)pq+(s-2)q+s-3,q=2(p-1).     

强阻尼波动方程的H1-Galerkin混合有限元超收敛分析

研究了强阻尼波动方程的H¹-Galerkin混合有限元方法的超收敛性. 借助于协调线性三角形元已有的分析估计式, 直接利用插值算子代替原始变量 u 的 Ritz 投影和应力变量 p 的 Ritz-Volterra 投影,对半离散和全离散格式, 得到了u在 H¹(Ω) 模和 p 在 H(div;Ω) 模意义下比以往文献高一阶的超逼近和超收敛结果.     

一个带非线性边界条件的强耦合抛物方程组的整体存在性与爆破

本文处理带非线性边界条件 u n=uα, v n=vβ ,(x ,t) ∈ Ω× (0 ,T)的抛物方程组ut =vpΔu ,vt=uqΔv ,(x ,t) ∈Ω× (0 ,T) ,其中Ω RN 为一个有界区域 ,p ,q>0和α ,β≥ 0为常数 .研究了上述问题正解的整体存在性和爆破 ,建立了整体存在和爆破的新标准 .证明了当max{p+β,q+α}≤ 1时正解 (u ,v)整体存在 ,当min{p+β ,q+α}>1且max{α ,β}<1时正解 (u ,v)在有限时刻爆破     

双曲空间中完备子流形的特征值估计

研究了$(n+p)$维双曲空间$\mathbb{H}^{n+p}$中完备非紧子流形的第一特征值的上界.特别地,证明了$\mathbb{H}^{n+p}$中具有平行平均曲率向量$H$和无迹第二基本形式有限$L^q(q\geq n)$范数的完备子流形的第一特征值不超过$\frac{(n-1)^2(1-|H|^2)}{4}$,和$\mathbb{H}^{n+1}(n\leq5)$中具有常平均曲率向量$H$和无迹第二基本形式有限$L^q(2(1-\sqrt{\frac{2}{n}})相似文献   

Exact controllability of the wave equation with Neumann boundary control

We consider the wave equation defined on a smooth bounded domainR ⁿ with boundary =₀₁, with ₀ possibly empty and ₁ nonempty and relatively open in . The control action is exercised in the Neumann boundary conditions only on ₁, while homogeneous boundary conditions of Dirichlet type are imposed on the complementary part ₀. We study by a direct method (i.e., without passing through uniform stabilization) the problem of exact controllability on some finite time interval [0,T] for initial data on some preassigned spaceZ=Z ₁ ×Z ₂ based on and with control functions in some preassigned space based on ₁ and [0,T]. We consider several choices of pairs [Z, ] of spaces, and others may be likewise studied by similar methods. Our main results are exact controllability results in the following cases: (i) and and , both under suitable geometrical conditions on the triplet {, ₀, ₁} expressed in terms of a general vector field; (iii)Z = L ² ()×[H ¹ ()] in the Neumann case ₀=Ø in the absence of geometrical conditions on , but with a special classV of controls, larger thanL ² (). The key technical issues are, in all cases, lower bounds on theL ² ( ₁)-norm of appropriate traces of the solution to the corresponding homogeneous problem. These are obtained by multiplier techniques.This paper was presented by the second named author at the IFIP WG 7.2 Conference on Boundary Control and Boundary Variations held at the Department de Mathématiques, Universite de Nice, France, June 10–13, 1986 (a preliminary version will appear as Lectures Notes in Control Sciences, edited by J. P. Zolezio [T4]); at the International Conference on Control of Distributed Parameter Systems, held at Vorau, Austria, July 6–12, 1986; at the Second Workshop on Control of Systems Governed by Partial Differential Equations, held at Val David, Quebec, Canada, October 5–9, 1986; and at the 26th IEEE Conference on Decision and Control, held at Los Angeles, CA, December 9–11, 1987 [LT6]. This research was partially supported by the National Science Foundation under Grant No. NSF 8301668 and by the Air Force Office of Scientific Research under Grant No. AFOSR-84-0365.     

高波数问题的超收敛性

 本文考虑求解Helmholtz方程的有限元方法的超逼近性质以及基于PPR后处理方法的超收敛性质.我们首先给出了矩形网格上的p-次元在收敛条件k（kh）^2p+1≤C₀下的有限元解和基于Lobatto点的有限元插值之间的超逼近以及重构的有限元梯度和精确解之间的超收敛分析.然后我们给出了四边形网格上的线性有限元方法的分析.这些估计都给出了与波数k和网格尺寸h的依赖关系.同时我们回顾了三角形网格上的线性有限元的超收敛结果.最后我们给出了数值实验并且结合Richardson外推进一步减少了误差.