Analysis of a continuous finite element method for H(\mathrm{curl},\mathrm{div})-elliptic interface problem

In this paper, we develop a continuous finite element method for the curlcurl-grad div vector second-order elliptic problem in a three-dimensional polyhedral domain occupied by discontinuous nonhomogeneous anisotropic materials. In spite of the fact that the curlcurl-grad div interface problem is closely related to the elliptic interface problem of vector Laplace operator type, the continuous finite element discretization of the standard variational problem of the former generally fails to give a correct solution even in the case of homogeneous media whenever the physical domain has reentrant corners and edges. To discretize the curlcurl-grad div interface problem by the continuous finite element method, we apply an element-local $L^2$ projector to the curl operator and a pseudo-local $L^2$ projector to the div operator, where the continuous Lagrange linear element enriched by suitable element and face bubbles may be employed. It is shown that the finite element problem retains the same coercivity property as the continuous problem. An error estimate $\mathcal{O }(h^r)$ in an energy norm is obtained between the analytical solution and the continuous finite element solution, where the analytical solution is in $\prod _{l=1}^L (H^r(\Omega _l))^3$ for some $r\in (1/2,1]$ due to the domain boundary reentrant corners and edges (e.g., nonconvex polyhedron) and due to the interfaces between the different material domains in $\Omega =\bigcup _{l=1}^L \Omega _l$ .     

Uniformly convergent H(div)-conforming rectangular elements for Darcy-Stokes problem

We consider a singular perturbation problem which describes 2D Darcy-Stokes flow.An H(div)-conforming rectangular element,DS-R14,is proposed and analyzed frst.This element has 14 degrees of freedom for velocity and is proved to be uniformly convergent with respect to perturbation constant.We then simplify this element to get another H(div)-conforming rectangular element,DS-R12,which has 12 degrees of freedom for velocity.The uniform convergence is also obtained for this element.Finally,we construct a de Rham complex corresponding to DS-R12 element.     

Multigrid in H (div) and H (curl)

Summary. We consider the solution of systems of linear algebraic equations which arise from the finite element discretization of variational problems posed in the Hilbert spaces and in three dimensions. We show that if appropriate finite element spaces and appropriate additive or multiplicative Schwarz smoothers are used, then the multigrid V-cycle is an efficient solver and preconditioner for the discrete operator. All results are uniform with respect to the mesh size, the number of mesh levels, and weights on the two terms in the inner products. Received June 12, 1998 / Revised version received March 12, 1999 / Published online January 27, 2000     

Volumes for $${mathrm{SL}}_N({mathbb {R}})$$ , the Selberg Integral and Random Lattices

Foundations of Computational Mathematics - There is a natural left and right invariant Haar measure associated with the matrix groups GL $${}_N(mathbb {R})$$ and SL $${}_N(mathbb {R})$$ due to...     

广义Witt代数$W(2, 1)$的投射不可分解模和Cartan不变理

研究了特征为2的代数闭域上广义Witt代数$W(2,\mathbf{1})$的 投射不可分解模, 给出了特征标高度$ht\chi\leq 0$的所有投射不可分解模同构类的代表元和Cartan不变量. 并且进一步讨论了既约包络代数$u(W(2,\mathbf{1}),\chi)$ 的表示型.     

Sparsity optimized high order finite element functions for H(div) on simplices

This paper deals with conforming high-order finite element discretizations of the vector-valued function space H(Div) in 2 and 3 dimensions. A new set of hierarchic basis functions on simplices with the following two main properties is introduced. When working with an affine, simplicial triangulation, first, the divergence of the basis functions is L ₂-orthogonal, and secondly, the L ₂-inner product of the interior basis functions is sparse with respect to the polynomial degree p. The construction relies on a tensor-product based construction with properly weighted Jacobi polynomials as well as an explicit splitting of the higher-order basis functions into solenoidal and non-solenoidal ones. The basis is suited for fast assembling strategies. The general proof of the sparsity result is done by the assistance of computer algebra software tools. Several numerical experiments document the proved sparsity patterns and practically achieved condition numbers for the parameter-dependent Div - Div problem. Even on curved elements a tremendous improvement in condition numbers is observed. The precomputed mass and stiffness matrix entries in general form are available online.     

Parameter-free H(div) preconditioning for a mixed finite element formulation of diffusion problems

^** Email: c.powell{at}manchester.ac.uk Mixed finite element formulations of generalised diffusion problemsyield linear systems with ill-conditioned, symmetric and indefinitecoefficient matrices. Preconditioners with optimal work complexitythat do not rely on artificial parameters are essential. Weimplement lowest order Raviart–Thomas elements and analysepractical issues associated with so-called ‘H(div) preconditioning’.Properties of the exact scheme are discussed in Powell &Silvester (2003, SIAM J. Matrix Anal. Appl., 25, 718–738).We extend the discussion, here, to practical implementation,the components of which are any available multilevel solverfor a weighted H(div) operator and a pressure mass matrix. Anew bound is established for the eigenvalue spectrum of thepreconditioned system matrix and extensive numerical resultsare presented.     

Multigrid methods for H(div) in three dimensions with nonoverlapping domain decomposition smoothers

We design and analyze V‐cycle multigrid methods for an H(div) problem discretized by the lowest‐order Raviart–Thomas hexahedral element. The smoothers in the multigrid methods involve nonoverlapping domain decomposition preconditioners that are based on substructuring. We prove uniform convergence of the V‐cycle methods on bounded convex hexahedral domains (rectangular boxes). Numerical experiments that support the theory are also presented.     

单圈图H(p,tK_{1,m})的Laplacian谱刻画

设图\,$H(p,tK_{1,m})$\,是一个顶点数为\,$p+mt$\,的连通单圈图,它是由圈\,$C_{p}$\,的依次相邻的\,$t(1\leq t\leq p)$\,个顶点、每一个顶点分别与星\,$K_{1,m}$\,的中心重合而得到的单圈图. 证明了单圈图\,$ H( p,p K_{1,4})$, $H(p,p K_{1,3})$, $H(p,(p-1)K_{1,3})$\,是由它们的\,Laplacian\,谱确定的,并证明了当\,$p$\,为偶数时,单圈图\,$H(p,$2K_{1,3})$, $H( p,(p-2) K_{1,3})$, $H(p,(p-3)K_{1,3})$\,也是由它们的\,Laplacian\,谱确定的.     

SUPG and grad‐div stabilized finite element methods for steady weakly compressible viscous flow

The extensive application of mathematical and computational methods has become an efficient and powerful approach to the investigation and solution of many problems and processes in fluid dynamics from qualitative as well as quantitative point of view. In this work a new class of advanced numerical approximation schemes to isothermal compressible viscous flow is presented. The schemes are based on an iteration between an Oseen like problem for the velocity and a hyperbolic transport equation for the density. Such schemes seem attractive for computations because they offer a reduction to simpler problems for which highly refined numerical methods either are known or can be built from existing approximation schemes to similar equations, and because of the guidance that can be drawn from an existence theory based on them. For the generalized Oseen subproblem a Taylor–Hood finite element method is proposed that is stabilized by a reduced SUPG and grad‐div technique (cf. [1, 4]) in the convection‐dominated case. Results of theoretical investigations and numerical studies are presented. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

( q,t)-analogues and $GL_{n}({\ma thbb{F}}_{ q})$

We start with a (q,t)-generalization of a binomial coefficient. It can be viewed as a polynomial in t that depends upon an integer q, with combinatorial interpretations when q is a positive integer, and algebraic interpretations when q is the order of a finite field. These (q,t)-binomial coefficients and their interpretations generalize further in two directions, one relating to column-strict tableaux and Macdonald’s “7 ^th variation” of Schur functions, the other relating to permutation statistics and Hilbert series from the invariant theory of GL_n(\mathbbF_q)GL_{n}({\mathbb{F}}_{q}) .     

Dilation operators and integral operators on amalgam space $$(L_{p},l_{q})$$(Lp,lq)

This paper establishes the Hardy–Littlewood–Pólya inequalities, the Hardy inequalities and the Hilbert inequalities on amalgam spaces. Moreover, it also gives the mapping properties of the Mellin convolutions, the Hadamard fractional integrals and the Hausdorff operators on amalgam spaces. We establish these properties by some estimates for the operator norms of the dilation operators on amalgam spaces.